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.

piRNA associates with immune diseases.

PIWI-interacting RNA (piRNA) is the most abundant small non-coding RNA in animal cells, typically 26-31 nucleotides in length and it binds with PIWI proteins, a subfamily of Argonaute proteins. Initially discovered in germ cells, piRNA is well known for its role in silencing transposons and maintaining genome integrity. However, piRNA is also present in somatic cells as well as in extracellular vesicles and exosomes. While piRNA has been extensively studied in various diseases, particular cancer, its function in immune diseases remains unclear. In this review, we summarize current research on piRNA in immune diseases. We first introduce the basic characteristics, biogenesis and functions of piRNA. Then, we review the association of piRNA with different types of immune diseases, including autoimmune diseases, immunodeficiency diseases, infectious diseases, and other immune-related diseases. piRNA is considered a promising biomarker for diseases, highlighting the need for further research into its potential mechanisms in disease pathogenesis.

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.

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.

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

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.

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.

Correction: Pharmacogenomic associations of adverse drug reactions in asthma: systematic review and research prioritization.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

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.

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.

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

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.

Gene expression changes in lymphoblastoid cell lines and primary B cells by dexamethasone.

BACKGROUND: Human Epstein-Barr virus-transformed lymphoblastoid cell lines (LCLs) have been thought to be a useful model system for pharmacogenomics studies. The purpose of this study was to determine the effect of Epstein-Barr virus transformation on gene expression changes by dexamethasone (Dex) in LCLs and primary B cells (PBCs) derived from the same individuals. PATIENTS AND METHODS: We prepared LCLs and purified PBCs from the same six male donors participating in the Childhood Asthma Management Program clinical trial, and compared mRNA profiles after 6 h incubation with Dex (10 mol/l) or sham buffer. We assessed differential expression and put the list of differentially expressed genes into the web interface of ConsensusPathDB to find the pathway-level interpretation of our genes specified. As a supplementary analysis, we looked at the expression of the Dex-regulated (inducing or repressing) genes in treatment-naive PBCs and LCLs (pre-Dex treatment) from the GSE30916 dataset. RESULTS: By hierarchical clustering, we found clustering of probes by cell types but not by individuals irrespective of Dex treatment. We observed that the Dex-regulated genes significantly overlapped in PBCs and LCLs. In addition, the expression of these genes showed significant correlations between treatment-naive PBCs and LCLs. Common genes showing significantly decreased expressions by the Dex treatment in both cells were enriched in immune responses and proinflammatory signaling pathways. CONCLUSION: Taken together, these results suggest the uses of LCLs are representative of the primary biologic effects of corticosteroids treatment.

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.

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.