Multi-omics in nasal epithelium reveals three axes of dysregulation for asthma risk in the African Diaspora populations.

Asthma has striking disparities across ancestral groups, but the molecular underpinning of these differences is poorly understood and minimally studied. A goal of the Consortium on Asthma among African-ancestry Populations in the Americas (CAAPA) is to understand multi-omic signatures of asthma focusing on populations of African ancestry. RNASeq and DNA methylation data are generated from nasal epithelium including cases (current asthma, N = 253) and controls (never-asthma, N = 283) from 7 different geographic sites to identify differentially expressed genes (DEGs) and gene networks. We identify 389 DEGs; the top DEG, FN1, was downregulated in cases (q = 3.26 × 10-9) and encodes fibronectin which plays a role in wound healing. The top three gene expression modules implicate networks related to immune response (CEACAM5; p = 9.62 × 10-16 and CPA3; p = 2.39 × 10-14) and wound healing (FN1; p = 7.63 × 10-9). Multi-omic analysis identifies FKBP5, a co-chaperone of glucocorticoid receptor signaling known to be involved in drug response in asthma, where the association between nasal epithelium gene expression is likely regulated by methylation and is associated with increased use of inhaled corticosteroids. This work reveals molecular dysregulation on three axes - increased Th2 inflammation, decreased capacity for wound healing, and impaired drug response - that may play a critical role in asthma within the African Diaspora.

Islet autoantibodies as precision diagnostic tools to characterize heterogeneity in type 1 diabetes: a systematic review.

BACKGROUND: Islet autoantibodies form the foundation for type 1 diabetes (T1D) diagnosis and staging, but heterogeneity exists in T1D development and presentation. We hypothesized that autoantibodies can identify heterogeneity before, at, and after T1D diagnosis, and in response to disease-modifying therapies. METHODS: We systematically reviewed PubMed and EMBASE databases (6/14/2022) assessing 10 years of original research examining relationships between autoantibodies and heterogeneity before, at, after diagnosis, and in response to disease-modifying therapies in individuals at-risk or within 1 year of T1D diagnosis. A critical appraisal checklist tool for cohort studies was modified and used for risk of bias assessment. RESULTS: Here we show that 152 studies that met extraction criteria most commonly characterized heterogeneity before diagnosis (91/152). Autoantibody type/target was most frequently examined, followed by autoantibody number. Recurring themes included correlations of autoantibody number, type, and titers with progression, differing phenotypes based on order of autoantibody seroconversion, and interactions with age and genetics. Only 44% specifically described autoantibody assay standardization program participation. CONCLUSIONS: Current evidence most strongly supports the application of autoantibody features to more precisely define T1D before diagnosis. Our findings support continued use of pre-clinical staging paradigms based on autoantibody number and suggest that additional autoantibody features, particularly in relation to age and genetic risk, could offer more precise stratification. To improve reproducibility and applicability of autoantibody-based precision medicine in T1D, we propose a methods checklist for islet autoantibody-based manuscripts which includes use of precision medicine MeSH terms and participation in autoantibody standardization workshops.

Islet autoantibodies are markers found in the blood when insulin-producing cells in the pancreas become damaged and can be used to predict future development of type 1 diabetes. We evaluated published literature to determine whether characteristics of islet antibodies (type, levels, numbers) could improve prediction and help understand differences in how individuals with type 1 diabetes respond to treatments. We found existing evidence shows that islet autoantibody type and number are most useful to predict disease progression before diagnosis. In addition, the age when islet autoantibodies first appear strongly influences rate of progression. These findings provide important information for patients and care providers on how islet autoantibodies can be used to understand future type 1 diabetes development and to identify individuals who have the potential to benefit from intervention or prevention therapy.

Disease-modifying therapies and features linked to treatment response in type 1 diabetes prevention: a systematic review.

BACKGROUND: Type 1 diabetes (T1D) results from immune-mediated destruction of insulin-producing beta cells. Prevention efforts have focused on immune modulation and supporting beta cell health before or around diagnosis; however, heterogeneity in disease progression and therapy response has limited translation to clinical practice, highlighting the need for precision medicine approaches to T1D disease modification. METHODS: To understand the state of knowledge in this area, we performed a systematic review of randomized-controlled trials with ≥50 participants cataloged in PubMed or Embase from the past 25 years testing T1D disease-modifying therapies and/or identifying features linked to treatment response, analyzing bias using a Cochrane-risk-of-bias instrument. RESULTS: We identify and summarize 75 manuscripts, 15 describing 11 prevention trials for individuals with increased risk for T1D, and 60 describing treatments aimed at preventing beta cell loss at disease onset. Seventeen interventions, mostly immunotherapies, show benefit compared to placebo (only two prior to T1D onset). Fifty-seven studies employ precision analyses to assess features linked to treatment response. Age, beta cell function measures, and immune phenotypes are most frequently tested. However, analyses are typically not prespecified, with inconsistent methods of reporting, and tend to report positive findings. CONCLUSIONS: While the quality of prevention and intervention trials is overall high, the low quality of precision analyses makes it difficult to draw meaningful conclusions that inform clinical practice. To facilitate precision medicine approaches to T1D prevention, considerations for future precision studies include the incorporation of uniform outcome measures, reproducible biomarkers, and prespecified, fully powered precision analyses into future trial design.

Type 1 diabetes (T1D) is a condition that results from the destruction of a type of cell in the pancreas that produces the hormone insulin, leading to lifelong dependence on insulin injections. T1D prevention remains a challenging goal, largely due to the immense variability in disease processes and progression. Therapies tested to date in medical research settings (clinical trials) work only in a subset of individuals, highlighting the need for more tailored prevention approaches. We reviewed clinical trials of therapies targeting the disease process in T1D. While the overall quality of trials was high, studies testing individual features affecting responses to treatments were low. This review reveals an important need to carefully plan high-quality analyses of features that affect treatment response in T1D, to ensure that tailored approaches may one day be applied to clinical practice.

Type 1 Diabetes Prevention: a systematic review of studies testing disease-modifying therapies and features linked to treatment response.

Background: Type 1 diabetes (T1D) results from immune-mediated destruction of insulin-producing beta cells. Efforts to prevent T1D have focused on modulating immune responses and supporting beta cell health; however, heterogeneity in disease progression and responses to therapies have made these efforts difficult to translate to clinical practice, highlighting the need for precision medicine approaches to T1D prevention. Methods: To understand the current state of knowledge regarding precision approaches to T1D prevention, we performed a systematic review of randomized-controlled trials from the past 25 years testing disease-modifying therapies in T1D and/or identifying features linked to treatment response, analyzing bias using a Cochrane-risk-of-bias instrument. Results: We identified 75 manuscripts, 15 describing 11 prevention trials for individuals with increased risk for T1D, and 60 describing treatments aimed at preventing beta cell loss in individuals at disease onset. Seventeen agents tested, mostly immunotherapies, showed benefit compared to placebo (only two prior to T1D onset). Fifty-seven studies employed precision analyses to assess features linked to treatment response. Age, measures of beta cell function and immune phenotypes were most frequently tested. However, analyses were typically not prespecified, with inconsistent methods reporting, and tended to report positive findings. Conclusions: While the quality of prevention and intervention trials was overall high, low quality of precision analyses made it difficult to draw meaningful conclusions that inform clinical practice. Thus, prespecified precision analyses should be incorporated into the design of future studies and reported in full to facilitate precision medicine approaches to T1D prevention. Plain Language Summary: Type 1 diabetes (T1D) results from the destruction of insulin-producing cells in the pancreas, necessitating lifelong insulin dependence. T1D prevention remains an elusive goal, largely due to immense variability in disease progression. Agents tested to date in clinical trials work in a subset of individuals, highlighting the need for precision medicine approaches to prevention. We systematically reviewed clinical trials of disease-modifying therapy in T1D. While age, measures of beta cell function, and immune phenotypes were most commonly identified as factors that influenced treatment response, the overall quality of these studies was low. This review reveals an important need to proactively design clinical trials with well-defined analyses to ensure that results can be interpreted and applied to clinical practice.

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.

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.

COVID-19 Mortality in the Colorado Center for Personalized Medicine Biobank.

Over 6.37 million people have died from COVID-19 worldwide, but factors influencing COVID-19-related mortality remain understudied. We aimed to describe and identify risk factors for COVID-19 mortality in the Colorado Center for Personalized Medicine (CCPM) Biobank using integrated data sources, including Electronic Health Records (EHRs). We calculated cause-specific mortality and case-fatality rates for COVID-19 and common pre-existing health conditions defined by diagnostic phecodes and encounters in EHRs. We performed multivariable logistic regression analyses of the association between each pre-existing condition and COVID-19 mortality. Of the 155,859 Biobank participants enrolled as of July 2022, 20,797 had been diagnosed with COVID-19. Of 5334 Biobank participants who had died, 190 were attributed to COVID-19. The case-fatality rate was 0.91% and the COVID-19 mortality rate was 122 per 100,000 persons. The odds of dying from COVID-19 were significantly increased among older men, and those with 14 of the 61 pre-existing conditions tested, including hypertensive chronic kidney disease (OR: 10.14, 95% CI: 5.48, 19.16) and type 2 diabetes with renal manifestations (OR: 5.59, 95% CI: 3.42, 8.97). Male patients who are older and have pre-existing kidney diseases may be at higher risk for death from COVID-19 and may require special care.

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.

COVID-19 Surveillance in the Biobank at the Colorado Center for Personalized Medicine: Observational Study.

BACKGROUND: Characterizing the experience and impact of the COVID-19 pandemic among various populations remains challenging due to the limitations inherent in common data sources, such as electronic health records (EHRs) or cross-sectional surveys. OBJECTIVE: This study aims to describe testing behaviors, symptoms, impact, vaccination status, and case ascertainment during the COVID-19 pandemic using integrated data sources. METHODS: In summer 2020 and 2021, we surveyed participants enrolled in the Biobank at the Colorado Center for Personalized Medicine (CCPM; N=180,599) about their experience with COVID-19. The prevalence of testing, symptoms, and impacts of COVID-19 on employment, family life, and physical and mental health were calculated overall and by demographic categories. Survey respondents who reported receiving a positive COVID-19 test result were considered a "confirmed case" of COVID-19. Using EHRs, we compared COVID-19 case ascertainment and characteristics in EHRs versus the survey. Positive cases were identified in EHRs using the International Statistical Classification of Diseases, 10th revision (ICD-10) diagnosis codes, health care encounter types, and encounter primary diagnoses. RESULTS: Of the 25,063 (13.9%) survey respondents, 10,661 (42.5%) had been tested for COVID-19, and of those, 1366 (12.8%) tested positive. Nearly half of those tested had symptoms or had been exposed to someone who was infected. Young adults (18-29 years) and Hispanics were more likely to have positive tests compared to older adults and persons of other racial/ethnic groups. Mental health (n=13,688, 54.6%) and family life (n=12,233, 48.8%) were most negatively affected by the pandemic and more so among younger groups and women; negative impacts on employment were more commonly reported among Black respondents. Of the 10,249 individuals who responded to vaccination questions from version 2 of the survey (summer 2021), 9770 (95.3%) had received the vaccine. After integration with EHR data up to the time of the survey completion, 1006 (4%) of the survey respondents had a discordant COVID-19 case status between EHRs and the survey. Using all longitudinal EHR and survey data, we identified 11,472 (6.4%) COVID-19-positive cases among Biobank participants. In comparison to COVID-19 cases identified through the survey, EHR-identified cases were younger and more likely to be Hispanic. CONCLUSIONS: We found that the COVID-19 pandemic has had far-reaching and varying effects among our Biobank participants. Integrated data assets, such as the Biobank at the CCPM, are key resources for population health monitoring in response to public health emergencies, such as the COVID-19 pandemic.

Lipid mediators are detectable in the nasal epithelium and differ by asthma status in female subjects.

BACKGROUND: Lipid mediators, bioactive products of polyunsaturated fatty acid metabolism, contribute to inflammation initiation and resolution in allergic diseases; however, their presence in lung-related biosamples has not been fully described. OBJECTIVE: We aimed to quantify lipid mediators in the nasal airway epithelium and characterize preliminary associations with asthma. METHODS: Using liquid chromatography-mass spectrometry, we conducted a pilot study to quantify 56 lipid mediators from nasal epithelial samples collected from 11 female participants of an outpatient asthma clinic and community controls (aged 30-55 years). We examined the presence of each compound using descriptive statistics to test whether lipid mediators could distinguish subjects with asthma (n = 8) from control subjects (n = 3) using linear regression and partial least squares discriminant analysis. RESULTS: Fifteen lipid mediators were detectable in all samples, including resolvin (Rv) D5 (RvD5), with the highest median concentrations (in pg/µg protein) of 13-HODE (126.481), 15-HETE (32.869), and 13-OxoODE (13.251). From linear regression adjusted for age, prostaglandin E2 (PGE2) had a trend (P < .1) for higher concentrations in patients with severe asthma compared to controls (mean difference, 0.95; 95% confidence interval, -0.04 to 1.95). Asthma patients had higher scores on principal component 3 compared to controls (mean difference, 2.42; 95% confidence interval, 0.89 to 3.96), which represented lower levels of proresolving 15-HEPE, 19,20-DiHDPA, RvD5, 14-HDHA, 17-HDHA, and 13-HOTrE. Most of these compounds were best at discriminating asthma cases from controls in partial least squares discriminant analysis. CONCLUSION: Lipid mediators are detectable in the nasal epithelium, and their levels distinguish asthma cases from controls.

Discovering metabolite quantitative trait loci in asthma using an isolated population.

BACKGROUND: Integration of metabolomics with genetics may advance understanding of disease pathogenesis but has been underused in asthma genetic studies. OBJECTIVE: We sought to discover new genetic effects in asthma and to characterize the molecular consequences of asthma genetic risk through integration with the metabolome in a homogeneous population. METHODS: From fasting serum samples collected on 348 Tangier Island residents, we quantified 2612 compounds using untargeted metabolomics. Genotyping was performed using Illumina's MEGA array imputed to the TOPMed reference panel. To prioritize metabolites for genome-wide association analysis, we performed a metabolome-wide association study with asthma, selecting asthma-associated metabolites with heritability q value less than 0.01 for genome-wide association analysis. We also tested the association between all metabolites and 8451 candidate asthma single nucleotide polymorphisms previously associated with asthma in the UK Biobank. We followed up significant associations by characterizing shared genetic signal for metabolites and asthma using colocalization analysis. For detailed Methods, please see this article's Online Repository at www.jacionline.org. RESULTS: A total of 60 metabolites were associated with asthma (P < .01), including 40 heritable metabolites tested in genome-wide association analysis. We observed a strong association peak for the endocannabinoid linoleoyl ethanolamide on chromosome 6 in VNN1 (P < 2.7 × 10-9). We found strong evidence (colocalization posterior probability >75%) for a shared causal variant between 3 metabolites and asthma, including the polyamine acisoga and variants in LPP, and derivative leukotriene B4 and intergenic variants in chr10p14. CONCLUSIONS: We identified novel metabolite quantitative trait loci with asthma associations. Identification and characterization of these genetically driven metabolites may provide insight into the functional consequences of genetic risk factors for asthma.

Epigenome-Wide Association Study of Infant Feeding and DNA Methylation in Infancy and Childhood in a Population at Increased Risk for Type 1 Diabetes.

We assessed associations between infant diet (e.g., breastfeeding and introduction to solid foods) and DNA methylation in infancy and childhood. We measured DNA methylation in peripheral blood collected in infancy (9-15 months of age) in 243 children; and in a subset of 50 children, we also measured methylation in childhood (6-9 years of age) to examine persistence, and at birth (in cord blood) to examine temporality. We performed multivariable linear regression of infant diet on the outcome of methylation using epigenome-wide and candidate site approaches. We identified six novel CpG sites associated with breastfeeding duration using an EWAS approach. One differentially methylated site presented directionally consistent associations with breastfeeding (cg00574958, CPT1A) in infancy and childhood but not at birth. Two differentially methylated sites in infancy (cg19693031, TXNIP; cg23307264, KHSRP) were associated with breastfeeding and were not present at birth; however, these associations did not persist into childhood. Associations between infant diet and methylation in infancy at three sites (cg22369607, AP001525.1; cg2409200, TBCD; cg27173510, PGBD5) were also present at birth, suggesting the influence of exposures other than infant diet. Infant diet exposures are associated with persistent methylation differences in CPT1A, which may be one mechanism behind infant diet's long-term health effects.

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.

Phospholipid Levels at Seroconversion Are Associated With Resolution of Persistent Islet Autoimmunity: The Diabetes Autoimmunity Study in the Young.

Reversion of islet autoimmunity (IA) may point to mechanisms that prevent IA progression. We followed 199 individuals who developed IA during the Diabetes Autoimmunity Study in the Young. Untargeted metabolomics was performed in serum samples following IA. Cox proportional hazards models were used to test whether the metabolites (2,487) predicted IA reversion: two or more consecutive visits negative for all autoantibodies. We conducted a principal components analysis (PCA) of the top metabolites; |hazard ratio (HR) >1.25| and nominal P < 0.01. Phosphatidylcholine (16:0_18:1(9Z)) was the strongest individual metabolite (HR per 1 SD 2.16, false discovery rate (FDR)-adjusted P = 0.0037). Enrichment analysis identified four clusters (FDR P < 0.10) characterized by an overabundance of sphingomyelin (d40:0), phosphatidylcholine (16:0_18:1(9Z)), phosphatidylcholine (30:0), and l-decanoylcarnitine. Overall, 63 metabolites met the criteria for inclusion in the PCA. PC1 (HR 1.4, P < 0.0001), PC2 (HR 0.85, P = 0.0185), and PC4 (HR 1.28, P = 0.0103) were associated with IA reversion. Given the potential influence of diet on the metabolome, we investigated whether nutrients were correlated with PCs. We identified 20 nutrients that were correlated with the PCs (P < 0.05). Total sugar intake was the top nutrient. Overall, we identified an association between phosphatidylcholine, sphingomyelin, and carnitine levels and reversion of IA.

Maternal food consumption during late pregnancy and offspring risk of islet autoimmunity and type 1 diabetes.

AIMS/HYPOTHESIS: We aimed to investigate the association between maternal consumption of gluten-containing foods and other selected foods during late pregnancy and offspring risk of islet autoimmunity (IA) and type 1 diabetes in The Environmental Determinants of Diabetes in the Young (TEDDY) study. METHODS: The TEDDY study recruited children at high genetic risk for type 1 diabetes at birth, and prospectively follows them for the development of IA and type 1 diabetes (n = 8556). A questionnaire on the mother's diet in late pregnancy was completed by 3-4 months postpartum. The maternal daily intake was estimated from a food frequency questionnaire for eight food groups: gluten-containing foods, non-gluten cereals, fresh milk, sour milk, cheese products, soy products, lean/medium-fat fish and fatty fish. For each food, we described the distribution of maternal intake among the four participating countries in the TEDDY study and tested the association of tertile of maternal food consumption with risk of IA and type 1 diabetes using forward selection time-to-event Cox regression. RESULTS: By 28 February 2019, 791 cases of IA and 328 cases of type 1 diabetes developed in TEDDY. There was no association between maternal late-pregnancy consumption of gluten-containing foods or any of the other selected foods and risk of IA, type 1 diabetes, insulin autoantibody-first IA or GAD autoantibody-first IA (all p ≥ 0.01). Maternal gluten-containing food consumption in late pregnancy was higher in Sweden (242 g/day), Germany (247 g/day) and Finland (221 g/day) than in the USA (199 g/day) (pairwise p < 0.05). CONCLUSIONS/INTERPRETATION: Maternal food consumption during late pregnancy was not associated with offspring risk for IA or type 1 diabetes. TRIAL REGISTRATION: ClinicalTrials.gov NCT00279318.

Predictors of oxylipins in a healthy pediatric population.

BACKGROUND: Oxylipins are formed from oxidation of omega-6 (n6) and omega-3 (n3) fatty acids (FAs). Evidence for inflammatory effects comes mostly from adults. METHODS: Oxylipins from n6 FA (27 n6-oxylipins) and n3 FA (12 n3-oxylipins) were measured through ultra-high-performance liquid chromatography-mass spectrometry (LC-MS/MS) in plasma from 111 children at risk of type 1 diabetes (age 1-17 years) studied longitudinally. Oxylipin precursor FAs (arachidonic acid, linoleic acid, alpha-linolenic acid, docosahexaenoic acid, eicosapentaenoic acid) were measured in red blood cell (RBC) membrane and plasma. Precursor FAs dietary intake was measured through food frequency questionnaire and environmental tobacco smoke (ETS) through questionnaires. Linear mixed models were used to test oxylipins with predictors. RESULTS: Age associated with 15 n6- and 6 n3-oxylipins; race/ethnicity associated with 3 n6- and 1 n3-oxylipins; sex associated with 2 n6-oxylipins. ETS associated with lipoxin-A4. Oxylipins associated with precursor FAs in plasma more often than RBC. RBC levels and dietary intake of precursor FAs more consistently associated with n3-oxylipins than with n6-oxylipins. CONCLUSIONS: In healthy children, oxylipin levels change with age. Oxylipins associated with precursor FAs more often in plasma than RBC or diet, suggesting that inflammatory regulation leading to FA release into plasma may also be a determinant of oxylipin generation. IMPACT: This is the first study to examine predictors of oxylipins in healthy children at risk of type 1 diabetes. In healthy children at risk of type 1 diabetes, many oxylipins change with age, and most oxylipins do not differ by sex or race/ethnicity. Environmental tobacco smoke exposure was associated with the presence of lipoxin A4. Omega-6- and omega-3-related oxylipin levels were consistently associated with their respective precursor fatty acid levels measured in the plasma. Proportionally more omega-3 compared to omega-6 oxylipins were associated with dietary intake and red blood cell membrane levels of the respective precursor fatty acid.

Prediction of the development of islet autoantibodies through integration of environmental, genetic, and metabolic markers.

BACKGROUND: The Environmental Determinants of the Diabetes in the Young (TEDDY) study has prospectively followed, from birth, children at increased genetic risk of type 1 diabetes. TEDDY has collected heterogenous data longitudinally to gain insights into the environmental and biological mechanisms driving the progression to persistent islet autoantibodies. METHODS: We developed a machine learning model to predict imminent transition to the development of persistent islet autoantibodies based on time-varying metabolomics data integrated with time-invariant risk factors (eg, gestational age). The machine learning was initiated with 221 potential features (85 genetic, 5 environmental, 131 metabolomic) and an ensemble-based feature evaluation was utilized to identify a small set of predictive features that can be interrogated to better understand the pathogenesis leading up to persistent islet autoimmunity. RESULTS: The final integrative machine learning model included 42 disparate features, returning a cross-validated receiver operating characteristic area under the curve (AUC) of 0.74 and an AUC of ~0.65 on an independent validation dataset. The model identified a principal set of 20 time-invariant markers, including 18 genetic markers (16 single nucleotide polymorphisms [SNPs] and two HLA-DR genotypes) and two demographic markers (gestational age and exposure to a prebiotic formula). Integration with the metabolome identified 22 supplemental metabolites and lipids, including adipic acid and ceramide d42:0, that predicted development of islet autoantibodies. CONCLUSIONS: The majority (86%) of metabolites that predicted development of islet autoantibodies belonged to three pathways: lipid oxidation, phospholipase A2 signaling, and pentose phosphate, suggesting that these metabolic processes may play a role in triggering islet autoimmunity.