- Research article
- Open Access
Differentially expressed non-coding RNAs induced by transmissible gastroenteritis virus potentially regulate inflammation and NF-κB pathway in porcine intestinal epithelial cell line
BMC Genomics volume 19, Article number: 747 (2018)
Transmissible gastroenteritis virus (TGEV) infection can activate NF-κB pathway in porcine intestinal epithelial cells and result in severe inflammation. Non-coding RNAs (ncRNAs) are not translated into proteins and play an important role in many biological and pathological processes such as inflammation, viral infection, and mitochondrial damage. However, whether ncRNAs participate in TGEV-induced inflammation in porcine intestinal epithelial cells is largely unknown.
In this study, the next-generation sequencing (NGS) technology was used to analyze the profiles of mRNAs, miRNAs, and circRNAs in Mock- and TGEV-infected intestinal porcine epithelial cell-jejunum 2 (IPEC-J2) cell line. A total of 523 mRNAs, 65 microRNAs (miRNAs), and 123 circular RNAs (circRNAs) were differentially expressed. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed differentially expressed mRNAs were linked to inflammation-related pathways, including NF-κB, Toll-like receptor, NOD-like receptor, Jak-STAT, TNF, and RIG-I-like receptor pathways. The interactions among mRNA, miRNA, and circRNA were analyzed. The data showed that ssc_circ_009380 and miR-22 might have interaction relationship. Dual-luciferase reporter assay confirmed that miR-22 directly bound to ssc_circ_009380. We also observed that overexpression of miR-22 led to a reduction of p-IκB-α and accumulation of p65 in nucleus in TGEV-infected IPEC-J2 cells. In contrast, inhibition of miR-22 had the opposite effects. Moreover, silencing of ssc_circ_009380 inhibited accumulation of p65 in nucleus and phosphorylation of IκB-α.
The data revealed that differentially expressed mRNAs and ncRNAs were primarily enriched in inflammation-related pathways and ssc_circ_009380 promoted activation of NF-κB pathway by binding miR-22 during TGEV-induced inflammation.
Transmissible gastroenteritis (TGE) is an acute viral disease, which is caused by TGEV infection and characterized by vomiting, dehydration, and severe diarrhea in pigs of all ages, especially less than 2-week-old suckling piglets. TGEV can impair porcine intestinal epithelial cells and trigger inflammatory response, mitochondrial injury, apoptosis, and complete mitophagy in IPEC-J2 . Not only TGEV, but also other coronavirus such as Severe acute respiratory syndrome (SARS), Middle east respiratory syndrome coronavirus (MERS), Porcine epidemic diarrhea virus (PEDV) can cause severe inflammation response [2,3,4,5], indicating that inflammation response is a common pathological process during coronaviruses infection.
Because of not coding, ncRNAs were previously known as junk RNAs. However, increasing evidence shows that ncRNAs play important regulatory roles in many biological and pathological processes, such as inflammation, viral infection, and apoptosis [6,7,8,9,10,11,12,13,14,15]. We previously reported that TGEV infection led to 21 differentially expressed miRNAs, among which miR-4331 inhibited transcription of TGEV gene7 via targeting CDCA7 , and miR-27b functioned as a negative regulator during TGEV-induced apoptosis . Nevertheless, little is known about whether circRNAs are involved in TGEV-induced inflammation in porcine intestinal epithelial cells.
Except for directly regulating genes, circRNAs may act as regulators by competing together with miRNAs to affect the stability of target RNAs or their translation. circRNAs harboring miRNA response elements (MREs) can function as miRNAs sponge to regulate gene expression. To investigate the effects of circRNAs during TGEV-induced inflammation, we performed RNA sequencing using NGS technology to test the alteration of mRNA, miRNA, and circRNA profiles during TGEV-induced inflammation in IPEC-J2 and conducted bioinformatics analysis of differentially expressed mRNAs, miRNAs, and circRNAs. We demonstrated that differentially expressed mRNAs, miRNAs, and circRNAs were mainly enriched in inflammation and immune response. ssc_circ_009380 attenuated TGEV-induced activation of NF-κB pathway via binding miR-22.
TGEV infection induces inflammatory response in IPEC-J2
IPEC-J2 cells were infected with TGEV at 1.0 MOI for 12 h, 24 h, and 36 h. The mRNA levels of IL-8, IL-6, and TNF-α were examined. The results displayed that mRNA levels of IL-8, IL-6, and TNF-α increased and reached to the peak at 24 h post infection (hpi) (Fig. 1).
Sequencing and analysis of ncRNAs and mRNAs of TGEV-infected IPEC-J2
The total RNA was extracted and sequenced using the NGS technology. In total, 523 differentially expressed mRNAs (including 462 up-regulated mRNAs and 61 down-regulated mRNAs), 65 differentially expressed miRNAs (including 46 up-regulated miRNAs and 19 down-regulated miRNAs) and 123 differentially expressed circRNAs (including 69 up-regulated circRNAs and 54 down-regulated circRNAs) were obtained (Fig. 2). In addition, there were 45 novel miRNAs in the differentially expressed miRNAs. More detailed information is respectively presented in Additional file 1: Tables S1–S3.
KEGG analysis of differentially expressed miRNAs, circRNAs, and mRNAs
We predicted the targets of differentially expressed miRNAs and source genes of differentially expressed circRNAs. 7571 target genes of miRNAs and 98 source genes of circRNAs were obtained (Additional file 1: Tables S4 and S5). These genes and differentially expressed mRNAs were respectively searched for functional enrichments by a KEGG database search (http://www.genome.jp/kegg/). KEGG analysis results indicated that differentially expressed mRNAs were mostly involved in Toll-like receptor, RIG-I-like receptor, TNF, NOD-like receptor, Jak-STAT, and NF-κB pathways (Fig. 3a). Differentially expressed miRNAs were primarily enriched in B cell receptor and Toll-like receptor pathways (Fig. 3b). The source genes of differentially expressed circRNAs were involved in RIG-I-like receptor, TNF, NOD-like receptor, and NF-κB pathways (Fig. 3c).
Analysis of ncRNA-mRNA regulatory network
To predict the interaction between ncRNAs and mRNAs during TGEV infection, the intersection of differentially expressed targets of miRNAs and mRNAs were collected. Then, 133 genes were obtained (Fig. 4a). The interaction network of the 133 genes was constructed (Fig. 4b) (More detailed information is shown in Additional file 1: Table S6). circRNA-miRNA-mRNA regulation network was generated (Fig. 5a and Additional file 1: Table S7). The mRNAs in circRNA-miRNA-mRNA regulatory network were searched against KEGG database for pathway enrichment. The results indicated that these mRNAs were most significantly involved in RIG-I-like receptor, TNF, NOD-like receptor, Toll-like receptor, and NF-κB pathways (Fig. 5b).
Validation of ncRNAs and mRNAs by qRT-PCR
Both divergent primers and convergent primers were designed and synthesized for identification of circular form. Head-to-tail splicing was validated by PCR. The results showed that ssc_circ_001964, ssc_circ_000470, ssc_circ_005884, and ssc_circ_009380 are circular (Fig. 6a). To validate the reliability of RNA sequencing, the expression levels of differentially expressed miRNAs (miR-218, mir-2413, mir-492, mir-7550, mir-2510, and miR-22), mRNAs (DDX58, CCL5, and IL-6), and circRNAs (ssc_circ_001964, ssc_circ_000470, ssc_circ_005884, and ssc_circ_009380) were measured by qRT-PCR. The results were consistent with that of RNA sequencing (Fig. 6b).
miR-22 directly binds to ssc_circ_009380
In the present study, we demonstrated that TGEV infection activated inflammation-related pathways and up-regulated miR-22. It was reported that miR-22 was involved in regulating inflammation and immune response [18, 19]. Therefore, we speculated that miR-22 might play a role in TGEV-induced inflammation response. Based on bioinformatics analysis, miR-22 might bind to ssc_circ_009380. To confirm that, the sequence of ssc_circ_009380 was amplified by PCR and cloned into 3′ UTR of Renilla luciferase in psiCHECK-2 to construct wild-type (WT) recombinant plasmid psi-ssc_circ_009380-WT. The binding sites of miR-22 in ssc_circ_009380 were mutated by point mutation to construct mutant recombinant plasmid psi-ssc_circ_009380-Mut (Fig. 7a and b). The recombinant plasmids were respectively co-transfected into IPEC-J2 cells together with miR-22 mimics (or miRNA mimics control, miR-22 inhibitors, miRNA inhibitors control). miR-22 was overexpressed using miR-22 mimics and inhibited using miR-22 inhibitors (Fig. 7c and d). The luciferase activities of psi-ssc_circ_009380-WT was decreased at 62% by miR-22 mimics and not affected by miR-22 inhibitors (Fig. 7e and f). However, the luciferase activities of psi-ssc_circ_009380-Mut was not affected by miR-22 mimics and inhibitors (Fig. 7e and f).The effects of miR-22 and ssc_circ_009380 on TGEV-induced activation of NF-κB in IPEC-J2.
To evaluate the effect of miR-22 on TGEV-induced NF-κB activation, IPEC-J2 cells were transfected with miR-22 mimics (or miRNA mimics control, miR-22 inhibitors, miRNA inhibitors control) and subsequently infected with TGEV at 1 MOI for 24 h. The miR-22 was overexpressed by miR-22 mimics and inhibited by miR-22 inhibitors (Fig. 8a and b). p-IκB-α and nucleic distribution of p65 were suppressed by miR-22 mimics and increased by miR-22 inhibitors, but not affected by miR-22 mimics and inhibitors in cytoplasm (Fig. 8d and e).
To evaluate the effect of ssc_circ_009380 on TGEV-induced activation of NF-κB, IPEC-J2 cells were transfected with siRNA of ssc_circ_009380 (siCirc009380) (or negative control) and subsequently infected with TGEV at 1 MOI for 24 h. The ssc_circ_009380 level was down-regulated by siCirc009380 (Fig. 8c). p-IκB-α and p65 in nucleus were decreased by siCirc009380 and not affected in cytoplasm (Fig. 8f).
In the present study, differentially expressed 523 mRNAs, 65 miRNAs, and 123 circRNAs were obtained. KEGG analysis showed that differentially expressed mRNAs and ncRNAs were primarily enriched in inflammation-related pathways. circRNA ssc_circ_009380 was identified as the sponge of miR-22 and enhanced activation of NF-κB pathway through binding miR-22 during TGEV infection.
During viral infection, host pathogen recognition receptors (PRRs), including Toll-like receptor, RIG-I-like receptor, NOD-like receptor, are responsible for recognizing viruses. Then, the recognition initiates a series of signaling cascades, including Jak-STAT, TNF, and NF-κB pathways to induce antiviral response . In this study, the 523 differentially expressed mRNAs are enriched in Toll-like receptors, RIG-I-like receptors, TNF, NOD-like receptors, Jak-STAT, and NF-κB signaling pathways. NF-κB pathway is a crucial pathway in inflammatory process and can be activated by Jak-STAT, TNF, Toll-like receptor, and RIG-I-like receptor pathways, which can be activated by TGEV infection [21,22,23,24,25,26,27]. RIG-I is a dsRNA helicase enzyme encoded by the DDX58 gene and recognizes viral double-stranded (ds) RNA to interfere with viral infection [28, 29]. Recently, it has been confirmed that TGEV infection significantly upregulates RIG-I and activates NF-κB pathway . Toll-like receptors (TLRs) are a class of proteins and play a key role in innate immune response by recognizing viral component . TLR3, a member of TLR family, has a fundamental role in recognizing dsRNA of viruses and activation of IRF3 and NF-κB pathway . In addition, other viruses, such as PSV and Influenza A virus, also give rise to activation of NF-κB pathway via recognition of viral RNA by TLR3 [32, 33]. Consistently, we find TGEV infection upregulates RIG-I and TLR3 and induces the activation of NF-κB pathway in IPEC-J2 cells, indicating TGEV may activate NF-κB pathway via RIG-I and TLR3 signaling. Jak-STAT signaling modulates transmission of signal from cell-membrane receptors to the nucleus and is essential for cytokines and growth factors, which cause critical cellular events such as immune responses, cell death, antiviral response . It is found that Jak-STAT pathway induced p65 phosphorylation and activation of NF-κB signaling during TGEV infection [26, 35]. In this study, TGEV infection activates the Jak-STAT signaling pathway in IPEC-J2 cells, consistent with previous studies. IL-6 is a key mediator in response to acute inflammation caused by viral infection, including TGEV, PEDV, Head strong violence (HSV), H1N1 influenza virus, and hepatitis B virus (HBV) [36,37,38,39]. Here, we have observed that TGEV infection up-regulates cytokines IL-6, TNF-α, IL-8, and CCL5. Extensive evidence support the conclusion [27, 40, 41]. IL-6 is both a proinflammatory cytokine and an anti-inflammatory cytokine. IL-6 induces Jak-STAT signaling pathway and subsequently activates NF-KB signaling by binding to IL6R. Moreover, activation of NF-KB can promote transcription of IL-6 [41, 42].
miRNAs and circular RNAs are crucial non-coding RNAs and implicated in various biological and pathological processes such as viral infection, immune responses, and inflammation. It was reported that cellular miRNA profiles were altered during TGEV infection in PK-15 cells and ST cells [16, 43]. In this study, miRNAs profile and circRNAs profile were changed during TGEV infection in IPEC-J2 cell line. In comparison with the miRNAs profile of ST cells and PK-15 cells infected with TGEV, we observed a distinct miRNAs profile, indicating miRNAs are tissue specific. miR-22 is involved in inhibition of myocardial ischemia-reperfusion injury by an anti-inflammation mechanism . We demonstrated that miR-22 was up-regulated by TGEV infection and had an anti-inflammation effect by targeting IL-6 and restraining DDX58, and CCL5, which is correlated with previous reports, but through different mechanisms. circRNAs are an abundant class of endogenous ncRNAs and function as a regulator via various mechanisms such as miRNA sponge, binding to RNA-binding protein, and regulating transcription and translation of gene [44,45,46]. We found that miR-22 could be directly attached to ssc_circ_9380 and mediate TGEV-induced inflammation response, suggesting ssc_circ_9380 acts as the sponge of miR-22.
In the present study, we performed the next generation sequencing to obtain the profiles of mRNAs, miRNAs, and circRNAs of IPEC-J2 cells during TGEV infection and investigated the effects of differentially expressed ncRNAs on TGEV induced inflammatory response. The data reveal that the differentially expressed mRNAs and ncRNAs are enriched in inflammation-related pathways and ssc_circ_009380 potentiate TGEV-induced activation of NF-κB pathway through attaching to miR-22.
Antibodies, cells, and virus
Histone H3 monoclonal antibody, phospho-I-κB monoclonal antibody, and NF-κB p65 (L8F6) mouse monoclonal antibody were purchased from Cell Signaling Technology (US). β-actin monoclonal antibody was purchased from Santa Cruz (US). Horseradish peroxidase (HRP)-conjugated secondary antibody was purchased from Pierce (US). Dylight594-conjugated secondary antibody was purchased from Genshare Biological (China). IPEC-J2 cell line was kindly gifted by Dr. Zhanyong Wei (Henan Agricultural University, China). Cells were cultured in Dulbecco’s Modified Eagle Medium (DMEM)/F-12/HAM (Thermo Fisher Scientific, US) supplemented with 100 IU of penicillin and 100 mg of streptomycin per ml, at 37 °C in an incubator with 5% CO2. The TGEV Shaanxi strain was separated from TGEV-infected piglets . miR-22 mimics, miRNA mimics control, miR-22 inhibitors, miRNA inhibitors control, siCirc009380, and negative control were synthesized by GenePharma (China) (The sequences were shown in Additional file 1: Table S8).
Strand-specific library construction and sequencing of mRNA and circRNA
IPEC-J2 cells were infected with TGEV at 1 MOI for 24 h (indicated by T1 and T2). Meanwhile, the mock infection (indicated by M1 and M2) was carried out. Total RNA was extracted with Trizol reagent (Invitrogen, Carlsbad, CA, US). Then rRNAs were removed to retain mRNAs and ncRNAs. The mRNAs and ncRNAs were fragmented into short fragments using fragmentation buffer and reversely transcribed into cDNA using M-MLV kit (Invitrogen, US). Second-strand cDNA were synthesized by DNA polymerase I, RNase H, dNTP (dUTP instead of dTTP) and buffer. Next, the cDNA fragments were purified with QiaQuick PCR extraction kit, end repaired, poly(A) added, and ligated to Illumina sequencing adapters. Then UNG (Uracil-N-Glycosylase) was used to digest the second-strand cDNA. The digested products were separated using agarose gel electrophoresis, amplified through PCR, and sequenced using Illumina HiSeqTM 2500 by Gene Denovo Biotechnology Co. (Guangzhou, China).
Alignment with reference genome
To get high quality clean reads, reads containing adapters, low quality reads, and rRNA reads were removed and then mapped to Sus scrofa reference genome (Sus scrofa 10.2) by TopHat2 (version 188.8.131.52) with default options.
The reconstruction of transcripts was carried out using Cufflinks (V2.2.1). The program reference annotation-based transcripts (RABT) was preferred. Cufflinks constructed faux reads according to reference to make up for the influence of low coverage sequencing. During the last step of assembly, all of the reassembled fragments were aligned with reference genes and then similar fragments were removed. Then we used Cuffmerge to merge transcripts from different replicates of a group into a comprehensive set of transcripts, and then merged the transcripts from multiple groups into a finally comprehensive set of transcripts for further downstream differential expression analysis.
Identification and annotations for novel transcripts
To identify the new transcripts, all of the reconstructed transcripts were aligned with reference genome and divided into twelve categories using Cuffcompare (V2.2.1). We used the following parameters to identify reliable novel transcripts: the length of transcript was longer than 200 bp and the exon number was more than 2. The novel transcripts were then searched against KEGG database for annotations (http://www.genome.jp/kegg/).
Quantification of transcript abundance
Transcripts abundance was quantified by RSEM (V1.2.8) and normalized to FPKM (Fragments Per Kilobase of transcript per Million mapped reads). The formula is shown as follow:
C, the number of fragments that are mapped to transcripts; N, the total number of fragments that are mapped to reference genes; L, the number of base pairs of transcript.
Identification and statistics of circRNAs
The Anchors reads were mapped to Sus scrofa reference genome, then subjected to find_circ (V1.0) to identify circRNAs. The identified circRNAs were subjected to statistical analysis of type, chromosome distribution, and length distribution.
Quantification of circRNA abundance
To quantify circRNAs, back-spliced junction reads were scaled to RPKM (Reads per kilobase of transcript per Million mapped reads). The formula is shown below:
C,the number of reads that were mapped to transcripts; N, the total number of reads that were mapped to reference genes; L, the number of base pairs of transcripts.
Library construction and sequencing of miRNAs
After total RNA was extracted using TRIzol, the RNA molecules in a size range of 18~ 30 nt were enriched by polyacrylamide gel electrophoresis (PAGE). Then the 3′ adapters were added. The 36~ 44-nt RNAs were enriched. The 5′ adapters were then ligated to the RNAs as well. The ligation products were reversely transcribed by PCR amplification. The 140~ 160-bp PCR products were enriched to generate a cDNA library and sequenced using Illumina HiSeqTM 2500 by Gene Denovo Biotechnology Co. (Guangzhou, China).
Alignment of miRNAs
To get clean tags, we removed low quality reads. All of the clean tags were aligned with small RNAs in GeneBank database (Release 209.0) and Rfam database (11.0) to identify and remove rRNA, scRNA, snoRNA, snRNA, and tRNA. Then, the rest of clean tags was aligned with Sus scrofa reference genome and searched against miRBase database (Release 21) to identify existed miRNAs. Unmapped miRNAs were aligned with other species. All of the unannotated tags were aligned with Sus scrofa reference genome. According to the genome positions and hairpin structures predicted by software Mireap (v0.2), the novel miRNA candidates were obaineds.
Quantification of miRNA abundance
The miRNA level was calculated and normalized to transcripts per million (TPM). The formula is show below:
Prediction of miRNAs targets
RNAhybrid (v2.1.2) + svm_light (v6.01), Miranda (v3.3a), and TargetScan (v7.0) were used to predict targets of miRNAs. The intersection of the targets was chosen as candidate targets of miRNAs.
Significance analysis of miRNAs, mRNAs, and circRNAs
The edgeR package (http://www.r-project.org/) was used to identify differentially expressed miRNAs, mRNAs, and circRNAs. A fold change ≥2 and ≤ 0.5, plus a P value < 0.05 were set as thresholds for significant differentially expressed miRNAs and circRNAs. mRNAs with a fold change ≥2 and ≤ 0.5, plus a false discovery rate (FDR) < 0.05, were identified as significant differentially mRNAs.
KEGG and interaction analysis of differentially expressed mRNAs, miRNAs, and circRNAs
KEGG database (http://www.genome.jp/kegg/) was used to annotate the pathways. The formula for P-value is shown below:
N,the number of all genes with KEGG annotation; N, the number of target genes (differentially expressed mRNAs, target genes of differentially expressed miRNAs, or source genes of differentially expressed circRNAs); M, the number of all genes in specific pathways. m; the number of target genes in specific pathways.
p-value ≤0.05 is set as a threshold. The interaction networks among miRNAs, circRNAs, and mRNAs were built and visualized using Cytoscape (v3.5.1) (http://www.cytoscape.org/).
Both divergent primers and convergent primers were designed to identify the circular form. Head-to-tail splicing was validated by PCR and sequencing after reverse transcription.
Quantification of miRNAs, circRNAs, and mRNAs using qRT-PCR
The total RNA of IPEC-J2 cells were isolated using TRIzol reagent and reversely transcribed using M-MLV reverse transcriptase (Invitrogen, US) according to the manufacturer′s instructions. For reverse transcription of circRNAs, the total RNA was treated with RNase R (Epicentre, US) to remove linear RNA and then was subjected to reverse transcription. qRT-PCR was performed on iQ5 real-time PCR System (Bio-Rad, US) . The primers are shown in Additional file 1: Table S9. Relative quantification of miRNAs was normalized to U6 using 2-ΔΔCt method. The relative quantification of mRNAs and circRNAs were normalized to GAPDH using the 2-ΔΔCt method.
The sequences of circRNAs and 3′ UTR of mRNAs were respectively amplified by PCR and cloned into vector psiCHECK-2 (Promega, US) (Primers are shown in Additional file 1: Table S9). The binding sites of miR-22 in targets and circRNA were mutated  (Additional file 1: Table S9). IPEC-J2 cells were co-transfected with 100 ng of plasmid and 100 nM of miR-22 mimics (or miRNA mimics control, miR-22 inhibitors, miRNA inhibitors control) using Lipofectamine 3000 (Invitrogen, US). At 48 h post transfection (hpt), the luciferase activities were tested using Dual-Glo Luciferase Assay System (Promega, US) following the manufacturer′s manual.
Western blot analysis
Cells were treated with RIPA lysis buffer containing phenylmethyl sulfonylfluoride (PMSF). Protein concentration was measured using BCA Protein Assay Reagent (Pierce, US). Proteins were separated on a 12% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and subsequently transferred onto a polyvinylidene difluoride (PVDF) membranes (Millipore, US). The PVDF membrane was blocked with 5% non-fat milk for 2 h at room temperature and then subsequently incubated with primary antibodies overnight at 4 °C and HRP-conjugated secondary antibody at room temperature for 1 h. Finally, the membrane was developed with enhanced chemiluminescence (ECL) (Promega, US).
The data are presented as the means ± SEM. Statistical significance were analyzed by unpaired Student′s t-test. p < 0.05 was defined as statistical significance.
Dulbecco’s Modified Eagle Medium (DMEM)/F-12/HAM
Fragments per kilobase of transcript per million mapped reads
Hepatitis B virus
Hour post transfection
Herpes simplex virus
Porcine intestinal epithelial cells-jejunum 2
Kyoto Encyclopedia of Genes and Genomes
Multiplicity of infection
MiRNA response elements
Pathogen recognition receptors
Retinoic acid-inducible gene I
Reads per kilobase of transcript per million mapped reads
Sodium dodecyl sulfate-polyacrylamide gel electrophoresis
Transmissible gastroenteritis virus
Transcripts per million
Zhu L, Mou C, Yang X, Lin J, Yang Q. Mitophagy in TGEV infection counteracts oxidative stress and apoptosis. Oncotarget. 2016;7(19):27122–41.
Zhu YG, Qu JM. Differential characteristics of the early stage of lung inflammation induced by SARS-CoV Nucleocapsid protein related to age in the mouse. Inflamm Res. 2009;58(6):312–20.
Oudit GY, Kassiri Z, Jiang C, Liu PP, Poutanen SM, Penninger JM, Butany J. SARS-coronavirus modulation of myocardial ACE2 expression and inflammation in patients with SARS. Eur J Clin Investig. 2009;39(7):618–25.
Houser KV, Broadbent AJ, Gretebeck L, Vogel L, Lamirande EW, Sutton T, Bock KW, Minai M, Orandle M, Moore IN, et al. Enhanced inflammation in New Zealand white rabbits when MERS-CoV reinfection occurs in the absence of neutralizing antibody. PLoS Pathog. 2017;13(8):e1006565.
Guo X, Hu H, Chen F, Li Z, Ye S, Cheng S, Zhang M, He Q. iTRAQ-based comparative proteomic analysis of Vero cells infected with virulent and CV777 vaccine strain-like strains of porcine epidemic diarrhea virus. J Proteome. 2016;130:65–75.
Mattick JS, Makunin IV: Non-coding RNA. Human molecular genetics 2006, 15 Spec No 1 :R17–R29.
Ratnik K, Viru L, Merits A. Control of the rescue and replication of Semliki Forest virus recombinants by the insertion of miRNA target sequences. PLoS One. 2013;8(9):e75802.
Bianchessi V, Badi I, Bertolotti M, Nigro P, D’Alessandra Y, Capogrossi MC, Zanobini M, Pompilio G, Raucci A, Lauri A. The mitochondrial lncRNA ASncmtRNA-2 is induced in aging and replicative senescence in endothelial cells. J Mol Cell Cardiol. 2015;81:62–70.
Dasgupta N, Peng Y, Tan Z, Ciraolo G, Wang D. Li R: miRNAs in mtDNA-less cell mitochondria. Cell death discovery. 2015;1:15004.
Ji J, Qin Y, Ren J, Lu C, Wang R, Dai X, Zhou R, Huang Z, Xu M, Chen M, et al. Mitochondria-related miR-141-3p contributes to mitochondrial dysfunction in HFD-induced obesity by inhibiting PTEN. Sci Rep. 2015;5:16262.
Zhou R, Wang R, Qin Y, Ji J, Xu M, Wu W, Chen M, Wu D, Song L, Shen H, et al. Mitochondria-related miR-151a-5p reduces cellular ATP production by targeting CYTB in asthenozoospermia. Sci Rep. 2015;5:17743.
Xie N, Liu G. ncRNA-regulated immune response and its role in inflammatory lung diseases. Am j physiol Lung cell molec physiol. 2015;309(10):L1076–87.
Grimaldi A, Zarone MR, Irace C, Zappavigna S, Lombardi A, Kawasaki H, Caraglia M, Misso G. Non-coding RNAs as a new dawn in tumor diagnosis. Semin Cell Dev Biol. 2018;78:37-50.
Fan C, Tang Y, Wang J, Xiong F, Guo C, Wang Y, Zhang S, Gong Z, Wei F, Yang L, et al. Role of long non-coding RNAs in glucose metabolism in cancer. Mol Cancer. 2017;16(1):130.
Han Li C, Chen Y. Small and long non-coding RNAs: novel targets in perspective Cancer therapy. Current genomics. 2015;16(5):319–26.
Song X, Zhao X, Huang Y, Xiang H, Zhang W, Tong D. Transmissible gastroenteritis virus (TGEV) infection alters the expression of cellular microRNA species that affect transcription of TGEV gene 7. Int J Biol Sci. 2015;11(8):913–22.
Zhao X, Song X, Bai X, Fei N, Huang Y, Zhao Z, Du Q, Zhang H, Zhang L, Tong D. miR-27b attenuates apoptosis induced by transmissible gastroenteritis virus (TGEV) infection via targeting runt-related transcription factor 1 (RUNX1). PeerJ. 2016;4:e1635.
Gidlof O, Sathanoori R, Magistri M, Faghihi MA, Wahlestedt C, Olde B, Erlinge D. Extracellular Uridine triphosphate and adenosine triphosphate attenuate endothelial inflammation through miR-22-mediated ICAM-1 inhibition. J Vasc Res. 2015;52(2):71–80.
Yang J, Fan Z, Yang J, Ding J, Yang C, Chen L. microRNA-22 attenuates myocardial ischemia-reperfusion injury via an anti-inflammatory mechanism in rats. Experimental and therapeutic medicine. 2016;12(5):3249–55.
Zhou Y, Wu W, Xie L, Wang D, Ke Q, Hou Z, Wu X, Fang Y, Chen H, Xiao S, et al. Cellular RNA helicase DDX1 is involved in transmissible gastroenteritis virus nsp14-induced interferon-Beta production. Front Immunol. 2017;8:940.
Hartman ZC, Poage GM, den Hollander P, Tsimelzon A, Hill J, Panupinthu N, Zhang Y, Mazumdar A, Hilsenbeck SG, Mills GB, et al. Growth of triple-negative breast cancer cells relies upon coordinate autocrine expression of the proinflammatory cytokines IL-6 and IL-8. Cancer Res. 2013;73(11):3470–80.
Kell AM, Gale M Jr. RIG-I in RNA virus recognition. Virology. 2015;479-480:110–21.
Zhang HX, Liu ZX, Sun YP, Zhu J, Lu SY, Liu XS, Huang QH, Xie YY, Zhu HB, Dang SY, et al. Rig-I regulates NF-kappaB activity through binding to Nf-kappab1 3’-UTR mRNA. Proc Natl Acad Sci U S A. 2013;110(16):6459–64.
Chen S, Zhu J, Chen G, Zuo S, Zhang J, Chen Z, Wang X, Li J, Liu Y, Wang P. 1,25-Dihydroxyvitamin D3 preserves intestinal epithelial barrier function from TNF-alpha induced injury via suppression of NF-kB p65 mediated MLCK-P-MLC signaling pathway. Biochem Biophys Res Commun. 2015;460(3):873–8.
De Simone V, Franze E, Ronchetti G, Colantoni A, Fantini MC, Di Fusco D, Sica GS, Sileri P, MacDonald TT, Pallone F, et al. Th17-type cytokines, IL-6 and TNF-alpha synergistically activate STAT3 and NF-kB to promote colorectal cancer cell growth. Oncogene. 2015;34(27):3493–503.
An K, Fang L, Luo R, Wang D, Xie L, Yang J, Chen H, Xiao S. Quantitative proteomic analysis reveals that transmissible gastroenteritis virus activates the JAK-STAT1 signaling pathway. J Proteome Res. 2014;13(12):5376–90.
Ding Z, An K, Xie L, Wu W, Zhang R, Wang D, Fang Y, Chen H, Xiao S, Fang L. Transmissible gastroenteritis virus infection induces NF-kappaB activation through RLR-mediated signaling. Virology. 2017;507:170–8.
Zhang Q, Wan H, Huang S, Zhang Y, Wang Y, Guo X, He P, Zhou M. Critical role of RIG-I-like receptors in inflammation in chronic obstructive pulmonary disease. Clin Respir J. 2016;10(1):22–31.
Zhu H, Lou F, Yin Q, Gao Y, Sun Y, Bai J, Xu Z, Liu Z, Cai W, Ke F, et al. RIG-I antiviral signaling drives interleukin-23 production and psoriasis-like skin disease. EMBO molecular medicine. 2017;9(5):589–604.
Kawai T, Akira S. Toll-like receptors and their crosstalk with other innate receptors in infection and immunity. Immunity. 2011;34(5):637–50.
Li C, Chen K, Kang H, Yan Y, Liu K, Guo C, Qi J, Yang K, Wang F, Guo L, et al. Double-stranded RNA released from damaged articular chondrocytes promotes cartilage degeneration via toll-like receptor 3-interleukin-33 pathway. Cell Death Dis. 2017;8(11):e3165.
Lan D, Tang C, Yue H, Sun H, Cui L, Hua X, Li J. Microarray analysis of differentially expressed transcripts in porcine intestinal epithelial cells (IPEC-J2) infected with porcine sapelovirus as a model to study innate immune responses to enteric viruses. Arch Virol. 2013;158(7):1467–75.
Coates BM, Staricha KL, Wiese KM, Ridge KM. Influenza a virus infection, innate immunity, and childhood. JAMA Pediatr. 2015;169(10):956–63.
Ahmad SF, Ansari MA, Zoheir KM, Bakheet SA, Korashy HM, Nadeem A, Ashour AE, Attia SM. Regulation of TNF-alpha and NF-kappaB activation through the JAK/STAT signaling pathway downstream of histamine 4 receptor in a rat model of LPS-induced joint inflammation. Immunobiology. 2015;220(7):889–98.
Yang CW, Lee YZ, Hsu HY, Shih C, Chao YS, Chang HY, Lee SJ. Targeting Coronaviral replication and cellular JAK2 mediated dominant NF-kappaB activation for comprehensive and ultimate inhibition of Coronaviral activity. Sci Rep. 2017;7(1):4105.
Aoki R, Kawamura T, Goshima F, Ogawa Y, Nakae S, Nakao A, Moriishi K, Nishiyama Y, Shimada S. Mast cells play a key role in host defense against herpes simplex virus infection through TNF-alpha and IL-6 production. J investigative dermatol. 2013;133(9):2170–9.
Chiaretti A, Pulitano S, Barone G, Ferrara P, Romano V, Capozzi D, Riccardi R. IL-1 beta and IL-6 upregulation in children with H1N1 influenza virus infection. Mediat Inflamm. 2013;2013:495848.
Quetier I, Brezillon N, Duriez M, Massinet H, Giang E, Ahodantin J, Lamant C, Brunelle MN, Soussan P, Kremsdorf D. Hepatitis B virus HBx protein impairs liver regeneration through enhanced expression of IL-6 in transgenic mice. J Hepatol. 2013;59(2):285–91.
Huan CC, Wang HX, Sheng XX, Wang R, Wang X, Liao Y, Liu QF, Tong GZ, Ding C, Fan HJ, et al. Porcine epidemic diarrhea virus nucleoprotein contributes to HMGB1 transcription and release by interacting with C/EBP-beta. Oncotarget. 2016;7(46):75064–80.
Ma R, Zhang Y, Liu H, Ning P. Proteome profile of swine testicular cells infected with porcine transmissible gastroenteritis coronavirus. PLoS One. 2014;9(10):e110647.
Fang XX, Jiang XL, Han XH, Peng YP, Qiu YH. Neuroprotection of interleukin-6 against NMDA-induced neurotoxicity is mediated by JAK/STAT3, MAPK/ERK, and PI3K/AKT signaling pathways. Cell Mol Neurobiol. 2013;33(2):241–51.
Hartman ZC, Poage GM, den Hollander P, Tsimelzon A, Hill J, Panupinthu N, Zhang Y, Mazumdar A, Hilsenbeck SG, Mills GB, et al. Growth of triple-negative breast cancer cells relies upon coordinate autocrine expression of the pro-inflammatory cytokines IL-6 and IL-8. Cancer research. 2013;73(11). https://doi.org/10.1158/0008-5472.CAN-1112-4524-T.
Liu X, Zhu L, Liao S, Xu Z, Zhou Y. The porcine microRNA transcriptome response to transmissible gastroenteritis virus infection. PLoS One. 2015;10(3):e0120377.
Fang Y, Wang X, Li W, Han J, Jin J, Su F, Zhang J, Huang W, Xiao F, Pan Q, et al. Screening of circular RNAs and validation of circANKRD36 associated with inflammation in patients with type 2 diabetes mellitus. Int J Mol Med. 2018;42(4):1865–74.
Zheng Q, Bao C, Guo W, Li S, Chen J, Chen B, Luo Y, Lyu D, Li Y, Shi G, et al. Circular RNA profiling reveals an abundant circHIPK3 that regulates cell growth by sponging multiple miRNAs. Nat Commun. 2016;7:11215.
Zhong Y, Du Y, Yang X, Mo Y, Fan C, Xiong F, Ren D, Ye X, Li C, Wang Y, et al. Circular RNAs function as ceRNAs to regulate and control human cancer progression. Mol Cancer. 2018;17(1):79.
Ding L, Chen GD, Xu XG, Tong DW. Isolation and identification of porcine transmissible gastroenteritis virus Shaanxi strain and sequence analysis of its. N gene Chin J Vet Med. 2011;47(10):09–11.
We thank Zhanyong Wei for gifting IPEC-J2 cell line.
This work was supported by National Natural Science Foundation of China (Grant No. 31472167), China Postdoctoral Science Foundation (Grant No. 2015 M570860), Key Research and Development Project in Shaanxi Province (Grant No. 2018ZDXM-NY-064), Central Project of Major Agricultural Technology Promotion Funds (Grant No. K3360217060), and Science and Technology Planning Project of Yangling demonstration zone (Grant No. 2018NY-35).
Availability of data and materials
All the sequencing data supporting the results in the current study have been deposited in NCBI with the accession number SRP128150 (https://www.ncbi.nlm.nih.gov/sra/SRP128150).
Consent to publish
Ethics approval and consent to participate
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. The detailed information of differentially expressed mRNAs. Table S2. The detailed information of differentially expressed miRNAs. Table S3. The detailed information of differentially expressed circRNAs. Table S4. The target genes of differentially expressed miRNAs. Table S5. Source genes of differentially expressed circRNAs. Table S6. The interactions between miRNAs and mRNAs. Table S7. The circRNA-miRNA-mRNA regulation network. Table S8. The sequences of miR-22 mimics and inhibitors. Table S9. The primers designed of qRT-PCR and PCR. (ZIP 282 kb)
About this article
Cite this article
Ma, X., Zhao, X., Zhang, Z. et al. Differentially expressed non-coding RNAs induced by transmissible gastroenteritis virus potentially regulate inflammation and NF-κB pathway in porcine intestinal epithelial cell line. BMC Genomics 19, 747 (2018). https://doi.org/10.1186/s12864-018-5128-5
- Circular RNA