A novel piwi-interacting RNA associates with type 2-high asthma phenotypes.

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.

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.

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

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

Novel genetic variants associated with inhaled corticosteroid treatment response in older adults with asthma.

INTRODUCTION: Older adults have the greatest burden of asthma and poorest outcomes. The pharmacogenetics of inhaled corticosteroid (ICS) treatment response is not well studied in older adults. METHODS: A genome-wide association study of ICS response was performed in asthmatics of European ancestry in Genetic Epidemiology Research on Adult Health and Aging (GERA) by fitting Cox proportional hazards regression models, followed by validation in the Mass General Brigham (MGB) Biobank and Rotterdam Study. ICS response was measured using two definitions in asthmatics on ICS treatment: (1) absence of oral corticosteroid (OCS) bursts using prescription records and (2) absence of asthma-related exacerbations using diagnosis codes. A fixed-effect meta-analysis was performed for each outcome. The validated single-nucleotide polymorphisms (SNPs) were functionally annotated to standard databases. RESULTS: In 5710 subjects in GERA, 676 subjects in MGB Biobank, and 465 subjects in the Rotterdam Study, four novel SNPs on chromosome six near PTCHD4 validated across all cohorts and met genome-wide significance on meta-analysis for the OCS burst outcome. In 4541 subjects in GERA and 505 subjects in MGB Biobank, 152 SNPs with p<5 × 10-5 were validated across these two cohorts for the asthma-related exacerbation outcome. The validated SNPs included methylation and expression quantitative trait loci for CPED1, CRADD and DST for the OCS burst outcome and GM2A, SNW1, CACNA1C, DPH1, and RPS10 for the asthma-related exacerbation outcome. CONCLUSIONS: Multiple novel SNPs associated with ICS response were identified in older adult asthmatics. Several SNPs annotated to genes previously associated with asthma and other airway or allergic diseases, including PTCHD4.

Metabo-Endotypes of Asthma Reveal Differences in Lung Function: Discovery and Validation in Two TOPMed Cohorts.

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

Intrauterine Smoke Exposure, microRNA Expression during Human Lung Development, and Childhood Asthma.

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.

Lung function, airway and peripheral basophils and eosinophils are associated with molecular pharmacogenomic endotypes of steroid response in severe asthma.

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, FEV1 change (Kendall tau-b R=-0.333 (-0.592 to -0.012)) and follow-up FEV1 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 FEV1 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.

Pharmacogenetics of inhaled corticosteroids and exacerbation risk in adults with asthma.

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.

Blood gene expression risk profiles and interstitial lung abnormalities: COPDGene and ECLIPSE cohort studies.

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.

Recent miRNA Research in Asthma.

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.

COPD-associated miR-145-5p is downregulated in early-decline FEV1 trajectories in childhood asthma.

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.

A polygenic risk score for asthma in a large racially diverse population.

BACKGROUND: Polygenic risk scores (PRSs) will have important utility for asthma and other chronic diseases as a tool for predicting disease incidence and subphenotypes. OBJECTIVE: We utilized findings from a large multiancestry GWAS of asthma to compute a PRS for asthma with relevance for racially diverse populations. METHODS: We derived two PRSs for asthma using a standard approach (based on genome-wide significant variants) and a lasso sum regression approach (allowing all genetic variants to potentially contribute). We used data from the racially diverse Kaiser Permanente GERA cohort (68 638 non-Hispanic Whites, 5874 Hispanics, 6870 Asians and 2760 Blacks). Race was self-reported by questionnaire. RESULTS: For the standard PRS, non-Hispanic Whites showed the highest odds ratio for a standard deviation increase in PRS for asthma (OR = 1.16 (95% CI 1.14-1.18)). The standard PRS was also associated with asthma in Hispanic (OR = 1.12 (95% CI 1.05-1.19)) and Asian (OR = 1.10 (95% CI 1.04-1.17)) subjects, with a trend towards increased risk in Blacks (OR = 1.05 (95% CI 0.97-1.15)). We detected an interaction by sex, with men showing a higher risk of asthma with an increase in PRS as compared to women. The lasso sum regression-derived PRS showed stronger associations with asthma in non-Hispanic White subjects (OR = 1.20 (95% CI 1.18-1.23)), Hispanics (OR = 1.17 (95% 1.10-1.26)), Asians (OR = 1.18 (95% CI 1.10-1.27)) and Blacks (OR = 1.10 (95% CI 0.99-1.22)). CONCLUSION: Polygenic risk scores across multiple racial/ethnic groups were associated with increased asthma risk, suggesting that PRSs have potential as a tool for predicting disease development.

Circulating MicroRNAs and Treatment Response in Childhood Asthma.

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 FEV1% (prebronchodilator FEV1 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 FEV1% 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.

Expression of SMARCD1 interacts with age in association with asthma control on inhaled corticosteroid therapy.

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.

Systems biology and in vitro validation identifies family with sequence similarity 129 member A (FAM129A) as an asthma steroid response modulator.

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.

Circulating microRNAs and prediction of asthma exacerbation in childhood asthma.

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 .

Plasma mitochondrial DNA and metabolomic alterations in severe critical illness.

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.

Factors influencing the infant gut microbiome at age 3-6 months: Findings from the ethnically diverse Vitamin D Antenatal Asthma Reduction Trial (VDAART).

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.

Genetic Association and Risk Scores in a Chronic Obstructive Pulmonary Disease Meta-analysis of 16,707 Subjects.

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.

Genetic control of gene expression at novel and established chronic obstructive pulmonary disease loci.

Genetic risk loci have been identified for a wide range of diseases through genome-wide association studies (GWAS), but the relevant functional mechanisms have been identified for only a small proportion of these GWAS-identified loci. By integrating results from the largest current GWAS of chronic obstructive disease (COPD) with expression quantitative trait locus (eQTL) analysis in whole blood and sputum from 121 subjects with COPD from the ECLIPSE Study, this analysis identifies loci that are simultaneously associated with COPD and the expression of nearby genes (COPD eQTLs). After integrative analysis, 19 COPD eQTLs were identified, including all four previously identified genome-wide significant loci near HHIP, FAM13A, and the 15q25 and 19q13 loci. For each COPD eQTL, fine mapping and colocalization analysis to identify causal shared eQTL and GWAS variants identified a subset of sites with moderate-to-strong evidence of harboring at least one shared variant responsible for both the eQTL and GWAS signals. Transcription factor binding site (TFBS) analysis confirms that multiple COPD eQTL lead SNPs disrupt TFBS, and enhancer enrichment analysis for loci with the strongest colocalization signals showed enrichment for blood-related cell types (CD3 and CD4+ T cells, lymphoblastoid cell lines). In summary, integrative eQTL and GWAS analysis confirms that genetic control of gene expression plays a key role in the genetic architecture of COPD and identifies specific blood-related cell types as likely participants in the functional pathway from GWAS-associated variant to disease phenotype.