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.

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 MicroRNAs associated with Bronchodilator Response in Childhood Asthma

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

Neutrophil extracellular trap production and CCL4L2 expression influence corticosteroid response in asthma

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.

Association of body anthropometry and obstructive sleep apnea in children: Variations observed in Hispanic children

Objectives

Obesity is a risk factor for obstructive sleep apnea (OSA) in children. Childhood obesity rates vary amongst different ethnic groups. Here the interaction of Hispanic ethnicity and obesity on OSA risk was evaluated.

Methods

Retrospective cross-sectional analysis of consecutive children undergoing polysomnography and anthropometry using bioelectrical impedance from 2017 to 2020. Demographics obtained from the medical chart. Children who had also undergone cardiometabolic testing were identified and the relationship of cardiometabolic markers with OSA and anthropometry was assessed.

Results

Data from 1217 children revealed Hispanic children were more likely to have moderate-severe OSA (36.0%) compared to Non-Hispanic children (26.5%), p < 0.001. Hispanic children had greater Body mass index (BMI), BMI percentile and percent body fat, p < 0.0001. In children that underwent cardiometabolic testing, Hispanic children had significantly greater serum alanine aminotransferase levels (ALT) levels. Following adjustment of age and sex, Hispanic ethnicity was not found to moderate the association of anthropometry with OSA, anthropometry with cardiometabolic markers, and OSA with cardiometabolic markers.

Conclusions

OSA was more likely in Hispanic children; this relationship was likely driven by obesity status rather than ethnicity. Among children undergoing cardiometabolic testing, Hispanic children were observed to have greater ALT concentrations however ethnicity did not impact the association of anthropometry and ALT or other cardiometabolic markers.

Systems Genomics Reveals microRNA Regulation of ICS Response in Childhood Asthma

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.

Pharmacogenetics of childhood uncontrolled asthma

INTRODUCTION

Asthma is a heterogeneous, multifactorial disease with multiple genetic and environmental risk factors playing a role in pathogenesis and therapeutic response. Understanding of pharmacogenetics can help with matching individualized treatments to specific genotypes of asthma to improve therapeutic outcomes especially in uncontrolled or severe asthma.

AREAS COVERED

In this review, we outline novel information about biology, pathways, and mechanisms related to interindividual variability in drug response (corticosteroids, bronchodilators, leukotriene modifiers, and biologics) for childhood asthma. We discuss candidate gene, genome-wide association studies and newer omics studies including epigenomics, transcriptomics, proteomics, and metabolomics as well as integrative genomics and systems biology methods related to childhood asthma. The articles were obtained after a series of searches, last updated November 2022, using database PubMed/CINAHL DB.

EXPERT OPINION

Implementation of pharmacogenetic algorithms can improve therapeutic targeting in children with asthma, particularly with severe or uncontrolled asthma who typically have challenges in clinical management and carry considerable financial burden. Future studies focusing on potential biomarkers both clinical and pharmacogenetic can help formulate a prognostic test for asthma treatment response that would represent true bench to bedside clinical implementation.

Researching COVID to enhance recovery (RECOVER) pediatric study protocol: Rationale, objectives and design

Importance

The prevalence, pathophysiology, and long-term outcomes of COVID-19 (post-acute sequelae of SARS-CoV-2 [PASC] or "Long COVID") in children and young adults remain unknown. Studies must address the urgent need to define PASC, its mechanisms, and potential treatment targets in children and young adults.

Observations

We describe the protocol for the Pediatric Observational Cohort Study of the NIH's RE searching COV ID to E nhance R ecovery (RECOVER) Initiative. RECOVER-Pediatrics is an observational meta-cohort study of caregiver-child pairs (birth through 17 years) and young adults (18 through 25 years), recruited from more than 100 sites across the US. This report focuses on two of five cohorts that comprise RECOVER-Pediatrics: 1) a de novo RECOVER prospective cohort of children and young adults with and without previous or current infection; and 2) an extant cohort derived from the Adolescent Brain Cognitive Development (ABCD) study ( n =10,000). The de novo cohort incorporates three tiers of data collection: 1) remote baseline assessments (Tier 1, n=6000); 2) longitudinal follow-up for up to 4 years (Tier 2, n=6000); and 3) a subset of participants, primarily the most severely affected by PASC, who will undergo deep phenotyping to explore PASC pathophysiology (Tier 3, n=600). Youth enrolled in the ABCD study participate in Tier 1. The pediatric protocol was developed as a collaborative partnership of investigators, patients, researchers, clinicians, community partners, and federal partners, intentionally promoting inclusivity and diversity. The protocol is adaptive to facilitate responses to emerging science.

Conclusions and Relevance

RECOVER-Pediatrics seeks to characterize the clinical course, underlying mechanisms, and long-term effects of PASC from birth through 25 years old. RECOVER-Pediatrics is designed to elucidate the epidemiology, four-year clinical course, and sociodemographic correlates of pediatric PASC. The data and biosamples will allow examination of mechanistic hypotheses and biomarkers, thus providing insights into potential therapeutic interventions.

Clinical Trialsgov Identifier

Clinical Trial Registration: http://www.clinicaltrials.gov . Unique identifier: NCT05172011.

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.

Antagonizing cholecystokinin A receptor in the lung attenuates obesity-induced airway hyperresponsiveness

Obesity increases asthma prevalence and severity. However, the underlying mechanisms are poorly understood, and consequently, therapeutic options for asthma patients with obesity remain limited. Here we report that cholecystokinin-a metabolic hormone best known for its role in signaling satiation and fat metabolism-is increased in the lungs of obese mice and that pharmacological blockade of cholecystokinin A receptor signaling reduces obesity-associated airway hyperresponsiveness. Activation of cholecystokinin A receptor by the hormone induces contraction of airway smooth muscle cells. In vivo, cholecystokinin level is elevated in the lungs of both genetically and diet-induced obese mice. Importantly, intranasal administration of cholecystokinin A receptor antagonists (proglumide and devazepide) suppresses the airway hyperresponsiveness in the obese mice. Together, our results reveal an unexpected role for cholecystokinin in the lung and support the repurposing of cholecystokinin A receptor antagonists as a potential therapy for asthma patients with obesity.

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.

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

Blood miRNAs Are Linked to Frequent Asthma Exacerbations in Childhood Asthma and Adult COPD

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.

Identification of ROBO2 as a Potential Locus Associated with Inhaled Corticosteroid Response in Childhood Asthma

Inhaled corticosteroids (ICS) are the most common asthma controller medication. An important contribution of genetic factors in ICS response has been evidenced. Here, we aimed to identify novel genetic markers involved in ICS response in asthma. A genome-wide association study (GWAS) of the change in lung function after 6 weeks of ICS treatment was performed in 166 asthma patients from the SLOVENIA study. Patients with an improvement in lung function ≥8% were considered as ICS responders. Suggestively associated variants (p-value ≤ 5 × 10⁻⁶) were evaluated in an independent study (n = 175). Validation of the association with asthma exacerbations despite ICS use was attempted in European (n = 2681) and admixed (n = 1347) populations. Variants previously associated with ICS response were also assessed for replication. As a result, the SNP rs1166980 from the ROBO2 gene was suggestively associated with the change in lung function (OR for G allele: 7.01, 95% CI: 3.29–14.93, p = 4.61 × 10⁻⁷), although this was not validated in CAMP. ROBO2 showed gene-level evidence of replication with asthma exacerbations despite ICS use in Europeans (minimum p-value = 1.44 × 10⁻⁵), but not in admixed individuals. The association of PDE10A-T with ICS response described by a previous study was validated. This study suggests that ROBO2 could be a potential novel locus for ICS response in Europeans.

Seasonal Variation in miR-328-3p and let-7d-3p Are Associated With Seasonal Allergies and Asthma Symptoms in Children

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.

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.

Genome-wide association study of asthma exacerbations despite inhaled corticosteroid use

RATIONALE

Substantial variability in response to asthma treatment with inhaled corticosteroids (ICS) has been described among individuals and populations, suggesting the contribution of genetic factors. Nonetheless, only a few genes have been identified to date. We aimed to identify genetic variants associated with asthma exacerbations despite ICS use in European children and young adults and to validate the findings in non-Europeans. Moreover, we explored whether a gene-set enrichment analysis could suggest potential novel asthma therapies.

METHODS

A genome-wide association study (GWAS) of asthma exacerbations was tested in 2681 children of European descent treated with ICS from eight studies. Suggestive association signals were followed up for replication in 538 European asthma patients. Further evaluation was performed in 1773 non-Europeans. Variants revealed by published GWAS were assessed for replication. Additionally, gene-set enrichment analysis focused on drugs was performed.

RESULTS

10 independent variants were associated with asthma exacerbations despite ICS treatment in the discovery phase (p≤5×10-6). Of those, one variant at the CACNA2D3-WNT5A locus was nominally replicated in Europeans (rs67026078; p=0.010), but this was not validated in non-European populations. Five other genes associated with ICS response in previous studies were replicated. Additionally, an enrichment of associations in genes regulated by trichostatin A treatment was found.

CONCLUSIONS

The intergenic region of CACNA2D3 and WNT5A was revealed as a novel locus for asthma exacerbations despite ICS treatment in European populations. Genes associated were related to trichostatin A, suggesting that this drug could regulate the molecular mechanisms involved in treatment response.

Circulating MicroRNA: Incident Asthma Prediction and Vitamin D Effect Modification

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.

Enhancing the prediction of childhood asthma remission: Integrating clinical factors with microRNAs

The novel integration of baseline clinical and microRNA variables significantly improves the long-term individualized prediction of childhood asthma remission by early adulthood compared to using clinical variables alone.

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.

Co-methylation analysis in lung tissue identifies pathways for fetal origins of COPD

COPD likely has developmental origins; however, the underlying molecular mechanisms are not fully identified. Investigation of lung tissue-specific epigenetic modifications such as DNA methylation using network approaches might facilitate insights linking in utero smoke (IUS) exposure and risk for COPD in adulthood.We performed genome-wide methylation profiling for adult lung DNA from 160 surgical samples and 78 fetal lung DNA samples isolated from discarded tissue at 8-18 weeks of gestation. Co-methylation networks were constructed to identify preserved modules that shared methylation patterns in fetal and adult lung tissues and associations with fetal IUS exposure, gestational age and COPD.Weighted correlation networks highlighted preserved and co-methylated modules for both fetal and adult lung data associated with fetal IUS exposure, COPD and lower adult lung function. These modules were significantly enriched for genes involved in embryonic organ development and specific inflammation-related pathways, including Hippo, phosphatidylinositol 3-kinase/protein kinase B (PI3K/AKT), Wnt, mitogen-activated protein kinase and transforming growth factor-β signalling. Gestational age-associated modules were remarkably preserved for COPD and lung function, and were also annotated to genes enriched for the Wnt and PI3K/AKT pathways.Epigenetic network perturbations in fetal lung tissue exposed to IUS and of early lung development recapitulated in adult lung tissue from ex-smokers with COPD. Overlapping fetal and adult lung tissue network modules highlighted putative disease pathways supportive of exposure-related and age-associated developmental origins of COPD.

Pharmacogenomic associations of adverse drug reactions in asthma: systematic review and research prioritisation

A systematic review of pharmacogenomic studies capturing adverse drug reactions (ADRs) related to asthma medications was undertaken, and a survey of Pharmacogenomics in Childhood Asthma (PiCA) consortia members was conducted. Studies were eligible if genetic polymorphisms were compared with suspected ADR(s) in a patient with asthma, as either a primary or secondary outcome. Five studies met the inclusion criteria. The ADRs and polymorphisms identified were change in lung function tests (rs1042713), adrenal suppression (rs591118), and decreased bone mineral density (rs6461639) and accretion (rs9896933, rs2074439). Two of these polymorphisms were replicated within the paper, but none had external replication. Priorities from PiCA consortia members (representing 15 institution in eight countries) for future studies were tachycardia (SABA/LABA), adrenal suppression/crisis and growth suppression (corticosteroids), sleep/behaviour disturbances (leukotriene receptor antagonists), and nausea and vomiting (theophylline). Future pharmacogenomic studies in asthma should collect relevant ADR data as well as markers of efficacy.

CASTER: Cross-Sectional Asthma STEroid Response Measurement

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.

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.

Real-Life Patterns of Exacerbations While on Inhaled Corticosteroids and Long-Acting Beta Agonists for Asthma over 15 Years

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.

COMPSRA: a COMprehensive Platform for Small RNA-Seq data Analysis

Small RNA-Seq is a common means to interrogate the small RNA’ome or the full spectrum of small RNAs (<200 nucleotide length) of a biological system. A pivotal problem in NGS based small RNA analysis is identifying and quantifying the small RNA’ome constituent components. For example, small RNAs in the circulatory system (circulating RNAs) are potential disease biomarkers and their function is being actively investigated. Most existing NGS data analysis tools focus on the microRNA component and a few other small RNA types like piRNA, snRNA and snoRNA. A comprehensive platform is needed to interrogate the full small RNA’ome, a prerequisite for down-stream data analysis. We present COMPSRA, a comprehensive modular stand-alone platform for identifying and quantifying small RNAs from small RNA sequencing data. COMPSRA contains prebuilt customizable standard RNA databases and sequence processing tools to enable turnkey basic small RNA analysis. We evaluated COMPSRA against comparable existing tools on small RNA sequencing data set from serum samples of 12 healthy human controls, and COMPSRA identified a greater diversity and abundance of small RNA molecules. COMPSRA is modular, stand-alone and integrates multiple customizable RNA databases and sequence processing tool and is distributed under the GNU General Public License free to non-commercial registered users at https://github.com/cougarlj/COMPSRA.

DNA methylation is associated with improvement in lung function on inhaled corticosteroids in pediatric asthmatics

Asthma is the most common chronic disease in children. Inhaled corticosteroids (ICS) are the first-line treatment for asthma control, but up to one-third of children have a poor treatment response. The mechanism of ICS resistance is poorly understood, and the role of DNA methylation in ICS treatment response is not known. We examined the association between peripheral blood DNA methylation and ICS treatment response in 152 pediatric persistent asthmatics from the Childhood Asthma Management Program. Response to ICS was measured by the percentage change in forced expiratory volume in 1 s (FEV1) 8 weeks after treatment initiation. The top CpG sites with a nominal P value less than 0.001 were correlated with gene expression using Pearson's and partial correlations. In 152 participants, mean±SD age was 9.8±2.0 years and median change in FEV1 after ICS initiation was 4.6% (interquartile range: 10.4%). A total of 545 CpG sites were differentially methylated (nominal P<0.05), and seven CpG sites had a nominal P value less than 0.001. Relative hypermethylation of cg20434811, cg02822723, cg14066280, cg27254601, and cg23913400 and relative hypomethylation of cg24937126 and cg24711626 were associated with an increase in FEV1 on ICS treatment. One CpG site was associated with gene expression. Relative hypermethylation of cg27254601 was associated with both an increase in FEV1 and BOLA2 expression (ρ=0.25, P=0.02). We identified a novel association between BOLA2 methylation, gene expression, and ICS response as measured by lung function. Pharmacoepigenetics has the potential to detect treatment sensitivity in persistent childhood asthma.

Genome-wide interaction study reveals age-dependent determinants of responsiveness to inhaled corticosteroids in individuals with asthma

While genome-wide association studies have identified genes involved in differential treatment responses to inhaled corticosteroids (ICS) in asthma, few studies have evaluated the potential effects of age in this context. A significant proportion of asthmatics experience exacerbations (hospitalizations and emergency department visits) during ICS treatment. We evaluated the interaction of genetic variation and age on ICS response (measured by the occurrence of exacerbations) through a genome-wide interaction study (GWIS) of 1,321 adult and child asthmatic patients of European ancestry. We identified 107 genome-wide suggestive (P<10-05) age-by-genotype interactions, two of which also met genome-wide significance (P<5x10-08) (rs34631960 [OR 2.3±1.6-3.3] in thrombospondin type 1 domain-containing protein 4 (THSD4) and rs2328386 [OR 0.5±0.3-0.7] in human immunodeficiency virus type I enhancer binding protein 2 (HIVEP2)) by joint analysis of GWIS results from discovery and replication populations. In addition to THSD4 and HIVEP2, age-by-genotype interactions also prioritized genes previously identified as asthma candidate genes, including DPP10, HDAC9, TBXAS1, FBXL7, and GSDMB/ORMDL3, as pharmacogenomic loci as well. This study is the first to link these genes to a pharmacogenetic trait for asthma.

Epigenome-wide meta-analysis of blood DNA methylation in newborns and children identifies numerous loci related to gestational age

BACKGROUND

Preterm birth and shorter duration of pregnancy are associated with increased morbidity in neonatal and later life. As the epigenome is known to have an important role during fetal development, we investigated associations between gestational age and blood DNA methylation in children.

METHODS

We performed meta-analysis of Illumina's HumanMethylation450-array associations between gestational age and cord blood DNA methylation in 3648 newborns from 17 cohorts without common pregnancy complications, induced delivery or caesarean section. We also explored associations of gestational age with DNA methylation measured at 4-18 years in additional pediatric cohorts. Follow-up analyses of DNA methylation and gene expression correlations were performed in cord blood. DNA methylation profiles were also explored in tissues relevant for gestational age health effects: fetal brain and lung.

RESULTS

We identified 8899 CpGs in cord blood that were associated with gestational age (range 27-42 weeks), at Bonferroni significance, P < 1.06 × 10- 7, of which 3343 were novel. These were annotated to 4966 genes. After restricting findings to at least three significant adjacent CpGs, we identified 1276 CpGs annotated to 325 genes. Results were generally consistent when analyses were restricted to term births. Cord blood findings tended not to persist into childhood and adolescence. Pathway analyses identified enrichment for biological processes critical to embryonic development. Follow-up of identified genes showed correlations between gestational age and DNA methylation levels in fetal brain and lung tissue, as well as correlation with expression levels.

CONCLUSIONS

We identified numerous CpGs differentially methylated in relation to gestational age at birth that appear to reflect fetal developmental processes across tissues. These findings may contribute to understanding mechanisms linking gestational age to health effects.

Pulmonary Neuroendocrine Cells Secrete γ-Aminobutyric Acid to Induce Goblet Cell Hyperplasia in Primate Models

Mucus overproduction is a major contributor to morbidity and mortality in asthma. Mucus overproduction is induced by orchestrated actions of multiple factors that include inflammatory cytokines and γ-aminobutyric acid (GABA). GABA is produced only by pulmonary neuroendocrine cells (PNECs) in the mouse lung. Recent studies in a neonatal mouse model of allergic inflammation have shown that PNECs play an essential role in mucus overproduction by GABA hypersecretion. Whether PNECs mediate dysregulated GABA signaling for mucus overproduction in asthma is unknown. In this study, we characterized the cellular source of GABA in the lungs of nonhuman primates and humans and assessed GABA secretion and signaling in primate disease models. We found that like in mice, PNECs were the major source of GABA in primate lungs. In addition, an infant nonhuman primate model of asthma exhibited an increase in GABA secretion. Furthermore, subjects with asthma had elevated levels of expression of a subset of GABA type α (GABAα) and type β (GABAβ) receptors in airway epithelium compared with those of healthy control subjects. Last, employing a normal human bronchial epithelial cell model of preinduced mucus overproduction, we showed pharmaceutical blockade of GABAα and GABAβ receptor signaling reversed the effect of IL-13 on MUC5AC gene expression and goblet cell proliferation. Together, our data demonstrate an evolutionarily conserved intraepithelial GABA signaling that, in concert with IL-13, plays an essential role in mucus overproduction. Our findings may offer new strategies to ameliorate mucus overproduction in patients with asthma by targeting PNEC secretion and GABA signaling.

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.

Genetic Signatures of Acute Asthma Exacerbation Related With Ineffective Response to Corticosteroid

Purpose

Acute exacerbation (AE) is an important domain of asthma management and may be related with ineffective response to corticosteroid. This study aimed to find mechanisms of AE using genome-wide gene expression profiles of blood cells from asthmatics and its perturbation by in vitro dexamethasone (Dex)-treatment.

Methods

We utilized lymphoblastoid B cells from 107 childhood asthmatics and peripheral blood mononuclear cells from 29 adult asthmatics who were treated with inhaled corticosteroids. We searched for a preserved co-expression gene module significantly associated with the AE rate in both cohorts and measured expression changes of genes belong to this module after Dex-treatment.

Results

We identified a preserved module composed of 77 genes. Among them, expressions of 2 genes (EIF2AK2 and NOL11) decreased significantly after Dex-treatment in both cohorts. EIF2AK2, a key gene acting antiviral defense mechanism, showed significantly higher expressions in asthmatics with AE. The protein repair pathway was enriched significantly in 64 genes which belong to the preserved module but showed no expression differences after Dex-treatment in both cohorts. Among them, MSRA and MSRB2 may play key roles by controlling oxidative stress.

Conclusions

Many genes belong to the AE rate-associated and preserved module identified in blood cells from childhood and adults asthmatics showed no expression changes after in vitro Dex-treatment. These findings suggest that we may need alternative treatment options to corticosteroids to prevent AE. EIF2AK2, MSRA and MSRB2 expressions on blood cells may help us select AE-susceptible asthmatics and adjust treatments to prevent AE.

A nuclease-polymerase chain reaction enables amplification of probes used for capture-based DNA target enrichment

DNA target enrichment via hybridization capture is a commonly adopted approach which remains expensive due in-part to using biotinylated-probe panels. Here we provide a novel isothermal amplification reaction to amplify rapidly existing probe panels without knowledge of the sequences involved, thereby decreasing a major portion of the overall sample preparation cost. The reaction employs two thermostable enzymes, BST-polymerase and duplex-specific nuclease DSN. DSN initiates random ‘nicks’ on double-stranded-DNA which enable BST to polymerize DNA by displacing the nicked-strand. Displaced strands re-hybridize and the process leads to an exponential chain-reaction generating biotinylated DNA fragments within minutes. When starting from single-stranded-DNA, DNA is first converted to double-stranded-DNA via terminal-deoxynucleotidyl-transferase (TdT) prior to initiation of BST–DSN reaction. Biotinylated probes generated by TdT–BST–DSN (TBD) reactions using panels of 33, 190 or 7186 DNA targets are used for hybrid-capture-based target enrichment from amplified circulating-DNA, followed by targeted re-sequencing. Polymerase-nuclease isothermal-chain-reactions generate random amplified probes with no apparent sequence dependence. One round of target-capture using TBD probes generates a modest on-target sequencing ratio, while two successive rounds of capture generate >80% on-target reads with good sequencing uniformity. TBD-reactions generate enough capture-probes to increase by approximately two to three orders-of-magnitude the target-enrichment experiments possible from an initial set of probes.

DNA methylation is associated with inhaled corticosteroid response in persistent childhood asthmatics

BACKGROUND

Response to inhaled corticosteroids is highly variable, and the association between DNA methylation and treatment response is not known.

OBJECTIVE

To examine the association between peripheral blood DNA methylation and inhaled corticosteroid response in children with persistent asthma.

METHODS

Epigenome-wide DNA methylation was analysed in individuals on inhaled corticosteroids in three independent and ethnically diverse cohorts-Childhood Asthma Management Program (CAMP); Children, Allergy, Milieu, Stockholm, Epidemiology (BAMSE); and Genetic Epidemiology of Asthma in Costa Rica Study (GACRS). Treatment response was evaluated using two definitions, the absence of emergency department visits and/or hospitalizations and the absence oral corticosteroid use while on inhaled corticosteroid therapy. CpG sites meeting nominal significance (P < 0.05) for each outcome were combined in a three-cohort meta-analysis with adjustment for multiple testing. DNA methylation was correlated with gene expression using Pearson and partial correlations.

RESULTS

In 154 subjects from CAMP, 72 from BAMSE, and 168 from GACRS, relative hypomethylation of cg00066816 (171 bases upstream of IL12B) was associated with the absence of emergency department visits and/or hospitalizations (Q = 0.03) in all cohorts and lower IL12B expression (ρ = 0.34, P = 0.01) in BAMSE. Relative hypermethylation of cg04256470 (688 bases upstream of CORT) was associated with the absence of oral corticosteroid use (Q = 0.04) in all cohorts and higher CORT expression (ρ = 0.20, P = 0.045) in CAMP.

CONCLUSION AND CLINICAL RELEVANCE

Differential DNA methylation of IL12B and CORT are associated with inhaled corticosteroid treatment response in persistent childhood asthmatics. Pharmaco-methylation can identify novel markers of treatment sensitivity in asthma.

Whole genome sequencing identifies high-impact variants in well-known pharmacogenomic genes

More than 1100 genetic loci have been correlated with drug response outcomes but disproportionately few have been translated into clinical practice. One explanation for the low rate of clinical implementation is that the majority of associated variants may be in linkage disequilibrium (LD) with the causal variants, which are often elusive. This study aims to identify and characterize likely causal variants within well-established pharmacogenomic genes using next-generation sequencing data from the 1000 Genomes Project. We identified 69,319 genetic variations within 160 pharmacogenomic genes, of which 8207 variants are in strong LD (r2>0.8) with known pharmacogenomic variants. Of the latter, eight are coding or structural variants predicted to have high impact, with 19 additional missense variants that are predicted to have moderate impact. In conclusion, we identified putatively functional variants within known pharmacogenomics loci that could account for the association signals and represent the missing causative variants underlying drug response phenotypes.

Large-scale, multiethnic genome-wide association study identifies novel loci contributing to asthma susceptibility in adults

In this large, multi-ethnic GWAS of asthma, we identified novel associations with potential functional relevance for asthma susceptibility in older adults of diverse racial backgrounds.

Remission of persistent childhood asthma: Early predictors of adult outcomes

BACKGROUND

Few data exist on the predictors of asthma remission by early adulthood in North America.

OBJECTIVE

The predictors of adult asthma remission were determined in a multiethnic population of patients with mild-to-moderate persistent childhood asthma.

METHODS

Asthma remission in early adulthood was measured by using 2 definitions: a clinical and a strict definition. Both included normal lung function and the absence of symptoms, exacerbations, and medication use. The strict definition also included normal airways responsiveness. Predictors were identified from 23 baseline measures by using multivariate logistic regression. The probability of remission was modeled by using decision tree analysis.

RESULTS

In 879 subjects the mean ± SD baseline age was 8.8 ± 2.1 years, 59.4% were male, and 68.7% were white. By adulthood, 229 (26.0%) of 879 participants were in clinical remission, and 111 (15.0%) of 741 participants were in strict remission. The degree of FEV₁/forced vital capacity (FVC) ratio impairment was the largest predictor of asthma remission. More than half of boys and two thirds of girls with baseline FEV₁/FVC ratios of 90% or greater were in remission at adulthood. Decreased airways responsiveness was also a predictor for both remission definitions (clinical remission odds ratio, 1.23 [95% CI, 1.09-1.39]; strict remission odds ratio, 1.52 [95% CI, 1.26-1.84]). The combination of normal FEV₁/FVC ratio, airways responsiveness, and serum eosinophil count at baseline yielded greater than 80% probability of remission by adulthood.

CONCLUSION

A considerable minority of patients with persistent childhood asthma will have disease remission by adulthood. Clinical prognostic indicators of asthma remission, including baseline lung function, can be seen from an early age.

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 .

Whole-Genome Sequencing of Pharmacogenetic Drug Response in Racially Diverse Children with Asthma

RATIONALE

Albuterol, a bronchodilator medication, is the first-line therapy for asthma worldwide. There are significant racial/ethnic differences in albuterol drug response.

OBJECTIVES

To identify genetic variants important for bronchodilator drug response (BDR) in racially diverse children.

METHODS

We performed the first whole-genome sequencing pharmacogenetics study from 1,441 children with asthma from the tails of the BDR distribution to identify genetic association with BDR.

MEASUREMENTS AND MAIN RESULTS

We identified population-specific and shared genetic variants associated with BDR, including genome-wide significant (P < 3.53 × 10-7) and suggestive (P < 7.06 × 10-6) loci near genes previously associated with lung capacity (DNAH5), immunity (NFKB1 and PLCB1), and β-adrenergic signaling (ADAMTS3 and COX18). Functional analyses of the BDR-associated SNP in NFKB1 revealed potential regulatory function in bronchial smooth muscle cells. The SNP is also an expression quantitative trait locus for a neighboring gene, SLC39A8. The lack of other asthma study populations with BDR and whole-genome sequencing data on minority children makes it impossible to perform replication of our rare variant associations. Minority underrepresentation also poses significant challenges to identify age-matched and population-matched cohorts of sufficient sample size for replication of our common variant findings.

CONCLUSIONS

The lack of minority data, despite a collaboration of eight universities and 13 individual laboratories, highlights the urgent need for a dedicated national effort to prioritize diversity in research. Our study expands the understanding of pharmacogenetic analyses in racially/ethnically diverse populations and advances the foundation for precision medicine in at-risk and understudied minority populations.

The phosphatidylinositide 3-kinase (PI3K) signaling pathway is a determinant of zileuton response in adults with asthma

Variable responsiveness to zileuton, a leukotriene antagonist used to treat asthma, may be due in part to genetic variation. While individual SNPs were previously associated with zileuton-related lung function changes, specific quantitative trait loci (QTLs) and biological pathways that may contribute have not been identified. In this study, we investigated the hypothesis that genetic variation within biological pathways is associated with zileuton response. We performed an integrative QTL mapping and pathway enrichment study to investigate data from a GWAS of zileuton response, in addition to mRNA expression profiles and leukotriene production data from lymphoblastoid cell lines (LCLs) (derived from asthmatics) that were treated with zileuton or ethanol (control). We identified 1060 QTLs jointly associated with zileuton-related differential LTB₄ production in LCLs and lung function change in patients taking zileuton, of which eight QTLs were also significantly associated with persistent LTB₄ production in LCLs following zileuton treatment (i.e., ‘poor’ responders). Four nominally significant trans-eQTLs were predicted to regulate three candidate genes (SELL, MTF2, and GAL), the expression of which was significantly reduced in LCLs following zileuton treatment. Gene and pathway enrichment analyses of QTL associations identified multiple genes and pathways, predominantly related to phosphatidyl inositol signaling via PI3K. We validated the PI3K pathway activation status in a subset of LCLs demonstrating variable zileuton-related LTB₄ production, and show that in contrast to LCLs that responded to zileuton, the PI3K pathway was activated in poor responder LCLs. Collectively, these findings demonstrate a role for the PIK3 pathway and its targets as important determinants of differential responsiveness to zileuton.

Circulating microRNAs and association with methacholine PC20 in the Childhood Asthma Management Program (CAMP) cohort

Introduction

Circulating microRNAs (miRNA) are promising biomarkers for human diseases. Our study hypothesizes that circulating miRNA would reveal candidate biomarkers related to airway hyperresponsiveness (AHR) and provide biologic insights into asthma epigenetic influences.

Methods

Serum samples obtained at randomization for 160 children in the Childhood Asthma Management Program were profiled using a TaqMan miRNA array set. The association of the isolated miRNA with methacholine PC₂₀ was assessed. Network and pathway analyses were performed. Functional validation of two significant miRNAs was performed in human airway smooth muscle cells (HASMs).

Results

Of 155 well-detected circulating miRNAs, eight were significantly associated with PC₂₀ with the strongest association with miR-296-5p. Pathway analysis revealed miR-16-5p as a network hub, and involvement of multiple miRNAs interacting with genes in the FoxO and Hippo signaling pathways by KEGG analysis. Functional validation of two miRNA in HASM showed effects on cell growth and diameter.

Conclusion

Reduced circulatory miRNA expression at baseline is associated with an increase in PC₂₀. These miRNA provide biologic insights into, and may serve as biomarkers of, asthma severity. miR-16-5p and -30d-5p regulate airway smooth muscle phenotypes critically involved in asthma pathogenesis, supporting a mechanistic link to these findings. Functional ASM phenotypes may be directly relevant to AHR.

Differential connectivity of gene regulatory networks distinguishes corticosteroid response in asthma

BACKGROUND

Variations in drug response between individuals have prevented us from achieving high drug efficacy in treating many complex diseases, including asthma. Genetics plays an important role in accounting for such interindividual variations in drug response. However, systematic approaches for addressing how genetic factors and their regulators determine variations in drug response in asthma treatment are lacking.

OBJECTIVE

We sought to identify key transcriptional regulators of corticosteroid response in asthma using a novel systems biology approach.

METHODS

We used Passing Attributes between Networks for Data Assimilations (PANDA) to construct the gene regulatory networks associated with good responders and poor responders to inhaled corticosteroids based on a subset of 145 white children with asthma who participated in the Childhood Asthma Management Cohort. PANDA uses gene expression profiles and published relationships among genes, transcription factors (TFs), and proteins to construct the directed networks of TFs and genes. We assessed the differential connectivity between the gene regulatory network of good responders versus that of poor responders.

RESULTS

When compared with poor responders, the network of good responders has differential connectivity and distinct ontologies (eg, proapoptosis enriched in network of good responders and antiapoptosis enriched in network of poor responders). Many of the key hubs identified in conjunction with clinical response are also cellular response hubs. Functional validation demonstrated abrogation of differences in corticosteroid-treated cell viability following siRNA knockdown of 2 TFs and differential downstream expression between good responders and poor responders.

CONCLUSIONS

We have identified and validated multiple TFs influencing asthma treatment response. Our results show that differential connectivity analysis can provide new insights into the heterogeneity of drug treatment effects.

Genetics and Genomics of Longitudinal Lung Function Patterns in Individuals with Asthma

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

Age, Sexual Dimorphism, and Disease Associations in the Developing Human Fetal Lung Transcriptome

The fetal origins of disease hypothesis suggests that variations in the course of prenatal lung development may affect life-long pulmonary function growth, decline, and pathobiology. Many studies support the existence of differences in the developing lung trajectory in males and females, and sex-specific differences in the prevalence of chronic lung diseases, such as asthma and bronchopulmonary dysplasia. The objectives of this study were to investigate the early developing fetal lung for transcriptomic correlates of postconception age (maturity) and sex, and their associations with chronic lung diseases. We analyzed whole-lung transcriptome profiles of 61 females and 78 males at 54-127 days postconception (dpc) from nonsmoking mothers using unsupervised principal component analysis and supervised linear regression models. We identified dominant transcriptomic correlates for postconception age and sex with corresponding gene sets that were enriched for developing lung structural and functional ontologies. We observed that the transcriptomic sex difference was not a uniform global time shift/lag, rather, lungs of males appear to be more mature than those of females before 96 dpc, and females appear to be more mature than males after 96 dpc. The age correlate gene set was consistently enriched for asthma and bronchopulmonary dysplasia genes, but the sex correlate gene sets were not. Despite sex differences in the developing fetal lung transcriptome, postconception age appears to be more dominant than sex in the effect of early fetal lung developments on disease risk during this early pseudoglandular phase of development.

Rationale and design of the multiethnic Pharmacogenomics in Childhood Asthma consortium

AIM

International collaboration is needed to enable large-scale pharmacogenomics studies in childhood asthma. Here, we describe the design of the Pharmacogenomics in Childhood Asthma (PiCA) consortium.

MATERIALS & METHODS

Investigators of each study participating in PiCA provided data on the study characteristics by answering an online questionnaire.

RESULTS

A total of 21 studies, including 14,227 children/young persons (58% male), from 12 different countries are currently enrolled in the PiCA consortium. Fifty six percent of the patients are Caucasians. In total, 7619 were inhaled corticosteroid users. Among patients from 13 studies with available data on asthma exacerbations, a third reported exacerbations despite inhaled corticosteroid use. In the future pharmacogenomics studies within the consortium, the pharmacogenomics analyses will be performed separately in each center and the results will be meta-analyzed.

CONCLUSION

PiCA is a valuable platform to perform pharmacogenetics studies within a multiethnic pediatric asthma population.