Conservation of Structure and Immune Antagonist Functions of Filoviral VP35 Homologs Present in Microbat Genomes

A novel approach to exploring the dark genome and its application to mapping of the vertebrate virus fossil record

BACKGROUND

Genomic regions that remain poorly understood, often referred to as the dark genome, contain a variety of functionally relevant and biologically informative features. These include endogenous viral elements (EVEs)-virus-derived sequences that can dramatically impact host biology and serve as a virus fossil record. In this study, we introduce a database-integrated genome screening (DIGS) approach to investigate the dark genome in silico, focusing on EVEs found within vertebrate genomes.

RESULTS

Using DIGS on 874 vertebrate genomes, we uncover approximately 1.1 million EVE sequences, with over 99% originating from endogenous retroviruses or transposable elements that contain EVE DNA. We show that the remaining 6038 sequences represent over a thousand distinct horizontal gene transfer events across 10 virus families, including some that have not previously been reported as EVEs. We explore the genomic and phylogenetic characteristics of non-retroviral EVEs and determine their rates of acquisition during vertebrate evolution. Our study uncovers novel virus diversity, broadens knowledge of virus distribution among vertebrate hosts, and provides new insights into the ecology and evolution of vertebrate viruses.

CONCLUSIONS

We comprehensively catalog and analyze EVEs within 874 vertebrate genomes, shedding light on the distribution, diversity, and long-term evolution of viruses and reveal their extensive impact on vertebrate genome evolution. Our results demonstrate the power of linking a relational database management system to a similarity search-based screening pipeline for in silico exploration of the dark genome.

Marburg Virus VP30 Is Required for Transcription Initiation at the Glycoprotein Gene

Marburg virus (MARV) is an enveloped, negative-sense RNA virus from the filovirus family that causes outbreaks of severe, frequently fatal illness in humans. Of the seven MARV proteins, the VP30 protein stands out because it is essential for viral growth but lacks a definitive function. Here, we used model MARV genome RNAs for one or two reporter genes and the MARV VP40, glycoprotein (GP), and VP24 genes to demonstrate that VP30 is dispensable for the transcription of some genes but critical for transcription reinitiation at the GP gene. This results in the loss of the expression of GP and downstream genes and the impaired production of infectious particles when VP30 is absent. Bicistronic minigenome assays demonstrate that the VP40 gene end/GP gene start junction specifically confers VP30 dependence. A region at the GP gene start site predicted to form a stem-loop contributes to VP30 dependence because the replacement of the GP stem-loop with corresponding sequences from the MARV VP35 gene relieves VP30 dependence. Finally, a Cys3-His zinc binding motif characteristic of filovirus VP30 proteins was demonstrated to be critical for reinitiation at GP. These findings address a long-standing gap in our understanding of MARV biology by defining a critical role for VP30 in MARV transcription. IMPORTANCE Marburg virus and Ebola virus encode VP30 proteins. While the role of VP30 in Ebola virus transcription has been well studied, the role of VP30 in the Marburg virus life cycle is not well understood. The work here demonstrates that different gene start sites within the Marburg viral genome have variable levels of dependence on Marburg virus VP30, with its expression being critical for transcription reinitiation at the GP gene start site. These findings address a long-standing question regarding Marburg virus VP30 function and further our understanding of how Marburg virus gene expression is regulated.

Development and Adoption of Genetically Engineered Plants for Virus Resistance: Advances, Opportunities and Challenges

Plant viruses cause yield losses to crops of agronomic and economic significance and are a challenge to the achievement of global food security. Although conventional plant breeding has played an important role in managing plant viral diseases, it will unlikely meet the challenges posed by the frequent emergence of novel and more virulent viral species or viral strains. Hence there is an urgent need to seek alternative strategies of virus control that can be more readily deployed to contain viral diseases. The discovery in the late 1980s that viral genes can be introduced into plants to engineer resistance to the cognate virus provided a new avenue for virus disease control. Subsequent advances in genomics and biotechnology have led to the refinement and expansion of genetic engineering (GE) strategies in crop improvement. Importantly, many of the drawbacks of conventional breeding, such as long lead times, inability or difficulty to cross fertilize, loss of desirable plant traits, are overcome by GE. Unfortunately, public skepticism towards genetically modified (GM) crops and other factors have dampened the early promise of GE efforts. These concerns are principally about the possible negative effects of transgenes to humans and animals, as well as to the environment. However, with regards to engineering for virus resistance, these risks are overstated given that most virus resistance engineering strategies involve transfer of viral genes or genomic segments to plants. These viral genomes are found in infected plant cells and have not been associated with any adverse effects in humans or animals. Thus, integrating antiviral genes of virus origin into plant genomes is hardly unnatural as suggested by GM crop skeptics. Moreover, advances in deep sequencing have resulted in the sequencing of large numbers of plant genomes and the revelation of widespread endogenization of viral genomes into plant genomes. This has raised the possibility that viral genome endogenization is part of an antiviral defense mechanism deployed by the plant during its evolutionary past. Thus, GM crops engineered for viral resistance would likely be acceptable to the public if regulatory policies were product-based (the North America regulatory model), as opposed to process-based. This review discusses some of the benefits to be gained from adopting GE for virus resistance, as well as the challenges that must be overcome to leverage this technology. Furthermore, regulatory policies impacting virus-resistant GM crops and some success cases of virus-resistant GM crops approved so far for cultivation are discussed.

Virus-like insertions with sequence signatures similar to those of endogenous nonretroviral RNA viruses in the human genome

Understanding the genetics and taxonomy of ancient viruses will give us great insights into not only the origin and evolution of viruses but also how viral infections played roles in our evolution. Endogenous viruses are remnants of ancient viral infections and are thought to retain the genetic characteristics of viruses from ancient times. In this study, we used machine learning of endogenous RNA virus sequence signatures to identify viruses in the human genome that have not been detected or are already extinct. Here, we show that the k-mer occurrence of ancient RNA viral sequences remains similar to that of extant RNA viral sequences and can be differentiated from that of other human genome sequences. Furthermore, using this characteristic, we screened RNA viral insertions in the human reference genome and found virus-like insertions with phylogenetic and evolutionary features indicative of an exogenous origin but lacking homology to previously identified sequences. Our analysis indicates that animal genomes still contain unknown virus-derived sequences and provides a glimpse into the diversity of the ancient virosphere.

Metagenomic Snapshots of Viral Components in Guinean Bats

To prevent the emergence of zoonotic infectious diseases and reduce their epidemic potential, we need to understand their origins in nature. Bats in the order Chiroptera are widely distributed worldwide and are natural reservoirs of prominent zoonotic viruses, including Nipah virus, Marburg virus, and possibly SARS-CoV-2. In this study, we applied unbiased metagenomic and metatranscriptomic approaches to decipher the virosphere of frugivorous and insectivorous bat species captured in Guéckédou, Guinea, the epicenter of the West African Ebola virus disease epidemic in 2013–2016. Our study provides a snapshot of the viral diversity present in these bat species, with several novel viruses reported for the first time in bats, as well as some bat viruses closely related to known human or animal pathogens. In addition, analysis of Mops condylurus genomic DNA samples revealed the presence of an Ebola virus nucleoprotein (NP)-derived pseudogene inserted in its genome. These findings provide insight into the evolutionary traits of several virus families in bats and add evidence that nonretroviral integrated RNA viruses (NIRVs) derived from filoviruses may be common in bat genomes.

Impact of Měnglà Virus Proteins on Human and Bat Innate Immune Pathways

Měnglà virus (MLAV), identified in Rousettus bats, is a phylogenetically distinct member of the family Filoviridae Because the filoviruses Ebola virus (EBOV) and Marburg virus (MARV) modulate host innate immunity, MLAV VP35, VP40, and VP24 proteins were compared with their EBOV and MARV homologs for innate immune pathway modulation. In human and Rousettus cells, MLAV VP35 behaved like EBOV and MARV VP35s, inhibiting virus-induced activation of the interferon beta (IFN-β) promoter and interferon regulatory factor 3 (IRF3) phosphorylation. MLAV VP35 also interacted with PACT, a host protein engaged by EBOV VP35 to inhibit RIG-I signaling. MLAV VP35 also inhibits PKR activation. MLAV VP40 was demonstrated to inhibit type I IFN-induced gene expression in human and bat cells. It blocked STAT1 tyrosine phosphorylation induced either by type I IFN or overexpressed Jak1, paralleling MARV VP40. MLAV VP40 also inhibited virus-induced IFN-β promoter activation, a property shared by MARV VP40 and EBOV VP24. A Jak kinase inhibitor did not recapitulate this inhibition in the absence of viral proteins. Therefore, inhibition of Jak-STAT signaling is insufficient to explain inhibition of IFN-β promoter activation. MLAV VP24 did not inhibit IFN-induced gene expression or bind karyopherin α proteins, properties of EBOV VP24. MLAV VP24 differed from MARV VP24 in that it failed to interact with Keap1 or activate an antioxidant response element reporter gene due to the absence of a Keap1-binding motif. These functional observations support a closer relationship of MLAV to MARV than to EBOV but also are consistent with MLAV belonging to a distinct genus.IMPORTANCE EBOV and MARV, members of the family Filoviridae, are highly pathogenic zoonotic viruses that cause severe disease in humans. Both viruses use several mechanisms to modulate the host innate immune response, and these likely contribute to the severity of disease. Here, we demonstrate that MLAV, a filovirus newly discovered in a bat, suppresses antiviral type I interferon responses in both human and bat cells. Inhibitory activities are possessed by MLAV VP35 and VP40, which parallels how MARV blocks IFN responses. However, whereas MARV activates cellular antioxidant responses through an interaction between its VP24 protein and host protein Keap1, MLAV VP24 lacks a Keap1-binding motif and fails to activate this cytoprotective response. These data indicate that MLAV possesses immune-suppressing functions that could facilitate human infection. They also support the placement of MLAV in a different genus than either EBOV or MARV.

Interactions among eukaryotes, retrotransposons and riboviruses: endogenous riboviral elements in eukaryotic genomes

Riboviruses are viruses that have RNA genomes and replicate only via RNA intermediates. Although they do not require a DNA phase for replication and do not encode reverse transcriptase, the presence of DNA forms of riboviral sequences in ribovirus-infected cells has been reported since the 1970s. Additionally, heritable ribovirus-derived sequences, called riboviral endogenous viral elements (EVEs), have been found in the genomes of many eukaryotes. These are now thought to be formed by the reverse transcription machineries of retrotransposons within eukaryotic genomes sometimes referred to as selfish elements. Surprisingly, some reverse-transcribed riboviral DNAs (including EVEs) provide physiological functions for their hosts, suggesting the occurrence of novel interactions among eukaryotic genomes, retrotransposons and riboviruses, and opening the door to new avenues of investigation. Here I review current knowledge on these triangular interactions, and discuss future directions in this field.

Evolutionary persistence of insect bunyavirus infection despite host acquisition and expression of the viral nucleoprotein gene

How insects combat RNA virus infection is a subject of intensive research owing to its importance in insect health, virus evolution, and disease transmission. In recent years, a pair of potentially linked phenomena have come to light as a result of this work-first, the pervasive production of viral DNA from exogenous nonretroviral RNA in infected individuals, and second, the widespread distribution of nonretroviral integrated RNA virus sequences (NIRVs) in the genomes of diverse eukaryotes. The evolutionary consequences of NIRVs for viruses are unclear and the field would benefit from studies of natural virus infections co-occurring with recent integrations, an exceedingly rare circumstance in the literature. Here, we provide evidence that a novel insect-infecting phasmavirus (Order Bunyavirales) has been persisting in a phantom midge host, Chaoborus americanus, for millions of years. Interestingly, the infection persists despite the host's acquisition (during the Pliocene), fixation, and expression of the viral nucleoprotein gene. We show that virus prevalence and geographic distribution are high and broad, comparable to the host-specific infections reported in other phantom midges. Short-read mapping analyses identified a lower abundance of the nucleoprotein-encoding genome segment in this virus relative to related viruses. Finally, the novel virus has facilitated the first substitution rate estimation for insect-infecting phasmaviruses. Over a period of approximately 16 million years, we find rates of (0.6 - 1.6) × 10^-7 substitutions per site per year in protein coding genes, extraordinarily low for negative-sense RNA viruses, but consistent with the few estimates produced over comparable evolutionary timescales.

Relevance of Ebola virus VP35 homo-dimerization on the type I interferon cascade inhibition

Ebola virus high lethality relies on its ability to efficiently bypass the host innate antiviral response, which senses the viral dsRNA through the RIG-I receptor and induces type I interferon α/β production. In the bypassing action, the Ebola virus protein VP35 plays a pivotal role at multiple levels of the RIG-I cascade, masking the viral 5′-triphosphorylated dsRNA from RIG-I, and interacting with other cascade components. The VP35 type I interferon inhibition is exerted by the C-terminal domain, while the N-terminal domain, containing a coiled-coil region, is primarily required for oligomerization. However, mutations at key VP35 residues L90/93/107A (VP35-3m) in the coiled-coil region were reported to affect oligomerization and reduce type I interferon antagonism, indicating a possible but unclear role of homo-oligomerization on VP35 interaction with the RIG-I pathway components. In this work, we investigated the VP35 dimerization thermodynamics and its contribution to type I interferon antagonism by computational and biological methods. Focusing on the coiled-coil region, we combined coarse-grained and all-atom simulations on wild type VP35 and VP35-3m homo-dimerization. According to our results, wild type VP35 coiled-coil is able to self-assemble into dimers, while VP35-3m coiled-coil shows poor propensity to even dimerize. Free-energy calculations confirmed the key role of L90, L93 and L107 in stabilizing the coiled-coil homo-dimeric structure. In vitro type I interferon antagonism studies, using full-length wild type VP35 and VP35-3m, revealed that VP35 homo-dimerization is an essential preliminary step for dsRNA binding, which appears to be the main factor of the VP35 RIG-I cascade inhibition, while it is not essential to block the other steps.