Putative endogenous filovirus VP35-like protein potentially functions as an IFN antagonist but not a polymerase cofactor

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.

Roles of Human Endogenous Retroviruses and Endogenous Virus-Like Elements in Cancer Development and Innate Immunity

Human endogenous retroviruses (HERVs) are remnants of ancient retroviral infections in the host genome. Although mutations and silencing mechanisms impair their original role in viral replication, HERVs are believed to play roles in various biological processes. Long interspersed nuclear elements (LINEs) are non-LTR retrotransposons that have a lifecycle resembling that of retroviruses. Although LINE expression is typically silenced in somatic cells, it also contributes to various biological processes. The aberrant expression of HERVs and LINEs is closely associated with the development of cancer and/or immunological diseases, suggesting that they are integrated into various pathways related to the diseases. HERVs/LINEs control gene expression depending on the context as promoter/enhancer elements. Some RNAs and proteins derived from HERVs/LINEs have oncogenic potential, whereas others stimulate innate immunity. Non-retroviral endogenous viral elements (nrEVEs) are a novel type of virus-like element in the genome. nrEVEs may also be involved in host immunity. This article provides a current understanding of how these elements impact cellular physiology in cancer development and innate immunity, and provides perspectives for future studies.

The Evolution and Expression Profiles of EC1 Gene Family during Development in Cotton

Fertilization is essential to sexual reproduction of flowering plants. EC1 (EGG CELL 1) proteins have a conserved cysteine spacer characteristic and play a crucial role in double fertilization process in many plant species. However, to date, the role of EC1 gene family in cotton is fully unknown. Hence, detailed bioinformatics analysis was explored to elucidate the biological mechanisms of EC1 gene family in cotton. In this study, we identified 66 genes in 10 plant species in which a total of 39 EC1 genes were detected from cotton genome. Phylogenetic analysis clustered the identified EC1 genes into three families (I-III) and all of them contain Prolamin-like domains. A good collinearity was observed in the synteny analysis of the orthologs from cotton genomes. Whole-genome duplication was determined to be one of the major impetuses for the expansion of the EC1 gene family during the process of evolution. qRT-PCR analysis showed that EC1 genes were highly expressed in reproductive tissues under multiple stresses, signifying their potential role in enhancing stress tolerance or responses. Additionally, gene interaction networks showed that EC1 genes may be involved in cell stress and response transcriptional regulator in the synergid cells and activate the expression of genes required for pollen tube guidance. Our results provide novel functional insights into the evolution and functional elucidation of EC1 gene family in cotton.

Functional Importance of Hydrophobic Patches on the Ebola Virus VP35 IFN-Inhibitory Domain

Viral protein 35 (VP35) of Ebola virus (EBOV) is a multifunctional protein that mainly acts as a viral polymerase cofactor and an interferon antagonist. VP35 interacts with the viral nucleoprotein (NP) and double-stranded RNA for viral RNA transcription/replication and inhibition of type I interferon (IFN) production, respectively. The C-terminal portion of VP35, which is termed the IFN-inhibitory domain (IID), is important for both functions. To further identify critical regions in this domain, we analyzed the physical properties of the surface of VP35 IID, focusing on hydrophobic patches, which are expected to be functional sites that are involved in interactions with other molecules. Based on the known structural information of VP35 IID, three hydrophobic patches were identified on its surface and their biological importance was investigated using minigenome and IFN-β promoter-reporter assays. Site-directed mutagenesis revealed that some of the amino acid substitutions that were predicted to disrupt the hydrophobicity of the patches significantly decreased the efficiency of viral genome replication/transcription due to reduced interaction with NP, suggesting that the hydrophobic patches might be critical for the formation of a replication complex through the interaction with NP. It was also found that the hydrophobic patches were involved in the IFN-inhibitory function of VP35. These results highlight the importance of hydrophobic patches on the surface of EBOV VP35 IID and also indicate that patch analysis is useful for the identification of amino acid residues that directly contribute to protein functions.

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.

Viral-derived DNA invasion and individual variation in an Indonesian population of large flying fox Pteropus vampyrus

Here, we performed next-generation sequencing (NGS) on six large flying foxes (Pteropus vampyrus) collected in Indonesia. Seventy-five virus species in the liver tissue of each specimen were listed. Viral homologous sequences in the bat genome were identified from the listed viruses. This finding provides collateral evidence of viral endogenization into the host genome. We found that two of the six specimens bore partial sequences that were homologous to the plant pathogens Geminiviridae and Luteoviridae. These sequences were absent in the P. vampyrus chromosomal sequences. Hence, plant viral homologous sequences were localized to the hepatocytes as extrachromosomal DNA fragments. Therefore, this suggests that the bat is a potential carrier or vector of plant viruses. The present investigation on wild animals offered novel perspectives on viral invasion, variation, and host interaction.

The Potential Role of Endogenous Viral Elements in the Evolution of Bats as Reservoirs for Zoonotic Viruses

Despite a small genome size, bats have comparable diversity of retroviral and non-retroviral endogenous sequences to other mammals. These include Class I and Class II retroviral sequences, foamy viruses, and deltaretroviruses, as well as filovirus, bornavirus, and parvovirus endogenous viral elements. Some of these endogenous viruses are sufficiently preserved in bat genomes to be expressed, with potential effects for host biology. It is clear that the bat immune system differs when compared with other mammals, yet the role that virus-derived endogenous elements may have played in the evolution of bat immunity is poorly understood. In this review, we discuss some of the bat-specific immune mechanisms that may have resulted in a virus-tolerant phenotype and link these to the long-standing virus-host coevolution that may have allowed a large diversity of endogenous retroviruses and other endogenous viral elements to colonize bat genomes. We also consider the possible effects of endogenization in the evolution of the bat immune system.

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.

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

Non-retroviral integrated RNA viral sequences (NIRVs) potentially encoding ∼280 amino acid homologs to filovirus VP35 proteins are present across the Myotis genus of bats. These are estimated to have been maintained for ∼18 million years, indicating their co-option. To address the reasons for co-option, 16 Myotis VP35s were characterized in comparison to VP35s from the extant filoviruses Ebola virus and Marburg virus, in which VP35s play critical roles in immune evasion and RNA synthesis. The Myotis VP35s demonstrated a conserved suppression of innate immune signaling, albeit with reduced potency, in either human or Myotis cells. Their attenuation reflects a lack of dsRNA binding that in the filoviral VP35s correlates with potent suppression of interferon responses. Despite divergent function, evolution has preserved in Myotis the structure of the filoviral VP35s, indicating that this structure is critical for co-opted function, possibly as a regulator of innate immune signaling.

[Virus-host coevolution: Endogenous RNA viral elements as pseudogenes]

RNA viruses do not need to take the form of DNAs, and RNAs alone complete their replication cycles. On the other hand, since the 1970s, it has been known that DNA fragments derived from RNA viruses can be detected in RNA virus-infected cells. Furthermore, in this decade, it has become clear that the eukaryotic genomes contain genetic sequences derived from non-retroviral RNA viruses. The DNA sequences derived from these RNA viruses are thought to be generatedby using a transposable mechanism of retrotransposon, such as LINE-1. Many endogenous RNA viral sequences are formed by the same mechanism as processed pseudogenes in eukaryotic cells, but the significance of the production of RNA viral "pseudogenes " in infected cells has not been elucidated. We have discovered endogenous bornavirus-like elements (EBLs), which derived from a negative-sense, single-stranded RNA virus, Bornaviruses, and have studied the evolution and function of EBLs in host animals. The analysis of EBLs provides us a clue to unravel the history of host-RNA virus coexistence. In this review, I overview about the function of endogenous RNA virus sequences, especially EBLs in mammalian genomes, and discuss the significance of endogenization of RNA viruses as viral pseudogenes in evolution.

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.

Cellular production of a counterfeit viral protein confers immunity to infection by a related virus

DNA copies of many non-retroviral RNA virus genes or portions thereof (NIRVs) are present in the nuclear genomes of many eukaryotes. These have often been preserved for millions of years of evolution, suggesting that they play an important cellular function. One possible function is resistance to infection by related viruses. In some cases, this appears to occur through the piRNA system, but in others by way of counterfeit viral proteins encoded by NIRVs. In the fungi, NIRVs may be as long as 1,400 uninterrupted codons. In one such case in the yeast Debaryomyces hansenii, one of these genes provides immunity to a related virus by virtue of expression of a counterfeit viral capsid protein, which interferes with assembly of viral capsids by negative complementation. The widespread occurrence of non-retroviral RNA virus genes in eukaryotes may reflect an underappreciated method of host resistance to infection. This work demonstrates for the first time that an endogenous host protein encoded by a gene that has been naturally acquired from a virus and fixed in a eukaryote can interfere with the replication of a related virus and do so by negative complementation.