Longitudinal changes in DNA methylation during the onset of islet autoimmunity differentiate between reversion versus progression of islet autoimmunity.

Background: Type 1 diabetes (T1D) is preceded by a heterogenous pre-clinical phase, islet autoimmunity (IA). We aimed to identify pre vs. post-IA seroconversion (SV) changes in DNAm that differed across three IA progression phenotypes, those who lose autoantibodies (reverters), progress to clinical T1D (progressors), or maintain autoantibody levels (maintainers). Methods: This epigenome-wide association study (EWAS) included longitudinal DNAm measurements in blood (Illumina 450K and EPIC) from participants in Diabetes Autoimmunity Study in the Young (DAISY) who developed IA, one or more islet autoantibodies on at least two consecutive visits. We compared reverters - individuals who sero-reverted, negative for all autoantibodies on at least two consecutive visits and did not develop T1D (n=41); maintainers - continued to test positive for autoantibodies but did not develop T1D (n=60); progressors - developed clinical T1D (n=42). DNAm data were measured before (pre-SV visit) and after IA (post-SV visit). Linear mixed models were used to test for differences in pre- vs post-SV changes in DNAm across the three groups. Linear mixed models were also used to test for group differences in average DNAm. Cell proportions, age, and sex were adjusted for in all models. Median follow-up across all participants was 15.5 yrs. (interquartile range (IQR): 10.8-18.7). Results: The median age at the pre-SV visit was 2.2 yrs. (IQR: 0.8-5.3) in progressors, compared to 6.0 yrs. (IQR: 1.3-8.4) in reverters, and 5.7 yrs. (IQR: 1.4-9.7) in maintainers. Median time between the visits was similar in reverters 1.4 yrs. (IQR: 1-1.9), maintainers 1.3 yrs. (IQR: 1.0-2.0), and progressors 1.8 yrs. (IQR: 1.0-2.0). Changes in DNAm, pre- vs post-SV, differed across the groups at one site (cg16066195) and 11 regions. Average DNAm (mean of pre- and post-SV) differed across 22 regions. Conclusion: Differentially changing DNAm regions were located in genomic areas related to beta cell function, immune cell differentiation, and immune cell function.

Pediatric and Adult Fracture Incidence: A Decreasing Trend With Increasing Hospital Admissions.

Large-scale studies examining fracture trends and epidemiological data are lacking. The purpose of this study was to evaluate the incidence of fractures presenting to US emergency departments using the National Electronic Injury Surveillance System. A total of 7,109,078 pediatric and 13,592,548 adult patients presenting to US emergency departments with a fracture between 2008 and 2017 were analyzed for patterns. Fractures accounted for 13.9% of pediatric injuries and 15% of adult injuries. Among children, fracture incidence was highest in the group 10 to 14 years old and most frequently involved the forearm (19.0%). Fracture incidence was highest in adults 80 years and older and most frequently involved the lower trunk (16.2%). On average, the rate of pediatric fractures decreased by 2.34% per year (95% CI, 0.25% increase to 4.88% decrease; P=.0757). Among adults, fracture incidence increased 0.33% per year (95% CI, 2.34% decrease to 2.85% increase; P=.7892). This change was significantly different between the pediatric and adult populations (P=.0152). There was an increase in the annual proportion of adults with fractures who were admitted (odds ratio per 1-year increase, 1.05; 95% CI, 1.03-1.07; P<.0001). There was no change in the proportion of pediatric patients with fractures who were admitted (odds ratio, 1.02; 95% CI, 0.99-1.05; P=.0606). The incidence of fractures decreased in pediatric patients yet was relatively stable in adult patients. Conversely, the proportion of patients with fractures who were admitted increased, particularly among adults. These findings may suggest that less severe fractures are presenting elsewhere, falsely inflating the observed rise in admissions. [Orthopedics. 2023;46(6):e369-e375.].

Investigating iron intake in risk of progression from islet autoimmunity to type 1 diabetes: The diabetes autoimmunity study in the young.

Background: Studies of the role of iron in the risk of type 1 diabetes (T1D) have been inconsistent. Given that iron generates reactive oxygen radicals, which can lead to oxidative damage and apoptosis in the beta cells of the pancreas, we examined whether iron intake was associated with the risk of progressing to T1D in individuals with islet autoimmunity (IA), the pre-clinical phase of T1D. Methods: DAISY is a prospective cohort following 2,547 children at increased risk for IA and progression to T1D. IA is defined as at least two consecutive serum samples positive for at least one autoantibody (insulin, GAD, IA-2, or ZnT8). We measured dietary intake at the time of IA seroconversion in 175 children with IA, and of these, 64 progressed to T1D. We used Cox regression to examine the association between energy-adjusted iron intake and progression to T1D, adjusting for HLA-DR3/4 genotype, race/ethnicity, age at seroconversion, presence of multiple autoantibodies at seroconversion, and multiple vitamin use. In addition, we tested whether this association was modified by vitamin C or calcium intake. Results: In children with IA, high iron intake (as defined as above the 75th percentile, > 20.3 mg/day) was associated with decreased risk of progression to T1D compared to moderate iron intake (as defined by the middle 25-75th percentiles, 12.7-20.3 mg/day) (adjusted hazard ratio (HR): 0.35; 95% confidence interval (CI): 0.15, 0.79). The association between iron intake and T1D was not modified by vitamin C nor calcium intake. In a sensitivity analysis, the removal of six children who had been diagnosed with celiac disease prior to IA seroconversion did not affect this association. Conclusion: Higher iron intake at the time of IA seroconversion is associated with a lower risk of progression to T1D, independent of multivitamin supplement use. Further research that includes plasma biomarkers of iron status is needed to investigate the relationship between iron and the risk of T1D.

Genome-wide analysis of oxylipins and oxylipin profiles in a pediatric population.

Background: Oxylipins are inflammatory biomarkers derived from omega-3 and-6 fatty acids implicated in inflammatory diseases but have not been studied in a genome-wide association study (GWAS). The aim of this study was to identify genetic loci associated with oxylipins and oxylipin profiles to identify biologic pathways and therapeutic targets for oxylipins. Methods: We conducted a GWAS of plasma oxylipins in 316 participants in the Diabetes Autoimmunity Study in the Young (DAISY). DNA samples were genotyped using the TEDDY-T1D Exome array, and additional variants were imputed using the Trans-Omics for Precision Medicine (TOPMed) multi-ancestry reference panel. Principal components analysis of 36 plasma oxylipins was used to capture oxylipin profiles. PC1 represented linoleic acid (LA)- and alpha-linolenic acid (ALA)-related oxylipins, and PC2 represented arachidonic acid (ARA)-related oxylipins. Oxylipin PC1, PC2, and the top five loading oxylipins from each PC were used as outcomes in the GWAS (genome-wide significance: p < 5×10-8). Results: The SNP rs143070873 was associated with (p < 5×10-8) the LA-related oxylipin 9-HODE, and rs6444933 (downstream of CLDN11) was associated with the LA-related oxylipin 13 S-HODE. A locus between MIR1302-7 and LOC100131146, rs10118380 and an intronic variant in TRPM3 were associated with the ARA-related oxylipin 11-HETE. These loci are involved in inflammatory signaling cascades and interact with PLA2, an initial step to oxylipin biosynthesis. Conclusion: Genetic loci involved in inflammation and oxylipin metabolism are associated with oxylipin levels.

Epidemiology of Physeal Fractures and Clinically Significant Growth Disturbances Affecting the Distal Tibia, Proximal Tibia, and Distal Femur: A Retrospective Cohort Study.

INTRODUCTION: Childhood fractures involving the physis potentially result in premature physeal closure that can lead to growth disturbances. Growth disturbances are challenging to treat with associated complications. Current literature focusing on physeal injuries to lower extremity long bones and risk factors for growth disturbance development is limited. The purpose of this study was to provide a review of growth disturbances among proximal tibial, distal tibial, and distal femoral physeal fractures. METHODS: Data were retrospectively collected from patients undergoing fracture treatment at a level I pediatric trauma center between 2008 and 2018. The study was limited to patients 0.5 to 18.9 years with a tibial or distal femoral physeal fracture, injury radiograph, and appropriate follow-up for determination of fracture healing. The cumulative incidence of clinically significant growth disturbance (CSGD) (a growth disturbance requiring subsequent physeal bar resection, osteotomy, and/or epiphysiodesis) was estimated, and descriptive statistics were used to summarize demographics and clinical characteristics among patients with and without CSGD. RESULTS: A total of 1,585 patients met the inclusion criteria. The incidence of CSGD was 5.0% (95% confidence interval, 3.8% to 6.6%). All cases of growth disturbance occurred within 2 years of initial injury. The risk of CSGD peaked at 10.2 years for males and 9.1 years for females. Complex fractures that required surgical treatment, distal femoral and proximal tibial fractures, age, and initial treatment at an outside hospital were significantly associated with an increased risk of a CSGD. DISCUSSION: All CSGDs occurred within 2 years of injury, indicating that these injuries should be followed for a period of at least 2 years. Patients with distal femoral or proximal tibial physeal fractures that undergo surgical treatment are at highest risk for developing a CSGD. LEVEL OF EVIDENCE: Level III Retrospective Cohort Study.

Rare Coding Variants in Patients with Non-Syndromic Vestibular Dysfunction.

Vertigo due to vestibular dysfunction is rare in children. The elucidation of its etiology will improve clinical management and the quality of life of patients. Genes for vestibular dysfunction were previously identified in patients with both hearing loss and vertigo. This study aimed to identify rare, coding variants in children with peripheral vertigo but no hearing loss, and in patients with potentially overlapping phenotypes, namely, Meniere's disease or idiopathic scoliosis. Rare variants were selected from the exome sequence data of 5 American children with vertigo, 226 Spanish patients with Meniere's disease, and 38 European-American probands with scoliosis. In children with vertigo, 17 variants were found in 15 genes involved in migraine, musculoskeletal phenotypes, and vestibular development. Three genes, OTOP1, HMX3, and LAMA2, have knockout mouse models for vestibular dysfunction. Moreover, HMX3 and LAMA2 were expressed in human vestibular tissues. Rare variants within ECM1, OTOP1, and OTOP2 were each identified in three adult patients with Meniere's disease. Additionally, an OTOP1 variant was identified in 11 adolescents with lateral semicircular canal asymmetry, 10 of whom have scoliosis. We hypothesize that peripheral vestibular dysfunction in children may be due to multiple rare variants within genes that are involved in the inner ear structure, migraine, and musculoskeletal disease.

Propensity scores as a novel method to guide sample allocation and minimize batch effects during the design of high throughput experiments.

BACKGROUND: We developed a novel approach to minimize batch effects when assigning samples to batches. Our algorithm selects a batch allocation, among all possible ways of assigning samples to batches, that minimizes differences in average propensity score between batches. This strategy was compared to randomization and stratified randomization in a case-control study (30 per group) with a covariate (case vs control, represented as ß1, set to be null) and two biologically relevant confounding variables (age, represented as ß2, and hemoglobin A1c (HbA1c), represented as ß3). Gene expression values were obtained from a publicly available dataset of expression data obtained from pancreas islet cells. Batch effects were simulated as twice the median biological variation across the gene expression dataset and were added to the publicly available dataset to simulate a batch effect condition. Bias was calculated as the absolute difference between observed betas under the batch allocation strategies and the true beta (no batch effects). Bias was also evaluated after adjustment for batch effects using ComBat as well as a linear regression model. In order to understand performance of our optimal allocation strategy under the alternative hypothesis, we also evaluated bias at a single gene associated with both age and HbA1c levels in the 'true' dataset (CAPN13 gene). RESULTS: Pre-batch correction, under the null hypothesis (ß1), maximum absolute bias and root mean square (RMS) of maximum absolute bias, were minimized using the optimal allocation strategy. Under the alternative hypothesis (ß2 and ß3 for the CAPN13 gene), maximum absolute bias and RMS of maximum absolute bias were also consistently lower using the optimal allocation strategy. ComBat and the regression batch adjustment methods performed well as the bias estimates moved towards the true values in all conditions under both the null and alternative hypotheses. Although the differences between methods were less pronounced following batch correction, estimates of bias (average and RMS) were consistently lower using the optimal allocation strategy under both the null and alternative hypotheses. CONCLUSIONS: Our algorithm provides an extremely flexible and effective method for assigning samples to batches by exploiting knowledge of covariates prior to sample allocation.

An Oxylipin-Related Nutrient Pattern and Risk of Type 1 Diabetes in the Diabetes Autoimmunity Study in the Young (DAISY).

Oxylipins, pro-inflammatory and pro-resolving lipid mediators, are associated with the risk of type 1 diabetes (T1D) and may be influenced by diet. This study aimed to develop a nutrient pattern related to oxylipin profiles and test their associations with the risk of T1D among youth. The nutrient patterns were developed with a reduced rank regression in a nested case-control study (n = 335) within the Diabetes Autoimmunity Study in the Young (DAISY), a longitudinal cohort of children at risk of T1D. The oxylipin profiles (adjusted for genetic predictors) were the response variables. The nutrient patterns were tested in the case-control study (n = 69 T1D cases, 69 controls), then validated in the DAISY cohort using a joint Cox proportional hazards model (n = 1933, including 81 T1D cases). The first nutrient pattern (NP1) was characterized by low beta cryptoxanthin, flavanone, vitamin C, total sugars and iron, and high lycopene, anthocyanidins, linoleic acid and sodium. After adjusting for T1D family history, the HLA genotype, sex and race/ethnicity, NP1 was associated with a lower risk of T1D in the nested case-control study (OR: 0.44, p = 0.0126). NP1 was not associated with the risk of T1D (HR: 0.54, p-value = 0.1829) in the full DAISY cohort. Future studies are needed to confirm the nested case-control findings and investigate the modifiable factors for oxylipins.

Prediction of Late Dysplasia Based on Ultrasound and Plain X-Ray at 6 Months.

BACKGROUND: Developmental dysplasia of the hip represents a spectrum of deformity. Residual dysplasia at 2 years of age is associated with an increased risk for osteoarthritis and functional limitations. We compared the prognostic value of 6-month imaging modalities and aimed to identify optimal diagnostic metrics for the prediction of residual dysplasia. METHODS: After IRB approval, patients who underwent Pavlik treatment between 2009 and 2018 with 2-year follow-up were identified. Sonographs [ultrasound (US)] and radiographs (x-ray) were obtained at 6-month and 2-year-old visits. Dysplasia at 2 years was defined as an acetabular index (AI) >24 degrees. Receiver operating characteristic curves were constructed to quantitatively compare the prognostic ability of US and x-ray-based measures at 6 months. Youden's index [(YI) (values range from 0 (poor test) to 1 (perfect test)] was used to evaluate existing cutoffs at 6 months of age (normal measurements: alpha angle (AA) ≥60 degrees, femoral head coverage (FHC) ≥50%, and AI <30 degrees) relative to newly proposed limits. RESULTS: Fifty-nine patients were included, of which 28.8% of patients (95% CI: 17.3 to 40.4%) had acetabular dysplasia at 2 years. After adjusting for sex, AA [Area under the Curve (AUC): 80] and AI (AUC: 79) at 6 months of age were better tests than FHC (AUC: 0.77). Current diagnostic cutoffs for AA (YI: 0.08), AI (YI: 0.0), and FHC (YI: 0.06) at 6 months had poor ability to predict dysplasia at 2 years. A composite test of all measures based on proposed cutoffs (AA ≥73 degrees, FHC > 62% and AI ≤24 degrees) was a better predictor of dysplasia at 2 years (Youden's index (YI): 0.63) than any single metric. CONCLUSIONS: The rate of residual dysplasia remains concerning. The 6-month x-ray and US both play a role in the ongoing management of the developmental dysplasia of the hip. The prediction of dysplasia is maximized when all metrics are considered collectively. Existing parameters were not accurate; We recommend the following cutoffs: AA ≥73 degrees, FHC > 62%, and AI ≤24 degrees. These cutoffs must be validated. LEVEL OF EVIDENCE: Prognostic Level II.

Application of Accepted Use Criteria for the Treatment of Developmental Dysplasia of the Hip Decreases the Number of Infants Treated With a Pavlik Harness.

BACKGROUND: This analysis examined how the application of the American Academy of Orthopedic Surgeons appropriate use criteria (AUC) for developmental dysplasia of the hip in infants would change treatment patterns and outcomes for Graf IIA hips at a single quaternary pediatric hospital. METHODS: After Institutional Review Board approval, patient medical records were reviewed and data were collected. Graf IIa hips were defined as alpha angle (AA) 50 to 59 degrees. AA and femoral head coverage (FHC) were measured from initial and 6-month ultrasounds and acetabular index (AI) was measured from radiographs at 6 months of age. Instability (positive Ortolani and Barlow tests) was noted. On the basis of the American Academy of Orthopedic Surgeons AUC for managing developmental dysplasia of the hip, hips were further categorized as normal (FHC ≥45%), borderline (FHC 35% to 44%), or dysplastic (FHC <35%). RESULTS: Overall, 13% (49/371) of Graf IIa hips (AA 50 to 59 degrees) were dysplastic (FHC <35%). Total 24% (89/371) were clinically unstable. Total 42% (37/89) of unstable Graf IIa hips were dysplastic. Only 4% of stable Graf IIa hips were dysplastic (12/282). Out of 371 Graf IIa hips, 256 were treated with Pavlik harness (n=250) or Rhino brace (n=6). Among stable, nondysplastic (SND) hips (those with normal and borderline FHC≥35%), 33% (52/158) were treated because of a more severe contralateral side. If the AUC had been applied, 67% (106/158) of SND Graf IIa hips would not have been treated. Among the n=162 hips that returned for a 6-month radiograph, there was no difference in AI in the 115 treated and 47 untreated SND hips (mean difference treatment vs. no treatment: -1.5, 95% CI, -3.1 to 0.2, P =0.0808). CONCLUSIONS: Using AUC recommendations, our center could reduce the number of SND Graf IIa hips we treat by 67%. Although 24% of Graf IIa hips were clinically unstable and 13% were dysplastic based on FHC, most Graf IIa hips had normal or borderline FHC per the AUC and may do well with observation and follow-up ultrasound at 12 weeks old. LEVEL OF EVIDENCE: Level III-diagnostic study.

DNA Methylation Near DLGAP2 May Mediate the Relationship between Family History of Type 1 Diabetes and Type 1 Diabetes Risk.

Given the differential risk of type 1 diabetes (T1D) in offspring of affected fathers versus affected mothers and our observation that T1D cases have differential DNA methylation near the imprinted DLGAP2 gene compared to controls, we examined whether methylation near DLGAP2 mediates the association between T1D family history and T1D risk. In a nested case-control study of 87 T1D cases and 87 controls from the Diabetes Autoimmunity Study in the Young, we conducted causal mediation analyses at 12 DLGAP2 region CpGs to decompose the effect of family history on T1D risk into indirect and direct effects. These effects were estimated from two regression models adjusted for the human leukocyte antigen DR3/4 genotype: a linear regression of family history on methylation (mediator model) and a logistic regression of family history and methylation on T1D (outcome model). For 8 of the 12 CpGs, we identified a significant interaction between T1D family history and methylation on T1D risk. Accounting for this interaction, we found that the increased risk of T1D for children with affected mothers compared to those with no family history was mediated through differences in methylation at two CpGs (cg27351978, cg00565786) in the DLGAP2 region, as demonstrated by a significant pure natural indirect effect (odds ratio (OR) = 1.98, 95% confidence interval (CI): 1.06-3.71) and nonsignificant total natural direct effect (OR = 1.65, 95% CI: 0.16-16.62) (for cg00565786). In contrast, the increased risk of T1D for children with an affected father or sibling was not explained by DNA methylation changes at these CpGs. Results were similar for cg27351978 and robust in sensitivity analyses. Lastly, we found that DNA methylation in the DLGAP2 region was associated (P<0:05) with gene expression of nearby protein-coding genes DLGAP2, ARHGEF10, ZNF596, and ERICH1. Results indicate that the maternal protective effect conferred through exposure to T1D in utero may operate through changes to DNA methylation that have functional downstream consequences.

Changes in the Coexpression of Innate Immunity Genes During Persistent Islet Autoimmunity Are Associated With Progression of Islet Autoimmunity: Diabetes Autoimmunity Study in the Young (DAISY).

Longitudinal changes in gene expression during islet autoimmunity (IA) may provide insight into biological processes that explain progression to type 1 diabetes (T1D). We identified individuals from Diabetes Autoimmunity Study in the Young (DAISY) who developed IA, autoantibodies present on two or more visits. Illumina's NovaSeq 6000 was used to quantify gene expression in whole blood. With linear mixed models we tested for changes in expression after IA that differed across individuals who progressed to T1D (progressors) (n = 25), reverted to an autoantibody-negative stage (reverters) (n = 47), or maintained IA positivity but did not develop T1D (maintainers) (n = 66). Weighted gene coexpression network analysis was used to identify coexpression modules. Gene Ontology pathway analysis of the top 150 differentially expressed genes (nominal P < 0.01) identified significantly enriched pathways including leukocyte activation involved in immune response, innate immune response, and regulation of immune response. We identified a module of 14 coexpressed genes with roles in the innate immunity. The hub gene, LTF, is known to have immunomodulatory properties. Another gene within the module, CAMP, is potentially relevant based on its role in promoting ß-cell survival in a murine model. Overall, results provide evidence of alterations in expression of innate immune genes prior to onset of T1D.

Lower extremity physeal fractures and post-traumatic surgical deformities; National Trauma Data Bank and institutional cohorts.

Background: Pediatric lower extremity physeal fractures carry a risk of developing deformities. Most epidemiological evidence is over 25 years old, single institution, and lacks follow-up, while recent studies report variable results. Understanding their epidemiology and deformity risk is important for patient counseling and follow-up. Methods: The National Trauma Data Bank (NTDB) from 2016 was queried to describe the modern epidemiology of physeal fractures. This was contrasted with our 10-year experience of surgically treated deformities. Basic descriptive statistics, Chi-square analysis, prevalence ratios and multivariable linear regression were used to interpret results. Results: The NTDB contained 22,048 non-physeal and 1,929 physeal fractures of the femur, tibia, and fibula. Physeal fracture prevalence rose after 8 years of age but decreased for girls 2 years sooner than boys. Salter Harris (SH) type 2 fractures predominated. Physeal fractures were more commonly associated with lower energy mechanisms of injury. Distal tibia fractures were more prevalent in the NTDB cohort, while distal femur and SH-1 fractures were more prevalent in the operative cohort. Over 10 years, only 52 (5.3%) of the deformity-correcting surgeries at our institution were for physeal fracture sequelae. Age at injury and intraarticular fractures were associated with shorter times from injury to deformity correction. Conclusion: Lower extremity physeal fractures are uncommon. Fracture pattern prevalence differs from an operative cohort. Proximal tibia physeal fractures appear to be an underappreciated source of deformity. The risk of developing deformity requiring operative intervention appears to be low and is generally treated within 2 years of initial injury.

An 8-Week Neuromuscular Training Program After Concussion Reduces 1-Year Subsequent Injury Risk: A Randomized Clinical Trial.

BACKGROUND: Existing data suggest that after concussion, athletes experience an increased risk of subsequent injury. Exploring methods that may reduce injury risk after successful postconcussion return to play may lead to new treatment approaches. PURPOSE: To examine the efficacy of a neuromuscular training (NMT) intervention on acute sports-related time-loss injury over the subsequent year relative to standard of care. STUDY DESIGN: Randomized clinical trial; Level of evidence, 1. METHODS: A total of 27 youth athletes were assessed initially postconcussion (median, 7 days postconcussion; interquartile range [IQR], 5-10) and after return-to-play clearance (median, 40 days postconcussion; IQR, 15-52). After return-to-play clearance, they were randomly assigned to NMT intervention (n = 11; mean ± SD age, 14.7 ± 1.7 years; 36% female) or standard of care (n = 16; mean ± SD age, 15.3 ± 1.8 years; 44% female). The intervention (duration, 8 weeks; frequency, 2 times per week) included guided strength exercises with landing stabilization focus. Standard of care received no recommendations. For the subsequent year, athletes prospectively completed a monthly log of sports-related injuries and organized sports competitions. RESULTS: During the first year after postconcussion return-to-play clearance, sports-related time-loss injuries were more common among standard of care relative to NMT intervention (75% [95% CI, 48%-93%] vs 36% [95% CI, 11%-69%]). After adjusting for age and sex, the hazard of subsequent injury in the standard-of-care group was 3.56 times (95% CI, 1.11-11.49; P = .0334) that of the NMT intervention group. Sports participation was similar between NMT intervention and standard of care during the year-long monitoring period (hours of organized sports per month; median, 12 [IQR, 2.6-32.1] vs 15.6 [IQR, 3.5-105.9]; P = .55). The age- and sex-adjusted incidence of injuries was 10.2 per 1000 competitive exposures (95% CI, 3.7-28.4) in the standard-of-care group as opposed to 3.4 per 1000 (95% CI, 0.9-13.4) in the NMT intervention group. After adjusting for age and sex, incidence of injuries was higher for standard of care vs NMT intervention (rate ratio, 2.96 [95% CI, 0.89-9.85]; P = .076). CONCLUSION: Although preliminary, our findings suggest that an NMT intervention initiated after return-to-play clearance may significantly reduce sports-related time-loss injuries over the subsequent year. REGISTRATION: NCT03917290 (ClinicalTrials.gov identifier).

How Do We Best Educate Our Patients' Caregivers? Comparing the Efficacy of Print Versus Media-based Education Materials in Peripheral Nerve Catheter and Pain Pump Education.

BACKGROUND: Media-based educational materials (EMs) are becoming prominent. The purpose of this study was to compare print versus media-based EMs given to caregivers of pediatric and young adult patients undergoing surgery. We aimed to see whether print or media-based EMs lead to greater caregiver satisfaction, comfort, and preparedness for outpatient peripheral nerve catheter and pain pump management. We also assessed caregiver preference for EM modality. HYPOTHESIS: We hypothesized that media-based EMs would demonstrate greater overall efficacy and thus generate higher caregiver preference. STUDY DESIGN: Randomized control trial. METHODS: After IRB approval, clinicaltrials.gov registration (17-0638), and informed consent, caregivers were randomized to either media or print-based EM groups. Caregivers reviewed their assigned EM and completed a standardized assessment of their comprehension. We assessed caregiver satisfaction, preparedness, and comfort level with the content on a 5-point Likert scale. On postoperative days 1 to 2, caregivers reported satisfaction, comfort, and preference for EM modality. An intent-to-treat analysis was used to compare the 2 groups. RESULTS: From our final cohort of 135 caregivers, we found no difference [P>0.05] in satisfaction, comfort level, level of preparedness, or discharge readiness scores between groups. After the caregivers were given both EMs, they were evenly split in their preference for print (49.6%) versus video (50.4%) based methods. CONCLUSIONS: We did not detect a significant difference in caregiver preference or feelings of preparedness between groups. Interestingly, a significant proportion of caregivers (25%) did not feel comfortable managing the peripheral nerve catheter and its pain pump at home. Future studies should work to improve caregiver comfort with educational content before patient discharge. CLINICAL RELEVANCE: Providers and institutions should feel comfortable providing both print and media-based patient and caregiver education. Caregiver education may be best suited based on caregiver preference of one EM modality versus the other. LEVEL OF EVIDENCE: Level I.

An effective processing pipeline for harmonizing DNA methylation data from Illumina's 450K and EPIC platforms for epidemiological studies.

OBJECTIVE: Illumina BeadChip arrays are commonly used to generate DNA methylation data for large epidemiological studies. Updates in technology over time create challenges for data harmonization within and between studies, many of which obtained data from the older 450K and newer EPIC platforms. The pre-processing pipeline for DNA methylation is not trivial, and influences the downstream analyses. Incorporating different platforms adds a new level of technical variability that has not yet been taken into account by recommended pipelines. Our study evaluated the performance of various tools on different versions of platform data harmonization at each step of pre-processing pipeline, including quality control (QC), normalization, batch effect adjustment, and genomic inflation. We illustrate our novel approach using 450K and EPIC data from the Diabetes Autoimmunity Study in the Young (DAISY) prospective cohort. RESULTS: We found normalization and probe filtering had the biggest effect on data harmonization. Employing a meta-analysis was an effective and easily executable method for accounting for platform variability. Correcting for genomic inflation also helped with harmonization. We present guidelines for studies seeking to harmonize data from the 450K and EPIC platforms, which includes the use of technical replicates for evaluating numerous pre-processing steps, and employing a meta-analysis.

A Mediation Approach to Discovering Causal Relationships between the Metabolome and DNA Methylation in Type 1 Diabetes.

Environmental factors including viruses, diet, and the metabolome have been linked with the appearance of islet autoimmunity (IA) that precedes development of type 1 diabetes (T1D). We measured global DNA methylation (DNAm) and untargeted metabolomics prior to IA and at the time of seroconversion to IA in 92 IA cases and 91 controls from the Diabetes Autoimmunity Study in the Young (DAISY). Causal mediation models were used to identify seven DNAm probe-metabolite pairs in which the metabolite measured at IA mediated the protective effect of the DNAm probe measured prior to IA against IA risk. These pairs included five DNAm probes mediated by histidine (a metabolite known to affect T1D risk), one probe (cg01604946) mediated by phostidyl choline p-32:0 or o-32:1, and one probe (cg00390143) mediated by sphingomyelin d34:2. The top 100 DNAm probes were over-represented in six reactome pathways at the FDR <0.1 level (q = 0.071), including transport of small molecules and inositol phosphate metabolism. While the causal pathways in our mediation models require further investigation to better understand the biological mechanisms, we identified seven methylation sites that may improve our understanding of epigenetic protection against T1D as mediated by the metabolome.

Severity of Idiopathic Scoliosis Is Associated with Differential Methylation: An Epigenome-Wide Association Study of Monozygotic Twins with Idiopathic Scoliosis.

Epigenetic mechanisms may contribute to idiopathic scoliosis (IS). We identified 8 monozygotic twin pairs with IS, 6 discordant (Cobb angle difference > 10°) and 2 concordant (Cobb angle difference ≤ 2°). Genome-wide methylation in blood was measured with the Infinium HumanMethylation EPIC Beadchip. We tested for differences in methylation and methylation variability between discordant twins and tested the association between methylation and curve severity in all twins. Differentially methylated region (DMR) analyses identified gene promoter regions. Methylation at cg12959265 (chr. 7 DPY19L1) was less variable in cases (false discovery rate (FDR) = 0.0791). We identified four probes (false discovery rate, FDR < 0.10); cg02477677 (chr. 17, RARA gene), cg12922161 (chr. 2 LOC150622 gene), cg08826461 (chr. 2), and cg16382077 (chr. 7) associated with curve severity. We identified 57 DMRs where hyper- or hypo-methylation was consistent across the region and 28 DMRs with a consistent association with curve severity. Among DMRs, 21 were correlated with bone methylation. Prioritization of regions based on methylation concordance in bone identified promoter regions for WNT10A (WNT signaling), NPY (regulator of bone and energy homeostasis), and others predicted to be relevant for bone formation/remodeling. These regions may aid in understanding the complex interplay between genetics, environment, and IS.

Single season re-injury risk after concussion and lower extremity injury among male, collision sport, high school athletes.

OBJECTIVES: Research suggests that a recent concussion increases subsequent lower extremity injury risk; however, data in high school athletes is limited. This study evaluates the association between concussion and subsequent injury risk among male, collision sport, high school athletes over a single season. DESIGN: Retrospective cohort study. METHODS: Data were obtained from 2005/06-2014/15 in the High School Reporting Information Online database. A two stage, multi-method matching process was used to identify athletes who suffered multiple injuries over a single athletic season. Demographics and injury characteristics were compared with Chi square and Student's t-tests. Multiple Cox Proportional Hazards regression analysis was used to test whether index injury type was associated with hazard of subsequent injury following return to play from index injury. RESULTS: Overall 1364 HS athletes sustained 2 injuries over a single athletic season (subsequent injury within 45 days of the index injury). Index injuries included 686 (50.2%) lower extremity injuries, 417 (30.6%) upper extremity injuries, and 261 (19.1%) concussions. Hazard of subsequent concussion was increased in the index concussion group relative to the index lower extremity injury group [hazard ratio (HR): 1.60, 95% CI: 1.15-2.23, p = 0.0052]. Hazard of a subsequent lower extremity injury was not significantly different for the index lower extremity injury group relative to the index concussion group [HR: 1.27, 95% CI: 0.98-1.65, p = 0.0728]. CONCLUSIONS: History of recent concussion or recent lower extremity injury are both risk factors for subsequent lower extremity injury in male, collision sport, high school athletes.

Inverse probability weighting is an effective method to address selection bias during the analysis of high dimensional data.

Omics studies frequently use samples collected during cohort studies. Conditioning on sample availability can cause selection bias if sample availability is nonrandom. Inverse probability weighting (IPW) is purported to reduce this bias. We evaluated IPW in an epigenome-wide analysis testing the association between DNA methylation (261,435 probes) and age in healthy adolescent subjects (n = 114). We simulated age and sex to be correlated with sample selection and then evaluated four conditions: complete population/no selection bias (all subjects), naïve selection bias (no adjustment), and IPW selection bias (selection bias with IPW adjustment). Assuming the complete population condition represented the "truth," we compared each condition to the complete population condition. Bias or difference in associations between age and methylation was reduced in the IPW condition versus the naïve condition. However, genomic inflation and type 1 error were higher in the IPW condition relative to the naïve condition. Postadjustment using bacon, type 1 error and inflation were similar across all conditions. Power was higher under the IPW condition compared with the naïve condition before and after inflation adjustment. IPW methods can reduce bias in genome-wide analyses. Genomic inflation is a potential concern that can be minimized using methods that adjust for inflation.