Skip to main content

Comparative analyses of long non-coding RNA profiles in vivo in cystic fibrosis lung airway and parenchyma tissues

Abstract

Background

Recent advances in the functional analyses of endogenous non-coding RNA (ncRNA) molecules, including long non-coding RNAs (LncRNAs), have provided a new perspective on the crucial roles of RNA in gene regulation. Consequently, LncRNA deregulation is a key factor in various diseases, including pulmonary disorders like Cystic Fibrosis (CF). CF is the most common life limiting recessive disease in the U.S., and is due to mutations in the CFTR gene. CF mutations, of which the most common is F508del-CFTR, prevents correct folding, trafficking and function of the mutant CFTR protein and is further manifested by the hyper-expression of pro-inflammatory cytokines and chemokines into the airway lumen leading to bronchiectasis and culminating in lung destruction.

Methods

Here we report a distinct LncRNA signature and corresponding mRNAs that distinguishes CF lung (airway and parenchyma) tissues from matched non-CF controls (n = 4 each group), generated by microarray specific for LncRNAs which includes corresponding mRNA expressions. In silico analyses of the cellular processes that are impacted by these LncRNAs was performed using Gene Ontology (GO). A selected subset of LncRNAs were validated by quantitative real-time PCR.

Results

We have identified 636 LncRNAs differentially expressed in CF airway epithelium and 1974 in CF lung parenchyma compared to matched non-CF controls (fold change ≥2, p < 0.05), majority of which (> 50%) are intergenic. Interestingly, 15 of these differentially expressed LncRNAs and 9 coding mRNAs are common to airway and parenchyma tissues. GO analyses indicates that signaling pathways and cell membrane functions are significantly affected by the alteration in LncRNA expressions in CF lung tissues. Seven of the differentially expressed LncRNAs, exhibit similar expression trends in CFBE41o- compared to control cells.

Conclusion

Understanding the mechanisms by which these LncRNAs regulate CF disease phenotype will help develop novel therapeutic targets for CF and related pulmonary diseases, such as COPD and Asthma.

Introduction

Cystic Fibrosis (CF) is the most common life limiting recessive disease in the U.S. and is due to mutations in the CFTR gene. CF mutations, the most common of which is F508del-CFTR, cause a massive pro-inflammatory phenotype in the lung arising from dys-regulated expression of inflammatory genes. Recently, endogenous non-coding RNA (ncRNA) molecules, including long non-coding RNAs (LncRNAs), have emerged as important therapeutic targets at the frontier of biomedical research. These LncRNAs coordinate with epigenetic factors to play a crucial role in the regulation of biological processes as well as in diseases. LncRNAs have recently emerged as novel epigenetic regulators of gene expression, including inflammatory genes. The recent FANTOM Atlas [1] has shown that many LncRNAs come from upstream enhancer elements and are functional.

LncRNAs are defined as non-coding RNAs (ncRNAs) that are transcribed by RNA polymerase II, and are at least 200 nucleotides in length [2]. LncRNAs do not have the ability to encode proteins, and they include all ncRNAs longer than 200 nucleotides (except rRNA and tRNA). These RNA molecules have provided a new perspective on the roles of RNAs in gene regulation [3, 4]. LncRNAs mostly originate within a 2-kb region surrounding the Transcription Start Site (TSS) of protein-coding genes and some originate from more distal (> 2 kb) unannotated regions. Thus LncRNAs are classified according to their position relative to protein-coding genes, and are divided into five classes [5, 6]: (i) intronic LncRNAs are located within an intron of a protein-coding gene in either direction; (ii) long intergenic ncRNAs (LincRNA) are separated by transcriptional units from protein-coding genes; (iii) bidirectional LncRNAs are transcribed in opposite directions in relation to the promoter of a protein-coding gene; (iv) antisense LncRNAs are transcribed across the exons of protein-coding genes from the opposite direction; and (v) transcribed pseudogene LncRNAs are transcribed from a gene without the ability to produce a protein.

Some LncRNAs are preferentially expressed in immune cells and play important roles in immune cell development [7]. LncRNAs are considered to be more species-, tissue- and developmental stage-specific than mRNAs [8]. Several studies have demonstrated that LncRNA deregulation has a role in various diseases [9, 10], including pulmonary disorders [11, 12]. One study indicates the role of LncRNAs in the regulation of mutant CFTR and its impact on the CF disease phenotype [13]. Suppression of LncRNA in intron 11 of CFTR, called BGAS (BG213071), or repression of its protein binding partners, has been shown to induce a 4-fold increase in mutant CFTR at the cell surface [13]. LncRNA expressions have been analyzed in CF bronchial epithelium [14] and in primary CF bronchial epithelial cells infected with Pseudomonas aeruginosa [15].

Here we have identified LncRNAs that are differentially expressed in lung airway and parenchyma tissues isolated from CF patients undergoing lung transplant compared to matched non-CF control tissues. Further analysis of the expression of selected subset of LncRNAs was performed by TaqMan-based qPCR assays. Seven LncRNAs exhibit similar expression trend in the CF epithelial cell line, CFBE41o-. The CF disease-specific LncRNA signature includes > 50% intergenic LncRNAs. Bioinformatic analyses of these differentially expressed CF-specific LncRNAs indicate their impact on CF-relevant cellular and biological processes, including cell membrane function as well as signaling pathways. Collectively, these LncRNAs are likely to provide new insights into epigenetic mechanisms that regulate CF disease phenotype. Understanding the associated mechanisms will ultimately help identify novel therapeutic targets for CF and related pulmonary disorders.

Materials and methods

Reagents

The following reagents were used: α-MEM (Sigma, M2279), 0.25% Trypsin-EDTA (Sigma, T3924), Fetal Bovine Serum (Millipore, ES-009-B), L-Glutamine (Millipore, TMS-002-C), Penicillin-Streptomycin solution (Millipore, TMS-AB2-C), Hygromycin B (Sigma, H0654-500MG) and miRVana kit (Ambion, AM1560).

Study populations and human specimens

Lung tissues from eight subjects undergoing lung transplant were obtained for this study from the University of Pittsburgh Cystic Fibrosis Research Center: 4 CF and 4 matched non-CFs, with mean ages of 29 ± 3.8 yr and 24.4 ± 6.4 yr, respectively. Bronchial epithelial and parenchyma tissues were collected from each of these individuals, under a University of Pittsburgh Institutional Review Board approved protocol.

The procedure for collection of bronchial epithelial and parenchymal tissues was similar to that described [16]. Following lung resection, the bronchial tree from the second to sixth generation was dissected and rinsed for 24 h at 4 °C in Eagle’s MEM/HEPES on a rocker to remove any blood and mucus. After 36–48 h exposure to 0.1% protease XIV and 1% DNase in EMEM/HEPES, bronchial epithelial samples were collected by gently scraping the luminal surface using a surgical blade. Samples were then treated with Accutase and passed through a 0.22um cell filter. Microscopic evaluation of these samples revealed that this approach isolates a pure population of bronchial epithelia cells without contamination from the underlying lamina propria, which was intact following this procedure. Distal lung parenchyma samples were obtained from a lobe showing no obvious pathology or inflammation. Samples collected in this manner were snap frozen in liquid nitrogen and stored at -80 °C until RNA extraction.

RNA profiling and statistical analyses of data

LncRNAs and mRNA profiling studies were performed by Arraystar, Inc. (Rockville, MD, USA) on RNA isolated from lung tissues. Total RNA from each sample was quantified by the NanoDrop Spectrophotometer ND-1000, and RNA integrity was assessed by agarose gel electrophoresis as well as by bioanalyzer (as part of Arraystar services). Agilent Feature Extraction software (version 11.0.1.1) was used to analyze the acquired array data. Quantile normalization and subsequent data processing were performed using the GeneSpring GX v12.1 software package (Agilent Technologies). Differentially expressed LncRNAs and mRNAs with statistical significance were identified through Fold Change filtering between two sample groups. GO analysis was applied to determine the roles that these differentially expressed mRNAs played in the identified biological pathways. Finally, Hierarchical Clustering was performed to show the distinguishable LncRNA and mRNA expression patterns among samples.

Cell culture

CFBE41o- (Millipore, SCC160) and CFBE41o-6.2WT-CFTR (Millipore, SCC151) cells, obtained from Millipore Sigma were used for these studies. Both the cell lines were grown in Collagen I coated T75 flask (Thermo Scientific, 132,707) in minimum Eagle’s medium and supplemented with 10% Fetal Bovine Serum, 2 mM L-glutamine, 1X Penicillin-Streptomycin and 300 μg/ml Hygromycin B (CFBE41o-6.2WT-CFTR). All cells were maintained in a humidified 5% CO2 incubator at 37 °C.

RT-PCR and LncRNA assays

Total RNA was isolated using the mirVana miRNA Isolation Kit, following the manufacturer’s instructions. After RNA extraction, RNA samples were reverse transcribed by the High Capacity cDNA Reverse Transcription Kit (Applied Biosystems, 4,368,813). Real-time quantification of individual LncRNAs in CF cell lines was performed with specific TaqMan assays designed using the custom TaqMan assay design tool available from Thermo Fisher. Real-time PCR data were normalized to the endogenous β-actin control. The relative fold changes of LncRNAs were analyzed using the 2-ΔΔCT method.

Statistical analysis

Statistical analysis was performed using Excel. Significance values (p ≤ 0.05) were determined by student’s t-test. Error bars on graphs represent SEM.

Results

LncRNA expression profiles in CF lung airway and parenchyma tissues

To understand the role of LncRNAs in CF lung disease we performed comprehensive analyses of LncRNAs and corresponding mRNA expression profiles in lung tissues, both airway and parenchyma, obtained from CF patients undergoing lung transplant. These data were compared with those obtained from age (ranging from 23 to 36 years) and sex-matched (4 male and 4 female) healthy controls (Table 1). The expression of LncRNAs and mRNAs were analyzed by microarray (Human LncRNA Array v3.0, Arraystar, Inc.), which is comprised of ~ 40,173 LncRNAs and ~ 20,730 mRNAs. We have identified 636 differentially expressed LncRNAs in CF airway tissues (Fig. 1a) and 1974 differentially expressed LncRNAs in CF lung parenchyma tissues (Fig. 1b) compared to matched respective healthy controls (n = 4 each group, fold change ≥2, p < 0.05). As depicted in Fig. 1, when all these LncRNAs are compared using a hierarchical clustering algorithm, the dendrogram clearly distinguishes CF tissues from matched controls. The corresponding mRNA expression profiles also clearly distinguish CF tissues from controls (Fig. 1a and b).

Table 1 Patient lung tissue data
Fig. 1
figure1

Analyses of LncRNA in CF lung tissues. The expressions of LncRNA and corresponding mRNAs significantly (p < 0.05) altered in CF lung tissues compared to matched non-CF control tissues (n = 4 each) were analyzed by LncRNA arrays: a airway and b parenchyma tissues. Density plots depict the variation in expression: X and Y axes represent averaged normalized values in each group (log2). Volcano plots depict the fold-changes: each dot represents a single LncRNA or mRNA and is colored black unless it was differentially expressed. Hierarchical clustering and the heat map indicate the expression level of the transcripts significantly altered in CF tissues compared to non-CF controls: red represents increased expression, while green represents reduced expression. (CF-A, CF lung airway tissue; HBE-A, control airway tissue; CF-P, CF lung parenchyma tissue; HBE-P, control lung parenchyma)

Interestingly, 15 LncRNAs exhibit similar expression trends in CF lung airways as well as in CF parenchyma tissues (Fig. 2). The Venn diagrams (Fig. 2a) indicate that of the total number of up-regulated LncRNAs in CF tissues, 15 LncRNAs (0.7%) are common to both CF airway and parenchyma tissues, while 380 LncRNAs (18.8%) are unique to CF airway and 1628 (80.5%) are unique to CF parenchyma tissues compared to matched controls. Similarly, amongst those LncRNAs that are down-regulated in CF tissues; coincidentally, 15 LncRNAs (2.7%) are common to both CF airway and parenchyma tissues, while 226 LncRNAs (40.4%) are unique to CF airway and 318 (56.9%) are unique to CF parenchyma tissues compared to matched controls. Additionally, the analyses of corresponding mRNAs in these tissues indicate that two mRNAs are up-regulated and seven mRNAs are down-regulated both in airway and parenchyma tissues (Fig. 2b).

Fig. 2
figure2

Comparative analyses and classification of LncRNAs significantly altered in CF airway and parenchyma tissues. Venn diagrams indicate the differentially expressed (a) LncRNAs and b mRNAs that are common to airway and parenchyma tissues: up-regulated and down-regulated in CF compared to control tissues. The LncRNAs and mRNAs that exhibit similar expression pattern in the CF airway and parenchyma tissues are listed in the adjacent tables. The two LncRNAs that were validated in cell line models are highlighted. c The pie charts indicate the classification of the aberrantly expressed LncRNAs. The intergenic LncRNAs account for majority (58–63%) of the differences in CF airway as well as lung parenchyma tissues compared to non-CF control tissues

Classification of the differentially expressed LncRNAs in CF lung tissues

Analyses of the differentially expressed LncRNAs indicate that majority of the differentially expressed LncRNAs belong to the class of intergenic LncRNAs. The classification type as well as respective fold changes of the LncRNAs that are significantly up-regulated or down-regulated in CF tissues are listed in Tables 2 and 3. Of the top 10 LncRNAs that are significantly up- or down-regulated in the CF airway tissues, EEF1E1-BLOC1S5 and ARHGAP19-SLIT1 are read through transcripts, uc.363+, uc.8- and uc.476 are ultra-conserved LncRNAs, while the remaining have not been characterized. A few of the top 10 LncRNAs significantly altered in the CF parenchyma tissues have been characterized, and include the two up-regulated natural anti-sense non-coding transcripts, MIR3617 and POT1-AS1, and two intergenic LncRNAs, LINC00176 and LINC01023.

Table 2 LncRNA expressions in CF airway tissues
Table 3 LncRNA expressions in CF parenchyma tissues

As depicted by the Venn diagram (Fig. 3a), the intergenic LncRNAs account for the majority (58–63%) of the differences in CF airway as well as CF lung parenchyma tissues compared to respective control tissues. LncRNAs that are in the natural antisense and intronic-antisense categories accounted for about 13–15% of the total, while the remaining categories that include exon-sense overlapping, intron-sense overlapping and bi-directional accounted for < 10% of the total LncRNAs.

Fig. 3
figure3

Analyses of cellular processes affected by differentially expressed LncRNAs in CF lung airway and parenchyma tissues. Gene Ontology (GO) analyses was used to evaluate the impact of the aberrantly altered LncRNAs in (a) up-regulated in CF airway, b down-regulated in CF airway, c up-regulated in CF lung parenchyma, and d down-regulated in CF lung parenchyma tissues compared to respective non-CF controls. The significantly impacted Biological Processes, Cellular Component and Molecular Functions are depicted by corresponding histograms that indicate the enrichment score

Analyses of LncRNA-mRNA networks

In silico analyses of the differentially expressed LncRNAs and mRNAs networks were performed with Gene Ontology (GO) analyses program. Our data indicate that CF relevant biological processes are associated with signaling pathways; cellular components that include the intracellular vesicle, and membrane systems; and molecular functions including transcriptional activator activity are impacted in CF lung tissues compared to matched control tissues. Consistent with the CF disease phenotype, the array data indicate cell membrane function, including ion transport, is significantly affected by the aberrant expression of LncRNAs in CF lung tissues (Fig. 3a-d).

The Biological Processes (BP) that are significantly up-regulated in CF airway tissues, include regulation of inflammatory molecules (viz., IL-6, Interferon-beta), response to LPS as well as Toll-like Receptor (TLR) signaling. Those BP that are down-regulated significantly in CF airways, include regulation of cytosolic calcium ion, membrane transport functions. The two key BP that are up-regulated in CF parenchyma tissues include regulation of Insulin-like Growth Factor signaling and G-protein couple receptor signaling. Moreover, cell-cell adhesion, response to cytokine, vesicle-mediated transport, oxidation-reduction reactions are reduced in the CF parenchyma tissues. The overall analyses of the Cellular Component (CC) on CF airways and parenchyma tissues indicate significant impact on membrane structure, vesicles, organelles (both intracellular and extracellular), exosomes, and trans-golgi network. Consistently, the overall Molecular Function (MF) impacted in the CF lung tissues include ion transport, membrane transport activity, ligand binding activity (viz. LPS), membrane receptor activity (viz. G-protein receptor binding).

Validation of selected LncRNAs in CF epithelial cells

Subsequent to comprehensive analyses of LncRNAs in CF lung tissues, we validated a selected subset of LncRNAs in cell culture models. We thus performed quantitative real-time PCR-based validation of selected LncRNAs in F508del-CFTR CF airway epithelial cell line, CFBE41o-, compared to WT-CFTR control cell line. Seven of the LncRNAs that are differentially expressed in CF lung tissues compared to matched non-CF control tissues exhibit similar expression trends in CF cell lines (Fig. 4). These include RN7SKP237 and ARHGAP19-SLIT, that are up-regulated, and AC003092.1and LOC101928516 that are down-regulated, in CF airway tissues (Fig. 4a); LINC01023 which is down-regulated in CF lung parenchyma tissues (Fig. 4b); and two that are common to CF airway and parenchyma tissues, LINC01558 is up-regulated, and LOC100133050 is down-regulated (Fig. 4c). The data indicate that the expression levels of all these LncRNAs were significantly (p < 0.05) different between the CF bronchial epithelial cell line, CFBE41o- as compared to that in control cells. Further studies are focused on functional analyses of these LncRNAs.

Fig. 4
figure4

Validation of LncRNAs. The levels of LncRNAs differentially expressed in (a) CF airway tissues, b CF parenchyma tissues and c those that exhibit similar trend in both CF airway and parenchyma tissues were analyzed (by TaqMan qPCR assays) in F508del-CFTR CF cell line CFBE410-, and in control cell line containing WT-CFTR

Discussion

LncRNAs have emerged as key regulators of cellular processes and have been found to be aberrantly expressed in various diseases. The role of LncRNAs in CF has not been studied extensively. Here we report a LncRNA signature that distinguishes F508del-CFTR CF lung airway and parenchyma tissues from matched controls (n = 4 each group). The CF airway tissues have 636 differentially expressed LncRNAs, while the CF lung parenchyma tissues exhibited 1974 differentially expressed LncRNAs, compared to respective matched control tissues. The majority (> 50%) of these LncRNAs belong to the class of intergenic LncRNAs (LincRNAs). Moreover, 15 of these LncRNAs are elevated in both types of CF lung tissues, airway as well as parenchyma tissues, and coincidentally 15 LncRNAs also exhibit reduced expression in both types of CF tissues. Interestingly, we also find mRNAs that exhibit a similar trend in the CF airway and parenchyma tissues. Though we were able to analyze a limited number of clinical samples (n = 4 per group), we find a distinct LncRNA signature for CF airway and parenchyma tissues.

Comprehensive analyses of LncRNAs in CF tissues have identified two read through non-coding transcripts, EEF1E1-BLOC1S5 and ARHGAP19-SLIT1, an antisense LncRNA RN7SKP237, and 3 ultra-conserved LncRNAs, uc.363+, uc.8- and uc.476, that are significantly up- or down-regulated in the CF airway tissues. While in the CF parenchyma tissues two natural anti-sense non-coding transcripts, MIR3617 and POT1-AS1, are up-regulated, and two intergenic LncRNAs, LINC00176 and LINC01023 exhibit opposite trend. While LINC00176 is a myc-target and is up-regulated, LINC01023 is known to be associated with IGF1R/Akt pathway and is down-regulated in the CF lung parenchyma tissues. The other significantly altered LncRNAs (see Tables 2 and 3) have not been characterized. To date there have been two reports directly focused on determining the expression of LncRNAs in CF [14, 15]. Both of these studies have identified LncRNA expression profiles in primary CF cells type different from those used in our study. McKiernan et al., [14] analyzed LncRNAs in cells isolated from bronchial brushings of CF patients and non-CF controls and have identified XIST and TLR8 to be differentially expressed in these samples. However, we did not find significant alteration in the expression of these LncRNAs in the CF lung tissues, which were isolated from CF patients undergoing lung transplant. The source of the samples could be contributing to these differences. Moreover, the study by Balloy et al. [15] also used different cell types, primary CF bronchial epithelial cells which were infected with Pseudomonas aeruginosa, and have identified the differential expression of distinct non-coding transcripts.

Analyses of the corresponding coding mRNA expressions yields interesting data. There is a total of 9 mRNAs which are commonly altered in expression in the CF airway and parenchyma tissues (Fig. 2b). The proteins coded by the two up-regulated mRNAs, RP11-105C20.2 and LOC90925, have not been characterized. However, the 7 mRNA transcripts that are down-regulated in the CF lung tissues (both airway and parenchyma) are interestingly associated with CF disease phenotype. GPNMB levels have been reported to be markedly increased in serum and circulating leukocytes from CF patients compared to healthy controls [17]. On the contrary we find decreased levels of the GPNMB transcript in CF lung tissues. Also, interestingly CFTR has been reported to regulate GnRH secretion and thereby regulate sexual maturation and infertility in CF women [18]. MARCH9, an E3 ubiquitin ligase, is down-regulated in both CF airway and parenchyma tissues. Its activity is similar to that of MARCH2, which promotes ubiquitination and subsequent lysosomal degradation of mature CFTR [19]. TPPP has been shown to reduce microtubule formation rates and decrease microtubule acetylation, replicating the CF cellular condition [20]. GLOD5 is predicted to interact with HDAC6, and interestingly depletion of HDAC6 has been shown to rescue CF disease phenotype in a CF mouse model [21]. KLK12 expression has been shown to be suppressed in CF cells compared to control cells 2 h post-infection with Pseudomonas aeruginosa [15]. TIMP4 restricts fibrosis by regulating ECM deposition (through TGF-beta signaling pathway) and also causes restriction of inflammation [22]. Consistently, TIMP4 expression is reduced in CF lung tissues.

In silico analyses of these differentially expressed LncRNAs and mRNAs by GO indicates impact on relevant CF disease-specific pathways and cellular processes. The biological processes impacted include Toll-like receptor signaling, oxidation-reduction process, regulation of IL-6 production, regulation of cytosolic Calcium ions, G-protein signaling, etc. Consistently, cellular processes and molecular functions impacted include, vesicle transport, membrane function, LPS-binding, cell adhesion, oxidoreductase activity, transcriptional and translation activity, etc. These bioinformatic analyses support the importance of the array data.

Subsequently, we analyzed the expression of the top 10 up-regulated and top-10 down-regulated LncRNAs in CF airway tissues compared to control tissues (see Table 2) in CF airway epithelial cell line using TaqMan assay specific for each of these LncRNAs. Due to lack of sufficient RNA samples from CF tissues and also since in vitro cell culture systems will be employed for further mechanistic studies, we have analyzed the expression of selected LncRNAs in CFBE41o- cell lines as well as in respective control cell line with WT-CFTR. The LncRNAs which exhibit similar expression trends in CF epithelial cell lines as compared with those observed in the microarray data obtained from CF lung tissues, include RN7SKP237, ARHGAP19-SLIT1, AC003092.1 and LOC101928516 differentially expressed in the CF airway tissues; LINC01023 which is down-regulated in CF lung parenchyma tissues; and LINC01558 and LOC100133050, that exhibit similar expression trend both in the CF airway and parenchyma tissues.

The biological functions of the LncRNAs, RN7SKP237 and ARHGAP19-SLIT1, are poorly understood. RN7SKP237 is a pseudogene, which belongs to the anti-sense family of LncRNAs, and is associated with the microsomal glutathione transferase gene (MGST2). Interestingly, MGST2 protein catalyzes the biogenesis of leukotriene C4, which is a potent pro-inflammatory mediator of the pathophysiology of CF [23]. ARHGAP19-SLIT1 is a read through transcript between neighboring Rho-GTPase-activating protein 19 and slit homolog1 (SLIT1) and is a target of non-sense-mediated decay (NMD). This LncRNA spans two genes, ARHGAP19 and SLIT1. Interestingly, the coding transcript ARHGAP19 is a negative regulator of Rho GTPases, which are involved in cell migration, proliferation, and differentiation, actin remodeling, and G1 cell cycle progression [24], again relevant to CF disease. AC003092.1 has been shown to induce apoptosis in Glioblastoma (GB) by increased sensitivity to chemotherapy and thereby a potential therapeutic target for GB [25]. LINC01023 is a regulator of the IGF1R/Akt pathway in glioma [26], again a signaling pathway that has been shown to regulate CF lung disease. LOC101928516, LINC01558 and LOC100133050 have not been characterized.

The differentially expressed LncRNAs in CF may play important roles in the pathophysiology of CF lung disease. The LncRNA expression profile and the corresponding mRNA expression implicate CF relevant biological processes as affected in the CF disease tissues compared to control tissues. Further studies are being directed towards understanding the role of these LncRNAs in the disease phenotype of CF, including lung inflammation and mutant CFTR function. These mechanisms will serve as paradigms for similar complex processes and are expected to lead to the development of novel therapeutic targets in CF and other pulmonary disorders, such as COPD and Asthma.

Availability of data and materials

The datasets generated and/or analyzed during the current study are not publicly available due [Ongoing research with these] but are available from the corresponding author on reasonable request.

Abbreviations

CF:

Cystic Fibrosis

CFTR:

Cystic Fibrosis Transmembrane Conductance Regulator

GO:

Gene Ontology

LincRNA:

long intergenic non-coding RNA

LncRNA:

long non-coding RNA

References

  1. 1.

    Hon CC, Ramilowski JA, Harshbarger J, Bertin N, Rackham OJ, Gough J, Denisenko E, Schmeier S, Poulsen TM, Severin J, et al. An atlas of human long non-coding RNAs with accurate 5′ ends. Nature. 2017;543:199–204.

    CAS  Article  Google Scholar 

  2. 2.

    Derrien T, Johnson R, Bussotti G, Tanzer A, Djebali S, Tilgner H, Guernec G, Martin D, Merkel A, Knowles DG, et al. The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression. Genome Res. 2012;22:1775–89.

    CAS  Article  Google Scholar 

  3. 3.

    Wang KC, Chang HY. Molecular mechanisms of long noncoding RNAs. Mol Cell. 2011;43:904–14.

    CAS  Article  Google Scholar 

  4. 4.

    Bonasio R, Shiekhattar R. Regulation of transcription by long noncoding RNAs. Annu Rev Genet. 2014;48:433–55.

    CAS  Article  Google Scholar 

  5. 5.

    Ma L, Bajic VB, Zhang Z. On the classification of long non-coding RNAs. RNA Biol. 2013;10:925–33.

    PubMed  Google Scholar 

  6. 6.

    St Laurent G, Wahlestedt C, Kapranov P. The landscape of long noncoding RNA classification. Trends Genet. 2015;31:239–51.

    CAS  Article  Google Scholar 

  7. 7.

    Collier SP, Henderson MA, Tossberg JT, Aune TM. Regulation of the Th1 genomic locus from Ifng through Tmevpg1 by T-bet. J Immunol. 2014;193:3959–65.

    CAS  Article  Google Scholar 

  8. 8.

    Djebali S, Davis CA, Merkel A, Dobin A, Lassmann T, Mortazavi A, Tanzer A, Lagarde J, Lin W, Schlesinger F, et al. Landscape of transcription in human cells. Nature. 2012;489:101–8.

    CAS  Article  Google Scholar 

  9. 9.

    Li J, Xuan Z, Liu C. Long non-coding RNAs and complex human diseases. Int J Mol Sci. 2013;14:18790–808.

    CAS  Article  Google Scholar 

  10. 10.

    Maass PG, Luft FC, Bahring S. Long non-coding RNA in health and disease. J Mol Med (Berl). 2014;92:337–46.

    CAS  Article  Google Scholar 

  11. 11.

    Mestdagh P, Vandesompele J, Brusselle G, Vermaelen K. Non-coding RNAs and respiratory disease. Thorax. 2015;70:388–90.

    Article  Google Scholar 

  12. 12.

    Booton R, Lindsay MA. Emerging role of MicroRNAs and long noncoding RNAs in respiratory disease. Chest. 2014;146:193–204.

    CAS  Article  Google Scholar 

  13. 13.

    Saayman SM, Ackley A, Burdach J, Clemson M, Gruenert DC, Tachikawa K, Chivukula P, Weinberg MS, Morris KV. Long non-coding RNA BGas regulates the cystic fibrosis Transmembrane conductance regulator. Mol Ther. 2016;24:1351–7.

    CAS  Article  Google Scholar 

  14. 14.

    McKiernan PJ, Molloy K, Cryan SA, McElvaney NG, Greene CM. Long noncoding RNA are aberrantly expressed in vivo in the cystic fibrosis bronchial epithelium. Int J Biochem Cell Biol. 2014;52:184–91.

    CAS  Article  Google Scholar 

  15. 15.

    Balloy V, Koshy R, Perra L, Corvol H, Chignard M, Guillot L, Scaria V. Bronchial epithelial cells from cystic fibrosis patients express a specific long non-coding RNA signature upon Pseudomonas aeruginosa infection. Front Cell Infect Microbiol. 2017;7:218.

    Article  Google Scholar 

  16. 16.

    Rojas M, Marie B, Vignaud JM, Martinet N, Siat J, Grosdidier G, Cascorbi I, Alexandrov K. High DNA damage by benzo [a] pyrene 7,8-diol-9,10-epoxide in bronchial epithelial cells from patients with lung cancer: comparison with lung parenchyma. Cancer Lett. 2004;207:157–63.

    CAS  Article  Google Scholar 

  17. 17.

    Bhattacharyya S, Feferman L, Sharma G, Tobacman JK. Increased GPNMB, phospho-ERK1/2, and MMP-9 in cystic fibrosis in association with reduced arylsulfatase B. Mol Genet Metab. 2018;124:168–75.

    CAS  Article  Google Scholar 

  18. 18.

    Hodges CA, Palmert MR, Drumm ML. Infertility in females with cystic fibrosis is multifactorial: evidence from mouse models. Endocrinology. 2008;149:2790–7.

    CAS  Article  Google Scholar 

  19. 19.

    Cheng J, Guggino W. Ubiquitination and degradation of CFTR by the E3 ubiquitin ligase MARCH2 through its association with adaptor proteins CAL and STX6. PLoS One. 2013;8:e68001.

    CAS  Article  Google Scholar 

  20. 20.

    Rymut SM, Ivy T, Corey DA, Cotton CU, Burgess JD, Kelley TJ. Role of exchange protein activated by cAMP 1 in regulating rates of microtubule formation in cystic fibrosis epithelial cells. Am J Respir Cell Mol Biol. 2015;53:853–62.

    CAS  Article  Google Scholar 

  21. 21.

    Rosenjack J, Hodges CA, Darrah RJ, Kelley TJ. HDAC6 depletion improves cystic fibrosis mouse airway responses to bacterial challenge. Sci Rep. 2019;9:10282.

    Article  Google Scholar 

  22. 22.

    Arpino V, Brock M, Gill SE. The role of TIMPs in regulation of extracellular matrix proteolysis. Matrix Biol. 2015;44-46:247–54.

    CAS  Article  Google Scholar 

  23. 23.

    Sampson AP, Spencer DA, Green CP, Piper PJ, Price JF. Leukotrienes in the sputum and urine of cystic fibrosis children. Br J Clin Pharmacol. 1990;30:861–9.

    CAS  Article  Google Scholar 

  24. 24.

    Lv L, Xu J, Zhao S, Chen C, Zhao X, Gu S, Ji C, Xie Y, Mao Y. Sequence analysis of a human RhoGAP domain-containing gene and characterization of its expression in human multiple tissues. DNA Seq. 2007;18:184–9.

    CAS  Article  Google Scholar 

  25. 25.

    Xu N, Liu B, Lian C, Doycheva DM, Fu Z, Liu Y, Zhou J, He Z, Yang Z, Huang Q, et al. Long noncoding RNA AC003092.1 promotes temozolomide chemosensitivity through miR-195/TFPI-2 signaling modulation in glioblastoma. Cell Death Dis. 2018;9:1139.

    Article  Google Scholar 

  26. 26.

    Yu M, Yu S, Gong W, Chen D, Guan J, Liu Y. Knockdown of linc01023 restrains glioma proliferation, migration and invasion by regulating IGF-1R/AKT pathway. J Cancer. 2019;10:2961–8.

    Article  Google Scholar 

Download references

Acknowledgements

The authors thank Dr. Michael M. Myerburg, Airway Cell and Tissues Core of the Research Development Program of the CF Foundation at the University of Pittsburgh for his assistance with lung tissue samples.

Disclaimer

The views expressed are those of the authors and do not reflect the official policy or position of the Uniformed Services University of the Health Sciences, the Department of the Defense, or the United States government.

Funding

This study was supported by NHLBI/USUHS-CHIRP funds 1I80VP000012 [RB]; NIH grants DK68196 and DK72506 [RAF]; and Cystic Fibrosis Foundation FRIZZE05XX [RAF].

Author information

Affiliations

Authors

Contributions

PK designed and performed the experiments, analyzed and interpreted the data, and contributed to writing the manuscript, CS performed experiments and analyzed data, KP designed the experiments, RAF designed the experiments and contributed to writing the manuscript, and RB designed the experiment, analyzed and interpreted the data, and was a major contributor in writing the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Roopa Biswas.

Ethics declarations

Ethics approval and consent to participate

The lung tissue samples were collected under a University of Pittsburgh Institutional Review Board approved protocol.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Kumar, P., Sen, C., Peters, K. et al. Comparative analyses of long non-coding RNA profiles in vivo in cystic fibrosis lung airway and parenchyma tissues. Respir Res 20, 284 (2019). https://doi.org/10.1186/s12931-019-1259-8

Download citation

Keywords

  • Cystic fibrosis
  • F508del-CFTR
  • long non-coding RNA
  • mRNA
  • lung disease