Endogenous retroviruses provide important insights into the deep history of this viral lineage. Endogenous foamy viruses are thought to be very rare and only a few cases have been identified to date. Here we report a novel endogenous foamy virus (CaEFV) within the genome of the Cape golden mole (Chrysochloris asiatica). The identification of CaEFV reveals the presence of foamy virus in the placental mammal superorder Afrotheria. Phylogenetic analyses place CaEFV basal to other foamy viruses of Eutherian origin, suggesting an ancient codivergence between foamy virus and placental mammals. These findings have implications for understanding the long-term evolution, diversity, and biology of retroviruses.

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.

Both the replication and transcription of the influenza virus are catalyzed by the viral polymerase complex. The polymerases of most avian influenza A viruses have poor performance in mammalian cells, which is considered to be one of the important species barriers. Pigs have been long considered as important intermediate hosts for interspecies transmission of the avian influenza virus, because of their susceptibility to infection with both avian and mammalian influenza viruses. However, the molecular basis of influenza polymerase adaptation in pigs remains largely unknown. ANP32A and ANP32B proteins have been identified as playing fundamental roles in influenza virus replication and host range determination. In this study, we found that swine ANP32A (swANP32A), unlike swine ANP32B or other mammalian ANP32A or B, shows stronger supporting activity to avian viral polymerase. Knockout of ANP32A in pig cells PK15 dramatically reduced avian influenza polymerase activity and viral infectivity, suggesting a unique feature of swANP32A in supporting avian influenza viral polymerase. This species-specific activity is mapped to two key sites, 106V and 156S, in swANP32A. Interestingly, the amino acid 106V is unique to pigs among all the vertebrate species studied, and when combined with 156S, exhibits positive epistasis in pigs. Mutation of 106V and 156S to the signature found in ANP32As from other mammalian species weakened the interaction between swANP32A and chicken viral polymerase, and reduced polymerase activity. Understanding the molecular basis of ANP32 proteins may help to discover new antiviral targets and design avian influenza resistant genome edited pigs.

Eukaryotes have occasionally acquired genetic material through horizontal gene transfer (HGT). However, little is known about the evolutionary and functional significance of such acquisitions. Lysozymes are ubiquitous enzymes that degrade bacterial cell walls. Here, we provide evidence that two subclasses of bivalves (Heterodonta and Palaeoheterodonta) acquired a lysozyme gene via HGT, building on earlier findings. Phylogenetic analyses place the bivalve lysozyme genes within the clade of bacteriophage lysozyme genes, indicating that the bivalves acquired the phage-type lysozyme genes from bacteriophages, either directly or through intermediate hosts. These bivalve lysozyme genes underwent dramatic structural changes after their co-option, including intron gain and fusion with other genes. Moreover, evidence suggests that recurrent gene duplication occurred in the bivalve lysozyme genes. Finally, we show the co-opted lysozymes exhibit a capacity for antibacterial action, potentially augmenting the immune function of related bivalves. This represents an intriguing evolutionary strategy in the eukaryote-microbe arms race, in which the genetic materials of bacteriophages are co-opted by eukaryotes, and then used by eukaryotes to combat bacteria, using a shared weapon against a common enemy.

The deep history and early diversification of retroviruses remains elusive, largely because few retroviruses have been characterized in vertebrates other than mammals and birds. Endogenous retroviruses (ERVs) documented past retroviral infections and thus provide 'molecular fossils' for studying the deep history of retroviruses. Here we perform a comprehensive phylogenomic analysis of ERVs within the genomes of 92 non-avian/mammalian vertebrates, including 72 fishes, 4 amphibians, and 16 reptiles. We find that ERVs are present in all the genomes of jawed vertebrates, revealing the ubiquitous presence of ERVs in jawed vertebrates. We identify a total of >8,000 ERVs and reconstruct ~450 complete or partial ERV genomes, which dramatically expands the phylogenetic diversity of retroviruses and suggests that the diversity of exogenous retroviruses might be much underestimated in non-avian/mammalian vertebrates. Phylogenetic analyses show that retroviruses cluster into five major groups with different host distributions, providing important insights into the classification and diversification of retroviruses. Moreover, we find retroviruses mainly underwent frequent host switches in non-avian/mammalian vertebrates, with exception of spumavirus-related viruses that codiverged with their ray-finned fish hosts. Interestingly, ray-finned fishes and turtles appear to serve as unappreciated hubs for the transmission of retroviruses. Finally, we find retroviruses underwent many independent water-land transmissions, indicating the water-land interface is not a strict barrier for retrovirus transmission. Our analyses provide unprecedented insights into and valuable resources for studying the diversification, key evolutionary transitions, and macroevolution of retroviruses.

Polydnaviruses are obligate mutualists of parasitoid wasps and are divided into two genera, Bracovirus and Ichnovirus. Bracoviruses are thought to originate from a single integration of an ancestral nudivirus into the ancestor of microgastroid complex ~100 million years ago. However, all the known nudiviruses are only distantly related to bracoviruses, and much remains obscure about the origin of bracoviruses. Here we employ a paleovirological method to screen endogenous nudivirus-like elements across arthropods. Interestingly, we identify many endogenous nudivirus-like elements within the genome of Eurytoma brunniventris, a species of the Chalcidoidea superfamily. Among them, we find 14 core gene sequences are likely to be derived from a betanudivirus (designated EbrENV-β), suggesting that betanudivirus has been circulating in parasitoid wasps. Phylogenomic analysis suggests that EbrENV-β is the known closest relative of bracoviruses. Synteny analyses show the order of core genes is not well conserved between EbrENV-β and nudiviruses, revealing the dynamic nature of the evolution of nudivirus genome structures. Our findings narrow down the evolutionary gap between bracoviruses and nudiviruses and provide novel insights into the origin and evolution of polydnaviruses.

Our knowledge of citrus viruses is largely skewed toward virus pathology in cultivated orchards. Little is known about the virus diversity in wild citrus species. Here, we used a metatranscriptomics approach to characterize the virus diversity in a wild citrus habitat within the proposed center of the origin of citrus plants. We discovered a total of 44 virus isolates that could be classified into species Citrus tristeza virus and putative species citrus associated ampelovirus 1, citrus associated ampelovirus 2, and citrus virus B within the family Closteroviridae, providing important information to explore the factors facilitating outbreaks of citrus viruses and the evolutionary history of the family Closteroviridae. We found that frequent horizontal gene transfer, gene duplication, and alteration of expression strategy have shaped the genome complexity and diversification of the family Closteroviridae. Recombination frequently occurred among distinct Closteroviridae members, thereby facilitating the evolution of Closteroviridae. Given the potential emergence of similar wild-citrus-originated novel viruses as pathogens, the need for surveillance of their pathogenic and epidemiological characteristics is of utmost priority for global citrus production.

Transposable elements (TEs) comprise a large proportion of the eukaryote genomes. Yet it remains poorly understood how TEs influence the fitness of the hosts carrying them. Here, we empirically test the impact of TEs on the host fitness in the fission yeast Schizosaccharomyces pombe. We find that two families of TEs (Tf1 and Tf2 elements), both of which belong to long terminal repeat retrotransposons, are highly polymorphic among individual S. pombe strains. Only 13 complete Tf2 elements are identified in S. pombe laboratory strain 972. These 13 Tf2 elements integrated into host genomes in very recent time and are segregating within the S. pombe population. Through knocking out each of the 13 Tf2 elements in S. pombe strain 972, we find Tf2 knockout does not affect the host fitness, and Tf2 elements do not alter the expression of nearby genes. Challenged by diverse forms of stress, the Tf2 knockout strains do not exhibit different growth rates from wild-type strain. Together, we conclude that segregating complete Tf2 elements insertions are largely neutral to host fitness in the fission yeast. Our study provides genome-wide empirical support for the selfish nature of TEs in fission yeast.

The origin and deep history of retroviruses remain mysterious and contentious, largely because the diversity of retroviruses is incompletely understood. Here, we report the discovery of lokiretroviruses, a novel major lineage of retroviruses, within the genomes of a wide range of vertebrates (at least 137 species), including lampreys, ray-finned fishes, lobe-finned fishes, amphibians, and reptiles. Lokiretroviruses share a similar genome architecture with known retroviruses, but display some unique features. Interestingly, lokiretrovirus Env proteins share detectable similarity with fusion glycoproteins of viruses within the Mononegavirales order, blurring the boundary between retroviruses and negative sense single-stranded RNA viruses. Phylogenetic analyses based on reverse transcriptase demonstrate that lokiretroviruses are sister to all the retroviruses sampled to date, providing a crucial nexus for studying the deep history of retroviruses. Comparing congruence between host and virus phylogenies suggests lokiretroviruses mainly underwent cross-species transmission. Moreover, we find that retroviruses replaced their ribonuclease H and integrase domains multiple times during their evolutionary course, revealing the importance of domain shuffling in the evolution of retroviruses. Overall, our findings greatly expand our views of the diversity of retroviruses, and provide novel insights into the origin and complex evolutionary history of retroviruses.

The origin and deep evolution of retroviruses remain largely unclear. It has been proposed that retroviruses might have originated from a Ty3/Gypsy retrotransposon, but all known Ty3/Gypsy retrotransposons are only distantly related to retroviruses. Retroviruses and some plant Athila/Tat elements (within Ty3/Gypsy retrotransposons) independently evolved a dual RNase H domain and an env/env-like gene. Here, we reported the discovery of a novel lineage of retrotransposons, designated Odin retrotransposons, in the genomes of eight sea anemones (order Actinaria) within the Cnidaria phylum. Odin retrotransposons exhibited unique genome features, encoding a dual RNase H domain (like retroviruses) but no env gene (like most Ty3/Gypsy retrotransposons). Phylogenetic analyses based on reverse transcriptase showed that Odin retrotransposons formed a sister group to lokiretroviruses, and lokiretroviruses and Odin retrotransposons together were sister to canonical retroviruses. Moreover, phylogenetic analyses based on RNase H and integrase also supported the hypothesis that Odin retrotransposons were sisters to lokiretroviruses. Lokiretroviruses and canonical retroviruses did not form a monophyletic group, indicating that lokiretroviruses and canonical retroviruses might represent two distinct virus families. Taken together, the discovery of Odin retrotransposons narrowed down the evolutionary gaps between retrotransposons and canonical retroviruses and lokiretroviruses. IMPORTANCE The origin of retroviruses remains largely unclear. In this study, we discovered a novel retrotransposon lineage, Odin retrotransposons, within the genomes of sea anemones (order Actinaria). In contrast to retroviruses and most retrotransposons, Odin retrotransposons encode a dual RNase H domain but no env gene. Phylogenetic analyses showed that Odin retrotransposons were sisters to lokiretroviruses, and lokiretroviruses and Odin retrotransposons were sisters to retroviruses, establishing an evolutionary framework to decipher the origin of retroviruses (canonical retroviruses and lokiretroviruses). Our results provided insights into the diversity and deep evolution of LTR retrotransposons closely related to retroviruses.

Little is known about the origin and long-term evolutionary mode of retroviruses. Retroviruses can integrate into their hosts' genomes, providing a molecular fossil record for studying their deep history. Here we report the discovery of an endogenous foamy virus-like element, which we designate 'coelacanth endogenous foamy-like virus' (CoeEFV), within the genome of the coelacanth (Latimeria chalumnae). Phylogenetic analyses place CoeEFV basal to all known foamy viruses, strongly suggesting an ancient ocean origin of this major retroviral lineage, which had previously been known to infect only land mammals. The discovery of CoeEFV reveals the presence of foamy-like viruses in species outside the Mammalia. We show that foamy-like viruses have likely codiverged with their vertebrate hosts for more than 407 million years and underwent an evolutionary transition from water to land with their vertebrate hosts. These findings suggest an ancient marine origin of retroviruses and have important implications in understanding foamy virus biology.

LTR retrotransposons comprise a major component of the genomes of eukaryotes. On occasion, retrotransposon genes can be recruited by their hosts for diverse functions, a process formally referred to as co-option. However, a comprehensive picture of LTR retrotransposon gag gene co-option in eukaryotes is still lacking, with several documented cases exclusively involving Ty3/Gypsy retrotransposons in animals. Here, we use a phylogenomic approach to systemically unearth co-option of retrotransposon gag genes above the family level of taxonomy in 2,011 eukaryotes, namely co-option occurring during the deep evolution of eukaryotes. We identify a total of 14 independent gag gene co-option events across more than 740 eukaryote families, eight of which have not been reported previously. Among these retrotransposon gag gene co-option events, nine, four, and one involve gag genes of Ty3/Gypsy, Ty1/Copia, and Bel-Pao retrotransposons, respectively. Seven, four, and three co-option events occurred in animals, plants, and fungi, respectively. Interestingly, two co-option events took place in the early evolution of angiosperms. Both selective pressure and gene expression analyses further support that these co-opted gag genes might perform diverse cellular functions in their hosts, and several co-opted gag genes might be subject to positive selection. Taken together, our results provide a comprehensive picture of LTR retrotransposon gag gene co-option events that occurred during the deep evolution of eukaryotes and suggest paucity of LTR retrotransposon gag gene co-option during the deep evolution of eukaryotes.

Discriminating self from nonself is fundamental to immunity. Yet, it remains largely elusive how the mechanisms of self and nonself discrimination originated. Sensing double-stranded RNA as nonself, the 2',5'-oligoadenylate synthetase (OAS)-ribonuclease L (RNase L) pathway represents a crucial component of innate immunity. Here, we combine phylogenomic and functional analyses to show that the functional OAS-RNase L pathway likely originated through tinkering with preexisting proteins before the rise of jawed vertebrates during or before the Silurian period (444 to 419 Mya). Multiple concerted losses of OAS and RNase L occurred during the evolution of jawed vertebrates, further supporting the ancient coupling between OAS and RNase L. Moreover, both OAS and RNase L genes evolved under episodic positive selection across jawed vertebrates, suggesting a long-running evolutionary arms race between the OAS-RNase L pathway and microbes. Our findings illuminate how an innate immune pathway originated via molecular tinkering.

Newcastle disease (ND), caused by ND virus (NDV), is one of the most serious illnesses of birds, particularly chickens, and has been one of the major causes of economic losses in the poultry industry. Live vaccines are widely used to prevent chicken from NDV all over the world. Given the implications that recombination has for RNA virus evolution, it is clearly important to determine the extent to which recombination plays a role in NDV evolution. In this study, we performed the phylogenetic and recombination analysis on complete NDV genomes. A natural multi-recombinant cockatoo/Indonesia/14698/90 (AY562985) was identified. Its two minor parental-like strains might be from the NDV vaccine lineage and anhinga/U.S.(Fl)/44083/93 lineage, respectively. Our study suggests that recombination plays a role in NDV evolution. Especially, the study also suggests that live vaccines have capacity to play roles in shaping NDV evolution by homologous recombination with circulating virus.

The ancestor of cetaceans underwent a macroevolutionary transition from land to water early in the Eocene Period >50 million years ago. However, little is known about how diverse retroviruses evolved during this shift from terrestrial to aquatic environments. Did retroviruses transition into water accompanying their hosts? Did retroviruses infect cetaceans through cross-species transmission after cetaceans invaded the aquatic environments? Endogenous retroviruses (ERVs) provide important molecular fossils for tracing the evolution of retroviruses during this macroevolutionary transition. Here, we use a phylogenomic approach to study the origin and evolution of ERVs in cetaceans. We identify a total of 8,724 ERVs within the genomes of 25 cetaceans, and phylogenetic analyses suggest these ERVs cluster into 315 independent lineages, each of which represents one or more independent endogenization events. We find that cetacean ERVs originated through two possible routes. 298 ERV lineages may derive from retrovirus endogenization that occurred before or during the transition from land to water of cetaceans, and most of these cetacean ERVs were reaching evolutionary dead-ends. 17 ERV lineages are likely to arise from independent retrovirus endogenization events that occurred after the split of mysticetes and odontocetes, indicating that diverse retroviruses infected cetaceans through cross-species transmission from non-cetacean mammals after the transition to aquatic life of cetaceans. Both integration time and synteny analyses support the recent or ongoing activity of multiple retroviral lineages in cetaceans, some of which proliferated into hundreds of copies within the host genomes. Although ERVs only recorded a proportion of past retroviral infections, our findings illuminate the complex evolution of retroviruses during one of the most marked macroevolutionary transitions in vertebrate history.

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has shown once again that coronavirus (CoV) in animals are potential sources for epidemics in humans. Porcine deltacoronavirus (PDCoV) is an emerging enteropathogen of swine with a worldwide distribution. Here, we implemented and described an approach to analyze the epidemiology of PDCoV following its emergence in the pig population. We performed an integrated analysis of full genome sequence data from 21 newly sequenced viruses, along with comprehensive epidemiological surveillance data collected globally over the last 15 years. We found four distinct phylogenetic lineages of PDCoV, which differ in their geographic circulation patterns. Interestingly, we identified more frequent intra- and interlineage recombination and higher virus genetic diversity in the Chinese lineages compared with the USA lineage where pigs are raised in different farming systems and ecological environments. Most recombination breakpoints are located in the ORF1ab gene rather than in genes encoding structural proteins. We also identified five amino acids under positive selection in the spike protein suggesting a role for adaptive evolution. According to structural mapping, three positively selected sites are located in the N-terminal domain of the S1 subunit, which is the most likely involved in binding to a carbohydrate receptor, whereas the other two are located in or near the fusion peptide of the S2 subunit and thus might affect membrane fusion. Finally, our phylogeographic investigations highlighted notable South-North transmission as well as frequent long-distance dispersal events in China that could implicate human-mediated transmission. Our findings provide new insights into the evolution and dispersal of PDCoV that contribute to our understanding of the critical factors involved in CoVs emergence.

Mitochondrial linear plasmids have been sporadically reported in fungi and plants. Yet, much remains obscure about the diversity, distribution, and evolution of mitochondrial linear plasmids. Here, through phylogenomic analyses across 7,163 cellular organisms (including 991 plants), we find that mitochondrial linear plasmids are widely present in land plants and fungi. Phylogenetic analyses indicate that plants are likely to have acquired mitochondrial linear plasmids horizontally from fungi before or during the conquest of terrestrial environments by plants. Gene content analyses show that mitochondrial linear plasmids harbor a highly dynamic and promiscuous repertoire of genes. Our study refines the understanding of the origin and evolution of mitochondrial linear plasmids.