Transcriptome analysis reveals the host response to Schmallenberg virus in bovine cells and antagonistic effects of the NSs protein

Transcriptome profiling highlights regulated biological processes and type III interferon antiviral responses upon Crimean-Congo hemorrhagic fever virus infection

Crimean-Congo hemorrhagic fever virus (CCHFV) is a biosafety level-4 (BSL-4) pathogen that causes Crimean-Congo hemorrhagic fever (CCHF) characterized by hemorrhagic manifestation, multiple organ failure and high mortality rate, posing great threat to public health. Despite the recently increasing research efforts on CCHFV, host cell responses associated with CCHFV infection remain to be further characterized. Here, to better understand the cellular response to CCHFV infection, we performed a transcriptomic analysis in human kidney HEK293 ​cells by high-throughput RNA sequencing (RNA-seq) technology. In total, 496 differentially expressed genes (DEGs), including 361 up-regulated and 135 down-regulated genes, were identified in CCHFV-infected cells. These regulated genes were mainly involved in host processes including defense response to virus, response to stress, regulation of viral process, immune response, metabolism, stimulus, apoptosis and protein catabolic process. Therein, a significant up-regulation of type III interferon (IFN) signaling pathway as well as endoplasmic reticulum (ER) stress response was especially remarkable. Subsequently, representative DEGs from these processes were well validated by RT-qPCR, confirming the RNA-seq results and the typical regulation of IFN responses and ER stress by CCHFV. Furthermore, we demonstrate that not only type I but also type III IFNs (even at low dosages) have substantial anti-CCHFV activities. Collectively, the data may provide new and comprehensive insights into the virus-host interactions and particularly highlights the potential role of type III IFNs in restricting CCHFV, which may help inform further mechanistic delineation of the viral infection and development of anti-CCHFV strategies.

Functional Variants Associated With CMPK2 and in ASB16 Influence Bovine Digital Dermatitis

Bovine digital dermatitis (BDD) is an infectious disease of the hoof in cattle with multifactorial etiology and a polygenic influence on susceptibility. With our study, we identified genomic regions with the impact on occurrence and development of BDD. We used 5,040 genotyped animals with phenotype information based on the M-stage system for genome-wide association. Significant associations for single-nucleotide polymorphisms were found near genes CMPK2 (chromosome 11) and ASB16 (chromosome 19) both being implicated in immunological processes. A sequence analysis of the chromosomal regions revealed rs208894039 and rs109521151 polymorphisms as having significant influence on susceptibility to the disease. Specific genotypes were significantly more likely to be affected by BDD and developed chronic lesions. Our study provides an insight into the genomic background for a genetic predisposition related to the pathogenesis of BDD. Results might be implemented in cattle-breeding programs and could pave the way for the establishment of a BDD prescreening test.

Host Cell Restriction Factors of Bunyaviruses and Viral Countermeasures

The Bunyavirales order comprises more than 500 viruses (generally defined as bunyaviruses) classified into 12 families. Some of these are highly pathogenic viruses infecting different hosts, including humans, mammals, reptiles, arthropods, birds, and/or plants. Host cell sensing of infection activates the innate immune system that aims at inhibiting viral replication and propagation. Upon recognition of pathogen-associated molecular patterns (PAMPs) by cellular pattern recognition receptors (PRRs), numerous signaling cascades are activated, leading to the production of interferons (IFNs). IFNs act in an autocrine and paracrine manner to establish an antiviral state by inducing the expression of hundreds of IFN-stimulated genes (ISGs). Some of these ISGs are known to restrict bunyavirus infection. Along with other constitutively expressed host cellular factors with antiviral activity, these proteins (hereafter referred to as “restriction factors”) target different steps of the viral cycle, including viral entry, genome transcription and replication, and virion egress. In reaction to this, bunyaviruses have developed strategies to circumvent this antiviral response, by avoiding cellular recognition of PAMPs, inhibiting IFN production or interfering with the IFN-mediated response. Herein, we review the current knowledge on host cellular factors that were shown to restrict infections by bunyaviruses. Moreover, we focus on the strategies developed by bunyaviruses in order to escape the antiviral state developed by the infected cells.

Interactions of Viral Proteins from Pathogenic and Low or Non-Pathogenic Orthohantaviruses with Human Type I Interferon Signaling

Rodent-borne orthohantaviruses are asymptomatic in their natural reservoir, but they can cause severe diseases in humans. Although an exacerbated immune response relates to hantaviral pathologies, orthohantaviruses have to antagonize the antiviral interferon (IFN) response to successfully propagate in infected cells. We studied interactions of structural and nonstructural (NSs) proteins of pathogenic Puumala (PUUV), low-pathogenic Tula (TULV), and non-pathogenic Prospect Hill (PHV) viruses, with human type I and III IFN (IFN-I and IFN-III) pathways. The NSs proteins of all three viruses inhibited the RIG-I-activated IFNβ promoter, while only the glycoprotein precursor (GPC) of PUUV, or its cleavage product Gn/Gc, and the nucleocapsid (N) of TULV inhibited it. Moreover, the GPC of both PUUV and TULV antagonized the promoter of IFN-stimulated responsive elements (ISRE). Different viral proteins could thus contribute to inhibition of IFNβ response in a viral context. While PUUV and TULV strains replicated similarly, whether expressing entire or truncated NSs proteins, only PUUV encoding a wild type NSs protein led to late IFN expression and activation of IFN-stimulated genes (ISG). This, together with the identification of particular domains of NSs proteins and different biological processes that are associated with cellular proteins in complex with NSs proteins, suggested that the activation of IFN-I is probably not the only antiviral pathway to be counteracted by orthohantaviruses and that NSs proteins could have multiple inhibitory functions.

Transcriptome analysis of senecavirus A-infected cells: Type I interferon is a critical anti-viral factor

Senecavirus A (SVA)-associated vesicular disease (SAVD) has extensively been present in the swine industry during the past years. The mechanisms of SVA-host interactions at the molecular level, subsequent to SVA infection, are unclear. We studied the gene expression profiles of LLC-PK1 cells, with or without SVA infection, for 6 h and 12 h using an RNA-seq technology. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed on differentially expressed genes (DEGs). Immune-related genes and pathways were significantly modified after SVA infection. To confirm the RNA-seq data, 28 important DEGs were selected for RT-qPCR assays. All DEGs exhibited expression patterns consistent with the RNA-seq results. Among them, type I IFNs (including IFN-α and IFN-β) showed the largest upregulation, followed by RSAD2, DDX58, MX1 and the 17 other DEGs. In contrary, ID2 and another 5 DEGs were down-regulated or unchanged. These results indicated that type I IFNs play a critical role in host immune responses against SVA infection at early stage, while other immune-regulated genes directly or indirectly participate in the host immune responses.

Transcriptome analysis of responses to bluetongue virus infection in Aedes albopictus cells

Background

Bluetongue virus (BTV) causes a disease among wild and domesticated ruminants which is not contagious, but which is transmitted by biting midges of the Culicoides species. BTV can induce an intense cytopathic effect (CPE) in mammalian cells after infection, although Culicoides- or mosquito-derived cell cultures cause non-lytic infection with BTV without CPE. However, little is known about the transcriptome changes in Aedes albopictus cells infected with BTV.

Methods

Transcriptome sequencing was used to identify the expression pattern of mRNA transcripts in A. albopictus cells infected with BTV, given the absence of the Culicoides genome sequence. Bioinformatics analyses were performed to examine the biological functions of the differentially expressed genes. Subsequently, quantitative reverse transcription–polymerase chain reaction was utilized to validate the sequencing data.

Results

In total, 51,850,205 raw reads were generated from the BTV infection group and 51,852,293 from the control group. A total of 5769 unigenes were common to both groups; only 779 unigenes existed exclusively in the infection group and 607 in the control group. In total, 380 differentially expressed genes were identified, 362 of which were up-regulated and 18 of which were down-regulated. Bioinformatics analyses revealed that the differentially expressed genes mainly participated in endocytosis, FoxO, MAPK, dorso-ventral axis formation, insulin resistance, Hippo, and JAK-STAT signaling pathways.

Conclusion

This study represents the first attempt to investigate transcriptome-wide dysregulation in A. albopictus cells infected with BTV. The understanding of BTV pathogenesis and virus–vector interaction will be improved by global transcriptome profiling.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1498-3) contains supplementary material, which is available to authorized users.

Transcriptome profiling in Rift Valley fever virus infected cells reveals modified transcriptional and alternative splicing programs

Rift Valley fever virus (RVFV) is a negative-sense RNA virus belonging to the Phenuiviridae family that infects both domestic livestock and humans. The NIAID has designated RVFV as a Category A priority emerging pathogen due to the devastating public health outcomes associated with epidemic outbreaks. However, there is no licensed treatment or vaccine approved for human use. Therefore it is of great interest to understand RVFV pathogenesis in infected hosts in order to facilitate creation of targeted therapies and treatment options. Here we provide insight into the host-pathogen interface in human HEK293 cells during RVFV MP-12 strain infection using high-throughput mRNA sequencing technology. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of differentially expressed genes showed robust innate immune and cytokine-mediated inflammatory pathway activation as well as alterations in pathways associated with fatty acid metabolism and extracellular matrix receptor signaling. We also analyzed the promoter regions of DEGs for patterns in transcription factor binding sites, and found several that are known to act synergistically to impact apoptosis, immunity, metabolism, and cell growth and differentiation. Lastly, we noted dramatic changes in host alternative splicing patterns in genes associated with mRNA decay and surveillance, RNA transport, and DNA repair. This study has improved our understanding of RVFV pathogenesis and has provided novel insight into pathways and signaling modules important for RVFV diagnostics and therapeutic development.

Comparison of Schmallenberg virus sequences isolated from mammal host and arthropod vector

Schmallenberg virus (SBV) is the member of Peribunyaviridae family, which comprises pathogens of importance for human and veterinary medicine. The virus is transmitted only between animals and mainly by biting midges of the genus Culicoides. This study was performed in order to determine SBV genetic diversity and elucidate the host–vector adaptation. All three viral segments were analysed for sequence variability and phylogenetic relations. The Polish SBV strains obtained from acute infections of cattle, congenital cases in sheep, and from Culicoides midges were sequenced using Sanger and next-generation sequencing (NGS) methods. The obtained sequences were genetically similar (99.2–100% identity) to the first-detected strain BH80/11—4 from German cattle. The sampling year and origin of Polish sequences had no effect on molecular diversity of SBV. Considering all analysed Polish as well as European sequences, ovine-derived sequences were the most variable, while the midge ones were more conserved and encompassed unique substitutions located mainly in nonstructural protein S. SBV sequences isolated from Culicoides are the first submitted to GenBank and reported.

Development of reverse genetics systems and investigation of host response antagonism and reassortment potential for Cache Valley and Kairi viruses, two emerging orthobunyaviruses of the Americas

Orthobunyaviruses such as Cache Valley virus (CVV) and Kairi virus (KRIV) are important animal pathogens. Periodic outbreaks of CVV have resulted in the significant loss of lambs on North American farms, whilst KRIV has mainly been detected in South and Central America with little overlap in geographical range. Vaccines or treatments for these viruses are unavailable. One approach to develop novel vaccine candidates is based on the use of reverse genetics to produce attenuated viruses that elicit immune responses but cannot revert to full virulence. The full genomes of both viruses were sequenced to obtain up to date genome sequence information. Following sequencing, minigenome systems and reverse genetics systems for both CVV and KRIV were developed. Both CVV and KRIV showed a wide in vitro cell host range, with BHK-21 cells a suitable host cell line for virus propagation and titration. To develop attenuated viruses, the open reading frames of the NSs proteins were disrupted. The recombinant viruses with no NSs protein expression induced the production of type I interferon (IFN), indicating that for both viruses NSs functions as an IFN antagonist and that such attenuated viruses could form the basis for attenuated viral vaccines. To assess the potential for reassortment between CVV and KRIV, which could be relevant during vaccination campaigns in areas of overlap, we attempted to produce M segment reassortants by reverse genetics. We were unable to obtain such viruses, suggesting that it is an unlikely event.

Alterations in the host transcriptome in vitro following Rift Valley fever virus infection

Rift Valley fever virus (RVFV) causes major outbreaks among livestock, characterized by "abortion storms" in which spontaneous abortion occurs in almost 100% of pregnant ruminants. Humans can also become infected with mild symptoms that can progress to more severe symptoms, such as hepatitis, encephalitis, and hemorrhagic fever. The goal of this study was to use RNA-sequencing (RNA-seq) to analyze the host transcriptome in response to RVFV infection. G2/M DNA damage checkpoint, ATM signaling, mitochondrial dysfunction, regulation of the antiviral response, and integrin-linked kinase (ILK) signaling were among the top altered canonical pathways with both the attenuated MP12 strain and the fully virulent ZH548 strain. Although several mRNA transcripts were highly upregulated, an increase at the protein level was not observed for the selected genes, which was at least partially due to the NSs dependent block in mRNA export. Inhibition of ILK signaling, which is involved in cell motility and cytoskeletal reorganization, resulted in reduced RVFV replication, indicating that this pathway is important for viral replication. Overall, this is the first global transcriptomic analysis of the human host response following RVFV infection, which could give insight into novel host responses that have not yet been explored.

Fetopathic effects of experimental Schmallenberg virus infection in pregnant goats

Schmallenberg virus (SBV) is an emerging virus responsible for congenital malformations in the offspring of domestic ruminants. It is speculated that infection of pregnant dams may also lead to a significant number of unrecognized fetal losses during the early period of gestation. To assess the pathogenic effects of SBV infection of goats in early pregnancy, we inoculated dams at day 28 or 42 of gestation and followed the animals until day 55 of gestation. Viremia in the absence of clinical signs was detected in all virus-inoculated goats. Fetal deaths were observed in several goats infected at day 28 or 42 of gestation and were invariably associated with the presence of viral genomic RNA in the affected fetuses. Among the viable fetuses, two displayed lesions in the central nervous system (porencephaly) in the presence of viral genome and antigen. All fetuses from goats infected at day 42 and the majority of fetuses from goats infected at day 28 of gestation contained viral genomic RNA. Viral genome was widely distributed in these fetuses and their respective placentas, and infectious virus could be isolated from several organs and placentomes of the viable fetuses. Our results show that fetuses of pregnant goats are susceptible to vertical SBV infection during early pregnancy spanning at least the period between day 28 and 42 of gestation. The outcomes of experimental SBV infection assessed at day 55 of gestation include fetal mortalities, viable fetuses displaying lesions of the central nervous system, as well as viable fetuses without any detectable lesion.

Nonstructural Protein NSs of Schmallenberg Virus Is Targeted to the Nucleolus and Induces Nucleolar Disorganization

Schmallenberg virus (SBV) was discovered in Germany in late 2011 and then spread rapidly to many European countries. SBV is an orthobunyavirus that causes abortion and congenital abnormalities in ruminants. A virus-encoded nonstructural protein, termed NSs, is a major virulence factor of SBV, and it is known to promote the degradation of Rpb1, a subunit of the RNA polymerase II (Pol II) complex, and therefore hampers global cellular transcription. In this study, we found that NSs is mainly localized in the nucleus of infected cells and specifically appears to target the nucleolus through a nucleolar localization signal (NoLS) localized between residues 33 and 51 of the protein. NSs colocalizes with nucleolar markers such as B23 (nucleophosmin) and fibrillarin. We observed that in SBV-infected cells, B23 undergoes a nucleolus-to-nucleoplasm redistribution, evocative of virus-induced nucleolar disruption. In contrast, the nucleolar pattern of B23 was unchanged upon infection with an SBV recombinant mutant with NSs lacking the NoLS motif (SBVΔNoLS). Interestingly, unlike wild-type SBV, the inhibitory activity of SBVΔNoLS toward RNA Pol II transcription is impaired. Overall, our results suggest that a putative link exists between NSs-induced nucleolar disruption and its inhibitory function on cellular transcription, which consequently precludes the cellular antiviral response and/or induces cell death.

IMPORTANCE Schmallenberg virus (SBV) is an emerging arbovirus of ruminants that spread in Europe between 2011 and 2013. SBV induces fetal abnormalities during gestation, with the central nervous system being one of the most affected organs. The virus-encoded NSs protein acts as a virulence factor by impairing host cell transcription. Here, we show that NSs contains a nucleolar localization signal (NoLS) and induces disorganization of the nucleolus. The NoLS motif in the SBV NSs is absolutely necessary for virus-induced inhibition of cellular transcription. To our knowledge, this is the first report of nucleolar functions for NSs within the Bunyaviridae family.

Transcriptome Analysis of HepG2 Cells Expressing ORF3 from Swine Hepatitis E Virus to Determine the Effects of ORF3 on Host Cells

Hepatitis E virus- (HEV-) mediated hepatitis has become a global public health problem. An important regulatory protein of HEV, ORF3, influences multiple signal pathways in host cells. In this study, to investigate the function of ORF3 from the swine form of HEV (SHEV), high-throughput RNA-Seq-based screening was performed to identify the differentially expressed genes in ORF3-expressing HepG2 cells. The results were validated with quantitative real-time PCR and gene ontology was employed to assign differentially expressed genes to functional categories. The results indicated that, in the established ORF3-expressing HepG2 cells, the mRNA levels of CLDN6, YLPM1, APOC3, NLRP1, SCARA3, FGA, FGG, FGB, and FREM1 were upregulated, whereas the mRNA levels of SLC2A3, DKK1, BPIFB2, and PTGR1 were downregulated. The deregulated expression of CLDN6 and FREM1 might contribute to changes in integral membrane protein and basement membrane protein expression, expression changes for NLRP1 might affect the apoptosis of HepG2 cells, and the altered expression of APOC3, SCARA3, and DKK1 may affect lipid metabolism in HepG2 cells. In conclusion, ORF3 plays a functional role in virus-cell interactions by affecting the expression of integral membrane protein and basement membrane proteins and by altering the process of apoptosis and lipid metabolism in host cells. These findings provide important insight into the pathogenic mechanism of HEV.