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

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 Guo et al., 2018) and with the potential to serve as a feed additive . 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 . 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 . 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 fulllength 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.

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).
The ORF of CRNG was predicted by bioinformatics tools in NCBI 1 ( 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 2 , it suggested that the gene CRNG was a non-coding RNA (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.

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.

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.

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, FIGURE 2 | Coding ability prediction of CRNG full-length and its reverse complementary strand. The score of coding ability was minus, explaining that CRNG was a non-coding gene.
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 FIGURE 3 | Relative expression levels of CRNG in different piglet tissues as measured by reverse transcription quantitative polymerase chain reaction (RT-qPCR). β-actin was used as a reference gene, liver was used as control. LncRNA CRNG was the highest in liver, high in duodenum and jejunum, low in kidney, adrenal gland, hypothalamus and pituitary, very low in other tissues. All experiments were performed in biological triplicates with three technical replicates. 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.

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 2 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 R 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 .

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 .

RACE and 5 RACE PCR
The 3 RACE and 5 RACE experiment was performed in full accordance with the classical RACE method of PCR Frontiers in Pharmacology | www.frontiersin.org 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 0 via the phanta max superfidelity DNA polymerase (Vazyme, China). And the inner PCR was carried out with gene specific primer CRNG2 and Q 1 via taq DNA Polymerase (Mg 2+ 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 0 and Q T via the phanta max superfidelity DNA polymerase (Vazyme, China), and the inner PCR was carried out with gene specific primer GSP3 and Q 1 via taq DNA Polymerase (Mg 2+ 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 , 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 -GGTGACGGCACGTTTGACC TCGTTGAGATGGTG-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.

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 antigenpresenting (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 immunityregulation 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 H 17 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 FIGURE 8 | Possible mechanisms of lncRNA regulating immunity and inflammation. Boosted lncRNA CRNG could increase FOXP3 and KITLG expression and inhibit ANPEP and STAT5A expression, thus promoting pro-inflammation cytokine and suppressing anti-inflammation cytokine and virus infection and production. Moreover, CRNG could enhance the expression of miR-451, thus protecting host from pathogens infection, and indirectly regulate the homeostatic of IL-6 levels. Each color line represents a signal pathway. ANPEP, Alanyl aminopeptidase, membrane; CRNG, Cyadox-related novel gene; FOXP3, Forkhead box protein 3; KITLG, Kit ligand; miR-451, microRNA 451, STAT5A: Signal transducer and activator of transcription 5A. (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 . miR-451 also was associated with autoimmune disease, which was increased in serum levels 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.

ETHICS STATEMENT
The current study was approved by the Ethical Committee of the Faculty of Veterinary Medicine (Huazhong Agricultural University).

AUTHOR CONTRIBUTIONS
QLu, LL, and AH generated, analyzed, and interpreted the data, and prepared the manuscript. LC, YZ, QLi, XW, YW, and ZL generated the format and grammar modification. ZY and MD generated the idea, designed the study, analyzed and interpreted the data, and wrote the manuscript.

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar. 2019.00539/full#supplementary-material DATA SHEET S1 | The result of comparing CRNG with NCBI database.
DATA SHEET S2 | The expressions of 78 genes upon the interfered CDRG via siRNA2, including 32 upregulated genes and 46 downregulated genes, show in excel.
DATA SHEET S3 | The silencing or overexpression efficiency after transfecting the cells with siRNA or CRNG-expression vector, the silencing efficiency of CRNG was showed in Figure 2 and CRNG induced by cyadox in PK-15 cells was shown in Figure 3. The siRNA sequences of CRNG were showed in Table 1.