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.

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.

Dysregulated STAT3 signaling and T cell immunometabolic dysfunction define a targetable, high mortality subphenotype of critically ill children.

Sepsis is the leading cause of death of hospitalized children worldwide. Despite the established link between immune dysregulation and mortality in pediatric sepsis, it remains unclear which host immune factors contribute causally to adverse sepsis outcomes. Identifying modifiable pathobiology is an essential first step to successful translation of biologic insights into precision therapeutics. We designed a prospective, longitudinal cohort study of 88 critically ill pediatric patients with multiple organ dysfunction syndrome (MODS), including patients with and without sepsis, to define subphenotypes associated with targetable mechanisms of immune dysregulation. We first assessed plasma proteomic profiles and identified shared features of immune dysregulation in MODS patients with and without sepsis. We then employed consensus clustering to define three subphenotypes based on protein expression at disease onset and identified a strong association between subphenotype and clinical outcome. We next identified differences in immune cell frequency and activation state by MODS subphenotype and determined the association between hyperinflammatory pathway activation and cellular immunophenotype. Using single cell transcriptomics, we demonstrated STAT3 hyperactivation in lymphocytes from the sickest MODS subgroup and then identified an association between STAT3 hyperactivation and T cell immunometabolic dysregulation. Finally, we compared proteomics findings between patients with MODS and patients with inborn errors of immunity that amplify cytokine signaling pathways to further assess the impact of STAT3 hyperactivation in the most severe patients with MODS. Overall, these results identify a potentially pathologic and targetable role for STAT3 hyperactivation in a subset of pediatric patients with MODS who have high severity of illness and poor prognosis.

Identification and transcriptomic assessment of latent profile pediatric septic shock phenotypes.

BACKGROUND: Sepsis poses a grave threat, especially among children, but treatments are limited owing to heterogeneity among patients. We sought to test the clinical and biological relevance of pediatric septic shock subclasses identified using reproducible approaches. METHODS: We performed latent profile analyses using clinical, laboratory, and biomarker data from a prospective multi-center pediatric septic shock observational cohort to derive phenotypes and trained a support vector machine model to assign phenotypes in an internal validation set. We established the clinical relevance of phenotypes and tested for their interaction with common sepsis treatments on patient outcomes. We conducted transcriptomic analyses to delineate phenotype-specific biology and inferred underlying cell subpopulations. Finally, we compared whether latent profile phenotypes overlapped with established gene-expression endotypes and compared survival among patients based on an integrated subclassification scheme. RESULTS: Among 1071 pediatric septic shock patients requiring vasoactive support on day 1 included, we identified two phenotypes which we designated as Phenotype 1 (19.5%) and Phenotype 2 (80.5%). Membership in Phenotype 1 was associated with ~ fourfold adjusted odds of complicated course relative to Phenotype 2. Patients belonging to Phenotype 1 were characterized by relatively higher Angiopoietin-2/Tie-2 ratio, Angiopoietin-2, soluble thrombomodulin (sTM), interleukin 8 (IL-8), and intercellular adhesion molecule 1 (ICAM-1) and lower Tie-2 and Angiopoietin-1 concentrations compared to Phenotype 2. We did not identify significant interactions between phenotypes, common treatments, and clinical outcomes. Transcriptomic analysis revealed overexpression of genes implicated in the innate immune response and driven primarily by developing neutrophils among patients designated as Phenotype 1. There was no statistically significant overlap between established gene-expression endotypes, reflective of the host adaptive response, and the newly derived phenotypes, reflective of the host innate response including microvascular endothelial dysfunction. However, an integrated subclassification scheme demonstrated varying survival probabilities when comparing patient endophenotypes. CONCLUSIONS: Our research underscores the reproducibility of latent profile analyses to identify pediatric septic shock phenotypes with high prognostic relevance. Pending validation, an integrated subclassification scheme, reflective of the different facets of the host response, holds promise to inform targeted intervention among those critically ill.

Urine metabolomics signature reveals novel determinants of adrenal suppression in children taking inhaled corticosteroids to control asthma symptoms.

BACKGROUND: Asthma is routinely treated with inhaled corticosteroids (ICS). Asthma patients on ICS are at increased risk of adrenal suppression, a potentially serious effect of long-term glucocorticoid exposure; however, this relationship is poorly understood. Therefore, this study aims to identify metabolite biomarkers related to adrenal suppression in asthma patients taking ICS. METHODS: A total of 571 urine metabolites from 200 children with asthma on ICS in the Pharmacogenetics of Adrenal Suppression with Inhaled Steroids (PASS) cohort were profiled. Samples were grouped by peak plasma cortisol measurement as adrenal sufficient (>350 nmol/L) or insufficient (≤350 nmol/L) (outcome). Regression and discriminant-based statistical models combined with network analyses were utilized to assess relationships between metabolites and the outcome. Finally, prioritized metabolites were validated using data from an ancillary study of the Childhood Asthma Management (CAMP) cohort with similar characteristics to PASS. RESULTS: Ninety metabolites were significantly associated with adrenal suppression, of which 57 also could discriminate adrenal status. While 26 metabolites (primarily steroids) were present at lower levels in the adrenal insufficient patients, 14 were significantly elevated in this group; the top metabolite, mannitol/sorbitol, was previously associated with asthma exacerbations. Network analyses identified unique clusters of metabolites related to steroids, fatty acid oxidation, and nucleoside metabolism, respectively. Four metabolites including urocanic acid, acetylcarnitine, uracil, and sorbitol were validated in CAMP cohort for adrenal suppression. CONCLUSIONS: Urinary metabolites differ among asthma patients on ICS, by adrenal status. While steroid metabolites were reduced in patients with poor adrenal function, our findings also implicate previously unreported metabolites involved in amino acid, lipid, and nucleoside metabolism.

The metabolic role of vitamin D in children's neurodevelopment: a network study.

Neurodevelopmental disorders are rapidly increasing in prevalence and have been linked to various environmental risk factors. Mounting evidence suggests a potential role of vitamin D in child neurodevelopment, though the causal mechanisms remain largely unknown. Here, we investigate how vitamin D deficiency affects children's communication development, particularly in relation to Autism Spectrum Disorder (ASD). We do so by developing an integrative network approach that combines metabolomic profiles, clinical traits, and neurodevelopmental data from a pediatric cohort. Our results show that low levels of vitamin D are associated with changes in the metabolic networks of tryptophan, linoleic, and fatty acid metabolism. These changes correlate with distinct ASD-related phenotypes, including delayed communication skills and respiratory dysfunctions. Additionally, our analysis suggests the kynurenine and serotonin sub-pathways may mediate the effect of vitamin D on early life communication development. Altogether, our findings provide metabolome-wide insights into the potential of vitamin D as a therapeutic option for ASD and other communication disorders.

Urinary eicosanoid levels in early life and risk of atopic disease in childhood.

BACKGROUND: Eicosanoids are lipid mediators including thromboxanes (TXs), prostaglandins (PGs), and leukotrienes with a pathophysiological role in established atopic disease. However, their role in the inception of disease is unclear. This study aimed to investigate the association between urinary eicosanoids in early life and development of atopic disease. METHODS: This study quantified the levels of 21 eicosanoids in urine from children from the COPSAC2010 (Copenhagen Prospective Studies on Asthma in Childhood 2010) (age 1 year, n = 450) and VDAART (Vitamin D Antenatal Asthma Reduction Trial) (age 3 years, n = 575) mother-child cohorts and analyzed the associations with development of wheeze/asthma, atopic dermatitis, and biomarkers of type-2 inflammation, applying false discovery rate of 5% (FDR5%) multiple testing correction. RESULTS: In both cohorts, analyses adjusted for environmental determinants showed that higher TXA2 eicosanoids in early life were associated with increased risk of developing atopic dermatitis (P < FDR5%) and type-2 inflammation (P < .05). In VDAART, lower PGE2 and PGI2 eicosanoids and higher isoprostanes were also associated with increased risk of atopic dermatitis (P < FDR5%). For wheeze/asthma, analyses in COPSAC2010 showed that lower isoprostanes and PGF2 eicosanoids and higher PGD2 eicosanoids at age 1 year associated with an increased risk at age 1-10 years (P < .05), whereas analyses in VDAART showed that lower PGE2 and higher TXA2 eicosanoids at age 3 years associated with an increased risk at 6 years (P < FDR5%). CONCLUSIONS: This study suggests that early life perturbations in the eicosanoid metabolism are present before the onset of atopic disease in childhood, which provides pathophysiological insight in the inception of atopic diseases.

Family history as the strongest predictor of aortic and peripheral aneurysms in patients with intracranial aneurysms.

OBJECTIVE: Intracranial aneurysms (IA) and aortic aneurysms (AA) are both abnormal dilations of arteries with familial predisposition and have been proposed to share co-prevalence and pathophysiology. Associations of IA and non-aortic peripheral aneurysms are less well-studied. The goal of the study was to understand the patterns of aortic and peripheral (extracranial) aneurysms in patients with IA, and risk factors associated with the development of these aneurysms. METHODS: 4701 patients were included in our retrospective analysis of all patients with intracranial aneurysms at our institution over the past 26 years. Patient demographics, comorbidities, and aneurysmal locations were analyzed. Univariate and multivariate analyses were performed to study associations with and without extracranial aneurysms. RESULTS: A total of 3.4% of patients (161 of 4701) with IA had at least one extracranial aneurysm. 2.8% had thoracic or abdominal aortic aneurysms. Age, male sex, hypertension, coronary artery disease, history of ischemic cerebral infarction, connective tissues disease, and family history of extracranial aneurysms in a 1st degree relative were associated with the presence of extracranial aneurysms and a higher number of extracranial aneurysms. In addition, family history of extracranial aneurysms in a second degree relative is associated with the presence of extracranial aneurysms and atrial fibrillation is associated with a higher number of extracranial aneurysms. CONCLUSION: Significant comorbidities are associated with extracranial aneurysms in patients with IA. Family history of extracranial aneurysms has the strongest association and suggests that IA patients with a family history of extracranial aneurysms may benefit from screening.

Consistent Multi-Omic Relationships Uncover Molecular Basis of Pediatric Asthma IgE Regulation.

Serum total immunoglobulin E levels (total IgE) capture the state of the immune system in relation to allergic sensitization. High levels are associated with airway obstruction and poor clinical outcomes in pediatric asthma. Inconsistent patient response to anti-IgE therapies motivates discovery of molecular mechanisms underlying serum IgE level differences in children with asthma. To uncover these mechanisms using complementary metabolomic and transcriptomic data, abundance levels of 529 named metabolites and expression levels of 22,772 genes were measured among children with asthma in the Childhood Asthma Management Program (CAMP, N=564) and the Genetic Epidemiology of Asthma in Costa Rica Study (GACRS, N=309) via the TOPMed initiative. Gene-metabolite associations dependent on IgE were identified within each cohort using multivariate linear models and were interpreted in a biochemical context using network topology, pathway and chemical enrichment, and representation within reactions. A total of 1,617 total IgE-dependent gene-metabolite associations from GACRS and 29,885 from CAMP met significance cutoffs. Of these, glycine and guanidinoacetic acid (GAA) were associated with the most genes in both cohorts, and the associations represented reactions central to glycine, serine, and threonine metabolism and arginine and proline metabolism. Pathway and chemical enrichment analysis further highlighted additional related pathways of interest. The results of this study suggest that GAA may modulate total IgE levels in two independent pediatric asthma cohorts with different characteristics, supporting the use of L-Arginine as a potential therapeutic for asthma exacerbation. Other potentially new targetable pathways are also uncovered.

Invited Commentary: Raising a High-Pressure Alarm About Pediatric Hypertension.

Hypertension is a common "silent killer" in adult medicine, but epidemiologic estimates of elevated blood pressure in children and adolescents are challenged by under-diagnosis and resultant low utilization of relevant administrative or billing codes. In the article by Horgan et al (Am J Epidemiol 2024), children and adolescents with hypertension and elevated blood pressure were identified using direct assessment of blood pressure measurements available in the electronic health record from both inpatient and outpatient visits ("clinical cohort") in comparison to diagnosis codes ("claims-based cohort"). The study population included 3.75 million pediatric healthcare visits available in the US Food and Drug Administration's Sentinel System. While the study applied a relatively novel methodology to interrogate available clinical data within the EHR to better understand the prevalence of pediatric hypertension and raised concern for a higher occurrence of hypertension among children and adolescents than previously realized using claims codes, the utility of the prevalence estimates may be limited by the potential for misclassification bias inherent in EHR data. However, these data raise important concerns about relaying solely on ICD-9-CM/ICD-10-CM codes to quantify the epidemiology of pediatric hypertension and highlight opportunities to address elevated blood pressure in children that could improve long-term cardiovascular health.

Diagnostic Identification of Acute Brain Dysfunction in Pediatric Sepsis and Septic Shock in the Electronic Health Record: A Comparison of Four Definitions in a Reference Dataset.

OBJECTIVES: Acute brain dysfunction (ABD) in pediatric sepsis has a prevalence of 20%, but can be difficult to identify. Our previously validated ABD computational phenotype (CPABD) used variables obtained from the electronic health record indicative of clinician concern for acute neurologic or behavioral change. We tested whether the CPABD has better diagnostic performance to identify confirmed ABD than other definitions using the Glasgow Coma Scale or delirium scores. DESIGN: Diagnostic testing in a curated cohort of pediatric sepsis/septic shock patients. SETTING: Quaternary freestanding children's hospital. SUBJECTS: The test dataset comprised 527 children with sepsis/septic shock managed between 2011 and 2021 with a prevalence (pretest probability) of confirmed ABD of 30% (159/527). MEASUREMENTS AND MAIN RESULTS: CPABD was based on use of neuroimaging, electroencephalogram, and/or administration of new antipsychotic medication. We compared the performance of the CPABD with three GCS/delirium-based definitions of ABD-Proulx et al, International Pediatric Sepsis Consensus Conference, and Pediatric Organ Dysfunction Information Update Mandate. The posttest probability of identifying ABD was highest in CPABD (0.84) compared with other definitions. CPABD also had the highest sensitivity (83%; 95% CI, 76-89%) and specificity (93%; 95% CI, 90-96%). The false discovery rate was lowest in CPABD (1-in-6) as was the false omission rate (1-in-14). Finally, the prevalence threshold for the definitions varied, with the CPABD being the definition closest to 20%. CONCLUSIONS: In our curated dataset of pediatric sepsis/septic shock, CPABD had favorable characteristics to identify confirmed ABD compared with GCS/delirium-based definitions. The CPABD can be used to further study the impact of ABD in studies using large electronic health datasets.

The Asthma Risk Gene, GSDMB, Promotes Mitochondrial DNA-induced ISGs Expression.

Released mitochondrial DNA (mtDNA) in cells activates cGAS-STING pathway, which induces expression of interferon-stimulated genes (ISGs) and thereby promotes inflammation, as frequently seen in asthmatic airways. However, whether the genetic determinant, Gasdermin B (GSDMB), the most replicated asthma risk gene, regulates this pathway remains unknown. We set out to determine whether and how GSDMB regulates mtDNA-activated cGAS-STING pathway and subsequent ISGs induction in human airway epithelial cells. Using qPCR, ELISA, native polyacrylamide gel electrophoresis, co-immunoprecipitation and immunofluorescence assays, we evaluated the regulation of GSDMB on cGAS-STING pathway in both BEAS-2B cells and primary normal human bronchial epithelial cells (nHBEs). mtDNA was extracted in plasma samples from human asthmatics and the correlation between mtDNA levels and eosinophil counts was analyzed. GSDMB is significantly associated with RANTES expression in asthmatic nasal epithelial brushing samples from the Genes-environments and Admixture in Latino Americans (GALA) II study. Over-expression of GSDMB promotes DNA-induced IFN and ISGs expression in bronchial epithelial BEAS-2B cells and nHBEs. Conversely, knockout of GSDMB led to weakened induction of interferon (IFNs) and ISGs in BEAS-2B cells. Mechanistically, GSDMB interacts with the C-terminus of STING, promoting the translocation of STING to Golgi, leading to the phosphorylation of IRF3 and induction of IFNs and ISGs. mtDNA copy number in serum from asthmatics was significantly correlated with blood eosinophil counts especially in male subjects. GSDMB promotes the activation of mtDNA and poly (dA:dT)-induced activation of cGAS-STING pathway in airway epithelial cells, leading to enhanced induction of ISGs.

Vitamin D constrains inflammation by modulating the expression of key genes on Chr17q12-21.1.

Vitamin D possesses immunomodulatory functions and vitamin D deficiency has been associated with the rise in chronic inflammatory diseases, including asthma (Litonjua and Weiss, 2007). Vitamin D supplementation studies do not provide insight into the molecular genetic mechanisms of vitamin D-mediated immunoregulation. Here, we provide evidence for vitamin D regulation of two human chromosomal loci, Chr17q12-21.1 and Chr17q21.2, reliably associated with autoimmune and chronic inflammatory diseases. We demonstrate increased vitamin D receptor (Vdr) expression in mouse lung CD4+ Th2 cells, differential expression of Chr17q12-21.1 and Chr17q21.2 genes in Th2 cells based on vitamin D status and identify the IL-2/Stat5 pathway as a target of vitamin D signaling. Vitamin D deficiency caused severe lung inflammation after allergen challenge in mice that was prevented by long-term prenatal vitamin D supplementation. Mechanistically, vitamin D induced the expression of the Ikzf3-encoded protein Aiolos to suppress IL-2 signaling and ameliorate cytokine production in Th2 cells. These translational findings demonstrate mechanisms for the immune protective effect of vitamin D in allergic lung inflammation with a strong molecular genetic link to the regulation of both Chr17q12-21.1 and Chr17q21.2 genes and suggest further functional studies and interventional strategies for long-term prevention of asthma and other autoimmune disorders.

Diagnosis and Severity Assessment of COPD Using a Novel Fast-Response Capnometer and Interpretable Machine Learning.

INTRODUCTION: Spirometry is the gold standard for COPD diagnosis and severity determination, but is technique-dependent, nonspecific, and requires administration by a trained healthcare professional. There is a need for a fast, reliable, and precise alternative diagnostic test. This study's aim was to use interpretable machine learning to diagnose COPD and assess severity using 75-second carbon dioxide (CO2) breath records captured with TidalSense's N-TidalTM capnometer. METHOD: For COPD diagnosis, machine learning algorithms were trained and evaluated on 294 COPD (including GOLD stages 1-4) and 705 non-COPD participants. A logistic regression model was also trained to distinguish GOLD 1 from GOLD 4 COPD with the output probability used as an index of severity. RESULTS: The best diagnostic model achieved an AUROC of 0.890, sensitivity of 0.771, specificity of 0.850 and positive predictive value (PPV) of 0.834. Evaluating performance on all test capnograms that were confidently ruled in or out yielded PPV of 0.930 and NPV of 0.890. The severity determination model yielded an AUROC of 0.980, sensitivity of 0.958, specificity of 0.961 and PPV of 0.958 in distinguishing GOLD 1 from GOLD 4. Output probabilities from the severity determination model produced a correlation of 0.71 with percentage predicted FEV1. CONCLUSION: The N-TidalTM device could be used alongside interpretable machine learning as an accurate, point-of-care diagnostic test for COPD, particularly in primary care as a rapid rule-in or rule-out test. N-TidalTM also could be effective in monitoring disease progression, providing a possible alternative to spirometry for disease monitoring.

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).

Machine Learning Prediction of Treatment Response to Inhaled Corticosteroids in Asthma.

BACKGROUND: Although inhaled corticosteroids (ICS) are the first-line therapy for patients with persistent asthma, many patients continue to have exacerbations. We developed machine learning models to predict the ICS response in patients with asthma. METHODS: The subjects included asthma patients of European ancestry (n = 1371; 448 children; 916 adults). A genome-wide association study was performed to identify the SNPs associated with ICS response. Using the SNPs identified, two machine learning models were developed to predict ICS response: (1) least absolute shrinkage and selection operator (LASSO) regression and (2) random forest. RESULTS: The LASSO regression model achieved an AUC of 0.71 (95% CI 0.67-0.76; sensitivity: 0.57; specificity: 0.75) in an independent test cohort, and the random forest model achieved an AUC of 0.74 (95% CI 0.70-0.78; sensitivity: 0.70; specificity: 0.68). The genes contributing to the prediction of ICS response included those associated with ICS responses in asthma (TPSAB1, FBXL16), asthma symptoms and severity (ABCA7, CNN2, PTRN3, and BSG/CD147), airway remodeling (ELANE, FSTL3), mucin production (GAL3ST), leukotriene synthesis (GPX4), allergic asthma (ZFPM1, SBNO2), and others. CONCLUSIONS: An accurate risk prediction of ICS response can be obtained using machine learning methods, with the potential to inform personalized treatment decisions. Further studies are needed to examine if the integration of richer phenotype data could improve risk prediction.

Circulating miRNAs associate with historical childhood asthma hospitalization in different serum vitamin D groups.

BACKGROUND: Vitamin D may help to alleviate asthma exacerbation because of its anti-inflammation effect, but the evidence is inconsistent in childhood asthma. MiRNAs are important mediators in asthma pathogenesis and also excellent non-invasive biomarkers. We hypothesized that circulating miRNAs are associated with asthma exacerbation and modified by vitamin D levels. METHODS: We sequenced baseline serum miRNAs from 461 participants in the Childhood Asthma Management Program (CAMP). Logistic regression was used to associate miRNA expression with asthma exacerbation through interaction analysis first and then stratified by vitamin D insufficient and sufficient groups. Microarray from lymphoblastoid B-cells (LCLs) treated by vitamin D or sham of 43 subjects in CAMP were used for validation in vitro. The function of miRNAs was associated with gene modules by weighted gene co-expression network analysis (WGCNA). RESULTS: We identified eleven miRNAs associated with asthma exacerbation with vitamin D effect modification. Of which, five were significant in vitamin D insufficient group and nine were significant in vitamin D sufficient group. Six miRNAs, including hsa-miR-143-3p, hsa-miR-192-5p, hsa-miR-151a-5p, hsa-miR-24-3p, hsa-miR-22-3p and hsa-miR-451a were significantly associated with gene modules of immune-related functions, implying miRNAs may mediate vitamin D effect on asthma exacerbation through immune pathways. In addition, hsa-miR-143-3p and hsa-miR-451a are potential predictors of childhood asthma exacerbation at different vitamin D levels. CONCLUSIONS: miRNAs are potential mediators of asthma exacerbation and their effects are directly impacted by vitamin D levels.

Gasdermin B, an asthma-susceptibility gene, promotes MAVS-TBK1 signalling and airway inflammation.

RATIONALE: Respiratory virus-induced inflammation is the leading cause of asthma exacerbation, frequently accompanied by induction of interferon-stimulated genes (ISGs). How asthma-susceptibility genes modulate cellular response upon viral infection by fine-tuning ISG induction and subsequent airway inflammation in genetically susceptible asthma patients remains largely unknown. OBJECTIVES: To decipher the functions of gasdermin B (encoded by GSDMB) in respiratory virus-induced lung inflammation. METHODS: In two independent cohorts, we analysed expression correlation between GSDMB and ISG s. In human bronchial epithelial cell line or primary bronchial epithelial cells, we generated GSDMB-overexpressing and GSDMB-deficient cells. A series of quantitative PCR, ELISA and co-immunoprecipitation assays were performed to determine the function and mechanism of GSDMB for ISG induction. We also generated a novel transgenic mouse line with inducible expression of human unique GSDMB gene in airway epithelial cells and infected the mice with respiratory syncytial virus to determine the role of GSDMB in respiratory syncytial virus-induced lung inflammation in vivo. RESULTS: GSDMB is one of the most significant asthma-susceptibility genes at 17q21 and acts as a novel RNA sensor, promoting mitochondrial antiviral-signalling protein (MAVS)-TANK binding kinase 1 (TBK1) signalling and subsequent inflammation. In airway epithelium, GSDMB is induced by respiratory viral infections. Expression of GSDMB and ISGs significantly correlated in respiratory epithelium from two independent asthma cohorts. Notably, inducible expression of human GSDMB in mouse airway epithelium led to enhanced ISGs induction and increased airway inflammation with mucus hypersecretion upon respiratory syncytial virus infection. CONCLUSIONS: GSDMB promotes ISGs expression and airway inflammation upon respiratory virus infection, thereby conferring asthma risk in risk allele carriers.

Joint clinical and molecular subtyping of COPD with variational autoencoders.

Chronic Obstructive Pulmonary Disease (COPD) is a complex, heterogeneous disease. Traditional subtyping methods generally focus on either the clinical manifestations or the molecular endotypes of the disease, resulting in classifications that do not fully capture the disease's complexity. Here, we bridge this gap by introducing a subtyping pipeline that integrates clinical and gene expression data with variational autoencoders. We apply this methodology to the COPDGene study, a large study of current and former smoking individuals with and without COPD. Our approach generates a set of vector embeddings, called Personalized Integrated Profiles (PIPs), that recapitulate the joint clinical and molecular state of the subjects in the study. Prediction experiments show that the PIPs have a predictive accuracy comparable to or better than other embedding approaches. Using trajectory learning approaches, we analyze the main trajectories of variation in the PIP space and identify five well-separated subtypes with distinct clinical phenotypes, expression signatures, and disease outcomes. Notably, these subtypes are more robust to data resampling compared to those identified using traditional clustering approaches. Overall, our findings provide new avenues to establish fine-grained associations between the clinical characteristics, molecular processes, and disease outcomes of COPD.

Lung function trajectories in a cohort of patients with moderate-to-severe asthma on mepolizumab, omalizumab, or dupilumab.

BACKGROUND: Lung function is an independent predictor of mortality. We evaluated the lung function trajectories of a cohort of patients with asthma receiving biologic therapy. METHODS: We identified 229 monoclonal antibody-naïve adult patients with moderate-to-severe asthma who initiated omalizumab, mepolizumab, or dupilumab between 2010 and 2022 in a large healthcare system in Boston, MA. Generalized additive mixed models were used to estimate the lung function trajectories during the 156 weeks following biologic initiation. Response was defined as an improvement in FEV1 or a decrease of ≤0.5% per year. The Kaplan-Meier estimator was used to assess time to no additional improvement in FEV1 in responders. All models were adjusted for age, sex, body mass index, smoking status, baseline exacerbation rate, and baseline blood eosinophil count. RESULTS: Eighty-eight patients initiated mepolizumab, 76 omalizumab, and 65 dupilumab. Baseline eosinophil count was highest in the mepolizumab group (405 cells/mcL) and lowest for omalizumab (250 cells/mcL). Both FEV1 and FVC improved in the mepolizumab group (FEV1 + 20 mL/year; FVC +43 mL/year). For omalizumab, there was an initial improvement in the first year followed by decline with an overall FEV1 loss of -44 mL/year and FVC -32 mL/year. For dupilumab, both FEV1 (+61 mL/year) and FVC (+74 mL/year) improved over time. Fifty percent of the mepolizumab group, 58% omalizumab, and 72% of dupilumab were responders. The median time to no additional FEV1 improvement in responders was 24 weeks for omalizumab, 48 weeks for mepolizumab, and 57 weeks for dupilumab. CONCLUSION: In this clinical cohort, mepolizumab, omalizumab, and dupilumab had beneficial effects on FEV1 and FVC with distinct post-initiation trajectories.