Mediators of Inflammation

Mediators of Inflammation / 2018 / Article
Special Issue

Airway Inflammatory/Immune Responses in COPD and Cystic Fibrosis

View this Special Issue

Research Article | Open Access

Volume 2018 |Article ID 7501851 | 14 pages | https://doi.org/10.1155/2018/7501851

Long Noncoding RNAs and mRNA Regulation in Peripheral Blood Mononuclear Cells of Patients with Chronic Obstructive Pulmonary Disease

Academic Editor: Catherine Greene
Received09 Aug 2017
Revised08 Nov 2017
Accepted27 Nov 2017
Published13 Mar 2018

Abstract

Background. Inflammation plays a pivotal role in the pathogenesis of chronic obstructive pulmonary disease (COPD). We evaluated the lncRNA and mRNA expression profile of peripheral blood mononuclear cells (PBMCs) from healthy nonsmokers, smokers without airflow limitation, and COPD patients. Methods. lncRNA and mRNA profiling of PBMCs from 17 smokers and 14 COPD subjects was detected by high-throughput microarray. The expression of dysregulated lncRNAs was validated by qPCR. The lncRNA targets in dysregulated mRNAs were predicted and the GO enrichment was analyzed. The regulatory role of lncRNA ENST00000502883.1 on CXCL16 expression and consequently the effect on PBMC recruitment were investigated by siRNA knockdown and chemotaxis analysis. Results. We identified 158 differentially expressed lncRNAs in PBMCs from COPD subjects compared with smokers. The dysregulated expression of 5 selected lncRNAs NR_026891.1 (FLJ10038), ENST00000502883.1 (RP11-499E18.1), HIT000648516, XR_429541.1, and ENST00000597550.1 (CTD-2245F17.3), was validated. The GO enrichment showed that leukocyte migration, immune response, and apoptosis are the main enriched processes that previously reported to be involved in the pathogenesis of COPD. The regulatory role of ENST00000502883.1 on CXCL16 expression and consequently the effect on PBMC recruitment was confirmed. Conclusion. This study may provide clues for further studies targeting lncRNAs to control inflammation in COPD.

1. Introduction

Chronic obstructive pulmonary disease (COPD) is characterized by progressive airflow obstruction that is only partly reversible, inflammation in the airways, and other systemic effects [1]. So far there is no satisfactory therapy to treat individuals once the disease is established. The mechanism of the disease generally involves aberrant, chronic inflammation coupled with the loss of lung structural cells due to heightened apoptosis.

Long noncoding RNAs (lncRNAs) are mRNA-like transcripts longer than 200 nucleotides without protein-coding functions [2]. According to the genomic proximity to protein-coding genes, lncRNA was generally classified into five types: (1) sense, (2) antisense, (3) bidirectional, (4) intronic, and (5) intergenic [3]. lncRNAs have been shown to play important roles in diverse biological and pathological processes and are dysregulated in various human diseases [4]. However, few studies concerning the dysregulation and regulatory role of lncRNA in lung diseases have been reported. Several hundred lncRNAs were identified in the developing mouse lung by sequencing polyadenylated RNAs in embryonic and adult lung tissue. NANCI and LL34 were found to regulate expression of hundreds of genes in mouse airway epithelial cell culture, but their role in lung diseases remains unclear [5]. A recent study showed that lncRNAs including RP11-46A10.4, LINC00883, BCYRN1, and LINC00882 act as miRNA “sponges” to regulate the growth of airway smooth muscle cells [6]. Bi et al. [7] analyzed lncRNA expression in the lung tissue of nonsmokers, smokers without COPD, and smokers with COPD; RNA44121|UCSC-2000-3182 and RNA43510|UCSC-1260-3754 were found the most over- and underregulated lncRNA, respectively.

Inflammation plays a pivotal role in the pathogenesis of COPD, where CD8+ T lymphocytes, neutrophils, and macrophages are the main type of immune cells of local inflammatory milieu of COPD [8]. Different immunoregulatory properties of T cells and monocytes have been demonstrated in COPD patients [9]. In the previous study, we have analyzed the expression profile of miRNAs as well as regulation network of dysregulated miRNAs and mRNAs in PBMCs of COPD patients [10]. Several studies on lncRNA expression of peripheral blood mononuclear cells (PBMC) in other diseases have been reported [11, 12]. However, the lncRNA expression profile of PBMCs in COPD patients remains undone.

In this study, we sought to determine if lncRNAs were differentially expressed in PBMC of patients with COPD and if lncRNA expression may be linked to dysregulated mRNA expression relevant to the pathogenesis of COPD. We analyzed lncRNA and mRNA expression profiles in PBMC from COPD patients versus nonsmokers and smokers without airflow limitation. The lncRNA targets in dysregulated mRNAs were predicted.

2. Methods

2.1. Subjects

Peripheral venous blood was taken in heparin-coated tubes from 20 healthy nonsmokers, 17 smokers without airflow limitation, and 14 COPD patients at the Department of Respiration, The First Affiliated Hospital, Xi’an Jiaotong University, Xi’an, China. The COPD patients were eligible for this study if they met the following criteria: age ≥ 50 and ≤76 years; smoking history (≥20 pack years); postbronchodilator FEV1 ≥ 25% of predicted value and postbronchodilator FEV1/forced vital capacity (FVC) ≤ 0.70; and no history of asthma, atopy (as assessed by an allergy skin prick test during screening) or any other active lung disease. Patients on home oxygen or with raised carbon dioxide tension (>44 mmHg), α1-antitrypsin deficiency, recent exacerbation (in the last 4 weeks), an uncontrolled medical condition, or hypersensitivity to inhaled corticosteroids and bronchodilators were not eligible for the study. All nonsmokers and smokers without airflow limitation met the following criteria: age ≥ 42 and ≤75 years, post-BD FEV1% predicted > 80, no diagnostic cancer, diabetes, cardiovascular disease and hypertension, no use of inhaled or oral corticosteroids in the previous 6 months, no atopy, and no respiratory tract infection 1 month prior to the study. Patient characteristics are in Table 1. The experimental procedures were carried out in accordance with the approved guidelines. All experimental protocols were approved by the Research Ethics Boards of The First Affiliated Hospital, Xi’an Jiaotong University (2015-015) and the informed consent was obtained from all subjects.


Healthy nonsmokersSmokers without COPDCOPD

Number201714
Age61 ± 1456 ± 1769 ± 8
Male/female12/817/013/1
Current/ex-smokers14/38/6
Post-BD FEV1% predicted101.0 ± 8.692.9 ± 9.029.0 ± 10.6
Gold stage
 I0
 II6
 III-IV8

Data are presented as mean ± SD. BD: bronchodilator; FEV1, forced expiratory volume in 1 s.
2.2. PBMC Isolation and RNA Extraction

PBMCs were isolated from venous blood by density gradient centrifugation over Ficoll-Paque PLUS reagent (GE Healthcare, Uppsala, Sweden) and suspended in QIAzol Lysis Reagent (Qiagen, Dusseldorf, Germany). Total RNA was extracted using miRNeasy Mini Kit (Qiagen) according to the manufacturer’s procedure. RNA integrity was determined by formaldehyde-denaturing gel electrophoresis. We followed the methods of the previous study [10].

2.3. lncRNA and mRNA Microarray

Equal amount of RNA sample from each nonsmokers (), smokers (), or COPD patients () was pooled, respectively, for lncRNA and mRNA profiling assay using Agilent Human lncRNA + mRNA Array v4.0 system by CapitalBiotech Company, Beijing, China. Each array contained a probe set comprising 40,000 human lncRNA transcripts and 34,000 human mRNAs. These lncRNA and mRNA target sequences were merged from multiple databases including Refseq, UCSC, H-InvDB, Human lincRNA catalog, NRED, lncRNAdb, and RNAdb, and 848 were from the Chen Ruisheng Lab (Institute of Biophysics, Chinese Academy of Science). Briefly, first-strand cDNAs were synthesized from 1 μg of total RNA using random or poly T primers carrying the T7 promoter sequence and the CbcScript II reverse transcriptase. Second strand cDNAs were then synthesized using RNaseH and DNA polymerase. The double-stranded cDNA was column-purified and used as a template to amplify cRNA by in vitro transcription reaction. The amplified cRNAs were purified and reverse transcribed into first-strand cDNA using CbcScript II reverse transcriptase. The second cDNA strand was then synthesized using Klenow enzyme, random primers, regular dNTP, and Cy3- or Cy5-labeled dCTP. The labeled cDNAs were hybridized with an array of analysis. Hybridized slides were then washed and scanned with Agilent Microarray Scanner System (G2565CA).

2.4. Data Analysis

The lncRNA and mRNA array data were analyzed for data summarization, normalization, and quality control using GeneSpring V11.5 software (Agilent). To select differentially expressed genes, we used threshold values of ≥2-fold change, and a Benjamini–Hochberg-corrected value of 0.05 performed on technically duplicated dots for each lncRNA. The data were Log2 transformed and median centered by genes using the Adjust Data function of CLUSTER 3.0 software. Further analysis was performed by hierarchical clustering with average linkages. Finally, we performed tree visualization using Java TreeView (Stanford University School of Medicine, Stanford, CA, USA).

2.5. Quantitative Reverse Transcription PCR Validation

Independent assays were performed using quantitative reverse transcription PCR (qRT-PCR) on all patient samples for individual lncRNA (NR_026891.1, ENST00000502883.1, HIT000648516, XR_429541.1, and ENST00000597550.1) and mRNAs (CXCL16, HMOX1, SLA2, and SIGLEC14) predicted to be regulated by lncRNAs. Total RNA was extracted using miRNeasy Mini Kit (Qiagen) according to the manufacturer’s procedure. Quality control and RNA concentrations were determined by spectrophotometer (IMPLEN, Munich, Germany). The reverse transcription was performed on 500 ng of total RNA by using the iScript™ cDNA Synthesis kit (Bio-Rad). The cDNA was then amplified by using the CFX Connect™ Real-Time PCR Detection System (Bio-Rad) with SYBR Green (Bio-Rad) and the primers listed in Table 2. In addition, the expression of the aforementioned lncRNA was detected on the isolated different cell types including CD4+ T cells, CD8+ T cells, CD14+ monocytes, and CD20+ B cells from PBMCs in some smokers and COPD patients by positive selection (Anti-PE MicroBeads UltraPure, Miltenyi Biotec, Teterow, Germany). Data were presented relative to β-actin for lncRNAs and mRNAs based on calculations of 2−ΔΔCt. Statistical significance was defined as as measured by the t-test using GraphPad Prism 5 software (GraphPad, San Diego, CA, USA).


GeneSequence (5-3)Direction

NR_026891.1GGACCACATCTCCTGTACCAForward
CCTCATGACGTGCACTTTACCReverse
ENST00000502883.1GTGTCCATGTAACTAACTCCTGGTForward
GTCAGCGGGAAGGAAGACAGReverse
HIT000648516ATAGAAGCGATCTACCCTCACAGForward
TGCTGGGCTCGTTCGTReverse
XR_429541.1TGGTCACTTTCCAGTTTCCACAForward
AATTCAGATCCCACATCAGCCTReverse
ENST00000597550.1GCCCTCCACCCCAGATTAACForward
CTCTTTCTTCTCTTCTGGACTTCTReverse
CXCL16CTGAGAGCTTACCATCGGTGTForward
TCAAGACAGCTCATCAATTCCTReverse
HMOX1GGCCAGCAACAAAGTGCAAGForward
ATGGCATAAAGCCCTACAGCAReverse
SLA2GACATCTGCTGCCTACTCAAGGForward
TGTGGCAGCTTCAGAAAACAGGReverse
SIGLEC14CTGCACAGTTGACAGCAACCForward
GTCTGGGAAGGATTGAGGGCReverse
β-ActinTACCTCATGAAGATCCTCACCForward
TTTCGTGGATGCCACAGGACReverse

2.6. lncRNA Target Prediction and Gene Ontology Analysis

The targets of differentially expressed lncRNAs were identified via cis or transregulatory effects. The validated differentially expressed lncRNAs were selected for target prediction. The genes transcribed within a 10 kb upstream or downstream of lncRNAs were considered as cis target genes. lncRNAs and potential cis target genes were paired and visualized using UCSC genome browser. The transtarget genes were chosen by BLAST software according to the impact of lncRNA binding on complete mRNA molecules. The gene ontology (GO) enrichment of dysregulated mRNAs predicted to be regulated by lncRNAs was analyzed by an integrated functional link enrichment of Gene Ontology or gene sets (http://lego.blueowl.cn/).

2.7. Stimulation of Cigarette Smoke Extract

The T lymphocyte cell line 6T-CEM (ATCC, Manassas, VA, USA) was stimulated with 5% cigarette smoke extract (CSE) [13] for 24 h, and the mRNA expression of ENST00000502883.1 and CXCL16 was examined by qRT-PCR.

2.8. lncRNA ENST00000502883.1 Knockdown by Small Interfering RNA (siRNA) Transfection

6T-CEM cells were cultured in 24-well plates with a 5~7105/ml density in RPMI1640 medium containing 10% FBS. The cells were then transfected with siRNA targeting lncRNA ENST00000502883.1 (sequence: 5-CAAACGUUCAUGUGAAAGATT-3; 5-UCUUUCACAUGAACGUUUGTT-3) and synthesized by GenePharma, Shanghai, China, using the siRNA transfection reagent (Roche, Penzberg, Germany) for 48 h. The scrambled siRNA (GenePharma) was used as a negative control. Six hours before the end of the experiment, 80 nM of PMA, 1 μg/mL of ionomycin and10 μg/mL of Brefeldin A was added to the culture. The efficiency of knockdown was examined by qRT-PCR. The expression of target gene CXCL16 was detected by qRT-PCR and Western blot.

2.9. Western Blot

The protein samples of 6T-CEM cells were loaded (5 μg) on a 10% acrylamide SDS-PAGE gel (Bio-Rad, Hercules, CA, USA) for protein separation, followed by transfer to PVDF membranes (Bio-Rad). The blots were then blocked with 1% BSA in 0.1% Tween 20/PBS for 1 h at room temperature and then incubated overnight at 4°C with antibodies specific for CXCL16 (PeproTech, Rocky Hill, CT). After washing with 0.1% Tween 20 in PBS, the membranes were incubated with a 1 : 3000 dilution of goat anti-rabbit IgG HRP (Calbiochem) in 1% solution of powdered milk in PBS/0.1% Tween 20. The membranes were exposed to ECL solution (Bio-Rad) and imaged by chemiluminescence (Clinx Science Instrument, Shanghai, China).

2.10. Chemotaxis Analysis

6T-CEM cells were transfected with lncRNA ENST00000502883.1 siRNA or scrambled siRNA for 48 h, the supernatant was collected for chemotaxis analysis. The PBMCs were isolated and suspended at a concentration of 106 cells/ml in chemotaxis buffer (RPMI 1640 containing 25 mM Hepes and 1% (v/w) endotoxin-free bovine serum albumin). The chemotaxis protocol was performed using a 48-well microchemotaxis Boyden chamber (Neuro Probe, Cabin John, MD) with 5 μm pore polycarbonate filters (Neuro Probe). The inferior wells were loaded with cell culture supernatants pretreated at 37°C for 30 min with neutralizing Ab against CXCL16 (R&D systems, Minneapolis, MN, USA) or goat IgG isotype control (R&D systems), chemotaxis buffer, and CXCL12 (PeproTech) at 10−7 M were used as negative and positive controls, respectively. 20 ng/ml of CXCL16 (PeproTech) pretreated with neutralizing Ab against CXCL16 (R&D systems) or goat IgG isotype control (R&D systems) was loaded into the inferior wells as well. The chemotaxis system was conducted for 2 h 30 min at 37°C in 5% CO2. Each condition was performed in triplicate. Cells having migrated through the filter were counted in the inferior well, and results were expressed as index of chemotaxis compared with chemotaxis buffer.

2.11. Statistical Analysis

Statistical analysis for expression of lncRNA and mRNA by qRT-PCR in PBMCs was performed by Mann-Whitney U test. The paired t-test was performed for cell culture experiments and chemotaxis analysis. Probability values of were considered significant. Data analysis was performed by using the GraphPad prism 5 software (GraphPad, San Diego, CA, USA).

3. Results

3.1. lncRNA Microarray

We firstly examined the lncRNA profiling in PBMCs from 20 nonsmokers, 17 smokers without airflow limitation, and 14 COPD patients. We compared the lncRNA expression between all paired groups (smokers versus nonsmokers, COPD versus nonsmokers, and COPD versus smokers) (Figures 1(a), 1(b), and 1(e)). Compared with nonsmokers, 27 lncRNAs were upregulated and 62 were downregulated in smokers (Table 3), while 165 lncRNAs were upregulated and 81 downregulated in COPD patients (Table 4). When comparing the lncRNA expression between smokers and COPD patients, there were 110 upregulated and 48 downregulated in COPD patients (Table 5).


lncRNAFold change value

Upregulation
ENST00000594469.14.1320.0049
TCONS_000176562.7080.0073
ENST00000452347.12.6100.0076
ENST00000445076.12.5230.0080
ENST00000517658.12.4950.0081
ENST00000379928.42.4880.0082
ENST00000425031.12.4690.0083
ENST00000571404.12.3980.0086
ENST00000417071.12.3920.0087
ENST00000447643.12.3720.0088
Downregulation
ENST00000602863.10.0040.0028
ENST00000445814.10.0180.0029
TCONS_000173430.0470.0031
ENST00000602587.10.0590.0031
ENST00000553269.10.1540.0040
uc.1730.1810.0044
ENST00000420213.10.2070.0047
TCONS_000021060.2200.0049
ENST00000502883.10.2210.0049
ENST00000434051.10.2260.0050


lncRNAFold change value

Upregulation
ENST00000594469.120.2230.0049
ENST00000416105.119.6040.0030
XR_428545.118.0630.0030
NR_103548.18.2890.0034
TCONS_000109847.5620.0035
TCONS_000104037.5220.0035
ENST00000416758.17.1130.0035
HIT0000646975.4780.0038
ENST00000513492.15.2920.0039
ENST00000607854.14.6220.0041
Downregulation
TCONS_000099620.1330.0041
ENST00000456917.10.1340.0041
ENST00000609385.10.1360.0041
ENST00000517983.10.1560.0045
ENST00000420213.10.1990.0053
TCONS_000053140.2090.0055
ENST00000584923.10.2120.0055
TCONS_000083600.2210.0057
ENST00000445814.10.2380.0062
ENST00000602813.10.2390.0062


lncRNAFold change value

Upregulation
ENST00000602863.179.1440.0028
ENST00000446595.119.7150.0030
TCONS_0001634017.9020.0030
XR_428545.114.4260.0031
TCONS_0001734313.4160.0031
ENST00000416105.113.0570.0032
ENST00000445814.112.9710.0032
TCONS_000092347.2670.0036
ENST00000434051.16.9200.0036
ENST00000416758.16.4760.0037
Downregulation
TCONS_000289040.0560.0032
ENST00000517983.10.0920.0035
ENST00000456917.10.1050.0037
ENST00000609385.10.1150.0038
TCONS_000053140.1870.0048
HIT000648524_020.2110.0053
TCONS_000099620.2160.0053
TCONS_000083600.2230.0055
XR_429946.10.2640.0065
ENST00000548760.20.2680.0066

3.2. mRNAs Microarry

We performed the parallel mRNA microarray on pooled RNA samples to compare the mRNA expression profiling in PBMCs from 20 nonsmokers, 17 smokers without airflow limitation, and 14 COPD patients. Eighty-two upregulated and 83 downregulated mRNAs were found in smokers compared with nonsmokers; 190 mRNA were upregulated and 156 downregulated in COPD patients compared with nonsmokers and 135 upregulated and 112 downregulated in COPD patients compared with smokers (Figures 1(c)1(e), part of the results was published in the previous study [10]). Table 6 shows the top 10 upregulated and 10 downregulated genes (part of results was published in the previous study [10]).


Gene symbolGene nameFold change
Smokers versus nonsmokers ( value)
Fold change
COPD versus nonsmokers
( value)
Fold change
COPD versus smokers
( value)
Function

Upregulation
CD177CD177 molecule3.48 (0.0056)78.70 (0.0029)22.59 (0.0030)Leukocyte migration
MUC17Mucin 17,cell surface associated1.08 (0.2270)23.22 (0.0030)21.54 (0.0030)Extracellular matrix constituent
IL1R2Interleukin 1 receptor, type II0.92 (0.7332)10.08 (0.0033)10.96 (0.0033)Decoy receptor, inhibits the activity of IL-1
SARDHSarcosine dehydrogenase0.19 (0.0045)1.83 (0.0100)9.52 (0.0033)Mitochondrial matrix
EGR3Early growth response 31.12 (0.1549)7.57 (0.0035)6.75 (0.0036)Positive regulation of endothelial cell proliferation
AREGAmphiregulin0.65 (0.0319)4.01 (0.0071)6.16 (0.0048)EGF family, promote the growth of normal epithelial cells
SLC6A2Solute carrier family 6 (neurotransmitter transporter, noradrenalin), member 20.39 (0.0093)1.84 (0.0099)4.70 (0.0042)Sodium symporter
TMEM167ATransmembrane protein 167A0.79 (0.8343)3.68 (0.0435)4.55 (0.0077)Golgi apparatus
KCNJ15Potassium inwardly-rectifying channel, subfamily J, member 150.78 (0.5713)3.45 (0.0049)4.42 (0.0043)Potassium channel activity
FCHO1FCH domain only 10.26 (0.0057)1.15 (0.0381)4.39 (0.0043)Clathrin-mediated endocytosis
Downregulation
IL1AInterleukin 1, alpha2.03 (0.0117)0.05 (0.0031)0.02 (0.0029)Immune response
IL6Interleukin 6 (interferon, beta 2)1.59 (0.0111)0.10 (0.0044)0.06 (0.0032)Proinflammatory and anti-inflammatory role
CXCL10Chemokine (C-X-C motif) ligand 101.06 (0.2735)0.07 (0.0033)0.07 (0.0033)Leukocyte chemotaxis
TNFTumor necrosis factor0.43 (0.0112)0.04 (0.0030)0.085 (0.0035)Inflammation, cause apoptosis
CCL20Chemokine (C-C motif) ligand 201.06 (0.2724)0.13 (0.0041)0.13 (0.0039)Lymphocytes chemotaxis
CCL4Chemokine (C-C motif) ligand 40.87 (0.3686)0.14 (0.0045)0.16 (0.0047)Leukocyte chemotaxis
CCL3L3Chemokine (C-C motif) ligand 3-like 30.89 (0.5755)0.17 (0.0066)0.19 (0.0063)Leukocyte chemotaxis
C9orf7Chromosome 9 open reading frame 71.40 (0.0126)0.27 (0.0070)0.19 (0.0049)Calcium channel activity
IL1RNInterleukin 1 receptor antagonist0.90 (0.5833)0.19 (0.0051)0.21 (0.0053)Inhibition of the activities of IL-1
RNF19BRing finger protein 19B2.21 (0.0099)0.53 (0.0362)0.24 (0.0059)Cytotoxic effects of natural killer (NK) cells

Table reproduced from Dang et al. [10], under the Creative Commons Attribution License/public domain.
3.3. Target Prediction

We further predicted the potential cis- and transtarget genes for the dysregulated lncRNAs between COPD patients and smokers within the dysregulated mRNAs. The predicted regulation network was shown in Figure 2. The GO enrichment analysis showed that the biological processes and molecular functions including leukocyte migration, immune response, and apoptosis are the main enriched GOs in the dysregulated mRNAs predicted to be regulated by lncRNAs (Table 7). According to the role of target mRNAs in the pathogenesis of COPD, the 5 lncRNAs NR_026891.1, ENST00000502883.1, HIT000648516, XR_429541.1, and ENST00000597550.1 were selected for the further validation (Tables 8 and 9).


GO valueEnrich scoreGene list

GO:0034612 response to tumor necrosis factor0.00027526.766CCL3, THBS1, CXCL16, and CASP8
GO:0002685 regulation of leukocyte migration0.0024122.613HMOX1, CCL3, and THBS1
GO:1,902,041 regulation of extrinsic apoptotic signaling pathway via death domain receptors0.0025241.244THBS1 and HMOX1
GO:2,001,236 regulation of extrinsic apoptotic signaling pathway0.0049415.411HMOX1, THBS1, and CASP8
GO:0002443 leukocyte-mediated immunity0.0049414.489SLA2, IGHG1–001, and CCL3
GO:0002757 immune response-activating signal transduction0.004949.249SLA2, C5AR1, IGHG1–001, CASP8, and BIRC3
GO:0002687 positive regulation of leukocyte migration0.0049423.266CCL3 and THBS1
GO:2,001,233 regulation of apoptotic signaling pathway0.004948.922HMOX1, SFPQ, THBS1, and CASP8
GO:0071356 cellular response to tumor necrosis factor0.0049420.39CCL3 and THBS1
GO:0010035 response to inorganic substance0.004948.309HMOX1, PPP1R15B, THBS1, and CASP8


lncRNAFold change
Smokers versus nonsmokers
Fold change
COPD versus nonsmokers
Fold change
COPD versus smokers

NR_026891.11.152.762.39
ENST00000502883.10.220.783.51
HIT0006485160.511.142.25
XR_429541.11.132.782.46
ENST00000597550.10.791.612.05


lncRNAGene symbolTarget geneMicroarray analyze fold changeFunction

NR_026891.1FLJ10038CXCL162.29Induced by the inflammatory cytokines IFN-gamma and TNF-alpha; cytokine-cytokine receptor interaction, organism-specific biosystem
ENST00000502883.1RP11-499E18.1CXCL162.29Induced by the inflammatory cytokines IFN-gamma and TNF-alpha; cytokine-cytokine receptor interaction, organism-specific biosystem
HIT000648516HMOX12.71Keap1-Nrf2 pathway, organism-specific biosystem
XR_429541.1SLA20.46Negative regulation of B cell activation; regulation of immune response
TCR signaling in naive CD4+ T cells
ENST00000597550.1CTD-2245F17.3SIGLEC142.57Innate immune system
Cell adhesion

3.4. qPCR Validation of Dysregulated lncRNAs and mRNAs in PBMCs

The upregulated expression of 5 selected lncRNA (NR_026891.1, ENST00000502883.1, HIT000648516, XR_429541.1, and ENST00000597550.1) in COPD patients compared with smokers was further validated by qRT-PCR (Figure 3(a)). In PBMCs of smokers, XR_429541.1 was the most highly expressed lncRNA. In PBMCs of COPD patients, ENST00000502883.1, HIT000648516, and XR_429541.1 had the higher expression (Figure 3(b)). To analyze the relationship between lncRNA expression and lung function of subjects, the correlation analysis between lncRNA expression and FEV1% predicted was performed. The significant negative relevance appeared in the expression of these lncRNAs (Figure 3(c)). However, there is no difference in lncRNA sets between ex and current smokers in the smokers and COPD groups (data not shown).

3.5. qPCR Validation of Dysregulated lncRNAs in Different Cell Types of PBMCs

PBMCs consist mainly of T lymphocytes, B lymphocytes, and monocytes. We therefore analyzed the expression of dysregulated lncRNAs in the isolated different cell types including CD4+ T cells, CD8+ T cells, CD14+ monocytes, and CD20+ B cells from PBMCs of smokers (Figure 4(a)) and COPD patients (Figure 4(b)). NR_026891.1 was highly expressed in monocytes and B cells of smokers, while it was expressed mainly in CD8+ T cells and B cells of COPD patients, suggesting CD8+ T cells contributed to the increased expression in COPD. ENST00000502883.1 was consistently expressed higher in B cells and CD4+ T cells. In COPD patients, the expression of HIT000648516 in monocytes was decreased compared with smokers, suggesting the expression in other cells was increased in COPD patients. XR_429541.1 was mainly expressed in CD4+ and CD8+ T cells in both smokers and COPD patients. CD8+ T cells mainly contributed to the increased expression of ENST00000597550.1 in COPD patients.

3.6. qPCR Validation of Dysregulated Predicted Target mRNAs in PBMCs

We further examined the expression of 4 predicted target genes known to be involved in the pathogenesis of COPD in PBMCs from smokers and COPD patients. Consistent with the results of microarray assay, the increased expression of HMOX1, SIGLEC14, and CXCL16 and decreased expression of SLA2 was observed in COPD patients compared with smokers (Figure 5).

3.7. Validation of ENST00000502883.1-Regulated CXCL16 Expression in T Cell Line

We examined the regulatory role of ENST00000502883.1 on CXCL16 expression. Firstly, 5% of cigarette smoke extract could significantly enhance the expression of ENST00000502883.1 and CXCL16 on T cells (Figure 6(a)). ENST00000502883.1 knockdown with siRNA consequently decreased the expression of CXCL16 at both mRNA (Figure 6(b)) and protein (Figure 6(c)) level, suggesting ENST00000502883.1 could regulate CXCL16 expression. We then analyzed the chemotactic effect of T cell supernatant on PBMCs. The pretreatment with CXCL16 neutralizing antibody significantly decreased the chemotactic effect of NC control culture medium, suggesting the CXCL16 secreted by T cells play a good biological chemotactic function. Meanwhile, the supernatant of T cells with siRNA transfection had the reduced effect compared with NC control, which indicted lncRNA ENST00000502883.1 could regulate CXCL16 expression and furtherl regulate the recruitment of inflammatory cells.

4. Discussion

Increasing evidence has confirmed lncRNAs to be one of the most important factors controlling gene expression [1416]. Therefore, for the first time, we evaluated the lncRNA and mRNA expression profile of PBMCs from healthy nonsmokers, smokers without airflow limitation, and COPD patients. We identified 158 differentially expressed lncRNAs in PBMCs from COPD subjects compared with smokers without airflow limitation. The dysregulated expression of 5 selected lncRNA including NR_026891.1 (FLJ10038), ENST00000502883.1 (RP11-499E18.1), HIT000648516, XR_429541.1, and ENST00000597550.1 (CTD-2245F17.3) were validated by qPCR. The GO enrichment analysis for the dysregulated mRNAs predicted to be regulated by lncRNAs showed that leukocyte migration, immune response, and apoptosis are the main enriched processes that previously reported to be involved in the pathogenesis of COPD. We further validated the dysregulation of target genes CXCL16, HMOX1, SLA2, and SIGLEC14 in PBMCs of COPD patients. It is noted that the predicted target genes discussed in the present study did not overlay the ones that we predicted in the study of dysregulated miRNAs of COPD patients [10]. This may suggest that the regulation profile of lncRNA and miRNA on target genes is different in PBMCs of COPD patients.

For the economic consideration, the equal amount of RNA sample from each smoker and COPD patients was pooled, respectively, for lncRNA profiling assay. And the expression of selected lncRNAs in each individual was further validated by qRT-PCR. The similar approach was previously used on lncRNA microarray assays [1719].

CXCL16 was reported as one of the systemic inflammatory markers for COPD in a large cohort of COPD patients and controls [20]. In this study, increased CXCL16 expression was also found in PBMCs of COPD patients versus smokers. CXCL16 is expressed by dendritic cells, macrophages, T cells, and B cells and can act as a chemoattractant for Th1 cells, which implies its relevance to COPD [21, 22]. CXCL16 was predicted to be regulated by NR_026891.1 (FLJ10038) and ENST00000502883.1 (RP11-499E18.1) through a transregulation manner. Few studies have been reported on the two lncRNAs except for the upregulation of NR_026891.1 (FLJ10038) in an experimental model of Alzheimer’s disease [23]. In this study, we confirmed the regulatory role of ENST00000502883.1 on CXCL16 expression and consequently the effect on PBMC recruitment.

The lncRNA HIT000648516 was predicted to target the HMOX1 gene by a sense manner. HMOX1 gene encodes heme oxygenase (HO)-1, which is induced under physiological conditions such as inflammation and oxidative stress and protects against inflammatory- and oxidant-mediated cellular injury [24]. HO-1 may also play a vital function in maintaining cellular homeostasis [25]. Moreover, the increased expression of HO-1 was found in sputum samples obtained at the onset of a severe COPD exacerbation [26] and in the peripheral blood monocytes during acute inflammatory illnesses of children [27]. Thus, the increased expression of HO-1 in PBMCs may reflect a protective response in COPD patients under the environment of inflammation and oxidative stress.

The expression of SLA2 gene was decreased in PBMCs of COPD patients compared with smokers, and lncRNA XR_429541.1 was predicted to regulate this gene by transregulation. The SLA2 gene encodes Src-like adaptor protein-2 (SLAP-2), which shares 36% sequence similarity with SLAP. SLAP-2 is predominantly expressed in hematopoietic cells and plays an inhibitory role in the activation of T cells [28]. In COPD, once activated, T cells are present in the lung and exert their effector functions by attracting other inflammatory cells like neutrophils and macrophages and enhancing their inflammatory functions [8]. Therefore, the downregulation of SLAP-2 in PBMCs seen in this study may be related with inflammation in COPD patients.

The SIGLEC14 gene expression was predicted to be regulated by ENST00000597550.1 (CTD-2245F17.3). Siglecs are a family of sialic acid-binding lectins expressed mainly on innate immune cells [29]. Siglec-14, a Siglec family member with an activating signaling property, is expressed on granulocytes and monocytes. Siglec-14 serum concentration can serve as a useful marker for COPD exacerbation susceptibility and consequential decline in pulmonary function [30]. Patients with COPD who are homozygous null for this allele had fewer inflammatory exacerbations than patients expressing the wild-type allele, which suggests that Siglec-14 may promote inflammatory sequelae caused by neutrophils [30]. Furthermore, inhaled corticosteroids may exert negative effects on treatment through increased Siglec-14 expression [31].

Overall, through lncRNA and mRNA expression profiling in nonsmokers, smokers, and COPD patients, we identified the dysregulated lncRNAs and mRNAs in PBMCs from COPD patients compared to smokers. We further analyzed the regulation network between lncRNAs and mRNAs, where the genes CXCL16, HMOX1, SLA2, and SIGLEC14 were predicted to be regulated by certain lncRNAs through sense or miRNA regulation. This study may provide clues for further studies targeting lncRNAs to control inflammation in COPD.

Conflicts of Interest

There is no competing financial interest.

Authors’ Contributions

Xiaoyan Qu isolated the PBMCs, performed qRT-PCR and cell culture, and analyzed the data. Xiaomin Dang and Dong Shang collected blood samples. Weijia Wang isolated part of the PBMCs. Ying Li performed part of cell culture. Dan Xu coordinated the collection of blood samples. Ying Chang designed the study and drafted the manuscript. All authors reviewed the manuscript.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (31501044), International Cooperation and Exchange Program of Science and Technology of Shaanxi Province (2017KW-044, 2016KW-026), and the Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry. This work was supported in part by Project 985 of Xi’an Jiaotong University. The authors acknowledge Dr. Yong-Ping Shao for his critical review of the manuscript.

References

  1. M. Decramer, W. Janssens, and M. Miravitlles, “Chronic obstructive pulmonary disease,” Lancet, vol. 379, no. 9823, pp. 1341–1351, 2012. View at: Publisher Site | Google Scholar
  2. B. Yan and Z. Wang, “Long noncoding RNA: its physiological and pathological roles,” DNA and Cell Biology, vol. 31, Supplement 1, pp. S-34–S-41, 2012. View at: Publisher Site | Google Scholar
  3. J. S. Mattick, “The genetic signatures of noncoding RNAs,” PLoS Genetics, vol. 5, no. 4, article e1000459, 2009. View at: Publisher Site | Google Scholar
  4. S. Lalevee and R. Feil, “Long noncoding RNAs in human disease: emerging mechanisms and therapeutic strategies,” Epigenomics, vol. 7, no. 6, pp. 877–879, 2015. View at: Publisher Site | Google Scholar
  5. M. J. Herriges, D. T. Swarr, M. P. Morley et al., “Long noncoding RNAs are spatially correlated with transcription factors and regulate lung development,” Genes & Development, vol. 28, no. 12, pp. 1363–1379, 2014. View at: Publisher Site | Google Scholar
  6. M. M. Perry, E. Tsitsiou, P. J. Austin et al., “Role of non-coding RNAs in maintaining primary airway smooth muscle cells,” Respiratory Research, vol. 15, no. 1, p. 58, 2014. View at: Publisher Site | Google Scholar
  7. H. Bi, J. Zhou, D. Wu et al., “Microarray analysis of long non-coding RNAs in COPD lung tissue,” Inflammation Research, vol. 64, no. 2, pp. 119–126, 2015. View at: Publisher Site | Google Scholar
  8. M. G. Cosio, J. Majo, and M. G. Cosio, “Inflammation of the airways and lung parenchyma in COPD: role of T cells,” Chest, vol. 121, no. 5, pp. 160S–165S, 2002. View at: Publisher Site | Google Scholar
  9. W. Stankiewicz, M. P. Dabrowski, A. Chcialowski, and T. Plusa, “Cellular and cytokine immunoregulation in patients with chronic obstructive pulmonary disease and bronchial asthma,” Mediators of Inflammation, vol. 11, no. 5, pp. 307–312, 2002. View at: Publisher Site | Google Scholar
  10. X. Dang, X. Qu, W. Wang et al., “Bioinformatic analysis of microRNA and mRNA regulation in peripheral blood mononuclear cells of patients with chronic obstructive pulmonary disease,” Respiratory Research, vol. 18, no. 1, p. 4, 2017, Epub 2017/01/07. View at: Publisher Site | Google Scholar
  11. J. Song, D. Kim, J. Han, Y. Kim, M. Lee, and E. J. P. B. M. C. Jin, “Exosome-derived Hotair is a critical regulator and potent marker for rheumatoid arthritis,” Clinical and Experimental Medicine, vol. 15, no. 1, pp. 121–126, 2015. View at: Publisher Site | Google Scholar
  12. W. Sui, Q. Yan, H. Li et al., “Genome-wide analysis of long noncoding RNA expression in peripheral blood mononuclear cells of uremia patients,” Journal of Nephrology, vol. 26, no. 4, pp. 731–738, 2013. View at: Publisher Site | Google Scholar
  13. Y. Chang, L. Al-Alwan, S. Audusseau et al., “Genetic deletion of IL-17A reduces cigarette smoke-induced inflammation and alveolar type II cell apoptosis,” American Journal of Physiology Lung Cellular and Molecular Physiology, vol. 306, no. 2, pp. L132–L143, 2014. View at: Publisher Site | Google Scholar
  14. A. N. Khachane and P. M. Harrison, “Mining mammalian transcript data for functional long non-coding RNAs,” PLoS One, vol. 5, no. 4, article e10316, 2010. View at: Publisher Site | Google Scholar
  15. Y. G. Chen, A. T. Satpathy, and H. Y. Chang, “Gene regulation in the immune system by long noncoding RNAs,” Nature Immunology, vol. 18, no. 9, pp. 962–972, 2017. View at: Publisher Site | Google Scholar
  16. L. Magagula, M. Gagliardi, J. Naidoo, and M. Mhlanga, “Lnc-ing inflammation to disease,” Biochemical Society Transactions, vol. 45, no. 4, pp. 953–962, 2017. View at: Publisher Site | Google Scholar
  17. B. Lee, J. Mazar, M. N. Aftab et al., “Long noncoding RNAs as putative biomarkers for prostate cancer detection,” The Journal of Molecular Diagnostics, vol. 16, no. 6, pp. 615–626, 2014. View at: Publisher Site | Google Scholar
  18. J. Carlevaro-Fita, A. Rahim, R. Guigo, L. A. Vardy, and R. Johnson, “Cytoplasmic long noncoding RNAs are frequently bound to and degraded at ribosomes in human cells,” RNA, vol. 22, no. 6, pp. 867–882, 2016. View at: Publisher Site | Google Scholar
  19. D. Scarano, R. Rao, and G. Corrado, “In silico identification and annotation of non-coding RNAs by RNA-seq and de novo assembly of the transcriptome of tomato fruits,” PLoS One, vol. 12, no. 2, article e0171504, 2017. View at: Publisher Site | Google Scholar
  20. T. M. Eagan, T. Ueland, P. D. Wagner et al., “Systemic inflammatory markers in COPD: results from the Bergen COPD cohort study,” The European Respiratory Journal, vol. 35, no. 3, pp. 540–548, 2010. View at: Publisher Site | Google Scholar
  21. L. E. Donnelly and P. J. Barnes, “Chemokine receptors as therapeutic targets in chronic obstructive pulmonary disease,” Trends in Pharmacological Sciences, vol. 27, no. 10, pp. 546–553, 2006. View at: Publisher Site | Google Scholar
  22. P. Shashkin, D. Simpson, V. Mishin, B. Chesnutt, and K. Ley, “Expression of CXCL16 in human T cells,” Arteriosclerosis, Thrombosis, and Vascular Biology, vol. 23, no. 1, pp. 148-149, 2003. View at: Publisher Site | Google Scholar
  23. F. J. Hoerndli, S. Pelech, A. Papassotiropoulos, and J. Gotz, “Abeta treatment and P301L tau expression in an Alzheimer’s disease tissue culture model act synergistically to promote aberrant cell cycle re-entry,” The European Journal of Neuroscience, vol. 26, no. 1, pp. 60–72, 2007. View at: Publisher Site | Google Scholar
  24. D. Morse and A. M. Choi, “Heme oxygenase-1: the “emerging molecule” has arrived,” American Journal of Respiratory Cell and Molecular Biology, vol. 27, no. 1, pp. 8–16, 2002. View at: Publisher Site | Google Scholar
  25. A. M. Choi and J. Alam, “Heme oxygenase-1: function, regulation, and implication of a novel stress-inducible protein in oxidant-induced lung injury,” American Journal of Respiratory Cell and Molecular Biology, vol. 15, no. 1, pp. 9–19, 1996. View at: Publisher Site | Google Scholar
  26. M. Tsoumakidou, N. Tzanakis, G. Chrysofakis, and N. M. Siafakas, “Nitrosative stress, heme oxygenase-1 expression and airway inflammation during severe exacerbations of COPD,” Chest, vol. 127, no. 6, pp. 1911–1918, 2005. View at: Publisher Site | Google Scholar
  27. A. Yachie, T. Toma, K. Mizuno et al., “Heme oxygenase-1 production by peripheral blood monocytes during acute inflammatory illnesses of children,” Experimental Biology and Medicine, vol. 228, no. 5, pp. 550–556, 2003. View at: Publisher Site | Google Scholar
  28. S. J. Holland, X. C. Liao, M. K. Mendenhall et al., “Functional cloning of Src-like adapter protein-2 (SLAP-2), a novel inhibitor of antigen receptor signaling,” The Journal of Experimental Medicine, vol. 194, no. 9, pp. 1263–1276, 2001. View at: Publisher Site | Google Scholar
  29. A. Varki and T. Angata, “Siglecs—the major subfamily of I-type lectins,” Glycobiology, vol. 16, no. 1, pp. 1R–27R, 2006. View at: Publisher Site | Google Scholar
  30. T. Angata, T. Ishii, T. Motegi et al., “Loss of Siglec-14 reduces the risk of chronic obstructive pulmonary disease exacerbation,” Cellular and Molecular Life Sciences, vol. 70, no. 17, pp. 3199–3210, 2013. View at: Publisher Site | Google Scholar
  31. P. Wielgat, R. M. Mroz, A. Stasiak-Barmuta et al., “Inhaled corticosteroids increase siglec-5/14 expression in sputum cells of COPD patients,” Advances in Experimental Medicine and Biology, vol. 839, pp. 1–5, 2015. View at: Publisher Site | Google Scholar

Copyright © 2018 Xiaoyan Qu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

2704 Views | 397 Downloads | 4 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.