Next Article in Journal
A Proteomic Approach to Uncover Neuroprotective Mechanisms of Oleocanthal against Oxidative Stress
Previous Article in Journal
Developing Hollow-Channel Gold Nanoflowers as Trimodal Intracellular Nanoprobes
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 19
Issue 8
10.3390/ijms19082328
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Transcriptomic Analysis of Immune-Related Gene Expression in Duck Embryo Fibroblasts Following Duck Tembusu Virus Infection

by
Guanliu Yu
1,2,3,
Yun Lin
1,2,3,
Yi Tang
1,2,3,* and
Youxiang Diao
1,2,3,*
1
College of Animal Science and Technology, Shandong Agricultural University, 61 Daizong Road, Tai’an 271018, Shandong, China
2
Shandong Provincial Key Laboratory of Animal Biotechnology and Disease Control and Prevention, Shandong Agricultural University, 61 Daizong Road, Tai’an 271018, Shandong, China
3
Shandong Provincial Engineering Technology Research Center of Animal Disease Control and Prevention, Shandong Agricultural University, 61 Daizong Road, Tai’an 271018, Shandong, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2018, 19(8), 2328; https://doi.org/10.3390/ijms19082328
Submission received: 30 June 2018 / Revised: 28 July 2018 / Accepted: 31 July 2018 / Published: 8 August 2018
(This article belongs to the Section Biochemistry)
Browse Figures
Versions Notes

Abstract

:
Duck is a major waterfowl species in China, providing high-economic benefit with a population of up to 20–30 billion per year. Ducks are commonly affected by severe diseases, including egg-drop syndrome caused by duck Tembusu virus (DTMUV). The immune mechanisms against DTMUV invasion and infection remain poorly understood. In this study, duck embryo fibroblasts (DEFs) were infected with DTMUV and harvested at 12 and 24 h post-infection (hpi), and their genomes were sequenced. In total, 911 (764 upregulated and 147 downregulated genes) and 3008 (1791 upregulated and 1217 downregulated) differentially expressed genes (DEGs) were identified at 12 and 24 hpi, respectively. Kyoto Encyclopedia of Genes and Genomes enrichment analysis revealed that DEGs were considerably enriched in immune-relevant pathways, including Toll-like receptor signaling pathway, Cytosolic DNA-sensing pathway, RIG-I-like receptor signaling pathway, Chemokine signaling pathway, NOD-like receptor signaling pathway, and Hematopoietic cell lineage at both time points. The key DEGs in immune system included those of the cytokines (IFN α2, IL-6, IL-8L, IL-12B, CCR7, CCL19, and CCL20), transcription factors or signaling molecules (IRF7, NF-κB, STAT1, TMEM173, and TNFAIP3), pattern recognition receptors (RIG-I and MDA5), and antigen-presenting proteins (CD44 and CD70). This suggests DTMUV infection induces strong proinflammatory/antiviral effects with enormous production of cytokines. However, these cytokines could not protect DEFs against viral attack. Our data revealed valuable transcriptional information regarding DTMUV-infected DEFs, thereby broadening our understanding of the immune response against DTMUV infection; this information might contribute in developing strategies for controlling the prevalence of DTMUV infection.

Graphical Abstract

1. Introduction

Hemorrhagic ovarian inflammation, also termed egg-drop syndrome, is primarily caused by duck Tembusu virus (DTMUV), and has resulted in significant economic losses in the poultry industry in Southeast Asia since 2010 [1]. DTMUV, similar to other flaviviruses, is a single-stranded, positive-sense RNA virus with an approximately 11 kb genome. The virus belongs to the Ntaya antigenic group of the family Flavirviridae [2]. It infects ducks in addition to geese [3,4], chickens [5], sparrows [6], mice [7,8], and pigeons [9]. More importantly, a wide spectrum of mammalian cells exhibit cytopathic effects (CPEs) following Tembusu virus infection, such as A549, BHK21, Hela, Vero, and SH-SY5Y [10], and also including HEK293 according to our recent study. These findings imply that the virus has expanded its host range and may pose a threat to mammals’ health [11].
Better understanding of the immune responses to DTMUV invasion and infection in host cells is essential for developing strategies to control infection by the virus. Once the virus invades host cells, the first line of cellular defense (i.e., host innate immune system) involves the activation of pathogen-associated molecular patterns by pattern recognition receptors (PRRs), including toll-like receptors, nucleotide oligomerization domain (NOD)-like receptors, and retinoic acid-inducible gene I (RIG-I)-like receptors. Subsequently, interferon (IFN)-related signaling pathways are activated, resulting in the induction of IFN-stimulated genes, which can inhibit viral replication [12]. Therefore, a broadened understanding of the host immune system will be beneficial for better comprehension the molecular mechanisms underlying pathogen-host interaction.
Notably, high-throughput RNA-sequencing (RNA-Seq) technology, which has the advantages of broad genome coverage, speed, low-cost, high precision, and unbiased quantification of transcript expression, has been widely used to explore the molecular mechanisms underlying pathogen-host interactions in recent years [13,14,15,16,17,18,19]. Transcriptome analysis via RNA-Seq can reveal all RNA transcripts in tissues or cells and further elucidate the expression of genes in specific tissues or cells in different physiological states and cellular environment [20].
To date, DTMUV-related studies have focused on pathogenicity [21,22], detection methods [23,24], vaccines [10,25,26], transmission routes [27], genetic evolution and epidemiology [1,28,29], and molecular structure [30,31]. Unfortunately, the immune mechanisms of DTMUV invasion and infection in host cells remain poorly understood. Therefore, to acquire transcriptome information and to study the interactions of DTMUV with host cells, duck embryo fibroblasts (DEFs) were infected with DTMUV, harvested at 12 and 24 h post-infection (hpi), and subjected to transcriptome sequencing. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were used to examine differentially expressed genes (DEGs) identified during different stages of infection. Our findings expand our understanding of the immune mechanisms underlying DTMUV infection.

2. Results

2.1. DTMUV Infection Induces Cytopathic Effects (CPEs) and Reduces Cell Viability

As shown in Figure 1, no CPEs were observed in mock-infected cells (Figure 1a,c,e). However, in cells infected with DTMUV, CPEs (e.g., cell shrinkage, rounding, and suspension) were observed as early as 24 hpi (Figure 1d). At 48 hpi, most cells detached from culture vessel (Figure 1f).
The WST-1 cell proliferation assay was performed to examine cell viability following viral infection. As shown in Figure 2, mock-infected DEFs displayed nearly 100% viability at 12 and 24 hpi. Conversely, the viabilities of DTMUV-infected DEFs were 72.05% and 45.82% at 12 and 24 hpi (Figure 2c,f), respectively, confirming the high efficiency of viral infection.

2.2. Transcriptome Sequencing Data

As shown in Table 1, the raw reads, clean reads, clean bases, Q30, and total mapped percentage of each sample were recorded for each library. For all libraries, the percentages of clean data, Q20, and Q30 were >98%, >97%, and >94%, respectively. These results illustrated the high quality of the sequencing data and ensured their suitability for the next step of the analysis. All sequencing data have been uploaded to the Sequence Read Archive of the National Center for Biotechnology Information (accession number SRP150572).
In addition, principal component analysis was conducted on the basis of the expression level of each sample. As shown in Figure 3a, mock-infected cells (i.e., 12D and 24D) were clearly distinct from those infected with DTMUV (i.e., 12S and 24S). For DTMUV-infected samples, the correlation coefficients were proportional to the time post-infection; in other words, samples at 12 hpi did not cluster with those at 24 hpi. Hence, the infection efficiency of samples was different at 12 and 24 hpi.

2.3. Identification and Analysis of DEGs

In this study, after DEFs were infected with DTMUV, DEGs between the mock- and DTMUV-infected groups were analyzed under the threshold of p-value < 0.05 & |log2(fold change)| > 1. As shown in Figure 3b, 911 (764 upregulated and 147 downregulated genes) and 3008 (1791 upregulated and 1217 downregulated genes) DEGs were identified at 12 and 24 hpi, respectively. The details of the DEGs are listed in Table S1.

2.4. GO and KEGG Enrichment Analysis

GO is a standardized system for categorizing genes, and it is often used to explore the roles of DEGs. As shown in Figure 4, DTMUV infection was mainly related to single-organism processes, cellular processes, biological regulation, cellular components, and membrane binding at 12 and 24 hpi. Besides, the details of the GO terms are included in Table S2.
Furthermore, to further explore the functions of DEGs post DTMUV infection in DEFs, KEGG enrichment analysis was also conducted. As shown in Figure 5, DEGs were mainly enriched in the categories of signal transduction, immune system, signaling molecules and interaction, endocrine system, nervous system, and cellular community/cell growth and death at 12 and 24 hpi. The number of enriched immune system pathways was higher than those for other pathways, excluding signal transduction. Moreover, the details of the KEGG terms are included in Table S3.

2.5. Signaling Pathways in the Immune System

In this study, the DEGs were considerably enriched in immune-relevant pathways, including the Toll-like receptor signaling pathway, Cytosolic DNA-sensing pathway, RIG-I-like receptor signaling pathway, Chemokine signaling pathway, NOD-like receptor signaling pathway, and Hematopoietic cell lineage, at 12 and 24 hpi (Figure 6). As shown in Figure 6, 95.7% (22/23) of DEGs were upregulated in Toll-like receptor pathways at 12 hpi compared with 81.3% (26/32) of DEGs at 24 hpi. In Cytosolic DNA-sensing pathways, 100% (13/13) of DEGs were upregulated at 12 hpi compared with 92.3% (12/13) at 24 hpi. In Chemokine signaling pathways, 100% (19/19) of DEGs were upregulated at 12 hpi compared with 77.4% (24/31) at 24 hpi. In NOD-like receptor signaling pathways, 100% (9/9) of DEGs were upregulated at 12 hpi compared with 76.9% (10/13) at 24 hpi. Regarding DEGs involved in the Hematopoietic cell lineage, a higher percentage of DEGs were upregulated at 12 hpi than at 24 hpi (90.9% (10/11) vs. 59.9% (13/22)). Cumulatively, these results indicated that stronger immune responses were induced by DTMUV infection at 12 hpi than at 24 hpi. Detailed information about the major immune-relevant enriched genes is shown in Table S4.

2.6. Expression of Key DEGs in Immune-Related Pathways

Genes in immune-related pathways with large fold changes in expression may play important roles in the process of host resistance to virus invasion. As shown in Figure 7, the top six DEGs in the Toll-like receptor signaling pathway were IFN stimulating factor 7 (IRF7), interleukin (IL)-12B, IFN-α2, IL-8L, LOC101794331 (T-lymphocyte activation antigen CD86-like), and signal transducer and activator of transcription 1 (STAT1). For the chemokine signaling pathway, the top six DEGs were CCL19, CCR7, IL-8L, LOC101795038 (fractalkine-like), STAT1, and CCL20. For the NOD-like receptor signaling pathway, the top six DEGs were TNFAIP3, IL-8L, IL-6, NFKBIA, NFKB1, and BIRC2. For the hematopoietic cell lineage, the top five DEGs were CD7, LOC101800410 (T-cell surface glycoprotein CD1b-3-like), IL-6, ANPEP, and CD44. For the RIG-I-like receptor signaling pathway, the top six DEGs were IRF7, IL-12B, IFN-α2, transmembrane protein 173 (TMEM173/STING), IL-8L, and IFIH1. For the Cytosolic DNA-sensing signaling pathway, the top six DEGs were IRF7, IFN-α2, TMEM173, IL-6, DDX58, and NFKBIA. After ranking DEGs by expression at 12 hpi, it was apparent that immune-related DEGs mainly included cytokines (IFN-α2, IL-6, IL-8L, IL-12B, CCR7, CCL19, and CCL20), transcription factors or signaling molecules (IRF7, nuclear factor κB, STAT1, TMEM173/STING, and TNFAIP3), PRRs (RIG-I/DDX58 and MDA5/IFIH1), and antigen-presenting proteins (CD44 and CD70) (Figure 8). Details of these DEGs at 12 and 24 hpi are summarized in Table S4.

2.7. Validation of the Transcriptomic Data by Quantitative Real-Time (qRT)-PCR

To further verify the expression profiles of DEGs generated by RNA-Seq, six DEGs in the Toll-like receptor pathway were selected randomly for qRT-PCR analysis. The expression levels of the DEGs are shown in Figure 9.
Overall, the expression patterns of all six DEGs obtained via qRT-PCR were consistent with those obtained via RNA-Seq; however, the relative expression levels were not completely identical (e.g., CCL4 in Figure 8b,). Furthermore, the correlation of the expression of the six DEGs between RNA-Seq and qRT-PCR was also analyzed using Pearson’s correlation analysis. The analysis revealed a Pearson’s correlation coefficient of 0.78 (p = 0.0026), confirming that the RNA-Seq data were reliable and accurate (Figure 9g).

2.8. Determination of Cytokine Levels in Culture Supernatant

Compared with the mock-infected group, the expression of IL-6, IL-8, IL-12, and IFN-α was significantly upregulated in the culture supernatant of the DTMUV-infected group at 12 and 24 hpi (p < 0.05 or p < 0.01, Figure 10). These data further confirmed the high expression of cytokines induced by DTMUV and reflected the reliability of the data generated by RNA-Seq.

3. Discussion

DTMUV is an important pathogen that endangers waterfowl, and it has caused loss of several billion dollars in the Chinese poultry industry since 2010. Therefore, it is of great significance to deepen our understanding of the molecular mechanisms underlying pathogen-host interactions to control the occurrence and prevalence of DTMUV infection. Over the past decade, RNA-Seq technology has emerged as a revolutionary and powerful tool for revealing the molecular expression profile and specific biological processes of disease development [32]. A well-known of molecular expression profile will contribute to clarity the molecular mechanisms between pathogen and host. However, to the best of our knowledge, no report has examined the molecular expression profile of the interaction between DTMUV and DEFs, excluding a transcriptome study of the spleens of DTMUV-infected goslings [4].
In this study, to better reflect the process by which DTMUV invades host cells, DEFs were selected as an in vitro model to investigate the molecular expression profile at the cellular level. It is noteworthy that the virus usually attacks immune components in host cells during the early or middle stages of infection, and a broadened understanding of the host immune system will help us better understand the molecular mechanisms underlying pathogen–host interaction [12]. For this reason, we primarily analyzed the effect of viral invasion on the host immune system at 12 and 24 hpi.
The innate immune response is crucial for fighting viral invasion during early-stage infection. Here, several cellular signaling pathways were identified to be involved in the response to DTMUV infection at 12 and 24 hpi, including Toll-like receptor signaling pathway, Cytosolic DNA-sensing pathway, RIG-I-like receptor signaling pathway, Chemokine signaling pathway, NOD-like receptor signaling pathway, and Hematopoietic cell lineage, and so on (Figure 6). Of the above signaling pathways, pattern recognition receptors (PRRs), including toll-like receptors (TLRs), RIG-I-like receptors, and NOD-like receptors, play a central role in recognizing viral RNA and initiating the innate immune response by activating type I IFN and inflammatory cytokines [4,33].
In this study, some key PRRs in the above mentioned pathways were recognized by DTMUV, such as RIG-I, and MDA5, which triggered the host innate immune response. As previously reported, MDA5 and RIG-I can recognize long RNA fragments or dsRNAs (not triphosphorylated) and 5′-triphosphorylated or short RNA fragments, respectively [5,34]. The two PRRs were also recognized by other flaviviruses, such as Dengue virus [35,36], Zika virus [34,37], Hepatitis virus [38], Japanese encephalitis virus [35,39], West nile virus [40,41], and Avian tembusu virus [5], indicating that DTMUV may trigger the same innate immune signaling pathway as the above flaviviruses. In this study, DTMUV was sensed by PRRs in DEFs, and various immune molecules were efficiently induced. The results were somewhat similar to those obtained in a previous study, which uncovered that IFN-α, IFN-β, IFN-γ, IL-1β, IL-2, IL-6, and IL-8 were upregulated in the spleens of DTMUV-infected ducks [33]. However, these cytokines could not inhibit DTMUV replication in DEFs. This may be the reason for the occurrence of “cytokine storm”, which is a phenomenon of multiple inflammatory cytokines (e.g., TNF-α, IL-2, IL-8, IL-12, IFN-α, IFN-β, IFN-γ, and MCP-1, etc.) overproduced for resisting pathogens. The overexpressed cytokines can also cause great harm to the host when they resist pathogen invasion [33,42].
It is well-known that the transcription of innate immune molecules is depends on the activation of several transcriptional factors [5,43], such as IRF, NF-κB, and STAT. IRF7, is activated by TLR-dependent signaling pathways and it regulates type I IFN responses [44]. NF-κB is a vital transcription factor that regulates genes responsible for a variety of immune response [42]. STAT1, is one of the most important members of the STAT protein family, which plays important roles in regulating cell growth, proliferation, differentiation, and apoptosis [45]. In addition, inactivation of NF-κB or disruption of IRF7 could significantly reduce the productions of type I and type III IFN in response to avian Tembusu virus infection in 293T cells [5]. The results of the current study illustrated that the three above mentioned transcriptional factors were significantly upregulated in DEFs at 12 and 24 hpi after DTMUV infection. Meanwhile, these highly expressed transcriptional factors provide favorable conditions for the high production of ILs, IFNs, and chemotactic factors.

4. Materials and Methods

4.1. Cell Culture and Virus Infection

Duck embryo fibroblasts (DEFs) were obtained from 9-day-old specific pathogen free (SPF) duck embryos (purchased from Harbin Veterinary Research Institute, Harbin, China) according to the manufacturer’s instructions [46], cultured in Dulbecco’s modified Eagle’s medium (DMEM/F-12 1:1) (01-172-1ACS, BI, Kibbutz, Beit Haemek, Israel) at pH 7.2, supplemented with 10% fetal bovine serum (FBS) (04-001-1ACS, BI, Kibbutz, Beit Haemek, Israel), 100 U/mL penicillin and 100 mg/mL streptomycin (P1400, Solarbio Science & Technology Co., Ltd., Beijing, China), and grown in 25 cm2 flask (07-8025, BIOLOGIX, Kennebunk, ME, USA) at 37 °C in a 5% CO2 cell incubator.
When DEFs reached 80–90% confluency, they were mock-infected or infected with DTMUV (SDSM strain, GenBank Accession No. KC333867.1, which was obtained from our laboratory, the Poultry Disease Lab of Shandong Agricultural University) at a multiplicity of infection (MOI) of 3. After 1.5 h of viral adsorption at 37 °C in a 5% CO2 incubator, the inoculum was replaced with maintenance medium (i.e., DMEM/F-12 with 2% of FBS), and then harvested at 12 and 24 hpi, respectively. Three flasks of cells were collected for each treatment at each time point.
Moreover, all animal experiments were performed according to the guidelines of the Committee on Ethics of Animals of Shandong and the appropriate biosecurity, and the Animal Care and Use Committee of Shandong Agricultural University approved the protocol (No. SDAUA-2018-165).

4.2. Cell Viability Detection by WST-1 Assay

In this study, in order to evaluate the infection efficiency of DTMUV in DEFs, a WST-1 cell proliferation assay kit (Beyotime, Shanghai, China) was chosen to detect cell viability, as per the manufacturer’s instructions [47]. In total, 2 × 103 DEFs were seeded in in 96-well plates and cultured in DMEM at pH 7.2, supplemented with 2% FBS maintenance medium. Wells containing culture medium alone were used as the control. After incubation at 37 °C in a 5% CO2 incubator for 5 h, DEFs were mock-infected or infected with DTMUV at a MOI of 0.01 and maintained for 12 or 24 h. Then, 10 μL of WST-1 solution were added to each well. After 3.5 h of incubation at 37 °C in a 5% CO2 incubator, the 96-well plate was examined using a microplate reader (BioTek Instruments, Winooski, VT, USA) at a wavelength of 450 nm. Cell viability was calculated as ((OD450 of DTMUV well − OD450 of blank control well)/(OD450 of PBS well − OD450 of blank control well)) × 100%

4.3. DEFs Genomic RNA Extraction, Library Construction and Illumina Sequencing

DEFs samples were collected at 12 and 24 hpi (three replicates each) and rapidly stored at −80 °C until further use. Total genomic RNA was extracted using Trizol reagent (Invitrogen, Waltham, MA, USA), as per the manufacturer’s protocols.
The RNA samples were treated with RNase-free DNase I to remove potential genomic DNA, purified using poly-T oligo-attached magnetic beads, and segregated into 200–300-bp fragments. The integrity and concentration of RNA samples were checked using an Agilent 2100 Bioanalyzer (Agilent, Palo Alto, CA, USA), high-quality samples (i.e., A260/A280 ≥ 1.5, A260/A280 ≥ 1.0, and RNA integrity ≥ 7) were used to synthesize first-strand cDNA using reverse transcriptase with 6-base random primers. Second-strand cDNA was synthesized using the first-strand cDNA as a template, and T bases were replaced with U to generate a chain-specific library. Further, library fragments were enriched by PCR amplification, and selection based on fragment size resulted in a 300–400-bp library. Next, 150 bp pair-end readers were performed on the library using an Illumina HiSeq 2000 sequencer at Personal Biotechnology Co., Ltd. (Personal Bio, Shanghai, China).

4.4. Transcriptome Data Analysis

Clean data obtained after RNA-Seq were modified by removing adaptors, poly-N sequences, and poor-quality data using Fast-QC software (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and Cut adapt software (version 1.2.1, http://code.google.com/p/cutadapt/) [48]. Then, the clean data were assembled de novo into unigenes using Trinity software (http://trinityrnaseq.sf.net) [49], and the key parameters (e.g., clean reads, Q20, and Q30) of clean data were also calculated. The assembled unigenes were then mapped to the Peking duck reference genome (duckbase.refseq.v4.fa, http://www.duckbase.org/Download) using TopHat2 software (version 2.1.1, Baltimore, MD, USA) [50,51]. The values of total mapped reads in different samples were also computed by HTSeq software (version 0.9.0, http://htseq.readthedocs.io/en/release_0.10.0/) [52]. Finally, differences in gens expression between mock- and DTMUV-infected groups were calculated and standardized by reads per kilobase of unigene per million mapped reads (FPKM) method [53].

4.5. Differential Expression Analysis

Differential expression analysis of data collected in different groups at 12 and 24 hpi was conducted using DESeq software (http://www-huber.embl.de/users/anders/DESeq) [52]. Genes with a p-values ≤ 0.05 & |log2(fold change)| ≥ 1 in DESeq analysis were considered as differentially expressed genes (DEGs). To assess the biological function of DEGs, Gene Ontology (GO) enrichment analysis was performed using the Blast2Go software (http://www.blast2go.org/) [54,55]. In this analysis, only categories with a p-values < 0.05 were considered significantly enriched in the network. The enrichment results were mainly classified as molecular function, biological process, and cellular component.
The online Kyoto Encyclopedia of Genes and Genomes (KEGG) automatic annotation (http://www.genome.jp/kegg/tool/map_pathway2.html) was used to reveal the location of the DEGs in the different pathways with default parameters [56]. p-values < 0.05, denoted statistical significance.

4.6. Validation of RNA-Seq Data by qRT-PCR

To verify the reliability of the RNA-Seq data, nine DEGs were randomly selected for quantitative real-time PCR (qRT-PCR) using the specific primers (Table S5). First-strand cDNA was synthesized using total RNA and a PrimeScript™ RT Reagent Kit with gDNA Eraser (TaKaRa, Dalian, China). The PCR assay was conducted using TransStrat Top Green qPCR SuperMix (TransGen Biotech Co., Ltd., Beijing, China) on a 7300 Real-Time PCR instrument (Applied Bio systems, Foster City, CA, USA). Each 20-μL reaction mixture contained 0.5 μL of each primer, 1 μL of DNA template, 10 μL of 2×Trans Strat® Top Green qPCR Super Mix, and 8.0 μL of ddH2O. The PCR amplification program followed the following protocol: 94 °C for 30 s, followed by 42 cycles at 94 °C for 5 s, 60 °C for 31 s, and dissociation stage. Each sample had three biological replicates, and the duck glyceraldehyde-3-phosphate dehydrogenase (GAPDH) gene served as the endogenous reference gene. The relative expression levels of DEGs were transformed using the 2−∆∆Ct method. All data are shown as the mean ± standard deviation of five replicates. Pearson’s correlation analysis was performed to compare RNA-Seq and qRT-PCR data using the SAS 9.1 software (SAS Institute, Inc., Cary, NC, USA)

4.7. Determination of Cytokine Levels in Culture Supernatants Using Enzyme-Linked Immunoassay (ELISA)

In this study, to evaluate the expression of cytokines at the protein level following DTMUV infection, the concentrations of IFN-α, IL-6, IL-8, and IL-12 in culture supernatants were measured using ELISA kits (Mlbio, Shanghai, China), as per the manufacturer’s protocols.

5. Conclusions

In summary, our data revealed the transcriptomic profile of DTMUV-infected DEFs for the first time. Upon DTMUV infection, key PRRs (RIG-I and MDA5) in the immune system were activated, inducing increased production of cytokines (IFN-α2, IL-6, IL-8L, IL-12B, CCR7, CCL19, and CCL20), transcription factors or signaling molecules (IRF7, NF-κB, STAT1, TMEM173/STING and TNFAIP3), and antigen-presenting proteins (CD44 and CD70). These findings broaden our understanding of the immune responses to DTMUV infection, which might contribute to the development of strategies to control the prevalence of DTMUV infection.

Reference

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/19/8/2328/s1.

Author Contributions

G.Y. completed most of the experiments with the help of Y.L.; Y.T. and Y.D. provided the materials, designed the experiments, and reviewed the manuscript. All authors have read and approved the manuscript for publication.

Funding

This research was funded by the National Natural Science Foundation of China (31472199); the National Key Research and Development Program of China (2016YFD0500107); the China Agriculture Research System (CARS-42-19); Science and Technology Development Plan of Shandong Province (2014GNC111023); Funds of Shandong “Double Tops” Program.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Lei, W.; Guo, X.; Fu, S.; Feng, Y.; Tao, X.; Gao, X.; Song, J.; Yang, Z.; Zhou, H.; Liang, G. The genetic characteristics and evolution of Tembusu virus. Vet. Microbiol. 2017, 201, 32–41. [Google Scholar] [CrossRef] [PubMed]
  2. Yu, G.; Yu, X.; Yang, G.; Tang, Y.; Diao, Y. A Novel Diagnostic Method to Detect Duck Tembusu Virus: A Colloidal Gold-Based Immunochromatographic Assay. Front. Microbiol. 2018, 9, 1001. [Google Scholar] [CrossRef] [PubMed]
  3. Ti, J.; Zhang, L.; Li, Z.; Zhao, D.; Zhang, Y.; Li, F.; Diao, Y. Effect of age and inoculation route on the infection of duck Tembusu virus in Goslings. Vet. Microbiol. 2015, 181, 190–197. [Google Scholar] [CrossRef] [PubMed]
  4. He, Y.; Wang, A.; Chen, S.; Wu, Z.; Zhang, J.; Wang, M.; Jia, R.; Zhu, D.; Liu, M.; Yang, Q. Differential immune-related gene expression in the spleens of duck Tembusu virus-infected goslings. Vet. Microbiol. 2017, 212, 39–47. [Google Scholar] [CrossRef] [PubMed]
  5. Fu, G.; Chen, C.; Huang, Y.; Cheng, L.; Fu, Q.; Wan, C.; Shi, S.; Chen, H.; Liu, W. Comparative analysis of transcriptional profiles of retinoic-acid-induced gene I-like receptors and interferons in seven tissues from ducks infected with avian Tembusu virus. Arch. Virol. 2016, 161, 11–18. [Google Scholar] [CrossRef] [PubMed]
  6. Tang, Y.; Diao, Y.; Yu, C.; Gao, X.; Ju, X.; Xue, C.; Liu, X.; Ge, P.; Qu, J.; Zhang, D. Characterization of a Tembusu virus isolated from naturally infected house sparrows (Passer domesticus) in Northern China. Transbound. Emerg. Dis. 2013, 60, 152–158. [Google Scholar] [CrossRef] [PubMed]
  7. Li, S.; Li, X.; Zhang, L.; Wang, Y.; Yu, X.; Tian, K.; Su, W.; Han, B.; Su, J. Duck Tembusu virus exhibits neurovirulence in BALB/c mice. Virol. J. 2013, 10, 260. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Ti, J.; Zhang, M.; Li, Z.; Li, X.; Diao, Y. Duck Tembusu Virus Exhibits Pathogenicity to Kunming Mice by Intracerebral Inoculation. Front. Microbiol. 2016, 7, 190. [Google Scholar] [CrossRef] [PubMed]
  9. Liu, P.; Lu, H.; Li, S.; Moureau, G.; Deng, Y.Q.; Wang, Y.; Zhang, L.; Jiang, T.; de Lamballerie, X.; Qin, C.F. Genomic and antigenic characterization of the newly emerging Chinese duck egg-drop syndrome flavivirus: Genomic comparison with Tembusu and Sitiawan viruses. J. Gen. Virol. 2012, 93, 2158–2170. [Google Scholar] [CrossRef] [PubMed]
  10. Wang, H.J.; Liu, L.; Li, X.F.; Ye, Q.; Deng, Y.Q.; Qin, E.D.; Qin, C.F. In vitro and in vivo characterization of chimeric duck Tembusu virus based on Japanese encephalitis live vaccine strain SA14-14-2. J. Gen. Virol. 2016, 97, 1551–1556. [Google Scholar] [CrossRef] [PubMed]
  11. Tang, Y.; Gao, X.; Diao, Y.; Feng, Q.; Chen, H.; Liu, X.; Ge, P.; Yu, C. Tembusu virus in human, China. Transbound. Emerg. Dis. 2013, 60, 193–196. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, A.; Liu, F.; Chen, S.; Wang, M.; Jia, R.; Zhu, D.; Liu, M.; Sun, K.; Wu, Y.; Chen, X.; et al. Transcriptome Analysis and Identification of Differentially Expressed Transcripts of Immune-Related Genes in Spleen of Gosling and Adult Goose. Int. J. Mol. Sci. 2015, 16, 22904–22926. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Cong, F.L.; Liu, X.; Han, Z.X.; Shao, Y.H.; Kong, X.G.; Liu, S.W. Transcriptome analysis of chicken kidney tissues following coronavirus avian infectious bronchitis virus infection. BMC Genom. 2013, 14, 743. [Google Scholar] [CrossRef] [PubMed]
  14. Jagya, N.; Varma, S.P.K.; Thakral, D.; Joshi, P.; Duragapal, H.; Panda, S.K. RNA-Seq Based Transcriptome Analysis of Hepatitis E Virus (HEV) and Hepatitis B Virus (HBV) Replicon Transfected Huh-7 Cells. PLoS ONE 2014, 9, e87835. [Google Scholar] [CrossRef] [PubMed]
  15. McLoughlin, K.E.; Nalpas, N.C.; Rue-Albrecht, K.; Browne, J.A.; Magee, D.A.; Killick, K.E.; Park, S.D.; Hokamp, K.; Meade, K.G.; O’Farrelly, C.; et al. RNA-seq Transcriptional Profiling of Peripheral Blood Leukocytes from Cattle Infected with Mycobacterium bovis. Front. Immunol. 2014, 5, 396. [Google Scholar] [CrossRef] [PubMed]
  16. Wang, H.K.; Li, S.C.; Qi, K.Z.; Xue, T.; Tu, J.; Zhou, X.H.; Liu, H.M. Application of RNA-Seq technology for screening differentially expressed genes in chicken small intestine cells affected by avian pathogenic Escherichia coli. Chin. J. Prev. Vet. Med. 2016, 38, 111–115. [Google Scholar]
  17. Menicucci, A.R.; Versteeg, K.; Woolsey, C.; Mire, C.E.; Geisbert, J.B.; Cross, R.W.; Agans, K.N.; Jankeel, A.; Geisbert, T.W.; Messaoudi, I. Transcriptome Analysis of Circulating Immune Cell Subsets Highlight the Role of Monocytes in Zaire Ebola Virus Makona Pathogenesis. Front. Immunol. 2017, 8, 1372. [Google Scholar] [CrossRef] [PubMed]
  18. Preusse, M.; Schughart, K.; Pessler, F. Host Genetic Background Strongly Affects Pulmonary microRNA Expression before and during Influenza A Virus Infection. Front. Immunol. 2017, 8. [Google Scholar] [CrossRef] [PubMed]
  19. Chen, G.; He, L.; Luo, L.; Huang, R.; Liao, L.; Li, Y.; Zhu, Z.; Wang, Y. Transcriptomics Sequencing Provides Insights into Understanding the Mechanism of Grass Carp Reovirus Infection. Int. J. Mol. Sci. 2018, 19, 488. [Google Scholar] [CrossRef] [PubMed]
  20. Yin, F.; Gao, Q.; Tang, B.; Sun, P.; Han, K.; Huang, W. Transcriptome and analysis on the complement and coagulation cascades pathway of large yellow croaker (Larimichthys crocea) to ciliate ectoparasite Cryptocaryon irritans infection. Fish Shellfish Immunol. 2016, 50, 127–141. [Google Scholar] [CrossRef] [PubMed]
  21. Li, N.; Lv, C.; Yue, R.; Shi, Y.; Wei, L.; Chai, T.; Liu, S. Effect of age on the pathogenesis of duck tembusu virus in Cherry Valley ducks. Front. Microbiol. 2015, 6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Shen, H.Q.; Lin, W.C.; Wang, Z.X.; Zhang, K.; Yan, Z.Q.; Zhou, Q.F.; Qin, J.P.; Xie, Q.M.; Bi, Y.Z.; Chen, F. Pathogenicity and genetic characterization of a duck Tembusu virus associated with egg-dropping in Muscovy ducks. Virus Res. 2016, 223, 52–56. [Google Scholar] [CrossRef] [PubMed]
  23. Tang, Y.; Yeh, Y.T.; Chen, H.; Yu, C.; Gao, X.; Diao, Y. Comparison of four molecular assays for the detection of Tembusu virus. Avian Pathol. 2015, 44, 379–385. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Wu, X.; Shi, Y.; Yan, D.; Li, X.; Yan, P.; Gao, X.; Zhang, Y.; Yu, L.; Ren, C.; Li, G.; et al. Development of a PCR-Based Reverse Genetics System for an Attenuated Duck Tembusu Virus Strain. PLoS ONE 2016, 11, e0156579. [Google Scholar] [CrossRef] [PubMed]
  25. Ma, T.; Liu, Y.; Cheng, J.; Liu, Y.; Fan, W.; Cheng, Z.; Niu, X.; Liu, J. Liposomes containing recombinant E protein vaccine against duck Tembusu virus in ducks. Vaccine 2016, 34, 2157–2163. [Google Scholar] [CrossRef] [PubMed]
  26. Zhang, L.; Li, Z.; Zhang, Q.; Sun, M.; Li, S.; Su, W.; Hu, X.; He, W.; Su, J. Efficacy assessment of an inactivated Tembusu virus vaccine candidate in ducks. Res. Vet. Sci. 2017, 110, 72–78. [Google Scholar] [CrossRef] [PubMed]
  27. Zhang, Y.; Li, X.; Chen, H.; Ti, J.; Yang, G.; Zhang, L.; Lu, Y.; Diao, Y. Evidence of possible vertical transmission of Tembusu virus in ducks. Vet. Microbiol. 2015, 179, 149–154. [Google Scholar] [CrossRef] [PubMed]
  28. Dai, L.; Li, Z.; Tao, P. Evolutionary analysis of Tembusu virus: Evidence for the emergence of a dominant genotype. Infect. Genet. Evol. 2015, 32, 124–129. [Google Scholar] [CrossRef] [PubMed]
  29. Zhou, X.; Zhang, T.; Song, D.; Huang, T.; Peng, Q.; Chen, Y.; Li, A.; Zhang, F.; Wu, Q.; Ye, Y.; Tang, Y. Whole-Genome Sequence of Duck Tembusu Virus Strain DTMUV/CH/2014, Isolated in China. Genome Announc. 2016, 4. [Google Scholar] [CrossRef] [PubMed]
  30. Han, K.; Zhao, D.; Liu, Y.; Huang, X.; Yang, J.; Liu, Q.; An, F.; Li, Y. Design and evaluation of a polytope construct with multiple B and T epitopes against Tembusu virus infection in ducks. Res. Vet. Sci. 2016, 104, 174–180. [Google Scholar] [CrossRef] [PubMed]
  31. Zhao, D.; Huang, X.; Liu, Y.; Han, K.; Zhang, J.; Yang, J.; Xie, X.; Li, Y. Domain I and II from newly emerging goose tembusu virus envelope protein functions as a dominant-negative inhibitor of virus infectivity. Res. Vet. Sci. 2015, 98, 121–126. [Google Scholar] [CrossRef] [PubMed]
  32. Wang, Z.; Gerstein, M.; Snyder, M. RNA-Seq: A revolutionary tool for transcriptomics. Nat. Rev. 2009, 10, 57–63. [Google Scholar] [CrossRef] [PubMed]
  33. Li, N.; Wang, Y.; Li, R.; Liu, J.; Zhang, J.; Cai, Y.; Liu, S.; Chai, T.; Wei, L. Immune responses of ducks infected with duck Tembusu virus. Front. Microbiol. 2015, 6, 425. [Google Scholar] [CrossRef] [PubMed]
  34. Chazal, M.; Beauclair, G.; Gracias, S.; Najburg, V.; Simon-Loriere, E.; Tangy, F.; Komarova, A.V.; Jouvenet, N. RIG-I Recognizes the 5’ Region of Dengue and Zika Virus Genomes. Cell Rep. 2018, 24, 320–328. [Google Scholar] [CrossRef] [PubMed]
  35. Loo, Y.M.; Fornek, J.; Crochet, N.; Bajwa, G.; Perwitasari, O.; Martinez-Sobrido, L.; Akira, S.; Gill, M.A.; Garcia-Sastre, A.; Katze, M.G.; et al. Distinct RIG-I and MDA5 signaling by RNA viruses in innate immunity. J. Virol. 2008, 82, 335–345. [Google Scholar] [CrossRef] [PubMed]
  36. Urcuqui-Inchima, S.; Cabrera, J.; Haenni, A.L. Interplay between dengue virus and Toll-like receptors, RIG-I/MDA5 and microRNAs: Implications for pathogenesis. Antivir. Res. 2017, 147, 47–57. [Google Scholar] [CrossRef] [PubMed]
  37. Khandia, R.; Munjal, A.; Dhama, K.; Karthik, K.; Tiwari, R.; Malik, Y.S.; Singh, R.K.; Chaicumpa, W. Modulation of Dengue/Zika Virus Pathogenicity by Antibody-Dependent Enhancement and Strategies to Protect Against Enhancement in Zika Virus Infection. Front. Immunol. 2018, 9, 597. [Google Scholar] [CrossRef] [PubMed]
  38. Zhu, C.; Xiao, F.; Hong, J.; Wang, K.; Liu, X.; Cai, D.; Fusco, D.N.; Zhao, L.; Jeong, S.W.; Brisac, C.; et al. EFTUD2 is a Novel Innate Immune Regulator Restricting Hepatitis C Virus Infection through the RIG-I/MDA5 Pathway. J. Virol. 2015, 89, 6608–6618. [Google Scholar] [CrossRef] [PubMed]
  39. Kato, H.; Takeuchi, O.; Sato, S.; Yoneyama, M.; Yamamto, M.; Matsui, K.; Uematsu, S.; Jung, A. Differential roles of MDA5 and RIG-I helicases in the recognition of RNA viruses. Nature 2006, 441, 101–105. [Google Scholar] [CrossRef] [PubMed]
  40. Fredericksen, B.L.; Keller, B.C.; Fornek, J.; Katze, M.G.; Gale, M.J. Establishment and maintenance of the innate antiviral response to West Nile Virus involves both RIG-I and MDA5 signaling through IPS-1. J. Virol. 2008, 82, 609–616. [Google Scholar] [CrossRef] [PubMed]
  41. Zhang, H.L.; Ye, H.Q.; Liu, S.Q.; Deng, C.L.; Li, X.D.; Shi, P.Y.; Zhang, B. West Nile Virus NS1 Antagonizes Interferon Beta Production by Targeting RIG-I and MDA5. J. Virol. 2017, 91. [Google Scholar] [CrossRef] [PubMed]
  42. Livolsi, A.; Busuttil, V.; Imbert, V.; Abraham, R.T.; Peyron, J.F. Tyrosine phosphorylation-dependent activation of NF-κB Requirement for p56 LCK and ZAP-70 protein tyrosine kinases. Eur. J. Biochem. 2001, 268, 1508–1515. [Google Scholar] [CrossRef] [PubMed]
  43. Xing, J.; Ni, L.; Wang, S.; Wang, K.; Lin, R.; Zheng, C. Herpes simplex virus 1-encoded tegument protein VP16 abrogates the production of beta interferon (IFN) by inhibiting NF-kappaB activation and blocking IFN regulatory factor 3 to recruit its coactivator CBP. J. Virol. 2013, 87, 9788–9801. [Google Scholar] [CrossRef] [PubMed]
  44. Shahsavandi, S.; Ebrahimi, M.M.; Mohammadi, A.; Zarrin Lebas, N. Impact of chicken-origin cells on adaptation of a low pathogenic influenza virus. Cytotechnology 2013, 65, 419–424. [Google Scholar] [CrossRef] [PubMed]
  45. Chen, L.J.; Liu, G.J. Progress and Prospect of STAT1 in Tumorigenesis. Chin. J. Cell Biol. 2017, 39, 1612–1628. [Google Scholar]
  46. Adams, S.; Xing, Z.; Li, J.; Mendoza, K.; Perez, D.; Reed, K.; Cardona, C. The effect of avian influenza virus NS1 allele on virus replication and innate gene expression in avian cells. Mol. Immunol. 2013, 56, 358–368. [Google Scholar] [CrossRef] [PubMed]
  47. Li, Y.; Xie, P.; Sun, M.; Xiang, B.; Kang, Y.; Gao, P.; Zhu, W.; Ning, Z.; Ren, T. S1PR1 expression correlates with inflammatory responses to Newcastle disease virus infection. Infect. Genet Evol. 2016, 37, 37–42. [Google Scholar] [CrossRef] [PubMed]
  48. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011, 17, 10–11. [Google Scholar] [CrossRef]
  49. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L. Trinity: Reconstructing a full-length transcriptome without a reference genome from RNA-Seq data. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef] [PubMed]
  50. Trapnell, C.; Pachter, L.; Salzberg, S.L. TopHat: Discovering splice junctions with RNA-Seq. Bioinformatics 2009, 25, 1105–1111. [Google Scholar] [CrossRef] [PubMed]
  51. Kim, D.; Pertea, G.; Trapnell, C.; Pimentel, H.; Kelley, R.; Salzberg, S.L. TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14, R36. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Anders, S.; Huber, W. Differential expression analysis for sequence count data. Genome Biol. 2010, 11, R106. [Google Scholar] [CrossRef] [PubMed]
  53. Mortazavi, A.; Williams, B.A.; McCue, K.; Schaeffer, L.; Wold, B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods 2008, 50, 621–628. [Google Scholar] [CrossRef] [PubMed]
  54. Conesa, A.; Gotz, S. Blast2GO: A comprehensive suite for functional analysis in plant genomics. Int. J. Plant Genom. 2008, 2008, 619832. [Google Scholar] [CrossRef] [PubMed]
  55. Bindea, G.; Mlecnik, B.; Hackl, H.; Charoentong, P.; Tosolini, M.; Kirilovsky, A.; Fridman, W.H.; Pages, F.; Trajanoski, Z.; Galon, J. ClueGO: A Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics 2009, 25, 1091–1093. [Google Scholar] [CrossRef] [PubMed]
  56. Kanehisa, M.; Araki, M.; Goto, S.; Hattori, M.; Hirakawa, M.; Itoh, M.; Katayama, T.; Kawashima, S.; Okuda, S.; Tokimatsu, T.; et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008, 36, D480–D484. [Google Scholar] [CrossRef] [PubMed]
Ijms 19 02328 g001 550
Figure 1. Cytopathic effects following infection of duck embryo fibroblasts with duck Tembusu virus at 12, 24, and 48 h post-infection. (a,c,e) no CPEs were observed in mock-infected cells. (b,d,f) the CPEs (e.g., cell shrinkage, rounding, and suspension) were found in cells infected with DTMUV. The blue scale bar represents 200 μm.
Figure 1. Cytopathic effects following infection of duck embryo fibroblasts with duck Tembusu virus at 12, 24, and 48 h post-infection. (a,c,e) no CPEs were observed in mock-infected cells. (b,d,f) the CPEs (e.g., cell shrinkage, rounding, and suspension) were found in cells infected with DTMUV. The blue scale bar represents 200 μm.
Ijms 19 02328 g001
Ijms 19 02328 g002 550
Figure 2. Viability of duck embryo fibroblasts at 12 (a–c) and 24 h (d–f) after infection with duck Tembusu virus (DTMUV). CPEs (cell shrinkage and rounding, red arrow) were observed in Figure 2b,e. Cell viability is shown as the mean ± SD of eight replicates.** p < 0.01.
Figure 2. Viability of duck embryo fibroblasts at 12 (a–c) and 24 h (d–f) after infection with duck Tembusu virus (DTMUV). CPEs (cell shrinkage and rounding, red arrow) were observed in Figure 2b,e. Cell viability is shown as the mean ± SD of eight replicates.** p < 0.01.
Ijms 19 02328 g002
Ijms 19 02328 g003 550
Figure 3. Principal component of samples (a) and visibility graph of differentially expressed genes (b) at 12 and 24 h post-infection (hpi). Note: (a): PC1 shows the differences among duck Tembusu virus (DTMUV)-infected samples; PC2 indicates differences between mock- and DTMUV-infected samples. 12D, mock-infected DEFs at 12 hpi; 24D, mock-infected DEFs at 24 hpi; 12S, DTMUV-infected DEFs at 12 hpi; 24S, DTMUV-infected DEFs at 24 hpi. (b): Upregulated differentially expressed genes (DEGs) are shown in red, whereas downregulated DEGs are shown in green. The green curve in the circle represents the common DEGs at 12 and 24 hpi.
Figure 3. Principal component of samples (a) and visibility graph of differentially expressed genes (b) at 12 and 24 h post-infection (hpi). Note: (a): PC1 shows the differences among duck Tembusu virus (DTMUV)-infected samples; PC2 indicates differences between mock- and DTMUV-infected samples. 12D, mock-infected DEFs at 12 hpi; 24D, mock-infected DEFs at 24 hpi; 12S, DTMUV-infected DEFs at 12 hpi; 24S, DTMUV-infected DEFs at 24 hpi. (b): Upregulated differentially expressed genes (DEGs) are shown in red, whereas downregulated DEGs are shown in green. The green curve in the circle represents the common DEGs at 12 and 24 hpi.
Ijms 19 02328 g003
Ijms 19 02328 g004 550
Figure 4. Gene Ontology enrichment of differentially expressed genes (a): 12 h post-infection (hpi); (b): 24 h post-infection 24 hpi.
Figure 4. Gene Ontology enrichment of differentially expressed genes (a): 12 h post-infection (hpi); (b): 24 h post-infection 24 hpi.
Ijms 19 02328 g004
Ijms 19 02328 g005 550
Figure 5. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis at 12 (a) and 24 h post-infection (b). The top six enriched KEGG terms at both the time points are marked in black bold font.
Figure 5. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis at 12 (a) and 24 h post-infection (b). The top six enriched KEGG terms at both the time points are marked in black bold font.
Ijms 19 02328 g005
Ijms 19 02328 g006 550
Figure 6. The immune system pathways enriched during duck Tembusu virus infections at 12 h (a) and 24 h (b) post-infection (hpi) (p < 0.05).
Figure 6. The immune system pathways enriched during duck Tembusu virus infections at 12 h (a) and 24 h (b) post-infection (hpi) (p < 0.05).
Ijms 19 02328 g006
Ijms 19 02328 g007 550
Figure 7. Differentially expressed genes (DEGs) in six common immune system-related signaling pathways at 12 and 24 h post-infection (hpi). (a) DEGs in Toll-like receptor signaling pathway at 12 and 24 hpi. (b) DEGs in Chemokine signaling pathway at 12 and 24 hpi. (c) DEGs in Nod-like receptor signaling pathway at 12 and 24 hpi. (d) DEGs in Hematopoietic cell lineage at 12 and 24 hpi. (e) DEGs in RIG-I-like receptor signaling pathway at 12 and 24 hpi. (f) DEGs in Cytosolic DNA-sensing signaling pathway at 12 and 24 hpi. The six DEGs with the largest fold changes in expression in each pathway are marked by black dots.
Figure 7. Differentially expressed genes (DEGs) in six common immune system-related signaling pathways at 12 and 24 h post-infection (hpi). (a) DEGs in Toll-like receptor signaling pathway at 12 and 24 hpi. (b) DEGs in Chemokine signaling pathway at 12 and 24 hpi. (c) DEGs in Nod-like receptor signaling pathway at 12 and 24 hpi. (d) DEGs in Hematopoietic cell lineage at 12 and 24 hpi. (e) DEGs in RIG-I-like receptor signaling pathway at 12 and 24 hpi. (f) DEGs in Cytosolic DNA-sensing signaling pathway at 12 and 24 hpi. The six DEGs with the largest fold changes in expression in each pathway are marked by black dots.
Ijms 19 02328 g007
Ijms 19 02328 g008 550
Figure 8. The top six differentially expressed genes in the Toll-like receptor signaling pathway, Chemokine signaling pathway, Nod-like receptor signaling pathway, Hematopoietic cell lineage, RIG-I-like receptor signaling pathway, and Cytosolic DNA-sensing signaling pathway. The colors in the grid represent the expression fold change as indicated by the color scale. Note: LOC101794331 (T-lymphocyte activation antigen CD86-like); LOC101795038 (fractalkine-like); LOC101800410 (T-cell surface glycoprotein CD1b-3-like).
Figure 8. The top six differentially expressed genes in the Toll-like receptor signaling pathway, Chemokine signaling pathway, Nod-like receptor signaling pathway, Hematopoietic cell lineage, RIG-I-like receptor signaling pathway, and Cytosolic DNA-sensing signaling pathway. The colors in the grid represent the expression fold change as indicated by the color scale. Note: LOC101794331 (T-lymphocyte activation antigen CD86-like); LOC101795038 (fractalkine-like); LOC101800410 (T-cell surface glycoprotein CD1b-3-like).
Ijms 19 02328 g008
Ijms 19 02328 g009 550
Figure 9. Confirmation of the transcriptome sequencing data by quantitative real-time (qRT)-PCR. Six differentially expressed genes (DEGs) involved in immune-related pathways were selected for qRT-PCR, and expression was estimated using the 2−∆∆CT method. The qRT-PCR data are presented as the mean ± standard deviation (n = 5) (a–f). Pearson’s correlation analyses of the expression of above six DEGs between Seq-RNA and qRT-PCR (g).
Figure 9. Confirmation of the transcriptome sequencing data by quantitative real-time (qRT)-PCR. Six differentially expressed genes (DEGs) involved in immune-related pathways were selected for qRT-PCR, and expression was estimated using the 2−∆∆CT method. The qRT-PCR data are presented as the mean ± standard deviation (n = 5) (a–f). Pearson’s correlation analyses of the expression of above six DEGs between Seq-RNA and qRT-PCR (g).
Ijms 19 02328 g009
Ijms 19 02328 g010 550
Figure 10. Concentrations of IL-6 (a); IL-8 (b); IL-12 (c); and IFN-α (d) in the culture supernatant of the mock- and duck Tembusu virus (DTMUV)-infected groups. The data are presented as the mean ± standard deviation (n = 6). * p < 0.05, ** p < 0.01.
Figure 10. Concentrations of IL-6 (a); IL-8 (b); IL-12 (c); and IFN-α (d) in the culture supernatant of the mock- and duck Tembusu virus (DTMUV)-infected groups. The data are presented as the mean ± standard deviation (n = 6). * p < 0.05, ** p < 0.01.
Ijms 19 02328 g010
Table
Table 1. Sequencing data generated by an Illumina HiSeqTM 2000 sequencer.
Table 1. Sequencing data generated by an Illumina HiSeqTM 2000 sequencer.
SampleRaw ReadsClean ReadsClean Data (%)Q20 (%)Q30 (%)Total Mapped (%)
12hD182,525,57681,938,15899.2897.7094.2976.51
12hD282,444,20481,830,30299.2597.6694.2477.84
12hD382,505,91281,871,54899.2397.6294.1576.91
12hS181,086,67680,390,53499.1497.7094.2869.68
12hS281,606,42281,000,12099.2597.7794.4161.35
12hS384,368,22283,680,92499.1897.7294.3065.83
24hD184,046,38283,435,38299.2897.6394.1573.81
24hD285,114,91284,487,66699.2697.5794.0371.89
24hD385,158,52084,518,93699.2497.5694.0373.36
24hS184,904,01683,892,83298.8097.7094.2446.64
24hS285,411,01084,413,20898.8397.8294.4645.81
24hS384,786,93083,734,22698.7597.6894.1745.10

Share and Cite

MDPI and ACS Style

Yu, G.; Lin, Y.; Tang, Y.; Diao, Y. Comparative Transcriptomic Analysis of Immune-Related Gene Expression in Duck Embryo Fibroblasts Following Duck Tembusu Virus Infection. Int. J. Mol. Sci. 2018, 19, 2328. https://doi.org/10.3390/ijms19082328

AMA Style

Yu G, Lin Y, Tang Y, Diao Y. Comparative Transcriptomic Analysis of Immune-Related Gene Expression in Duck Embryo Fibroblasts Following Duck Tembusu Virus Infection. International Journal of Molecular Sciences. 2018; 19(8):2328. https://doi.org/10.3390/ijms19082328

Chicago/Turabian Style

Yu, Guanliu, Yun Lin, Yi Tang, and Youxiang Diao. 2018. "Comparative Transcriptomic Analysis of Immune-Related Gene Expression in Duck Embryo Fibroblasts Following Duck Tembusu Virus Infection" International Journal of Molecular Sciences 19, no. 8: 2328. https://doi.org/10.3390/ijms19082328

APA Style

Yu, G., Lin, Y., Tang, Y., & Diao, Y. (2018). Comparative Transcriptomic Analysis of Immune-Related Gene Expression in Duck Embryo Fibroblasts Following Duck Tembusu Virus Infection. International Journal of Molecular Sciences, 19(8), 2328. https://doi.org/10.3390/ijms19082328

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop