The role of cytokines in acute and chronic postsurgical pain after major musculoskeletal surgeries in a quaternary pediatric center.

STUDY OBJECTIVE: To determine if baseline cytokines/chemokines and their changes over postoperative days 0-2 (POD0-2) predict acute and chronic postsurgical pain (CPSP) after major surgery. DESIGN: Prospective, observational, longitudinal nested study. SETTING: University-affiliated quaternary children's hospital. PATIENTS: Subjects (≥8 years old) with idiopathic scoliosis undergoing spine fusion or pectus excavatum undergoing Nuss procedure. MEASUREMENTS: Demographics, surgical, psychosocial measures, pain scores, and opioid use over POD0-2 were collected. Cytokine concentrations were analyzed in serial blood samples collected before and up to two weeks after surgery, using Luminex bead arrays. After data preparation, relationships between pre- and post-surgical cytokine concentrations with acute (% time in moderate-severe pain over POD0-2) and chronic (pain score > 3/10 beyond 3 months post-surgery) post-surgical pain were analyzed using univariable and multivariable regression analyses with adjustment for covariates and mixed effects models were used to associate longitudinal cytokine concentrations with pain outcomes. MAIN RESULTS: Analyses included 3,164 repeated measures of 16 cytokines/chemokines from 112 subjects (median age 15.3, IQR 13.5-17.0, 54.5 % female, 59.8 % pectus). Acute postsurgical pain was associated with higher baseline concentrations of GM-CSF (ß = 0.95, SE 0.31; p = 0.003), IL-1ß (ß = 0.84, SE 0.36; p = 0.02), IL-2 (ß = 0.78, SE 0.34; p = 0.03), and IL-12 p70 (ß = 0.88, SE 0.40; p = 0.03) and longitudinal postoperative elevations in GM-CSF (ß = 1.38, SE 0.57; p = 0.03), IFNγ (ß = 1.36, SE 0.6; p = 0.03), IL-1ß (ß = 1.25, SE 0.59; p = 0.03), IL-7 (ß = 1.65, SE 0.7; p = 0.02), and IL-12 p70 (ß = 1.17, SE 0.58; p = 0.04). In contrast, CPSP was associated with lower baseline concentration of IL-8 (ß = -0.39, SE 0.17; p = 0.02), and the risk of developing CPSP was elevated in patients with lower longitudinal postoperative concentrations of IL-6 (ß = -0.57, SE 0.26; p = 0.03), IL-8 (ß = -0.68, SE 0.24; p = 0.006), and IL-13 (ß = -0.48, SE 0.22; p = 0.03). Covariates female (vs. male) sex and surgery type (pectus surgery vs. spine) were associated with higher odds for CPSP in baseline adjusted cytokine-CPSP association models for IL-2, IL-4, IL-5, IL-6, IL-8, IL-10, TNFα, and IL-8, IL-10, respectively. CONCLUSION: We identified pro-inflammatory cytokine profiles associated with higher risk of acute postoperative pain. Interestingly, pleiotropic cytokine IL-6, chemokine IL-8 (which promotes neutrophil infiltration and monocyte differentiation), and monocyte-released anti-inflammatory cytokine IL-13, were associated with lower CPSP risk. Our results suggest heterogenous outcomes of cytokine/chemokine signaling that can both promote and protect against post-surgical pain. These may serve as predictive and prognostic biomarkers of pain outcomes following surgery.

The Role of Cytokines in Acute and Chronic Postsurgical Pain in Pediatric Patients after Major Musculoskeletal Surgeries.

Study Objective: To determine if baseline cytokines and their changes over postoperative days 0-2 (POD0-2) predict acute and chronic postsurgical pain (CPSP) after major surgery. Design: Prospective, observational, longitudinal nested study. Setting: University-affiliated quaternary children's hospital. Patients: Subjects (≥8 years old) with idiopathic scoliosis undergoing spine fusion or pectus excavatum undergoing Nuss procedure. Measurements: Demographics, surgical, psychosocial measures, pain scores, and opioid use over POD0-2 were collected. Cytokine concentrations were analyzed in serial blood samples collected before and after (up to two weeks) surgery, using Luminex bead arrays. After data preparation, relationships between pre- and post-surgical cytokine concentrations with acute (% time in moderate-severe pain over POD0-2) and chronic (pain score>3/10 beyond 3 months post-surgery) pain were analyzed. After adjusting for covariates, univariate/multivariate regression analyses were conducted to associate baseline cytokine concentrations with postoperative pain, and mixed effects models were used to associate longitudinal cytokine concentrations with pain outcomes. Main Results: Analyses included 3,164 measures of 16 cytokines from 112 subjects (median age 15.3, IQR 13.5-17.0, 54.5% female, 59.8% pectus). Acute postsurgical pain was associated with higher baseline concentrations of GM-CSF (ß=0.95, SE 0.31; p=.003), IL-1ß (ß=0.84, SE 0.36; p=.02), IL-2 (ß=0.78, SE 0.34; p=.03), and IL-12 p70 (ß=0.88, SE 0.40; p=.03) and longitudinal postoperative elevations in GM-CSF (ß=1.38, SE 0.57; p=.03), IFNγ (ß=1.36, SE 0.6; p=.03), IL-1ß (ß=1.25, SE 0.59; p=.03), IL-7 (ß=1.65, SE 0.7, p=.02), and IL-12 p70 (ß=1.17, SE 0.58; p=.04). In contrast, CPSP was associated with lower baseline concentration of IL-8 (ß= -0.39, SE 0.17; p=.02), and the risk of developing CPSP was elevated in patients with lower longitudinal postoperative concentrations of IL-6 (ß= -0.57, SE 0.26; p=.03), IL-8 (ß= -0.68, SE 0.24; p=.006), and IL-13 (ß= -0.48, SE 0.22; p=.03). Furthermore, higher odds for CPSP were found for females (vs. males) for IL-2, IL-4, IL-5, IL-6, IL-8, IL-10, and TNFα, and for pectus (vs. spine) surgery for IL-8 and IL-10. Conclusion: We identified pro-inflammatory cytokines associated with increased acute postoperative pain and anti-inflammatory cytokines associated with lower CPSP risk, with potential to serve as predictive and prognostic biomarkers.

Systems Biology Guided Gene Enrichment Approaches Improve Prediction of Chronic Post-surgical Pain After Spine Fusion.

OBJECTIVES: Incorporation of genetic factors in psychosocial/perioperative models for predicting chronic postsurgical pain (CPSP) is key for personalization of analgesia. However, single variant associations with CPSP have small effect sizes, making polygenic risk assessment important. Unfortunately, pediatric CPSP studies are not sufficiently powered for unbiased genome wide association (GWAS). We previously leveraged systems biology to identify candidate genes associated with CPSP. The goal of this study was to use systems biology prioritized gene enrichment to generate polygenic risk scores (PRS) for improved prediction of CPSP in a prospectively enrolled clinical cohort. METHODS: In a prospectively recruited cohort of 171 adolescents (14.5 ± 1.8 years, 75.4% female) undergoing spine fusion, we collected data about anesthesia/surgical factors, childhood anxiety sensitivity (CASI), acute pain/opioid use, pain outcomes 6-12 months post-surgery and blood (for DNA extraction/genotyping). We previously prioritized candidate genes using computational approaches based on similarity for functional annotations with a literature-derived "training set." In this study, we tested ranked deciles of 1336 prioritized genes for increased representation of variants associated with CPSP, compared to 10,000 randomly selected control sets. Penalized regression (LASSO) was used to select final variants from enriched variant sets for calculation of PRS. PRS incorporated regression models were compared with previously published non-genetic models for predictive accuracy. RESULTS: Incidence of CPSP in the prospective cohort was 40.4%. 33,104 case and 252,590 control variants were included for association analyses. The smallest gene set enriched for CPSP had 80/1010 variants associated with CPSP (p < 0.05), significantly higher than in 10,000 randomly selected control sets (p = 0.0004). LASSO selected 20 variants for calculating weighted PRS. Model adjusted for covariates including PRS had AUROC of 0.96 (95% CI: 0.92-0.99) for CPSP prediction, compared to 0.70 (95% CI: 0.59-0.82) for non-genetic model (p < 0.001). Odds ratios and positive regression coefficients for the final model were internally validated using bootstrapping: PRS [OR 1.98 (95% CI: 1.21-3.22); ß 0.68 (95% CI: 0.19-0.74)] and CASI [OR 1.33 (95% CI: 1.03-1.72); ß 0.29 (0.03-0.38)]. DISCUSSION: Systems biology guided PRS improved predictive accuracy of CPSP risk in a pediatric cohort. They have potential to serve as biomarkers to guide risk stratification and tailored prevention. Findings highlight systems biology approaches for deriving PRS for phenotypes in cohorts less amenable to large scale GWAS.

NAT1 genetic variation increases asthma risk in children with secondhand smoke exposure.

OBJECTIVE: We previously reported that children exposed to secondhand smoke (SHS) that carried variants in the NAT1 gene had over two-fold higher hair cotinine levels. Our objective was to determine if NAT1 polymorphisms confer increased risk for developing asthma in children exposed to SHS. METHODS: White participants in the Cincinnati Childhood Allergy and Air Pollution Study (n = 359) were genotyped for 10 NAT1 variants. Smoke exposure was defined by hair cotinine and parental report. Asthma was objectively assessed by spirometry and methacholine challenge. Findings were replicated in the Genomic Control Cohort (n = 638). RESULTS: Significant associations between 5 NAT1 variants and asthma were observed in the CCAAPS exposed group compared to none in the unexposed group. There was a significant interaction between NAT1 rs13253389 and rs4921581 with smoke exposure (p = 0.02, p = 0.01) and hair cotinine level (p = 0.048, p = 0.042). Children wildtype for rs4921581 had increasing asthma risk with increasing hair cotinine level, whereas those carrying the NAT1 minor allele had an increased risk of asthma regardless of cotinine level. In the GCC, 13 NAT1 variants were associated with asthma in the smoke-exposed group, compared to 0 in the unexposed group, demonstrating gene-level replication. CONCLUSIONS: Variation in the NAT1 gene modifies asthma risk in children exposed to secondhand-smoke. To our knowledge, this is the first report of a gene-environment interaction between NAT1 variants, smoke exposure, cotinine levels, and pediatric asthma. NAT1 genotype may have clinical utility as a biomarker of increased asthma risk in children exposed to smoke.

Hydromorphone population pharmacokinetics in pediatric surgical patients.

BACKGROUND: Hydromorphone is an opioid agonist used for pediatric analgesia. Due to lack of data, pediatric dosing (based on adult pharmacokinetic models) is not optimal. AIM: This study characterizes hydromorphone population pharmacokinetics in pediatric surgical patients. METHODS: In this prospective observational study, 34 children (4-18 years, bodyweight 23-89.6 kg) received multiple intravenous hydromorphone boluses followed by postoperative hydromorphone patient-controlled analgesia. Arterial blood samples were collected before and at 3, 10, 30, and 90 (and few samples at 1350) minutes after the first dose. Hydromorphone concentrations were measured by validated LC-MS/MS assay. Nonlinear mixed-effects modeling was used for pharmacokinetic model development. The final population pharmacokinetic model was evaluated by visual predictive check and bootstrap analysis. Monte Carlo simulations based on the final pharmacokinetic model determined optimal patient-controlled analgesia parameters to achieve a target of 20 ng/mL (as the median effective analgesic concentration), using minimum effective analgesic concentration of 4 ng/mL as a proxy for patient-controlled analgesia dose demand, and not exceeding the defined safe upper threshold of 40 ng/mL. RESULTS: Hydromorphone pharmacokinetic profiles were adequately described by a two-compartmental model with first-order elimination. Bodyweight was found to be a significant covariate for hydromorphone clearance. Allometrically scaledpharmacokinetic parameter estimates (per 70 kg), systemic clearance (0.748 L/min), volume of distribution (33 L), peripheral clearance (1.57 L/min), and peripheral volume of distribution (146 L) were similar to reported adult parameter estimates. Sex, race, age, and type of surgery were not identified as significant covariates. To identify optimal patient-controlled analgesia dosing parameters, we simulated several initial loading doses, demand doses, and lockout intervals. Our simulations support an initial patient-controlled analgesia loading dose of 15 µg/kg followed by a demand dose of 6 µg/kg with lockout intervals of 20 minutes. CONCLUSIONS: After intravenous hydromorphone, plasma pharmacokinetic profiles in children undergoing different surgeries were well described by a two-compartment population allometric pharmacokinetic model using bodyweight as the size descriptor. Model informed simulations identified patient-controlled analgesia parameters to inform initial settings, with adjustments as needed based on observed individual effects.

Prevalence of Noonan spectrum disorders in a pediatric population with valvar pulmonary stenosis.

OBJECTIVE: To evaluate the prevalence of Noonan spectrum disorders (NSD) in a pediatric population with valvar pulmonary stenosis (vPS) and identify the clinical characteristics that differentiate those with NSD from those without NSD. DESIGN: A retrospective chart review of 204 patients diagnosed with vPS between 9/1/2012 and 12/1/2016 at a pediatric medical center was performed. The quantitative features of vPS, genetic diagnosis information, and phenotypic characteristics of Noonan syndrome were collected. Chi-square test, Fisher's exact test, t test, Wilcoxon rank-sum test, and ANOVA were used for comparisons among the groups. Logistic regression was used to test for the association between the clinical characteristics and the presence of NSD. RESULTS: Syndromic diagnoses were made in 10% of the children with vPS, with NSD accounting for 6%. Hypertrophic cardiomyopathy (P < .0001), short stature (P < .0001), developmental delay (P < .0001), ophthalmological abnormalities (P < .0001), pectus carinatum/excavatum (P = .01), neurological abnormalities (P = .022), and aortic stenosis (P = .031) were present more often in individuals with NSD compared to nonsyndromic vPS. A logistic regression analysis showed a 4.8-fold increase in odds for NSD for each additional characteristic (P < .0001). CONCLUSIONS: At least 6% of the children with vPS have an underlying NSD. Individuals with vPS and NSD were significantly more likely to have additional features known to be associated with NSD than those with vPS without NSD. We conclude that vPS in the presence of one or more significant characteristics should prompt referral for genetic evaluation as a guide to ascertain patients at risk for NSD while optimizing the use of clinical genetics evaluation and potential genetic testing.

Posterior fossa anomalies diagnosed with fetal MRI: associated anomalies and neurodevelopmental outcomes.

OBJECTIVE: The purpose of this study was to describe the relationship between intracranial and extracranial anomalies and neurodevelopmental outcome for fetuses diagnosed with a posterior fossa anomaly (PFA) on fetal MRI. METHODS: Cases of Dandy-Walker malformation, vermian hypogenesis/hypoplasia, and mega cisterna magna (MCM) were identified through the Fetal Care Center of Cincinnati between January 2004 and December 2010. Parental interview and retrospective chart review were used to assess neurodevelopmental outcome. RESULTS: Posterior fossa anomalies were identified in 59 fetuses; 9 with Dandy-Walker malformation, 36 with vermian hypogenesis/hypoplasia, and 14 with MCM. Cases with isolated PFAs (14/59) had better outcomes than those with additional anomalies (p = 0.00016), with isolated cases of MCM all being neurodevelopmentally normal. Cases with additional intracranial anomalies had a worse outcome than those without intracranial anomalies (p = 0.00017). The presence of extracranial anomalies increased the likelihood of having a poor outcome (p = 0.00014) as did the identification of an abnormal brainstem (p = 0.00018). CONCLUSION: Intracranial and extracranial anomalies were good predictors of neurodevelopmental outcome in this study. The prognosis was poor for individuals with an abnormal brainstem, whereas those with isolated MCM had normal neurodevelopmental outcome.

Genomic organization and expression analysis of the murine Fam3c gene.

Previously, we identified FAM3C as a candidate gene for autosomal recessive nonsyndromic hearing loss locus 17 (DFNB17). This gene has since been found to be a member of a cytokine-like gene family, but its function has not been determined. The purpose of this study was thus to elucidate the gene structure and pattern of expression, providing information that might allow a hypothesis to be developed about FAM3C function of in the inner ear. To do this we analyzed its mouse ortholog, Fam3c. Fam3c was found to be ubiquitously expressed in all analyzed tissues, and had two major transcript variants presumed to result from an alternative use of two distinct polyadenylation signals. In situ hybridization experiments revealed a predominant Fam3c pattern of expression in the nonsensory epithelium of the growing semicircular canals at embryonic day (E) 15.5. This expression pattern resembles the known pattern of the Nkx5 homeobox genes. Analysis of the Fam3c promoter region demonstrated a putative Nkx5.1 binding site. Based on our findings, we hypothesize that Fam3c may be a downstream target gene for the Nkx5.1 transcription factor, and may thus be involved in cell differentiation and proliferation during inner ear embryogenesis. Additionally, analyses of putative amino acid sequences of FAM3C orthologous proteins showed that their primary and secondary structures and overall topology were highly conserved. Further study is underway to determine the role of FAM3C in inner ear development.

Genotypic and phenotypic correlations of DFNB1-related hearing impairment in the Midwestern United States.

OBJECTIVE: To determine the genotypic and phenotypic correlations of hearing impairment (HI) in a midwestern US population related to autosomal recessive nonsyndromic hearing loss locus 1 (DFNB1). DESIGN: A retrospective review. SETTING: Tertiary care children's hospital. PATIENTS: A total of 160 consecutive children diagnosed with idiopathic sensorineural hearing loss. MAIN OUTCOME MEASURES: GJB2 genotype and audiometric phenotype. RESULTS: The prevalence of subjects with HI having biallelic GJB2-related mutations was 15.3% (24/157). Of these 24 patients, 9 (38%) were homozygous 35delG, 6 (25%) had other biallelic nonsense mutations, and 9 (38%) had a missense mutation of at least 1 allele. The allelic prevalence of 35delG was 8.6% (27/314) in the study population and 48% (23/48) in the DFNB1 group. The M34T allele mutation was next most prevalent at 2.2% (7/314) in the study population and 10% (5/48) in the DFNB1 group. Severe to profound HI occurred in 59% of DFNB1 subjects. Genotypes with biallelic nonsense mutations had a high risk of severe to profound HI (88%). DFNB1-related HI was usually bilateral, symmetric, nonprogressive, and had flat audiograms. However, asymmetric HI (22%), sloping audiograms (26%), and even borderline-normal hearing in 1 ear was observed, and these were associated with the presence of at least 1 missense mutation. Two novel mutations, K15T and L90V, were identified. A subject presenting to our clinic with severe to profound HI had a 40% risk of biallelic GJB2 mutation. CONCLUSIONS: Our population represents a consecutively enrolled clinic population with sensorineural hearing loss. In our DFNB1-related HI cohort, the 35delG mutation and severe to profound HI rates were lower than previously reported. Our missense mutation and M34T allelic prevalence rates were higher than expected and were associated with a less severe hearing loss. The presence of biallelic nonsense mutations was associated with severe to profound hearing loss in nearly 90% of cases. Mild asymmetric HI and sloping audiograms were more often associated with missense mutations.