Airway Epithelial Gene Expression Research Articles

Dysfunctional neutrophil type 1 interferon responses in preschool children with recurrent wheezing and IL-4–mediated aeroallergen sensitization

BackgroundThe innate mechanisms associated with viral exacerbations in preschool children with recurrent wheezing are not understood. ObjectiveWe assessed differential gene expression in blood neutrophils from preschool children with recurrent wheezing, stratified by aeroallergen sensitization, at baseline and after exposure to polyinosinic-polycytidylic acid (poly(I:C)). We also examined whether poly(I:C)-stimulated blood neutrophils influenced airway epithelial gene expression. MethodsBlood neutrophils were purified and cultured overnight with poly(I:C) and underwent next-generation sequencing with Reactome pathway analysis. Primary human small airway epithelial cells were treated with poly(I:C)-treated neutrophil culture supernatants and were analyzed for Type 1 interferon gene expression with a targeted array. Symptoms and exacerbations were assessed in participants over 12 months. Results436 genes were differently expressed in neutrophils from children with versus without aeroallergen sensitization at baseline, with significant downregulation of Type 1 interferons. These Type 1 interferons were significantly upregulated in sensitized children after poly(I:C) stimulation. Confirmatory experiments demonstrated similar upregulation of Type 1 interferons in IL-4 treated neutrophils stimulated with poly(I:C). Poly(I:C)-treated neutrophil supernatants from children with aeroallergen sensitization also induced a Type 1 interferon response in epithelial cells. Children with aeroallergen sensitization also had higher symptom scores during exacerbations and these symptom differences persisted for three days after prednisolone treatment. ConclusionsType 1 interferon responses are dysregulated in preschool children with aeroallergen sensitization, which is in turn associated with exacerbation severity. Given the importance of Type 1 interferon signaling in viral resolution, additional studies of neutrophil Type 1 interferon responses are needed in this population.

Clonal Somatic Mutations in Chronic Lung Diseases Are Associated with Reduced Lung Function.

Rationale: Constantly exposed to the external environment and mutagens such as tobacco smoke, human lungs have one of the highest somatic mutation rates among all human organs. However, the relationship of these mutations to lung disease and function is not known. Objectives: To identify the prevalence and significance of clonal somatic mutations in chronic lung diseases. Methods: We analyzed the clonal somatic mutations from 1,251 samples of normal and diseased noncancerous lung tissue RNA sequencing with paired whole-genome sequencing from the Lung Tissue Research Consortium. We examined the associations of somatic mutations with lung function, disease status, and computationally deconvoluted cell types in two of the most common diseases represented in our dataset, chronic obstructive pulmonary disease (COPD; 29%) and idiopathic pulmonary fibrosis (IPF; 13%). Measurements and Main Results: Clonal somatic mutational burden was associated with reduced lung function in both COPD and IPF. We identified an increased prevalence of clonal somatic mutations in individuals with IPF compared with normal control subjects and individuals with COPD independent of age and smoking status. IPF clonal somatic mutations were enriched in disease-related and airway epithelial-expressed genes such as MUC5B in IPF. Patients who were MUC5B risk variant carriers had increased odds of developing somatic mutations of MUC5B that were explained by increased expression of MUC5B. Conclusions: Our identification of an increased prevalence of clonal somatic mutation in diseased lung that correlates with airway epithelial gene expression and disease severity highlights for the first time the role of somatic mutational processes in lung disease genetics.

Nepalese indoor cookstove smoke extracts alter human airway epithelial gene expression, DNA methylation and hydroxymethylation

Household air pollution caused by inefficient cooking practices causes 4 million deaths a year worldwide. In Nepal, 86% of the rural population use solid fuels for cooking. Over 25% of premature deaths associated with air pollution are respiratory in nature. Here we aimed to identify molecular signatures of different cookstove and fuel type exposures in human airway epithelial cells, to understand the mechanisms mediating cook stove smoke induced lung disease. Primary human airway epithelial cells in submerged culture were exposed to traditional cook stove (TCS), improved cook stove (ICS) and liquefied petroleum gas (LPG) stove smoke extracts. Changes to gene expression, DNA methylation and hydroxymethylation were measured by bulk RNA sequencing and HumanMethylationEPIC BeadChip following oxidative bisulphite conversion, respectively. TCS smoke extract alone reproducibly caused changes in the expression of 52 genes enriched for oxidative stress pathways. TCS, ICS and LPG smoke extract exposures were associated with distinct changes to DNA methylation and hydroxymethylation. A subset of TCS induced genes were associated with differentially methylated and/or hydroxymethylated CpGs sites, and enriched for the ferroptosis pathway and the upstream regulator NFE2L2. DNA methylation and hydroxymethylation changes not associated with a concurrent change in gene expression, were linked to biological processes and molecular pathways important to airway health, including neutrophil function, transforming growth factor beta signalling, GTPase activity, and cell junction organisation. Our data identified differential impacts of TCS, ICS and LPG cook stove smoke on the human airway epithelium transcriptome, DNA methylome and hydroxymethylome and provide further insight into the association between indoor air pollution exposure and chronic lung disease mechanisms.

Cis- and trans-eQTM analysis reveals novel epigenetic and transcriptomic immune markers of atopic asthma in airway epithelium

Expression quantitative trait methylation (eQTM) analyses uncover associations between DNA methylation markers and gene expression. Most eQTM analyses of complex diseases have focused on cis-eQTM pairs (within 1 megabase). This study sought to identify cis- and trans-methylation markers associated with gene expression in airway epithelium from youth with and without atopic asthma. In this study, the investigators conducted both cis- and trans-eQTM analyses in nasal (airway) epithelial samples from 158 Puerto Rican youth with atopic asthma and 100 control subjects without atopy or asthma. The investigators then attempted to replicate their findings in nasal epithelial samples from 2 studies of children, while also examining whether their results in nasal epithelium overlap with those from an eQTM analysis in white blood cells from the Puerto Rican subjects. This study identified 9,108 cis-eQTM pairs and 2,131,500 trans-eQTM pairs. Trans-associations were significantly enriched for transcription factor and microRNA target genes. Furthermore, significant cytosine-phosphate-guanine sites (CpGs) were differentially methylated in atopic asthma and significant genes were enriched for genes differentially expressed in atopic asthma. In this study, 50.7% to 62.6% of cis- and trans-eQTM pairs identified in Puerto Rican youth were replicated in 2 smaller cohorts at false discovery rate-adjusted P< .1. Replicated genes in the trans-eQTM analysis included biologically plausible asthma-susceptibility genes (eg, HDC, NLRP3, ITGAE, CDH26, and CST1) and are enriched in immune pathways. Studying both cis- and trans-epigenetic regulation of airway epithelial gene expression can identify potential causal and regulatory pathways or networks for childhood asthma. Trans-eQTM CpGs may regulate gene expression in airway epithelium through effects on transcription factor and microRNA target genes.

WS18.06 The impact of cytomegalovirus on airway epithelial gene expression

Introduction: The pathogenicity of cytomegalovirus (CMV) has predominantly been recognized in immunosuppressed populations & newborns. However, a recent study in persons with CF observed that CMV seropositive individuals died from lung disease 8 years sooner than those who were seronegative. Further, CMV DNA was identified in those individuals’ sputum (Parkins 2019, ERJ). This observation aligns with increasing evidence associating CMV infection with worse outcomes in many inflammatory conditions, and a higher overall mortality. To date, models of CMV infection of human bronchial epithelial (HBE) cells have not been established, limiting our ability to study its impact in CF airways.

Transcriptomic Response of Primary Human Bronchial Cells to Repeated Exposures of Cigarette and ENDS Preparations.

Cigarette smoke deregulates several biological pathways by modulating gene expression in airway epithelial cells and altering the physiology of the airway epithelium. The effects of repeated exposures of electronic cigarette delivery systems (ENDS) on gene expression in airway epithelium are relatively unknown. In order to assess the effect of repeated exposures of ENDS, primary normal human bronchial epithelial (NHBE) cells grown at air-liquid interface (ALI) were exposed to cigarette and ENDS preparations daily for 10 days. Cigarette smoke preparations significantly altered gene expression in a dose-dependent manner compared to vehicle control, including genes linked to oxidative stress, xenobiotic metabolism, cancer pathways, epithelial-mesenchymal transition, fatty acid metabolism, degradation of collagen and extracellular matrix, O-glycosylation, and chemokines/cytokines, which are known pathways found to be altered in smokers. Conversely, ENDS preparations had minimal effect on transcriptional pathways. This study revealed that a sub-chronic exposure of primary NHBE cultures to cigarette and ENDS preparations differentially regulated genes and canonical pathways, with minimal effect observed with ENDS preparations compared to cigarette preparations. This study also demonstrates the versatility of primary NHBE cultures at ALI to evaluate repeat-dose exposures of tobacco products.

Potentiation of long-acting \u03b22-agonist and glucocorticoid responses in human airway epithelial cells by modulation of intracellular cAMP

IntroductionOver 300 million people in the world live with asthma, resulting in 500,000 annual global deaths with future increases expected. It is estimated that around 50–80% of asthma exacerbations are due to viral infections. Currently, a combination of long-acting beta agonists (LABA) for bronchodilation and glucocorticoids (GCS) to control lung inflammation represent the dominant strategy for the management of asthma, however, it is still sub-optimal in 35–50% of moderate-severe asthmatics resulting in persistent lung inflammation, impairment of lung function, and risk of mortality. Mechanistically, LABA/GCS combination therapy results in synergistic efficacy mediated by intracellular cyclic adenosine monophosphate (cAMP).HypothesisIncreasing intracellular cAMP during LABA/GCS combination therapy via inhibiting phosphodiesterase 4 (PDE4) and/or blocking the export of cAMP by ATP Binding Cassette Transporter C4 (ABCC4), will potentiate anti-inflammatory responses of mainstay LABA/GCS therapy.MethodsExpression and localization experiments were performed using in situ hybridization and immunohistochemistry in human lung tissue from healthy subjects, while confirmatory transcript and protein expression analyses were performed in primary human airway epithelial cells and cell lines. Intervention experiments were performed on the human airway epithelial cell line, HBEC-6KT, by pre-treatment with combinations of LABA/GCS with PDE4 and/or ABCC4 inhibitors followed by Poly I:C or imiquimod challenge as a model for viral stimuli. Cytokine readouts for IL-6, IL-8, CXCL10/IP-10, and CCL5/RANTES were quantified by ELISA.ResultsUsing archived human lung and human airway epithelial cells, ABCC4 gene and protein expression were confirmed in vitro and in situ. LABA/GCS attenuation of Poly I:C or imiquimod-induced IL-6 and IL-8 were potentiated with ABCC4 and PDE4 inhibition, which was greater when ABCC4 and PDE4 inhibition was combined. Modulation of cAMP levels had no impact on LABA/GCS modulation of Poly I:C-induced CXCL10/IP-10 or CCL5/RANTES.ConclusionModulation of intracellular cAMP levels by PDE4 or ABCC4 inhibition potentiates LABA/GCS efficacy in human airway epithelial cells challenged with viral stimuli. The data suggest further exploration of the value of adding cAMP modulators to mainstay LABA/GCS therapy in asthma for potentiated anti-inflammatory efficacy.

Integrative Analysis Reveals a miRNA-mRNA Regulatory Network and Potential Causative Agents in the Asthmatic Airway Epithelium.

BackgroundDuring asthma progression, the intricate molecular networks, including microRNA (miRNA) transcriptional regulation in airway epithelium, remain largely undefined. The abnormal expression of miRNAs in asthmatic airway epithelium is a recent and fast-growing area in developing diagnostic and therapeutic targets for asthma.Material and MethodsAnalyses were conducted to compare airway epithelial miRNAs and gene expression between patients with asthma and healthy subjects from three datasets (two for miRNAs expression profiles and one for gene expression profile). The interactions network between differentially expressed (DE)-miRNAs and mRNAs was further identified for functional analysis. To verify the involvement and functions of all the identified miRNAs in asthma, we constructed two cellular models of asthma. The most promising causal miRNA candidate, miR-1246, was examined in an in vitro system to explore its targets and roles in asthma pathophysiology.ResultsThrough integrative analysis, we obtained six miRNAs with 31 validated target genes in airway epithelium associated with asthma. Next, we confirmed that these miRNAs were all associated with asthma progression by in vitro functional experiments. They may participate in eosinophilic inflammation (miR-92b-3p, miR-1246, miR-197-3p, and miR-124-5p) or remodeling (miR-197-3p, miR-193a-5p, miR-1246, and miR-92b-3p). Additionally, some other non-screened valuable miRNAs were also examined and identified (miR-21-5p and miR-19b-3p), and some detected in blood correlated with the disease status. Furthermore, we found that miR-1246 could directly target POSTN and influence epithelial-to-mesenchymal transition and fibrosis in airway epithelial cells.ConclusionWe constructed a preliminary epithelial regulatory network in asthma based on six identified miRNAs and their valuable target genes. Candidate factors in the biological miRNA-mRNA network in airway epithelium may provide further information on the pathogenesis of asthma. Strikingly, among all screened miRNAs, miR-1246, which could interact with POSTN may have multifunctional effects in the course of asthma and be a promising agent for asthma treatment and molecular subtyping.

Improving lung cancer risk stratification leveraging whole transcriptome RNA sequencing and machine learning across multiple cohorts

BackgroundBronchoscopy for suspected lung cancer has low diagnostic sensitivity, rendering many inconclusive results. The Bronchial Genomic Classifier (BGC) was developed to help with patient management by identifying those with low risk of lung cancer when bronchoscopy is inconclusive. The BGC was trained and validated on patients in the Airway Epithelial Gene Expression in the Diagnosis of Lung Cancer (AEGIS) trials. A modern patient cohort, the BGC Registry, showed differences in key clinical factors from the AEGIS cohorts, with less smoking history, smaller nodules and older age. Additionally, we discovered interfering factors (inhaled medication and sample collection timing) that impacted gene expressions and potentially disguised genomic cancer signals.MethodsIn this study, we leveraged multiple cohorts and next generation sequencing technology to develop a robust Genomic Sequencing Classifier (GSC). To address demographic composition shift and interfering factors, we synergized three algorithmic strategies: 1) ensemble of clinical dominant and genomic dominant models; 2) development of hierarchical regression models where the main effects from clinical variables were regressed out prior to the genomic impact being fitted in the model; and 3) targeted placement of genomic and clinical interaction terms to stabilize the effect of interfering factors. The final GSC model uses 1232 genes and four clinical covariates – age, pack-years, inhaled medication use, and specimen collection timing.ResultsIn the validation set (N = 412), the GSC down-classified low and intermediate pre-test risk subjects to very low and low post-test risk with a specificity of 45% (95% CI 37–53%) and a sensitivity of 91% (95%CI 81–97%), resulting in a negative predictive value of 95% (95% CI 89–98%). Twelve percent of intermediate pre-test risk subjects were up-classified to high post-test risk with a positive predictive value of 65% (95%CI 44–82%), and 27% of high pre-test risk subjects were up-classified to very high post-test risk with a positive predictive value of 91% (95% CI 78–97%).ConclusionsThe GSC overcame the impact of interfering factors and achieved consistent performance across multiple cohorts. It demonstrated diagnostic accuracy in both down- and up-classification of cancer risk, providing physicians actionable information for many patients with inconclusive bronchoscopy.

Intra-tracheal delivery of AAV6 vectors results in sustained transduction in murine lungs without genomic integration.

Despite the progress made in AAV-based gene therapy targeting different organ systems, lung-targeted gene therapy using AAV vectors has not been effective, mostly due to the poor transduction and un-sustained gene expression in airway epithelium. Furthermore, concerns over possible harmful insertional mutagenesis seen in other cell types, particularly hepatocytes, raised a question about AAV safety. In this study, we evaluate the long-term persistence of this vector in mouse lungs and any possible harmful integration of these vectors into the host genome. AAV6 vectors expressing reporter gene (firefly luciferase) were delivered to the lungs of C57BL/6 mice through intra-tracheal intubation. Despite the large variation among individual animals, most animals had high and sustained luciferase activity with a peak from 2 to 3 weeks post-transduction before a significant decline between 15 and 19 weeks post-transduction. More importantly, even after its decline, most animals maintained detectable luciferase expression for 150 days or more, which was confirmed by post-necropsy qPCR analysis of luciferase gene expression. At the termination point of experiments, an average of one copy of AAV expression cassette per mouse genome was detected. We also found that partial overlaps between the AAV6 expression cassette and the mouse genome were distributed broadly with no apparent systematic preference in any mouse chromosomal map location. In summary, our data suggest that AAV6 mediated long-term gene expression in the lungs with no evidence of genomic integration, and thus, any insertional mutagenesis.

Expression Quantitative Trait Methylation Analysis Reveals Methylomic Associations With Gene Expression in Childhood Asthma

Nasal (airway) epithelial methylation profiles have been associated with asthma, but the effects of such profiles on expression of distant cis-genes are largely unknown. To identify genes whose expression is associated with proximal and distal CpG probes (within 1 Mb), and to assess whether and how such genes are differentially expressed in atopic asthma. Genome-wide expression quantitative trait methylation (eQTM) analysis in nasal epithelium from Puerto Rican subjects (aged 9-20 years) with (n= 219) and without (n= 236) asthma. After the eQTM analysis, a Gene Ontology Enrichment analysis was conducted for the top 500 eQTM genes, and mediation analyses were performed to identify paths from DNA methylation to atopic asthma through gene expression. Asthma was defined as physician-diagnosed asthma and wheeze in the previous year, and atopy was defined as at least one positive IgE to allergens. Atopic asthma was defined as the presence of both atopy and asthma. We identified 16,867 significant methylation-gene expression pairs (false-discovery rate-adjusted P< .01) in nasal epithelium from study participants. Most eQTM methylation probes were distant (average distance, ∼378 kb) from their target genes, and also more likely to be located in enhancer regions of their target genes in lung tissue than control probes. The top 500 eQTM genes were enriched in pathways for immune processes and epithelial integrity and were more likely to have been previously identified as differentially expressed in atopic asthma. In a mediation analysis, we identified 5,934 paths through which methylation markers could affect atopic asthma through gene expression in nasal epithelium. Previous epigenome-wide association studies of asthma have estimated the effects of DNA methylation markers on expression of nearby genes in airway epithelium. Our findings suggest that distant epigenetic regulation of gene expression in airway epithelium plays a role in atopic asthma.

Gene Expression Alterations in the Bronchial Epithelium of e-Cigarette Users

Although e-cigarette (ECIG) use has increased in the United States, their potential health effects remain uncertain. Understanding the effects of tobacco cigarette (TCIG) smoke on bronchial airway epithelial gene expression have previously provided insights into tobacco-related disease pathogenesis. Identifying the impact of ECIGs on airway gene expression could provide insights into their potential long-term health effects. We sought to compare the bronchial airway gene-expression profiles of former TCIG smokers now using ECIGs with the profiles of former and current TCIG smokers. We performed gene-expression profiling of bronchial epithelial cells collected from current TCIG smokers (n= 9), current ECIG users who are former TCIG smokers (n= 15), and former TCIG smokers (n= 21). We then compared our findings with previous studies of the effects of TCIG use on bronchial epithelium, as well an invitro model of ECIG exposure. Among 3,165 genes whose expression varied between the three study groups (q< 0.05), we identified 468 genes altered in ECIG users relative to former smokers (P< .05). Seventy-nine of these genes were up- or down-regulated concordantly among ECIG and TCIG users. We did not detect ECIG-associated gene-expression changes in known pathways associated with TCIG usage. Genes downregulated in ECIG users are enriched among the genes most downregulated by exposure of airway epithelium to ECIG vapor invitro. ECIGs induce both distinct and shared patterns of gene expression relative to TCIGs in the bronchial airway epithelium. The concordance of the genes altered in ECIG users and in the invitro study suggests that genes altered in ECIG users are likely to be changed as the direct effect of ECIG exposure.

MEG3 is increased in idiopathic pulmonary fibrosis and regulates epithelial cell differentiation

Idiopathic pulmonary fibrosis (IPF) is a chronic interstitial lung disease causing fibrotic remodeling of the peripheral lung, leading to respiratory failure. Peripheral pulmonary epithelial cells lose normal alveolar epithelial gene expression patterns and variably express genes associated with diverse conducting airway epithelial cells, including basal cells. Single-cell RNA sequencing of pulmonary epithelial cells isolated from IPF lung tissue demonstrated altered expression of LncRNAs, including increased MEG3. MEG3 RNA was highly expressed in subsets of the atypical IPF epithelial cells and correlated with conducting airway epithelial gene expression patterns. Expression of MEG3 in human pulmonary epithelial cell lines increased basal cell-associated RNAs, including TP63, KRT14, STAT3, and YAP1, and enhanced cell migration, consistent with a role for MEG3 in regulating basal cell identity. MEG3 reduced expression of TP73, SOX2, and Notch-associated RNAs HES1 and HEY1, in primary human bronchial epithelial cells, demonstrating a role for MEG3 in the inhibition of genes influencing basal cell differentiation into club, ciliated, or goblet cells. MEG3 induced basal cell genes and suppressed genes associated with terminal differentiation of airway cells, supporting a role for MEG3 in regulation of basal progenitor cell functions, which may contribute to tissue remodeling in IPF.

Abstract 2531: Diagnostic evaluation of indeterminate pulmonary nodules via rna-seq of bronchial epithelium

Abstract RATIONALE: Diagnosis of the estimated 1.5 million indeterminate pulmonary nodules found annually on chest CT in the US poses a significant clinical challenge. Prior work from our group has identified alterations in the bronchial epithelial cell microarray-derived transcriptome which are associated with lung cancer and which can serve as a clinically useful biomarker for the evaluation of patients undergoing bronchoscopy for suspect lung cancer. Given that RNA sequencing has a broader measurement range and allows for unbiased detection of novel transcripts, we evaluated the ability of bronchial epithelial gene-expression from RNA sequencing to identify patients with lung cancer among individuals with indeterminate pulmonary nodules. METHODS: Bronchial epithelial brushings were collected from current and former smokers (n=93; age&amp;gt;45, pack-year&amp;gt;20) undergoing diagnostic work up for indeterminate pulmonary nodules (7-30 mm in diameter) at 4 military and 7 VA hospitals within the DECAMP consortium. Patients were followed clinically for up to two years after sample collection until a final diagnosis of lung cancer (n=52) or benign disease (n=41) was made. We performed total RNA sequencing using the Illumina TruSeq Stranded Total RNA Sample Preparation kit. Individual samples were sequenced to generate paired-end reads per sample. Reads were aligned to the human assembly Genome Reference Consortium Human Build 37 (GRCh37) in Ensembl using Spliced Transcripts Alignment to a Reference (STAR) and summarized into count data via RNA-Seq by Expectation Maximization (RSEM). RESULTS: We found no differences in age, gender, or smoking status between patients with malignant vs. benign nodules. Using a generalized linear model accounting for smoking status, we identified 34 differentially expressed genes (DEGs, FDR q &amp;lt; 0.2) from the bronchial epithelium associated with lung cancer. Of these 34 DEGs, 28 were up-regulated in individuals with lung cancer and enriched for pathways related to the response to hypoxia and p53. The 6 down-regulated DEGs were enriched for immune cell activation, supporting an immunosuppressive field effect of lung cancer. Concordant enrichment of the previous microarray-based bronchial gene-expression classifier was observed using Gene Set Enrichment Analysis (GSEA, FDR q &amp;lt; 0.05), supporting the cross-platform and cross-cohort validity of the microarray-based classifier. CONCLUSION: Taken together, findings from this pilot study suggest that there are distinct lung-cancer associated airway epithelial gene-expression alterations detected by RNA-sequencing which replicate a previously published microarray-based signature. Additional alterations identified by our RNAseq studies may be leveraged to refine and further develop an RNA-sequencing based airway biomarker for the early detection of lung cancer in high-risk smokers. Citation Format: Xingyi Shi, Ehab Billatos, Jiarui Zhang, Jennifer Beane, Elizabeth Moses, Gang Liu, Christopher Stevenson, Marc E. Lenburg, Avrum Spira. Diagnostic evaluation of indeterminate pulmonary nodules via rna-seq of bronchial epithelium [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 2531.

Impact of acute exposure to cigarette smoke on airway gene expression

Understanding effects of acute smoke exposure (ASE) on airway epithelial gene expression and their relationship with the effects of chronic smoke exposure may provide biological insights into the development of smoking-related respiratory diseases. Bronchial airway epithelial cell brushings were collected from 63 individuals without recent cigarette smoke exposure and before and 24 h after smoking three cigarettes. RNA from these samples was profiled on Affymetrix Human Gene 1.0 ST microarrays. We identified 91 genes differentially expressed 24 h after ASE (false discovery rate < 0.25). ASE induced genes involved in xenobiotic metabolism, oxidative stress, and inflammation and repressed genes related to cilium morphogenesis and cell cycle. While many genes altered by ASE are altered similarly in chronic smokers, metallothionein genes are induced by ASE and suppressed in chronic smokers. Metallothioneins are also suppressed in current and former smokers with lung cancer relative to those without lung cancer. Acute exposure to as little as three cigarettes and chronic smoking induce largely concordant changes in airway epithelial gene expression. Differences in short-term and long-term effects of smoking on metallothionein expression and their relationship to lung cancer requires further study given these enzymes' role in the oxidative stress response.

Stability of gene expression by primary bronchial epithelial cells over increasing passage number

BackgroundAn increasing number of studies using primary human bronchial epithelial cells (BECs) have reported intrinsic differences in the expression of several genes between cells from asthmatic and non-asthmatic donors. The stability of gene expression by primary BECs with increasing cell passage number has not been well characterized.MethodsTo determine if expression by primary BECs from asthmatic and non-asthmatic children of selected genes associated with airway remodeling, innate immune response, immunomodulatory factors, and markers of differentiated airway epithelium, are stable over increasing cell passage number, we studied gene expression patterns in passages 1, 2, 3, 4, and 5 BECs from asthmatic (n = 6) and healthy (n = 6) subjects that were differentiated at an air-liquid interface. RNA was harvested from BECs and RT-PCR was performed for TGFβ1, TGFβ2, activin A, FSTL3, MUC5AC, TSLP, IL-33, CXCL10, IFIH1, p63, KT5, TUBB4A, TJP1, OCLN, and FOXJ1.ResultsExpression of TGFβ1, TGFβ2, activin A, FSTL3, MUC5AC, CXCL10, IFIH1, p63, KT5, TUBB4A, TJP1, OCLN, and FOXJ1 by primary BECs from asthmatic and healthy children was stable with no significant differences between passages 1, 2 and 3; however, gene expression at cell passages 4 and 5 was significantly greater and more variable compared to passage 1 BECs for many of these genes. IL-33 and FOXJ1 expression was also stable between passages 1 through 3, however, expression at passages 4 and 5 was significantly lower than by passage 1 BECs. TSLP, p63, and KRT5 expression was stable across BEC passages 1 through 5 for both asthmatic and healthy BECs.ConclusionsThese observations illustrate the importance of using BECs from passage ≤3 when studying gene expression by asthmatic and non-asthmatic primary BECs and characterizing the expression pattern across increasing cell passage number for each new gene studied, as beyond passage 3 genes expressed by primary BECs appear to less accurately model in vivo airway epithelial gene expression.

Meta-analysis of airway epithelium gene expression in asthma.

Differential gene expression in the airway epithelium of patients with asthma versus controls has been reported in several studies. However, there is no consensus on which genes are reproducibly affected in asthma. We sought to identify a consensus list of differentially expressed genes (DEGs) using a meta-analysis approach.We identified eight studies with data that met defined inclusion criteria. These studies comprised 355 cases and 193 controls and involved sampling either bronchial or nasal epithelium. We conducted study-level analyses, followed by a meta-analysis. Likewise, we applied a meta-analysis framework to the results of study-level pathway enrichment.We identified 1273 DEGs, 431 of which had not been identified in previous studies. 450 DEGs exhibited large effect sizes and were robust to study population differences in age, sex, race/ethnicity, medication use, smoking status and exacerbations. The magnitude of differential expression of these 450 genes was highly similar in bronchial and nasal airway epithelia. Meta-analysis of pathway enrichment revealed a number of consistently dysregulated biological pathways, including putative transcriptional and post-transcriptional regulators.In total, we identified a set of genes that is consistently dysregulated in asthma, that links to known and novel biological pathways, and that will inform asthma subtype identification.

Mandatory role of HMGA1 in human airway epithelial normal differentiation and post-injury regeneration.

Due to high levels of expression in aggressive tumors, high mobility group AT-hook 1 (HMGA1) has recently attracted attention as a potential anti-tumor target. However, HMGA1 is also expressed in normal somatic progenitor cells, raising the question: how might systemic anti-HMGA1 therapies affect the structure and function of normal tissue differentiation? In the present study, RNA sequencing data demonstrated HMGA1 is highly expressed in human airway basal stem/progenitor cells (BC), but decreases with BC differentiation in air-liquid interface cultures (ALI). BC collected from nonsmokers, healthy smokers, and smokers with chronic obstructive pulmonary disease (COPD) displayed a range of HMGA1 expression levels. Low initial expression levels of HMGA1 in BC were associated with decreased ability to maintain a differentiated ALI epithelium. HMGA1 down-regulation in BC diminished BC proliferation, suppressed gene expression related to normal proliferation and differentiation, decreased airway epithelial resistance, suppressed junctional and cell polarity gene expression, and delayed wound closure of airway epithelium following injury. Furthermore, silencing of HMGA1 in airway BC in ALI increased the expression of genes associated with airway remodeling in COPD including squamous, epithelial-mesenchymal transition (EMT), and inflammatory genes. Together, the data suggests HMGA1 plays a central role in normal airway differentiation, and thus caution should be used to monitor airway epithelial structure and function in the context of systemic HMGA1-targeted therapies.

Meta-analysis of airway epithelium gene expression in asthma

Meta-analysis of airway epithelium gene expression in asthma

Airway Mucosal Host Defense Is Key to Genomic Regulation of Cystic Fibrosis Lung Disease Severity.

The severity of cystic fibrosis (CF) lung disease varies widely, even for Phe508del homozygotes. Heritability studies show that more than 50% of the variability reflects non-cystic fibrosis transmembrane conductance regulator (CFTR) genetic variation; however, the full extent of the pertinent genetic variation is not known. We sought to identify novel CF disease-modifying mechanisms using an integrated approach based on analyzing "in vivo" CF airway epithelial gene expression complemented with genome-wide association study (GWAS) data. Nasal mucosal RNA from 134 patients with CF was used for RNA sequencing. We tested for associations of transcriptomic (gene expression) data with a quantitative phenotype of CF lung disease severity. Pathway analysis of CF GWAS data (n = 5,659 patients) was performed to identify novel pathways and assess the concordance of genomic and transcriptomic data. Association of gene expression with previously identified CF GWAS risk alleles was also tested. Significant evidence of heritable gene expression was identified. Gene expression pathways relevant to airway mucosal host defense were significantly associated with CF lung disease severity, including viral infection, inflammation/inflammatory signaling, lipid metabolism, apoptosis, ion transport, Phe508del CFTR processing, and innate immune responses, including HLA (human leukocyte antigen) genes. Ion transport and CFTR processing pathways, as well as HLA genes, were identified across differential gene expression and GWAS signals. Transcriptomic analyses of CF airway epithelia, coupled to genomic (GWAS) analyses, highlight the role of heritable host defense variation in determining the pathophysiology of CF lung disease. The identification of these pathways provides opportunities to pursue targeted interventions to improve CF lung health.