Annual (2023) taxonomic update of RNA-directed RNA polymerase-encoding negative-sense RNA viruses (realm Riboviria: kingdom Orthornavirae: phylum Negarnaviricota)

In April 2023, following the annual International Committee on Taxonomy of Viruses (ICTV) ratification vote on newly proposed taxa, the phylum Negarnaviricota was amended and emended. The phylum was expanded by one new family, 14 new genera, and 140 new species. Two genera and 538 species were renamed. One species was moved, and four were abolished. This article presents the updated taxonomy of Negarnaviricota as now accepted by the ICTV.

The toxins of vertically transmitted Spiroplasma

Vertically transmitted (VT) microbial symbionts play a vital role in the evolution of their insect hosts. A longstanding question in symbiont research is what genes help promote long-term stability of vertically transmitted lifestyles. Symbiont success in insect hosts is due in part to expression of beneficial or manipulative phenotypes that favor symbiont persistence in host populations. In Spiroplasma, these phenotypes have been linked to toxin and virulence domains among a few related strains. However, these domains also appear frequently in phylogenetically distant Spiroplasma, and little is known about their distribution across the Spiroplasma genus. In this study, we present the complete genome sequence of the Spiroplasma symbiont of Drosophila atripex, a non-manipulating member of the Ixodetis clade of Spiroplasma, for which genomic data are still limited. We perform a genus-wide comparative analysis of toxin domains implicated in defensive and reproductive phenotypes. From 12 VT and 31 non-VT Spiroplasma genomes, ribosome-inactivating proteins (RIPs), OTU-like cysteine proteases (OTUs), ankyrins, and ETX/MTX2 domains show high propensity for VT Spiroplasma compared to non-VT Spiroplasma. Specifically, OTU and ankyrin domains can be found only in VT-Spiroplasma, and RIP domains are found in all VT Spiroplasma and three non-VT Spiroplasma. These domains are frequently associated with Spiroplasma plasmids, suggesting a possible mechanism for dispersal and maintenance among heritable strains. Searching insect genome assemblies available on public databases uncovered uncharacterized Spiroplasma genomes from which we identified several spaid-like genes encoding RIP, OTU, and ankyrin domains, suggesting functional interactions among those domain types. Our results suggest a conserved core of symbiont domains play an important role in the evolution and persistence of VT Spiroplasma in insects.

2022 taxonomic update of phylum Negarnaviricota (Riboviria: Orthornavirae), including the large orders Bunyavirales and Mononegavirales

In March 2022, following the annual International Committee on Taxonomy of Viruses (ICTV) ratification vote on newly proposed taxa, the phylum Negarnaviricota was amended and emended. The phylum was expanded by two new families (bunyaviral Discoviridae and Tulasviridae), 41 new genera, and 98 new species. Three hundred forty-nine species were renamed and/or moved. The accidentally misspelled names of seven species were corrected. This article presents the updated taxonomy of Negarnaviricota as now accepted by the ICTV.

Highly transmissible cytoplasmic incompatibility by the extracellular insect symbiont Spiroplasma

Cytoplasmic incompatibility (CI) is a form of reproductive manipulation caused by maternally inherited endosymbionts infecting arthropods, like Wolbachia, whereby matings between infected males and uninfected females produce few or no offspring. We report the discovery of a new CI symbiont, a strain of Spiroplasma causing CI in the parasitoid wasp Lariophagus distinguendus. Its extracellular occurrence enabled us to establish CI in uninfected adult insects by transferring Spiroplasma-infected hemolymph. We sequenced the CI-Spiroplasma genome and did not find any homologues of any of the cif genes discovered to cause CI in Wolbachia, suggesting independent evolution of CI. Instead, the genome contains other potential CI-causing candidate genes, such as homologues of high-mobility group (HMG) box proteins that are crucial in eukaryotic development but rare in bacterial genomes. Spiroplasma’s extracellular nature and broad host range encompassing medically and agriculturally important arthropods make it a promising tool to study CI and its applications.

•

A Spiroplasma strain induces cytoplasmic incompatibility (CI) in a parasitoid wasp

•

Spiroplasma and CI are efficiently transferred to naive hosts by hemolymph injection

•

The CI-Spiroplasma genome contains potential CI-causing genes, but no Wolbachia cifs

•

Absence of cif genes in the Spiroplasma genome suggests an independent origin of CI

Evolution and Diversity of Inherited Viruses in the Nearctic Phantom Midge, Chaoborus americanus

Inherited mutualists, parasites, and commensals occupy one of the most intimate ecological niches available to invertebrate-associated microbes. How this transmission environment influences microbial evolution is increasingly understood for inherited bacterial symbionts, but in viruses, research on the prevalence of vertical transmission and its effects on viral lineages is still maturing. The evolutionary stability of this strategy remains difficult to assess, although phylogenetic evidence of frequent host shifts and selective sweeps have been interpreted as strategies favoring parasite persistence. In this study, we describe and investigate a natural insect system in which species-wide sweeps have been restricted by the isolation of host populations. Previous work identified evidence of pronounced mitochondrial genetic structure among North American populations of the phantom midge, Chaoborus americanus. Here we take advantage of the geographical isolation in this species to investigate the diversity and persistence of its inherited virome. We identify eight novel RNA viruses from six families and use small RNA sequencing in reproductive tissues to provide evidence of vertical transmission. We report region-specific virus strains that mirror the continental phylogeography of the host, demonstrating that members of the inherited virome have independently persisted in parallel host lineages since they last shared a common ancestor in the Mid Pleistocene. We find that the small interfering RNA pathway, a frontline of antiviral defense in insects, targets members of this inherited virome. Finally, our results suggest that the Piwi-mediated RNA silencing pathway is unlikely to function as a general antiviral defense in Chaoborus, in contrast to its role in some mosquitoes. However, we also report that the PIWI-interacting RNA pathway generates abundant piRNAs from endogenous viral elements closely related to actively infecting inherited viruses, potentially helping to explain idiosyncratic patterns of virus-specific Piwi targeting in this insect.

2021 Taxonomic update of phylum Negarnaviricota (Riboviria: Orthornavirae), including the large orders Bunyavirales and Mononegavirales

In March 2021, following the annual International Committee on Taxonomy of Viruses (ICTV) ratification vote on newly proposed taxa, the phylum Negarnaviricota was amended and emended. The phylum was expanded by four families (Aliusviridae, Crepuscuviridae, Myriaviridae, and Natareviridae), three subfamilies (Alpharhabdovirinae, Betarhabdovirinae, and Gammarhabdovirinae), 42 genera, and 200 species. Thirty-nine species were renamed and/or moved and seven species were abolished. This article presents the updated taxonomy of Negarnaviricota as now accepted by the ICTV.

2020 taxonomic update for phylum Negarnaviricota (Riboviria: Orthornavirae), including the large orders Bunyavirales and Mononegavirales

In March 2020, following the annual International Committee on Taxonomy of Viruses (ICTV) ratification vote on newly proposed taxa, the phylum Negarnaviricota was amended and emended. At the genus rank, 20 new genera were added, two were deleted, one was moved, and three were renamed. At the species rank, 160 species were added, four were deleted, ten were moved and renamed, and 30 species were renamed. This article presents the updated taxonomy of Negarnaviricota as now accepted by the ICTV.

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.

Taxonomy of the order Bunyavirales: update 2019

In February 2019, following the annual taxon ratification vote, the order Bunyavirales was amended by creation of two new families, four new subfamilies, 11 new genera and 77 new species, merging of two species, and deletion of one species. This article presents the updated taxonomy of the order Bunyavirales now accepted by the International Committee on Taxonomy of Viruses (ICTV).

The defensive Spiroplasma

Defensive microbes are of great interest for their roles in arthropod health, disease transmission, and biocontrol efforts. Obligate bacterial passengers of arthropods, such as Spiroplasma, confer protection against the natural enemies of their hosts to improve their own fitness. Although known for less than a decade, Spiroplasma's defensive reach extends to diverse parasites, both microbial and multicellular. We provide an overview of known defensive phenotypes against nematodes, parasitoid wasps, and fungi, and highlight recent studies supporting the role of Spiroplasma-encoded ribosome-inactivating proteins in protection. With cellular features well-suited for life in the hemolymph, broad distribution among invertebrate hosts, and the capacity to repeatedly evolve vertical transmission, Spiroplasma may be uniquely equipped to form intimate, defensive associations to combat extracellular parasites. Along with insights into defensive mechanisms, recent significant advances have been made in male-killing - a phenotype with interesting evolutionary ties to defense. Finally, we look forward to an exciting decade using the genetic tools of Drosophila, and the rapidly-advancing tractability of Spiroplasma itself, to better understand mechanisms and evolution in defensive symbiosis.

Taxonomy of the order Bunyavirales: second update 2018

In October 2018, the order Bunyavirales was amended by inclusion of the family Arenaviridae, abolishment of three families, creation of three new families, 19 new genera, and 14 new species, and renaming of three genera and 22 species. This article presents the updated taxonomy of the order Bunyavirales as now accepted by the International Committee on Taxonomy of Viruses (ICTV).

Toxin and Genome Evolution in a Drosophila Defensive Symbiosis

Defenses conferred by microbial symbionts play a vital role in the health and fitness of their animal hosts. An important outstanding question in the study of defensive symbiosis is what determines long term stability and effectiveness against diverse natural enemies. In this study, we combine genome and transcriptome sequencing, symbiont transfection and parasite protection experiments, and toxin activity assays to examine the evolution of the defensive symbiosis between Drosophila flies and their vertically transmitted Spiroplasma bacterial symbionts, focusing in particular on ribosome-inactivating proteins (RIPs), symbiont-encoded toxins that have been implicated in protection against both parasitic wasps and nematodes. Although many strains of Spiroplasma, including the male-killing symbiont (sMel) of Drosophila melanogaster, protect against parasitic wasps, only the strain (sNeo) that infects the mycophagous fly Drosophila neotestacea appears to protect against parasitic nematodes. We find that RIP repertoire is a major differentiating factor between strains that do and do not offer nematode protection, and that sMel RIPs do not show activity against nematode ribosomes in vivo. We also discovered a strain of Spiroplasma infecting a mycophagous phorid fly, Megaselia nigra. Although both the host and its Spiroplasma are distantly related to D. neotestacea and its symbiont, genome sequencing revealed that the M. nigra symbiont encodes abundant and diverse RIPs, including plasmid-encoded toxins that are closely related to the RIPs in sNeo. Our results suggest that distantly related Spiroplasma RIP toxins may perform specialized functions with regard to parasite specificity and suggest an important role for horizontal gene transfer in the emergence of novel defensive phenotypes.

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.

An exceptional family: Ophiocordyceps-allied fungus dominates the microbiome of soft scale insects (Hemiptera: Sternorrhyncha: Coccidae)

Hemipteran insects of the suborder Sternorrhyncha are plant sap feeders, where each family is obligately associated with a specific bacterial endosymbiont that produces essential nutrients lacking in the sap. Coccidae (soft scale insects) is the only major sternorrhynchan family in which obligate symbiont(s) have not been identified. We studied the microbiota in seven species from this family from Israel, Spain and Cyprus, by high-throughput sequencing of ribosomal genes, and found that no specific bacterium was prevalent and abundant in all the tested species. In contrast, an Ophiocordyceps-allied fungus sp.-a lineage widely known as entomopathogenic-was highly prevalent. All individuals of all the tested species carried this fungus. Phylogenetic analyses showed that the Ophiocordyceps-allied fungus from the coccids is closely related to fungi described from other hemipterans, and they appear to be monophyletic, although the phylogenies of the Ophiocordyceps-allied fungi and their hosts do not appear to be congruent. Microscopic observations show that the fungal cells are lemon-shaped, are distributed throughout the host's body and are present in the eggs, suggesting vertical transmission. Taken together, the results suggest that the Ophiocordyceps-allied fungus may be a primary symbiont of Coccidae-a major evolutionary shift from bacteria to fungi in the Sternorrhyncha, and an important example of fungal evolutionary lifestyle switch.

Generality of toxins in defensive symbiosis: Ribosome-inactivating proteins and defense against parasitic wasps in Drosophila

While it has become increasingly clear that multicellular organisms often harbor microbial symbionts that protect their hosts against natural enemies, the mechanistic underpinnings underlying most defensive symbioses are largely unknown. Spiroplasma bacteria are widespread associates of terrestrial arthropods, and include strains that protect diverse Drosophila flies against parasitic wasps and nematodes. Recent work implicated a ribosome-inactivating protein (RIP) encoded by Spiroplasma, and related to Shiga-like toxins in enterohemorrhagic Escherichia coli, in defense against a virulent parasitic nematode in the woodland fly, Drosophila neotestacea. Here we test the generality of RIP-mediated protection by examining whether Spiroplasma RIPs also play a role in wasp protection, in D. melanogaster and D. neotestacea. We find strong evidence for a major role of RIPs, with ribosomal RNA (rRNA) from the larval endoparasitic wasps, Leptopilina heterotoma and Leptopilina boulardi, exhibiting the hallmarks of RIP activity. In Spiroplasma-containing hosts, parasitic wasp ribosomes show abundant site-specific depurination in the α-sarcin/ricin loop of the 28S rRNA, with depurination occurring soon after wasp eggs hatch inside fly larvae. Interestingly, we found that the pupal ectoparasitic wasp, Pachycrepoideus vindemmiae, escapes protection by Spiroplasma, and its ribosomes do not show high levels of depurination. We also show that fly ribosomes show little evidence of targeting by RIPs. Finally, we find that the genome of D. neotestacea’s defensive Spiroplasma encodes a diverse repertoire of RIP genes, which are differ in abundance. This work suggests that specificity of defensive symbionts against different natural enemies may be driven by the evolution of toxin repertoires, and that toxin diversity may play a role in shaping host-symbiont-enemy interactions.

Nearly all insects harbor bacterial partners. These microbes can be defensive symbionts, protecting their hosts against parasites and pathogens, and in some cases, may defend against more than one enemy, presenting opportunity to study the evolution of specificity underlying defensive interactions. What factors determine specificity and generality of defense? Can symbiont-encoded effector molecules act generally against enemies without causing harm to the host? We show here that symbiont-encoded ribosome-inactivating toxins, previously implicated in protection of a Drosophila fruit fly against its nematode parasite, are also implicated in defending flies against parasitic wasps. We quantify activity of the toxin as the proportion of ribosomes depurinated and, importantly, find that susceptible wasps are attacked early in development and are strongly affected, while hosts and a resistant wasp are not. We also show that not one, but a family of toxins is maintained and expressed by symbionts. Together, our findings implicate toxin diversity as a factor contributing to the evolution of specificity in a symbiont-mediated defense.

Unexpected differences in the population genetics of phasmavirids (Bunyavirales) from subarctic ponds

Little is known of the evolution of RNA viruses in aquatic systems. Here, we assess the genetic connectivity of two bunyaviruses (Kigluaik phantom orthophasmavirus or KIGV and Nome phantom orthophasmavirus or NOMV) with zooplanktonic hosts from subarctic ponds. We expected weak genetic structure among populations as the hosts (phantom midges) have a terrestrial winged dispersal stage. To test whether their respective viruses mirror this structure, we collected and analyzed population datasets from 21 subarctic freshwater ponds and obtained sequences from all four genes in the viral genomes. Prevalence averaged 66 per cent for 514 host specimens and was not significantly different between recently formed thaw ponds and glacial ponds. Unexpectedly, KIGV from older ponds showed pronounced haplotype divergence with little evidence of genetic connectivity. However, KIGV populations from recent thaw ponds appeared to be represented by a closely related haplotype group, perhaps indicating a genotypic dispersal bias. Unlike KIGV, NOMV had modest structure and diversity in recently formed thaw ponds. For each virus, we found elevated genetic diversity relative to the host, but similar population structures to the host. Our results suggest that non-random processes such as virus-host interactions, genotypic bias, and habitat effects differ among polar aquatic RNA viruses.

Bacterial associates of Hyalesthes obsoletus (Hemiptera: Cixiidae), the insect vector of bois noir disease, with a focus on cultivable bacteria

The planthopper Hyalesthes obsoletus (Hemiptera: Cixiidae) is an important vector of phytoplasma diseases in grapevine. In the current study, the bacterial community compositions of symbionts of this insect were examined. Two dominant bacterial lineages were identified by mass sequencing: the obligate symbiont Candidatus Sulcia, and a facultative symbiont that is closely related to Pectobacterium sp. and to BEV, a cultivable symbiont of another phytoplasma vector, the leafhopper Euscelidius variegatus. In addition, one bacterium was successfully isolated in this study - a member of the family Xanthomonadaceae that is most closely related to the genus Dyella. This Dyella-like bacterium (DLB) was detected by FISH analysis in H. obsoletus guts but not ovaries, and its prevalence in H. obsoletus increased during the fall, suggesting that it was acquired by the host through feeding. We found that DLB inhibits Spiroplasma melliferum, a cultivable relative of phytoplasma, suggesting that it is a potential candidate for biological control against phytoplasma in grapevines.

Evidence that ebolaviruses and cuevaviruses have been diverging from marburgviruses since the Miocene

An understanding of the timescale of evolution is critical for comparative virology but remains elusive for many RNA viruses. Age estimates based on mutation rates can severely underestimate divergences for ancient viral genes that are evolving under strong purifying selection. Paleoviral dating, however, can provide minimum age estimates for ancient divergence, but few orthologous paleoviruses are known within clades of extant viruses. For example, ebolaviruses and marburgviruses are well-studied mammalian pathogens, but their comparative biology is difficult to interpret because the existing estimates of divergence are controversial. Here we provide evidence that paleoviral elements of two genes (ebolavirus-like VP35 and NP) in cricetid rodent genomes originated after the divergence of ebolaviruses and cuevaviruses from marburgviruses. We provide evidence of orthology by identifying common paleoviral insertion sites among the rodent genomes. Our findings indicate that ebolaviruses and cuevaviruses have been diverging from marburgviruses since the early Miocene.

Selectively maintained paleoviruses in Holarctic water fleas reveal an ancient origin for phleboviruses

The ecological model, Daphnia pulex (Cladocera: Daphniidae), is broadly distributed in Holarctic freshwater habitats and has been the subject of multidisciplinary study for over half a century, but never has a natural RNA virus infection been reported in daphnids. Here we report on a group of paleoviruses related to RNA dependent RNA polymerase in the genome of D. pulex. Phylogenetic analysis suggests that these paleoviruses are derived from a viral lineage within the genus Phlebovirus. Comparison of the genomic sequences flanking individual paleoviruses reveal that some are orthologous viral insertions having been present in the common ancestor of the D. pulex species complex, which is millions of years old. Still, we detected some sites that have the signature of purifying selection. In contrast, other paleoviruses in this group seem to be unique to specific host lineages and even contain undisrupted open reading frames, suggesting either more recent acquisition, or selective maintenance.

Virus-host co-evolution under a modified nuclear genetic code

Among eukaryotes with modified nuclear genetic codes, viruses are unknown. However, here we provide evidence of an RNA virus that infects a fungal host (Scheffersomyces segobiensis) with a derived nuclear genetic code where CUG codes for serine. The genomic architecture and phylogeny are consistent with infection by a double-stranded RNA virus of the genus Totivirus. We provide evidence of past or present infection with totiviruses in five species of yeasts with modified genetic codes. All but one of the CUG codons in the viral genome have been eliminated, suggesting that avoidance of the modified codon was important to viral adaptation. Our mass spectroscopy analysis indicates that a congener of the host species has co-opted and expresses a capsid gene from totiviruses as a cellular protein. Viral avoidance of the host's modified codon and host co-option of a protein from totiviruses suggest that RNA viruses co-evolved with yeasts that underwent a major evolutionary transition from the standard genetic code.

Phylogeny, integration and expression of sigma virus-like genes in Drosophila

The recent and surprising discovery of widespread NIRVs (non-retroviral integrated RNA viruses) has highlighted the importance of genomic interactions between non-retroviral RNA viruses and their eukaryotic hosts. Among the viruses with integrated representatives are the rhabdoviruses, a family of negative sense single-stranded RNA viruses. We identify sigma virus-like NIRVs of Drosophila spp. that represent unique cases where NIRVs are closely related to exogenous RNA viruses in a model host organism. We have used a combination of bioinformatics and laboratory methods to explore the evolution and expression of sigma virus-like NIRVs in Drosophila. Recent integrations in Drosophila provide a promising experimental system to study functionality of NIRVs. Moreover, the genomic architecture of recent NIRVs provides an unusual evolutionary window on the integration mechanism. For example, we found that a sigma virus-like polymerase associated protein (P) gene appears to have been integrated by template switching of the blastopia-like LTR retrotransposon. The sigma virus P-like NIRV is present in multiple retroelement fused open reading frames on the X and 3R chromosomes of Drosophila yakuba - the X-linked copy is transcribed to produce an RNA product in adult flies. We present the first account of sigma virus-like NIRVs and the first example of NIRV expression in a model animal system, and therefore provide a platform for further study of the possible functions of NIRVs in animal hosts.

Evolutionary maintenance of filovirus-like genes in bat genomes

BACKGROUND

Little is known of the biological significance and evolutionary maintenance of integrated non-retroviral RNA virus genes in eukaryotic host genomes. Here, we isolated novel filovirus-like genes from bat genomes and tested for evolutionary maintenance. We also estimated the age of filovirus VP35-like gene integrations and tested the phylogenetic hypotheses that there is a eutherian mammal clade and a marsupial/ebolavirus/Marburgvirus dichotomy for filoviruses.

RESULTS

We detected homologous copies of VP35-like and NP-like gene integrations in both Old World and New World species of Myotis (bats). We also detected previously unknown VP35-like genes in rodents that are positionally homologous. Comprehensive phylogenetic estimates for filovirus NP-like and VP35-like loci support two main clades with a marsupial and a rodent grouping within the ebolavirus/Lloviu virus/Marburgvirus clade. The concordance of VP35-like, NP-like and mitochondrial gene trees with the expected species tree supports the notion that the copies we examined are orthologs that predate the global spread and radiation of the genus Myotis. Parametric simulations were consistent with selective maintenance for the open reading frame (ORF) of VP35-like genes in Myotis. The ORF of the filovirus-like VP35 gene has been maintained in bat genomes for an estimated 13. 4 MY. ORFs were disrupted for the NP-like genes in Myotis. Likelihood ratio tests revealed that a model that accommodates positive selection is a significantly better fit to the data than a model that does not allow for positive selection for VP35-like sequences. Moreover, site-by-site analysis of selection using two methods indicated at least 25 sites in the VP35-like alignment are under positive selection in Myotis.

CONCLUSIONS

Our results indicate that filovirus-like elements have significance beyond genomic imprints of prior infection. That is, there appears to be, or have been, functionally maintained copies of such genes in mammals. "Living fossils" of filoviruses appear to be selectively maintained in a diverse mammalian genus (Myotis).