Search

Background: MicroRNAs (miRNAs) are crucial post-transcriptional regulators involved in inflammatory diseases, such as asthma. Poor lung function and airflow issues in childhood are linked to the development of chronic obstructive pulmonary disease (COPD) in adulthood., Methods: We analyzed small RNA-Seq data from 365 peripheral whole blood samples from the Genetics of Asthma in Costa Rica Study (GACRS) for association with airflow levels measured by FEV1/FVC. Differentially expressed (DE) miRNAs were identified using DESeq2 in R, adjusting for covariates and applying a 10% false discovery rate (FDR). The analysis included 361 samples and 649 miRNAs. The two DE miRNAs were further tested for association with airflow obstruction in a study of adult former smokers with and without COPD., Results: We found 1 upregulated and 1 downregulated miRNA in participants with airflow below the threshold compared to those above it. In the adult study, the same miRNAs were upregulated and downregulated in individuals with FEV1/FVC < 0.7 versus those with FEV1/FVC > 0.7, showing suggestive statistical evidence. The target genes of these miRNAs were enriched for PI3K-Akt, Hippo, WNT, MAPK, and focal adhesion pathways., Conclusions: Two differentially expressed miRNAs were associated with airflow levels in children with asthma and airflow obstruction in adults with COPD. This suggests that shared genetic regulatory systems may influence childhood airflow and contribute to adulthood airflow obstruction., Competing Interests: Declarations. Institutional review board statement: This study was conducted in accordance with the Declaration of Helsinki and was approved by the Institutional Review Board of Brigham and Women’s Hospital, IRB# 2017P001799, 7/29/2020. Additional approval for the Costa Rica cohort was obtained from the Hospital Nacional de Niños (San José, Costa Rica) and Brigham and Women’s Hospital (Boston, MA, USA). COPDGene Institutional Review Board (IRB) approval was obtained at each of the participating study centers before study initiation. Consent for publication: Not applicable. Clinical trial number: Not applicable. Competing interests: The authors declare no competing interests., (© 2025. The Author(s).)

Introduction: Uncontrolled or severe asthma results in symptomatic usage of short-acting ß2-agonist (SABA) usage. MicroRNAs (miRNAs) are posttranslational regulators that can influence asthma biology. This study aims to identify miRNAs that are associated with increased SABA usage., Methods: Small RNA sequenced from blood serum of 1,132 children with asthma aged 6 to 14 years in the Genetics of Asthma in Costa Rica Study was used for this analysis. Logistic regression identified miRNAs in patients who required increased SABA usage. These miRNA were validated for association with SABA-induced bronchodilator responsiveness (BDR). Gene target pathway analysis was performed on validated miRNAs., Results: Twenty-one miRNAs were significantly associated with increased SABA usage with OR ranging from 0.87 to 1.23. Two miRNAs, miR-378a-3p and miR-144-3p, had odds ratio 1.14 (1-1.29, p = .05) and 1.11 (1.01-1.22, p = .035), respectively, for increased SABA usage and were also significantly associated with bronchodilator response. Furthermore, a linear regression analysis involving these miRNA and bronchodilator response revealed that increased miR-378a-3p correlated with decreased BDR and increased expression of miR-144-3p correlated with improving pulmonary function with bronchodilators. In gene target Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, the erythroblastosis viral oncogene signaling pathway had among one of the highest fold enrichment and p-value., Conclusion: Increased expression of miR-378a-3p and miR-144-3p was seen in this patient population who required increased SABA usage. There were different bronchodilatory effects seen in these two miRNAs, suggesting different potential mechanisms underlying increased SABA usage., (© 2024 The Author(s). Pediatric Pulmonology published by Wiley Periodicals LLC.)

Background: Bronchodilator response (BDR) is a measure of improvement in airway smooth muscle tone, inhibition of liquid accumulation and mucus section into the lumen in response to short-acting beta-2 agonists that varies among asthmatic patients. MicroRNAs (miRNAs) are well-known post-translational regulators. Identifying miRNAs associated with BDR could lead to a better understanding of the underlying complex pathophysiology., Objective: The purpose of this study is to identify circulating miRNAs associated with bronchodilator response in asthma and decipher possible mechanism of bronchodilator response variation., Methods: We used available small RNA sequencing on blood serum from 1,134 asthmatic children aged 6 to 14 years who participated in the Genetics of Asthma in Costa Rica Study (GACRS). We filtered the participants into the highest and lowest bronchodilator response (BDR) quartiles and used DeSeq2 to identify miRNAs with differential expression (DE) in high (N = 277) vs. low (N = 278) BDR group. Replication was carried out in the Leukotriene modifier Or Corticosteroids or Corticosteroid-Salmeterol trial (LOCCS), an adult asthma cohort. The putative target genes of DE miRNAs were identified, and pathway enrichment analysis was performed., Results: We identified 10 down-regulated miRNAs having odds ratios (OR) between 0.37 and 0.76 for a doubling of miRNA counts and one up-regulated miRNA (OR = 2.26) between high and low BDR group. These were assessed for replication in the LOCCS cohort, where two miRNAs (miR-200b-3p and miR-1246) were associated. Further, functional annotation of 11 DE miRNAs were performed as well as of two replicated miRs. Target genes of these miRs were enriched in regulation of cholesterol biosynthesis by SREBPs, ESR-mediated signaling, G1/S transition, RHO GTPase cycle, and signaling by TGFB family pathways., Conclusion: MiRNAs miR-1246 and miR-200b-3p are associated with both childhood and adult asthma BDR. Our findings add to the growing body of evidence that miRNAs play a significant role in the difference of asthma treatment response among patients as it points to genomic regulatory machinery underlying difference in bronchodilator response among patients., Trial Registration: LOCCS cohort [ClinicalTrials.gov number NCT00156819, Registration date 20050912], GACRS cohort [ClinicalTrials.gov number NCT00021840]., (© 2024. The Author(s).)

The mechanisms underlying severe allergic asthma are complex and unknown, meaning it is a challenge to provide the most appropriate treatment. This study aimed to identify novel biomarkers for stratifying allergic asthmatic patients according to severity, and to uncover the biological mechanisms that lead to the development of the severe uncontrolled phenotype. By using miRNA PCR panels, we analyzed the expression of 752 miRNAs in serum samples from control subjects ( n = 15) and mild ( n = 11) and severe uncontrolled ( n = 10) allergic asthmatic patients. We identified 40 differentially expressed miRNAs between severe uncontrolled and mild allergic asthmatic patients. Functional enrichment analysis revealed signatures related to inflammation, angiogenesis, lipid metabolism and mRNA regulation. A random forest classifier trained with DE miRNAs achieved a high accuracy of 97% for severe uncontrolled patient stratification. Validation of the identified biomarkers was performed on a subset of allergic asthmatic patients from the CAMP cohort at Brigham and Women's Hospital, Harvard Medical School. Four of these miRNAs (hsa-miR-99b-5p, hsa-miR-451a, hsa-miR-326 and hsa-miR-505-3p) were validated, pointing towards their potential as biomarkers for stratifying allergic asthmatic patients by severity and providing insights into severe uncontrolled asthma molecular pathways.

MicroRNAs (miRNAs) are recognized as key regulatory factors in numerous human diseases, with the same miRNA often involved in several diseases simultaneously or being identified as a biomarker for dozens of separate diseases. While of evident biological importance, miRNA pleiotropy remains poorly understood, and quantifying this could greatly aid in understanding the broader role miRNAs play in health and disease. To this end, we introduce miRAIDD (miRNA Artificial Intelligence Disease Database), a comprehensive database of human miRNA-disease causal associations constructed using large language models (LLM). Through this endeavor, we provide two entirely novel contributions: 1) we systematically quantify miRNA pleiotropy, a property of evident translational importance; and 2) describe biological and bioinformatic characteristics of miRNAs which lead to increased pleiotropy. Further, we provide our code, database, and experience using AI LLMs to the broader research community., Competing Interests: Competing Interests The authors declare that they have no competing interests.

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., (© 2024 The Authors. Immunity, Inflammation and Disease published by John Wiley & Sons Ltd.)

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., (© 2024. The Author(s).)

Background: Piwi-interacting RNAs (piRNAs), comprising the largest noncoding RNA group, regulate transcriptional processes. Whether piRNAs are associated with type 2 (T2)-high asthma is unknown., Objective: We sought to investigate the association between piRNAs and T2-high asthma in childhood asthma., Methods: We sequenced plasma samples from 462 subjects in the Childhood Asthma Management Program (CAMP) as the discovery cohort and 1165 subjects in the Genetics of Asthma in Costa Rica Study (GACRS) as a replication cohort. Sequencing reads were filtered first, and piRNA reads were annotated and normalized. Linear regression was used for the association analysis of piRNAs and peripheral blood eosinophil count, total serum IgE level, and long-term asthma exacerbation in children with asthma. Mediation analysis was performed to investigate the effect direction. We then ascertained if the circulating piRNAs were present in asthmatic airway epithelial cells in a Gene Expression Omnibus (GEO; www.ncbi.nlm.nih.gov/geo) public data set., Results: Fifteen piRNAs were significantly associated with eosinophil count in CAMP (P ≤ .05), and 3 were successfully replicated in GACRS. Eleven piRNAs were associated with total IgE in CAMP, and one of these was replicated in GACRS. All 22 significant piRNAs were identified in epithelial cells in vitro, and 6 of these were differentially expressed between subjects with asthma and healthy controls. Fourteen piRNAs were associated with long-term asthma exacerbation, and effect of piRNAs on long-term asthma exacerbation are mediated through eosinophil count and serum IgE level., Conclusion: piRNAs are associated with peripheral blood eosinophils and total serum IgE in childhood asthma and may play important roles in T2-high asthma., (Copyright © 2023 American Academy of Allergy, Asthma & Immunology. All rights reserved.)

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 (PRS Asthma ) and spirometry (FEV 1 and FEV 1 /forced vital capacity; PRS spiro ) 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 PRS Asthma and a previously published PRS spiro 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 PRS Asthma (odds ratio [OR] per SD, 1.16 [95% CI, 1.14-1.19]) and PRS spiro (OR per SD, 1.19 [95% CI, 1.17-1.22]) both predicted asthma, whereas the PRS spiro predicted COPD (OR per SD, 1.25 [95% CI, 1.21-1.30]). However, results differed by cohort. The PRS spiro was not associated with COPD in GERA and Mass General Brigham Biobank. In the Genetic Epidemiology of COPD study, the PRS Asthma (OR per SD: Whites, 1.3; African Americans, 1.2) and PRS spiro (OR per SD: Whites, 2.2; African Americans, 1.6) were both associated with ACO. In GERA, the PRS Asthma was associated with asthma exacerbations (OR, 1.18) in Whites; the PRS spiro 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., (Published by Elsevier Inc.)

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., Competing Interests: Conflicts of interest. This work was funded in part by TruDiagnostic. JLS is a scientific advisor to Precision Inc and Ahara Inc. MM and VNG have filed patents on measuring cellular aging .STW receives royalties from UpToDate and is on the Board of Histolix, a digital pathology company.

Rationale: Bronchodilator response (BDR) is a measure of improvement in airway smooth muscle tone, inhibition of liquid accumulation and mucus section into the lumen in response to short-acting beta-2 agonists that varies among asthmatic patients. MicroRNAs (miRNAs) are well-known post-translational regulators. Identifying miRNAs associated with BDR could lead to a better understanding of the underlying complex pathophysiology., Objective: The purpose of this study is to identify circulating miRNAs associated with bronchodilator response in asthma and decipher possible mechanism of bronchodilator response variation., Methods: We used available small RNA sequencing on blood serum from 1,134 asthmatic children aged 6 to 14 years who participated in the Genetics of Asthma in Costa Rica Study (GACRS). We filtered the participants into high and low bronchodilator response (BDR) quartiles and used DeSeq2 to identify miRNAs with differential expression (DE) in high (N= 277) vs low (N= 278) BDR group. Replication was carried out in the Leukotriene modifier Or Corticosteroids or Corticosteroid-Salmeterol trial (LOCCS), an adult asthma cohort. The putative target genes of DE miRNAs were identified, and pathway enrichment analysis was performed., Results: We identified 10 down-regulated miRNAs having odds ratios (OR) between 0.37 and 0.76 for a doubling of miRNA counts and one up-regulated miRNA (OR=2.26) between high and low BDR group. These were assessed for replication in the LOCCS cohort, where two miRNAs (miR-200b-3p and miR-1246) were associated. Further, functional annotation of 11 DE miRNAs were performed as well as of two replicated miRs. Target genes of these miRs were enriched in regulation of cholesterol biosynthesis by SREBPs, ESR-mediated signaling, G1/S transition, RHO GTPase cycle, and signaling by TGFB family pathways., Conclusion: MiRNAs miR-1246 and miR-200b-3p are associated with both childhood and adult asthma BDR. Our findings add to the growing body of evidence that miRNAs play a significant role in the difference of asthma treatment response among patients as it points to genomic regulatory machinery underlying difference in bronchodilator response among patients., Trial Registration: LOCCS cohort [ClinicalTrials.gov number: NCT00156819], GACRS cohort [ClinicalTrials.gov number: NCT00021840]., Competing Interests: Competing interests The authors declare no conflict of interest in relation to this manuscript.

The association between neutrophil extracellular traps (NETs) and response to inhaled corticosteroids (ICS) in asthma is unclear. To better understand this relationship, we analyzed the blood transcriptomes from children with controlled and uncontrolled asthma in the Taiwanese Consortium of Childhood Asthma Study using weighted gene coexpression network analysis and pathway enrichment methods. We identified 298 uncontrolled asthma-specific differentially expressed genes and one gene module associated with neutrophil-mediated immunity, highlighting a potential role for neutrophils in uncontrolled asthma. We also found that NET abundance was associated with nonresponse to ICS in patients. In a neutrophilic airway inflammation murine model, steroid treatment could not suppress neutrophilic inflammation and airway hyperreactivity. However, NET disruption with deoxyribonuclease I (DNase I) efficiently inhibited airway hyperreactivity and inflammation. Using neutrophil-specific transcriptomic profiles, we found that CCL4L2 was associated with ICS nonresponse in asthma, which was validated in human and murine lung tissue. CCL4L2 expression was also negatively correlated with pulmonary function change after ICS treatment. In summary, steroids fail to suppress neutrophilic airway inflammation, highlighting the potential need to use alternative therapies such as leukotriene receptor antagonists or DNase I that target the neutrophil-associated phenotype. Furthermore, these results highlight CCL4L2 as a potential therapeutic target for individuals with asthma refractory to ICS.

Background: Asthmatic patients' responses to inhaled corticosteroids (ICS) are variable and difficult to quantify. We have previously defined a Cross-sectional Asthma STEroid Response (CASTER) measure of ICS response. MicroRNAs (miRNAs) have shown strong effects on asthma and inflammatory processes., Objective: The purpose of this study was to identify key associations between circulating miRNAs and ICS response in childhood asthma., Methods: Small RNA sequencing in peripheral blood serum from 580 children with asthma on ICS treatment from The Genetics of Asthma in Costa Rica Study (GACRS) was used to identify miRNAs associated with ICS response using generalized linear models. Replication was conducted in children on ICS from the Childhood Asthma Management Program (CAMP) cohort. The association between replicated miRNAs and the transcriptome of lymphoblastoid cell lines in response to a glucocorticoid was assessed., Results: The association study on the GACRS cohort identified 36 miRNAs associated with ICS response at 10% false discovery rate (FDR), three of which (miR-28-5p, miR-339-3p, and miR-432-5p) were in the same direction of effect and significant in the CAMP replication cohort. In addition, in vitro steroid response lymphoblastoid gene expression analysis revealed 22 dexamethasone responsive genes were significantly associated with three replicated miRNAs. Furthermore, Weighted Gene Co-expression Network Analysis (WGCNA) revealed a significant association between miR-339-3p and two modules (black and magenta) of genes associated with immune response and inflammation pathways., Conclusion: This study highlighted significant association between circulating miRNAs miR-28-5p, miR-339-3p, and miR-432-5p and ICS response. miR-339-3p may be involved in immune dysregulation, which leads to a poor response to ICS treatment.

Intrauterine smoke (IUS) exposure during early childhood has been associated with a number of negative health consequences, including reduced lung function and asthma susceptibility. The biological mechanisms underlying these associations have not been established. MicroRNAs regulate the expression of numerous genes involved in lung development. Thus, investigation of the impact of IUS on miRNA expression during human lung development may elucidate the impact of IUS on post-natal respiratory outcomes. We sought to investigate the effect of IUS exposure on miRNA expression during early lung development. We hypothesized that miRNA-mRNA networks are dysregulated by IUS during human lung development and that these miRNAs may be associated with future risk of asthma and allergy. Human fetal lung samples from a prenatal tissue retrieval program were tested for differential miRNA expression with IUS exposure (measured using placental cotinine concentration). RNA was extracted and miRNA-sequencing was performed. We performed differential expression using IUS exposure, with covariate adjustment. We also considered the above model with an additional sex-by-IUS interaction term, allowing IUS effects to differ by male and female samples. Using paired gene expression profiles, we created sex-stratified miRNA-mRNA correlation networks predictive of IUS using DIABLO. We additionally evaluated whether miRNAs were associated with asthma and allergy outcomes in a cohort of childhood asthma. We profiled pseudoglandular lung miRNA in n = 298 samples, 139 (47%) of which had evidence of IUS exposure. Of 515 miRNAs, 25 were significantly associated with intrauterine smoke exposure (q-value < 0.10). The IUS associated miRNAs were correlated with well-known asthma genes (e.g., ORM1-Like Protein 3, ORDML3 ) and enriched in disease-relevant pathways (oxidative stress). Eleven IUS-miRNAs were also correlated with clinical measures (e.g., Immunoglobulin E andlungfunction) in children with asthma, further supporting their likely disease relevance. Lastly, we found substantial differences in IUS effects by sex, finding 95 significant IUS-miRNAs in male samples, but only four miRNAs in female samples. The miRNA-mRNA correlation networks were predictive of IUS (AUC = 0.78 in males and 0.86 in females) and suggested that IUS-miRNAs are involved in regulation of disease-relevant genes (e.g., A disintegrin and metalloproteinase domain 19 ( ADAM19) , LBH regulator of WNT signaling ( LBH) ) and sex hormone signaling (Coactivator associated methyltransferase 1( CARM1) ). Our study demonstrated differential expression of miRNAs by IUS during early prenatal human lung development, which may be modified by sex. Based on their gene targets and correlation to clinical asthma and atopy outcomes, these IUS-miRNAs may be relevant for subsequent allergy and asthma risk. Our study provides insight into the impact of IUS in human fetal lung transcriptional networks and on the developmental origins of asthma and allergic disorders.

Purpose of Review: The study of microRNA in asthma has revealed a vibrant new level of gene regulation underlying asthma pathology. Several miRNAs have been shown to be important in asthma, influencing various biological mechanisms which lead to asthma pathology and symptoms. In addition, miRNAs have been proposed as biomarkers of asthma affection status, asthma severity, and asthma treatment response. We review all recent asthma-miRNA work, while also presenting comprehensive tables of all miRNA results related to asthma., Recent Findings: We here reviewed 63 recent studies published reporting asthma and miRNA research, and an additional 14 reviews of the same. We summarized the information for both adult and childhood asthma, as well as research on miRNAs in asthma-COPD overlap syndrome (ACOs), and virus-induced asthma exacerbations. We attempted to present a comprehensive collection of recently published asthma-associated miRNAs as well as tables of all published asthma-related miRNA results., (© 2022. The Author(s).)

MiRNAs have been shown to play a powerful regulatory role in the progression of serious diseases, including cancer, Alzheimer's, and others, raising the possibility of new miRNA-based therapies for these conditions. Current experimental methods, such as differential expression analysis, can discover disease-associated miRNAs, yet many of these miRNAs play no functional role in disease progression. Interventional experiments used to discover disease causal miRNAs can be time consuming and costly. We present DisiMiR: a novel computational method that predicts pathogenic miRNAs by inferring biological characteristics of pathogenicity, including network influence and evolutionary conservation. DisiMiR separates disease causal miRNAs from merely disease-associated miRNAs, and was accurate in four diseases: breast cancer (0.826 AUC), Alzheimer's (0.794 AUC), gastric cancer (0.853 AUC), and hepatocellular cancer (0.957 AUC). Additionally, DisiMiR can generate hypotheses effectively: 78.4% of its false positives that are mentioned in the literature have been confirmed to be causal through recently published research. In this work, we show that DisiMiR is a powerful tool that can be used to efficiently and flexibly to predict pathogenic miRNAs in an expression dataset, for the further elucidation of disease mechanisms, and the potential identification of novel drug targets.

Background: Interstitial lung abnormalities (ILA) are radiologic findings that may progress to idiopathic pulmonary fibrosis (IPF). Blood gene expression profiles can predict IPF mortality, but whether these same genes associate with ILA and ILA outcomes is unknown. This study evaluated if a previously described blood gene expression profile associated with IPF mortality is associated with ILA and all-cause mortality., Methods: In COPDGene and ECLIPSE study participants with visual scoring of ILA and gene expression data, we evaluated the association of a previously described IPF mortality score with ILA and mortality. We also trained a new ILA score, derived using genes from the IPF score, in a subset of COPDGene. We tested the association with ILA and mortality on the remainder of COPDGene and ECLIPSE., Results: In 1469 COPDGene (training n = 734; testing n = 735) and 571 ECLIPSE participants, the IPF score was not associated with ILA or mortality. However, an ILA score derived from IPF score genes was associated with ILA (meta-analysis of test datasets OR 1.4 [95% CI: 1.2-1.6]) and mortality (HR 1.25 [95% CI: 1.12-1.41]). Six of the 11 genes in the ILA score had discordant directions of effects compared to the IPF score. The ILA score partially mediated the effects of age on mortality (11.8% proportion mediated)., Conclusions: An ILA gene expression score, derived from IPF mortality-associated genes, identified genes with concordant and discordant effects on IPF mortality and ILA. These results suggest shared, and unique biologic processes, amongst those with ILA, IPF, aging, and death., (© 2022. The Author(s).)

Introduction: Asthma is a complex disease with heterogeneous expression/severity. There is growing interest in defining asthma endotypes consistently associated with different responses to therapy, focusing on type 2 inflammation (Th2) as a key pathological mechanism. Current asthma endotypes are defined primarily by clinical/laboratory criteria. Each endotype is likely characterised by distinct molecular mechanisms that identify optimal therapies., Methods: We applied unsupervised (without a priori clinical criteria) principal component analysis on sputum airway cells RNA-sequencing transcriptomic data from 19 asthmatics from the Severe Asthma Research Program at baseline and 6-8 weeks follow-up after a 40 mg dose of intramuscular corticosteroids. We investigated principal components PC1, PC3 for association with 55 clinical variables., Results: PC3 was associated with baseline Th2 clinical features including blood (rank-sum p=0.0082) and airway (rank-sum p=0.0024) eosinophilia, FEV 1 change (Kendall tau-b R=-0.333 (-0.592 to -0.012)) and follow-up FEV 1 albuterol response (Kendall tau-b R=0.392 (0.079 to 0.634)). PC1 with blood basophlia (rank-sum p=0.0191). The top 5% genes contributing to PC1, PC3 were enriched for distinct immune system/inflammation ontologies suggesting distinct subject-specific clusters of transcriptomic response to corticosteroids. PC3 association with FEV 1 change was reproduced in silico in a comparable independent 14-subject (baseline, 8 weeks after daily inhaled corticosteroids (ICS)) airway epithelial cells microRNAome dataset., Conclusions: Transcriptomic PCs from this unsupervised methodology define molecular pharmacogenomic endotypes that may yield novel biology underlying different subject-specific responses to corticosteroid therapy in asthma, and optimal personalised asthma care. Top contributing genes to these PCs may suggest new therapeutic targets., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2022. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.)

MicroRNAs have been independently associated with asthma and COPD; however, it is unclear if microRNA associations will overlap when evaluating retrospective acute exacerbations. Objective: We hypothesized that peripheral blood microRNAs would be associated with retrospective acute asthma exacerbations in a pediatric asthma cohort and that such associations may also be relevant to acute COPD exacerbations. Methods: We conducted small-RNA sequencing on 374 whole-blood samples from children with asthma ages 6-14 years who participated in the Genetics of Asthma in Costa Rica Study (GACRS) and 450 current and former adult smokers with and without COPD who participated in the COPDGene study. Measurements and Main Results: After QC, we had 351 samples and 649 microRNAs for Differential Expression (DE) analysis between the frequent ( n = 183) and no or infrequent exacerbation ( n = 168) groups in GACRS. Fifteen upregulated miRs had odds ratios (OR) between 1.22 and 1.59 for a doubling of miR counts, while five downregulated miRs had ORs between 0.57 and 0.8. These were assessed for generalization in COPDGene, where three of the upregulated miRs (miR-532-3p, miR-296-5p, and miR-766-3p) and two of the downregulated miRs (miR-7-5p and miR-451b) replicated. Pathway enrichment analysis showed MAPK and PI3K-Akt signaling pathways were strongly enriched for target genes of DE miRNAs and miRNAs generalizing to COPD exacerbations, as well as infection response pathways to various pathogens. Conclusion: miRs (451b; 7-5p; 532-3p; 296-5p and 766-3p) associated with both childhood asthma and adult COPD exacerbations may play a vital role in airflow obstruction and exacerbations and point to shared genomic regulatory machinery underlying exacerbations in both diseases.

Rationale: Current guidelines do not sufficiently capture the heterogeneous nature of asthma; a more detailed molecular classification is needed. Metabolomics represents a novel and compelling approach to derive asthma endotypes (i.e., subtypes defined by functional and/or pathobiological mechanisms). Objectives: To validate metabolomic-driven endotypes of asthma and explore their underlying biology. Methods: In the Genetics of Asthma in Costa Rica Study (GACRS), untargeted metabolomic profiling, similarity network fusion, and spectral clustering was used to identify metabo-endotypes of asthma, and differences in asthma-relevant phenotypes across these metabo-endotypes were explored. The metabo-endotypes were recapitulated in the Childhood Asthma Management Program (CAMP), and clinical differences were determined. Metabolomic drivers of metabo-endotype membership were investigated by meta-analyzing findings from GACRS and CAMP. Measurements and Main Results: Five metabo-endotypes were identified in GACRS with significant differences in asthma-relevant phenotypes, including prebronchodilator (p-ANOVA = 8.3 × 10 -5 ) and postbronchodilator (p-ANOVA = 1.8 × 10 -5 ) FEV 1 /FVC. These differences were validated in the recapitulated metabo-endotypes in CAMP. Cholesterol esters, trigylcerides, and fatty acids were among the most important drivers of metabo-endotype membership. The findings suggest dysregulation of pulmonary surfactant homeostasis may play a role in asthma severity. Conclusions: Clinically meaningful endotypes may be derived and validated using metabolomic data. Interrogating the drivers of these metabo-endotypes has the potential to help understand their pathophysiology.

Background: Inhaled corticosteroids (ICS) are a cornerstone of asthma treatment. However, their efficacy is characterized by wide variability in individual responses., Objective: We investigated the association between genetic variants and risk of exacerbations in adults with asthma and how this association is affected by ICS treatment., Methods: We investigated the pharmacogenetic effect of 10 single nucleotide polymorphisms (SNPs) selected from the literature, including SNPs previously associated with response to ICS (assessed by change in lung function or exacerbations) and novel asthma risk alleles involved in inflammatory pathways, within all adults with asthma from the Dutch population-based Rotterdam study with replication in the American GERA cohort. The interaction effects of the SNPs with ICS on the incidence of asthma exacerbations were assessed using hurdle models adjusting for age, sex, BMI, smoking and treatment step according to the GINA guidelines. Haplotype analyses were also conducted for the SNPs located on the same chromosome., Results: rs242941 (CRHR1) homozygotes for the minor allele (A) showed a significant, replicated increased risk for frequent exacerbations (RR = 6.11, P < 0.005). In contrast, rs1134481 T allele within TBXT (chromosome 6, member of a family associated with embryonic lung development) showed better response with ICS. rs37973 G allele (GLCCI1) showed a significantly poorer response on ICS within the discovery cohort, which was also significant but in the opposite direction in the replication cohort., Conclusion: rs242941 in CRHR1 was associated with poor ICS response. Conversely, TBXT variants were associated with improved ICS response. These associations may reveal specific endotypes, potentially allowing prediction of exacerbation risk and ICS response., (© 2021 The Authors. Clinical & Experimental Allergy published by John Wiley & Sons Ltd.)

Metabolomic indicators of asthma treatment responses have yet to be identified. In this study, we aimed to uncover plasma metabolomic profiles associated with asthma exacerbations while on inhaled corticosteroid (ICS) treatment. We determined whether these profiles change with age from adolescence to adulthood. We utilized data from 170 individuals with asthma on ICS from the Mass General Brigham Biobank to identify plasma metabolites associated with asthma exacerbations while on ICS and examined potential effect modification of metabolite-exacerbation associations by age. We used liquid chromatography-high-resolution mass spectrometry-based metabolomic profiling. Sex-stratified analyses were also performed for the significant associations. The age range of the participating individuals was 13-43 years with a mean age of 33.5 years. Of the 783 endogenous metabolites tested, eight demonstrated significant associations with exacerbation after correction for multiple comparisons and adjusting for potential confounders (Bonferroni p value < 6.2 × 10 -4 ). Potential effect modification by sex was detected for fatty acid metabolites, with males showing a greater reduction in their metabolite levels with ICS exacerbation. Thirty-eight metabolites showed suggestive interactions with age on exacerbation (nominal p -value < 0.05). Our findings demonstrate that plasma metabolomic profiles differ for individuals who experience asthma exacerbations while on ICS. The differentiating metabolites may serve as biomarkers of ICS response and may highlight metabolic pathways underlying ICS response variability.

Objective: MicroRNAs (miRs) are small non-coding RNA molecules of around 18-22 nucleotides that are key regulators of many biologic processes, particularly inflammation. The purpose of this study was to determine the association of circulating miRs from asthmatic children with seasonal variation in allergic inflammation and asthma symptoms., Methods: We used available small RNA sequencing on blood serum from 398 children with mild-to-moderate asthma from the Childhood Asthma Management Program. We used seasonal asthma symptom data at the study baseline and allergen affection status from baseline skin prick tests as primary outcomes. We identified differentially expressed (DE) miRs between pairs of seasons using DESeq2. Regression analysis was used to identify associations between allergy status to specific seasonal allergens and DE miRs in 4 seasons and between seasonal asthma symptom data and DE miRs. We performed pathway enrichment analysis for target genes of the DE miRs using DAVID., Results: After quality control, 398 samples underwent differential analysis between the 4 seasons. We found 52 unique miRs from a total of 81 DE miRs across seasons. Further investigation of the association between these miRs and sensitization to seasonal allergens using skin prick tests revealed that 26 unique miRs from a total of 38 miRs were significantly associated with a same-season allergen. Comparison between seasonal asthma symptom data revealed that 2 of these 26 miRs also had significant associations with asthma symptoms in the same seasons: miR-328-3p ( P < 0.03) and let-7d-3p ( P < 0.05). Enrichment analysis showed that the most enriched pathway clusters were Rap1, Ras, and MAPK signaling pathways., Conclusion: Our results show seasonal variation in miR-328-3p and let-7d-3p are significantly associated with seasonal asthma symptoms and seasonal allergies. These indicate a potentially protective role for let-7d-3p and a deleterious role for miR-328-3p in asthmatics sensitized to mulberry. Further work will determine whether these miRs are drivers or results of the allergic response., Competing Interests: There are no financial or other issue that might lead to conflict of interest., (Copyright © 2021 The Korean Academy of Asthma, Allergy and Clinical Immunology · The Korean Academy of Pediatric Allergy and Respiratory Disease.)

Background: Many microRNAs (miRNAs) have been associated with asthma and chronic obstructive pulmonary disease (COPD). Longitudinal lung function growth trajectories of children with asthma-normal growth, reduced growth (RG), early decline (ED), and RG with an ED (RGED)-have been observed, with RG and RGED associated with adverse outcomes, including COPD., Objective: Our aim was to determine whether circulating miRNAs from an early age in children with asthma would be prognostic of reduced lung function growth patterns over the next 16 years., Methods: We performed small RNA sequencing on sera from 492 children aged 5 to 12 years with mild-to-moderate asthma from the CAMP clinical trial, who were subsequently followed for 12 to 16 years. miRNAs were assessed for differential expression between previously assigned lung function growth patterns., Results: We had 448 samples and 259 miRNAs for differential analysis. In a comparison of the normal and the most severe group (ie, normal growth compared with RGED), we found 1 strongly dysregulated miRNA, hsa-miR-145-5p (P < 8.01E-05). This miR was downregulated in both ED groups (ie, ED and RGED). We verified that miR-145-5p was strongly associated with airway smooth muscle cell growth in vitro., Conclusion: Our results showed that miR-145-5p is associated with the ED patterns of lung function growth leading to COPD in children with asthma and additionally increases airway smooth muscle cell proliferation. This represents a significant extension of our understanding of the role of miR-145-5p in COPD and suggests that reduced expression of miR-145-5p is a risk factor for ED of long-term lung function., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Genome-wide association studies (GWAS) of response to asthma medications have primarily focused on Caucasian populations, with findings that may not be generalizable to minority populations. We derived a polygenic risk score (PRS) for response to albuterol as measured by bronchodilator response (BDR), and examined the PRS in a cohort of Hispanic school-aged children with asthma. We leveraged a published GWAS of BDR to identify relevant genetic variants, and ranked the top variants according to their Combined Annotation Dependent Depletion (CADD) scores. Variants with CADD scores greater than 10 were used to compute the PRS. Once we derived the PRS, we determined the association of the PRS with BDR in a cohort of Hispanic children with asthma (the Genetics of Asthma in Costa Rica Study (GACRS)) in adjusted linear regression models. Mean BDR in GACRS participants was5.6% with a standard deviation of 10.2%. We observed a 0.63% decrease in BDR in response to albuterol for a standard deviation increase in the PRS ( p = 0.05). We also observed decreased odds of a BDR response at or above the 12% threshold for a one standard deviation increase in the PRS (OR = 0.80 (95% CI 0.67 to 0.95)). Our findings show that combining variants from a pharmacogenetic GWAS into a PRS may be useful for predicting medication response in asthma.

Of children with recurrent wheezing in early childhood, approximately half go on to develop asthma. MicroRNAs have been described as excellent non-invasive biomarkers due to their prognostic utility. We hypothesized that circulating microRNAs can predict incident asthma and that that prediction might be modified by vitamin D. We selected 75 participants with recurrent wheezing at 3 years old from the Vitamin D Antenatal Asthma Reduction Trial (VDAART). Plasma samples were collected at age 3 and sequenced for small RNA-Seq. The read counts were normalized and filtered by depth and coverage. Logistic regression was employed to associate miRNAs at age 3 with asthma status at age 5. While the overall effect of miRNA on asthma occurrence was weak, we identified 38 miRNAs with a significant interaction effect with vitamin D and 32 miRNAs with a significant main effect in the high vitamin D treatment group in VDAART. We validated the VDAART results in Project Viva for both the main effect and interaction effect. Meta-analysis was performed on both cohorts to obtain the combined effect and a logistic regression model was used to predict incident asthma at age 7 in Project Viva. Of the 23 overlapped miRNAs in the stratified and interaction analysis above, 9 miRNAs were replicated in Project Viva with strong effect size and remained in the meta-analysis of the two populations. The target genes of the 9 miRNAs were enriched for asthma-related Kyoto Encyclopedia of Genes and Genomes (KEGG) signaling pathways. Using logistic regression, microRNA hsa-miR-574-5p had a good prognostic ability for incident asthma prognosis with an area under the receiver operating characteristic (AUROC) of 0.83. In conclusion, miRNAs appear to be good biomarkers of incident asthma, but only when vitamin D level is considered.

Background: COPD is a leading cause of mortality., Research Question: We hypothesized that applying machine learning to clinical and quantitative CT imaging features would improve mortality prediction in COPD., Study Design and Methods: We selected 30 clinical, spirometric, and imaging features as inputs for a random survival forest. We used top features in a Cox regression to create a machine learning mortality prediction (MLMP) in COPD model and also assessed the performance of other statistical and machine learning models. We trained the models in subjects with moderate to severe COPD from a subset of subjects in Genetic Epidemiology of COPD (COPDGene) and tested prediction performance in the remainder of individuals with moderate to severe COPD in COPDGene and Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints (ECLIPSE). We compared our model with the BMI, airflow obstruction, dyspnea, exercise capacity (BODE) index; BODE modifications; and the age, dyspnea, and airflow obstruction index., Results: We included 2,632 participants from COPDGene and 1,268 participants from ECLIPSE. The top predictors of mortality were 6-min walk distance, FEV 1 % predicted, and age. The top imaging predictor was pulmonary artery-to-aorta ratio. The MLMP-COPD model resulted in a C index ≥ 0.7 in both COPDGene and ECLIPSE (6.4- and 7.2-year median follow-ups, respectively), significantly better than all tested mortality indexes (P < .05). The MLMP-COPD model had fewer predictors but similar performance to that of other models. The group with the highest BODE scores (7-10) had 64% mortality, whereas the highest mortality group defined by the MLMP-COPD model had 77% mortality (P = .012)., Interpretation: An MLMP-COPD model outperformed four existing models for predicting all-cause mortality across two COPD cohorts. Performance of machine learning was similar to that of traditional statistical methods. The model is available online at: https://cdnm.shinyapps.io/cgmortalityapp/., (Copyright © 2020 American College of Chest Physicians. Published by Elsevier Inc. All rights reserved.)

Asthma patient response to inhaled corticosteroids (ICS) is variable and difficult to quantify. We aimed to define a measure of steroid response suitable for pharmacogenetic research in longitudinal and cross-sectional cohorts. Using longitudinal data from the Childhood Asthma Management Program (CAMP) asthma cohort, we defined the Cross-sectional Asthma STEroid Response (CASTER) measure in cross-sectional data. We then applied this to cross-sectional slices of four independent asthma cohorts: The Improving Asthma Control Trial (IMPACT), the Salmeterol or Corticosteroids Study (SOCS), the Pediatric Asthma Controller Trial (PACT), and the Genetics of Asthma in Costa Rica Study (GACRS). CASTER achieved high accuracy on the childhood asthma cohorts: GACRS, PACT, and also on cross-sectional data from CAMP (AUCs 82%, 71%, 63%, respectively). This demonstrates that select cross-sectional clinical information is sufficient to identify good and poor responders to ICS treatment in childhood asthma. Thus, CASTER represents a major improvement in the usability and applicability of steroid response measures in asthma research.

Rationale: Inhaled corticosteroids (ICS) are key treatments for controlling asthma and preventing asthma attacks. However, the responsiveness to ICS varies among individuals. MicroRNAs (miRNAs) have been lauded for their prognostic utility. Objectives: We hypothesized that circulating miRNAs obtained at baseline/prerandomization in the Childhood Asthma Management Program (CAMP) could serve as biomarkers and biologic mediators of ICS clinical response over the 4-year clinical trial period. Methods: We selected baseline serum samples from 462 CAMP subjects subsequently randomized to either ICS (budesonide) or placebo. Samples underwent small RNA sequencing, and read counts were normalized and filtered by depth and coverage. Linear regression was used to associate miRNAs with change in FEV 1 % (prebronchodilator FEV 1 as a percent predicted) over the 4-year treatment period in both main effects and interaction models. We validated the function of the top associated miRNAs by luciferase reporter assays of glucocorticoid-mediated transrepression and predicted response to ICS through logistic regression models. Measurements and Main Results: We identified 7 miRNAs significantly associated with FEV 1 % change ( P  ≤ 0.05) and 15 miRNAs with significant interaction ( P  ≤ 0.05) to ICS versus placebo treatments. We selected three miRNAs for functional validation, of which hsa-miR-155-5p and hsa-miR-532-5p were significantly associated with changes in dexamethasone-induced transrepression of NF-κB. Combined, these two miRNAs were predictive of ICS response over the course of the clinical trial, with an area under the receiver operating characteristic curve of 0.86. Conclusions: We identified two functional circulating miRNAs predictive of asthma ICS treatment response over time.

Background: Genetic factors influence chronic obstructive pulmonary disease (COPD) risk, but the individual variants that have been identified have small effects. We hypothesised that a polygenic risk score using additional variants would predict COPD and associated phenotypes., Methods: We constructed a polygenic risk score using a genome-wide association study of lung function (FEV 1 and FEV 1 /forced vital capacity [FVC]) from the UK Biobank and SpiroMeta. We tested this polygenic risk score in nine cohorts of multiple ethnicities for an association with moderate-to-severe COPD (defined as FEV 1 /FVC <0·7 and FEV 1 <80% of predicted). Associations were tested using logistic regression models, adjusting for age, sex, height, smoking pack-years, and principal components of genetic ancestry. We assessed predictive performance of models by area under the curve. In a subset of studies, we also studied quantitative and qualitative CT imaging phenotypes that reflect parenchymal and airway pathology, and patterns of reduced lung growth., Findings: The polygenic risk score was associated with COPD in European (odds ratio [OR] per SD 1·81 [95% CI 1·74-1·88] and non-European (1·42 [1·34-1·51]) populations. Compared with the first decile, the tenth decile of the polygenic risk score was associated with COPD, with an OR of 7·99 (6·56-9·72) in European ancestry and 4·83 (3·45-6·77) in non-European ancestry cohorts. The polygenic risk score was superior to previously described genetic risk scores and, when combined with clinical risk factors (ie, age, sex, and smoking pack-years), showed improved prediction for COPD compared with a model comprising clinical risk factors alone (AUC 0·80 [0·79-0·81] vs 0·76 [0·75-0·76]). The polygenic risk score was associated with CT imaging phenotypes, including wall area percent, quantitative and qualitative measures of emphysema, local histogram emphysema patterns, and destructive emphysema subtypes. The polygenic risk score was associated with a reduced lung growth pattern., Interpretation: A risk score comprised of genetic variants can identify a small subset of individuals at markedly increased risk for moderate-to-severe COPD, emphysema subtypes associated with cigarette smoking, and patterns of reduced lung growth., Funding: US National Institutes of Health, Wellcome Trust., (Copyright © 2020 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY 4.0 license. Published by Elsevier Ltd.. All rights reserved.)

Asthma affects more than 300 million people in the world, costs over $80 billion annually in the United States, and is efficaciously treated with inhaled corticosteroids (ICS). To our knowledge, no studies have examined the real-world effectiveness of ICS, including the combination therapy consisting of ICS and long-acting beta agonists (LABAs), and patterns of use over a 15-year time period. We used data from the Kaiser Permanente Northern California multi-ethnic Genetic Epidemiology Research on Adult Health and Aging (GERA) Cohort which comprises longitudinal electronic health record data of over 100,000 people. Data included longitudinal asthma-related events, such as ambulatory office visits, hospitalizations, emergency department (ED) visits, and fills of ICS and ICS-LABA combination. Asthma exacerbations were defined as an asthma-related ED visit, hospitalization, or oral corticosteroid (OCS) burst. We used an expected-value approach to determine ICS and ICS-LABA coverage over exacerbation events. We compared rates of exacerbation of subjects on ICS or ICS-LABAs to their own rates of exacerbation when off controller medications. We found ICS-LABA therapy had significant effects, reducing all types of exacerbations per day by a factor of 1.76 (95% CI (1.06, 2.93), p = 0.03) and, specifically, bursts per day by a factor of 1.91 (95% CI (1.04, 3.53), p = 0.037). In conclusion, ICS-LABA therapy was significantly associated with fewer asthma-related exacerbations in a large population of individuals with asthma who were followed for 15 years.

Background: Asthma is known to display different phenotypes across the life-course, suggesting that age related changes are particularly relevant to understanding asthma pathogenesis and remission. We have previously demonstrated that a lung function phenotype associated with asthma, bronchodilator response, is reduced with age, at rate of 0.24 percent per year. Methods: In this study, we interrogated the serum metabolome, to determine whether circulating metabolites mediate age-related changes in bronchodilator response (BDR) for individuals with asthma. We used data on 295 participants from the follow-up phase of the CAMP clinical trial (age 12.2-25.9 years; mean BDR of 8%, standard deviation 7%). Using a counterfactual framework, we analyzed over 500 pareto-scaled metabolites using mediation analysis to identify indirect effects of age through potential metabolite mediators. Results: There was a significant indirect effect of age on BDR through 4 plasmalogens (C36:1 PC and related metabolites) (Indirect Effect Beta = -0.001, p = 0.006). Conclusions: Our findings suggest that plasmalogens may contribute to age-related asthma phenotypes, and may also serve as potential pharmacologic targets for enhancement of lung function in individuals with asthma. Trial Registration: This work uses data from the previous clinical trial of asthma, the Childhood Asthma Management Program (CAMP), registered at ClinicalTrials.gov, # NCT00000575., (Copyright © 2020 Sordillo, Lutz, Kelly, McGeachie, Dahlin, Tantisira, Clish, Lasky-Su and Wu.)

Background: Global gene expression levels are known to be highly dependent upon gross demographic features including age, yet identification of age-related genomic indicators has yet to be comprehensively undertaken in a disease and treatment-specific context., Methods: We used gene expression data from CD4+ lymphocytes in the Asthma BioRepository for Integrative Genomic Exploration (Asthma BRIDGE), an open-access collection of subjects participating in genetic studies of asthma with available gene expression data. Replication population participants were Puerto Rico islanders recruited as part of the ongoing Genes environments & Admixture in Latino Americans (GALA II), who provided nasal brushings for transcript sequencing. The main outcome measure was chronic asthma control as derived by questionnaires. Genomic associations were performed using regression of chronic asthma control score on gene expression with age in years as a covariate, including a multiplicative interaction term for gene expression times age., Results: The SMARCD1 gene (SWI/SNF-related matrix-associated actin-dependent regulator of chromatin subfamily D member 1) interacted with age to influence chronic asthma control on inhaled corticosteroids, with a doubling of expression leading to an increase of 1.3 units of chronic asthma control per year (95% CI [0.86, 1.74], p = 6 × 10 - 9 ), suggesting worsening asthma control with increasing age. This result replicated in GALA II (p = 3.8 × 10 - 8 ). Cellular assays confirmed the role of SMARCD1 in glucocorticoid response in airway epithelial cells., Conclusion: Focusing on age-dependent factors may help identify novel indicators of asthma medication response. Age appears to modulate the effect of SMARCD1 on asthma control with inhaled corticosteroids.

The role of metabolism in modifying age-related differential responses to asthma medications is insufficiently understood. The objective of this study was to determine the role of the metabolome in modifying the effect of age on bronchodilator response (BDR) in individuals with asthma. We used longitudinal measures of BDR and plasma metabolomic profiling in 565 children with asthma from the Childhood Asthma Management Program (CAMP) to identify age by metabolite interactions on BDR. The mean ages at the three studied time-points across 16 years of follow-up in CAMP were 8.8, 12.8, and 16.8 years; the mean BDRs were 11%, 9% and 8%, respectively. Of 501 identified metabolites, 39 (7.8%) demonstrated a significant interaction with age on BDR ( p -value < 0.05). We were able to validate two significant interactions in 320 children with asthma from the Genetics of Asthma in Costa Rica Study; 2-hydroxyglutarate, a compound involved in butanoate metabolism (interaction; CAMP: β = -0.004, p = 1.8 × 10 -4 ; GACRS: β = -0.015, p = 0.018 ), and a cholesterol ester; CE C18:1 (CAMP: β = 0.005, p = 0.006; GACRS: β = 0.023, p = 0.041 ) Five additional metabolites had a p -value < 0.1 in GACRS, including Gammaminobutyric acid (GABA), C16:0 CE, C20:4 CE, C18.0 CE and ribothymidine. These findings suggest Cholesterol esters and GABA may modify the estimated effect of age on bronchodilator response.

Background: Cell-free plasma mitochondrial DNA (mtDNA) levels are associated with endothelial dysfunction and differential outcomes in critical illness. A substantial alteration in metabolic homeostasis is commonly observed in severe critical illness. We hypothesized that metabolic profiles significantly differ between critically ill patients relative to their level of plasma mtDNA., Methods: We performed a metabolomic study with biorepository plasma samples collected from 73 adults with systemic inflammatory response syndrome or sepsis at a single academic medical center. Patients were treated in a 20-bed medical ICU between 2008 and 2010. To identify key metabolites and metabolic pathways related to plasma NADH dehydrogenase 1 (ND1) mtDNA levels in critical illness, we first generated metabolomic data using gas and liquid chromatography-mass spectroscopy. We performed fold change analysis and volcano plot visualization based on false discovery rate-adjusted p values to evaluate the distribution of individual metabolite concentrations relative to ND1 mtDNA levels. We followed this by performing orthogonal partial least squares discriminant analysis to identify individual metabolites that discriminated ND1 mtDNA groups. We then interrogated the entire metabolomic profile using pathway overrepresentation analysis to identify groups of metabolite pathways that were different relative to ND1 mtDNA levels., Results: Metabolomic profiles significantly differed in critically ill patients with ND1 mtDNA levels ≥ 3200 copies/μl plasma relative to those with an ND1 mtDNA level < 3200 copies/μl plasma. Several analytical strategies showed that patients with ND1 mtDNA levels ≥ 3200 copies/μl plasma had significant decreases in glycerophosphocholines and increases in short-chain acylcarnitines., Conclusions: Differential metabolic profiles during critical illness are associated with cell-free plasma ND1 mtDNA levels that are indicative of cell damage. Elevated plasma ND1 mtDNA levels are associated with decreases in glycerophosphocholines and increases in short-chain acylcarnitines that reflect phospholipid metabolism dysregulation and decreased mitochondrial function, respectively.

Background: Variation in response to the most commonly used class of asthma controller medication, inhaled corticosteroids, presents a serious challenge in asthma management, particularly for steroid-resistant patients with little or no response to treatment., Objective: We applied a systems biology approach to primary clinical and genomic data to identify and validate genes that modulate steroid response in asthmatic children., Methods: We selected 104 inhaled corticosteroid-treated asthmatic non-Hispanic white children and determined a steroid responsiveness endophenotype (SRE) using observations of 6 clinical measures over 4 years. We modeled each subject's cellular steroid response using data from a previously published study of immortalized lymphoblastoid cell lines under dexamethasone (DEX) and sham treatment. We integrated SRE with immortalized lymphoblastoid cell line DEX responses and genotypes to build a genome-scale network using the Reverse Engineering, Forward Simulation modeling framework, identifying 7 genes modulating SRE., Results: Three of these genes were functionally validated by using a stable nuclear factor κ-light-chain-enhancer of activated B cells luciferase reporter in A549 human lung epithelial cells, IL-1β cytokine stimulation, and DEX treatment. By using small interfering RNA transfection, knockdown of family with sequence similarity 129 member A (FAM129A) produced a reduction in steroid treatment response (P < .001)., Conclusion: With this systems-based approach, we have shown that FAM129A is associated with variation in clinical asthma steroid responsiveness and that FAM129A modulates steroid responsiveness in lung epithelial cells., (Copyright © 2018 American Academy of Allergy, Asthma & Immunology. Published by Elsevier Inc. All rights reserved.)

To explore novel methods for the analysis of metabolomics data, we compared the ability of Partial Least Squares Discriminant Analysis (PLS-DA) and Bayesian networks (BN) to build predictive plasma metabolite models of age three asthma status in 411 three year olds ( n = 59 cases and 352 controls) from the Vitamin D Antenatal Asthma Reduction Trial (VDAART) study. The standard PLS-DA approach had impressive accuracy for the prediction of age three asthma with an Area Under the Curve Convex Hull (AUCCH) of 81%. However, a permutation test indicated the possibility of overfitting. In contrast, a predictive Bayesian network including 42 metabolites had a significantly higher AUCCH of 92.1% ( p for difference < 0.001), with no evidence that this accuracy was due to overfitting. Both models provided biologically informative insights into asthma; in particular, a role for dysregulated arginine metabolism and several exogenous metabolites that deserve further investigation as potential causative agents. As the BN model outperformed the PLS-DA model in both accuracy and decreased risk of overfitting, it may therefore represent a viable alternative to typical analytical approaches for the investigation of metabolomics data.

Background: Single omic analyses have provided some insight into the basis of lung function in children with asthma, but the underlying biologic pathways are still poorly understood., Methods: Weighted gene coexpression network analysis (WGCNA) was used to identify modules of coregulated gene transcripts and metabolites in blood among 325 children with asthma from the Genetic Epidemiology of Asthma in Costa Rica study. The biology of modules associated with lung function as measured by FEV 1 , the FEV 1 /FVC ratio, bronchodilator response, and airway responsiveness to methacholine was explored. Significantly correlated gene-metabolite module pairs were then identified, and their constituent features were analyzed for biologic pathway enrichments., Results: WGCNA clustered 25,060 gene probes and 8,185 metabolite features into eight gene modules and eight metabolite modules, where four and six, respectively, were associated with lung function (P ≤ .05). The gene modules were enriched for immune, mitotic, and metabolic processes and asthma-associated microRNA targets. The metabolite modules were enriched for lipid and amino acid metabolism. Integration of correlated gene-metabolite modules expanded the single omic findings, linking the FEV 1 /FVC ratio with ORMDL3 and dysregulated lipid metabolism. This finding was replicated in an independent population., Conclusions: The results of this hypothesis-generating study suggest a mechanistic basis for multiple asthma genes, including ORMDL3, and a role for lipid metabolism. They demonstrate that integrating multiple omic technologies may provide a more informative picture of asthmatic lung function biology than single omic analyses., (Copyright © 2018 American College of Chest Physicians. Published by Elsevier Inc. All rights reserved.)

It is well known that depression is associated with asthma symptoms. We assessed the combined effects of genetic factors and depression on asthma symptom severity using Bayesian network (BN) analysis. The common 100 top-ranked single-nucleotide polymorphisms (SNPs) were obtained from two genome-wide association studies of symptom severity in two childhood asthmatics trials (CAMP (Childhood Asthma Management Program) and CARE (Childhood Asthma Research and Education)). Using SNPs plus five discretized variables (depression, anxiety, age, sex, and race), we performed BN analysis in 529 CAMP subjects. We identified two nodes (depression and rs4672619 mapping to ERBB4 (Erb-B2 receptor tyrosine kinase 4)) that were within the Markov neighborhood of the symptom node in the network and then evaluated the interactive effects of depressive status and rs4672619 genotypes on asthma symptom severity. In childhood asthmatics with homozygous reference alleles, severe depression was related to less severe symptoms. However, in childhood asthmatics with heterozygous alleles and homozygous variant alleles, depression and symptom severity showed a positive correlation (interaction permutation P value = 0.019). We then tried to evaluate whether the interactive effects that we found were sustained in another independent cohort of elderly asthmatics. Contrary to the findings from childhood asthmatics, elderly asthmatics with homozygous reference alleles showed a positive correlation between depression and symptom severity, and elderly asthmatics with heterozygous alleles and homozygous variant alleles showed a negative correlation (interaction permutation P value = 0.003). In conclusion, we have identified a novel SNP, rs4672619, that shows interactive effects with depression on asthma symptom severity in childhood and elderly asthmatics in opposite directions.

Background: Circulating microRNAs have shown promise as non-invasive biomarkers and predictors of disease activity. Prior asthma studies using clinical, biochemical and genomic data have not shown excellent prediction of exacerbation. We hypothesized that a panel of circulating microRNAs in a pediatric asthma cohort combined with an exacerbation clinical score might predict exacerbation better than the latter alone., Methods: Serum samples from 153 children at randomization in the Childhood Asthma Management Program were profiled for 754 microRNAs. Data dichotomized for asthma exacerbation one year after randomization to inhaled corticosteroid treatment were used for binary logistic regression with miRNA expressions and exacerbation clinical score., Results: 12 of 125 well-detected circulating microRNAs had significant odd ratios for exacerbation with miR-206 being most significant. Each doubling of expression of the 12 microRNA corresponded to a 25-67% increase in exacerbation risk. Stepwise logistic regression yielded a 3-microRNA model (miR-146b, miR-206 and miR-720) that, combined with the exacerbation clinical score, had excellent predictive power with a 0.81 AUROC. These 3 microRNAs were involved in NF-kβ and GSK3/AKT pathways., Conclusions: This combined circulating microRNA-clinical score model predicted exacerbation in asthmatic subjects on inhaled corticosteroids better than each constituent feature alone., Trial Registration: ClinicalTrials.gov Identifier: NCT00000575 .

Background: Metabolic homeostasis is substantially disrupted in critical illness. Given the pleiotropic effects of vitamin D, we hypothesized that metabolic profiles differ between critically ill patients relative to their vitamin D status., Methods: We performed a metabolomics study on biorepository samples collected from a single academic medical center on 65 adults with systemic inflammatory response syndrome or sepsis treated in a 20-bed medical ICU between 2008 and 2010. To identify key metabolites and metabolic pathways related to vitamin D status in critical illness, we first generated metabolomic data using gas and liquid chromatography mass spectroscopy. We followed this by partial least squares-discriminant analysis to identify individual metabolites that were significant. We then interrogated the entire metabolomics profile using metabolite set enrichment analysis to identify groups of metabolites and pathways that were differentiates of vitamin D status. Finally we performed logistic regression to construct a network model of chemical-protein target interactions important in vitamin D status., Results: Metabolomic profiles significantly differed in critically ill patients with 25(OH)D ≤ 15 ng/ml relative to those with levels >15 ng/ml. In particular, increased 1,5-anhydroglucitol, tryptophan betaine, and 3-hydroxyoctanoate as well as decreased 2-arachidonoyl-glycerophosphocholine and N-6-trimethyllysine were strong predictors of 25(OH)D >15 ng/ml. The combination of these five metabolites led to an area under the curve for discrimination for 25(OH)D > 15 ng/ml of 0.82 (95% CI 0.71-0.93). The metabolite pathways related to glutathione metabolism and glutamate metabolism are significantly enriched with regard to vitamin D status., Conclusion: Vitamin D status is associated with differential metabolic profiles during critical illness. Glutathione and glutamate pathway metabolism, which play principal roles in redox regulation and immunomodulation, respectively, were significantly altered with vitamin D status.

The heritability of chronic obstructive pulmonary disease (COPD) cannot be fully explained by recognized genetic risk factors identified as achieving genome-wide significance. In addition, the combined contribution of genetic variation to COPD risk has not been fully explored. We sought to determine: (1) whether studies of variants from previous studies of COPD or lung function in a larger sample could identify additional associated variants, particularly for severe COPD; and (2) the impact of genetic risk scores on COPD. We genotyped 3,346 single-nucleotide polymorphisms (SNPs) in 2,588 cases (1,803 severe COPD) and 1,782 control subjects from four cohorts, and performed association testing with COPD, combining these results with existing genotyping data from 6,633 cases (3,497 severe COPD) and 5,704 control subjects. In addition, we developed genetic risk scores from SNPs associated with lung function and COPD and tested their discriminatory power for COPD-related measures. We identified significant associations between SNPs near PPIC (P = 1.28 × 10 -8 ) and PPP4R4/SERPINA1 (P = 1.01 × 10 -8 ) and severe COPD; the latter association may be driven by recognized variants in SERPINA1. Genetic risk scores based on SNPs previously associated with COPD and lung function had a modest ability to discriminate COPD (area under the curve, ∼0.6), and accounted for a mean 0.9-1.9% lower forced expiratory volume in 1 second percent predicted for each additional risk allele. In a large genetic association analysis, we identified associations with severe COPD near PPIC and SERPINA1. A risk score based on combining genetic variants had modest, but significant, effects on risk of COPD and lung function.

Purpose of Review: This review will catalog the many recent longitudinal studies that have investigated the relationship between asthma and lung function, or the persistence and trajectories of lung function deficits., Recent Findings: Recent work has reported on 50-year follow-ups of some prominent population cohorts. A history of asthma confers a 10-30-fold risk of chronic obstructive pulmonary disease. Individuals reaching a reduced maximum growth of forced expiratory volume in 1 s in early adulthood are at risk for early or more severe chronic obstructive pulmonary disease (COPD)., Summary: Taken together, there is a wealth of overlapping cohort studies of lung function, asthma and COPD. These show that asthma is associated with reduced lung function, which may start in infancy or prenatally, persists through childhood and adulthood and predisposes for early or more severe COPD.

Background: The gut microbiome in infancy influences immune system maturation, and may have an important impact on allergic disease risk., Objective: We sought to determine how prenatal and early life factors impact the gut microbiome in a relatively large, ethnically diverse study population of infants at age 3 to 6 months, who were enrolled in Vitamin D Antenatal Asthma Reduction Trial, a clinical trial of vitamin D supplementation in pregnancy to prevent asthma and allergies in offspring., Methods: We performed 16S rRNA gene sequencing on 333 infants' stool samples. Microbial diversity was computed using the Shannon index. Factor analysis applied to the top 25 most abundant taxa revealed 4 underlying bacterial coabundance groups; the first dominated by Firmicutes (Lachnospiraceae/Clostridiales), the second by Proteobacteria (Klebsiella/Enterobacter), the third by Bacteriodetes, and the fourth by Veillonella. Scores for coabundance groups were used as outcomes in regression models, with prenatal/birth and demographic characteristics as independent predictors. Multivariate analysis, using all microbial community members, was also conducted., Results: White race/ethnicity was associated with lower diversity but higher Bacteroidetes coabundance scores. C-section birth was associated with higher diversity, but decreased Bacteroidetes coabundance scores. Firmicutes scores were higher for infants born by C-section. Breast-fed infants had lower proportions of Clostridiales. Cord blood vitamin D was linked to increased Lachnobacterium, but decreased Lactococcus., Conclusions: The findings presented here suggest that race, mode of delivery, breast-feeding, and cord blood vitamin D levels are associated with infant gut microbiome composition, with possible long-term implications for immune system modulation and asthma/allergic disease incidence., (Copyright © 2016 American Academy of Allergy, Asthma & Immunology. Published by Elsevier Inc. All rights reserved.)

Asthma is a complex disease well-suited to metabolomic profiling, both for the development of novel biomarkers and for the improved understanding of pathophysiology. In this review, we summarize the 21 existing metabolomic studies of asthma in humans, all of which reported significant findings and concluded that individual metabolites and metabolomic profiles measured in exhaled breath condensate, urine, plasma, and serum could identify people with asthma and asthma phenotypes with high discriminatory ability. There was considerable consistency across the studies in terms of the reported biomarkers, regardless of biospecimen, profiling technology, and population age. In particular, acetate, adenosine, alanine, hippurate, succinate, threonine, and trans-aconitate, and pathways relating to hypoxia response, oxidative stress, immunity, inflammation, lipid metabolism and the tricarboxylic acid cycle were all identified as significant in at least two studies. There were also a number of nonreplicated results; however, the literature is not yet sufficiently developed to determine whether these represent spurious findings or reflect the substantial heterogeneity and limited statistical power in the studies and their methods to date. This review highlights the need for additional asthma metabolomic studies to explore these issues, and, further, the need for standardized methods in the way these studies are conducted. We conclude by discussing the potential of translation of these metabolomic findings into clinically useful biomarkers and the crucial role that integrated omics is likely to play in this endeavor., (Copyright © 2016 American College of Chest Physicians. Published by Elsevier Inc. All rights reserved.)

Introduction: Preeclampsia is a leading cause of maternal and fetal mortality worldwide, yet its exact pathogenesis remains elusive., Objectives: This study, nested within the Vitamin D Antenatal Asthma Reduction Trial (VDAART), aimed to develop integrated omics models of preeclampsia that have utility in both prediction and in the elucidation of underlying biological mechanisms., Methods: Metabolomic profiling was performed on first trimester plasma samples of 47 pregnant women from VDAART who subsequently developed preeclampsia and 62 controls with healthy pregnancies, using liquid-chromatography tandem mass-spectrometry. Metabolomic profiles were generated based on logistic regression models and assessed using Received Operator Characteristic Curve analysis. These profiles were compared to profiles from generated using third trimester samples. The first trimester metabolite profile was then integrated with a pre-existing transcriptomic profile using network methods., Results: In total, 72 (0.9%) metabolite features were associated (p<0.01) with preeclampsia after adjustment for maternal age, race, and gestational age. These features had moderate to good discriminatory ability; in ROC curve analyses a summary score based on these features displayed an area under the curve (AUC) of 0.794 (95%CI 0.700, 0.888). This profile retained the ability to distinguish preeclamptic from healthy pregnancies in the third trimester (AUC:0.762 (95% CI 0.663, 0.860)). Additionally, metabolite set enrichment analysis identified common pathways, including glycerophospholipid metabolism, at the two time-points. Integration with the transcriptomic signature refined these results suggesting a particular role for lipid imbalance, immune function and the circulatory system., Conclusions: These findings suggest it is possible to develop a predictive metabolomic profile of preeclampsia. This profile is characterized by changes in lipid and amino acid metabolism and dysregulation of immune response and can be refined through interaction with transcriptomic data. However validation in larger and more diverse populations is required., Competing Interests: Conflict of Interest: All authors declare that they have no conflicts of interest.

Rationale: Patterns of longitudinal lung function growth and decline in childhood asthma have been shown to be important in determining risk for future respiratory ailments including chronic airway obstruction and chronic obstructive pulmonary disease., Objectives: To determine the genetic underpinnings of lung function patterns in subjects with childhood asthma., Methods: We performed a genome-wide association study of 581 non-Hispanic white individuals with asthma that were previously classified by patterns of lung function growth and decline (normal growth, normal growth with early decline, reduced growth, and reduced growth with early decline). The strongest association was also measured in two additional cohorts: a small asthma cohort and a large chronic obstructive pulmonary disease metaanalysis cohort. Interaction between the genomic region encompassing the most strongly associated single-nucleotide polymorphism and nearby genes was assessed by two chromosome conformation capture assays., Measurements and Main Results: An intergenic single-nucleotide polymorphism (rs4445257) on chromosome 8 was strongly associated with the normal growth with early decline pattern compared with all other pattern groups (P = 6.7 × 10 -9 ; odds ratio, 2.8; 95% confidence interval, 2.0-4.0); replication analysis suggested this variant had opposite effects in normal growth with early decline and reduced growth with early decline pattern groups. Chromosome conformation capture experiments indicated a chromatin interaction between rs4445257 and the promoter of the distal CSMD3 gene., Conclusions: Early decline in lung function after normal growth is associated with a genetic polymorphism that may also protect against early decline in reduced growth groups. Clinical trial registered with www.clinicaltrials.gov (NCT00000575).

Introduction: While whole genome prediction (WGP) methods have recently demonstrated successes in the prediction of complex genetic diseases, they have not yet been applied to asthma and related phenotypes. Longitudinal patterns of lung function differ between asthmatics, but these phenotypes have not been assessed for heritability or predictive ability. Herein, we assess the heritability and genetic predictability of asthma-related phenotypes., Methods: We applied several WGP methods to a well-phenotyped cohort of 832 children with mild-to-moderate asthma from CAMP. We assessed narrow-sense heritability and predictability for airway hyperresponsiveness, serum immunoglobulin E, blood eosinophil count, pre- and post-bronchodilator forced expiratory volume in 1 sec (FEV 1 ), bronchodilator response, steroid responsiveness, and longitudinal patterns of lung function (normal growth, reduced growth, early decline, and their combinations). Prediction accuracy was evaluated using a training/testing set split of the cohort., Results: We found that longitudinal lung function phenotypes demonstrated significant narrow-sense heritability (reduced growth, 95%; normal growth with early decline, 55%). These same phenotypes also showed significant polygenic prediction (areas under the curve [AUCs] 56% to 62%). Including additional demographic covariates in the models increased prediction 4-8%, with reduced growth increasing from 62% to 66% AUC. We found that prediction with a genomic relatedness matrix was improved by filtering available SNPs based on chromatin evidence, and this result extended across cohorts., Conclusions: Longitudinal reduced lung function growth displayed extremely high heritability. All phenotypes with significant heritability showed significant polygenic prediction. Using SNP-prioritization increased prediction across cohorts. WGP methods show promise in predicting asthma-related heritable traits.