Research article | Open | Published:
Transcriptome analysis reveals the host response to Schmallenberg virus in bovine cells and antagonistic effects of the NSs protein
BMC Genomicsvolume 16, Article number: 324 (2015)
Schmallenberg virus (SBV) is a member of the Orthobunyavirus genus (Bunyaviridae family) causing malformations and abortions in ruminants. Although, as for other members of this family/genus, the non-structural protein NSs has been shown to be an interferon antagonist, very little is known regarding the overall inhibitory effects and targets of orthobunyavirus NSs proteins on host gene expression during infection. Therefore, using RNA-seq this study describes changes to the transcriptome of primary bovine cells following infection with Schmallenberg virus (SBV) or with a mutant lacking the non-structural protein NSs (SBVdelNSs) providing a detailed comparison of the effect of NSs expression on the host cell.
The sequence reads from all samples (uninfected cells, SBV and SBVdelNSs) assembled well to the bovine host reference genome (on average 87.43% of the reads). During infection with SBVdelNSs, 649 genes were differentially expressed compared to uninfected cells (78.7% upregulated) and many of these were known antiviral and IFN-stimulated genes. On the other hand, only nine genes were differentially expressed in SBV infected cells compared to uninfected control cells, demonstrating the strong inhibitory effect of NSs on cellular gene expression. However, the majority of the genes that were expressed during SBV infection are involved in restriction of viral replication and spread indicating that SBV does not completely manage to shutdown the host antiviral response.
In this study we show the effects of SBV NSs on the transcriptome of infected cells as well as the cellular response to wild type SBV. Although NSs is very efficient in shutting down genes of the host innate response, a number of possible antiviral factors were identified. Thus the data from this study can serve as a base for more detailed mechanistic studies of SBV and other orthobunyaviruses.
Schmallenberg virus (SBV) is a negative-sense, single-stranded, segmented, RNA virus belonging to the Bunyaviridae family, within the genus Orthobunyavirus. The virus was discovered 2011 in connection with outbreaks of diarrhoea, reduced milk production and fever in cattle in Germany and the Netherlands  and has since the initial discovery rapidly spread to many European countries . Although infection is often connected with milder symptoms it can cause foetal abortions and malformations in ruminants [1,3]. SBV was the first orthobunyavirus detected in Europe, but the genus consists of over 170 viruses grouped into 18 serogroups and many of these are important pathogens. In humans, this includes for example La Crosse virus (LACV) causing encephalitis and Ngari Virus causing haermorrhagic fever . Animal pathogens, apart from SBV, includes Akabane virus (AKBV) and Cache Valley virus (CVV) that also infect ruminants and can cause abortions and congenital malformations.
The SBV genome consists of three segments (L, large; M, medium; and S, small) coding for four structural proteins (N, nucleoprotein; Gn and Gc, glycoproteins; L, RNA-dependent RNA polymerase) and two non-structural proteins (NSm and NSs). Through the production of recombinant viruses lacking NSs (LACVdelNSs, AKBVdelNSs and Bunyamwera(BUNV)delNSs), NSs has been shown to be a major virulence factor for orthobunyaviruses as these deletant viruses are attenuated and show reduced growth rate compared to wild type virus [5-7]. SBV NSs has also been shown to be able to counteract host antiviral responses. Recombinant SBV lacking ability to express NSs (SBVdelNSs) is a potent inducer of type I IFN while SBV does not induce the IFN response after infection [8,9]. Apart from interfering with the IFN pathway, SBV NSs has also been shown to induce degradation of the RPB1 subunit of RNA polymerase II and to potentially have a pro-apoptotic role . However, little is known about the overall inhibitory effects and targets of orthobunyavirus NSs proteins on gene expression in infected cells. Therefore, in this study we analysed the changes in the transcriptome of primary bovine cells following infection with either wild type SBV or SBVdelNSs (5, 6).
Results and discussion
Viral infection and sequencing output – quality and reference assembly
We carried out RNA-seq analysis in infected and mock-infected primary bovine cells, in order to use cells with an intact cell-autonomous innate immune system. Bovine primary fibroblasts were infected with either SBV, SBVdelNSs or mock-infected and RNA was extracted at 16 h p.i.. Three biological replicates were assesses for each sample and cells were also infected in parallel for immunofluorescence (Figure 1A) and Western blot (Figure 1B) in order to confirm that the bovine primary cells were infected by SBV and SBVdelNSs at 16 h. Both the SBV and the SBVdelNSs infected cells were positive for the virus with an estimated infection rate at about 30%.
The samples for transcriptome analysis were sequenced on Illumina MiSeq. On average, 26.5 million reads/sample were generated with a phred quality value of 30 or higher and 87.43% of the sequence reads assembled against the Bos taurus genome (Ensembl Btau_4.0) using TopHat2 . mRNA enrichment was carried out prior to sequencing and bunyaviruses in general (including SBV) lack a poly-A tail in their genome and mRNAs, consequently no assembly against the virus genome was performed.
Differential expression analysis
Cuffdiff2  was used to identify differentially expressed (DE) genes (genes with a +/− 2-fold change or more and with p ≤ 0.05 were considered significant). 651 DE genes were identified and most were found in the SBVdelNSs infected cells (Figure 2A). Hence, the majority of the DE genes are as a direct or indirect result of the loss of NSs. Most (78.7%) of the DE genes affected by the loss of NSs were upregulated and fold differences ranged from 12.7-fold to the 1-fold cut-off value (on a log2 scale). The fold changes for down-regulated genes (21.3%) were more subtle, ranging from 2.72-fold to the 1-fold cut-off (log2 scale). The DE sequence analysis was validated for 10 genes using Sybr green realtime PCR (Figure 2B); the fold changes for SBVdelNSs compared to the SBV infected cells were all significant (p ≤ 0.05) and corresponded to the sequencing data.
Ingenuity Pathway Analysis (IPA) (http://www.ingenuity.com/products/ipa) showed that the DE genes are, to a great extent, involved in pathways associated with host antiviral responses (Table 1), such as type I IFN-signalling and IFN-dependent gene expression, as well as pattern recognition. This is also evident when extracting the top 30 most upregulated genes (Table 2) as most of these have antiviral functions. The major molecules involved in viral RNA recognition including DDX58 (RIG-I), TRIM25, IFIH1 (MDA5), PKR, TLR3 were highly up-regulated in SBVdelNSs infected cells (Figure 3A), as are those involved in antigen presentation to CD8+ T-lymphocytes including MHC I α/β and TAP 1/2 (Figure 3B). As a consequence of activation of the viral RNA recognition pathways many interferon stimulated genes such OAS1/2, MX1 and several guanylate binding proteins, were upregulated [13-16]. Several interleukins were found among the DE transcripts, for example IL-8 and its downstream molecules (e.g. VCAM-1, ICAM-1 and Cox2) involved in inflammation, activation of angiogenesis and leucocyte infiltration/activation (Figure 3B).
In SBV-infected cells only nine DE genes (RSAD2, ISG15, OAS1, OAS2, IFIT2, MX1, GNAL, RPS3A and MDFI) were identified (eight upregulated and one down-regulated) (Figure 2A and Table 1). The majority of these genes are involved in antiviral responses. All up-regulated genes in SBV infected cells, except guanine nucleotide binding protein (G protein), alpha activating activity polypeptide, olfactory type (GNAL) and ribosomal protein S3A (RPS3A), were also upregulated following SBVdelNSs infection, although to a much higher extent, indicating that NSs may not be able to completely shutdown all genes of the host antiviral response. Viperin (RSAD2) can be induced by a number of different viruses such as sindbis virus, Japanese encephalitis virus and lassa fever virus (LASV) either dependently and independently of IFN . For example, Zapata JC et al.  showed that LASV strongly induce viperin early in infection (PBMC) while the attenuated ML29 has a weaker and delayed viperin induction suggesting that ML29 has a mechanism to interfere with the host-response signalling pathways. Antiviral activity of RSAD2 has been shown against for example BUNV, chikungunya virus, influenza virus and dengue virus through different mechanisms such as inhibition of viral replication and budding/egress [13-16,19]. MX1 has also been shown to have antiviral activity against a number of different RNA viruses such influenza and bunyaviruses . MX1 is significantly upregulated in foetuses infected with CVV and it has been proposed to be involved in the clearance of the virus . During LACV infection MX1 can bind to the nucleocapsid and inhibit viral replication . Also, ISG15 is upregulated during foetal infection with CVV . The OAS proteins performs their antiviral activity through the activation of RNAseL which leads to degradation of cellular and viral RNA  and OAS1, but not OAS2, is reported to contribute to a slight inhibition of BUNV and BUNVdelNSs viruses . Hence although very few genes were up regulated during SBV infection most of the upregulated genes have all been previously shown to have antiviral activity for other similar viruses indicating that these genes, individually or collectively, may have an important role in the antiviral response and the consequent restriction of SBV infection. However, it should be noted that we only investigated host responses at one time point, and it is possible that at an earlier or later time point during the viral replication cycle, more DE genes may be found, as shown for dengue virus for example [24,25]. As not all primary cells were infected by SBV or SBVdelNSs it is possible that weak host responses to SBV might be masked by transcripts from uninfected cells. This may also explain why the DE analysis of SBV infected cells did not show a general host transcription shutdown as described , although mRNA turnover and cell types used in this and other studies may be further reasons. Only one transcript (MDFI – MyoD family inhibitor) was down-regulated in SBV infected cells compared to uninfected cells. This gene was also down-regulated in SBVdelNSs-infected cells, indicating that this change is not due to the NSs protein alone. This is a transcription factor that negatively regulates myogenic proteins. Knockout studies of the murine homolog, inhibitory of myogenic family (l-mfa), have shown that lack of l-mfa can lead to embryonic lethality and placental defects as well as skeletal patterning defects . Although, the high viral levels of SBV seen in brain tissue are likely to cause the teratogenic effects seen during SBV infection the down-regulation of MDFI in SBV and SBVdelNSs infected cells is potentially of interest in relation to foetal abnormalities linked to SBV infection. Future studies will determine whether SBV and SBVdelINSs affect expression of MDFI in fetal brain tissue.
We have investigated the effect of SBV infection on the host transcriptome and more specifically the effects induced by lack of NSs expression. Our results show the response to orthobunyavirus infection in a relevant host cell system and that NSs is very efficient in shutting down the immune response of the host but that despite this a number of known antiviral proteins (such as viperin, MX1, OAS1/2 etc.) are still induced during SBV infection. The data from this study serves to identify possible antiviral factors and serve as a basis for more detailed mechanistic studies of SBV and other orthobunyaviruses as well as the mechanism(s) of action of their NSs protein.
Bovine fibroblast cells were isolated from cow aortas by collagenase treatment using a method adapted described earlier . Aortas were harvested from killed animals sourced from an abattoir and incubated at room temperature for 3 h in Dulbecco’s modified Eagles medium (DMEM) (Life Technologies) supplemented with 5 % foetal bovine serum (FBS), 100 U/ml penicillin 100 μg/ml streptomycin (P/S) and 2.5 μg/ml amphotericin B. After incubation, the aortas were cleaned, opened longitudinally and placed, intima layer down, into collagenase (2 mg/ml, DMEM) for 60 min at 37°C. After incubation the cells were isolated by collection of the collagenase supernatant and scraping of the aorta wall that was in contact with the collagenase. The supernatant containing the cells was centrifuged at 1000 rpm for 5 mins. The cells were then resuspended, seeded in 12-well plates and maintained at 37°C and 5% CO2 in DMEM with 20% FBS, 100 U/ml penicillin 100 μg/ml Streptomycin, 2.5 μg/ml amphotericin B. The population was initially a mixture of endothelial cells and fibroblasts but frequent passaging for approximately 10 days allowed the fibroblasts to outgrow the endothelial cells and form a pure population of fibroblasts.
Viruses and infections
Rescued wild type SBV as well as a SBV lacking NSs described in Elliot et al.  were used throughout all the experiments.
Bovine fibroblasts were seeded in 24-well plates and at a confluence of 80% the cells were infected with SBV or SBVdelNSs at a MOI of 5 for 16 h. Mock-infected cells were included in parallel as control. For each condition three individual replicates were set up. Viral infection was confirmed by immunofluorescence and western blot targeting the N-protein of the virus. For western blot γ-tubulin (Sigma) was used as an internal control.
RNA was extracted from infected and uninfected cells using a combination of Trizol (Life Technology) and RNeasy Mini kit (Qiagen). In short, the cell media was discarded from the wells and the cells were mixed and homogenised in 750 μl Trizol. Chloroform (150 μl) was added and the sample mixed prior to a short incubation at room temperature and a centrifugation step at +4°C for 15 minutes. The upper aqueous phase was transferred to a new tube and mixed with 70% EtOH and then transferred to a RNeasy Mini Spin Column. The sample was bound to the column through centrifugation and then washed once with RW1 buffer and twice with RPE-buffer before eluted with 30 μl RNase free water. The RNA was stored at −80°C until further use.
Library preparation and MiSeq sequencing
The TruSeq stranded mRNA sample preparation kit (Illumina) was used to enrich samples for mRNA and construct libraries for sequencing. Single end datasets with a read length of 150 nucleotides were generated from each sample library on an Illumina MiSeq. Each sample (SBV, SBVdelNSs and uninfected control) was set up in triplicate, and sequenced on a separate MiSeq run, with 9 runs being carried out in total. The sequences from the MiSeq runs have been deposited in the European Nucleotide Archive (ENA) and can be accessed through the study accession number PRJEB9007.
Sequence quality and assembly
The FastQC software (http://www.bioinformatics.babraham.ac.uk/projects/fastqc) was used to check the RNA-Seq reads quality in order to remove low quality reads (Q > 30). The reads that passed the quality check were assembled against the Bos taurus genome (Ensembl Btau_4.0) using TopHat2 .
Differential expression and pathway analysis
Cuffdiff2 within the Cufflink package (v2.2.1)  was used to identify differentially expressed (DE) genes and genes with a +/− 2-fold change or more and with p ≤ 0.05 were considered significant and used in the further analysis. Ingenuity Pathway Analysis (IPA) (http://www.ingenuity.com/products/ipa) was applied for the functional annotation and pathway analysis.
Realtime PCR verification
150 ng RNA were converted into cDNA using superscript (Life Technologies) with 150 ng oligo(dT) (Promega) according to the manufacturer’s instructions. Fast SYBR Green Master mix (Life Technologies) was run according to the manufacturer’s protocol, primers are available on request. GAPDH was used as a housekeeping gene and the expression fold-changes was calculated using the 2-∆∆Ct method.
La Crosse virus
Cache Valley virus
Hoffmann B, Scheuch M, Hoper D, Jungblut R, Holsteg M, Schirrmeier H, et al. Novel orthobunyavirus in Cattle, Europe, 2011. Emerg Infect Dis. 2012;18(3):469–72.
Authority EFS. “Schmallenberg virus”: analysis of the epidemiological data (November 2012). In: Supporting Publications, vol. EN-360. 2012. p. 22.
van den Brom R, Luttikholt SJ, Lievaart-Peterson K, Peperkamp NH, Mars MH, van der Poel WH, et al. Epizootic of ovine congenital malformations associated with Schmallenberg virus infection. Tijdschr Diergeneeskd. 2012;137(2):106–11.
Elliott RM. Orthobunyaviruses: recent genetic and structural insights. Nat Rev Microbiol. 2014;12(10):673–85.
Blakqori G, Delhaye S, Habjan M, Blair CD, Sanchez-Vargas I, Olson KE, et al. La Crosse bunyavirus nonstructural protein NSs serves to suppress the type I interferon system of mammalian hosts. J Virol. 2007;81(10):4991–9.
Bridgen A, Weber F, Fazakerley JK, Elliott RM. Bunyamwera bunyavirus nonstructural protein NSs is a nonessential gene product that contributes to viral pathogenesis. Proc Natl Acad Sci U S A. 2001;98(2):664–9.
Ogawa Y, Sugiura K, Kato K, Tohya Y, Akashi H. Rescue of Akabane virus (family Bunyaviridae) entirely from cloned cDNAs by using RNA polymerase I. J Gen Virol. 2007;88(Pt 12):3385–90.
Elliott RM, Blakqori G, van Knippenberg IC, Koudriakova E, Li P, McLees A, et al. Establishment of a reverse genetics system for Schmallenberg virus, a newly emerged orthobunyavirus in Europe. J Gen Virol. 2013;94(Pt 4):851–9.
Varela M, Schnettler E, Caporale M, Murgia C, Barry G, McFarlane M, et al. Schmallenberg virus pathogenesis, tropism and interaction with the innate immune system of the host. PLoS Pathog. 2013;9(1), e1003133.
Barry G, Varela M, Ratinier M, Blomstrom AL, Caporale M, Seehusen F, et al. NSs protein of Schmallenberg virus counteracts the antiviral response of the cell by inhibiting its transcriptional machinery. J Gen Virol. 2014;95(Pt 8):1640–6.
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.
Trapnell C, Hendrickson DG, Sauvageau M, Goff L, Rinn JL, Pachter L. Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat Biotechnol. 2013;31(1):46–53.
Seo JY, Yaneva R, Cresswell P. Viperin: a multifunctional, interferon-inducible protein that regulates virus replication. Cell Host Microbe. 2011;10(6):534–9.
Helbig KJ, Carr JM, Calvert JK, Wati S, Clarke JN, Eyre NS, et al. Viperin is induced following dengue virus type-2 (DENV-2) infection and has anti-viral actions requiring the C-terminal end of viperin. PLoS Negl Trop Dis. 2013;7(4), e2178.
Teng TS, Foo SS, Simamarta D, Lum FM, Teo TH, Lulla A, et al. Viperin restricts chikungunya virus replication and pathology. J Clin Invest. 2012;122(12):4447–60.
Wang X, Hinson ER, Cresswell P. The interferon-inducible protein viperin inhibits influenza virus release by perturbing lipid rafts. Cell Host Microbe. 2007;2(2):96–105.
Helbig KJ, Beard MR. The role of viperin in the innate antiviral response. J Mol Biol. 2014;426(6):1210–9.
Zapata JC, Carrion Jr R, Patterson JL, Crasta O, Zhang Y, Mani S, et al. Transcriptome analysis of human peripheral blood mononuclear cells exposed to Lassa virus and to the attenuated Mopeia/Lassa reassortant 29 (ML29), a vaccine candidate. PLoS Negl Trop Dis. 2013;7(9), e2406.
Carlton-Smith C, Elliott RM. Viperin, MTAP44, and protein kinase R contribute to the interferon-induced inhibition of Bunyamwera Orthobunyavirus replication. J Virol. 2012;86(21):11548–57.
Haller O, Staeheli P, Kochs G. Interferon-induced Mx proteins in antiviral host defense. Biochimie. 2007;89(6–7):812–8.
Rodrigues Hoffmann A, Dorniak P, Filant J, Dunlap KA, Bazer FW, de la Concha-Bermejillo A, et al. Ovine fetal immune response to Cache Valley virus infection. J Virol. 2013;87(10):5586–92.
Kochs G, Janzen C, Hohenberg H, Haller O. Antivirally active MxA protein sequesters La Crosse virus nucleocapsid protein into perinuclear complexes. Proc Natl Acad Sci U S A. 2002;99(5):3153–8.
Sadler AJ, Williams BR. Interferon-inducible antiviral effectors. Nat Rev Immunol. 2008;8(7):559–68.
Fink J, Gu F, Ling L, Tolfvenstam T, Olfat F, Chin KC, et al. Host gene expression profiling of dengue virus infection in cell lines and patients. PLoS Negl Trop Dis. 2007;1(2), e86.
Sessions OM, Tan Y, Goh KC, Liu Y, Tan P, Rozen S, et al. Host cell transcriptome profile during wild-type and attenuated dengue virus infection. PLoS Negl Trop Dis. 2013;7(3), e2107.
Kraut N, Snider L, Chen CM, Tapscott SJ, Groudine M. Requirement of the mouse I-mfa gene for placental development and skeletal patterning. EMBO J. 1998;17(21):6276–88.
Gillespie JI, Greenwell JR, Johnson C. Agonist-induced fluctuations in cytoplasmic calcium in primary cultures of bovine endothelial cells. Exp Physiol. 1991;76(5):667–76.
This work was supported by the Swedish Research Council (Grant no. 2012–6555) (A-LB) and the UK Medical Research Council (AK).
The authors declare that they have no competing interests.
ALB and JKS carried out the experimental work, apart from the library preparation and MiSeq sequencing that was performed by GW and MB. GB isolated and grew the bovine primary cells and RME contributed with the viruses used in this study. The bioinformatic analysis (quality, assambly, DE analysis) was perfomed by QG, the IPA analysis was carried out by QG and ALB. ES and MMF were involved in experimental discussions. ALB and AK concieved the study, designed the experiments and wrote the initial manuscrip draft. QG, GB, GW, JKS, MB, MMF, ES, RME and MP contributed to the draft and revising of the final manuscript. All authors read and approved the final manuscript.