Uncovering the Repertoire of Endogenous Flaviviral Elements in Aedes Mosquito Genomes

Endogenous viral elements derived from nonretroviral RNA viruses have been described in various animal genomes. Whether they have a biological function, such as host immune protection against related viruses, is a field of intense study. Here, we investigated the repertoire of endogenous flaviviral elements (EFVEs) in Aedes mosquitoes, the vectors of arboviruses such as dengue and chikungunya viruses. Previous studies identified three EFVEs from Aedes albopictus cell lines and one from Aedes aegypti cell lines. However, an in-depth characterization of EFVEs in wild-type mosquito populations and individual mosquitoes in vivo has not been performed. We detected the full-length DNA sequence of the previously described EFVEs and their respective transcripts in several A. albopictus and A. aegypti populations from geographically distinct areas. However, EFVE-derived proteins were not detected by mass spectrometry. Using deep sequencing, we detected the production of PIWI-interacting RNA-like small RNAs, in an antisense orientation, targeting the EFVEs and their flanking regions in vivo. The EFVEs were integrated in repetitive regions of the mosquito genomes, and their flanking sequences varied among mosquito populations. We bioinformatically predicted several new EFVEs from a Vietnamese A. albopictus population and observed variation in the occurrence of those elements among mosquitoes. Phylogenetic analysis of an A. aegypti EFVE suggested that it integrated prior to the global expansion of the species and subsequently diverged among and within populations. The findings of this study together reveal the substantial structural and nucleotide diversity of flaviviral integrations in Aedes genomes. Unraveling this diversity will help to elucidate the potential biological function of these EFVEs.

IMPORTANCE Endogenous viral elements (EVEs) are whole or partial viral sequences integrated in host genomes. Interestingly, some EVEs have important functions for host fitness and antiviral defense. Because mosquitoes also have EVEs in their genomes, characterizing these EVEs is a prerequisite for their potential use to manipulate the mosquito antiviral response. In the study described here, we focused on EVEs related to the Flavivirus genus, to which dengue and Zika viruses belong, in individual Aedes mosquitoes from geographically distinct areas. We show the existence in vivo of flaviviral EVEs previously identified in mosquito cell lines, and we detected new ones. We show that EVEs have evolved differently in each mosquito population. They produce transcripts and small RNAs but not proteins, suggesting a function at the RNA level. Our study uncovers the diverse repertoire of flaviviral EVEs in Aedes mosquito populations and contributes to an understanding of their role in the host antiviral system.

Virus Latency and the Impact on Plants

Plant viruses are thought to be essentially harmful to the lives of their cultivated crop hosts. In most cases studied, the interaction between viruses and cultivated crop plants negatively affects host morphology and physiology, thereby resulting in disease. Native wild/non-cultivated plants are often latently infected with viruses without any clear symptoms. Although seemingly non-harmful, these viruses pose a threat to cultivated crops because they can be transmitted by vectors and cause disease. Reports are accumulating on infections with latent plant viruses that do not cause disease but rather seem to be beneficial to the lives of wild host plants. In a few cases, viral latency involves the integration of full-length genome copies into the host genome that, in response to environmental stress or during certain developmental stages of host plants, can become activated to generate and replicate episomal copies, a transition from latency to reactivation and causation of disease development. The interaction between viruses and host plants may also lead to the integration of partial-length segments of viral DNA genomes or copy DNA of viral RNA genome sequences into the host genome. Transcripts derived from such integrated viral elements (EVEs) may be beneficial to host plants, for example, by conferring levels of virus resistance and/or causing persistence/latency of viral infections. Studies on viral latency in wild host plants might help us to understand and elucidate the underlying mechanisms of latency and provide insights into the raison d’être for viruses in the lives of plants.

Endogenous hepadnaviruses, bornaviruses and circoviruses in snakes

We report the discovery of endogenous viral elements (EVEs) from Hepadnaviridae, Bornaviridae and Circoviridae in the speckled rattlesnake, Crotalus mitchellii, the first viperid snake for which a draft whole genome sequence assembly is available. Analysis of the draft assembly reveals genome fragments from the three virus families were inserted into the genome of this snake over the past 50 Myr. Cross-species PCR screening of orthologous loci and computational scanning of the python and king cobra genomes reveals that circoviruses integrated most recently (within the last approx. 10 Myr), whereas bornaviruses and hepadnaviruses integrated at least approximately 13 and approximately 50 Ma, respectively. This is, to our knowledge, the first report of circo-, borna- and hepadnaviruses in snakes and the first characterization of non-retroviral EVEs in non-avian reptiles. Our study provides a window into the historical dynamics of viruses in these host lineages and shows that their evolution involved multiple host-switches between mammals and reptiles.

The Interplay Between Viruses and RNAi Pathways in Insects

As an overarching immune mechanism, RNA interference (RNAi) displays pathogen specificity and memory via different pathways. The small interfering RNA (siRNA) pathway is the primary antiviral defense mechanism against RNA viruses of insects and plays a lesser role in defense against DNA viruses. Reflecting the pivotal role of the siRNA pathway in virus selection, different virus families have independently evolved unique strategies to counter this host response, including protein-mediated, decoy RNA-based, and microRNA-based strategies. In this review, we outline the interplay between insect viruses and the different pathways of the RNAi antiviral response; describe practical application of these interactions for improved expression systems and for pest and disease management; and highlight research avenues for advancement of the field.

piRNAs derived from ancient viral processed pseudogenes as transgenerational sequence-specific immune memory in mammals

Endogenous bornavirus-like nucleoprotein elements (EBLNs) are sequences within vertebrate genomes derived from reverse transcription and integration of ancient bornaviral nucleoprotein mRNA via the host retrotransposon machinery. While species with EBLNs appear relatively resistant to bornaviral disease, the nature of this association is unclear. We hypothesized that EBLNs could give rise to antiviral interfering RNA in the form of PIWI-interacting RNAs (piRNAs), a class of small RNA known to silence transposons but not exogenous viruses. We found that in both rodents and primates, which acquired their EBLNs independently some 25-40 million years ago, EBLNs are present within piRNA-generating regions of the genome far more often than expected by chance alone (ℙ = 8 × 10(-3)-6 × 10(-8)). Three of the seven human EBLNs fall within annotated piRNA clusters and two marmoset EBLNs give rise to bona fide piRNAs. In both rats and mice, at least two of the five EBLNs give rise to abundant piRNAs in the male gonad. While no EBLNs are syntenic between rodent and primate, some of the piRNA clusters containing EBLNs are; thus we deduce that EBLNs were integrated into existing piRNA clusters. All true piRNAs derived from EBLNs are antisense relative to the proposed ancient bornaviral nucleoprotein mRNA. These observations are consistent with a role for EBLN-derived piRNA-like RNAs in interfering with ancient bornaviral infection. They raise the hypothesis that retrotransposon-dependent virus-to-host gene flow could engender RNA-mediated, sequence-specific antiviral immune memory in metazoans analogous to the CRISPR/Cas system in prokaryotes.

Natural Variation in Resistance to Virus Infection in Dipteran Insects

The power and ease of Drosophila genetics and the medical relevance of mosquito-transmitted viruses have made dipterans important model organisms in antiviral immunology. Studies of virus-host interactions at the molecular and population levels have illuminated determinants of resistance to virus infection. Here, we review the sources and nature of variation in antiviral immunity and virus susceptibility in model dipteran insects, specifically the fruit fly Drosophila melanogaster and vector mosquitoes of the genera Aedes and Culex. We first discuss antiviral immune mechanisms and describe the virus-specificity of these responses. In the following sections, we review genetic and microbiota-dependent variation in antiviral immunity. In the final sections, we explore less well-studied sources of variation, including abiotic factors, sexual dimorphism, infection history, and endogenous viral elements. We borrow from work on other pathogen types and non-dipteran species when it parallels or complements studies in dipterans. Understanding natural variation in virus-host interactions may lead to the identification of novel restriction factors and immune mechanisms and shed light on the molecular determinants of vector competence.

An Ancient Lineage of Highly Divergent Parvoviruses Infects both Vertebrate and Invertebrate Hosts

Chapparvoviruses (ChPVs) comprise a divergent, recently identified group of parvoviruses (family Parvoviridae), associated with nephropathy in immunocompromised laboratory mice and with prevalence in deep sequencing results of livestock showing diarrhea. Here, we investigate the biological and evolutionary characteristics of ChPVs via comparative in silico analyses, incorporating sequences derived from endogenous parvoviral elements (EPVs) as well as exogenous parvoviruses. We show that ChPVs are an ancient lineage within the Parvoviridae, clustering separately from members of both currently established subfamilies. Consistent with this, they exhibit a number of characteristic features, including several putative auxiliary protein-encoding genes, and capsid proteins with no sequence-level homology to those of other parvoviruses. Homology modeling indicates the absence of a β-A strand, normally part of the luminal side of the parvoviral capsid protein core. Our findings demonstrate that the ChPV lineage infects an exceptionally broad range of host species, including both vertebrates and invertebrates. Furthermore, we observe that ChPVs found in fish are more closely related to those from invertebrates than they are to those of amniote vertebrates. This suggests that transmission between distantly related host species may have occurred in the past and that the Parvoviridae family can no longer be divided based on host affiliation.

ViralRecall-A Flexible Command-Line Tool for the Detection of Giant Virus Signatures in 'Omic Data

Giant viruses are widespread in the biosphere and play important roles in biogeochemical cycling and host genome evolution. Also known as nucleo-cytoplasmic large DNA viruses (NCLDVs), these eukaryotic viruses harbor the largest and most complex viral genomes known. Studies have shown that NCLDVs are frequently abundant in metagenomic datasets, and that sequences derived from these viruses can also be found endogenized in diverse eukaryotic genomes. The accurate detection of sequences derived from NCLDVs is therefore of great importance, but this task is challenging owing to both the high level of sequence divergence between NCLDV families and the extraordinarily high diversity of genes encoded in their genomes, including some encoding for metabolic or translation-related functions that are typically found only in cellular lineages. Here, we present ViralRecall, a bioinformatic tool for the identification of NCLDV signatures in ‘omic data. This tool leverages a library of giant virus orthologous groups (GVOGs) to identify sequences that bear signatures of NCLDVs. We demonstrate that this tool can effectively identify NCLDV sequences with high sensitivity and specificity. Moreover, we show that it can be useful both for removing contaminating sequences in metagenome-assembled viral genomes as well as the identification of eukaryotic genomic loci that derived from NCLDV. ViralRecall is written in Python 3.5 and is freely available on GitHub: https://github.com/faylward/viralrecall.

Population genomics in the arboviral vector Aedes aegypti reveals the genomic architecture and evolution of endogenous viral elements

Horizontal gene transfer from viruses to eukaryotic cells is a pervasive phenomenon. Somatic viral integrations are linked to persistent viral infection whereas integrations into germline cells are maintained in host genomes by vertical transmission and may be co-opted for host functions. In the arboviral vector Aedes aegypti, an endogenous viral element from a nonretroviral RNA virus (nrEVE) was shown to produce PIWI-interacting RNAs (piRNAs) to limit infection with a cognate virus. Thus, nrEVEs may constitute a heritable, sequence-specific mechanism for antiviral immunity, analogous to piRNA-mediated silencing of transposable elements. Here, we combine population genomics and evolutionary approaches to analyse the genomic architecture of nrEVEs in A. aegypti. We conducted a genome-wide screen for adaptive nrEVEs and searched for novel population-specific nrEVEs in the genomes of 80 individual wild-caught mosquitoes from five geographical populations. We show a dynamic landscape of nrEVEs in mosquito genomes and identified five novel nrEVEs derived from two currently circulating viruses, providing evidence of the environmental-dependent modification of a piRNA cluster. Overall, our results show that virus endogenization events are complex with only a few nrEVEs contributing to adaptive evolution in A. aegypti.

A single unidirectional piRNA cluster similar to the flamenco locus is the major source of EVE-derived transcription and small RNAs in Aedes aegypti mosquitoes

Endogenous viral elements (EVEs) are found in many eukaryotic genomes. Despite considerable knowledge about genomic elements such as transposons (TEs) and retroviruses, we still lack information about nonretroviral EVEs. Aedes aegypti mosquitoes have a highly repetitive genome that is covered with EVEs. Here, we identified 129 nonretroviral EVEs in the AaegL5 version of the A. aegypti genome. These EVEs were significantly associated with TEs and preferentially located in repeat-rich clusters within intergenic regions. Genome-wide transcriptome analysis showed that most EVEs generated transcripts although only around 1.4% were sense RNAs. The majority of EVE transcription was antisense and correlated with the generation of EVE-derived small RNAs. A single genomic cluster of EVEs located in a 143 kb repetitive region in chromosome 2 contributed with 42% of antisense transcription and 45% of small RNAs derived from viral elements. This region was enriched for TE-EVE hybrids organized in the same coding strand. These generated a single long antisense transcript that correlated with the generation of phased primary PIWI-interacting RNAs (piRNAs). The putative promoter of this region had a conserved binding site for the transcription factor Cubitus interruptus, a key regulator of the flamenco locus in Drosophila melanogaster Here, we have identified a single unidirectional piRNA cluster in the A. aegypti genome that is the major source of EVE transcription fueling the generation of antisense small RNAs in mosquitoes. We propose that this region is a flamenco-like locus in A. aegypti due to its relatedness to the major unidirectional piRNA cluster in Drosophila melanogaster.

Assessing the Diversity of Endogenous Viruses Throughout Ant Genomes

Endogenous viral elements (EVEs) can play a significant role in the evolution of their hosts and have been identified in animals, plants, and fungi. Additionally, EVEs potentially provide an important snapshot of the evolutionary frequency of viral infection. The purpose of this study is to take a comparative host-centered approach to EVE discovery in ant genomes to better understand the relationship of EVEs to their ant hosts. Using a comprehensive bioinformatic pipeline, we screened all nineteen published ant genomes for EVEs. Once the EVEs were identified, we assessed their phylogenetic relationships to other closely related exogenous viruses. A diverse group of EVEs were discovered in all screened ant host genomes and in many cases are similar to previously identified exogenous viruses. EVEs similar to ssRNA viral proteins are the most common viral lineage throughout the ant hosts, which is potentially due to more chronic infection or more effective endogenization of certain ssRNA viruses in ants. In addition, both EVEs similar to viral glycoproteins and retrovirus-derived proteins are also abundant throughout ant genomes, suggesting their tendency to endogenize. Several of these newly discovered EVEs are found to be potentially functional within the genome. The discovery and analysis of EVEs is essential in beginning to understand viral–ant interactions over evolutionary time.

Nudivirus Remnants in the Genomes of Arthropods

Endogenous viral elements (EVEs), derived from all major types of viruses, have been discovered in many eukaryotic genomes, representing "fossil records" of past viral infections. The endogenization of nudiviruses has been reported in several insects, leading to the question of whether genomic integration is a common phenomenon for these viruses. In this study, genomic assemblies of insects and other arthropods were analyzed to identify endogenous sequences related to Nudiviridae. A total of 359 nudivirus-like genes were identified in 43 species belonging to different groups; however, none of these genes were detected in the known hosts of nudiviruses. A large proportion of the putative EVEs identified in this study encode intact open reading frames or are transcribed as mRNAs, suggesting that they result from recent endogenization of nudiviruses. Phylogenetic analyses of the identified EVEs and inspections of their flanking regions indicated that integration of nudiviruses has occurred recurrently during the evolution of arthropods. This is the first report of a comprehensive screening for nudivirus-derived EVEs in arthropod genomes. The results of this study demonstrated that a large variety of arthropods, especially hemipteran and hymenopteran insects, have previously been or are still infected by nudiviruses. These findings have greatly extended the host range of Nudiviridae and provide new insights into viral diversity, evolution, and host-virus interactions.

Potential Links between Hepadnavirus and Bornavirus Sequences in the Host Genome and Cancer

Various viruses leave their sequences in the host genomes during infection. Such events occur mainly in retrovirus infection but also sometimes in DNA and non-retroviral RNA virus infections. If viral sequences are integrated into the genomes of germ line cells, the sequences can become inherited as endogenous viral elements (EVEs). The integration events of viral sequences may have oncogenic potential. Because proviral integrations of some retroviruses and/or reactivation of endogenous retroviruses are closely linked to cancers, viral insertions related to non-retroviral viruses also possibly contribute to cancer development. This article focuses on genomic viral sequences derived from two non-retroviral viruses, whose endogenization is already reported, and discusses their possible contributions to cancer. Viral insertions of hepatitis B virus play roles in the development of hepatocellular carcinoma. Endogenous bornavirus-like elements, the only non-retroviral RNA virus-related EVEs found in the human genome, may also be involved in cancer formation. In addition, the possible contribution of the interactions between viruses and retrotransposons, which seem to be a major driving force for generating EVEs related to non-retroviral RNA viruses, to cancers will be discussed. Future studies regarding the possible links described here may open a new avenue for the development of novel therapeutics for tumor virus-related cancers and/or provide novel insights into EVE functions.

A Mutualistic Poxvirus Exhibits Convergent Evolution with Other Heritable Viruses in Parasitoid Wasps

For insects known as parasitoid wasps, successful development as a parasite results in the death of the host insect. As a result of this lethal interaction, wasps and their hosts have coevolved strategies to gain an advantage in this evolutionary arms race. Although normally considered to be strict pathogens, some viruses have established persistent infections within parasitoid wasp lineages and are beneficial to wasps during parasitism. Heritable associations between viruses and parasitoid wasps have evolved independently multiple times, but most of these systems remain largely understudied with respect to viral origin, transmission and replication strategies of the virus, and interactions between the virus and host insects. Here, we report a detailed characterization of Diachasmimorpha longicaudata entomopoxvirus (DlEPV), a poxvirus found within the venom gland of Diachasmimorpha longicaudata wasps. Our results show that DlEPV exhibits similar but distinct transmission and replication dynamics compared to those of other parasitoid viral elements, including vertical transmission of the virus within wasps, as well as virus replication in both female wasps and fruit fly hosts. Functional assays demonstrate that DlEPV is highly virulent within fly hosts, and wasps without DlEPV have severely reduced parasitism success compared to those with a typical viral load. Taken together, the data presented in this study illustrate a novel case of beneficial virus evolution, in which a virus of unique origin has undergone convergent evolution with other viral elements associated with parasitoid wasps to provide an analogous function throughout parasitism.IMPORTANCE Viruses are generally considered to be disease-causing agents, but several instances of beneficial viral elements have been identified in insects called parasitoid wasps. These virus-derived entities are passed on through wasp generations and enhance the success of the wasps' parasitic life cycle. Many parasitoid-virus partnerships studied to date exhibit common features among independent cases of this phenomenon, including a mother-to-offspring route of virus transmission, a restricted time and location for virus replication, and a positive effect of virus activity on wasp survival. Our characterization of Diachasmimorpha longicaudata entomopoxvirus (DlEPV), a poxvirus found in Diachasmimorpha longicaudata parasitoid wasps, represents a novel example of beneficial virus evolution. Here, we show that DlEPV exhibits functional similarities to known parasitoid viral elements that support its comparable role during parasitism. Our results also demonstrate unique differences that suggest DlEPV is more autonomous than other long-term viral associations described in parasitoid wasps.

Patchy DNA forms of the Zika virus RNA genome are generated following infection in mosquito cell cultures and in mosquitoes

Zika virus (ZIKV) is a mosquito-borne flavivirus and has historically been reported to cause mild symptomatic diseases during human infections. More recently, the explosion of microcephaly among infants born to ZIKV-infected women has made ZIKV a global public health concern. While ZIKV causes acute human diseases, infections of vector mosquitoes are basically non-pathogenic, allowing persistent infections and conferring lifelong ability to transmit the virus. Recent studies have revealed that DNA forms of arboviral RNA genomes play a significant role in viral persistence in mosquitoes. We have initiated experiments to determine whether ZIKV generates viral DNA (vDNA) forms following infection in mosquitoes. Here we show that vDNAs are generated following ZIKV infection both in mosquito cell cultures and in its primary vector Aedes aegypti. vDNA formation is more extensive in RNA interference (RNAi)-deficient Aedes albopictus-derived C6/36 cells compared to RNAi-proficient mosquito cells. In addition, vDNAs are generated via multiple template-switching events.

Endogenous Retrovirus 3 - History, Physiology, and Pathology

Endogenous viral elements (EVE) seem to be present in all eukaryotic genomes. The composition of EVE varies between different species. The endogenous retrovirus 3 (ERV3) is one of these elements that is present only in humans and other Catarrhini. Conservation of ERV3 in most of the investigated Catarrhini and the expression pattern in normal tissues suggest a putative physiological role of ERV3. On the other hand, ERV3 has been implicated in the pathogenesis of auto-immunity and cancer. In the present review we summarize knowledge about this interesting EVE. We propose the model that expression of ERV3 (and probably other EVE loci) under pathological conditions might be part of a metazoan SOS response.

Deducing the Role of Virus Genome-Derived PIWI-Associated RNAs in the Mosquito-Arbovirus Arms Race

The P-element-induced wimpy testis (PIWI)-associated RNA (piRNA) pathway is known for its role in the protection of genome integrity in the germline of Drosophila melanogaster by silencing transposable elements. The piRNAs that target transposons originate from piRNA clusters in transposon-rich regions of the Drosophila genome and are processed by three PIWI family proteins. In Aedes aegypti and Aedes albopictus mosquitoes, which are two of the most important vectors of arthropod-borne viruses (arboviruses), the number of PIWI family genes has expanded and some are expressed in somatic, as well as germline, tissues. These discoveries have led to active research to explore the possible expanded functional roles of the piRNA pathway in vector mosquitoes. Virus genome-derived piRNAs (which will be referred to as (virus name) vpiRNAs) have been demonstrated in Aedes spp. cultured cells and mosquitoes after infection by arthropod-borne alpha-, bunya-, and flaviviruses. However, although Culex quinquefasciatus also is an important arbovirus vector and has an expansion of PIWI family genes, vpiRNAs have seldom been documented in this genus after arbovirus infection. Generation of complementary DNA (cDNA) fragments from RNA genomes of alpha-, bunya-, and flaviviruses (viral-derived cDNAs, vDNAs) has been demonstrated in cultured Aedes spp. cells and mosquitoes, and endogenous viral elements (EVEs), cDNA fragments of non-retroviral RNA virus genomes, are found more abundantly in genomes of Ae. aegypti and Ae. albopictus than other vector mosquitoes. These observations have led to speculation that vDNAs are integrated into vector genomes to form EVEs, which serve as templates for the transcription of antiviral vpiRNA precursors. However, no EVEs derived from alphavirus genomes have been demonstrated in genomes of any vector mosquito. In addition, although EVEs have been shown to be a source of piRNAs, the preponderance of EVEs described in Aedes spp. vectors are more closely related to the genomes of persistently infecting insect-specific viruses than to acutely infecting arboviruses. Furthermore, the signature patterns of the “ping-pong” amplification cycle that maintains transposon-targeting piRNAs in Drosophila are also evident in alphavirus and bunyavirus vpiRNAs, but not in vpiRNAs of flaviviruses. These divergent observations have rendered deciphering the mechanism(s) of biogenesis and potential role of vpiRNAs in the mosquito–arbovirus arms race difficult, and the focus of this review will be to assemble major findings regarding vpiRNAs and antiviral immunity in the important arbovirus vectors from Aedes and Culex genera.

Shrimp Parvovirus Circular DNA Fragments Arise From Both Endogenous Viral Elements and the Infecting Virus

Some insects use endogenous reverse transcriptase (RT) to make variable viral copy DNA (vcDNA) fragments from viral RNA in linear (lvcDNA) and circular (cvcDNA) forms. The latter form is easy to extract selectively. The vcDNA produces small interfering RNA (siRNA) variants that inhibit viral replication via the RNA interference (RNAi) pathway. The vcDNA is also autonomously inserted into the host genome as endogenous viral elements (EVE) that can also result in RNAi. We hypothesized that similar mechanisms occurred in shrimp. We used the insect methods to extract circular viral copy DNA (cvcDNA) from the giant tiger shrimp (Penaeus monodon) infected with a virus originally named infectious hypodermal and hematopoietic necrosis virus (IHHNV). Simultaneous injection of the extracted cvcDNA plus IHHNV into whiteleg shrimp (Penaeus vannamei) resulted in a significant reduction in IHHNV replication when compared to shrimp injected with IHHNV only. Next generation sequencing (NGS) revealed that the extract contained a mixture of two general IHHNV-cvcDNA types. One showed 98 to 99% sequence identity to GenBank record AF218266 from an extant type of infectious IHHNV. The other type showed 98% sequence identity to GenBank record DQ228358, an EVE formerly called non-infectious IHHNV. The startling discovery that EVE could also give rise to cvcDNA revealed that cvcDNA provided an easy means to identify and characterize EVE in shrimp and perhaps other organisms. These studies open the way for identification, characterization and use of protective cvcDNA as a potential shrimp vaccine and as a tool to identify, characterize and select naturally protective EVE to improve shrimp tolerance to homologous viruses in breeding programs.

Endogenous Caulimovirids: Fossils, Zombies, and Living in Plant Genomes

The Caulimoviridae is a family of double-stranded DNA viruses that infect plants. The genomes of most vascular plants contain endogenous caulimovirids (ECVs), a class of repetitive DNA elements that is abundant in some plant genomes, resulting from the integration of viral DNA in the chromosomes of germline cells during episodes of infection that have sometimes occurred millions of years ago. In this review, we reflect on 25 years of research on ECVs that has shown that members of the Caulimoviridae have occupied an unprecedented range of ecological niches over time and shed light on their diversity and macroevolution. We highlight gaps in knowledge and prospects of future research fueled by increased access to plant genome sequence data and new tools for genome annotation for addressing the extent, impact, and role of ECVs on plant biology and the origin and evolutionary trajectories of the Caulimoviridae.

Protein Structure-Guided Hidden Markov Models (HMMs) as A Powerful Method in the Detection of Ancestral Endogenous Viral Elements

It has been believed for a long time that the transfer and fixation of genetic material from RNA viruses to eukaryote genomes is very unlikely. However, during the last decade, there have been several cases in which “virus-to-host” gene transfer from various viral families into various eukaryotic phyla have been described. These transfers have been identified by sequence similarity, which may disappear very quickly, especially in the case of RNA viruses. However, compared to sequences, protein structure is known to be more conserved. Applying protein structure-guided protein domain-specific Hidden Markov Models, we detected homologues of the Virgaviridae capsid protein in Schizophora flies. Further data analysis supported “virus-to-host” transfer into Schizophora ancestors as a single transfer event. This transfer was not identifiable by BLAST or by other methods we applied. Our data show that structure-guided Hidden Markov Models should be used to detect ancestral virus-to-host transfers.

Modeling of Mutational Events in the Evolution of Viruses

Diverse studies of viral evolution have led to the recognition that the evolutionary rates of viral taxa observed are dependent on the time scale being investigated—with short-term studies giving fast substitution rates, and orders of magnitude lower rates for deep calibrations. Although each of these factors may contribute to this time dependent rate phenomenon, a more fundamental cause should be considered. We sought to test computationally whether the basic phenomena of virus evolution (mutation, replication, and selection) can explain the relationships between the evolutionary and phylogenetic distances. We tested, by computational inference, the hypothesis that the phylogenetic distances between the pairs of sequences are functions of the evolutionary path lengths between them. A Basic simulation revealed that the relationship between simulated genetic and mutational distances is non-linear, and can be consistent with different rates of nucleotide substitution at different depths of branches in phylogenetic trees.

Molecular fossils reveal ancient associations of dsDNA viruses with several phyla of fungi

Little is known about the infections of double-stranded DNA (dsDNA) viruses in fungi. Here, we use a paleovirological method to systematically identify the footprints of past dsDNA virus infections within the fungal genomes. We uncover two distinct groups of endogenous nucleocytoplasmic large DNA viruses (NCLDVs) in at least seven fungal phyla (accounting for about a third of known fungal phyla), revealing an unprecedented diversity of dsDNA viruses in fungi. Interestingly, one fungal dsDNA virus lineage infecting six fungal phyla is closely related to the giant virus Pithovirus, suggesting giant virus relatives might widely infect fungi. Co-speciation analyses indicate fungal NCLDVs mainly evolved through cross-species transmission. Taken together, our findings provide novel insights into the diversity and evolution of NCLDVs in fungi.

First Evidence of Past and Present Interactions between Viruses and the Black Soldier Fly, Hermetia illucens

Black soldier flies (BSFs, Hermetia illucens) are becoming a prominent research model encouraged by the insect as food and feed and waste bioconversion industries. Insect mass-rearing facilities are at risk from the spread of viruses, but so far, none have been described in BSFs. To fill this knowledge gap, a bioinformatic approach was undertaken to discover viruses specifically associated with BSFs. First, BSF genomes were screened for the presence of endogenous viral elements (EVEs). This led to the discovery and mapping of seven orthologous EVEs integrated into three BSF genomes originating from five viral families. Secondly, a virus discovery pipeline was used to screen BSF transcriptomes. This led to detecting a new exogenous totivirus that we named hermetia illucens totivirus 1 (HiTV1). Phylogenetic analyses showed this virus belongs to a clade of insect-specific totiviruses and is closely related to the largest EVE located on chromosome 1 of the BSF genome. Lastly, this EVE was found to express a small transcript in some BSFs infected by HiTV1. Altogether, this data mining study showed that far from being unscathed from viruses, BSFs bear traces of past interactions with several viral families and of present interactions with the exogenous HiTV1.

Current trends in RNA virus detection through metatranscriptome sequencing data

With advances in sequencing technology, metatranscriptome sequencing from a variety of environmental and biological sources has revealed the existence of various previously unknown RNA viruses. This review presents recent major RNA virome studies sampled from invertebrate and vertebrate species as well as aquatic environments. In particular, we focus on severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) and related RNA virus identification through metatranscriptome sequencing analyses. Recently developed bioinformatics software and databases for RNA virus identification are introduced. A relationship between newly identified RNA viruses and endogenous viral elements in host genomes is also discussed.

The Metagenomic Analysis of Viral Diversity in Colorado Potato Beetle Public NGS Data

The Colorado potato beetle (CPB) is one of the most serious insect pests due to its high ecological plasticity and ability to rapidly develop resistance to insecticides. The use of biological insecticides based on viruses is a promising approach to control insect pests, but the information on viruses which infect leaf feeding beetles is scarce. We performed a metagenomic analysis of 297 CPB genomic and transcriptomic samples from the public National Center for Biotechnology Information Sequence Read Archive (NCBI SRA) database. The reads that were not aligned to the reference genome were assembled with metaSPAdes, and 13314 selected contigs were analyzed with BLAST tools. The contigs and non-aligned reads were also analyzed with Kraken2 software. A total of 3137 virus-positive contigs were attributed to different viruses belonging to 6 types, 17 orders, and 32 families, matching over 97 viral species. The annotated sequences can be divided into several groups: those that are homologous to genetic sequences of insect viruses (Adintoviridae, Ascoviridae, Baculoviridae, Dicistroviridae, Chuviridae, Hytrosaviridae, Iflaviridae, Iridoviridae, Nimaviridae, Nudiviridae, Phasmaviridae, Picornaviridae, Polydnaviriformidae, Xinmoviridae etc.), plant viruses (Betaflexiviridae, Bromoviridae, Kitaviridae, Potyviridae), and endogenous retroviral elements (Retroviridae, Metaviridae). Additionally, the full-length genomes and near-full length genome sequences of several viruses were assembled. We also found sequences belonging to Bracoviriform viruses and, for the first time, experimentally validated the presence of bracoviral genetic fragments in the CPB genome. Our work represents the first attempt to discover the viral genetic material in CPB samples, and we hope that further studies will help to identify new viruses to extend the arsenal of biopesticides against CPB.