integrated analysis of lncrna and mrna transcriptomes
TRANSCRIPT
Research ArticleIntegrated Analysis of lncRNA and mRNA TranscriptomesReveals New Regulators of Ubiquitination and the ImmuneResponse in Silica-Induced Pulmonary Fibrosis
Yao Zhou,1 Li He,1 Xiao-Dan Liu,2 Hua Guan,2 Ying Li,3
Rui-Xue Huang ,1 and Ping-Kun Zhou 2,4
1Department of Occupational and Environmental Health, Xiangya School of Public Heath, Central South University, Changsha,Hunan Province 410078, China
2Department of Radiation Biology, Beijing Key Laboratory for Radiobiology, Beijing Institute of Radiation Medicine, AMMS,Beijing 100850, China
3Hunan Prevention and Treatment Center for Occupational Diseases, Changsha, China4Institute for Chemical Carcinogenesis, State Key Laboratory of Respiratory, School of Public Health, Guangzhou Medical University,Guangzhou 511436, China
Correspondence should be addressed to Rui-Xue Huang; [email protected] and Ping-Kun Zhou; [email protected]
Yao Zhou and Li He contributed equally to this work.
Received 30 October 2018; Revised 19 December 2018; Accepted 26 December 2018; Published 13 January 2019
Academic Editor: Jinsong Ren
Copyright © 2019 Yao Zhou et al.This is an open access article distributed under theCreative Commons Attribution License, whichpermits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Objectives. As an epigenetic player, long noncoding RNAs (LncRNAs) have been reported to participate in multiple biologicalprocesses; however, their biological functions in silica-induced pulmonary fibrosis (SIPF) occurrence and development remainincompletely understood.Methods. Five case/control pairs were used to perform integrated transcriptomes analysis of lncRNA andmRNA. Prediction of lncRNA and mRNA functions was aided by the Gene Ontology (GO) and Kyoto Encyclopedia of Genes andGenomes (KEGG) databases. Additionally, we constructed a coexpression network of lncRNAs and mRNAs to identify targets ofregulation. Results. In total, 1069 differentially expressed mRNAs and 366 lncRNAs were identified with the changes more than 2times (p<0.05), of which 351 downregulatedmRNA and 31 downregulated lncRNAwere <0.5 (p<0.05) and those of 718 upregulatedmRNAs and 335 upregulated lncRNA were >2 (p<0.05). The levels of 10 lncRNAs were measured via qRT-PCR; the results wereconsistent with the microarray data. Four genes named of FEM1B, TRIM39, TRIM32, and KLHL15 were enriched significantly withubiquitination and immune response. Cytokine-cytokine receptor interaction was the most significantly enriched KEGG pathwayin both mRNAs and lncRNAs. The coexpression network revealed that a single lncRNA can interact with multiple mRNAs, andvice versa. Conclusions. lncRNA and mRNA expression were aberrant in patients with SIPF compared to controls, indicating thatdifferentially expressed lncRNAs and mRNAs may play critical roles in SIPF development. Our study affords new insights into themolecular mechanisms of SIPF and identifies potential biomarkers and targets for SIPF diagnosis and treatment.
1. Introduction
Evidence has been presented that silica dust acute andchronic exposures will trigger an inflammatory cascade,such as macrophages and alveolar epithelial cells proliferatedsuccessively, followed by inflammatory cascade; the fibroticdevelopment will be ignited, which eventually will proceedto determine the silicosis [1, 2]. Silicosis is an irreversible
and incurable pulmonary disorder; even when patients areno longer exposed to silica, the fibrosis remains progressive[3, 4]. Silica dust is widely distributed in workplaces inwhich drilling, grinding, and hammering activities occur[5, 6]. China has the highest silicosis burden worldwide, withmore than 600,000 cases recorded over the past 30 years.Furthermore, the number of new cases is increasing, withover 24,000 deaths annually [7]. In South African gold mines,
HindawiBioMed Research InternationalVolume 2019, Article ID 6305065, 11 pageshttps://doi.org/10.1155/2019/6305065
2 BioMed Research International
the incidence of silicosis is 13–25% in long-term miners[8]. Silicosis is also prevalent in developing countries. TheCDC states that approximately 2 million US workers arecurrently exposed to silica; from 2001 to 2010, 1,437 deathswere attributable, inwhole or in part, to silicosis; the youngestpatient was aged only 19 years [9]. In the UK, over 3.2 millionpersons are occupationally exposed to silica [10]. Many pre-ventative efforts have been made, including dust control andthe development of personal protective equipment; in somedeveloped countries, the incidence has thus steadily declinedover the past few decades, but new cases or outbreaks aresporadically reported [11]. In developing countries, silica-induced silicosis associated with progressive lung fibrosisstill remains a major concern [6, 12]. No effective therapyis available, and the molecular details of fibrotic progressionremain unclear [13]. Early biomarkers are urgently required;these would aid in implementation of preventative measuresand allow for early diagnosis, intervention, and treatment.
Some biomarkers aiding in early diagnosis have beenidentified and can be divided into exposure, effect, andsusceptibility biomarkers. Serum KL-6, MMP-2, and SP-D are potential biomarkers of silicosis, with levels beingsignificantly higher in cases than controls (p<0.05) [14].HO-1, the enzyme catabolizing heme to bilirubin, is a rec-ognized biomarker [15]. Silica increases the ceruloplasminlevel, which thus serves as a diagnostic biomarker [16]. Theserum CC16 level indicates the risk of silica exposure toxicity[17]. However, current biomarkers studies regarding lungfibrosis induced by silica mostly focus on the genes level,seldom associated with transcriptome changes, which is anemerging concern which have been reported to be potentialbiomarkers in the development of lung fibrosis. Recently,epigenetic research has shown that the changes inmicroRNA(miRNA)/long noncoding RNA (lncRNA) expression mayaffect the stability and translation of genes involved in silica-induced lung fibrosis. MiRNAs are short noncoding RNAsthat bind to targeting mRNAs in a sequence-specific manner[18]. For example, miR-489 was increased significantly in ananimal model of lung fibrosis, suppressing TGF-𝛽1 synthesis(and where TGF-𝛽1 is a transcription factor precipitatingthe inflammation associated with lung fibrosis [19]). MiR-449a inhibits fibrosis by targeting Bcl2, in turn regulatingautophagy [20]. miR-19a is downregulated in the early stagesof fibrosis and was proposed as a biomarker for earlydiagnosis [21]. The miR-146a level in bronchoalveolar lavagefluid aids in early diagnosis [22]. lncRNA research onlycommenced in the last decade; lncRNAs are long noncodingRNAs [23] modulating a series of important biological pro-cesses, e.g., metabolism and apoptosis, and serve also as use-ful biomarkers of various diseases, including silica-inducedlung fibrosis. Sun et al. showed that the lncRNAs suc.77 and2700086A05Rik regulated the progression of lung fibrosis[24]. lncRNA-ATB, first identified in 2014, was activated byTGF-𝛽 and induced lung fibrosis [25]. Thus, the potentialutility of lncRNAs as fibrosis biomarkers requires attention,especially the investigation on human cases of silica asso-ciated fibrosis. However, such work in the context of silica-induced lung fibrosis is limited. Here, we collected peripheralblood samples from patients with silica-induced lung fibrosis
followed by transcriptome analysis using a microarray; wethen employed bioinformatics tools to identify changes inlncRNA levels in patients with progressive lung fibrosis; weidentified the top 10 most-affected lncRNAs and furtheranalyzed the interaction networks of lncRNAs and mRNAs.Our study may improve our mechanistic understanding ofthe silicosis and identifies novel diagnostic and therapeuticbiomarkers.
2. Methods and Materials
2.1. Participants. We enrolled 10 subjects: 5 with phase Ilung fibrosis (cases) and 5 healthy subjects who worked inthe same industry but were not exposed to silica (controls).The inclusion criteria were (1) an occupational history ofsilica exposure and a fibrosis diagnosis based on the Chi-nese Pneumoconiosis Diagnosis Standards (GBZ70-2013),as agreed by at least three occupational physicians; (2)lack of complications such as tuberculosis, lung cancer, oran infection; and (3) voluntary agreement to participate.The mean age of the cases was 55±6.2 years and that ofcontrolswas 51±7.9 years; all participantsweremale.The silicadust concentration in the workplace was 6.43±0.77 mg/m3,exceeding the exposure limit of 0.2–1 mg/m3 (GBZ2-2002).At least 5 mL of peripheral blood was collected from eachsubject into a PAX gene Blood RNA tube (BD Biosciences,San Jose, CA, USA), held at room temperature (22–25∘C)for at least 2 h to completely lyse all cells, placed at –70∘Cand immediately delivered to the Beijing Kangpusen BiotechCompany for microarray analysis. The study was performedin accordance with all relevant tenets of the Declarationof Helsinki and was approved by the Ethics Committee ofXiangya School of Public Health (Approval no. XYGW-2018-11). All participants provided written informed consent.
2.2. Total RNA Isolation. Total RNA was isolated usingthe mirVanaRNA Isolation Kit following the manufacturer’sinstructions. Briefly, a 1.25 volume of absolute ethanol wasadded to each blood sample followed by mixing; the mixturewas passed through a filter and the filter was washed threetimes with 700 𝜇L of miRNA wash solution and 500 𝜇Lgeneral wash solution; then, RNA was eluted into 100 𝜇L ofelution solution at 95∘C and stored at –70∘C until analyzed.RNA purity and concentration were assessed using a QIA-GEN RNeasy Mini Kit and by deriving a spectrophotometricratio with the aid of a NanoDrop-2000 instrument, respec-tively (Thermo-Scientific, Waltham, MA, USA).
2.3. Microarray Analysis. In brief, microarray analysis wasperformed as follows: RNA/-poly-A-RNA-control mixtureswere prepared by mixing test RNA samples (50–500 ng), adiluted poly-A-RNA control, and nuclease-free water; 4 𝜇Lof first-strand buffer mix and 1 𝜇L of first-strand enzymesolution were then mixed and added, and first-strand wascDNA synthesized; this was followed by second-strand syn-thesis. Purified cDNA and second-cycle single-strand cDNAwere subjected to WT cartridge hybridization, washed, andscanned using a DNA Microarray Scanner (Agilent, SantaClara, CA, USA). Expression Console software (ver. 1.4.1)
BioMed Research International 3
was used to adjust the raw data background and standardizemicroarray data.
2.4. Differential mRNA and lncRNA Expression and Clus-ter Analysis. Transcriptome Analysis Console (version 3.1,Affymetrix, Santa Clara, CA, USA) software was employed toexplore the differential expression of mRNAs and lncRNAsbetween cases and controls; we selected RNAs exhibiting thegreatest fold differenceswith reference to the p values. Clusteranalysis [26] was used to identify genes exhibiting similarbiological functions.
2.5. Gene Ontology and Pathway Enrichment Analysis. GeneOntology (GO) and pathway enrichment analysis were per-formed to predict biological processes, cellular components,and molecular functions affected in disease. We used acontrolled gene vocabulary (http://www.geneontology.org).The Kyoto Encyclopedia of Genes and Genomes (KEGG)database was employed for pathway enrichment analysis.
2.6. lncRNA-mRNA Coexpression Network. To explore inter-actions among differentially expressed lncRNAs and mRNAsin cases and controls, we prepared coexpression networksusing only differentially expressed mRNAs with Pearsoncorrelation coefficients ≥0.95 (p<0.01) and generated visualdata. We included five essential lncRNAs and all essen-tial mRNAs in a network assessment. Only differentiallyexpressed mRNAs with Pearson correlation coefficients ≥0.90 (p<0.01) were used to construct the lncRNA-mRNAnetwork and generate visual representations.
2.7. Quantitative Real-Time PCR (RT-PCR) Validation. TotalRNA was isolated and reverse-transcribed into cDNA usinga PEXBIO RNA kit according to the manufacturer’s protocol.qPCR was performed with the aid of Super Real PreMix Plus(SYBRGreen, Tiangen, Beijing, China) on a CFX96TMReal-Time System (Bio-Rad, Hercules, CA, USA). GAPDH servedas the internal control and fold change calculations weremade using the 2–Δ Δ CT method. The primer sequences areshown in Table S1. Differences between mRNA and lncRNAexpression levels were compared using Student’s two-tailedt-test. A p value <0.01 was considered to reflect significanceafter false discovery rate correction for multiple testing.
2.8. Statistical Analysis. R software (R Development CoreTeam, Vienna, Austria) was used to compare the expressionlevels of mRNAs and lncRNAs, for cluster analysis, and todraw volcano maps. One-way ANOVA was employed toexplore the significance of differences in mRNA and lncRNAlevels after log
2transformation. The fold changes between
cases and controls were calculated; changes >2 or <0.5 wereconsidered to indicate differential expression; p<0.05 wastaken to indicate significance.
3. Results
3.1. Comprehensive Transcriptome Analysis Reveals Differen-tially Expressed mRNAs and lncRNAs. Through microarray
analysis, 1069 differentially expressed mRNAswere identifiedwith the changes more than 2 times (p<0.05), of which one-way ANOVA showed that the fold changes of 351 downregu-lated mRNAwere< 0.5 (p<0.05) and those of 718 upregulatedmRNAs were > 2 (p<0.05). The most prominently upreg-ulated mRNA was TC0700007850.hg.1 (fold change, 8.00;p<0.05).The top 25 up- and downregulated mRNAs are listedin Table 1. Meanwhile, 366 differentially expressed lncRNAswere identified with the changes more than 2 times (p<0.05),of which one-way ANOVA showed that the fold changes of31 downregulated lncRNAs were <0.5 (p<0.05) and those of335 upregulated lncRNAs were >2 (p<0.05).Themost promi-nently downregulated lncRNA was TC1300009322.hg.1, witha fold change of 0.395 (p<0.05). The most prominentlyupregulated lncRNA was TC1700010761.hg.1, with a foldchange of 6.70 (p<0.05). The top 25 up- and downregulatedlncRNAs are listed in Table 2.
Unsupervised cluster analysis of the mRNAs of the 10samples, in terms of differentially expressed genes, yieldedthe results shown in Figure 1(a). A volcano plot of thesemRNAs based on p values and fold changes is shown inFigure 1(b). Figure 1(b) shows that more mRNAs were up-than downregulated. Unsupervised cluster analysis of thelncRNAs of the 10 samples is shown in Figure 1(c). Highlyexpressed lncRNAs are shown in red and those expressed atlow levels are in green. A volcano plot based on p values andfold changes is shown in Figure 1(d). Analysis of the normal-ized mRNA expression levels revealed significantly differentmRNAs patterns in case group and control group (718 up-and 351 downregulated) along with more mRNAs were up-than downregulated (Figures 1(a) and 1(b)). The upregulatedmRNAs with the greatest fold change over 5, which are indescending rank, were LOYBOY, DEEZOY, SNEYSPORBY,NYSPSWBY, ZORSKABU, TOYPYBU, and WUSLAWBUTO.Meanwhile, the downregulated mRNAswith the greatest foldchange <0.35, which are in descending rank, were CPA3,FCER1A, and CCR5. Analysis of the normalized lncRNAexpression levels revealed significantly different lncRNAspatterns in case group and control group (335 up- and 31downregulated) along with more lncRNAs were up- thandownregulated (Figures 1(c) and 1(d)). The upregulatedlncRNAs with the greatest fold change over 4, which arein descending rank, were RP11-138I1.4, RP4-620F22.2, andRPL41P3. Meanwhile, the downregulated lncRNAs with thegreatest fold change <0.45, which are in descending rank,were RP11-13A1.3, LOC286437, RP4-539M6.21, RP11-333J10.2,and FTX. As expected, the mRNAs and lncRNAs expressionlevels in Figures 1(a)–1(d) were observed to be consistent withcase group and control group.
3.2. Comprehensive Transcriptome Functional Analysis ofDifferentially Expressed mRNAs and lncRNAs Reveals NewRegulators of Ubiquitination and the Immune Response. Tohighlight the biological functions and signaling pathwayschange in differentially expressed mRNAs and lncRNAs,we performed GO analysis and KEGG pathway assessmenton mRNAs and lncRNAs, respectively (Figure 2). FormRNAs, in the main GO biological process affectedwas sensory perception of smell and interleukin-12
4 BioMed Research International
Table 1: The top 25 upregulated and top 25 downregulated mRNAs from case group compared with control group.
Upregulated mRNAs Downregulated mRNAsAccession Symbol Fold p value Accession Symbol Fold p valueTC0X00007799.hg.1 loyboy 5.78 0.012 TC0300009122.hg.1 CPA3 0.30 0.025TC1300008271.hg.1 deezoy 5.54 0.025 TC0100010252.hg.1 FCER1A 0.31 0.022TC0300010812.hg.1 sneysporby 5.44 0.038 TC0300007257.hg.1 CCR5 0.33 0.002TC0300006637.hg.1 nyspawby 5.31 0.016 TC1700011444.hg.1 Y RNA 0.35 0.029TC0600009711.hg.1 zorskabu 5.31 0.019 TC0500009637.hg.1 morflaw 0.35 0.003TC0100009141.hg.1 toypybu 5.10 0.016 TC0300012816.hg.1 P2RY12 0.36 0.0005TC0700012732.hg.1 wuslawbu 5.01 0.017 TC1700010782.hg.1 titari 0.36 0.007TC0100008816.hg.1 IFI44 4.87 0.024 TC1200009872.hg.1 KLRB1 0.36 0.006TC1700009890.hg.1 RP11-138I1.4 4.60 0.005 TC1500009459.hg.1 gozo 0.37 0.002TC0100017874.hg.1 slokubo 4.59 0.034 TC0400008033.hg.1 sagee 0.38 0.0009TC1800008840.hg.1 glorlo 4.46 0.016 TC0500012614.hg.1 slorskarby 0.38 0.0108TC0100011221.hg.1 snawkobu 4.46 0.013 TC0600012379.hg.1 ELOVL4 0.38 0.009TC0500012458.hg.1 doysheeby 4.42 0.012 TC0200012982.hg.1 veyklar 0.39 0.0107TC0100008972.hg.1 RP4-620F22.2 4.37 0.021 TC0300007255.hg.1 CCR3 0.39 0.007TC1400006879.hg.1 snovaw 4.35 0.048 TC0100016431.hg.1 TNFSF4 0.39 0.008TC0500010789.hg.1 seystyby 4.32 0.037 TC0600013284.hg.1 nimime 0.39 0.009TC2200007220.hg.1 RPL41P3 4.29 0.013 TC1200011466.hg.1 RP11-13A1.3 0.40 0.009TC0900011875.hg.1 RNU6ATAC 4.26 0.012 TC0400010745.hg.1 CLOCK 0.40 0.0001TC1100011736.hg.1 Gorbla 4.22 0.021 TC0100008797.hg.1 AK5 0.40 0.013TC2000008618.hg.1 LOC100287 4.20 0.048 TC0X00006723.hg.1 CDKL5 0.40 0.003TC1600008087.hg.1 RP11-2K6.1 4.14 0.037 TC0300010302.hg.1 kleyjubo 0.40 0.0004TC1900011146.hg.1 poyjoy 4.08 0.002 TC0800008361.hg.1 vokey 0.40 0.0006TC0600011130.hg.1 HIST1H4C 4.07 0.013 TC0X00010443.hg.1 LOC286437 0.40 0.0107TC0X00007364.hg.1 SNORA11 3.99 0.021 TC1800006915.hg.1 HRH4 0.40 0.007TC0500011210.hg.1 deychoby 3.92 0.013 TC0700010552.hg.1 snarcheeby 0.40 0.008“---” presents that no symbols are available in the microarray databases, but they were actually presenting dysregulation in the case group.
secretion, while in the affected GO cellular componentswere principally the immune response including secretorydimeric IgA immunoglobulin complex and monomericIgA immunoglobulin complex, while the GO molecularfunctions affected were olfactory receptor activity, cytokinebinding, and, in particular, thiol-dependent ubiquitinylhydrolase activity, ubiquitinyl hydrolase activity, andubiquitin-like protein-specific protease activity (Figure 2(a)).Four genes named FEM1B, TRIM39, TRIM32, and KLHL15were enriched significantly with ubiquitination and immuneresponse. KEGG pathway analysis emphasized olfactorytransduction, chemokine signaling, cytokine-cytokinereceptor interaction, maturity-onset diabetes of the young,and bile secretion (Figure 2(b)). GO and KEGG analysesof differentially expressed lncRNAs are shown in Figure 3.Differentially expressed lncRNAs were predicted to havethe following main functions: keratinization, keratinocytedifferentiation, and epidermal cell differentiation in thebiological process category; intermediate filament andthe GO molecular components in the cellular componentcategory; and G-protein-coupled amine peptide receptoractivity, hormonal activity, transcription factor activity, andRNA polymerase/distal enhancer sequence-specific bindingin the molecular function category (Figure 3(a)). KEGGpathway analysis of all differentially expressed lncRNAs
significant at p<0.05 revealed that the affected pathwaysincluded olfactory transduction, neuroactive ligand-receptor interaction, estrogen signaling, the cholinergicsynapse, renin secretion, and cytokine-cytokine receptorinteraction (Figure 3(b)). In order to explore individuallncRNA function in the development of lung fibrosis, wesubjected the following 10 lncRNAs to individual GO andKEGG analysis: 18S F, RPL41P3, RN7SL541P, RP11.537H15.3,RN7SL293P, RP4-620F22.2, Metazoa-SRP, RP11-609D21.3,RP11-138I1.4, and RN7SL783P. GO analysis revealed mostof these selected individual lncRNAs of involvement inchemokine binding, CoA-transferase activity, and G-protein-coupled serotonin receptor activity, all of whichare associated with SIPF development (Figure S1-S5);intriguingly, some lncRNAs were predicted to have immuneresponse such as RP11.537H15.3 affected immunoglobulinproduction and adaptive immune response. The KEGG datashowed that these selected individual lncRNAs participatedin multiple pathways, including the herpes simplex infectionpathway, a hepatitis C pathway involving inflammation, andthe PPAR signaling pathway.
3.3. Verification of Expression Changes in lncRNAs. To verifythe differential expression of lncRNAs, we quantified theexpression levels of 10 lncRNAs in 20 cases and 20 controls,
BioMed Research International 5
Table 2: The top 25 upregulated and top 25 downregulated lncRNAs from case group compared with control group.
Upregulated lncRNAs Downregulated lncRNAsAccession Symbol Fold p value Accession Symbol Fold p valueTC1700009890.hg.1 RP11-138I1.4 4.603 0.013 TC1200011466.hg.1 RP11-13A1.3 0.400 0.0094TC0100008972.hg.1 RP4-620F22.2 4.37 0.0024 TC0X00010443.hg.1 LOC286437 0.404 0.0012TC2200007220.hg.1 RPL41P3 4.298 0.013 TC2200008440.hg.1 RP4-539M6.21 0.418 0.007TC1700009599.hg.1 RP11-609D21.3 3.66 0.012 TC1700007659.hg.1 RP11-333J10.2 0.434 0.00009TC1500009338.hg.1 RP11-605F22.1 3.26 0.000 TC0X00010064.hg.1 FTX 0.445 0.0026TC1900006725.hg.1 CTC-518P12.6 3.18 0.04 TC0800011834.hg.1 ASAP1-IT1 0.458 0.0095TC1000007371.hg.1 AP001610.5 3.12 0.02 TC1100006449.hg.1 RP11-326C3.10 0.471 0.041TC1600011549.hg.1 skeyglerbu 3.07 0.02 TC1700011172.hg.1 RP11-257O5.2 0.471 0.0015TC1200011218.hg.1 RP11-368L12.1 3.05 0.00 TC0800010416.hg.1 RP11-11C20.3 0.479 0.0020TC0500006875.hg.1 RP11-1143G9.4 2.94 0.00 TC0200015562.hg.1 AC083900.1 0.483 0.0005TC1700011014.hg.1 CTD-2139B152 2.76 0.01 TC0500013315.hg.1 CTD-2306M10.1 0.486 0.013TC0700011233.hg.1 rerlee 2.74 0.01 TC0300008280.hg.1 LINC01215 0.487 0.003TC1700011336.hg.1 rasleebu 2.70 0.01 TC1500009462.hg.1 RP11-519C12.1 0.487 0.0004TC1600011191.hg.1 RP11-264B14.2 2.65 0.03 TC1200011369.hg.1 RP11-530C5.4 0.494 0.0006TC0100009018.hg.1 GS1-21A4.1 2.58 0.02 TC0100007717.hg.1 RP11-415J8.7 0.49 0.0006TC0600012967.hg.1 fawnoybu 2.57 0.00 TC1200006889.hg.1 RP11-180M15.7 0.49 0.005TC0900011401.hg.1 RP1-93H18.6 2.55 0.00 TC0100018423.hg.1 LINC01355 0.50 0.002TC0100017733.hg.1 werteebo 2.50 0.00 TC1200010249.hg.1 RP11-996F15.6 0.51 0.005TC1400008210.hg.1 TSNAX 2.49 0.01 TC1600008532.hg.1 RP11-319G9.3 0.51 0.001TC1000010277.hg.1 RP11-638I2.8 2.49 0.007 TC1100011979.hg.1 RP11-178H8.7 0.51 0.001TC0300009362.hg.1 RN7SL398P 2.46 0.000 TC0300013876.hg.1 EPHB1 0.51 0.003TC0X00009457.hg.1 RP11-3P17.5 2.45 0.04 TC1300007758.hg.1 RP11-10E18.7 0.51 0.0001TC14 GL000009v2 random00006439.hg.1 RN7SL15P 2.43 0.000 TC0100017904.hg.1 CHRM3-AS2 0.519 0.005
TC1100012004.hg.1 RP11-685N10.1TC1100009942.hg.1 RP11-435B5.6 2.43 0.02 TC1400008411.hg.1 RP11-73M18.6 0.21 0.005TC1700008351.hg.1 TRIM22 2.42 0.02“---“ presents that no symbols are available in the microarray databases, but they were actually presenting dysregulation in the case group.
18S F, RPL41P3, RN7SL541P, RP11.537H15.3, RN7SL293P,RP4-620F22.2, Metazoa-SRP, RP11-609D21.3, RP11-138I1.4,and RN7SL783P, via quantitative RT-PCR. The expressionlevels of cases and controls differed significantly (Figure 1(e)),consistent with the microarray data.
3.4. LncRNA-mRNA Network Analysis. To analyze thelncRNA-mRNA coexpression interaction network, wesought genes that were significantly enriched in both casesand controls. A total of 2,053 such genes were identified, andlncRNA-mRNA coexpression networks were constructed;this identified five lncRNAs (RPL41P3, RP11-609D21.3, RP4-620F22.2, AP001610.5, and RN7SL783P) that determinedthe essential differences between cases and controls (FiguresS6). These lncRNAs and their coexpressed mRNAs, mayindicate how silica exerts biological effects and may alsoserve as useful therapeutic targets. For example, AP001610.5was coexpressed with the following mRNAs/genes: nuraro,IFIT3, TAGLN3, bawma, rasuru, IFIT1, CT45A9, nymawbo,and CT45A8. RN7SL783P was coexpressed with many moremRNAs/genes, including pulmonary fibrosis-associatedgenes such as POYNEY and PODUBU. RP4-620F22.2,
RPL41P3, and RP11-605F22.1 were also associated with manymRNAs (Table S2).
4. Discussion
Microarray and bioinformatics techniques now allow us todetect changes in lncRNA expression; most work to date hasfocused on miRNAs [27–33]. lncRNAs are longer sequenceand have more complex functions than miRNAs [34]. Forthe first time, we herein identified mRNAs and lncRNAs thatwere abnormally expressed in human lung fibrosis induced bysilica exposure, which exhibits complex pathological changes;no effective therapy is available. Biomarkers discovered viamRNAs and lncRNAs analysis might facilitate diagnosis andtherapy [35, 36].
Differentially expressed mRNAs and lncRNAs have beenreported to be associated with development of lung fibrosis.Laurent et al. explored the cultured lung fibroblast transcrip-tome of idiopathic pulmonary fibrosis; mRNAs encodingLIMS2, ASB1, andHHAT were upregulated and those encod-ing TRANK1, IFIT1, and SLC15A3 were downregulated [37].Biswas et al. found that, in cystic fibrosis, CTFR expression
6 BioMed Research International
silic
a-1
norm
al-1
silic
a-2
norm
al-2
silic
a-3
norm
al-3
silic
a-4
norm
al-4
silic
a-5
norm
al-5
−2 2−1 0
low C-1
color cay
(a)
log2 (Fold Change)−3 −2 −1 0
0
2
4
6
8
1 2 3
−log10
(P_v
alue
)
(b)
silic
a-1
norm
al-1
silic
a-2
norm
al-2
silic
a-3
norm
al-3
silic
a-4
norm
al-4
silic
a-5
norm
al-5
−2 2−1 0
low C-1
color cay
(c)
log2 (Fold Change)−4 −3 −2 −1 0 1 2 43
0
2
4
6
8
−log10
(P_v
alue
)
(d)
Fold
expr
essio
n ch
ange
s
7
6
5
4
3
2
1
0
18S
RN7SL783P
RN7SL541P
RN7SL293P
RPL41P3
Metazo
a-SRP
RP11-537H15.3
RP11-609D21.3
RP11-138I1.4
RP4-620F22.2
(e)
Figure 1: Differentially expressed mRNAs and lncRNAs among 10 samples. (a) Unsupervised cluster analysis of included differentiallyexpressed mRNAs. Red rectangles mean the clustered upregulated mRNAs, while green rectangles mean clustered downregulated mRNAs.(b) A volcano plot of these mRNAs based on p values and fold changes. Red dots mean the upregulated mRNAs, while green dots meandownregulated mRNAs. (c) Unsupervised cluster analysis of included differentially expressed lncRNAs. Red rectangles mean the clusteredupregulated lncRNAs, while green rectangles mean clustered downregulated lncRNAs. (d) A volcano plot of these mRNAs based on p valuesand fold changes. Red dots mean the upregulated lncRNAs, while green dots mean downregulated lncRNAs. (e) Differential expression oflncRNAs was quantified using quantitative RT-PCR.
BioMed Research International 7
Num
bers
of g
enes
20
15
10
5
0
Biological ProcessCellular ComponentMolecular Function
Gene Function Classification (GO)dete
ction
of ch
emica
l stim
ulus invo
lved in
senso
ry perc
eptio
n
detecti
on of
chem
ical st
imulus i
nvolve
d in se
nsory
percep
tion of
smell
senso
ry perc
eptio
n of sm
ell
interleu
kin-12 se
cretio
n
regulat
ion of intra
cellular
mRNA lo
caliza
tion
positive
regu
lation of in
tracel
lular m
RNA local
izatio
n
exter
nal sid
e plas
ma mem
brane
side o
f mem
brane
dimeri
c IgA
immunoglo
bulin co
mplex
secret
ory dim
eric I
gA im
munoglobulin
complex
IgA im
munoglobulin co
mplex
IgA im
munoglobulin co
mplex, ci
rculat
ing
monomeric I
gA im
munoglobulin
complex
polymeri
c IgA
immunoglo
bulin co
mplex
secret
ory IgA
immunoglo
bulin co
mplex
olfacto
ry rec
eptor a
ctivit
y
cytokine b
inding
C-C ch
emokine b
inding
chem
okine binding
G-protei
n coupled
chem
oattrac
tant r
eceptor a
ctivit
y
chem
okine rece
ptor acti
vity
C-C ch
emokine r
eceptor a
ctivit
y
thiol-dep
enden
t ubiquiti
nyl hyd
rolas
e acti
vity
ubiquitinyl
hydro
lase a
ctivit
y
ubiquitin-li
ke pro
tein-sp
ecific p
rotea
se act
ivity
(a)
Num
bers
of g
enes
15
10
5
0
Gene Function Classification (KEGG)
Olfacto
ry tr
ansd
uctio
n
Bile
secre
tion
Cyto
kine-c
ytokin
e rec
epto
r int
erac
tion
Matu
rity o
nset
diab
etes o
f the
youn
g
Chem
okin
e sign
aling
path
way
(b)
Figure 2: GO analysis and KEGG pathway assessment on mRNAs. (a) GO analysis on mRNAs. (b) KEGG pathway analysis on mRNAs.
was regulated at the mRNA level [38]. Certain mRNAsassociated with lung inflammation (including that encodingIL-8) were appeared to be a regulator of other miRNAsbiogenesis in lung epithelial cells, showing they were poten-tial candidates for anti-inflammatory therapeutics for cysticfibrosis [39, 40]. We found that mRNAs encoding LOYBOY,DEEZOY, and SNEYSPORBY were upregulated, and thoseencoding CPA3, FCER1A, and CCR5 were downregulated;these mRNAsmay be involved in silica-induced lung fibrosis.Although lncRNAs may be potential therapeutic targets infibrosis [41], the only relevant report is that of Sun et al.;in the paraquat-induced mouse lung fibrosis model, 513lncRNAs were upregulated and 204 were downregulated,which affected cell differentiation, epithelial morphogenesis,and immune response, all of which are closely associatedwith the epidermal-mesenchymal transition (EMT) [24].Our study identified 366 differentially expressed lncRNAs inhuman sample in which the numbers were lower than Sun etal. study. This might be due to the species variation.
Go and KEGG analyses are commonly used to pre-dict biological functions and/or signaling pathways affectedby mRNAs and lncRNAs. On GO analysis of significantchanges in mRNAs (fold change >2; p<0.05), in termsof biological proses, cellular components, and molecularfunction, the most-affected processes were interleukin-12
secretion, regulation of intracellular mRNA localization,effects on the dimeric IgA immunoglobulin complex, C-C chemokine binding, inflammation response, immuneresponse, ubiquitin-like protein-specific protease activity,and thiol-dependent ubiquitinyl hydrolase activity. Here, theresults which the inflammation response affected by mRNAsare consistent with the previous report; however, mRNAsaffecting the ubiquitination and immune response are newinsights into the mRNAs function in the lung fibrosis. Thissuggests mRNAsmight participate in lung fibrosis process byregulating immune response and ubiquitination. Ubiquitinis a highly conserved low-molecular-weight protein andubiquitination represents a very common posttranslationalmodification. Extracellular ubiquitin regulates the immuneresponse and exhibits anti-inflammatory and neuroprotec-tive activities. Zhou et al. reported that silica activatedmacrophages and increased circRNA HECTD1 levels viaubiquitination [42]. Tsubouchi et al. suggested that, in ableomycin-induced lung fibrosis mouse model, azithromycin(AZM; a second-generation antibacterial macrolide) inhib-ited autophagy via ubiquitination of NOX4 [43]. Nan et al.found that ubiquitination of the transcription factors Smad2and Smad3 was changed via the action of ubiquitin carboxyl-terminal hydrolase-L5 during the pathogenesis of idiopathicpulmonary fibrosis [44]. However, the role of ubiquitination
8 BioMed Research International
Gene Function Classification (GO)N
umbe
rs o
f gen
es
80
60
40
20
0
Biological ProcessCellular ComponentMolecular Function
detecti
on of
chem
ical st
imulus i
nvolve
d in se
nsory
percep
tion
detecti
on of
chem
ical st
imulus i
nvolve
d in se
nsory
percep
tion of
smell
senso
ry perc
eptio
n of sm
ell
kerat
inizatio
n
Keratin
ocyte
differe
ntiatio
n
epiderm
al cel
l diff
erentia
tion
epiderm
is dev
elopmen
t
skin dev
elopmen
t
neuro
peptid
e sign
aling p
athway
Spon
glotro
7phob
last la
yer d
evelo
pment
phospholi
pase C
-activ
ating G
-pro
tein co
upled re
ceptor
sign
aling p
athway
Inter
med
iate fi
lamen
t
Inter
med
iate fi
lamen
t cyto
skele
ton
kerat
in filamen
t
comfied
enve
lope
olfac
tory r
eceptor
activ
ity
odoran
t binding
G-Pro
tein co
upled pep
tide r
eceptor
activ
ity
peptid
e rece
ptor ac
tivity
hormon
e acti
vity
G-pro
tein co
upled am
ine rece
ptor ac
tivity
transcr
iption fa
ctor a
ctivit
y RNA poly
merase
II dist
al en
hancer
sequ
ence-
speci
fic binding
C-C ch
emok
ine rece
ptor ac
tivity
volta
ge-ga
ted pota
ssium ch
annel
activ
ity
potassi
um chan
nel ac
tivity
protei
n kinase
inhibi
tor ac
tivity
transcr
iptionala
ctiva
tor ac
tivity
RNA polymera
se II
transcr
iption re
gulat
ory r
egion
sequ
ence-
speci
fic binding
(a)
Num
bers
of g
enes
80
60
40
20
0
Gene Function Classification (KEGG)
Olfac
tory
tran
sduc
tion
Neur
oacti
ve lig
and-
rece
ptor
inter
actio
nTa
ste tr
ansd
uctio
n
Estro
gen
signa
ling p
athwa
yCh
olin
ergic
syna
pse
Cyto
kine
-cyt
okin
e rec
epto
r int
erac
tion
Reni
n se
creti
onHe
patit
is C
Syste
mic
lupu
s ery
them
atosu
s
Chem
okin
e sign
aling
path
way
Long
-term
pot
entia
tion
(b)
Figure 3: GO analysis and KEGG pathway assessment on lncRNAs. (a) GO analysis on lncRNAs. (b) KEGG pathway analysis on lncRNAs.
in the development of SIPF remains poorly understood andfurther studies are thus required. GO analysis indicatedthat both mRNAs and lncRNAs participate in the immuneresponse in lung fibrosis. Balloy V. et al. assessed transcrip-tomes of lncRNA between cystic fibrosis and found a specificcystic fibrosis signature of lncRNA expression was associatedwith immune response of patients [45]. A study by SookoianS. et al. found that lncRNA MALAT1 participated in thedevelopment of liver fibrosis [46]. However, studies regardingthe underlying mechanism of lncRNAs participating in thedevelopment of silica-induced lung fibrosis remain furtherinvestigated.
Liu et al. created an lncRNA/mRNA coexpression net-work for human pulmonary epithelial cells; 24 mRNAswere coexpressed and interacted with 33 lncRNAs; eachmRNA interacted with a few lncRNAs and each lncRNAwith a few mRNAs [47]. We included RPL41P3, RP11-609D21.39, RP4-620F22.1, AP001610.5, and RN7SL783P in
a coexpression network with related mRNAs. For RPL41P3,GO analysis revealed involvement in chemokine binding,CoA-transferase activity, and G-protein-coupled serotoninreceptor activity, all of which are associated with SIPF devel-opment (Figure S1); lncRNA-RP11-609D21.3 was principallyassociated with carbohydrate transporter, endoribonuclease,and ribonuclease A activities, and with spliceosome activ-ity, ferroptosis, and Kaposi’s sarcoma-associated herpesvirusinfection (Figure S2). The lncRNARP4-620F22.2 was asso-ciated with sphinganine-1-phosphate metabolism, penetra-tion of the zona pellucida, diol and sphingoid metabolism,sperm fusion to the egg plasma membrane, and viralRNA genome replication. Cell component analysis indi-cated that this lncRNA affected the exterior and marginalplasma membrane; the molecular functions involved weredihydrosphingosine-1-phosphate phosphatase, sphingosine-1-phosphate phosphatase, and hyaluronoglucosaminidase, aswell as sequence-specific single-stranded DNA binding and
BioMed Research International 9
hexosaminidase (Figure S3). The biological processes inwhich Metazoa-SRP lncRNA was involved included type Iinterferon signaling pathway, the cellular response to typeI interferon, negative regulation of helicase activity, theviral defense response and its negative regulation, intra-cellular viral protein transport, intracellular symbiont pro-tein transport, and cellular responses to interferon-𝛼 andexogenous dsRNA.The KEGG data showed that the lncRNAparticipated in multiple pathways, including the herpessimplex infection pathway, a hepatitis C pathway involvinginflammation, and the PPAR signaling pathway (FigureS4). The lncRNA RN7SL783P molecular functions includedmelanocortin receptor binding and adenylate cyclase andphosphorus-oxygen lyase. The KEGG data showed that thelncRNAwas active inmelanogenesis, estrogen signaling, gas-tric acid secretion, GnRH signaling, circadian entrainment,aldosterone synthesis, and secretion, as well as inflammatoryregulation of TRP channels, glucagon signaling, vascularsmooth muscle contraction, oocyte meiosis apelin signaling,adrenergic signaling in cardiomyocytes, oxytocin signaling,cGMP-PKG signaling, calcium signaling, cAMP signaling,rap1 signaling, and phototransduction (Figure S5).
Our work provides an initial understanding of the rolesplayed by lncRNAs in SIPF.The lncRNA levels determined viaquantitative reverse transcription (qRT-PCR)were consistentwith the microarray data. Our results suggest that (1) silicaexposure modulates the expression of both lncRNAs andgenes and (2) some silica-induced genes may be involved inubiquitination, and some lncRNAs may affect immune cellinteractions. The mechanisms by which lncRNA and mRNAcoexpression mediates SIPF require further study.
However, our study had certain limitations. First, thelncRNAs measured via qRT-PCR constituted only a smallproportion of those identified via microarray; some lncRNAsmay have been missed such that further studies are required.Second, we recruited only five cases; studies with largersamples are required. However, we are the first to use humanperipheral blood samples to identify lncRNAs andmRNAs inpatients with SIPF.
5. Conclusions
Gene and lncRNA expression profiles were dysregulated inthe peripheral blood of SIPF patients, suggesting future diag-nostic and therapeutic biomarkers. LncRNA dysregulationmay play an important role in disease development, given itsplethora of effects on posttranscriptional mRNA regulation.However, further biology function work is required.
Data Availability
The data used to support the findings of this study areincluded within the supplementary information file.
Additional Points
Key Messages. What is already known about this subject?Silica exposure in the workplace is the potential risk factor for
workers, which has been identified to be the cause of silicosis,leading to damage workers health, with many complicationsincluding tuberculosis. Preliminary studies suggest that earlydiscovery and diagnosis of biomarkers would benefit fromproviding interventions on improving patients’ life quality.What are the new findings? To the best of our knowledge,this study focusing on the uncovering of the transcriptomessignature of patients with silicosis is firstly reported. Previoussimilar studies were conducted using animal model or cellmodel rather than using peripheral blood samples. In total,1069 differentially expressed mRNAs and 366 lncRNAs wereidentified. Four genes named FEM1B, TRIM39, TRIM32, andKLHL15 were enriched significantly with ubiquitination andimmune response. How might it impact clinical practicein the foreseeable future? This study indicated that theselncRNAs andmRNAsmay play a critical biological role in theprogression and development of silica-induced silicosis andare thus candidate therapeutic targets and potential diagnosisbiomarkers.Meanwhile, this transcriptomes studymight pro-videmore baseline understanding themolecularmechanismsof silicosis which could pave the way for a further functionalstudy of underlying regulatory mechanisms.
Ethical Approval
The study was performed in accordance with all relevanttenets of the Declaration of Helsinki and was approved bythe Ethics Committee of Xiangya School of Public Health(Approval no. XYGW-2018-11).
Consent
This work is published under the standard license to publishagreement. All participants provided written informed con-sent.
Conflicts of Interest
The authors declare no conflicts of interest.
Authors’ Contributions
Ping-Kun Zhou and Rui-Xue Huang contributed study con-cept and critical design and Xiao-Dan Liu, Li He, Yao Zhou,and Hua Guan conducted the cell experiments. Hua Guan,Li He, Yao Zhou, and Xiao-Dan Liu acquired, analyzed, andinterpreted data. Hua Guan, Li He, Yao Zhou, and Xiao-DanLiu facilitated quality assurance of data. Hua Guan, Xiao-Dan Liu, Li He, Yao Zhou, and Ping-Kun Zhou fulfill theinitial manuscript and Ping-Kun Zhou and Rui-Xue Huangcritically reviewed and revised the final manuscript. Allauthors were involved in the current study and contributedto this study. Yao Zhou and Li He contributed equally to thestudy.
Acknowledgments
This study is supported by grants from National Key BasicResearch Program (973 Program) ofMOST, China (Grant no.
10 BioMed Research International
2015CB910601), the National Natural Science Foundations ofChina (Grant nos. 81530085, U1432248, and U1803124), andthe Fundamental Research Funds for the Central Universitiesof Central South University (2018zzts849).
Supplementary Materials
Figure S1. GO analysis and KEGG pathway assessment onlncRNA RP11-138I1.4. A. GO analysis on lncRNA RP11-138I1.4. B. KEGG pathway analysis on lncRNA RP11-138I1.4.Figure S2. GO analysis on lncRNA RPL41p3. Figure S3.GO analysis on 18S F. Figure S4. GO analysis and KEGGpathway assessment on lncRNA AP001610.5. A. GO anal-ysis on lncRNA AP001610.5. B. KEGG pathway analysison lncRNA AP001610.5. Figure S5. GO analysis and KEGGpathway assessment on lncRNA RN7SL783P. A. GO anal-ysis on lncRNA RN7SL783P. B. KEGG pathway analysison lncRNA RN7SL783P. Figure S6. lncRNA-mRNA coex-pression network. Figure S7. lncRNA-mRNA coexpressionnetwork of AP001610.5, RP4-620F22.2, RP11-605F22.1, andRP11-609D21.3, respectively. Table S1. Primer sequences usedin Quantitative Real-Time PCR (qRT-PCR). Table S2. Fivepathways and the corresponding affected genes in case groupcompared with control group. (Supplementary Materials)
References
[1] K. M. Pollard, “Silica, silicosis, and autoimmunity,” Frontiers inImmunology, vol. 7, no. 97, 2016.
[2] J. Li, W. Yao, J. Y. Hou et al. et al., “Crystalline silica promotesrat fibrocyte differentiation in vitro, and fibrocytes participatein silicosis in vivo,” Biomedical and Environmental Sciences, vol.30, no. 9, pp. 649–660, 2017.
[3] A. D. Siribaddana, K.Wickramasekera,W.M. Palipana et al., “Astudy on silicosis among employees of a silica processing factoryin the Central Province of Sri Lanka,” �e Ceylon MedicalJournal, vol. 61, no. 1, pp. 6–10, 2016.
[4] D. F. Cao, Y. Jin, and B. Shang, “The analysis of the silicosismisdiagnosis in the silica bricks industry,” Chinese Journal ofIndustrial Hygiene and Occupational Diseases, vol. 34, no. 3, pp.211-212, 2016.
[5] P.-A. Rosental, D. Rosner, and P. D. Blanc, “From silicosis tosilica hazards: An experiment in medicine, history, and thesocial sciences,”American Journal of IndustrialMedicine, vol. 58,pp. 3–5, 2015.
[6] F. Cainelli, M. N. Tanko, and S. Vento, “Silica exposure andsilicosis: Action is needed now,” Southern Medical Journal, vol.103, no. 11, p. 1078, 2010.
[7] H.Ning, Y. Zhou, Z. Zhou, S. Cheng, andR.Huang, “Challengesto improving occupational health in China,” Occupational andEnvironmental Medicine, vol. 74, no. 12, pp. 924-925, 2017.
[8] R. Ehrlich, A. Montgomery, P. Akugizibwe, and G. Gonsalves,“Public health implications of changing patterns of recruitmentinto the South African mining industry, 1973-2012: a databaseanalysis,” BMC Public Health, vol. 18, no. 1, Article ID 93, 2017.
[9] N. L. Bragazzi, G. Dini, A. Toletone, F. Brigo, and P. Durando,“Infodemiological data concerning silicosis in the USA in theperiod 2004–2010 correlating with real-world statistical data,”Data in Brief, vol. 10, pp. 457–464, 2017.
[10] T. Kauppinen, J. Toikkanen, D. Pedersen et al. et al., “Occu-pational exposure to carcinogens in the european union,”Occupational and Environmental Medicine, vol. 57, no. 1, pp. 10–18, 2000.
[11] C. M. Barber, D. Fishwick, M. J. Seed, M. Carder, and M.van Tongeren, “Artificial stone-associated silicosis in the UK,”Occupational and Environmental Medicine, vol. 75, no. 7, pp.541.1-541, 2018.
[12] X. Wang and E. Yano, “Pulmonary dysfunction in silica-exposed workers: A relationship to radiographic signs of silico-sis and emphysema,” American Journal of Industrial Medicine,vol. 36, no. 2, pp. 299–306, 1999.
[13] G. K. Kulkarni, “Prevention and control of silicosis: A nationalchallenge,” Indian Journal of Occupational and EnvironmentalMedicine, vol. 11, no. 3, pp. 95-96, 2007.
[14] C. Xue, N. Wu, X. Li, M. Qiu, X. Du, and Q. Ye, “Serumconcentrations of Krebs von den Lungen-6, surfactant proteinD, and matrix metalloproteinase-2 as diagnostic biomarkersin patients with asbestosis and silicosis: A case-control study,”BMC Pulmonary Medicine, vol. 17, no. 1, Article ID 144, 2017.
[15] A. M. K. Choi and J. Alam, “Heme oxygenase-1: function,regulation, and implication of a novel stress-inducible protein inoxidant-induced lung injury,” American Journal of RespiratoryCell and Molecular Biology, vol. 15, no. 1, pp. 9–19, 1996.
[16] S. Fati and L. Aloia, “The behavior of ceruloplasmin and ofserum copper in pulmonary silicosis and in other diseases ofthe respiratory tract,” FoliaMedica FoliaMedica, vol. 44, pp. 731–751, 1961.
[17] J. Pandey and D. Agarwal, “Biomarkers: A potential prognostictool for silicosis,” Indian Journal of Occupational and Environ-mental Medicine, vol. 16, no. 3, pp. 101–107, 2012.
[18] X. Zhang, J. Hou, and M. Lu, “Regulation of hepatitis Bvirus replication by epigenetic mechanisms and microRNAs,”Frontiers in Genetics, vol. 4, Article ID 202, 2013.
[19] S. X. Jin, Q. Y. Wu, W. W. Yan, and C. H. Ni, “Therapeuticeffect of miR-489 in a mouse model of silica-induced maturedpulmonary fibrosis,” Chinese Journal of Industrial Hygiene andOccupational Diseases, vol. 35, no. 5, pp. 337–341, 2017.
[20] R. Han, X. Ji, R. Rong et al., “MiR-449a regulates autophagyto inhibit silica-induced pulmonary fibrosis through targetingBcl2,” Journal of Molecular Medicine, vol. 94, no. 11, pp. 1267–1279, 2016.
[21] Z. Yang, Q. Li, S. Yao et al., “Down-Regulation of miR-19a as aBiomarker for Early Detection of Silicosis,” Anatomical Record,vol. 299, no. 9, pp. 1300–1307, 2016.
[22] Y. Zhang, D. Zhou, F. Wang et al., “Bronchoalveolar lavage fluidmicroRNA-146a a biomarker of disease severity and pulmonaryfunction in patients with silicosis,” Journal of Occupational andEnvironmental Medicine, vol. 58, no. 5, pp. e177–e182, 2016.
[23] V. A. Moran, R. J. Perera, and A. M. Khalil, “Emerging func-tional and mechanistic paradigms of mammalian long non-coding RNAs,” Nucleic Acids Research, vol. 40, no. 14, pp. 6391–6400, 2012.
[24] H. Sun, J. Chen, W. Qian et al., “Integrated long non-coding RNA analyses identify novel regulators of epithelial-mesenchymal transition in the mouse model of pulmonaryfibrosis,” Journal of Cellular andMolecularMedicine, vol. 20, no.7, pp. 1234–1246, 2016.
[25] J.-H. Yuan, F. Yang, F. Wang et al., “A Long Noncoding RNAActivated by TGF-𝛽 promotes the invasion-metastasis cascadein hepatocellular carcinoma,”Cancer Cell, vol. 25, no. 5, pp. 666–681, 2014.
BioMed Research International 11
[26] H. Chen, K. Zhang, J. Lu, G. Wu, H. Yang, and K. Chen,“Comprehensive analysis ofmRNA-lncRNA co-expression pro-file revealing crucial role of imprinted gene clusterDLK1-MEG3in chordoma,”Oncotarget , vol. 8, no. 68, pp. 112623–112635, 2017.
[27] P. J. McKiernan, O. Cunningham, C. M. Greene, and S. A.Cryan, “Targeting miRNA-based medicines to cystic fibrosisairway epithelial cells using nanotechnology,” InternationalJournal of Nanomedicine, vol. 8, pp. 3907–3915, 2013.
[28] I. K. Oglesby, R. Agrawal, M. A. Mall, N. G. McElvaney, andC. M. Greene, “miRNA-221 is elevated in cystic fibrosis airwayepithelial cells and regulates expression ofATF6,”Molecular andCellular Pediatrics, vol. 2, no. 1, 2015.
[29] I. K. Oglesby and P. J. McKiernan, “MiRNA expression in cysticfibrosis bronchial epithelial cells,”Methods inMolecularBiology,vol. 1509, pp. 57–69, 2017.
[30] B. D. Tong, M. Y. Xiao, and J. X. Zeng, “MiRNA-21 promotesfibrosis in orbital fibroblasts from thyroid-associated ophthal-mopathy,”Molecular Vision, vol. 21, pp. 324–334, 2015.
[31] R. Huang, T. Yu, Y. Li, and J. Hu, “Upregulated has-miR-4516 as a potential biomarker for early diagnosis of dust-induced pulmonary fibrosis in patients with pneumoconiosis,”Toxicology Research, vol. 7, no. 3, pp. 415–422, 2018.
[32] L. Mo, M. Song, Q. Huang et al., “Exosome-packaged miR-1246 contributes to bystander DNA damage by targeting LIG4,”British Journal of Cancer, vol. 119, no. 4, pp. 492–502, 2018.
[33] X. Liu, D. Xie, Y. Wang, H. Guan, R. Huang, and P. Zhou,“Integrated analysis of lncRNA–mRNAco-expression networksin the 𝛼-particle induced carcinogenesis of human branchialepithelial cells,” International Journal of Radiation Biology, pp.1–12, 2018.
[34] J. Qin, H. Ning, Y. Zhou, Y. Hu, L. Yang, and R. Huang,“LncRNA MIR31HG overexpression serves as poor prognosticbiomarker and promotes cells proliferation in lung adenocarci-noma,” Biomedicine & Pharmacotherapy, vol. 99, pp. 363–368,2018.
[35] P. K. Zhou and R. X. Huang, “Targeting of the respiratorychain by toxicants: beyond the toxicities to mitochondrialmorphology,” Toxicology Research, vol. 7, no. 6, pp. 1008–1011,2018.
[36] J. Qin,H.Ning, Y. Zhou et al., “LncRNAUc.173 is a keymoleculefor the regulation of lead-induced renal tubular epithelial cellapoptosis,” Biomedicine & Pharmacotherapy, vol. 100, pp. 101–107, 2018.
[37] L. Plantier, H. Renaud, R. Respaud, S. Marchand-Adam, andB. Crestani, “Transcriptome of cultured lung fibroblasts inidiopathic pulmonary fibrosis:Meta-analysis of publically avail-able microarray datasets reveals repression of inflammationand immunity pathways,” International Journal of MolecularSciences, vol. 17, no. 12, 2016.
[38] R. Biswas, P. Kumar, and H. B. Pollard, “Regulation of mRNAturnover in cystic fibrosis lung disease,”Wiley InterdisciplinaryReviews: RNA, vol. 8, no. 4, 2017.
[39] S. Bhattacharyya, P. Kumar, M. Tsuchiya, A. Bhattacharyya,and R. Biswas, “Regulation of miR-155 biogenesis in cysticfibrosis lung epithelial cells: Antagonistic role of two mRNA-destabilizing proteins, KSRP and TTP,” Biochemical and Bio-physical Research Communications, vol. 433, no. 4, pp. 484–488,2013.
[40] S. Bhattacharyya, U. Gutti, J. Mercado, C. Moore, H. B. Pol-lard, and R. Biswas, “MAPK signaling pathways regulate IL-8mRNA stability and IL-8 protein expression in cystic fibrosis
lung epithelial cell lines,” American Journal of Physiology-LungCellular and Molecular Physiology, vol. 300, no. 1, pp. L81–L87,2011.
[41] X. Jiang and F. Zhang, “Long noncoding RNA: A new contrib-utor and potential therapeutic target in fibrosis,” Epigenomics,vol. 9, no. 9, pp. 1233–1241, 2017.
[42] Z. Zhou, R. Jiang, X. Yang et al., “circRNA mediates silica-induced macrophage activation via HECTD1/ZC3H12A-dependent ubiquitination,” �eranostics, vol. 8, no. 2, pp.575–592, 2018.
[43] K. Tsubouchi, J. Araya, S. Minagawa et al., “Azithromycin atten-uates myofibroblast differentiation and lung fibrosis develop-ment through proteasomal degradation of NOX4,” Autophagy,vol. 13, no. 8, pp. 1420–1434, 2017.
[44] L. Nan, A. M. Jacko, J. Tan et al., “Ubiquitin carboxyl-terminalhydrolase-L5 promotes TGF𝛽-1 signaling by de-ubiquitinatingand stabilizing Smad2/Smad3 in pulmonary fibrosis,” ScientificReports, vol. 6, Article ID 33116, 2016.
[45] V. Balloy, R. Koshy, L. Perra et al., “Bronchial epithelial cellsfrom cystic fibrosis patients express a specific long non-codingRNA signature uponPseudomonas aeruginosa infection,” Fron-tiers in Cellular and InfectionMicrobiology, vol. 7, Article ID 218,2017.
[46] S. Sookoian, D. Flichman, M. E. Garaycoechea, J. San Martino,G. O. Castano, and C. J. Pirola, “Metastasis-associated lungadenocarcinoma transcript 1 as a common molecular driver inthe pathogenesis of nonalcoholic steatohepatitis and chronicimmune-mediated liver damage,”Hepatology Communications,vol. 2, no. 6, pp. 654–665, 2018.
[47] H. Liu, X. Zhao, J. Xiang et al., “Interaction network ofcoexpressedmRNA, miRNA, and lncRNA activated by TGF-𝛽1regulates EMT in human pulmonary epithelial cell,”MolecularMedicine Reports, vol. 16, no. 6, pp. 8045–8054, 2017.
Hindawiwww.hindawi.com
International Journal of
Volume 2018
Zoology
Hindawiwww.hindawi.com Volume 2018
Anatomy Research International
PeptidesInternational Journal of
Hindawiwww.hindawi.com Volume 2018
Hindawiwww.hindawi.com Volume 2018
Journal of Parasitology Research
GenomicsInternational Journal of
Hindawiwww.hindawi.com Volume 2018
Hindawi Publishing Corporation http://www.hindawi.com Volume 2013Hindawiwww.hindawi.com
The Scientific World Journal
Volume 2018
Hindawiwww.hindawi.com Volume 2018
BioinformaticsAdvances in
Marine BiologyJournal of
Hindawiwww.hindawi.com Volume 2018
Hindawiwww.hindawi.com Volume 2018
Neuroscience Journal
Hindawiwww.hindawi.com Volume 2018
BioMed Research International
Cell BiologyInternational Journal of
Hindawiwww.hindawi.com Volume 2018
Hindawiwww.hindawi.com Volume 2018
Biochemistry Research International
ArchaeaHindawiwww.hindawi.com Volume 2018
Hindawiwww.hindawi.com Volume 2018
Genetics Research International
Hindawiwww.hindawi.com Volume 2018
Advances in
Virolog y Stem Cells International
Hindawiwww.hindawi.com Volume 2018
Hindawiwww.hindawi.com Volume 2018
Enzyme Research
Hindawiwww.hindawi.com Volume 2018
International Journal of
MicrobiologyHindawiwww.hindawi.com
Nucleic AcidsJournal of
Volume 2018
Submit your manuscripts atwww.hindawi.com