# Bronchial Epithelial Cells from Cystic Fibrosis Patients Express a Specific Long Non-coding RNA Signature upon Pseudomonas aeruginosa Infection

Viviane Balloy, Remya Koshy, Lea Perra, Harriet Corvol, Michel Chignard, Loïc Guillot, Vinod Scaria

<figure><img src="image_1.png" /></figure>

## To cite this version:

Viviane Balloy, Remya Koshy, Lea Perra, Harriet Corvol, Michel Chignard, et al., Bronchial Epithelial Cells from Cystic Fibrosis Patients Express a Specific Long Non-coding RNA Signature upon Pseudomonas aeruginosa Infection. Frontiers in Cellular and Infection Microbiology, 2017, 7, pp.218. 10.3389/fcimb.2017.00218. [hal-01541639](https://hal.sorbonne-universite.fr/hal-01541639v1)

HAL Id: hal-01541639

https://hal.sorbonne-universite.fr/hal-01541639v1

Submitted on 19 Jun 2017

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_3.png" /></figure>

Check for
updates

# Bronchial Epithelial Cells from Cystic Fibrosis Patients Express a Specific Long Non-coding RNA Signature upon *Pseudomonas aeruginosa* Infection

Viviane Balloy<fnref n="1" />, Remya Koshy², Lea Perra<fnref n="1" />, Harriet Corvol<fnref n="1" />,<fnref n="3" />, Michel Chignard<fnref n="1" />, Loïc Guillot<fnref n="1" />*† and Vinod Scaria², 4*†

*Pseudomonas aeruginosa* (Pa) is the leading cause of chronic lung infection in Cystic Fibrosis (CF) patients. It is well recognized that CF epithelial cells fail to develop an appropriate response to infection, allowing bacterial colonization and a chronic inflammatory response. Since long non-coding RNAs (lncRNAs), are known to play a key role in regulating mammalian innate immune response, we hypothesized that CF cells exposed to Pa could express a specific lncRNA signature responsible of the maladaptive CF response. We analyzed transcriptomic datasets to compare the expression profiles of lncRNAs in primary CF and non-CF epithelial cells infected with Pa at 0, 2, 4, and 6 h of infection. Our analysis identified temporal expression signatures of 25, 73, 15, and 26 lncRNA transcripts differentially expressed at 0, 2, 4, and 6 h post-infection respectively, between CF and non-CF cells. In addition, we identified profiles specific to CF and non-CF cells. The differential expression of two candidate lncRNAs were independently validated using real-time PCR. We identified a specific CF signature of lncRNA expression in a context of Pa infection that could potentially play a role in the maladaptive immune response of CF patients.

**Keywords:** cystic fibrosis, lung, epithelium, *Pseudomonas aeruginosa*, lncRNA

## OPEN ACCESS

**Edited by:**

Silvia Mercedes Uriarte,
University of Louisville, United States

**Reviewed by:**

Gregory Anderson,
Indiana University - Purdue University
Indianapolis, United States

Gill Diamond,
University of Florida, United States

*Correspondence:**

Loïc Guillot
loic.guillot@inserm.fr
Vinod Scaria
vinods@igib.in

**Received:** 02 February 2017

**Accepted:** 11 May 2017

**Published:** 29 May 2017

**Citation:**

Balloy V, Koshy R, Perra L, Corvol H, Chignard M, Guillot L and Scaria V (2017) Bronchial Epithelial Cells from Cystic Fibrosis Patients Express a Specific Long Non-coding RNA Signature upon Pseudomonas aeruginosa Infection. *Front. Cell. Infect. Microbiol.* 7:218. doi: 10.3389/fcimb.2017.00218

## INTRODUCTION

Cystic fibrosis (CF) is an autosomal recessive genetic disorder affecting one in 3,500 live births, in Caucasian populations. This disease is due to mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene encoding the transporter channel for chloride and bicarbonate ions in epithelial cells. In airways, a deficient CFTR results in excessive secretion of abnormally thick and viscous mucus with impairment of innate immune host defenses. This results in chronic infection and inflammation, which together acts to further deteriorate the lung function (Palmer and Whiteley, 2015). *Pseudomonas aeruginosa* (Pa) is the leading cause of chronic lung infections and is responsible for significant morbidity and mortality in CF patients. Chronic Pa infection is characterized by repeated episodes of pulmonary exacerbations, during which lung function declines abruptly (Elborn, 2016).

<footnote n="1">¹ Sorbonne Universités, UPMC Univ Paris 06, INSERM, Centre de Recherche Saint-Antoine (CRSA), Paris, France, ² GN Ramachandran Knowledge Center for Genome Informatics, CSIR Institute of Genomics and Integrative Biology, Delhi, India,</footnote>

<footnote n="3">³ Pneumologie Pédiatrique, AP-HP, Hôpital Trousseau, Paris, France, ⁴ CSIR Institute of Genomics and Integrative Biology, Academy of Scientific and Innovative Research, Delhi, India</footnote>

<footnote>† These authors have contributed
equally to this work.</footnote>
Innate immunity system provides immediate defense against microbial invasion by recruiting immune cells and driving inflammation signals to sites of infection. Response of pulmonary epithelial cells to external stimuli requires tight regulation of genes involved in the innate immune response. This fine-tuned regulation, at the transcriptional and post-transcriptional levels, is necessary to prevent detrimental effects of uncontrolled cell activation. In a previous study, we compared transcriptomes of CF and non-CF bronchial epithelial cells during *Pa* infection (Balloy et al., 2015). We observed different profiles with genes up- and down-regulated in CF compared to non-CF cells that could illuminate the inflammatory response or be involved in the fighting of pathogens as well as in the decreased protection of tissue integrity.

Non-coding RNAs (ncRNAs) have recently emerged as important regulators of gene expression. Long ncRNAs (lncRNAs) larger than 200 nucleotides, represent the largest class of the ncRNA molecules. Although the current strongest evidence supports their role in cancer (Tsai et al., 2011), lncRNAs play key roles in regulating mammalian innate immune responses as well (Sigdel et al., 2015). Depending on their localization, they act on gene expression through interactions with DNA, RNA or proteins (Rinn and Chang, 2012) through various mechanisms including chromatin remodeling, epigenetic regulation, transcription, mRNA splicing, RNA decay and enhancer functions (Li et al., 2016).

It has been previously shown that lncRNA expression profile is altered in inflammatory lung diseases such as chronic obstructive pulmonary diseases (COPD), asthma and acute lung injury (ALI) (Xie and Liu, 2015). In CF, the description and the role of lncRNAs is beginning to be unraveled. However, *CFTR* expression and chloride ion function was recently shown to be regulated by lncRNA (McKiernan et al., 2014). A total of 1,063 lncRNAs with differential expression were identified from CF bronchial brush sampling compared to non-CF (Saayman et al., 2016). In particular, XIST and TLR8-AS1 were confirmed to be respectively up- and down-regulated in CF compared to non-CF samples. However, in this study, the influence of infection was not addressed, although there are increased evidences of lncRNA involvement in infectious diseases (Scaria and Pasha, 2012). Reports have revealed roles for lncRNAs as regulators of antimicrobial functions such as the lncRNA NEAT1 which directly regulates HIV replication (Zhang et al., 2013) or the lncRNA NRAV that functions as a negative regulator of the host antiviral immunity by repressing the expression of IFN-stimulated genes through strict control of the transcription rate (Ouyang et al., 2014).

We therefore investigated the differential expression of lncRNAs in CF and non-CF epithelial cells during infection by *Pa* that may partly mimic and possibly explain the differential mRNA profiles between the two cell types observed previously (Balloy et al., 2015). In the present study, we re-analyzed the RNA-seq data with this perspective and evaluated the expression profiles of lncRNAs in infected CF and non-CF epithelial cells at the aim of identifying a specific CF signature for lncRNA expression in the presence of *Pa* infection.

# METHODS

## Dataset

Raw datasets publicly available (European Nucleotide Archive (ENA), primary accession number PRJEB9292) were from a transcriptional analysis designed to compare the molecular determinants of CF and non-CF epithelial cells to *Pa* infection. The PAK strain of *Pa* (given by Stephen Lowry) is relevant for pulmonary infections in the context of CF Like most of clinical strains responsible for CF infection, it expresses a full complement of virulence factors, including pili, flagella, type II secreted enzymes, type III secreted exo-enzymes S, T, and Y, exotoxin A, elastase and phospholipase. The data were generated from human bronchial epithelial cells from four CF patients homozygous for the p.F508del mutation and four healthy donors (non-CF), each infected with *Pa* at 0.25 Multiplicity of Infection (MOI) for 0, 2, 4, and 6 h. The data with an average of over 150 million reads with each data set were previously published (Balloy et al., 2015).

## Reads Alignment, Mapping, and Assembly

Read quality was checked with FastQC program, a quality control tool for high throughput sequence data (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and reads were further aligned to the hg38 build of the human genome (UCSC Genome Browser). The splice-aware aligner Tophat (Mortazavi et al., 2008) (v2.0.11) which is integrated with Bowtie2, a fast and sensitive gapped read aligner, was used for aligning the sequence reads. One dataset was excluded from further analysis after observing only less than 85% of read alignment. Datasets, raw reads and the proportion of aligned reads are summarized in Supplementary Table 1. The expression levels across the transcripts and genes were based on GENCODE v23 (http://www.GENCODEgenes.org/). Briefly, this annotation set comprised of a total of 198,619 transcripts in the human genome, including 79,795 protein-coding and 118,824 non-coding transcripts, among which 27,817 lncRNAs. The expression levels were computed as fragments per kilobase of exon per million (FPKM) units and the differential expression between the different sets were computed using *cuffdiff*.

## Identifying Potential Differentially Expressed lncRNA Transcripts

We used a nominal cutoff of 1 FPKM to filter genes and transcripts. The Supplementary Table 2 lists out the number of genes and transcripts covered by each sample within the range of FPKM ≥ 1, 5, and 10. The differentially expressed genes and transcripts were filtered using stringent cutoffs. We considered lncRNAs with FPKM ≥ 1 with a log2-fold change ≤ -1 or ≥1 and p ≤ 0.05 as differentially expressed.

We performed two orthologous analyses across datasets. We first determined genes and transcripts which were differentially up- or down-regulated at each infection time points between CF and non-CF epithelial cells. Then, we examined the lncRNA transcripts differentially expressed between the 0 h time point and 2, 4, and 6 h post-infection in both CF and non-CF groups analyzed separately. From this, we sorted out the differentially
expressed lncRNA transcripts which were found only in CF or in non-CF cells. The overlapped transcripts of the first set comprising the differentially expressed lncRNA transcripts of CF/non-CF analysis and second set comprising the differentially expressed transcripts at 2, 4, and 6 h post-infection, identified only in CF cells, and gave us the potential lncRNA transcripts with potential for association with *Pa* infection in CF epithelial cells.

## Validation of Expression Level Using qRT-PCR

Quantitative PCR was performed using an ABI StepOnePlus™ (Applied Biosystems, Carlsbad, CA, USA) and TaqMan technology. TaqMan probes used (Applied Biosystems) were XIST (Hs01079824), MEG9 (Hs01593046), BLACAT1 (Hs03839366), and GAPDH (Hs02758991). Relative quantifications of the lncRNA level were carried out using the 2⁻⁰⁻⁰ΔΔCt method, and normalized with respect to GAPDH and for the expression levels to respective non-infected cells (T0). Each sample was assessed in triplicates to ensure experiment quality.

## Statistical Analysis

Data were described as mean ± SEM. Between-group differences were tested using the paired *t*-test. Values of *p* lower than 0.05 were considered significant; in the figures, statistically significant differences with *p* < 0.05 (*), *p* < 0.01 (**), and *p* < 0.001 (***) are indicated.

## Correlation Analysis of Gene Expression Data

Correlation analysis was carried out using Pearson's algorithm with the help of R package; rcorr. We observed the mRNA transcripts positively or negatively correlated with the overlapped lncRNAs from two analyses. Transcripts with an empirical r-score greater than 0.85 were considered to be positively regulated whereas transcripts with r-score less than -0.85 were considered to be negatively regulated and *p*-value less than 0.01 was considered as significant. We additionally performed enrichment analysis for the protein coding genes co-regulated with the lncRNAs using DAVID functional annotation tool.

# RESULTS

## Identification of Temporal Infection Signatures from RNA-Seq Data

For identifying differentially expressed genes and transcripts in CF compared to non-CF epithelial cells, we used *cuffdiff* and the reference gene and transcript annotations from the Gencode version 23 for guiding the assembly. The *cuffdiff* output provided FPKM values, log2-fold change value of CF upon non-CF for every time points studied. The numbers of genes differentially expressed at 0, 2, 4, and 6 h are found to be 297, 958, 281, and 439 respectively and a total of 336, 958, 336, and 518 transcripts respectively (Table 1). Their FPKM values and fold changes are summarized in Supplementary Table 3. The lncRNA annotation files were separately retrieved from the Gencode database to extract all the significantly differentially expressed lncRNA transcripts. We identified a total of 25, 73, 15, and 26 differentially expressed lncRNAs, corresponding to the 0, 2, 4, and 6 h time points. The unique number of lncRNA transcripts differentially expressed throughout the dataset was found to be 108. These latter are shown in a Venn diagram (Figure 1) to visualize shared lncRNAs differentially expressed between CF

<figure><img src="image_4.png" /><figcaption>FIGURE 1 | Venn Diagram of differentially expressed lncRNA transcripts in CF vs. non-CF cells at each time point of infection.</figcaption></figure>

<table><caption>TABLE 1 | Summary of genes and transcripts differentially expressed in CF vs. non-CF cells at each time point of infection.</caption><thead><tr><th>Differentially expressed genes/transcripts</th><th>0 h (CF/non-CF)</th><th>2 h (CF/non-CF)</th><th>4 h (CF/non-CF)</th><th>6 h (CF/non-CF)</th></tr></thead><tbody><tr><td>Genes</td><td>297</td><td>958</td><td>281</td><td>439</td></tr><tr><td>Transcripts</td><td>336</td><td>958</td><td>336</td><td>518</td></tr><tr><td>lncRNA transcripts</td><td>25</td><td>73</td><td>15</td><td>26</td></tr><tr><td>lncRNA genes</td><td>21</td><td>73</td><td>15</td><td>22</td></tr><tr><td>Unique lncRNA transcripts</td><td></td><td></td><td>108</td><td></td></tr><tr><td>Unique lncRNA genes</td><td></td><td></td><td>104</td><td></td></tr></tbody></table>
and non-CF infected cells. Out of the 108 lncRNA transcripts, a total of 12 lncRNAs are differentially expressed in at least 2 time points. Table 2 lists these lncRNAs, their cognate gene locus and functions as annotated from literature. Three of the 12 lncRNAs have been previously shown to be expressed in the human lung. The lncRNA RP11-44F21.5 has been previously suggested to be down-regulated in squamous cell carcinoma of the lung in the following analysis: https://genevisible.com/perturbations/HS/Gene%20Symbol/RP11-44F21.5. The two other lncRNAs,

LINC00704 and LINC00992, have been shown to be expressed in lungs, kidneys and salivary glands (Stelzer et al., 2016).

## Highlighting Distinct Specific Infection Signatures in Non-CF and CF Cells

To identify the lncRNA transcripts specific of Pa infection, we used *cuffdiff* to compare the expression values for each time point against 0 h for CF and non-CF datasets separately. We retrieved the differentially expressed transcripts using

<table><caption>TABLE 2 | Summary of lncRNA genes differentially expressed in CF vs. non-CF cells for at least 2 of the time points of infection.</caption><thead><tr><th>Conditions differentially expressed</th><th>IncRNA</th><th>Ensembl ID</th><th>Function</th></tr></thead><tbody><tr><td></td><td>CTD-2619J13.13</td><td>ENST00000596379.1<br>ENSG00000268307.1</td><td>CTD-2619J13.13 is a processed transcript type of long coding RNA. It belongs to chromosome 19 on the forward strand comprising 459 bp.</td></tr><tr><td></td><td>LINC00862</td><td>ENST00000367356.1<br>ENSG00000203721.5</td><td>LINC00862 is a long Intergenic Non-Protein Coding RNA 862. It is also known as Small Integral Membrane Protein 16 or C1orf98. It belongs to chromosome 12, comprising 31,811 bp on reverse strand.</td></tr><tr><td></td><td>CMB9-55F22.1</td><td>ENST00000623846.1<br>ENSG00000279672.1</td><td>CMB9-55F22.1 is a known pseudogene which is annotated as long coding RNA. It is found to have expressed in acute monoblastic/monocytic leukemia cell lines. It belongs to chromosome 11, encompassing 1,139 bp on forward strand</td></tr><tr><td></td><td>POLR2J4</td><td>ENST00000427076.5<br>ENSG00000214783.9</td><td>POLR2J4 is expanded as Polymerase (RNA) II (DNA Directed) Polypeptide J4. It is a pseudogene, which is annotated as long coding RNA. It is also known as RPB11-Phi. It is found to be on chromosome 12 carrying 78,300 bps on reverse strand. It is studied to be expressed in majority all tissue types.</td></tr><tr><td></td><td>WFDC21P</td><td>ENST00000566140.5<br>ENSG00000261040.6</td><td>WFDC21P is expanded as WAP Four-Disulfide Core Domain 21 which is a pseudogene type of long coding RNA. It is also known as IncDC. The STAT3-Binding Long Non-coding RNA Inc-DC Controls Human Dendritic Cell Differentiation. It belongs to chromosome 17 containing 8,320 bp on reverse strand</td></tr><tr><td></td><td>RP11-442H21.2</td><td>ENST00000491934.2<br>ENSG00000269926.1</td><td>RP11-442H21.2 is a known antisense long non-coding RNA. It is found to be expressed in bone marrow and cord blood. It belongs to chromosome 10 comprising 1,066 bp on reverse strand.</td></tr><tr><td></td><td>RP11-890B15.2</td><td>ENST00000533381.2<br>ENSG00000254842.6</td><td>RP11-890B15.2 is a known long intergenic non-coding RNA. It has found to have some regulatory activities like CTCF binding, promoter seq binding. It belongs to chromosome 11 comprising 2,203 bps on reverse strand.</td></tr><tr><td></td><td>RP11-445P17.8</td><td>ENST00000449457.1<br>ENSG00000224034.1</td><td>RP11-445P17.8 is a known lincRNA. It has found to be expressed in gallbladder and acting as an enhancer. It belongs to chromosome 10, 5,204 bp long on reverse strand.</td></tr><tr><td></td><td>RP11-44F21.5</td><td>ENST00000561705.1<br>ENSG00000260265.1</td><td>RP11-44F21.5 is a known lincRNA. It has found to be expressed in bladder and brain. The lncRNA is studied to be downregulated in pancreatic cancer studies and lung squamous carcinoma with normal non-cancerous tissues while found upregulated between allergic house dust mite and non-allergic dust mite. This belongs to chromosome 4 comprising 3,016 bps.</td></tr><tr><td></td><td>RP11-495P10.8</td><td>ENST00000434245.2<br>ENSG00000231196.3</td><td>RP11-495P10.8 is a long intergenic non-coding RNA which has shown role in epigenetic modifications in liver and brain cells. The overall mean expression of RP11-495P10.8 from different datasets is 0.20 RPM, ranging from as low as 0.00 RPM to 7.29 RPM. It is studied to be expressed more in erythrocytic leukemia. It belongs to chromosome 1 spanning 4,971 bp on reverse strand.</td></tr><tr><td></td><td>LINC00704</td><td>ENST00000430998.6<br>ENSG00000231298.6</td><td>LINC00704 is an RNA Gene, and is affiliated with the non-coding RNA class. It is found to have expressed in internal organs like lungs and kidneys and secretory tissue like salivary and pituitary gland. It belongs to chromosome 10 comprising 27,970 bps.</td></tr><tr><td></td><td>LINC00992</td><td>ENST00000504107.1<br>ENSG00000248663.6</td><td>LINC00992 expanded as Long Intergenic Non-Protein Coding RNA 992 is an RNA Gene, and is affiliated with the non-coding RNA class. It is also known as CTC-504A5.1. It is studied to have expressed in internal organs like lungs and kidneys and secretory tissue like salivary and pituitary gland. It belongs to chromosome 5 spanning the size 164,236 bps on forward strand.</td></tr></tbody></table>
<table><caption>TABLE 3 | Genes and lncRNAs differentially expressed in CF and non-CF cells at different time points of infection compared to the non-infected status (0 h).</caption><thead><tr><th rowspan="2">Differentially expressed genes/transcripts</th><th colspan="2">2 h vs. 0 h</th><th colspan="2">4 h vs. 0 h</th><th colspan="2">6 h vs. 0 h</th></tr><tr><th>CF</th><th>Non-CF</th><th>CF</th><th>Non-CF</th><th>CF</th><th>Non-CF</th></tr></thead><tbody><tr><td>Genes</td><td>72</td><td>587</td><td>99</td><td>75</td><td>272</td><td>140</td></tr><tr><td>Transcripts</td><td>72</td><td>603</td><td>102</td><td>76</td><td>305</td><td>146</td></tr><tr><td>lncRNA transcripts</td><td>2</td><td>53</td><td>1</td><td>7</td><td>17</td><td>5</td></tr><tr><td>lncRNA Genes</td><td>2</td><td>53</td><td>1</td><td>7</td><td>17</td><td>5</td></tr><tr><td>Unique lncRNA transcripts</td><td colspan="6">80</td></tr><tr><td>Unique lncRNA genes</td><td colspan="6">79</td></tr></tbody></table>

the same filters as mentioned previously. The total number of genes or transcripts, differentially expressed were found comparatively less in CF than in non-CF samples at 2 h vs. 0 h (Table 3). We further analyzed the common lncRNA transcripts differentially expressed at each time point in CF and non-CF datasets. Except for 2 up-regulated lncRNAs at 6 h vs. 0 h, we obtained two distinct profiles of lncRNAs as none of lncRNAs differentially expressed during infection in non-CF are common to those differentially expressed in CF cells (Figure 2). The FPKM values and fold changes of the up- and down-regulated transcripts are summarized in Supplementary Tables 4, 5, respectively.

For finding out the lncRNA signature specific to infected CF cells in the one hand and to non-CF cells in the other hand, we summarized overlaps between differentially expressed (both up and down) transcripts in the CF sets (Figure 3A) and in the non-CF sets (Figure 3B) respectively for each of the respective time points against non-infected time point (0 h). We found that 62 lncRNA transcripts (62 up-regulated) are specifically regulated in *Pa* infected non-CF epithelial cells, whereas 17 lncRNA transcripts (2 up- and 15 down-regulated) are specific of the CF epithelial cell response to the infection. The FPKM values and fold changes of the up-regulated transcripts in non-CF cells and the up- and down-regulated in CF cells are summarized in Supplementary Tables 6–8, respectively.

Among the lncRNA transcripts specific signature of infected CF cells, we selected two of them already studied in the literature, MEG9 and BLACAT1, for validation of their expression by qPCR. Thus, the pattern of expression found by *in silico* analysis for MEG9 and BLACAT1 was confirmed with a significant down regulation at 2 and 6 h post-infection, respectively (Figure 4). One lncRNA transcript, XIST, identified in the infected non-CF signature and already described by Greene's team comparing CF and non-CF samples, was measured by qPCR. Validation of its expression was not as clear as for MEG9 and BLACAT1. Indeed, we observed varying levels of expression between epithelial cells derived from male and female, independently of CFTR genotype, i.e., very low or no expression in male cells (Ct around 36 or undetermined) whereas high expression in female cells (Ct around 18.5) (data not shown).

## Proteins Coding Transcripts Correlated with BLACAT1 and MEG9 Expression

The Gene Ontology analysis did not reveal protein-coding genes significantly correlated with BLACAT1 expression, but it identified protein coding transcripts positively and negatively correlated with lncRNA MEG9. Among the genes positively correlated, we observed 62.3% of proteins engaged in protein binding, 18.5% in positive regulation of cellular process and around 10% in extracellular matrix. The genes negatively correlated (17.4%) with lncRNA MEG9 were involved in structural components of ribosome (Table 4).

## DISCUSSION

The airway epithelium is the first site in contact with inhaled pathogens. Its response to bacterial invasions activates innate immune mechanisms through the involvement of various receptors. In CF patients, this normal process is disturbed leading to a decrease of the bacterial clearance resulting in chronic infection in lungs.

Over the years, lncRNAs have been recognized as major regulators of multiple cellular processes among which the regulation of immune effectors. They can interact with transcription factors which is the case of the lncRNA NRON that inhibits the translocation of NFAT in the nucleus by sequestering it in the cytosol and so negatively regulating T-cell activation (Willingham et al., 2005). Three lncRNAs (murine NeST, human THRIL, and NEAT1) have been previously shown to regulate the innate immune response by modulating the transcription of cytokines such as IFN-γ, TNF-α, and IL-8 (Cullen, 2013; Gomez et al., 2013; Imamura et al., 2014; Li et al., 2014). It is also known that toll-like receptor (TLR) signaling can activate endogenous feedback-regulation networks to limit the potentially damaging effects of an excessive inflammation. This is the case with the lncRNA THRIL which helps in down-regulation to restrain TLR-induced gene activation (Li et al., 2014).

lncRNA expression is tissue-specific indicating a tight regulation that could be affected by *Pa* infection in CF patients. For this reason, we decided to explore the expression of lncRNAs in the infectious context of respiratory epithelial cells that are key cells involved in CF. CF vs. non-CF infected cells study brought 108 unique lncRNAs as differentially up- or down-regulated with 12 lncRNAs which are expressed in at least two time points. Among them, the lncRNAs, LINC00862, and CTD-2619J13, are differentially expressed at all four time points in spite of bacterial invasion. As such, they could be two potential signatures for CF. Most of the 108 unique lncRNAs was differentially expressed at 2 h, almost immediately after bacterial attack. According to Ensembl and lncRNAb, some of these differentially expressed lncRNAs are involved in regulatory activities namely CTCF binding, promoter sequence binding transcripts and few were found differentially expressed in different tissues such as in lungs and in erythrocytes.

To pinpoint a specific signature related to infection in the CF vs. non-CF context, we took in account lncRNA transcripts specific to *Pa* infection expressed only in CF in one hand and in
<figure><img src="image_5.png" /><figcaption>**FIGURE 2** | Overlaps between lncRNA transcripts differentially expressed at different time points in CF and non-CF datasets. Red circles represent down-regulated transcripts and green circles represent up-regulated transcripts. (A) 2 h in comparison with 0 h time point, (B) 4 h in comparison with 0 h time point, (C) 6 h in comparison with 0 h time point.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**FIGURE 3** | Summary of overlaps between differentially expressed transcripts only in the CF set (A) and only in non-CF set (B) at each of the respective time points against non-infected time point (0 h). Red circles represent down-regulated transcripts and green circles represent up-regulated transcripts. No transcripts were found down-regulated exclusively in non-CF set.</figcaption></figure>

non-CF cells in other hand from each time point (i.e., 2, 4, and 6 h) and compared with their expression in uninfected cells.

Numbers of lncRNAs differentially expressed for CF and non-CF cells are totally different with 17 and 62 lncRNAs, respectively. Among the lncRNAs belonging to the specific signature of *Pa* infected non-CF cells; we observed a well annotated lncRNA, XIST. As described by Greene's team ⁹, we detected the lncRNA XIST more expressed in CF than in non-CF cells at the non-infected status (Fold change >4). Moreover, we observed that XIST expression is up-regulated in non-CF cells during infection at 4 and 6 h (FC > 2). Nonetheless, the validation of XIST by qPCR does not fit with the *in silico* analysis. Indeed, we observed differences of expression between cells from male and female donors but no differences of expression between CF and non-CF epithelial cells. This result is consistent with the literature as XIST is a lncRNA required for the establishment of

X-chromosome silencing in placental mammals. X inactivation is an early developmental process in mammalian females that transcriptionally silences one of the pair of X chromosomes, thus providing dosage equivalence between males and females (Furlan and Rougeulle, 2016). Due to the differences in the gender of donors, taken individually, values did not show the variations averaged *in silico* analysis. Other transcripts, ENST00000452120.6 and ENST00000398461.5 were found in the CF cells specific signatures. These transcripts are annotated as MEG3 (Maternally Expressed 3) which is a maternally imprinted gene with all spliced transcripts annotated as lncRNAs and is a well-known lncRNA tumor suppressor. MEG3 has been previously reported to regulate the TGF-β pathway genes by forming RNA-DNA triplex structures (Mondal et al., 2015). This pathway is frequently regulated during tissue injury and repair and the up-regulation of TGF-β is involved in lung fibrogenic
<figure><img src="image_7.png" /><figcaption>**FIGURE 4** | Validation by qRT-PCR of BLACAT1 and MEG9 expressions in bronchial primary epithelial CF cells infected with *Pa*. RNA were purified in CF cells collected from 4 different patients and infected with *Pa*. LncRNA levels were measured at 0 (white bar) and 2 h (black bar) post-infection for MEG9 and 0 (white bar) and 6 h (black bar) post-infection for BLACAT1. The statistical analysis consisted in a paired *t*-test. \*P < 0.05 and \*\*P < 0.001.</figcaption></figure>

<table><caption>TABLE 4 | DAVID's GO terms for protein coding transcripts positively and negatively correlated with the lncRNA MEG9.</caption><thead><tr><th>Category</th><th>GO term</th><th>Count</th><th>Percentage</th><th>P-value</th><th>Benjamini corrected p-value</th></tr></thead><tbody><tr><td colspan="6"><strong>GENES POSITIVELY CORRELATED IN EXPRESSION</strong></td></tr><tr><td rowspan="2">Biological process</td><td>Regulation of localization</td><td>19</td><td>9.2</td><td>3.0E-4</td><td>4.6E-2</td></tr><tr><td>Positive regulation of cellular process</td><td>38</td><td>18.5</td><td>7.0E-4</td><td>5.30E-02</td></tr><tr><td rowspan="2">Cellular components</td><td>Extracellular matrix part</td><td>10</td><td>4.8</td><td>1.3E-5</td><td>4.4E-4</td></tr><tr><td>Extracellular matrix</td><td>14</td><td>6.8</td><td>3.1E-4</td><td>5.3E-3</td></tr><tr><td>Molecular function</td><td>Protein binding</td><td>129</td><td>62.3</td><td>1.60E-4</td><td>7.80E-3</td></tr><tr><td colspan="6"><strong>GENES NEGATIVELY CORRELATED IN EXPRESSION</strong></td></tr><tr><td>Molecular function</td><td>Structural constituent of ribosome</td><td>4</td><td>17.4</td><td>1.9E-3</td><td>3.10E-2</td></tr></tbody></table>

diseases. The expression of the lncRNA MEG3 in signature could be associated with the absence of injury after invasion of *Pa* in non-CF. A recent study in mice characterized MEG3 as a novel pulmonary inflammatory regulator of bacterial infection through miR-138 (Li and Wu, 2016).

The highlighting of a specific signature in CF epithelial cells in response to *Pa* infection revealed 17 lncRNA transcripts. Among them, 4 lncRNA transcripts MEG9, BLACAT1, RP11-477I4.4, and RP11-1334A24.5 overlapped with the 108 unique lncRNAs which are differentially expressed between CF and non-CF cells. RT-PCR gave a similar pattern of expression of MEG9 and BLACAT1 as from *in silico* analysis and validate the different signatures obtained with the computational analysis. ENST00000626538.1, otherwise known as lncRNA BLACAT1 (Bladder Cancer Associated Transcript 1), was found down-regulated in CF patients. Until now, BLACAT1, firstly characterized in bladder cancer, was also described to be involved in the development of gastric cancer. It has been shown that its depletion decreased the gastric cancer cell proliferation, motility and invasion (Hu et al., 2015). Since shedding and abnormal repair are characteristics of CF epithelial cells, it would be interesting to study whether their altered expression of BLACAT1 influences their repair, proliferation and migration. In our study, analysis did not reveal protein-coding gene significantly correlated with BLACAT1 expression. However,

as shown previously BLACAT1 could interact directly with other molecular partners than transcripts such as the proteins of the polycomb repressive complex 2 (PRC2) (He et al., 2013).

Concerning MEG9 (Maternally Expressed 9), a lncRNA belonging to the MEG family, it has been described to be up-regulated upon hypoxia in human endothelial cells and in a model of mouse ischemia (Voellenkle et al., 2016). In our study, MEG9 was found dramatically down-regulated 2 h post-infection compared to the non-infected status (0 h) in CF cells (Supplementary Table 8). From our correlation analysis of gene expression, MEG9 was found co-regulated positively with 226 and negatively with 25 protein coding genes. Gene ontology analysis showed that MEG9 was found positively co-regulated mostly with proteins engaged in protein binding (62.3%), in positive regulation of cellular process (18.5%) and few of them in extracellular matrix (Table 4), all processes involved in inflammation of the lung in CF (Cohen-Cumberknoh et al., 2013).

Our study showed that *Pa* infection altered the expression of lncRNAs in CF respiratory epithelial cells which could have a potential role in the innate immune system and thus could contribute to the maladaptive immune response observed in patients. We have experimentally validated the two potential transcripts namely MEG9 and BLACAT1 expression, which
were found differentially expressed in both CF/non-CF and infected/non-infected analyses. Further investigations either *in vitro* (respiratory epithelial cells in culture) or *in vivo* (study of conserved LncRNA (MEGs) in mice model of pulmonary infection) are necessary to understand the role of these lncRNAs in the innate immune response. The present study may contribute to provide further insight into the biological functions and molecular mechanisms of lncRNAs regulation in CF in which Pa infection plays a crucial role.

## AUTHOR CONTRIBUTIONS

Conceived and designed the experiments: LG, VB, and MC. Performed the biological experiments: LG, VB. Analyzed the data: RK, LG, VB, and VS. Wrote the paper: RK, VB, LG, MC, LP, HC, and VS.

## REFERENCES

Balloy, V., Varet, H., Dillies, M.-A., Proux, C., Jagla, B., Coppée, J.-Y., et al. (2015). Normal and cystic fibrosis human bronchial epithelial cells infected with *Pseudomonas aeruginosa* exhibit distinct gene activation patterns. *PLoS ONE* 10:e0140979. doi: 10.1371/journal.pone.0140979

Cohen-Cymberknoh, M., Kerem, E., Ferkol, T., and Elizur, A. (2013). Airway inflammation in cystic fibrosis: molecular mechanisms and clinical implications. *Thorax* 68, 1157-1162. doi: 10.1136/thoraxjnl-2013-203204

Cullen, B. R. (2013). Making a NeST for a persistent virus. *Cell Host Microbe* 13, 241-242. doi: 10.1016/j.chom.2013.02.016

Elborn, J. S. (2016). *Cystic Fibrosis*. Available online at: www.thelancet.com, 388.

Furlan, G., and Rougeulle, C. (2016). Function and evolution of the long noncoding RNA circuitry orchestrating X-chromosome inactivation in mammals. *Wiley Interdiscip. Rev. RNA* 7, 702-722. doi: 10.1002/wrn a.1359

Gomez, J. A., Wapinski, O. L., Yang, Y. W., Bureau, J.-F., Gopinath, S., Monack, D. M., et al. (2013). The NeST long ncRNA controls microbial susceptibility and epigenetic activation of the interferon-γ locus. *Cell* 152, 743-754. doi: 10.1016/j.cell.2013.01.015

He, W., Cai, Q., Sun, F., Zhong, G., Wang, P., Liu, H., et al. (2013). linc-UBC1 physically associates with polycomb repressive complex 2 (PRC2) and acts as a negative prognostic factor for lymph node metastasis and survival in bladder cancer. *Biochim. Biophys. Acta* 1832, 1528-1537. doi: 10.1016/j.bbadis.2013.05.010

Hu, Y., Pan, J., Wang, Y., Li, L., and Huang, Y. (2015). Long noncoding RNA linc-UBC1 is negative prognostic factor and exhibits tumor pro-oncogenic activity in gastric cancer. *Int. J. Clin. Exp. Pathol.* 8, 594-600.

Imamura, K., Imamachi, N., Akizuki, G., Kumakura, M., Kawaguchi, A., Nagata, K., et al. (2014). Long noncoding RNA NEAT1-dependent SFPQ relocation from promoter region to paraspeckle mediates IL8 expression upon immune stimuli. *Mol. Cell* 53, 393-406. doi: 10.1016/j.molcel.2014.01.009

Li, R., and Wu, M. (2016). Long noncoding RNA MEG3-4 regulates IL-1β during pulmonary bacterial infection by repressing microRNA-138. *J. Immunol.* 196(1 Suppl.), 127-19.

Li, R., Zhu, H., and Luo, Y. (2016). Understanding the functions of long non-coding RNAs through their higher-order structures. *Int. J. Mol. Sci.* 17:702. doi: 10.3390/ijms17050702

Li, Z., Chao, T.-C., Chang, K.-Y., Lin, N., Patil, V. S., Shimizu, C., et al. (2014). The long noncoding RNA THRIL regulates TNF expression through its interaction with hnRNPL. *Proc. Natl. Acad. Sci. U.S.A.* 111, 1002-1007. doi: 10.1073/pnas.1313768111

## ACKNOWLEDGMENTS

This work was supported by a grant from two French non-profit cystic fibrosis organizations, Vaincre la mucoviscidose and Association Gregory Lemarchal. The analysis was funded by CSIR, India through Grant BSC0123(GENCODE-C). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Authors also acknowledge Dr. Srinivasan Ramachandran for editorial help.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fcimb.2017.00218/full#supplementary-material

McKiernan, P. J., Molloy, K., Cryan, S. A., McElvaney, N. G., and Greene, C. M. (2014). Long noncoding RNA are aberrantly expressed *in vivo* in the cystic fibrosis bronchial epithelium. *Int. J. Biochem. Cell Biol.* 52, 184-191. doi: 10.1016/j.biocel.2014.02.022

Mondal, T., Subhash, S., Vaid, R., Enroth, S., Uday, S., Reinius, B., et al. (2015). MEG3 long noncoding RNA regulates the TGF-β pathway genes through formation of RNA-DNA triplex structures. *Nat. Commun.* 6:7743. doi: 10.1038/ncomms8743

Mortazavi, A., Williams, B. A., McCue, K., Schaeffer, L., Wold, B., Salzberg, S. L., et al. (2008). Mapping and quantifying mammalian transcriptomes by RNA-Seq. *Nat. Methods* 5, 621-628. doi: 10.1038/nmeth.1226

Ouyang, J., Zhu, X., Chen, Y., Wei, H., Chen, Q., Chi, X., et al. (2014). NRAV, a long noncoding RNA, modulates antiviral responses through suppression of interferon-stimulated gene transcription. *Cell Host Microbe* 16, 616-626. doi: 10.1016/j.chom.2014.10.001

Palmer, G. C., and Whiteley, M. (2015). Metabolism and pathogenicity of *Pseudomonas aeruginosa* infections in the lungs of individuals with cystic fibrosis. *Microbiol. Spectr.* 3. doi: 10.1128/microbiolspec.MBP-0003-2014

Rinn, J. L., and Chang, H. Y. (2012). Genome regulation by long noncoding RNAs. *Annu. Rev. Biochem.* 81, 145-166. doi: 10.1146/annurev-biochem-051410-092902

Saayman, S. M., Ackley, A., Burdach, J., Clemson, M., Gruenert, D. C., Tachikawa, K., et al. (2016). Long non-coding RNA BGAs regulates the cystic fibrosis transmembrane conductance regulator. *Mol. Ther.* 24, 1351-1357. doi: 10.1038/mt.2016.112

Scaria, V., and Pasha, A. (2012). Long non-coding RNAs in infection biology. *Front. Genet.* 3:308. doi: 10.3389/fgene.2012.00308

Sigdel, K. R., Cheng, A., Wang, Y., Duan, L., and Zhang, Y. (2015). The emerging functions of long noncoding RNA in immune cells: autoimmune diseases. *J. Immunol. Res.* 2015, 1-9. doi: 10.1155/2015/848790

Stelzer, G., Rosen, N., Plaschkes, I., Zimmerman, S., Twik, M., Fishilevich, S., et al. (2016). The genecards suite: from gene data mining to disease genome sequence analyses. *Curr. Protoc. Bioinform.* 54, 1.30.1-1.30.33. doi: 10.1002/cpbi.5

Tsai, M.-C., Spitale, R. C., and Chang, H. Y. (2011). Long intergenic noncoding RNAs: new links in cancer progression. *Cancer Res.* 71, 3-7. doi: 10.1158/0008-5472.CAN-10-2483

Voellenkle, C., Garcia-Manteiga, J. M., Pedrotti, S., Perfetti, A., De Toma, I., Da Silva, D., et al. (2016). Implication of long noncoding RNAs in the endothelial cell response to hypoxia revealed by RNA-sequencing. *Sci. Rep.* 6:24141. doi: 10.1038/srep24141
Willingham, A. T., Orth, A. P., Batalov, S., Peters, E. C., Wen, B. G., Aza-Blanc, P., et al. (2005). A strategy for probing the function of noncoding RNAs finds a repressor of NFAT. *Science* 309, 1570–1573. doi: 10.1126/science.1115901

Xie, N., and Liu, G. (2015). ncRNA-regulated immune response and its role in inflammatory lung diseases. *Am. J. Physiol. Lung Cell. Mol. Physiol.* 309, L1076–L1087. doi: 10.1152/ajplung.00286.2015

Zhang, Q., Chen, C.-Y., Yedavalli, V. S., and Jeang, K.-T. (2013). NEAT1 long noncoding RNA and paraspeckle bodies modulate HIV-1 posttranscriptional expression. *MBio* 4:e00596-12. doi: 10.1128/mbio.000596-12

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2017 Balloy, Koshy, Perra, Corvol, Chignard, Guillot and Scaria. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.