Metabolomics of IgE-Mediated Food Allergy and Oral Immunotherapy Outcomes based on Metabolomic Profiling.

Background: The immunometabolic mechanisms underlying variable responses to oral immunotherapy (OIT) in patients with IgE-mediated food allergy are unknown. Objective: To identify novel pathways associated with tolerance in food allergy, we used metabolomic profiling to find pathways important for food allergy in multi-ethnic cohorts and responses to OIT. Methods: Untargeted plasma metabolomics data were generated from the VDAART healthy infant cohort (N=384), a Costa Rican cohort of children with asthma (N=1040), and a peanut OIT trial (N=20) evaluating sustained unresponsiveness (SU, protection that lasts after therapy) versus transient desensitization (TD, protection that ends immediately afterwards). Generalized linear regression modeling and pathway enrichment analysis identified metabolites associated with food allergy and OIT outcomes. Results: Compared with unaffected children, those with food allergy were more likely to have metabolomic profiles with altered histidines and increased bile acids. Eicosanoids (e.g., arachidonic acid derivatives) (q=2.4×10 -20 ) and linoleic acid derivatives (q=3.8×10 -5 ) pathways decreased over time on OIT. Comparing SU versus TD revealed differing concentrations of bile acids (q=4.1×10 -8 ), eicosanoids (q=7.9×10 -7 ), and histidine pathways (q=0.015). In particular, the bile acid lithocholate (4.97[1.93,16.14], p=0.0027), the eicosanoid leukotriene B4 (3.21[1.38,8.38], p=0.01), and the histidine metabolite urocanic acid (22.13[3.98,194.67], p=0.0015) were higher in SU. Conclusions: We observed distinct profiles of bile acids, histidines, and eicosanoids that vary among patients with food allergy, over time on OIT and between SU and TD. Participants with SU had higher levels of metabolites such as lithocholate and urocanic acid, which have immunomodulatory roles in key T-cell subsets, suggesting potential mechanisms of tolerance in immunotherapy. Key Messages: - Compared with unaffected controls, children with food allergy demonstrated higher levels of bile acids and distinct histidine/urocanic acid profiles, suggesting a potential role of these metabolites in food allergy. - In participants receiving oral immunotherapy for food allergy, those who were able to maintain tolerance-even after stopping therapyhad lower overall levels of bile acid and histidine metabolites, with the exception of lithocholic acid and urocanic acid, two metabolites that have roles in T cell differentiation that may increase the likelihood of remission in immunotherapy. Capsule summary: This is the first study of plasma metabolomic profiles of responses to OIT in individuals with IgE-mediated food allergy. Identification of immunomodulatory metabolites in allergic tolerance may help identify mechanisms of tolerance and guide future therapeutic development.

Network Analysis Reveals Protein Modules Associated with Childhood Respiratory Diseases.

Background: The first year of life is a period of rapid immune development that can impact health trajectories and the risk of developing respiratory-related diseases, such as asthma, recurrent infections, and eczema. However, the biology underlying subsequent disease development remains unknown. Methods: Using weighted gene correlation network analysis (WGCNA), we derived modules of highly correlated immune-related proteins in plasma samples from children at age 1 year (N=294) from the Vitamin D Antenatal Asthma Reduction Trial (VDAART). We applied regression analyses to assess relationships between protein modules and development of childhood respiratory diseases up to age 6 years. We then characterized genomic, environmental, and metabolomic factors associated with modules. Results: WGCNA identified four protein modules at age 1 year associated with incidence of childhood asthma and/or recurrent wheeze (Padj range: 0.02-0.03), respiratory infections (Padj range: 6.3×10-9-2.9×10-6), and eczema (Padj=0.01) by age 6 years; three modules were associated with at least one environmental exposure (Padj range: 2.8×10-10-0.03) and disrupted metabolomic pathway(s) (Padj range: 2.8×10-6-0.04). No genome-wide SNPs were identified as significant genetic risk factors for any protein module. Relationships between protein modules with clinical, environmental, and 'omic factors were temporally sensitive and could not be recapitulated in protein profiles at age 6 years. Conclusion: These findings suggested protein profiles as early as age 1 year predicted development of respiratory-related diseases through age 6 and were associated with changes in pathways related to amino acid and energy metabolism. These may inform new strategies to identify vulnerable individuals based on immune protein profiling.

A genome-wide association study of mass spectrometry proteomics using the Seer Proteograph platform.

Genome-wide association studies (GWAS) with proteomics are essential tools for drug discovery. To date, most studies have used affinity proteomics platforms, which have limited discovery to protein panels covered by the available affinity binders. Furthermore, it is not clear to which extent protein epitope changing variants interfere with the detection of protein quantitative trait loci (pQTLs). Mass spectrometry-based (MS) proteomics can overcome some of these limitations. Here we report a GWAS using the MS-based Seer Proteograph™ platform with blood samples from a discovery cohort of 1,260 American participants and a replication in 325 individuals from Asia, with diverse ethnic backgrounds. We analysed 1,980 proteins quantified in at least 80% of the samples, out of 5,753 proteins quantified across the discovery cohort. We identified 252 and replicated 90 pQTLs, where 30 of the replicated pQTLs have not been reported before. We further investigated 200 of the strongest associated cis-pQTLs previously identified using the SOMAscan and the Olink platforms and found that up to one third of the affinity proteomics pQTLs may be affected by epitope effects, while another third were confirmed by MS proteomics to be consistent with the hypothesis that genetic variants induce changes in protein expression. The present study demonstrates the complementarity of the different proteomics approaches and reports pQTLs not accessible to affinity proteomics, suggesting that many more pQTLs remain to be discovered using MS-based platforms.

Metabolomic-derived endotypes of age-related macular degeneration (AMD): a step towards identification of disease subgroups.

Age-related macular degeneration (AMD) is a leading cause of blindness worldwide, with a complex pathophysiology and phenotypic diversity. Here, we apply Similarity Network Fusion (SNF) to cluster AMD patients into putative metabolomics-derived endotypes. Using a discovery cohort of 163 AMD patients from Boston, US, and a validation cohort of 214 patients from Coimbra, Portugal, we identified four distinct metabolomics-derived endotypes with varying retinal structural and functional characteristics, confirmed across both cohorts. Patients clustered into Endotype 1 exhibited a milder form of AMD and were characterized by low levels of amino acids in specific metabolic pathways. Meanwhile, patients clustered into both Endotype 3 and 4 were associated with more severe AMD and exhibited low levels of fatty acid metabolites and elevated levels of sphingomyelins and fatty acid metabolites, respectively. These preliminary findings indicate that metabolomics-derived endotyping may offer a refined strategy for categorizing AMD patients based on their specific pathophysiological underpinnings, rather than relying solely on traditional observational clinical indicators.

Quantifying the stochastic component of epigenetic aging.

DNA methylation clocks can accurately estimate chronological age and, to some extent, also biological age, yet the process by which age-associated DNA methylation (DNAm) changes are acquired appears to be quasi-stochastic, raising a fundamental question: how much of an epigenetic clock's predictive accuracy could be explained by a stochastic process of DNAm change? Here, using DNAm data from sorted immune cells, we build realistic simulation models, subsequently demonstrating in over 22,770 sorted and whole-blood samples from 25 independent cohorts that approximately 66-75% of the accuracy underpinning Horvath's clock could be driven by a stochastic process. This fraction increases to 90% for the more accurate Zhang's clock, but is lower (63%) for the PhenoAge clock, suggesting that biological aging is reflected by nonstochastic processes. Confirming this, we demonstrate that Horvath's age acceleration in males and PhenoAge's age acceleration in severe coronavirus disease 2019 cases and smokers are not driven by an increased rate of stochastic change but by nonstochastic processes. These results significantly deepen our understanding and interpretation of epigenetic clocks.

Metabolite signatures associated with microRNA miR-143-3p serve as drivers of poor lung function trajectories in childhood asthma.

BACKGROUND: Lung function trajectories (LFTs) have been shown to be an important measure of long-term health in asthma. While there is a growing body of metabolomic studies on asthma status and other phenotypes, there are no prospective studies of the relationship between metabolomics and LFTs or their genomic determinants. METHODS: We utilized ordinal logistic regression to identify plasma metabolite principal components associated with four previously-published LFTs in children from the Childhood Asthma Management Program (CAMP) (n = 660). The top significant metabolite principal component (PCLF) was evaluated in an independent cross-sectional child cohort, the Genetic Epidemiology of Asthma in Costa Rica Study (GACRS) (n = 1151) and evaluated for association with spirometric measures. Using meta-analysis of CAMP and GACRS, we identified associations between PCLF and microRNA, and SNPs in their target genes. Statistical significance was determined using an false discovery rate-adjusted Q-value. FINDINGS: The top metabolite principal component, PCLF, was significantly associated with better LFTs after multiple-testing correction (Q-value = 0.03). PCLF is composed of the urea cycle, caffeine, corticosteroid, carnitine, and potential microbial (secondary bile acid, tryptophan, linoleate, histidine metabolism) metabolites. Higher levels of PCLF were also associated with increases in lung function measures and decreased circulating neutrophil percentage in both CAMP and GACRS. PCLF was also significantly associated with microRNA miR-143-3p, and SNPs in three miR-143-3p target genes; CCZ1 (P-value = 2.6 × 10-5), SLC8A1 (P-value = 3.9 × 10-5); and TENM4 (P-value = 4.9 × 10-5). INTERPRETATION: This study reveals associations between metabolites, miR-143-3p and LFTs in children with asthma, offering insights into asthma physiology and possible interventions to enhance lung function and long-term health. FUNDING: Molecular data for CAMP and GACRS via the Trans-Omics in Precision Medicine (TOPMed) program was supported by the National Heart, Lung, and Blood Institute (NHLBI).

Metabolomic data presents challenges for epidemiological meta-analysis: a case study of childhood body mass index from the ECHO consortium.

INTRODUCTION: Meta-analyses across diverse independent studies provide improved confidence in results. However, within the context of metabolomic epidemiology, meta-analysis investigations are complicated by differences in study design, data acquisition, and other factors that may impact reproducibility. OBJECTIVE: The objective of this study was to identify maternal blood metabolites during pregnancy (> 24 gestational weeks) related to offspring body mass index (BMI) at age two years through a meta-analysis framework. METHODS: We used adjusted linear regression summary statistics from three cohorts (total N = 1012 mother-child pairs) participating in the NIH Environmental influences on Child Health Outcomes (ECHO) Program. We applied a random-effects meta-analysis framework to regression results and adjusted by false discovery rate (FDR) using the Benjamini-Hochberg procedure. RESULTS: Only 20 metabolites were detected in all three cohorts, with an additional 127 metabolites detected in two of three cohorts. Of these 147, 6 maternal metabolites were nominally associated (P < 0.05) with offspring BMI z-scores at age 2 years in a meta-analytic framework including at least two studies: arabinose (Coefmeta = 0.40 [95% CI 0.10,0.70], Pmeta = 9.7 × 10-3), guanidinoacetate (Coefmeta = - 0.28 [- 0.54, - 0.02], Pmeta = 0.033), 3-ureidopropionate (Coefmeta = 0.22 [0.017,0.41], Pmeta = 0.033), 1-methylhistidine (Coefmeta = - 0.18 [- 0.33, - 0.04], Pmeta = 0.011), serine (Coefmeta = - 0.18 [- 0.36, - 0.01], Pmeta = 0.034), and lysine (Coefmeta = - 0.16 [- 0.32, - 0.01], Pmeta = 0.044). No associations were robust to multiple testing correction. CONCLUSIONS: Despite including three cohorts with large sample sizes (N > 100), we failed to identify significant metabolite associations after FDR correction. Our investigation demonstrates difficulties in applying epidemiological meta-analysis to clinical metabolomics, emphasizes challenges to reproducibility, and highlights the need for standardized best practices in metabolomic epidemiology.

Sphingolipid classes and the interrelationship with pediatric asthma and asthma risk factors.

BACKGROUND: While dysregulated sphingolipid metabolism has been associated with risk of childhood asthma, the specific sphingolipid classes and/or mechanisms driving this relationship remain unclear. We aimed to understand the multifaceted role between sphingolipids and other established asthma risk factors that complicate this relationship. METHODS: We performed targeted LC-MS/MS-based quantification of 77 sphingolipids in plasma from 997 children aged 6 years from two independent cohorts (VDAART and COPSAC2010 ). We examined associations of circulatory sphingolipids with childhood asthma, lung function, and three asthma risk factors: functional SNPs in ORMDL3, low vitamin D levels, and reduced gut microbial maturity. Given racial differences between these cohorts, association analyses were performed separately and then meta-analyzed together. RESULTS: We observed elevations in circulatory sphingolipids with asthma phenotypes and risk factors; however, there were differential associations of sphingolipid classes with clinical outcomes and/or risk factors. While elevations from metabolites involved in ceramide recycling and catabolic pathways were associated with asthma and worse lung function [meta p-value range: 1.863E-04 to 2.24E-3], increased ceramide levels were associated with asthma risk factors [meta p-value range: 7.75E-5 to .013], but not asthma. Further investigation identified that some ceramides acted as mediators while some interacted with risk factors in the associations with asthma outcomes. CONCLUSION: This study demonstrates the differential role that sphingolipid subclasses may play in asthma and its risk factors. While overall elevations in sphingolipids appeared to be deleterious overall; elevations in ceramides were uniquely associated with increases in asthma risk factors only; while elevations in asthma phenotypes were associated with recycling sphingolipids. Modification of asthma risk factors may play an important role in regulating sphingolipid homeostasis via ceramides to affect asthma. Further function work may validate the observed associations.

OMICmAge: An integrative multi-omics approach to quantify biological age with electronic medical records.

Biological aging is a multifactorial process involving complex interactions of cellular and biochemical processes that is reflected in omic profiles. Using common clinical laboratory measures in ~30,000 individuals from the MGB-Biobank, we developed a robust, predictive biological aging phenotype, EMRAge, that balances clinical biomarkers with overall mortality risk and can be broadly recapitulated across EMRs. We then applied elastic-net regression to model EMRAge with DNA-methylation (DNAm) and multiple omics, generating DNAmEMRAge and OMICmAge, respectively. Both biomarkers demonstrated strong associations with chronic diseases and mortality that outperform current biomarkers across our discovery (MGB-ABC, n=3,451) and validation (TruDiagnostic, n=12,666) cohorts. Through the use of epigenetic biomarker proxies, OMICmAge has the unique advantage of expanding the predictive search space to include epigenomic, proteomic, metabolomic, and clinical data while distilling this in a measure with DNAm alone, providing opportunities to identify clinically-relevant interconnections central to the aging process.

Development and psychometric validation of the Pandemic-Related Traumatic Stress Scale for children and adults.

To assess the public health impact of the COVID-19 pandemic on mental health, investigators from the National Institutes of Health Environmental influences on Child Health Outcomes (ECHO) research program developed the Pandemic-Related Traumatic Stress Scale (PTSS). Based on the Diagnostic and Statistical Manual of Mental Disorders, 5th Edition (DSM-5) acute stress disorder symptom criteria, the PTSS is designed for adolescent (13-21 years) and adult self-report and caregiver-report on 3-12-year-olds. To evaluate psychometric properties, we used PTSS data collected between April 2020 and August 2021 from non-pregnant adult caregivers (n = 11,483), pregnant/postpartum individuals (n = 1,656), adolescents (n = 1,795), and caregivers reporting on 3-12-year-olds (n = 2,896). We used Mokken scale analysis to examine unidimensionality and reliability, Pearson correlations to evaluate relationships with other relevant variables, and analyses of variance to identify regional, age, and sex differences. Mokken analysis resulted in a moderately strong, unidimensional scale that retained nine of the original 10 items. We detected small to moderate positive associations with depression, anxiety, and general stress, and negative associations with life satisfaction. Adult caregivers had the highest PTSS scores, followed by adolescents, pregnant/postpartum individuals, and children. Caregivers of younger children, females, and older youth had higher PTSS scores compared to caregivers of older children, males, and younger youth, respectively. (PsycInfo Database Record (c) 2023 APA, all rights reserved).

Polygenic risk scores identify heterogeneity in asthma and chronic obstructive pulmonary disease.

BACKGROUND: Asthma and chronic obstructive pulmonary disease (COPD) have distinct and overlapping genetic and clinical features. OBJECTIVE: We sought to test the hypothesis that polygenic risk scores (PRSs) for asthma (PRSAsthma) and spirometry (FEV1 and FEV1/forced vital capacity; PRSspiro) would demonstrate differential associations with asthma, COPD, and asthma-COPD overlap (ACO). METHODS: We developed and tested 2 asthma PRSs and applied the higher performing PRSAsthma and a previously published PRSspiro to research (Genetic Epidemiology of COPD study and Childhood Asthma Management Program, with spirometry) and electronic health record-based (Mass General Brigham Biobank and Genetic Epidemiology Research on Adult Health and Aging [GERA]) studies. We assessed the association of PRSs with COPD and asthma using modified random-effects and binary-effects meta-analyses, and ACO and asthma exacerbations in specific cohorts. Models were adjusted for confounders and genetic ancestry. RESULTS: In meta-analyses of 102,477 participants, the PRSAsthma (odds ratio [OR] per SD, 1.16 [95% CI, 1.14-1.19]) and PRSspiro (OR per SD, 1.19 [95% CI, 1.17-1.22]) both predicted asthma, whereas the PRSspiro predicted COPD (OR per SD, 1.25 [95% CI, 1.21-1.30]). However, results differed by cohort. The PRSspiro was not associated with COPD in GERA and Mass General Brigham Biobank. In the Genetic Epidemiology of COPD study, the PRSAsthma (OR per SD: Whites, 1.3; African Americans, 1.2) and PRSspiro (OR per SD: Whites, 2.2; African Americans, 1.6) were both associated with ACO. In GERA, the PRSAsthma was associated with asthma exacerbations (OR, 1.18) in Whites; the PRSspiro was associated with asthma exacerbations in White, LatinX, and East Asian participants. CONCLUSIONS: PRSs for asthma and spirometry are both associated with ACO and asthma exacerbations. Genetic prediction performance differs in research versus electronic health record-based cohorts.

Metabolomics and Self-Reported Depression, Anxiety, and Phobic Symptoms in the VA Normative Aging Study.

Traditional approaches to understanding metabolomics in mental illness have focused on investigating a single disorder or comparisons between diagnoses, but a growing body of evidence suggests substantial mechanistic overlap in mental disorders that could be reflected by the metabolome. In this study, we investigated associations between global plasma metabolites and abnormal scores on the depression, anxiety, and phobic anxiety subscales of the Brief Symptom Inventory (BSI) among 405 older males who participated in the Normative Aging Study (NAS). Our analysis revealed overlapping and distinct metabolites associated with each mental health dimension subscale and four metabolites belonging to xenobiotic, carbohydrate, and amino acid classes that were consistently associated across all three symptom dimension subscales. Furthermore, three of these four metabolites demonstrated a higher degree of alteration in men who reported poor scores in all three dimensions compared to men with poor scores in only one, suggesting the potential for shared underlying biology but a differing degree of perturbation when depression and anxiety symptoms co-occur. Our findings implicate pathways of interest relevant to the overlap of mental health conditions in aging veterans and could represent clinically translatable targets underlying poor mental health in this high-risk population.

A meta-analysis of immune-cell fractions at high resolution reveals novel associations with common phenotypes and health outcomes.

BACKGROUND: Changes in cell-type composition of tissues are associated with a wide range of diseases and environmental risk factors and may be causally implicated in disease development and progression. However, these shifts in cell-type fractions are often of a low magnitude, or involve similar cell subtypes, making their reliable identification challenging. DNA methylation profiling in a tissue like blood is a promising approach to discover shifts in cell-type abundance, yet studies have only been performed at a relatively low cellular resolution and in isolation, limiting their power to detect shifts in tissue composition. METHODS: Here we derive a DNA methylation reference matrix for 12 immune-cell types in human blood and extensively validate it with flow-cytometric count data and in whole-genome bisulfite sequencing data of sorted cells. Using this reference matrix, we perform a directional Stouffer and fixed effects meta-analysis comprising 23,053 blood samples from 22 different cohorts, to comprehensively map associations between the 12 immune-cell fractions and common phenotypes. In a separate cohort of 4386 blood samples, we assess associations between immune-cell fractions and health outcomes. RESULTS: Our meta-analysis reveals many associations of cell-type fractions with age, sex, smoking and obesity, many of which we validate with single-cell RNA sequencing. We discover that naïve and regulatory T-cell subsets are higher in women compared to men, while the reverse is true for monocyte, natural killer, basophil, and eosinophil fractions. Decreased natural killer counts associated with smoking, obesity, and stress levels, while an increased count correlates with exercise and sleep. Analysis of health outcomes revealed that increased naïve CD4 + T-cell and N-cell fractions associated with a reduced risk of all-cause mortality independently of all major epidemiological risk factors and baseline co-morbidity. A machine learning predictor built only with immune-cell fractions achieved a C-index value for all-cause mortality of 0.69 (95%CI 0.67-0.72), which increased to 0.83 (0.80-0.86) upon inclusion of epidemiological risk factors and baseline co-morbidity. CONCLUSIONS: This work contributes an extensively validated high-resolution DNAm reference matrix for blood, which is made freely available, and uses it to generate a comprehensive map of associations between immune-cell fractions and common phenotypes, including health outcomes.

Elevated third trimester corticosteroid levels are associated with fewer offspring infections.

Respiratory infections are a leading cause of morbidity and mortality in early life, and recurrent infections increase the risk of developing chronic diseases. The maternal environment during pregnancy can impact offspring health, but the factors leading to increased infection proneness have not been well characterized during this period. Steroids have been implicated in respiratory health outcomes and may similarly influence infection susceptibility. Our objective was to describe relationships between maternal steroid levels and offspring infection proneness. Using adjusted Poisson regression models, we evaluated associations between sixteen androgenic and corticosteroid metabolites during pregnancy and offspring respiratory infection incidence across two pre-birth cohorts (N = 774 in VDAART and N = 729 in COPSAC). Steroid metabolites were measured in plasma samples from pregnant mothers across all trimesters of pregnancy by ultrahigh-performance-liquid-chromatography/mass-spectrometry. We conducted further inquiry into associations of steroids with related respiratory outcomes: asthma and lung function spirometry. Higher plasma corticosteroid levels in the third trimester of pregnancy were associated with lower incidence of offspring respiratory infections (P = 4.45 × 10-7 to 0.002) and improved lung function metrics (P = 0.020-0.036). Elevated maternal androgens were generally associated with increased offspring respiratory infections and worse lung function, with some associations demonstrating nominal significance at P < 0.05, but these trends were inconsistent across individual androgens. Increased maternal plasma corticosteroid levels in the late second and third trimesters were associated with lower infections and better lung function in offspring, which may represent a potential avenue for intervention through corticosteroid supplementation in late pregnancy to reduce offspring respiratory infection susceptibility in early life.Clinical Trial Registry information: VDAART and COPSAC were originally conducted as clinical trials; VDAART: ClinicalTrials.gov identifier NCT00920621; COPSAC: ClinicalTrials.gov identifier NCT00798226.

Phenotypically driven subgroups of ASD display distinct metabolomic profiles.

Autism Spectrum Disorder (ASD) is a heterogeneous condition that includes a broad range of characteristics and associated comorbidities; however, the biology underlying the variability in phenotypes is not well understood. As ASD impacts approximately 1 in 100 children globally, there is an urgent need to better understand the biological mechanisms that contribute to features of ASD. In this study, we leveraged rich phenotypic and diagnostic information related to ASD in 2001 individuals aged 4 to 17 years from the Simons Simplex Collection to derive phenotypically driven subgroups and investigate their respective metabolomes. We performed hierarchical clustering on 40 phenotypes spanning four ASD clinical domains, resulting in three subgroups with distinct phenotype patterns. Using global plasma metabolomic profiling generated by ultrahigh-performance liquid chromatography mass spectrometry, we characterized the metabolome of individuals in each subgroup to interrogate underlying biology related to the subgroups. Subgroup 1 included children with the least maladaptive behavioral traits (N = 862); global decreases in lipid metabolites and concomitant increases in amino acid and nucleotide pathways were observed for children in this subgroup. Subgroup 2 included children with the highest degree of challenges across all phenotype domains (N = 631), and their metabolome profiles demonstrated aberrant metabolism of membrane lipids and increases in lipid oxidation products. Subgroup 3 included children with maladaptive behaviors and co-occurring conditions that showed the highest IQ scores (N = 508); these individuals had increases in sphingolipid metabolites and fatty acid byproducts. Overall, these findings indicated distinct metabolic patterns within ASD subgroups, which may reflect the biological mechanisms giving rise to specific patterns of ASD characteristics. Our results may have important clinical applications relevant to personalized medicine approaches towards managing ASD symptoms.

Revealing the importance of prenatal gut microbiome in offspring neurodevelopment in humans.

BACKGROUND: It has been widely recognized that a critical time window for neurodevelopment occurs in early life and the host's gut microbiome plays an important role in neurodevelopment. Following recent demonstrations that the maternal prenatal gut microbiome influences offspring brain development in murine models, we aim to explore whether the critical time window for the association between the gut microbiome and neurodevelopment is prenatal or postnatal for human. METHODS: Here we leverage a large-scale human study and compare the associations between the gut microbiota and metabolites from mothers during pregnancy and their children with the children's neurodevelopment. Specifically, using multinomial regression integrated in Songbird, we assessed the discriminating power of the maternal prenatal and child gut microbiome for children's neurodevelopment at early life as measured by the Ages & Stages Questionnaires (ASQ). FINDINGS: We show that the maternal prenatal gut microbiome is more relevant than the children's gut microbiome to the children's neurodevelopment in the first year of life (maximum Q2 = 0.212 and 0.096 separately using the taxa at the class level). Moreover, we found that Fusobacteriia is more associated with high fine motor skills in ASQ in the maternal prenatal gut microbiota but become more associated with low fine motor skills in the infant gut microbiota (rank = 0.084 and -0.047 separately), suggesting the roles of the same taxa with respect to neurodevelopment can be opposite at the two stages of fetal neurodevelopment. INTERPRETATION: These findings shed light, especially in terms of timing, on potential therapeutic interventions to prevent neurodevelopmental disorders. FUNDING: This work was supported by the National Institutes of Health (grant numbers: R01AI141529, R01HD093761, RF1AG067744, UH3OD023268, U19AI095219, U01HL089856, R01HL141826, K08HL148178, K01HL146980), and the Charles A. King Trust Postdoctoral Fellowship.

Reduced Steroid Metabolites Identify Infection-Prone Children in Two Independent Pre-Birth Cohorts.

Recurrent respiratory infections are a leading cause of morbidity and mortality in early life, but there is no broadly accepted means to identify infection-prone children during this highly vulnerable period. In this study, we investigated associations between steroid metabolites and incident respiratory infections in two pre-birth cohorts to identify novel metabolomic signatures of early infection proneness. Children from the Vitamin D Antenatal Asthma Reduction Trial and the Copenhagen Prospective Studies on Asthma in Childhood were included, and profiling was performed on plasma samples collected at ages 1 and 6 years. Both cohorts recorded incidence of lower respiratory infections, upper respiratory infections, ear infections, and colds. Poisson regression analysis assessed the associations between 18 steroid metabolites and the total number of respiratory infections that occurred in offspring during follow-up. We found that steroid metabolites across androgenic, corticosteroid, pregnenolone, and progestin classes were reduced in children that suffered more infections, and these patterns persisted at age 6 years, generally reflecting consistency in direction of effect and significance. Our analysis suggested steroid metabolite measurement may be useful in screening for infection proneness during this critical developmental period. Future studies should clinically evaluate their potential utility as a clinical screening tool.

Circulating N-formylmethionine and metabolic shift in critical illness: a multicohort metabolomics study.

BACKGROUND: Cell stress promotes degradation of mitochondria which release danger-associated molecular patterns that are catabolized to N-formylmethionine. We hypothesized that in critically ill adults, the response to N-formylmethionine is associated with increases in metabolomic shift-related metabolites and increases in 28-day mortality. METHODS: We performed metabolomics analyses on plasma from the 428-subject Correction of Vitamin D Deficiency in Critically Ill Patients trial (VITdAL-ICU) cohort and the 90-subject Brigham and Women's Hospital Registry of Critical Illness (RoCI) cohort. In the VITdAL-ICU cohort, we analyzed 983 metabolites at Intensive Care Unit (ICU) admission, day 3, and 7. In the RoCI cohort, we analyzed 411 metabolites at ICU admission. The association between N-formylmethionine and mortality was determined by adjusted logistic regression. The relationship between individual metabolites and N-formylmethionine abundance was assessed with false discovery rate correction via linear regression, linear mixed-effects, and Gaussian graphical models. RESULTS: Patients with the top quartile of N-formylmethionine abundance at ICU admission had a significantly higher adjusted odds of 28-day mortality in the VITdAL-ICU (OR, 2.4; 95%CI 1.5-4.0; P = 0.001) and RoCI cohorts (OR, 5.1; 95%CI 1.4-18.7; P = 0.015). Adjusted linear regression shows that with increases in N-formylmethionine abundance at ICU admission, 55 metabolites have significant differences common to both the VITdAL-ICU and RoCI cohorts. With increased N-formylmethionine abundance, both cohorts had elevations in individual short-chain acylcarnitine, branched chain amino acid, kynurenine pathway, and pentose phosphate pathway metabolites. CONCLUSIONS: The results indicate that circulating N-formylmethionine promotes a metabolic shift with heightened mortality that involves incomplete mitochondrial fatty acid oxidation, increased branched chain amino acid metabolism, and activation of the pentose phosphate pathway.

Effect of perfluoroalkyl exposure in pregnancy and infancy on intrauterine and childhood growth and anthropometry. Sub study from COPSAC2010 birth cohort.

BACKGROUND: Perfluoroalkyl substances PFOS and PFOA are persistent and bioaccumulative exogenous chemicals in the human body with a range of suspected negative health effects. It is hypothesised that exposure during prenatal and early postnatal life might have particularly detrimental effects on intrauterine and childhood growth. In a Danish longitudinal mother-child cohort we investigate effect of PFOS and PFOA in pregnancy and infancy on intrauterine and childhood growth and anthropometry. METHODS: COPSAC2010 is an ongoing population based mother-child cohort of 738 pregnant women and their children followed from 24 week gestation with longitudinal deep clinical phenotyping until age 10 years. In this observational cohort sub study plasma PFOS and PFOA concentrations were semi-quantified by untargeted metabolomics in the mothers at week 24 and 1 week postpartum and in the children at ages 6 and 18 months and calibrated using a targeted pipeline. We examined associations to intrauterine and childhood growth and anthropometry, including interactions with child sex. Untargeted and targeted blood metabolomics profiles were integrated to investigate underlying mechanisms. FINDINGS: Pregnancy plasma PFOA concentrations were associated with lower birth size -0.19 [-0.33; -0.05] BMI z-score per 1-ng/mL and increased childhood height (z-scored) at age 6: 0.18 [0.05; 0.31], but there was no association between childs' own infancy plasma PFOA concentration and height. Pregnancy plasma PFOS concentrations were also associated with lower birth BMI (-0.04 [-0.08; -0.01]), but in childhood pregnancy plasma PFOS concentration interacted with child sex on BMI and fat percentage at 6 years with negative associations in girls and positive in boys. The effect of maternal plasma PFOS concentration on lower girl BMI was borderline mediated through increasing child plasma lactosyl-ceramide levels (p-mediation=0.08). Similarly the effect of maternal plasma PFOS concentration on higher boy fat percentage was borderline mediated through increasing child plasma lactosyl-ceramide levels (p-mediation=0.07). Infancy concentrations of plasma PFOS associated with lower height in childhood, -0.06 z-score at age 6 [-0.19; -0.03]. INTERPRETATION: Higher PFOS and PFOA plasma concentrations during pregnancy had detrimental effects on fetal growth. The effects on childhood growth were not similar as PFOA increased child height, opposite of PFOS in multipollutant models suggesting a differing fetal programming effect. Sex specific growth effects were borderline mediated through an altered lactosyl-ceramide metabolism, proposing a possible mechanism of PFOS that has long-lasting health consequences in this observational study. FUNDING: All funding received by COPSAC are listed on www.copsac.com. The Lundbeck Foundation (Grant no R16-A1694); The Novo Nordic Foundation (Grant nos NNF20OC0061029, NNF170C0025014, NNF180C0031764) The Ministry of Health (Grant no 903516); Danish Council for Strategic Research (Grant no 0603-00280B) and The Capital Region Research Foundation have provided core support to the COPSAC research center. Effort from JALS is supported by R01HL123915, R01HL141826, and R01HL155742 from NIH/NHLBI. CEW was supported by the Swedish Heart Lung Foundation (HLF 20180290, HLF 20200693). BC has received funding for this project from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No. 946228). The funding agencies did not have any role in design and conduct of the study; collection, management, and interpretation of the data; or preparation, review, or approval of the manuscript.

Genome-wide association study in minority children with asthma implicates DNAH5 in bronchodilator responsiveness.

Variability in response to short-acting ß2-agonists (e.g., albuterol) among patients with asthma from diverse racial/ethnic groups may contribute to asthma disparities. We sought to identify genetic variants associated with bronchodilator response (BDR) to identify potential mechanisms of drug response and risk factors for worse asthma outcomes. Genome-wide association studies of bronchodilator response (BDR) were performed using TOPMed Whole Genome Sequencing data of the Asthma Translational Genomic Collaboration (ATGC), which corresponded to 1136 Puerto Rican, 656 Mexican and 4337 African American patients with asthma. With the population-specific GWAS results, a trans-ethnic meta-analysis was performed to identify BDR-associated variants shared across the three populations. Replication analysis was carried out in three pediatric asthma cohorts, including CAMP (Childhood Asthma Management Program; n = 560), GACRS (Genetics of Asthma in Costa Rica Study; n = 967) and HPR (Hartford-Puerto Rico; n = 417). A genome-wide significant locus (rs35661809; P = 3.61 × 10-8) in LINC02220, a non-coding RNA gene, was identified in Puerto Ricans. While this region was devoid of protein-coding genes, capture Hi-C data showed a distal interaction with the promoter of the DNAH5 gene in lung tissue. In replication analysis, the GACRS cohort yielded a nominal association (1-tailed P < 0.05). No genetic variant was associated with BDR at the genome-wide significant threshold in Mexicans and African Americans. Our findings help inform genetic underpinnings of BDR for understudied minority patients with asthma, but the limited availability of genetic data for racial/ethnic minority children with asthma remains a paramount challenge.