Molecular Characterization and Biological Function of a Novel LncRNA CRNG in Swine

Abstract

Our previous study has showed that a novel gene is differentially expressed in the liver of cyadox-fed piglets, but its sequence and function are unknown. Here, rapid amplification of cDNA ends (RACE) and bioinformatics analysis showed that the novel gene is 953 bp without protein-coding ability and locates in chromosome 11. Hence, we identified the novel gene as long non-coding RNA (lncRNA) and named it cyadox-related novel gene (CRNG). Fluorescence in situ hybridization (FISH) showed that CRNG mainly distributes in cytoplasm. Moreover, microarray assay in combination with CRNG interference and overexpression showed that the differential genes such as ANPEP, KITLG, STAT5A, FOXP3, miR-451, IL-2, IL-10, IL-6, and TNF-α are mainly involved in viral and pathogens infection and the immune-inflammatory responses in PK-15 cells. This work reveals that CRNG might play a role in preventing the host from being infected by pathogens and viruses and exerting immune regulatory effects in the cytoplasm, which may be involved in prophylaxis of cyadox in piglets.

Introduction

Recent advances in high throughput sequencing technology have led to markedly expanding our knowledge of transcriptome, i.e., LncRNAs, which are longer than 200 nucleotides and unable to be translated into proteins (Li et al., 2014; Wu et al., 2014). The role of lncRNAs may far exceed people’s expectations, and increasing amount of evidence suggests that lncRNAs may play critical roles in control of a wide array of cellular functions, for example viral infection (Lemler et al., 2017), immune cell homeostasis and function (Chen et al., 2017; Mowel et al., 2018), inflammation (Chew et al., 2018), forcing us to radically attach importance to lncRNAs. In particular, attention is now shifting toward a type of RNA called LncRNAs, which is one of the most poorly understood, yet most common (Dinger et al., 2008; Kwok and Tay, 2017; Li et al., 2018).

Cyadox is a novel derivative of quinoxaline-1,4-dioxides with antibacterial and growth promotion effects (Cui et al., 2018; Guo et al., 2018) and with the potential to serve as a feed additive (Liu et al., 2018). A previous study in our lab has shown that eight differential genes were found through mRNA differential display technology in the liver of cyadox-fed piglets, including insulin-like growth factor-1, epidermal growth factor, poly ADP-ribose polymerase, the defender against apoptotic death 1, complement component 3, transketolase, sus scrofa zinc finger CCHC domain containing 3 and a novel gene (Yu et al., 2018). Additionally, the sequence of the novel gene is matched with a predicted sequence, sus scrofa uncharacterized LOC100626416, in nation center for biotechnology information, which was constantly updated in 2015, 2017, and 2018, and has not been experimentally demonstrated. Moreover, the novel gene was demonstrated without highly homologous protein sequences in protein data bank (Yu et al., 2018). Hence, we tentatively named the gene as CRNG. In the previous study, it was found that CRNG is related to the NF-kB, P38, TGF-β, JNK, PI3K, and JAK-STAT signaling pathway in primary cultured pig hepatocytes exposure to cyadox (Guo et al., 2018), suggesting that CRNG might be very important in cyadox action. Consequently, it is necessary to illustrate the full length, structural characteristics and biological function of CRNG, and further explain the role of the CRNG in the cyadox-mediated functional effects.

Considerable research efforts have been devoted to study the basic characterizations of the unknown gene via using the technology of RACE (Dieffenbach et al., 2003) to clone the full-length of novel gene, and bioinformatics, CPC (Kong et al., 2007) and ORF finder (Garcia et al., 2015), to analyze the ability of ending protein of novel genes. Additionally, RNA FISH was used to investigate the subcellular localization of lncRNAs to further elucidate the mechanisms and functions of lncRNAs (Huang et al., 2017; Das et al., 2018). Moreover, the combination of microarray and RT-qPCR can better analyze the properties and functions of genes (Shi et al., 2018).

This study aims at illustrating the characterization, function of CRNG in swine. The study on cyadox-related gene CRNG will help to provide a new sight on the pharmacological mechanism of cyadox. Our study showed that CRNG is a non-coding RNA mainly distributed in liver, followed by the jejunum and duodenum, and again the kidney of swine, and the cytoplasm of PK-15 cells. Microarray and RT-qPCR reveal important biological functions of CRNG, such as regulation of inflammation, pathogen infection and antiviral immunity, which provides a new view to better explain the development and application of cyadox and relevant immune mechanisms of CRNG in swine.

Results

Molecular Characteristics of LncRNA CRNG

To explore the biological functions of CRNG, 953 bp of the full-length cDNA sequence of CRNG was obtained by 5′ and 3′ RACE (Figure 1A). Using BLAST searches for porcine HTGS database we obtained a partially match (GenBank NC-010453.5, 76514352–76516982) with the obtained porcine CRNG cDNA, which revealed that the genomic sequence of swine CRNG located in chromosome 11. Then CRNG sequence was searched in the whole genome of pigs, and the cDNA of which consists of exon 1 (76513882–76514572), exon 2 (76516721–76516981) and one intron (76514573–76516720). However, the exon 1 obtained by 5′ RACE is not full-match to the genomic sequence, which was verified by sequencing. It means that the CRNG sequence is one base “C” more than the genome sequence at 76514142 and base “C” replaces base “T” at 76514197 (Supplementary Data Sheet S1), which was further determined to be 953 bp by DNA amplification and sequencing (Figure 1B).

FIGURE 1

The ORF of CRNG was predicted by bioinformatics tools in NCBI¹ (Table 2). The prediction results showed that the ORF has no coding capacity in porcine genome. Therefore, we also assessed the protein-coding potential of CRNG and its reverse complementary in a Support Vector Machine-based classifier: CPC², it suggested that the gene CRNG was a non-coding RNA (Figure 2).

FIGURE 2

Transcriptional Levels of CRNG in Different Tissues of Swine

The transcriptional levels of CRNG were evaluated relative to the endogenous β-actin mRNA levels in heart, liver, spleen, lung, kidney, duodenum, jejunum, ileum, longissimus dorsi muscle, thymus, hypothalamus, pituitary, bone marrow, adrenal gland, cecum, colon, rectum of swine by RT-qPCR. The result demonstrated that mRNA level of CRNG was the highest in liver, high in duodenum and jejunum, low in kidney, adrenal gland, hypothalamus and pituitary, very low in other tissues (Figure 3). In the pretest, the growth rate of the transfected IPEC-J2 cells slowed down, and there was a tendency to shed death, and the effect after transfection could not be detected. Therefore, PK-15 cells were finally selected as the cell line for the transfection experiment.

FIGURE 3

Subcellular Localization of CRNG

To gain insights into the precise mechanism of CRNG, we first examine its subcellular localization because the function of lncRNA depends on its subcellular distribution (Kopp and Mendell, 2018). FISH was used to observe the distribution of CRNG in PK-15 cells. As observed in Figure 4, CRNG predominantly is distributed in the cytoplasmic region of PK-15 cells. This finding provided the evidence that CRNG might act as endogenous sponge RNA to influence the expression of its target mRNA in cytoplasm.

FIGURE 4

Global Differential Gene Expression in PK-15 Cells Interfering CRNG

To gain further insight into how the manipulation of CRNG might modulate the function and biological process in PK-15 cells, we next determined the effects of siRNA-CRNG and siRNA-NC in PK-15 cells for 24 h by performing RNA-sequencing analysis. The statistical analysis showed that absence of the CRNG significantly affect the expression of 78 genes, including 32 upregulated genes and 46 downregulated genes. These gene sets were related to immune response pathways including interferon signaling, the cellular response to cytokine production, and viral process (Figure 5 and Supplementary Data Sheet S2). Additionally, the differential genes were individually annotated and elucidated the main biological process via GO-based enrichment analysis, including the cellular component, biological process and molecular function. Intriguingly, the analysis data of GO biological process confirmed that several biological pathways related to viral production, cytokine production, and inflammation were affected (Figure 6). These finding suggest that CRNG is an important part of immune and inflammation process.

FIGURE 5

FIGURE 6

Validation of Differentially Expressed Genes in PK-15 Cells Overexpressing and Interfering CRNG

To verify the reliability of the microarray assay results, several genes related to inflammation, infection and antiviral immunity, such as ANPEP, KITLG, STAT5A, FOXP3, and miR-451 were selected. Data indicated that CRNG significantly increase KITLG, FOXP3 and miR-451, and decrease ANPEP and STAT5A mRNA expression (Figure 7), suggesting that CRNG may be necessary for the modulation of antiviral immunity, inflammation and pathogen infection. Additionally, according to the GO analysis, several genes associated with inflammation, such as IL-2, IL-10, IL-6, TNF-α, were examined using RT-qPCR in PK-15 cells via overexpressing and interfering CRNG, respectively. The transcriptional levels of IL-2, IL-10, and IL-6 were decreased and increased in boosted-CRNG PK-15 cells, respectively, while that of TNF-α was unchanged (Figure 7). Intriguingly, in opposition to boosted-CRNG, siRNA-CRNG boosted IL-2 mRNA level and decreased IL-6 and TNF-α expression, while IL-10 was induced with no significance (Figure 7). These data indicate that CRNG directly modulate inflammatory response. Taken together, these results suggested that lncRNA CRNG might be a vital regulator of inflammation, pathogen infection and antiviral immunity.

FIGURE 7

Materials and Methods

Cell Line and Culture Conditions

The cryopreserved porcine kidney cell line PK-15 was purchased from American Type Culture Collection (ATCC, Rockville, MD, United States). The cryopreserved PK-15 was quickly thawed at 37°C in a water bath and centrifuged at low speed (50 g) for 5 min to remove DMSO, then suspended and cultured in Dulbecco’s modified Eagle’s medium/high glucose (HyClone, United States) containing 10% fetal bovine serum (Gibco, United States), 2 mM L-glutamine, 100 μg/mL streptomycin and 100 UI/mL penicillin. Cells were maintained in T-25 flasks at 37°C in a humidified 5% CO₂ atmosphere with daily exchanges of fresh culture medium. When cells approached 100% confluence, cells were rinsed with PBS twice, and incubated with fresh 0.25% trypsin for 3 to 5 min, then sub-cultured at a ratio of 1:3 or 1:4 into new flasks. Cells were passaged three times prior to use in experiments.

RT-qPCR

The RNA of PK-15 cells were lysed with Trizol^® Reagent (Invitrogen, United States) according to the instructions. The total RNA contaminating DNA was degraded by treating each sample with RNase-Free DNase Set (Sangon Biotech, China). After purification, total RNA was quantified by optical density (Q3000, Quawell, United States) and the quality was evaluated by gel electrophoresis. The first cDNA was synthesized by HiScript II 1st strand cDNA synthesis kit (Vazyme, China). All primers of target genes were commercially synthesized by Tianyi huiyuan biotechnology company, China (Table 1). The conditions of RT-qPCR were conducted by ChamQ SYBR Color qPCR Master Mix (Vazyme, China) according to the manufacturer’s instructions. The mRNA and miRNA levels were normalized against the amount of the housekeeping gene transcript β-actin and U6, respectively. And the relative expression was calculated by the 2^-ΔΔCt method (Wang X. et al., 2012).

Expression Profile of CRNG Gene

Seventeen porcine tissues, i.e., heart, liver, spleen, lung, kidney, duodenum, jejunum, ileum, longissimus dorsi muscle, thymus, hypothalamus, pituitary, bone marrow, adrenal gland, cecum, colon, and rectum, were collected from three Landrace × Large White crossbred barrows in age from 5 to 6 weeks from the pig farm of Huazhong Agricultural University. The current study was approved by the Ethical Committee of the Faculty of Veterinary Medicine (Huazhong Agricultural University). All the tissue RNAs were extracted completely according to the EZ-10 total RNA mini-preps kit instructions (Sangon Biotech, China). β-actin gene was used as an internal control. The expression profile and internal control primers were designed according to the CRNG sequence (Table 1). Analysis of relative gene expression was conducted using RT-qPCR and the 2^-ΔΔCt method (Wang X. et al., 2012).

3′ RACE and 5′ RACE PCR

The 3′ RACE and 5′ RACE experiment was performed in full accordance with the classical RACE method of PCR Primer: a laboratory manual (Dieffenbach et al., 2003). For 3′ RACE, the first cDNA was synthesized by reverse transcriptase M-MLV (RNase H-; Takara, Japan) with Q_T, then the ribonuclease H (Beyotime, China) was utilized to remove the RNA from cDNA. Simultaneously the cDNA was purified by SanPrep column PCR product purification kit (Sangon Biotech, China). The outer PCR was carried out with gene specific primer CRNG1 and Q₀ via the phanta max super-fidelity DNA polymerase (Vazyme, China). And the inner PCR was carried out with gene specific primer CRNG2 and Q₁ via taq DNA Polymerase (Mg²⁺ Plus Buffer; Vazyme, China). For 5′ RACE, the first cDNA was synthesized by HiScript II 1st strand cDNA synthesis kit (Vazyme, China) with gene specific primer GSP1, then the ribonuclease H (Beyotime, China) was utilized to remove the RNA from cDNA. And the poly “A” tail was added at the end of cDNA by dATP and terminal deoxynucleotidyl transferase (Takara, Japan). Then the cDNA was purified by SanPrep column PCR product purification kit (Sangon Biotech, China). The outer PCR was carried out with gene specific primer GSP2, Q₀ and Q_T via the phanta max super-fidelity DNA polymerase (Vazyme, China), and the inner PCR was carried out with gene specific primer GSP3 and Q₁ via taq DNA Polymerase (Mg²⁺ Plus Buffer; Vazyme, China). All the PCR products were analyzed by 1% agarose gel and sequenced by TA cloning (Tianyi huiyuan, China). All the primers were summarized in Table 1.

LncRNA – CRNG Identification

The full-length cDNA of CRNG was obtained via RACE. Simultaneously, for identifying the characteristic of CRNG, genome location and coding potential were analyzed via NCBI Genome and ORF Finder (Zhang et al., 2012), CPC (Kong et al., 2007), respectively. Additionally, the ORF finder can search novel DNA sequence for potential protein encoding fragments, then verify predicted protein via SMART BLAST or regular BLASTP. No homologous protein sequence of the result of ORF Finder and the score of CPC < 0 were used to assess the coding potential of transcripts. Then combine the results of the above coding ability analysis software to determine the coding ability of lncRNA.

RNA FISH

RNA FISH was performed based on the protocol in the Regulatory Non-Coding RNAs (Carmichael, 2015) with some adjustments. A FITC-Oligo nucleotide probe was designed and synthesized by Shanghai Bogoo Bio-Technique Co., Ltd., for the target gene CRNG: 5′-GGTGACGGCACGTTTGACCTCGTTGAGATGGTG-3′. The specificity of the probe sequence was detected in the NCBI database. Cultured cells were seeded in the 6-wells of chamber slides and exactly fixed by paraformaldehyde. The slides were pretreated for hybridization by a 0.01 % proteinase K digestion (20 min, 37°C) in 0.01 M HCl, followed by a short wash in 0.1 M Glycine Irrigation. Then, the slides were fixed in 4% formaldehyde for 10 min at room temperature. The slides were washed in PBS three times. The slides were washed with acetic anhydride pH = 8.0 (acetylation, reduce the background) at room temperature, 5 min, 2 times. The slides were washed in PBS two times, and then washed in 5 × SSC (pH = 7.5), 1 min, 2 times. The slides were covered with the hybridization solution and incubated in a humidified chamber for 1 h at 65°C. The hybridization mixture contained 500 ng/mL of FITC-probe cover the slides in hybridization instrument at 62°C, 72 h. The slides were washed in 2 × SSC (pH = 7.5), at room temperature, 1 min, 1 time and washed in PBS two times. The slides were stained with DAPI, mounted with anti-quenching agent, and then observed under a fluorescent microscope (IX-71, Olympus, Japan). RNA specificity was confirmed by destruction of signals when samples were pretreated with RNase.

Overexpression and Interference of LncRNA-CRNG

The over-expression plasmid pcDNA3.0-EGFP-CRNG was constructed by Tianyi huiyuan biotechnology company, China. siRNA2 for lncRNA-CRNG (sense 5′-GGACUCUCUUUGACCCUUUTT-3′; antisense 5′-AAAGGGUCAAAGAGAGUCCTT-3′) and a negative control siNC (sense 5′-UUCUCCGAACGGUCACGUTT-3′; antisense 5′-ACGUGACACGUUCGGAGAATT-3′) were obtained from GenePharma, Shanghai, China. The PK-15 cell transfection of the over-expression plasmid was conducted by ExFect^®2000 Transfection Reagent (Vazyme, China) according to the manufacturer’s instructions. The PK-15 cell transfection of siRNA was conducted by Lipofectamine^TM 2000 Transfection Reagent (Invitrogen, United States) according to the manufacturer’s instructions. The efficiency of overexpression or interference was determined by RT-qPCR. Moreover, the siRNA sequence of CRNG and the silencing or overexpression efficiency after transfecting the cells with siRNA or CRNG-expression vector were shown in Supplementary Data Sheet S3.

Microarray Analysis

The genechip was performed by Oebiotech Company (Shanghai, China). Total RNA of siRNA was quantified by the NanoDrop ND-2000 (Thermo Scientific) and the RNA integrity was assessed using Agilent Bioanalyzer 2100 (Agilent Technologies). The sample labeling, microarray hybridization and washing were performed based on the manufacturer’s standard protocols. Briefly, total RNAs were transcribed to double strand cDNAs and then synthesized cRNAs. Next, 2nd cycle cDNAs were synthesized from cRNAs. Followed fragmentation and biotin labeling, the 2nd cycle cDNAs were hybridized onto the microarray. After washing and staining, the arrays were scanned by the Affymetrix Scanner 3000 (Affymetrix). Affymetrix GeneChip Command Console (version 4.0, Affymetrix) software was used to extract raw data. Next, Expression Console (version1.3.1, Affymetrix) software offered RMA normalization for gene. Then the gene expression analysis was carried out.

Statistical Analysis

All data were presented as mean ± SD, and statistical analyses were analyzed by PASW Statistics 18 software. Only p < 0.05 was considered significant. “^∗” indicates p< 0.05; “^∗∗” indicates p< 0.01. GeneSpring software (version 13.1; Agilent Technologies) was employed to finish the basic analysis. Differentially expressed genes were then identified through fold change as well as P-value calculated with t-test. The threshold set for up- and down-regulated genes was a fold change > = 1.5 and a P-value < = 0.05. Afterward, GO analysis was applied to determine the roles of these differentially expressed mRNAs played in these GO terms. Finally, Hierarchical Clustering was performed to display the distinguishable genes’ expression pattern among samples.

Discussion

Long non-coding RNAs, non-coding transcripts of more than 200 nucleotides, are becoming key regulators for the expression and biological processes of some genes (Esteller, 2011). Based on previous study, we firstly clone the full-length sequence of CRNG, which consists of 953 bp without protein-coding ability. Moreover, the distribution of tissue expression profiles suggested that CRNG may play an important role in regulating liver and small intestine development and maintaining their functions. Our results indicate that CRNG is a long-non-coding RNA and might involve in various biological process in the development of pigs. Additionally, RNA FISH revealed that CRNG was predominantly located in the cytoplasm, indicating CRNG might serve as a competing endogenous RNA to sponge miRNAs and restore mRNA translation (Wang et al., 2014), or form complexes with diverse structural and regulatory functions to control mRNA turnover, translation, protein stability, sponging of cytosolic factors, and modulation of signaling pathways (Noh et al., 2018).

Accumulating studies have demonstrated the importance of lncRNAs in the regulation of immune and inflammatory responses (Atianand and Fitzgerald, 2014; Fitzgerald and Caffrey, 2014; Elling et al., 2016). One of the functional studies revealed that lincRNA-Cox2 broadly regulated the expression of a large number of immune genes and inflammatory response genes including pro-inflammatory cytokines and other inflammatory mediators (Carpenter et al., 2013). Previous study illustrated that after exposure to cyadox at the final concentration of 2 μM for 0.5, 1, 2, 4, and 8 h, CRNG mRNA expression would significantly increase in porcine primary hepatocytes compared with that of the control at 1 h. Moreover, the signaling pathways also showed that inflammation-related signaling pathways such as JNK and NF-κB, were involved in the regulation of CRNG (Guo et al., 2018). CRNG mRNA could also be activated by NF-κB, JAK-STAT and JNK pathways in cyadox-treated PK-15 cells (data not published). Present study showed that lncRNA CRNG predominantly associated with viral reproduction, response to virus, pathogen infection and the immune-inflammatory responses through regulating the genes of inflammation, pathogen infection and antiviral immunity, such as ANPEP, KITLG, STAT5A, FOXP3, miR-451 and immune and inflammatory factors, such as IL-2, IL-10, IL-6, TNF-α. Our study revealed that the CRNG has strongly relationship with immune, inflammation and pathogen infection.

ANPEP (also known as CD13) is a membrane metalloprotease consisting of 150 kDa, and an ectoenzyme of zinc-dependent aminopeptidases (Reinhold et al., 2007; Nefla et al., 2015). Recent studies have illustrated that ANPEP could involve in many biological process, such as the regulation of antigen-presenting (Thomas et al., 1994), immunoregulation (Biton et al., 2006), intestinal cholesterol absorption (Kramer et al., 2005). The compound of human CD13 and antibody could inhibit infection and block binding of HCMV virions to susceptible cells (Soderberg et al., 1993). In present study, we found that lncRNA CRNG could negatively regulate the expression of ANPEP, in addition, GO biological process had shown that lncRNA CRNG may involve in virus-host interaction. APN/CD13 could also serve as a modulator to T cell and target of tissue-specific autoimmunity in the CNS (Biton et al., 2006). CD13 served as a negative regulator for activation of mast cells in vitro and in vivo (Metz and Maurer, 2007; Ghosh et al., 2012). The lack of CD13 would activate the inflammation pathways, such as IL-6 and TNF-α, in a Fc𝜀RI-dependent manner (Zotz et al., 2016). From our data, it is clear that the absence of CRNG increased CD13 expression and inhibited the expression of pro-inflammatory factor, IL-6 and TNF-α, suggesting that CRNG may depend on CD13 signaling pathway involved in the regulation of inflammatory cytokines. These results indicated that lncRNA CRNG could inhibit the expression of CD13, thus activating inflammatory immune pathway and preventing the virus infection.

Kit ligand also called mast cell growth factor SCF or c-kit ligand, which is the main survival and developmental factor for mast cells (Ceponis et al., 1998). The combination of hyper-IL-6, SCF and GM-CSF could promote the differentiation of dendritic cells, thus stimulating the resting T cells against the processed antigen (Bernhard et al., 2000). Furthermore, IL-6 and SCF could partly promote the maturation of human cultured mast cells, and other factors may be involved in this process (Matsushima et al., 2000). In this paper, IL-6 and SCF had the same up-regulated trends, when lncRNA CRNG was overexpressed. It demonstrated that lncRNA CRNG may oppose antigen through the network of IL-6 and SCF.

The cellular responses to cytokines, such as SCF and interferons, depend on prior activation of the JAK/STAT signaling pathway. STAT5A showed a strong connection with immunity and virus production in primary CD T cells (Selliah et al., 2006). Notably, SCF not only participates in the regulation of inflammatory factors but also regulates the immunity-regulation through the activity of STAT5, which was reduced in G-CSF or SCF stimulated PNH clone cells (Ding et al., 2012).Additionally, IL-2 also could activate STAT5 and then inhibit the binding of STAT3 to IL17 locus, thus determining the extent of T_H17 cell generation (Yang et al., 2011). Splenocytes from STAT5A-knockout mice showed the partial impairment in IL-2-induced proliferation and defected proliferation in T-lymphocyte and NK-cell (Nakajima et al., 1997; Moriggl et al., 1999). Present study revealed that boosted lncRNA CRNG could decrease the levels of IL-2 and STAT5A, and increase the level of KITLG and IL-6. A reasonable inference can be drawn that lncRNA CRNG could inhibit the activity of STAT5A, and thereby promote inflammation release and reduce virus production, thus revealing the pharmacological activity of cyadox.

FOXP3 is a transcriptional regulators, a member of the forkhead/winged-helix family (Triulzi et al., 2013), implicated in the regulation of the development and inhibitory function of T-cells (Treg), serving as a mediators for self-tolerance, immune homeostasis and various inflammatory responses (Bin Dhuban et al., 2014). The mutation of FOXP3 would lead to the loss of immune homeostasis in mice and humans (Gambineri et al., 2003; Kasprowicz et al., 2003). Moreover, ectopic expression of FOXP3 repressed the production of IL-2 in conventional T cells, through interacting physically with AML1 (Ono et al., 2007), thus controlling the physiological and pathological immune responses mediated by T cell. Additionally, Foxp3+ cells might have a role in the pathogenesis of active periodontal lesions through repressing the expression of TGF-β1 and IL-10 (Scott et al., 2009). They are in agreement with our results performed in PK-15 cells demonstrating that lncRNA CRNG could increase IL-6,TNF-α, FOXP3 mRNA expression, and inhibit IL-10 mRNA expression, thus suppressing viral reproduction and preventing immune hyperactivity.

Evidence that host and pathogens interact through miRNA pathways was documented in mammalian infectious diseases (Roberts et al., 2011; Liang et al., 2013). Wang et al. illustrated that miR-451 could markedly reduce the var gene of parasite virulence factor P. falciparum erythrocyte membrane protein-1 (Wang et al., 2017). miR-451 also was associated with autoimmune disease, which was increased in serum levels (Wang H. et al., 2012; Yamada et al., 2014). Notably, increased levels of IL-6 and IFN-β can positively regulate miR-451 expression, and miR-451 expression ultimately leads to a decrease in IL-6 expression, which in turn buffers IL-6 expression (Rosenberger et al., 2012). Our study also showed that lncRNA CRNG could enhance miR-451 expression, thus protecting host from pathogens infection. Meanwhile, lncRNA CRNG can also increase the expression level of IL-6, but whether CRNG is involved in the regulation of IL-6 through miR-451 remains to be further studied.

Collectively, the molecular characterization and biological function of the novel gene lncRNA CRNG was initially analyzed for the first time. A reasonable conjecture is concluded that lncRNA CRNG could participate in several viral reproductions, responding to virus, infection and immune-inflammation through regulating the genes expression of inflammation and antiviral immunity, including ANPEP, KITLG, STAT5A, FOXP3, miR-451 and immune and inflammatory factors, IL-2, IL-10, IL-6, TNF-α (Figure 8). However, the mechanisms of lncRNA CRNG coupled with those related genes should be further investigated, in order to provide the alternative strategy for preventing the virus replication and further reveal the molecular mechanism of cyadox.

FIGURE 8

References

• 1

  AtianandM. K.FitzgeraldK. A. (2014). Long non-coding RNAs and control of gene expression in the immune system.Trends Mol. Med.20623–631. 10.1016/j.molmed.2014.09.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2

  BernhardH.LohmannM.BattenW. Y.MetzgerJ.LohrH. F.PeschelC.et al (2000). The gp130-stimulating designer cytokine hyper-IL-6 promotes the expansion of human hematopoietic progenitor cells capable to differentiate into functional dendritic cells.Exp. Hematol.28365–372.

  • Pubmed Abstract
  • Google Scholar

• 3

  Bin DhubanK.KorneteM.MasonE. S.PiccirilloC. A. (2014). Functional dynamics of Foxp3(+) regulatory T cells in mice and humans.Immunol. Rev.259140–158. 10.1111/imr.12168

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4

  BitonA.BankU.TagerM.AnsorgeS.ReinholdD.LendeckelU.et al (2006). Dipeptidyl peptidase IV (DP IV, CD26) and aminopeptidase N (APN, CD13) as regulators of T cell function and targets of immunotherapy in CNS inflammation.Adv. Exp. Med. Biol.575177–186. 10.1007/0-387-32824-6_19

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5

  CarmichaelG. G. (2015). Regulatory Non-Coding RNAs. Methods in Molecular Biology.New York, NY: Humana Press.

  • Google Scholar

• 6

  CarpenterS.AielloD.AtianandM. K.RicciE. P.GandhiP.HallL. L.et al (2013). A long noncoding RNA mediates both activation and repression of immune response genes.Science341789–792. 10.1126/science.1240925

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7

  CeponisA.KonttinenY. T.TakagiM.XuJ. W.SorsaT.Matucci-CerinicM.et al (1998). Expression of stem cell factor (SCF) and SCF receptor (c-kit) in synovial membrane in arthritis: correlation with synovial mast cell hyperplasia and inflammation.J. Rheumatol.252304–2314.

  • Pubmed Abstract
  • Google Scholar

• 8

  ChenY. G.SatpathyA. T.ChangH. Y. (2017). Gene regulation in the immune system by long noncoding RNAs.Nat. Immunol.18962–972. 10.1038/ni.3771

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9

  ChewC. L.ConosS. A.UnalB.TergaonkarV. (2018). Noncoding RNAs: master Regulators Of Inflammatory Signaling.Trends Mol. Med.2466–84. 10.1016/j.molmed.2017.11.003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10

  CuiL.XingD.HuangD.LiD.LuQ.WangX.et al (2018). Signaling pathways involved in the expression of SZNF and the target genes binding with SZNF related to cyadox.Biomed. Pharmacother.1081879–1893. 10.1016/j.biopha.2018.09.141

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11

  DasS.ZhangE.SenapatiP.AmaramV.ReddyM. A.StapletonK.et al (2018). A novel angiotensin II-induced long noncoding rna giver regulates oxidative stress, inflammation, and proliferation in vascular smooth muscle cells.Circ. Res.1231298–1312. 10.1161/circresaha.118.313207

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12

  DieffenbachC. W.DvekslerG. S.DieffenbachC. W.DvekslerG. S. (2003). PCR Primer: A Laboratory Manual.New York, NY: Cold Spring Harbor Laboratory Press.

  • Google Scholar

• 13

  DingS. X.FuR.WangH. L.ZhangT.RuanE. B.QuW.et al (2012). Expression of phosphorylated STAT5 in bone marrow hematopoietic stem cells of patients with paroxysmal nocturnal hemoglobinuria before and after in vitro G-CSF or SCF stimulation.Zhonghua Yi Xue Za Zhi92956–959.

  • Google Scholar

• 14

  DingerM. E.AmaralP. P.MercerT. R.PangK. C.BruceS. J.GardinerB. B.et al (2008). Long noncoding RNAs in mouse embryonic stem cell pluripotency and differentiation.Genome Res.181433–1445. 10.1101/gr.078378.108

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15

  EllingR.ChanJ.FitzgeraldK. A. (2016). Emerging role of long noncoding RNAs as regulators of innate immune cell development and inflammatory gene expression.Eur. J. Immunol.46504–512. 10.1002/eji.201444558

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16

  EstellerM. (2011). Non-coding RNAs in human disease.Nat. Rev. Genet.12861–874. 10.1038/nrg3074

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17

  FitzgeraldK. A.CaffreyD. R. (2014). Long noncoding RNAs in innate and adaptive immunity.Curr. Opin. Immunol.26140–146. 10.1016/j.coi.2013.12.001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  GambineriE.TorgersonT. R.OchsH. D. (2003). Immune dysregulation, polyendocrinopathy, enteropathy, and X-linked inheritance (IPEX), a syndrome of systemic autoimmunity caused by mutations of FOXP3 a critical regulator of T-cell homeostasis.Curr. Opin. Rheumatol.15430–435.

  • Pubmed Abstract
  • Google Scholar

• 19

  GarciaR. A.PereiraM. R.MaesterT. C.De Macedo LemosE. G. (2015). Investigation, expression, and molecular modeling of ORF2 a metagenomic lipolytic enzyme.Appl. Biochem. Biotechnol.1753875–3887. 10.1007/s12010-015-1556-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20

  GhoshM.McauliffeB.SubramaniJ.BasuS.ShapiroL. H. (2012). CD13 regulates dendritic cell cross-presentation and T cell responses by inhibiting receptor-mediated antigen uptake.J. Immunol.1885489–5499. 10.4049/jimmunol.1103490

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21

  GuoJ.CuiL.LuQ.ZhangY.LiuQ.WangX.et al (2018). Cyadox regulates the transcription of different genes by activation of the PI3K signaling pathway in porcine primary hepatocytes.J. Cell. Biochem.10.1002/jcb.28037 [Epub ahead of print].

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22

  HuangY.ZhangJ.HouL.WangG.LiuH.ZhangR.et al (2017). LncRNA AK023391 promotes tumorigenesis and invasion of gastric cancer through activation of the PI3K/Akt signaling pathway.J. Exp. Clin. Cancer Res.36:194. 10.1186/s13046-017-0666-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23

  KasprowiczD. J.SmallwoodP. S.TyznikA. J.ZieglerS. F. (2003). Scurfin (FoxP3) controls T-dependent immune responses in vivo through regulation of CD4+ T cell effector function.J. Immunol.1711216–1223.

  • Pubmed Abstract
  • Google Scholar

• 24

  KongL.ZhangY.YeZ. Q.LiuX. Q.ZhaoS. Q.WeiL.et al (2007). CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine.Nucleic Acids Res.35W345–W349. 10.1093/nar/gkm391

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25

  KoppF.MendellJ. T. (2018). Functional classification and experimental dissection of long noncoding RNAs.Cell172393–407. 10.1016/j.cell.2018.01.011

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26

  KramerW.GirbigF.CorsieroD.PfenningerA.FrickW.JahneG.et al (2005). Aminopeptidase N (CD13) is a molecular target of the cholesterol absorption inhibitor ezetimibe in the enterocyte brush border membrane.J. Biol. Chem.2801306–1320. 10.1074/jbc.M406309200

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27

  KwokZ. H.TayY. (2017). Long noncoding RNAs: lincs between human health and disease.Biochem. Soc. Trans.45805–812. 10.1042/bst20160376

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28

  LemlerD. J.BrochuH. N.YangF.HarrellE. A.PengX. (2017). Elucidating the role of host long non-coding rna during viral infection: challenges and paths forward.Vaccines5:37. 10.3390/vaccines5040037

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29

  LiX.WuZ.FuX.HanW. (2014). lncRNAs: insights into their function and mechanics in underlying disorders.Mutat. Res. Rev. Mutat. Res.7621–21. 10.1016/j.mrrev.2014.04.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30

  LiZ.LiX.JiangC.QianW.TseG.ChanM. T. V.et al (2018). Long non-coding RNAs in rheumatoid arthritis.Cell Prolif.51:e12404. 10.1111/cpr.12404

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31

  LiangH.ZenK.ZhangJ.ZhangC. Y.ChenX. (2013). New roles for microRNAs in cross-species communication.RNA Biol.10367–370. 10.4161/rna.23663

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32

  LiuQ.LeiZ.ZhouK.YuH.LiuS.SunQ.et al (2018). N-O reduction and ROS-Mediated AKT/FOXO1 and AKT/P53 pathways are involved in growth promotion and cytotoxicity of cyadox.Chem. Res. Toxicol.311219–1229. 10.1021/acs.chemrestox.8b00194

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33

  MatsushimaY.IshikawaO.KurosawaM.MiyachiY. (2000). Stem cell factor and IL-6 do not promote complete maturation of human cultured mast cells from umbilical cord blood cells: an ultrastructural study.J. Dermatol. Sci.244–13.

  • Google Scholar

• 34

  MetzM.MaurerM. (2007). Mast cells–key effector cells in immune responses.Trends Immunol.28234–241. 10.1016/j.it.2007.03.003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  MorigglR.TophamD. J.TeglundS.SexlV.MckayC.WangD.et al (1999). Stat5 is required for IL-2-induced cell cycle progression of peripheral T cells.Immunity10249–259. 10.1016/S1074-7613(00)80025-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36

  MowelW. K.KotzinJ. J.MccrightS. J.NealV. D.Henao-MejiaJ. (2018). Control of immune cell homeostasis and function by lncRNAs.Trends Immunol.3955–69. 10.1016/j.it.2017.08.009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37

  NakajimaH.LiuX.-W.Wynshaw-BorisA.RosenthalL. A.ImadaK.FinbloomD. S.et al (1997). An indirect effect of stat5a in IL-2–induced proliferation: a critical role for Stat5a in IL-2–mediated IL-2 receptor α chain Induction.Immunity7691–701. 10.1016/S1074-7613(00)80389-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38

  NeflaM.SudreL.DenatG.PriamS.Andre-LerouxG.BerenbaumF.et al (2015). The pro-inflammatory cytokine 14-3-3epsilon is a ligand of CD13 in cartilage.J. Cell. Sci.1283250–3262. 10.1242/jcs.169573

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39

  NohJ. H.KimK. M.MccluskyW. G.AbdelmohsenK.GorospeM. (2018). Cytoplasmic functions of long noncoding RNAs.Wiley Interdiscip. Rev. RNA9:e1471. 10.1002/wrna.1471

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40

  OnoM.YaguchiH.OhkuraN.KitabayashiI.NagamuraY.NomuraT.et al (2007). Foxp3 controls regulatory T-cell function by interacting with AML1/Runx1.Nature446685–689. 10.1038/nature05673

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  ReinholdD.BitonA.GoihlA.PieperS.LendeckelU.FaustJ.et al (2007). Dual inhibition of dipeptidyl peptidase IV and aminopeptidase N suppresses inflammatory immune responses.Ann. N.Y. Acad. Sci.1110402–409. 10.1196/annals.1423.042

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42

  RobertsA. P.LewisA. P.JoplingC. L. (2011). The role of microRNAs in viral infection.Prog. Mol. Biol. Transl. Sci.102101–139. 10.1016/b978-0-12-415795-8.00002-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43

  RosenbergerC. M.PodyminoginR. L.NavarroG.ZhaoG. W.AskovichP. S.WeissM. J.et al (2012). miR-451 regulates dendritic cell cytokine responses to influenza infection.J. Immunol.1895965–5975. 10.4049/jimmunol.1201437

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44

  ScottM. E.MaY.KuzmichL.MoscickiA. B. (2009). Diminished IFN-gamma and IL-10 and elevated Foxp3 mRNA expression in the cervix are associated with CIN 2 or 3.Int. J. Cancer1241379–1383. 10.1002/ijc.24117

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45

  SelliahN.ZhangM.DesimoneD.KimH.BrunnerM.IttenbachR. F.et al (2006). The gammac-cytokine regulated transcription factor, STAT5 increases HIV-1 production in primary CD4 T cells.Virology344283–291. 10.1016/j.virol.2005.09.063

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46

  ShiL.ZhangN.LiuH.ZhaoL.LiuJ.WanJ.et al (2018). Lysyl oxidase inhibition via beta-aminoproprionitrile hampers human umbilical vein endothelial cell angiogenesis and migration in vitro.Mol. Med. Rep.175029–5036. 10.3892/mmr.2018.8508

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47

  SoderbergC.GiugniT. D.ZaiaJ. A.LarssonS.WahlbergJ. M.MollerE. (1993). CD13 (human aminopeptidase N) mediates human cytomegalovirus infection.J. Virol.676576–6585.

  • Pubmed Abstract
  • Google Scholar

• 48

  ThomasR.DavisL. S.LipskyP. E. (1994). Rheumatoid synovium is enriched in mature antigen-presenting dendritic cells.J. Immunol.1522613–2623.

  • Pubmed Abstract
  • Google Scholar

• 49

  TriulziT.TagliabueE.BalsariA.CasaliniP. (2013). FOXP3 expression in tumor cells and implications for cancer progression.J. Cell Physiol.22830–35. 10.1002/jcp.24125

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50

  WangH.PengW.OuyangX.LiW.DaiY. (2012). Circulating microRNAs as candidate biomarkers in patients with systemic lupus erythematosus.Transl. Res.160198–206. 10.1016/j.trsl.2012.04.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 51

  WangX.LiuQ.IhsanA.HuangL.DaiM.HaoH.et al (2012). JAK/STAT pathway plays a critical role in the proinflammatory gene expression and apoptosis of RAW264.7 cells induced by trichothecenes as DON and T-2 toxin.Toxicol. Sci.127412–424. 10.1093/toxsci/kfs106

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52

  WangP.XueY.HanY.LinL.WuC.XuS.et al (2014). The STAT3-binding long noncoding RNA lnc-DC controls human dendritic cell differentiation.Science344310–313. 10.1126/science.1251456

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53

  WangZ.XiJ.HaoX.DengW.LiuJ.WeiC.et al (2017). Red blood cells release microparticles containing human argonaute 2 and miRNAs to target genes of Plasmodium falciparum.Emerg. Microb. Infect6:e75. 10.1038/emi.2017.63

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54

  WuZ.LiuX.LiuL.DengH.ZhangJ.XuQ.et al (2014). Regulation of lncRNA expression.Cell Mol. Biol. Lett.19561–575. 10.2478/s11658-014-0212-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55

  YamadaH.ItohM.HiratsukaI.HashimotoS. (2014). Circulating microRNAs in autoimmune thyroid diseases.Clin. Endocrinol.81276–281. 10.1111/cen.12432

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56

  YangX. P.GhoreschiK.Steward-TharpS. M.Rodriguez-CanalesJ.ZhuJ.GraingerJ. R.et al (2011). Opposing regulation of the locus encoding IL-17 through direct, reciprocal actions of STAT3 and STAT5.Nat. Immunol.12247–254. 10.1038/ni.1995

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57

  YuR.ZhangY.LuQ.CuiL.WangY.WangX.et al (2018). Differentially expressed genes in response to cyadox in swine liver analyzed by DDRT-PCR.Res. Vet. Sci.11872–78. 10.1016/j.rvsc.2018.01.014

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58

  ZhangJ.ZhouY.WuY.MaL.FanY.KangX.et al (2012). Isolation and characterization of a novel noncoding RNA from nickel-induced lung cancer.Biol. Trace Elem. Res.150258–263. 10.1007/s12011-012-9460-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 59

  ZotzJ. S.WolbingF.LassnigC.KauffmannM.SchulteU.KolbA.et al (2016). CD13/aminopeptidase N is a negative regulator of mast cell activation.FASEB J.302225–2235. 10.1096/fj.201600278

  • Pubmed Abstract
  • CrossRef
  • Google Scholar