- Research article
- Open Access
Transcriptome analysis of avian reovirus-mediated changes in gene expression of normal chicken fibroblast DF-1 cells
BMC Genomicsvolume 18, Article number: 911 (2017)
Avian reovirus (ARV) is an important poultry pathogen that can cause immunosuppression. In this study, RNA-Seq technology was applied to investigate the transcriptome-wide changes of DF-1 cells upon ARV infection at the middle stage.
Total RNA of ARV-infected or mock-infected samples at 10 and 18 h post infection (hpi) was extracted to build RNA-Seq datasets. Analysis of the sequencing data revealed that the expressions of numerous genes were altered, and a panel of differentially expressed genes were confirmed with RT-qPCR. At 10 hpi, 104 genes were down-regulated and 64 were up-regulated, while the expressions of 47 genes were increased and only one was down-regulated, which may play a role in retinoic acid biosynthesis, at 18 hpi in the ARV-infected cells. The similar profiles of up-regulated genes between the two groups of infected cells suggest that ARV infection activated a prolonged antiviral response of host cells. Alternative splicing analysis found no significantly changed events altered by ARV infection.
Overall, the differential expression profile presented in this study can be used to expand our understanding of the comprehensive interactions between ARV and the host cells, and may be helpful for us to reveal the pathogenic mechanism on the molecular level.
Avian reovirus (ARV) is member of the Orthoreovirus genus that has recently been classified into the Spinareovirinae subfamily, which is one of two subfamilies in the Reoviridae family . ARV is an important pathogen of birds and has been impacting poultry for nearly 60 years since it was first detected in 1957 [2, 3], and it is still prevalent in poultry until now, causing considerable economical loss in the global poultry industry [4,5,6]. Efficient and simple detecting methods may be helpful to control and prevent ARV infection . Horizontal transmission is the main route of infection, with infrequent egg transmission . Though ARV was found to be ubiquitous in poultry flocks, several strains could cause severe diseases . These pathogenic strains can cause tenosynovitis individually , and additionally usually cause mixed infections together with other pathogens, such as chicken anemia agent [10, 11]. It has been demonstrated that ARV can replicate in macrophages and cause immunosuppression [8, 12].
The pathogenicity and epizootiology of ARV have been well studied, but the pathogenesis at the molecular level is poorly understood. An excellent review on the structural and biological characteristics of ARV was published 10 years ago . Though many researchers have done brilliant work to reveal the pathogenesis of ARV infection at the molecular level in recent years, several major questions raised in the review remain unresolved. Previous studies showed that σC and P10 can induce apoptosis in different ways [13, 14], and a subsequent study correlated ARV-induced apoptosis with tissue injury . Another study demonstrated that ARV can induce autophagy to promote viral titer . Subsequent studies revealed the connection between ARV-induced autophagy and apoptosis [17, 18]. It was also demonstrated that ARV disrupts many cellular pathways, regulating protein translation, cell proliferation, and cell metabolism [19,20,21,22]. However, these results are scattered and hard to reconcile. Some studies applied proteomic analysis and microarray analysis to get a comparatively integrated data set [23, 24]. However, these methods have several disadvantages compared with RNA-Seq. RNA-Seq now provides a way to investigate virus-mediated changes on the transcriptome of host cells, with high accuracy and low background [25, 26]. Additionally, it provides information on alternative splicing events, analyzing single nucleotide polymorphism, and predicting novel transcripts [27,28,29]. In this study, we tried to build a complete expression profile of ARV-mediated changes at the transcriptional level using RNA-Seq to unveil the complex interactions between ARV and host cells.
Cell culture and virus inoculation
Chicken embryonic fibroblast cell line DF-1 (CRL-12203, ATCC) cells were cultured with high glucose (4.5 g D-Glucose/L) Dulbecco’s Modified Eagle Medium (HG-DMEM) (Basal Media, Shanghai, China) supplemented with 10% (v/v) fetal bovine serum (FBS) (Gibco, Shanghai, China) at 37 °C and 5% CO2. The ARV strain GX/2010/1, causing severe tenosynovitis and enteritis, was isolated by our lab and propagated in chicken embryo fibroblasts (CEF) cells, and was reported to trigger autophagy in host cells to promote virus production . The sequences of this strain are available in the GenBank database under the accession numbers KJ476699−KJ476708. The sixth generation of the purified virus was used in this study and the median tissue culture infective dose (TCID50) per milliliter (ml) of the virus was determined by the Reed-Muench method in CEF cells .
One day before virus inoculation, approximately 2 × 106 DF-1 cells were seeded into 75 cm2 flasks (Corning, ME, USA). When monolayer was complete (approximately 7 × 106 cells), culture medium was discarded and the cells were rinsed with phosphate buffered saline once. The purified virus was diluted to 10 multiplicity of infection (MOI) per 5 mL with HG-DMEM and applied into each flask of the ARV-infected group and an equal volume of HG-DMEM was added in the mock-infected group. After incubated at 37 °C for 1.5 h, the medium was changed to HG-DMEM supplemented with 2% FBS. Then the cells were continued to be incubated at 37 °C for 2, 10, 18 and 24 h (Fig. 1a).
Total RNA extraction and cDNA library construction
At specified hours post infection (hpi), the medium was discarded and total RNA was extracted from triplicate samples of uninfected or ARV infected groups using the Ultrapure RNA Kit (CWBIO, Beijing, China) according to manufacturer’s protocol. RNA degradation and contamination was monitored on 1% agarose gels. RNA purity was checked using the NanoPhotometer® spectrophotometer (IMPLEN, CA, USA). RNA concentration was measured using Qubit® RNA Assay Kit in Qubit® 2.0 Flurometer (Life Technologies, CA, USA). RNA integrity was assessed using the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, CA, USA).
A total amount of 3 μg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using NEBNext® Ultra™ RNA Library Prep Kit for Illumina® (NEB, USA) following the manufacturer’s recommendations, and index codes were added to attribute sequences to each sample.
Clustering and sequencing
The clustering of the index-coded samples was performed on a cBot Cluster Generation System using TruSeq PE Cluster Kit v3-cBot-HS (Illumia) according to the manufacturer’s instructions. After cluster generation, the library preparations were sequenced on an Illumina Hiseq platform and 150 bp paired-end reads were generated. Raw reads of fastq format were first processed through custom written Perl scripts. At the same time, Q20, Q30, and GC content were calculated. All of the downstream analyses were based on clean, high quality data.
Reads mapping and quantification of gene expression level
The index of the chicken reference genome (Ensembl, Galgal4, updated 11-2015) was built using Bowtie v2.2.3  and paired-end clean reads were aligned to the reference genome using TopHat v2.0.12 . And HTSeq v0.6.1 was used to count the number of reads mapped to each gene . Then, the expression level of each gene was calculated by the expected Fragments Per Kilobase of transcript per Million fragments mapped (FPKM) .
Differential expression analysis
Differential expression analysis of two groups was performed using the DESeq R package (1.18.0) . DESeq provides statistical routines for determining differential expression in digital gene expression data using a model based on a negative binomial distribution. The resulting P-values were adjusted using the Benjamini and Hochberg’s approach for controlling the false discovery rate . Genes with an adjusted P-value < 0.05 found by DESeq were considered to be differentially expressed. Additionally, KOBAS 2.0 software was used to test the statistical enrichment of differentially expressed genes (DEGs) in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways .
Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was carried out based on the basic rules of the MIQE guidelines . Briefly, 5 μg of total RNA (described above) was reverse transcribed using M-MLV reverse transcriptase (Transgen, Beijing, China) with a random hexamer primer (Genscript, Nanjing, China). The mixtures were diluted 1:10 with nuclease free water and then used as templates for qPCR. The qPCR analysis was performed using AceQ® qPCR SYBR® Green Master Mix (Vazyme, Nanjing, China) with 250 nM forward and reverse primers (Additional file 1). The reaction was carried out using LightCycler® Nano (Roche) with the following cycling conditions: an initial denaturation at 95 °C for 600 s followed by 45 cycles of 95 °C for 10 s and 60 °C for 30 s. Fold change was determined by the 2-△△Ct method .
Novel transcripts and alternative splicing prediction
The Cufflinks v2.1.1 Reference Annotation Based Transcript (RABT) assembly method was used to construct and identify both known and novel transcripts from TopHat alignment results . Alternative splicing (AS) events were classified to five major types by the software rMATS (Multivariate Analysis of Transcript Splicing) v3.2.5 . The number of AS events in each sample was estimated separately. Because the chicken genome has been recently updated , the differentially expressed novel transcripts were retrieved by the BLAST tool on National Center for Biotechnology Information (NCBI).
ARV infection of DF-1 cells and viral replication dynamics
To further study the molecular mechanism of ARV infection, DF-1 cells were infected with the virus for different time points at 10 MOI. The high dosage of virus was used to overcome the influence of uninfected cells [43, 44]. Infection and mock-infection were performed in biological triplicate for each time point and total RNA was extracted from both groups. The replication of the viral genome was determined by RT-qPCR and the fold change of M3 and S1 showed similar trends (Fig. 1c). Cytopathic effects could be seen at 18 hpi (Fig. 1b). To obtain an obviously changed transcriptome profile and minimize the influence of cell death and lysis, data at 2 hpi and 24 hpi were discarded and the remaining two groups were analyzed by RNA-Seq. One uninfected 10 hpi sample was lost, leaving a final total of 11 samples that were sequenced.
After an overall quality review, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads, and then the cDNA library was constructed, with quality assessment. After cluster generation, the library preparations were sequenced on an Illumina Hiseq platform and 150 bp paired-end reads were generated. The sequence run of each sample yielded at least 42 million clean reads and the lowest value of the reads possessing a Q-score > 20 was 95%, and the bottom line of the reads with a Q-score > 30 was 88% (Table 1). These results meet the requirements that more than 10 million reads are needed to construct a high quality eukaryotic transcriptome for discovering new genes and transcripts . Importantly, all samples had between 71.81% and 77.45% of total reads mapped to the chicken reference genome, and the percentage of uniquely mapped reads was between 70.75% and 76.26% (Table 1).
Then, the mapped data were used to predict the novel transcripts and analyze the five major types of AS events, including SE (skipped exon), A5SS (alternative 5′ splice site), A3SS (alternative 3′ splice site), MXE (mutually exclusive exon), and RI (retained intron). The predicted novel genes were further analyzed together with the known genes. No significantly changed AS events were found between each group. The mapped data was normalized by calculating the FPKM and the distribution of mean FPKM per gene was found to be uniform between the four conditions (Fig. 2a). The correlation of gene expression levels between all of the samples was investigated using the squared Pearson correlation coefficient (R2), and the minimum value was 0.979 (Fig. 2b). These results indicate that the expression levels of different genes or groups of genes are comparable, suggesting that the treatment is repeatable and has little variation. Therefore, the accuracy of the subsequent analysis of differentially expressed genes is likely to be high. Because the transcripts of the ARV genome do not have a poly-A 3′ tail , there are no reads can be used to indicate the replication of ARV.
Differentially expressed genes upon ARV infection
To further investigate the differential expression patterns in DF-1 cells between infected and mock-infected samples, the normalized gene expression level data were analyzed by DESeq. The resulting P-values were adjusted after correction for multiple testing and the DEGs were defined by having adjusted P-values (padj) < 0.05. The infected and mock-infected were compared with each other and the outline of the DEGs can be seen in Fig. 3a. Though there were significant alterations in the pairwise comparisons of different time points in the infected samples (Fig. 3b), the corresponding mock-infected samples also had big differences (Fig. 3c). This interference might result from cell culture and so should be discarded in the future studies. All of the DEGs were clustered and the results exhibited a clear time-dependent change in gene expression (Additional file 2).
The distinct effect on gene expression upon ARV infection was carefully examined. Compared with mock infected controls, 168 changes in the transcriptome, with 64 up-regulated and 104 down-regulated DEGs, were observed in response to ARV infection at 10 hpi (Fig. 3d, Additional file 3). Interestingly, only 47 up-regulated DEGs, with 31 genes in accordance with the 10 hpi group, and a single down-regulated novel gene were identified at 18 hpi (Fig. 3e, Additional file 3). These novel genes among the DEGs were retrieved by the BLAST tool. In addition, the pathway enrichment result can be seen at additional file 4. The DEGs with a fold change larger than 2 are listed in Table 2 with UniProtKB Keywords annotation , as a 2-fold threshold is commonly used to indicate biological significance .
Four up-regulated DEGs and two down-regulated DEGs were selected to be validated with RT-qPCR. The results show a similar pattern of ARV-mediated changes as was seen in the DEG analysis of RNA-Seq data (Fig. 4).
ARV is one of the major pathogens that can cause immunosuppression in poultry . Though the pathology and some molecular characteristics of ARV have been well studied [2, 8, 48], there are only a few reports that can be used to help us understand the molecular basis of ARV infection. In this study, the DEG listed in Table 2 show that DF-1 cells exerted a prolonged antiviral response upon ARV infection. The up-regulation of RSAD2 (radical S-adenosyl methionine domain containing 2), IFIT5 (interferon induced protein with tetratricopeptide repeats 5), OASL (2′-5′-oligoadenylate synthetase-like), ISG12(2) (interferon-stimulated genes) and Mx (myxovirus resistance) have been reported in the infection of infectious bursal disease virus (IBDV), which is another important pathogen similar to ARV but can cause much higher mortality and much more serious immunosuppression . EPSTI1 (epithelial-stromal interaction 1) does not have a Gene Ontology (GO) annotation. A recent report indicated that EPSTI1 plays a key role in IL-28A (interferon-λ2) mediated antiviral activity . These interferon- (IFN-) induced genes (ISGs) with high expression levels reflect the stimulation of IFNs. Even though no significant elevation of expression levels of IFN genes were identified in this study, the up-regulation of TLR3 (Toll-like receptor 3), MYD88 (Myeloid differentiation factor 88), IRF1 (IFN regulatory factor 1), and IRF3 (IFN regulatory factor 3) were found. TLR3 plays key roles in detecting virus-derived dsRNA and the TLR3 genes are polymorphic among different chicken breeds [50, 51]. In addition to TLR-induced pathways, members of the RLR family (retinoic acid inducible gene-I like receptor) constitute another TLR-independent anti-virus system. In our results, DHX58 (DEXH-box helicase 58, also known as LGP2, laboratory of genetics and physiology 2, or RLR3, RIG-I like receptor 3), IFIH1 (IFN-induced helicase C domain-containing protein 1, also known as MDA5, melanoma differentiation-associated protein 5), TRIM25 (Tripartite motif-containing protein 25), and a predicted TRIM25-like gene were found to be up-regulated at 10 hpi or at both 10 and 18 hpi time points. Chickens lack RIG-I (retinoic acid-inducible gene I), but the function of sensing viral infections can be performed by LGP2 and MDA5, which can interact with MAVS (mitochondrial antiviral signaling protein) or STING (stimulator of IFN genes) to stimulate the expression of IFNs [52,53,54].
DF-1 cells construct an antiviral environment through the expression of ISGs, including EIF2AK2 (Eukaryotic translation initiation factor 2-alpha kinase 2, also known as PKR, protein kinase RNA-activated). PKR is IFN-induced dsRNA-dependent enzymes. Active PKR can catalyze Ser-51 phosphorylation of the alpha subunit of EIF2, resulting in inhibition of protein synthesis at the initiation step of translation . However, a previous report demonstrated that σA, an ARV encoded dsRNA binding protein, can block the activation of PKR and restore translation. In that report, the inhibition of vaccinia virus replication might reflect mechanisms other than OAS and PKR to be responsible for the antiviral effects . Interestingly, another eukaryotic translation initiation factor 2-alpha kinase, EIF2AK3, was found to be down-regulated in our results. EIF2AK3, also known as PERK (PKR-like endoplasmic reticulum-resident kinase), is one of eIF2α kinases regulating gene expression in the unfolded protein response (UPR) and in amino acid starved cells . Protein synthesis can be inhibited during viral infection due to ER stress triggered by UPR, and different viruses may adapt different strategies to interfere with the activity of PERK [58, 59]. The depressed expression of PERK may reflect that ARV can impair the stress response and activate protein translation in DF-1 cells. The regulation of the host cell translation system ensures efficient replication of ARV. There is also a hypothesis that gene expression of ARV is mainly regulated at the translational level, rather than transcriptional level . The replication level of viral genome determined by RT-qPCR in our results is consistent with this hypothesis.
ARV was initially detected from the clinical case of tenosynovitis, and a direct link between the virus and disease had been conclusively demonstrated [2, 8]. Though this virus has been studied for many years, the molecular pathogenesis of the disease remains unclear. In our results, the elevated expression of a gene, WNT9a (also known as Wnt14), was observed, which might play a key role in the development of the disease. A previous report identified that Wnt14 plays a pivotal role in initiating synovial joint formation in the chick limb, but the researchers were unable to determine the specific pathway that is responsible for transducing the Wnt14 signal in joint formation . Later, studies demonstrated that the Wnt/β-catenin signaling pathway is necessary and sufficient to induce early steps of synovial joint formation . Subsequently, a precise expression pattern of various Wnts was analyzed during chick wing development . Continued expression of Wnt14 in the mature joint might be good for the maintenance of joint integrity and was presumed to play a role in the etiology of rheumatoid arthritis in humans . ARV can replicate and perhaps be persistent at hock joint of chicken . The up-regulation of Wnt14, combined with the induction of apoptosis , may be responsible for ARV-induced joint damage and more severe tendon rupture.
In conclusion, our results show that ARV infection stimulates a prolonged antiviral response in host cells and interferes with cell growth and cell death pathways. Our results also provide information that may be helpful to further investigate the pathogenesis of ARV infection. Combined with previous studies, we can begin to piece together the interactions between ARV and host cells (Fig. 5). However, the details of these interactions need to be further investigated in future studies.
Chicken embryo fibroblasts
Differentially expressed genes
Fetal bovine serum
Fragments Per Kilobase of transcript per Million fragments mapped
High glucose Dulbecco’s Modified Eagle Medium
Hours post infection
Kyoto Encyclopedia of Genes and Genomes
Multiplicity of infection
Reverse transcription-quantitative polymerase chain reaction
- TCID50 :
Median tissue culture infective dose
King AM, Adams MJ, Lefkowitz EJ, Carstens EB. Virus taxonomy: ninth report of the international committee on taxonomy of viruses. London: Academic; 2012.
Van DHL. The history of avian reovirus. Avian Dis. 2000;44(3):638–41.
Van DHL. Viral arthritis/tenosynovitis: a review. Avian Pathol. 1977;6(4):271–84.
Woźniakowski G, Samorek-Salamonowicz E, Gaweł A. Occurrence of reovirus infection in Muscovy ducks (Cairina Moschata) in south western Poland. Pol J Vet Sci. 2015;17(2):299–305.
Ayalew LE, Gupta A, Fricke J, Ahmed KA, Popowich S, Lockerbie B, Tikoo SK, Ojkic D, Gomis S. Phenotypic, genotypic and antigenic characterization of emerging avian reoviruses isolated from clinical cases of arthritis in broilers in Saskatchewan, Canada. Sci Rep. 2017;7(1):3565.
Nham EG, Pearl DL, Slavic D, Ouckama R, Ojkic D, Guerin MT. Flock-level prevalence, geographical distribution, and seasonal variation of avian reovirus among broiler flocks in Ontario. Can Vet J. 2017;58(8):828–34.
Woźniakowski G, Niczyporuk JS, Samorek-Salamonowicz E, Gaweł A. The development and evaluation of cross-priming amplification (CPA) for the detection of avian reovirus (ARV). J Appl Microbiol. 2015;118(2):528–36.
Schat KA, Skinner MA. Avian immunosuppressive diseases and immunoevasion. In: Davison F, Kaspers B, Schat KA, editors. Avian immunology. 2nd ed. London: Academic; 2014. p. 275–97.
Benavente J, Martinez-Costas J. Avian reovirus: structure and biology. Virus Res. 2007;123(2):105–19.
Engstrom BE. Blue wing disease of chickens: isolation of avian reovirus and chicken anaemia agent. Avian Pathol. 1988;17(1):23–32.
Engstrom BE, Fossum O, Luthman M. Blue wing disease of chickens: experimental infection with a Swedish isolate of chicken anaemia agent and an avian reovirus. Avian Pathol. 1988;17(1):33–50.
Mills JN, Wilcox GE. Replication of four antigenic types of avian reovirus in subpopulations of chicken leukocytes. Avian Pathol. 1993;22(2):353–61.
Shih WL, Hsu HW, Liao MH, Lee LH, Liu HJ. Avian reovirus σC protein induces apoptosis in cultured cells. Virology. 2004;321(1):65–74.
Salsman J, Top D, Boutilier J, Duncan R. Extensive syncytium formation mediated by the reovirus FAST proteins triggers apoptosis-induced membrane instability. J Virol. 2005;79(13):8090–100.
Lin HY, Chuang ST, Chen YT, Shih WL, Chang CD, Liu HJ. Avian reovirus-induced apoptosis related to tissue injury. Avian Pathol. 2007;36(2):155–9.
Meng S, Jiang K, Zhang X, Zhang M, Zhou Z, Hu M, Yang R, Sun C, Avian WY. Reovirus triggers autophagy in primary chicken fibroblast cells and Vero cells to promote virus production. Arch Virol. 2012;157(4):661–8.
Lin PY, Chang CD, Chen YC, Shih WL. RhoA/ROCK1 regulates avian Reovirus S1133-induced switch from autophagy to apoptosis. BMC Vet Res. 2015;11(1):1–12.
Duan S, Cheng J, Li C, Yu L, Zhang X. Autophagy inhibitors reduce avian-reovirus-mediated apoptosis in cultured cells and in chicken embryos. Arch Virol. 2015;160(7):1679–85.
Ji WT, Wang L, Lin RC, Huang WR, Liu HJ. Avian reovirus influences phosphorylation of several factors involved in host protein translation including eukaryotic translation elongation factor 2 (eEF2) in Vero cells. Biochem Bioph Res Co. 2009;384(3):301–5.
Liu HJ, Lin PY, Lee JW, Hsu HY, Shih WL. Retardation of cell growth by avian reovirus p17 through the activation of p53 pathway. Biochem Bioph Res Co. 2005;336(2):709–15.
Huang WR, Wang YC, Chi PI, Wang L, Wang CY, Lin CH, Liu HJ. Cell entry of avian reovirus follows a caveolin-1-mediated and dynamin-2-dependent endocytic pathway that requires activation of p38 mitogen-activated protein kinase (MAPK) and Src signaling pathways as well as microtubules and small GTPase Rab5 protein. J Biol Chem. 2011;286(35):30780–94.
Chi PI, Huang WR, Lai IH, Cheng CY, Liu HJ. The p17 nonstructural protein of avian reovirus triggers autophagy enhancing virus replication via activation of phosphatase and tensin deleted on chromosome 10 (PTEN) and AMP-activated protein kinase (AMPK), as well as dsRNA-dependent protein kinase (PKR). J Biol Chem. 2013;288(5):3571–84.
Lin PY, Liu HJ, Chang CD, Chang CI, Hsu JL, Liao MH, Lee JW, Shih WL. Avian reovirus S1133-induced DNA damage signaling and subsequent apoptosis in cultured cells and in chickens. Arch Virol. 2011;156(11):1917–29.
Chen WT, YL W, Chen T, Cheng CS, Chan HL, Chou HC, Chen YW, Yin HS. Proteomics analysis of the DF-1 chicken fibroblasts infected with avian reovirus strain S1133. PLoS One. 2014;9(3):e92154.
Kogenaru S, Yan Q, Guo Y, Wang N. RNA-seq and microarray complement each other in transcriptome profiling. BMC Genomics. 2012;13(1):629.
Liu JJ, Sturrock RN, Benton R. Transcriptome analysis of Pinus Monticola primary needles by RNA-seq provides novel insight into host resistance to Cronartium Ribicola. BMC Genomics. 2013;14(1):884.
Zhao C, Cees W, De WPJGM, Tang D, Theo VDL. RNA-Seq analysis reveals new gene models and alternative splicing in the fungal pathogen Fusarium graminearum. BMC Genomics. 2013;14(1):21.
Mcqueen CM, Whitfieldcargile CM, Konganti K, Blodgett GP, Dindot SV, Cohen ND. TRPM2 SNP genotype previously associated with susceptibility to Rhodococcus equi pneumonia in quarter horse foals displays differential gene expression identified using RNA-Seq. BMC Genomics. 2016;17(1):993.
Sheynkman GM, Johnson JE, Jagtap PD, Shortreed MR, Onsongo G, Frey BL, Griffin TJ, Smith LM. Using galaxy-P to leverage RNA-Seq for the discovery of novel protein variations. BMC Genomics. 2014;15(1):1–9.
Reed LJ, Muench H. A simple method of estimating fifty percent endpoints. Am J Epidemiol. 1938;27(3):493–9.
Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357–9.
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14(4):R36.
Anders S, Pyl PT, Huber W. HTSeq-a python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and abundance estimation from RNA-Seq reveals thousands of new transcripts and switching among isoforms. Nat Biotechnol. 2010;28(5):511–5.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc B. 1995;57(1):289–300.
Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, Kong L, Gao G, Li CY, Wei LKOBAS. 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011;39:W316–22.
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL, et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55(4):611–22.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25(4):402–8.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and cufflinks. Nat Protoc. 2012;7(3):562–78.
Shen S, Park JW, Lu ZX, Lin L, Henry MD, Wu YN, Zhou Q, Xing Y. rMATS: robust and flexible detection of differential alternative splicing from replicate RNA-Seq data. P Natl Acad Sci USA. 2014;111(51):E5593.
Genome database on National Center for Biotechnology Information. https://www.ncbi.nlm.nih.gov/genome/?term=gallus%20gallus. Accessed 30 May 2017.
Hui RK, Leung FC. Differential expression profile of chicken embryo fibroblast DF-1 cells infected with cell-adapted infectious Bursal disease virus. PLoS One. 2015;10(6):e0111771.
Sessions OM, Tan Y, Goh KC, Liu Y, Tan P, Rozen S, Ooi EE. Host cell transcriptome profile during wild-type and attenuated dengue virus infection. PLoS Negl Trop Dis. 2013;7(3):e2107.
Chen S, Wang A, Sun L, Liu F, Wang M, Jia R, Zhu D, Liu M, Yang Q, Wu Y, et al. Immune-Related Gene Expression Patterns in GPV- or H9N2- Infected Goose Spleens. Int J Mol Sci. 2016;17(12):E1990.
Universal Protein Resource databases. The European Bioinformatics Institute (EMBL-EBI), the SIB Swiss Institute of Bioinformatics and the Protein Information Resource (PIR). http://www.uniprot.org/. Accessed 30 May 2017.
Lamontagne J, Mell JC, Bouchard MJ. Transcriptome-wide analysis of hepatitis B virus-mediated changes to normal Hepatocyte gene expression. PLoS Pathog. 2016;12(2):e1005438.
Jones RC. Avian reovirus infection. Rev Sci Tech OIE. 2000;19(2):614–25.
Meng X, Yang D, Yu R, Zhu H. EPSTI1 is involved in IL-28A-mediated inhibition of HCV infection. Mediat Inflamm. 2015;2015:716315.
Karpala AJ, Lowenthal JW, Bean AG. Activation of the TLR3 pathway regulates IFNbeta production in chickens. Dev Comp Immunol. 2008;32(4):435–44.
Ruan W, An J, Wu Y. Polymorphisms of chicken TLR3 and 7 in different breeds. PLoS One. 2015;10(3):e0119967.
Liniger M, Summerfield A, Ruggli N. MDA5 can be exploited as efficacious genetic adjuvant for DNA vaccination against lethal H5N1 influenza virus infection in chickens. PLoS One. 2012;7(12):e49952.
Liniger M, Summerfield A, Zimmer G, Mccullough KC, Ruggli N. Chicken cells sense influenza a virus infection through MDA5 and CARDIF signaling involving LGP2. J Virol. 2012;86(2):705–17.
Cheng Y, Sun Y, Wang H, Yan Y, Ding C, Sun J, Chicken STING. Mediates activation of the IFN gene independently of the RIG-I gene. J Immunol. 2015;195(8):3922–36.
Zykova T, Zhu F, Bode AM, Zhang Y, Dong Z. ERK2 and RSK2 mediate phosphorylation of PKR (Thr451) and PKR directly catalyzes the phosphorylation of eIF2α at Ser51. Cancer Res. 2006;66(Suppl 8):1005.
Martínezcostas J, Gonzálezlópez C, Vakharia VN, Benavente J. Possible involvement of the double-stranded RNA-binding Core protein ςA in the resistance of avian Reovirus to interferon. J Virol. 2000;74(3):1124–31.
Harding HP, Novoa I, Zhang Y, Zeng H, Wek R, Schapira M, Ron D. Regulated translation initiation controls stress-induced gene expression in mammalian cells. Mol Cell. 2000;6(5):1099–108.
Pavio N, Romano PR, Graczyk TM, Feinstone SM, Taylor DR. Protein synthesis and endoplasmic reticulum stress can be modulated by the hepatitis C virus envelope protein E2 through the eukaryotic initiation factor 2α Kinase PERK. J Virol. 2003;77(6):3578–85.
Minakshi R, Padhan K, Rani M, Khan N, Ahmad F, Jameel S. The SARS Coronavirus 3a protein causes endoplasmic reticulum stress and induces ligand-independent downregulation of the type 1 interferon receptor. PLoS One. 2009;4(12):e8342.
Hartmann C, Tabin CJ. Wnt-14 plays a pivotal role in inducing synovial joint formation in the developing appendicular skeleton. Cell. 2001;104(3):341–51.
Guo X, Day TF, Jiang X, Garrett-Beal L, Topol L, Yang Y. Wnt/beta-catenin signaling is sufficient and necessary for synovial joint formation. Genes Dev. 2004;18(19):2404–17.
Loganathan PG, Nimmagadda S, Huang R, Scaal M, Christ B. Comparative analysis of the expression patterns of Wnts during chick limb development. Histochem Cell Biol. 2005;123(2):195–201.
Thanks to Huaisheng Wu for the assistance in performing RNA-Seq analysis.
This work was supported by grants from the National Natural Science Foundation of China (31272576), the China Agriculture Research System (CARS-41-K08), and the Priority Academic Program Development of Jiangsu Higher Education Institutions.
Availability of data and materials
The raw data of the RNA-Seq datasets was deposited in Sequence Read Archive (SRA) under the accession number SRP104704.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
RT-qPCR Primers. All RT-qPCR primers for detection of the replication of ARV and the DEGs. (DOCX 46 kb)
Cluster analysis of differentially expressed genes. Heatmap of the DEGs across all datasets based on log10 (FPKM + 1). (PDF 283 kb)
Differentially expressed genes between sample groups. Sheet 1. The DEGs between mock- or ARV-infected samples at 10 hpi. Sheet 2. The DEGs between mock- or ARV-infected samples at 18 hpi. (XLSX 74 kb)
KEGG pathway enrichment result. Sheet 1. The KEGG enrichment result of the DEGs between mock- or ARV-infected samples at 10hpi. Sheet 2. The KEGG enrichment result of the DEGs between mock- or ARV-infected samples at 18hpi. (XLSX 60 kb)