Skip to main content

Analysis of long non-coding RNAs in neonatal piglets at different stages of porcine deltacoronavirus infection



PDCoV (Porcine Deltacoronavirus) is a novel porcine coronavirus that causes intestinal necrosis of piglets, thinning of the intestinal wall and severe villus atrophy in the small intestine. PDCoV is a highly contagious infectious disease characterized by diarrhea, dehydration and vomiting. It has been reported that lncRNA has a significant effect on viral replication and increased or decreased virulence. At present, there is almost no research on lncRNA related to PDCoV infection. With the development of the research, a large number of lncRNAs related to PDCoV infection have been discovered. Identifying the role of these lncRNAs in the infection process facilitates the screening of diagnostically significant biomarkers.


Using high throughput sequencing to screen differentially expressed long non-coding RNA (lncRNA) during PDCoV infection, we identified 99, 41 and 33 differentially expressed lncRNAs in the early, middle and late stages of infection, respectively. These lncRNAs were involved in glycolysis / gluconeogenesis, histidine metabolism and pentose and Chloroalkane and chloroalkene degradation pathway. We obtained expression data of miRNAs, lncRNAs and mRNAs during PDCoV infection and constructed and investigated an interaction network. The qRT-PCR validation results of 6 differentially expressed lncRNAs were consistent with RNA-Seq results.


This study is the first to examine differentially expressed lncRNAs after PDCoV infection of piglets. These results can provide new insights into PDCoV infection and antiviral strategies.


Porcine deltacoronavirus (PDCoV) belongs to the newly identified deltacoronavirus genus in the Coronaviridae family. This family possesses single-stranded positive-strand RNA genomes contained in an nucleocapside, so PDCoV is an enveloped virus. Intestinal infections in pigs result in diarrhea, vomiting and dehydration and are highly lethal. In 2009, 17 cases of swine diarrhea were positive for PDCoV and were linked to the deltacoronavirus [1]. The virus appeared again in 2014 causing large-scale outbreaks of pig diarrhea in the United States [2]. A retrospective study in 2015 identified PDCoV from samples collected from four Chinese provinces between 2004 and 2014 indicated an infection rate of 6.5% [3]. Subsequently, Thailand, Laos, Vietnam and the Philippines successively reported the existence of PDCoV [4, 5].

PDCoV has a broad host range and interacts with the catalytic domain of receptor aminopeptidase N (APN) in cat, chicken and human cells through the S1 domain B of the S protein, even though PDCoV infections in animals other than pigs have not been reported [6]. Ectopic expression of porcine aminopeptidase N (pAPN) in non-susceptible cells converted them to PDCoV-susceptible although pAPN is not a critical receptor but is an important accessory factor for infection [7]. PDCoV infection also affects the host’s innate immunity. PDCoV infection inhibited Interferon-beta (IFN-β) production by inhibiting nuclear factor-κB (NF-κB) and interferon regulation factor 3(IFN3) activation and inhibiting IFN-β promoter activation in the RIG-I pathway by suppressing IFN3 [8].

Long non-coding RNAs (lncRNA) are RNAs > 200 nucleotides in length with no protein encoding function. These RNAs play key roles in numerous biological functions as well as human disease. LncRNAs are involved as regulatory factors and are involved in epigenetic and transcriptional and post-transcriptional regulation. The functions of mammalian lncRNAs have been linked to metabolism, apoptosis and cell proliferation as well as immune responses including regulation of pattern recognition receptors and inflammation [9,10,11,12].

LncRNAs play crucial roles in host-virus interactions that include immune responsiveness and viral infections can induce their expression. For example, the differential expression of > 4800 lncRNAs were found in rhabdomyosarcoma cells after Enterovirus 71 infection [13]. In addition, > 3000 differentially expressed genes (DEG) were associated with lncRNAs during infectious salmon anaemia (ISA) virus infection and most were regulated in response to infection [14]. A total of 1236 lncRNA transcripts were differentially regulated at different stages of bovine viral diarrhea virus infected Madin Darby bovine kidney (MDBK) cells [15].

In the current study, we analyzed lncRNAs that were differentially expressed after PDCoV infection of piglets to examine whether lncRNAs are involved in PDCoV infection.


Determination of tissue viral loads

We plaque-purified the virus that we used for swine testicular cells (ST) cell infections and the virus with titer 1.0 × 108.6 TCID50/ml was used for experimental infection of piglets. We found distinct differences in intestinal pathology between the experimental and the control groups at 2, 4 and 11 day post-infection (dpi) (Fig. 1). Piglet tissue samples taken at 2, 4 and 11 dpi were aseptically collected and the PDCoV RNA viral loads in each organ were determined by qPCR. The virus was found primarily in the small intestine. At mid-infection, viral loads in all organs except the kidney reached a maximum for all three periods. The heart, spleen, kidney, duodenum, jejunum, ileum, mesenteric lymph nodes and tonsils were all positive for PDCoV. At the late stage of infection, all organs except the ileum were PDCoV-negative (Fig. 2). The ileum was the most serious site of infection, so we used these samples for high throughput sequencing.

Fig. 1

Abdominal anatomy of pigs at different infection stages of PDCoV infection. a, Experimental group at 2dpi; b, Control group at 2dpi; c, Experimental group at 4dpi; d, Control group at 4dpi; e, Experimental group at 11dpi; f, Control group at 11dpi

Fig. 2

qRT-PCR results of tissue viral load after PDCoV infection in piglets

Identification of lncRNAs in PDCoV -infected piglets by RNA-Seq

The samples we used for RNA-Seq were taken in the pre-onset (2 dpi), mid-onset (4 dpi) and late-onset (11 dpi). During the experiment, we chose three piglets from eighteen experimental piglets during the corresponding infection period for dissection, and the control group also dissected three piglets. The samples we used for RNA-seq were obtained from six piglets dissected at each stage. We analyzed the RNA-Seq data for 18 samples and each contained between 53,173,880 and 94,027,536 raw data points. After filtering, 52,830,072 to 93,138,526 clean data points were obtained for each sample and after assembly resulted in 2130 novel lncRNAs. In this group there were 173 differentially expressed lncRNAs identified in the three infection phases at a fold change > 2 and P < 0.05. The early, middle and late periods included 99 (21 up, 78 down), 41 (9 up, 32 down) and 33 (10 up, 23 down) lncRNAs, respectively (Table 1).

Table 1 Quantitative analysis of differential expression of lncRNA between experimental group and control group in different periods

Target gene prediction and lncRNA analysis

LncRNAs are regulators of protein coding genes that lie near their genomic locations. We analyzed these RNAs for all protein coding genes within 100 kb of the lncRNA as potential cis-regulatory targets (Additional file 1).

Gene ontology (GO) analysis of DEGs indicated that 7, 7 and 9 terms were significantly enriched (P < 0.05) during early, middle and late infection periods, respectively. Molecular functions in the early period involved ferric-chelate reductase and oxidoreductase activity and metal ion oxidization. The middle period had terms related to oxidoreductase and prenyltransferase activity as well as positive regulation of activated T cell proliferation and catalytic activity. The late period terms were related to oxidoreductase activity, chromosome organization involved in the meiotic cell cycle and small molecule binding. Target genes related to oxidoreductase activity were found in all three periods (Additional file 2).

To elucidate the role of lncRNAs in the virus-host relationship, we performed Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of lncRNA target genes during each challenge period. We found that 6 pathways in the early period, 11 pathways in the middle and 14 pathways in the late period that were significantly enriched (P < 0.05). The temporal expression pattern of these pathways was related to both anabolic and catabolic functions. Anabolic functions included ascorbate and aldarate metabolism, histidine, β-alanine, glyoxylate and dicarboxylate metabolism. Catabolic pathways involved metabolism of limonene, pinene, chloroalkane and chloroalkene (Fig. 3a, b and c).

Fig. 3

KEGG pathway enrichment analysis of target genes. a, The enrichment pathway of genes in the early infection stage; b, The enrichment pathway of genes in the middle infection stage; c, The enrichment pathway of genes in the late infection stage. Dot diameter is proportional to the number of differential genes and color depth is proportional to significance. Abscissa, enrichment ratio; ordinate, the different pathways. (The signal pathway shown in the figure that represents in the top 20 p-value ranking)

lncRNA-miRNA-mRNA interaction network analysis

CeRNAs (competitive endogenous RNA) are a class of functionally defined RNAs possessing miRNA binding sites that can compete with miRNA and inhibit their regulation of target genes [16]. Some lncRNAs in the network we constructed function as ceRNAs and we were able to construct a network based on these RNAs. The network was composed of 11 lncRNA nodes, 33 mRNA nodes and 10 miRNA nodes. The expression level of lncRNA is significantly up-regulated or down-regulated, and its target miRNA is down-regulated or up-regulated. The target gene expression level of miRNA is consistent with lncRNA. This suggests that there may be significant competitive RNA during infection. During the early infection stage, a total of 3 differentially expressed miRNAs were targeted by 4 lncRNAs and 9 mRNAs in the network. MSTRG.20658.3 and MSTRG.24537.4 can competitively bind to ssc-miR-194b-3p with mRNA. MSTRG.16902.19 and MSTRG.43450.4 can competitively bind to ssc-miR-885-5p and MSTRG.16902.19 can also competitively bind to ssc-miR-7857-3p (Fig. 4a). In the middle period, a total of 6 differentially expressed miRNAs were targeted by 5 lncRNAs and 18 mRNAs in the network. We found that MSTRG.45362.10 was capable of competitive binding to both ssc-miR-490-5p and ssc-miR-216, and MSTRG.31817.2 was capable of competitive binding to ssc-miR-215, ssc-miR-874 and ssc-miR-184. MSTRG.35510.2, MSTRG.22938.2 and MSTRG.53204.2 can competitively bind to ssc-miR-133b, ssc-miR-184, ssc-miR-874, respectively (Fig. 4b). In the late stage, a total of 2 differentially expressed miRNAs were targeted by 2 lncRNAs and 6 mRNAs in the network. MSTRG.31541.4 was capable of competitive binding to ssc-miR-196b-5p, and MSTRG.41555.2 can competitively bind to ssc-miR-9785-5p (Fig. 4c). The biological significance of ceRNA network in PDCoV infection was also reflected by topological structures including hubs and nodes as well as direct connections.

Fig. 4

The lncRNA-miRNA-mRNA interaction network. a, pre-infection; b, middle stage and c, late stage of infection. Triangle, lncRNA; rectangle, miRNA; circle, mRNA. Red, up-regulation; green, down-regulation

Validation of differentially expressed lncRNAs

The RNA-Seq analysis indicated that the expression level of many lncRNAs were significantly changed in the experimental group compared with the control group. In the mid-infection period, we randomly selected 6 lncRNAs (MSTRG.35510.2, MSTRG.4915.1, MSTRG.16902.37, MSTRG.45362.10, MSTRG.16541.3, MSTRG.45362.9) for expression level verification. Our validation results were consistent with the RNA-Seq results, so the results of RNA-Seq were credible. (Fig. 5).

Fig. 5

Verification of RNA-Seq results. Six differentially expressed lncRNAs in the middle stages of infection were randomly selected for real-time qRT-PCR experiments to verify whether their expression levels were consistent with the RNA-Seq results. a, the expression verification result of six lncRNAs (MSTRG.35510.2, MSTRG.4915.1, MSTRG.16902.37, MSTRG.45362.10, MSTRG.16541.3, MSTRG.45362.9); b, the RNA-Seq results of lncRNAs (MSTRG.35510.2, MSTRG.4915.1, MSTRG.16902.37, MSTRG.45362.10, MSTRG.16541.3, MSTRG.45362.9). (n = 3; *p < 0.05, **p < 0.01, ***p < 0.001)


PDCoV infects pigs of all ages, but primarily causes diarrhea in newborn pigs. Clinically, PDCoV infection is similar to porcine intestinal coronavirus, but PDCoV infection has a wider tissue tropism and can be detected in organs other than the digestive tract. This suggests a complex pathogenic mechanism for this virus and an in-depth understanding of its pathogenic and immune mechanisms is necessary for infection control.

Current lncRNA research has focused on human medicine including cardiovascular disease and cancer [17,18,19]. In livestock and poultry, examination of lncRNAs is in its infancy and existing research has focused on muscle, bone and embryonic development as well as fat metabolism [20,21,22]. Additionally, these studies have focused on lncRNA regulation of protein-coding genes.

Recent studies have shown that viral infections can induce lncRNAs to promote or inhibit viral responses. The lncRNA NEAT1 can up-regulate anti-HIV factors during infection and promote human immunodeficiency virus 1(HIV-1) replication [23]. The lncRNA ACOD1 enhances the replication of multiple viruses in both mouse and human cells [24]. However, an examination of lncRNA expression during PDCoV infections was lacking. The current study is the first to use comprehensive deep-sequencing technology that implicates lncRNAs in the response to PDCoV infection in pigs.

Our RNA-Seq data can assist in understanding the mechanism of action of differentially expressed lncRNAs at different stages of PDCoV infection. We established a lncRNA gene library that was generated during the early, middle and late stages of PDCoV infection. We identified 173 differentially expressed lncRNAs and 2130 novel lncRNAs. The greatest number of differentially expressed lncRNAs were found during the early stage of infection (2 dpi). The number of down regulated was >up-regulated lncRNAs. In addition, we found lncRNA MSTRG.18455 was significantly down regulated (− 7-fold) and its target gene IGF1 was significantly enriched. The insulin-like growth factor 1 (IGF1) is a member of the growth and development promoting signaling system and the main determinant of animal growth [25, 26].

Porcine enteroviruses can enter the digestive system through the mouth and subsequently attach to the intestinal villi. This causes villus atrophy resulting in diarrhea, dehydration, vomiting and weight loss. Our results from target mRNA pathway analysis revealed that pre-challenge target mRNAs of lncRNA were enriched for the signaling pathway of glyoxylate and dicarboxylate metabolism, limonene and pinene degradation, chloroalkane and chloroalkene degradation as well as glycolysis / gluconeogenesis. During the middle stages of infection, limonene and pinene degradation, glycolysis / gluconeogenesis, ascorbate and aldarate metabolism, histidine metabolism signaling pathway were prominent and late-stage target mRNAs were primarily concentrated on histidine metabolism, glyoxylate and dicarboxylate metabolism, beta-Alanine metabolism. It can be seen that some differentially expressed lncRNAs are related to the metabolism of organisms. The GO analysis indicated that at the middle stage of infection, these target genes were significantly enriched in the biological process category including regulation of activated T cell proliferation. Our future work will include determining lncRNA mechanisms of T cell activation and proliferation. This will assist in determining how the immune system is compromised by PDCoV infection in pigs.

In this study, we constructed a lncRNA-miRNA-mRNA interaction network containing a ceRNA network. Interestingly, both lncRNA and mRNA in the ceRNA network were negatively correlated with miRNA. Previous studies have confirmed that ceRNAs can act as miRNA ‘sponges’ and this is especially important for cancer and tumor diseases [27, 28]. By analyzing the miRNAs in the ceRNA network we found that the target miRNAs ssc-miR-885-5p, ssc-miR-490-5p, ssc-miR-196b-5p and ssc-miR-133b of lncRNAs were all up regulated. MiR-885-5p is a direct regulator of the IGF1 receptor that in combination with p73 regulate emergence of aggressive cancer stem-like features [29]. Decreased expression of miR-133b is associated with poor survival and increased metastasis in colorectal cancer [30]. miR-490-5p inhibits cell proliferation, migration and invasion but miR-490-5p can promote apoptosis of Human hepatoma (Hep3B) cells by inhibiting Roundabout Guidance Receptor 1 (ROBO1) [31]. MiR-196b-5p overexpression may also be associated with a risk of conversion of myelodysplastic syndrome (MDS) to acute myeloid leukemia (AML) [32]. The target miRNAs ssc-miR-194b-3p, ssc-miR-184, ssc-miR-215, ssc-miR-874 were all down regulated. miR-874 and miR-215 can act as a tumor suppressor [33, 34]. Therefore, this network plays a role in the expression of immune-related genes and the stimulation of host immune responses.

Our research provides a scientific reference for the lncRNAs that regulate PDCoV replication that can assist studies of ceRNA and PDCoV infection. This work can also aid the development of effective drugs and genetic engineering of pigs with PDCoV resistance.


In this study, we provide the first analysis of differentially expressed lncRNAs after PDCoV infection in piglets. we constructed a lncRNA-miRNA-mRNA interaction network. This study provides insights into the relationships between lncRNAs and PDCoV immune modulation. Future studies will address lncRNA functions in immune escape used by PDCoV.


Cell culture and virus

ST cells used for experimental PDCoV infections in the study were kept in the Poultry Laboratory of the College of Animal Science, South China Agricultural University and cultured at 37 °C in a humidified 5% CO2 atmosphere in Dulbecco’s modified Eagle’s medium (DMEM, HyClone, Logan, UT, USA) supplemented with 10% fetal bovine serum (FBS, Hyclone) [35]. ST cells and culture media were checked to ensure the absence of PDCoV, porcine epidemic diarrhea virus (PEDV), transmissible gastroenteritis virus (TGEV) by gel electrophoresis of RT-PCR products. The virus strain PDCoV-CHN-GD16–05 (GenBank 74 accession no.KY363868.1) [36] was isolated and preserved by our laboratory from the watery diarrhea feces of nursing piglets in Guangdong Province, China. ST monolayers at 80% confluency were infected at a multiplicity of infection (MOI) of 1.5. The cells were cultured in serum-free medium at 37 °C for 1 h and the medium was replaced with fresh culture medium containing 2% FBS. Infected cells were collected 36 h post-infection (hpi).

Determination of viral growth in tissues

The tissues we have taken include heart, liver, spleen, lung, kidney, duodenum, jejunum, ileum, mesenteric lymph nodes (MLN), inguinal lymph nodes and tonsils. Absolute quantification of viral RNA is done by real-time qRT-PCR. using the PCR primer pair 5′-TGGCTGATCCTCGCATCATGG-3’and 5′-GAGCGCATCCTTAAGTCT CTC-3′. One-step RT-PCR reactions were performed using an ABI PRISM 7500 (Applied Biosystems, Foster City, CA, USA) instrument according to the TAKARA company ‘PrimeScript™ RT reagent Kit with gDNA Eraser’ operating instructions. We used SYBR Green qPCR Super Mix (Invitrogen) for real-time qRT-PCR experiments based on amplification conditions. The following steps were used: denaturation at 95 °C for 10 mins and 40 cycles at 95 °C for 15 s, 60 °C for 30 s, 72 °C for 30 s and finally a melting curve.

Piglet challenge experiments

The piglets were purchased from a farm in Huanong Wen’s Co., Ltd. Thirty 5-day-old piglets lacking any overt signs of infection were selected after observation for 24 h to ensure they were not exposed to stress before infection. The animals were infected with 5 mL of a virus solution containing 1.0 × 108.6 TCID50/mL administered orally (18 piglets), and 12 control piglets were administered the same volume of DMEM. At the end of the experiment, the experimental animals showed weight loss, loss of appetite, and diarrhea, so they were euthanized. Experimental animal euthanasia method is based on the experimental animal management and practical technical manual. Intravenous injection of sodium pentobarbital at a dose of 90–100 mg/kg. Corpses were put into the septic tank, and they were fermented and used for fertilizer.

RNA-Seq analysis

Three piglets were euthanized and dissected at the early (2 dpi), middle stage (4 dpi) and the late stages of infection period (11 dpi). Intestinal tissues with the most obvious lesions at each stage were collected in triplicate, frozen in liquid nitrogen and stored at − 80 °C. For controls, intestinal tissue was taken from control animals at similar time points.

Total RNA was extracted using Trizol reagent (Invitrogen) and purified using the Qiagen RNeasy Mini Kit (Qiagen, Valencia, CA, USA) according to the manufacturer’s instructions. RNA concentration and integrity were measured using the Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). We prepared RNA sequencing libraries from small intestine samples and performed 150-bp paired-end sequencing using the Illumina HiSeq platform. RNA sequencing libraries were prepared from 2 μg of total RNA using the TruSeq Kit (Illumina, San Diego, CA, USA) with the following modification. Instead of purifying poly-A RNA using poly-dT primer beads, we removed ribosomal RNA using the Ribo-Zero rRNA Removal Kit (Illumina). All other steps were performed according to the manufacturer’s protocol.

RNA-Seq libraries were quality control analyzed and the average insert size was 200 to 300 bp. The library was sequenced using a Hiseq platform (Illumina) at Shanghai Personal Biotechnology (Shanghai, China). Raw data was filtered and high-quality data was obtained. At the same time, the quality of the clean data was checked using the Q20, Q30 and GC content.

The Stringtie algorithm ( was used to analyze expression levels of lncRNA transcripts. The coding potential of candidate lncRNAs was analyzed to obtain high-confidence lncRNA using the Coding Potential Calculator (, Coding-noncoding Index and Pfamscan (

LncRNAs are primarily located near genes with coding functions and the function of the associated lncRNA can be approximated by determining target gene function within 100 kb. GO ( and KEGG ( databases were used to determine functions and metabolic pathways of DEGs. GO terms and KEGG pathways with Q-values ≤0.05 were considered significantly enriched.

Construction of lncRNA-miRNA-mRNA interaction network

LncRNAs, miRNAs and mRNAs that were differentially expressed between the experimental and control group tissues were chosen for analysis. The differential expression of miRNAs and lncRNAs was identified using standard selection criteria at P < 0.05 and fold change > 2. We used miRanda and psRobot software to predict target miRNA for lncRNA and target mRNA for miRNA ( and Target mRNAs were separately analyzed. The ceRNA network was generated from the interaction network and used Pearson correlation coefficient (PCC) to calculate the correlation between mRNA and lncRNA and chose pairs with positive correlations (PCC > 0.99 and P < 0.05). Visualization of the lncRNA-miRNA-mRNA interaction network was constructed using Cytoscape software (version 3.0;

Screening and qRT-PCR validation of differentially expressed lncRNAs

Six lncRNAs (MSTRG.35510.2, MSTRG.4915.1, MSTRG.16902.37, MSTRG.45362.10, MSTRG.16541.3, MSTRG.45362.9) were randomly selected from the middle of the infection for real-time qRT-PCR validation, with glyceraldehyde-3-phosphate dehydrogenase (GAPDH) used as an endogenous control. Total RNA was extracted from uninfected and infected cells using a total RNA extraction kit (Tiangen Biotech, Beijing, China). cDNA was synthesized using a reverse transcription kit (Fisher Scientific, Pittsburg, PA, USA) with 2 μg of total RNA according to the manufacturer’s instructions. The primers were designed using Primer 5.0 software (Table 2). One-step real-time RT-PCR reactions were performed using a One Step PrimeScript RT-PCR kit (Takara) on an ABI PRISM 7500 instrument. The following steps were used: reverse transcription at 42 °C for 10 min, denaturation at 95 °C for 10 s and 40 cycles at 95 °C for 5 s, 55 °C for 20 s, 72 °C for 10 s and finally a melting curve. Each reaction was performed in triplicate. LncRNA expression levels were calculated based on the 2 -ΔΔCT method [37]. Expression levels were normalized to those of GAPDH.

Table 2 DNA primers used for qRT-PCR

Statistical analysis

RNA-Seq data was analyzed using statistical R and expression levels of lncRNAs were compared using the paired sample t-test. Data were expressed as the mean ± standard deviation from at least three independent experiments. SPSS 17.0 software package (SPSS, Chicago, IL, USA) was used to analyze the qRT-PCR data. P < 0.05 was considered statistically significant. Statistical differences between the control and PDCoV infected cells were analyzed using the Student’s t-test.



Acute myeloid leukemia


Aminopeptidase N


Competitive endogenous RNA


Differentially expressed genes


Dulbecco’s modified Eagle’s medium


days post-infection


Fetal bovine serum


Glyceraldehyde-3-phosphate dehydrogenase


Gene ontology


Human hepatoma


Human immunodeficiency virus


hours post-infection


Interferon regulation factor 3




Insulin-like growth factor 1


Infectious salmon anaemia


Kyoto Encyclopedia of Genes and Genomes


Long non-coding RNA


Madin Darby bovine kidney


Myelodysplastic syndrome


A multiplicity of infection


Nuclear factor-κB


Porcine aminopeptidase N


Pearson correlation coefficient


Porcine deltacoronavirus


Porcine epidemic diarrhea virus


Roundabout Guidance Receptor 1


Swine testicular cells


Transmissible gastroenteritis virus


  1. 1.

    Woo PCY, Lau SKP, Lam CSF. Discovery of seven novel mammalian and avian coronaviruses in the genus Deltacoronavirus supports bat coronaviruses as the gene source of Alphacoronavirus and Betacoronavirus and avian coronaviruses as the gene source of Gammacoronavirus and Deltacoronavir. J Virol. 2012;86(7):3995–4008.

    CAS  Article  Google Scholar 

  2. 2.

    Wang L, Byrum B, Zhang Y. Detection and genetic characterization of deltacoronavirus in pigs, Ohio, USA, 2014. Emerg Infect Dis. 2014;20(7):1227.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Dong N, Fang L, Zeng S, Sun Q, Chen H, Xiao S. Porcine Deltacoronavirus in mainland China. Emerg Infect Dis. 2015;21(12):2254–5.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Lorsirigool A, Saeng-Chuto K, Temeeyasen G, Madapong A, Tripipat T, Wegner M, et al. The first detection and full-length genome sequence of porcine deltacoronavirus isolated in Lao PDR. Arch Virol. 2016;161(10):2909–11.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Madapong A, Saeng-Chuto K, Lorsirigool A, Temeeyasen G, Srijangwad A, Tripipat T, et al. Complete genome sequence of porcine Deltacoronavirus isolated in Thailand in 2015. Genome Announc. 2016;4(3):e408–16.

    Article  Google Scholar 

  6. 6.

    Li W, Hulswit RJG, Kenney SP, Widjaja I, Jung K, Alhamo MA, et al. Broad receptor engagement of an emerging global coronavirus may potentiate its diverse cross-species transmissibility. Pnas. 2018;115(22):E5135–43.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Zhu X, Liu S, Wang X, Luo Z, Shi Y, Wang D, et al. Contribution of porcine aminopeptidase N to porcine deltacoronavirus infection. Emerging Microbes & Infections. 2018;7(1):65.

    CAS  Article  Google Scholar 

  8. 8.

    Luo J, Fang L, Dong N, Fang P, Ding Z, Wang D, et al. Porcine deltacoronavirus (PDCoV) infection suppresses RIG-I mediated interferon-β production. Virology. 2016;495:10–7.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Zhu M, Liu J, Xiao J, Yang L, Cai M, Shen H, et al. Lnc-mg is a long non-coding RNA that promotes myogenesis. Nat Commun. 2017;8:14718.

    Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Rapicavoli NA, Qu K, Zhang J, Mikhail M, Laberge RM, Chang HY. A mammalian pseudogene lncRNA at the interface of inflammation and anti-inflammatory therapeutics. Elife. 2013;2(2):e00762.

    Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Xiong Y, Yuan J, Zhang C, Zhu Y, Kuang X, Lan L, et al. The STAT3-regulated long non-coding RNA Lethe promote the HCV replication. Biomed Pharmacother. 2015;72:165–71.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Landeras-Bueno S, Ortín J. Regulation of influenza virus infection by long noncoding RNAs. Virus Res. 2016;212:78–84.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Yin Z, Guan D, Fan Q, Su J, Zheng W, Ma W, et al. lncRNA expression signatures in response to enterovirus 71 infection. Biochemical & Biophysical Research Communications. 2013;430(2):629299 633.

    CAS  Article  Google Scholar 

  14. 14.

    Boltaña S, Valenzuela-Miranda D, Aguilar A, Mackenzie S, Gallardo-Escárate C. Long noncoding RNAs (lncRNAs) dynamics evidence immunomodulation during ISAV-infected Atlantic salmon (Salmo salar). Sci Rep. 2016;6:22698.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Ma Q, Li L, Tang Y, Fu Q, Liu S, Hu S, et al. Analyses of long non-coding RNAs and mRNA profiling through RNA sequencing of MDBK cells at different stages of bovine viral diarrhea virus infection. Res Vet Sci. 2017;115:508–16.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Fang L, Du WW, Yang X, Chen K, Ghanekar A, Levy G, et al. Versican 3′-untranslated region (3′-UTR) functions as a ceRNA in inducing the development of hepatocellular carcinoma by regulating miRNA activity. FASEB J. 2013;27:907–19.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Ma MZ, Chu BF, Zhang Y, Weng MZ, Qin YY, Gong W, et al. Long non-coding RNA CCAT1 promotes gallbladder cancer development via negative modulation of miRNA-218-5p. Cell Death Dis. 2015;6(1):e1583.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Xia T, Liao Q, Jiang X, Shao Y, Xiao B, Xi Y, et al. Long noncoding RNA associated-competing endogenous RNAs in gastric cancer. Sci Rep. 2014;4:6088.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Yu X, Tang W, Yang Y, Tang L, Dai R, Pu B, et al. Long noncoding RNA NKILA enhances the anti-cancer effects of baicalein in hepatocellular carcinoma via the regulation of NF-κB signaling. Chem Biol Interact. 2018;285:48–58.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Li T, Wang S, Wu R, Zhou X, Zhu D, Zhang Y. Identification of long non-protein coding RNAs in chicken skeletal muscle using next generation sequencing. Genomics. 2012;99(5):292–8.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Roeszler KN, Itman C, Sinclair AH, Smith CA. The long non-coding RNA, MHM, plays a role in chicken embryonic development, including gonadogenesis. Dev Biol 2012; 318 366(2):317–326. doi:

    CAS  Article  Google Scholar 

  22. 22.

    Divoux A, Karastergiou K, Xie H, Guo W, Perera RJ, Fried SK, et al. Identification of a novel lncRNA in gluteal adipose tissue and evidence for its positive effect on preadipocyte differentiation. Obesity. 2014;22(8):1781–5.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Zhang Q, Chen CY, Yedavalli VS, Jeang KT. NEAT1 long noncoding RNA and Paraspeckle bodies modulate HIV-1 posttranscriptional expression. Mbio. 2013;4(1):00596–12.

    CAS  Article  Google Scholar 

  24. 24.

    Wang P, Xu J, Wang Y, Cao X. An interferon-independent lncRNA promotes viral replication by modulating cellular metabolism. Science. 2017;358(6366):1051–5.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Butler AA, LeRoith D. Tissue-specific versus generalized gene targeting of the igf1 and igf1r genes and their roles in insulin-like growth factor physiology. Endocrinology. 2001;142(5):1685–8.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Stratikopoulos E, Szabolcs M, Dragatsis I, Klinakis A, Efstratiadis A. The hormonal action of IGF1 in postnatal mouse growth. Proc Natl Acad Sci U S A. 2008;105(49):19378–83.

    Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Zhang Y, Li Y, Wang Q, Zhang X, Wang D, Tang H, et al. Identification of an lncRNA-miRNA-mRNA interaction mechanism in breast cancer based on bioinformatic analysis. Mol Med Rep. 2017;16(4):5113–20.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Jin J, Chu Z, Ma P, Meng Y, Yang Y. Long non-coding RNA SPRY4-IT1 promotes proliferation and invasion by acting as a ceRNA of miR-101-3p in colorectal cancer cells. Tumour Biol. 2017;39(7):1010428317716250.

    Article  PubMed  Google Scholar 

  29. 29.

    Meier C, Hardtstock P, Joost S, Alla V, Pützer BM. p73 and IGF1R regulate emergence of aggressive Cancer stem-like features via miR-885-5p control. Cancer Res. 2016;76(2):197–205.

    Article  PubMed  Google Scholar 

  30. 30.

    Hu G, Chen D, Li X, Yang K, Wang H, Wu W. miR-133b regulates the MET proto-oncogene and inhibits the growth of colorectal cancer cells in vitro and in vivo. Cancer Biol Ther. 2010;10(2):190–7.

    Article  Google Scholar 

  31. 31.

    Chen W, Ye L, Wen D, Chen F. MiR-490-5p inhibits hepatocellular carcinoma cell proliferation, migration and invasion by directly regulating ROBO1. Pathol Oncol Res. 2017.

    Article  Google Scholar 

  32. 32.

    Wen J, Huang Y, Li H, Zhang X, Cheng P, Deng D, et al. Over-expression of miR-196b-5p is significantly associated with the progression of myelodysplastic syndrome. Int J Hematol. 2017;105(6):777–83.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Zhang X, Tang J, Zhi X, Xie K, Wang W, Li Z, et al. miR-874 functions as a tumor suppressor by inhibiting angiogenesis through STAT3/VEGF-A pathway in gastric cancer. Oncotarget. 2015;6(3):1605–17.

    Article  PubMed  PubMed Central  Google Scholar 

  34. 34.

    White NM, Khella HW, Grigull J, Adzovic S, Youssef YM, Honey RJ, et al. MiRNA profiling in metastatic renal cell carcinoma reveals a tumour-suppressor effect for miR-215. Br J Cancer. 2011;105(11):1741–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Jung K, Hu H, Saif LJ. Porcine deltacoronavirus infection: etiology, cell culture for virus isolation and propagation, molecular epidemiology and pathogenesis. Virus Research Volume. December 2016;226(2):50–9.

    CAS  Article  Google Scholar 

  36. 36.

    Mai K, Li D, Wu J, Wu Z, Cheng J, He L, et al. Complete genome sequences of two porcine Deltacoronavirus strains, CHN-GD16-03 and CHN-GD16-05, isolated in southern China, 2016. Genome Announc. 2018;6(4). Pii: e01545-17).

  37. 37.

    Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008;3(6):1101–8.

    CAS  Article  Google Scholar 

Download references


we would like to acknowledge Guangdong Wen’s Foodstuffs Group Co., Ltd. China, for providing us with piglets’ tissue samples.


This work was supported by the National Key Research and Development Program of China (No. 2016YFD0501304). The funders did not play any role in the design, conclusions or interpretation of the study.

Availability of data and materials

The data and material used and analyzed during the current study are available from the corresponding author on reasonable request. The cells used in the study were kept in the Poultry Laboratory of the College of Animal Science, South China Agricultural University and the Piglets purchased from a farm in Huanong Wen’s Co., Ltd. Do not require any administrative or ethical approval to use them.

Author information




Conceptualization, JYM; Data curation, ZHZ; Formal analysis, RTW; Funding acquisition, JYM; Investigation, XYT and YYZ; Project administration, TL and YS; Resources, JYM; Supervision, TL; Validation, YQC and XYT; Visualization, XYT; Writing – original draft, XYT. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jingyun Ma.

Ethics declarations

Ethics approval

This study was carried out in accordance with the recommendations of National Standards for Laboratory Animals of the People’s Republic of China (GB149258–2010). The protocol was approved by Animal Research Committees of South China Agricultural University. Pigs used for the study were handled in accordance with good animal practices required by the Animal Ethics Procedures and Guidelines of the People’s Republic of China.

Consent for publication

“Not applicable”.

Competing interests

The authors report no conflicts of interest. The authors themselves are responsible for the content and writing of the paper.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

The target gene of lncRNAs. (XLSX 238 kb)

Additional file 2:

Gene ontology (GO) analysis of target genes. (XLSX 10 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Tang, X., Lan, T., Wu, R. et al. Analysis of long non-coding RNAs in neonatal piglets at different stages of porcine deltacoronavirus infection. BMC Vet Res 15, 111 (2019).

Download citation


  • Long non-coding RNA (lncRNA)
  • Porcine deltacoronavirus infection
  • Functional enrichment
  • High throughput sequencing
  • Interaction network