Low frequency of paleoviral infiltration across the avian phylogeny

Characterization of Caulimovirid-like Sequences from Upland Cotton (Gossypium hirsutum L.) Exhibiting Terminal Abortion in Georgia, USA

In this study, we investigated the potential involvement of endogenous viral elements (EVEs) in the development of apical tissue necrosis, resulting in the terminal abortion of upland cotton (Gossypium hirsutum L.) in Georgia. The high-throughput sequence analysis of symptomatic and asymptomatic plant tissue samples revealed near-complete EVE-Georgia (EVE-GA) sequences closely related to caulimoviruses. The analysis of EVE-GA's putative open reading frames (ORFs) compared to cotton virus A and endogenous cotton pararetroviral elements (eCPRVE) revealed their similarity in putative ORFs 1-4. However, in the ORF 5 and ORF 6 encoding putative coat protein and reverse transcriptase, respectively, the sequences from EVE-GA have stop codons similar to eCPRVE sequences from Mississippi. In silico mining of the cotton genome database using EVE-GA as a query uncovered near-complete viral sequence insertions in the genomes of G. hirsutum species (~7 kb) but partial in G. tomentosum (~5.3 kb) and G. mustelinum (~5.1 kb) species. Furthermore, cotton EVEs' episomal forms and messenger RNA (mRNA) transcripts were detected in both symptomatic and asymptomatic plants collected from cotton fields. No significant yield difference was observed between symptomatic and asymptomatic plants of the two varieties evaluated in the experimental plot. Additionally, EVEs were also detected in cotton seeds and seedlings. This study emphasizes the need for future research on EVE sequences, their coding capacity, and any potential role in host immunity or pathogenicity.

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.

Adaptive expansion of ERVK solo-LTRs is associated with Passeriformes speciation events

Endogenous retroviruses (ERVs) are ancient retroviral remnants integrated in host genomes, and commonly deleted through unequal homologous recombination, leaving solitary long terminal repeats (solo-LTRs). This study, analysing the genomes of 362 bird species and their reptilian and mammalian outgroups, reveals an unusually higher level of solo-LTRs formation in birds, indicating evolutionary forces might have purged ERVs during evolution. Strikingly in the order Passeriformes, and especially the parvorder Passerida, endogenous retrovirus K (ERVK) solo-LTRs showed bursts of formation and recurrent accumulations coinciding with speciation events over past 22 million years. Moreover, our results indicate that the ongoing expansion of ERVK solo-LTRs in these bird species, marked by high transcriptional activity of ERVK retroviral genes in reproductive organs, caused variation of solo-LTRs between individual zebra finches. We experimentally demonstrated that cis-regulatory activity of recently evolved ERVK solo-LTRs may significantly increase the expression level of ITGA2 in the brain of zebra finches compared to chickens. These findings suggest that ERVK solo-LTRs expansion may introduce novel genomic sequences acting as cis-regulatory elements and contribute to adaptive evolution. Overall, our results underscore that the residual sequences of ancient retroviruses could influence the adaptive diversification of species by regulating host gene expression.

Different Host-Endogenous Retrovirus Relationships between Mammals and Birds Reflected in Genome-Wide Evolutionary Interaction Patterns

Mammals and birds differ largely in their average endogenous retrovirus loads, namely the proportion of endogenous retrovirus in the genome. The host-endogenous retrovirus relationships, including conflict and co-option, have been hypothesized among the causes of this difference. However, there has not been studies about the genomic evolutionary signal of constant host-endogenous retrovirus interactions in a long-term scale and how such interactions could lead to the endogenous retrovirus load difference. Through a phylogeny-controlled correlation analysis on ∼5,000 genes between the dN/dS ratio of each gene and the load of endogenous retrovirus in 12 mammals and 21 birds, separately, we detected genes that may have evolved in association with endogenous retrovirus loads. Birds have a higher proportion of genes with strong correlation between dN/dS and the endogenous retrovirus load than mammals. Strong evidence of association is found between the dN/dS of the coding gene for leucine-rich repeat-containing protein 23 and endogenous retrovirus load in birds. Gene set enrichment analysis shows that gene silencing rather than immunity and DNA recombination may have a larger contribution to the association between dN/dS and the endogenous retrovirus load for both mammals and birds. The above results together showing different evolutionary patterns between bird and mammal genes can partially explain the apparently lower endogenous retrovirus loads of birds, while gene silencing may be a universal mechanism that plays a remarkable role in the evolutionary interaction between the host and endogenous retrovirus. In summary, our study presents signals that the host genes might have driven or responded to endogenous retrovirus load changes in long-term evolution.

A Diverse Virome Is Identified in Parasitic Flatworms of Domestic Animals in Xinjiang, China

Parasitic flatworms infect diverse vertebrates and are major threats to animal and even human health; however, little is known about the virome of these lower life forms. Using viral metagenomic sequencing, we characterized the virome of the parasitic flatworms collected from major domestic animals, including Dicrocoelium lanceatum and Taenia hydatigena, Echinococcus granulosus sensu stricto and Echinococcus multilocularis. Seven and three different viruses were discovered from D. lanceatum and T. hydatigena, respectively, and no viral sequences were found in adult tapeworms and protoscoleces of E. granulosus sensu stricto and E. multilocularis. Two out of the five parasitic flatworm species carry viruses, showing a host specificity of these viruses. These viruses belong to the Parvoviridae, Circoviridae, unclassified circular, Rep-encoding single-stranded (CRESS) DNA virus, Rhabdoviridae, Endornaviridae, and unclassified RNA viruses. The presence of multiple highly divergent RNA viruses, especially those that cluster with viruses found in marine animals, implies a deep evolutionary history of parasite-associated viruses. In addition, we found viruses with high identity to common pathogens in dogs, including canine circovirus and canine parvovirus 2. The presence of these viruses in the parasites implies that they may infect parasitic flatworms but does not completely exclude the possibility of contamination from host intestinal contents. Furthermore, we demonstrated that certain viruses, such as CRESS DNA virus may integrate into the genome of their host. Our results expand the knowledge of viral diversity in parasites of important domestic animals, highlighting the need for further investigations of their prevalence among other parasites of key animals. IMPORTANCE Characterizing the virome of parasites is important for unveiling the viral diversity, evolution, and ecology and will help to understand the "Russian doll" pattern among viruses, parasites, and host animals. Our data indicate that diverse viruses are present in specific parasitic flatworms, including viruses that may have an ancient evolutionary history and viruses currently circulating in parasite-infected host animals. These data also raise the question of whether parasitic flatworms acquire and/or carry some viruses that may have transmission potential to animals. In addition, through the study of virus-parasite-host interactions, including the influence of viral infection on the life cycle of the parasite, as well as its fitness and pathogenicity to the host, we could find new strategies to prevent and control parasitic diseases.

Dinosaurs: Comparative Cytogenomics of Their Reptile Cousins and Avian Descendants

Reptiles known as dinosaurs pervade scientific and popular culture, while interest in their genomics has increased since the 1990s. Birds (part of the crown group Reptilia) are living theropod dinosaurs. Chromosome-level genome assemblies cannot be made from long-extinct biological material, but dinosaur genome organization can be inferred through comparative genomics of related extant species. Most reptiles apart from crocodilians have both macro- and microchromosomes; comparative genomics involving molecular cytogenetics and bioinformatics has established chromosomal relationships between many species. The capacity of dinosaurs to survive multiple extinction events is now well established, and birds now have more species in comparison with any other terrestrial vertebrate. This may be due, in part, to their karyotypic features, including a distinctive karyotype of around n = 40 (~10 macro and 30 microchromosomes). Similarity in genome organization in distantly related species suggests that the common avian ancestor had a similar karyotype to e.g., the chicken/emu/zebra finch. The close karyotypic similarity to the soft-shelled turtle (n = 33) suggests that this basic pattern was mostly established before the Testudine-Archosaur divergence, ~255 MYA. That is, dinosaurs most likely had similar karyotypes and their extensive phenotypic variation may have been mediated by increased random chromosome segregation and genetic recombination, which is inherently higher in karyotypes with more and smaller chromosomes.

Comparative analysis reveals the long-term coevolutionary history of parvoviruses and vertebrates

Parvoviruses (family Parvoviridae) are small DNA viruses that cause numerous diseases of medical, veterinary, and agricultural significance and have important applications in gene and anticancer therapy. DNA sequences derived from ancient parvoviruses are common in animal genomes and analysis of these endogenous parvoviral elements (EPVs) has demonstrated that the family, which includes twelve vertebrate-specific genera, arose in the distant evolutionary past. So far, however, such "paleovirological" analysis has only provided glimpses into the biology of ancient parvoviruses and their long-term evolutionary interactions with hosts. Here, we comprehensively map EPV diversity in 752 published vertebrate genomes, revealing defining aspects of ecology and evolution within individual parvovirus genera. We identify 364 distinct EPV sequences and show these represent approximately 200 unique germline incorporation events, involving at least five distinct parvovirus genera, which took place at points throughout the Cenozoic Era. We use the spatiotemporal and host range calibrations provided by these sequences to infer defining aspects of long-term evolution within individual parvovirus genera, including mammalian vicariance for genus Protoparvovirus, and interclass transmission for genus Dependoparvovirus. Moreover, our findings support a model of virus evolution in which the long-term cocirculation of multiple parvovirus genera in vertebrates reflects the adaptation of each viral genus to fill a distinct ecological niche. Our findings show that efforts to develop parvoviruses as therapeutic tools can be approached from a rational foundation based on comparative evolutionary analysis. To support this, we published our data in the form of an open, extensible, and cross-platform database designed to facilitate the wider utilisation of evolution-related domain knowledge in parvovirus research.

Ineffective silencing of transposable elements on an avian W Chromosome

One of the defining features of transposable elements (TEs) is their ability to move to new locations in the host genome. To minimise the potentially deleterious effects of de novo TE insertions, hosts have evolved several mechanisms to control TE activity, including recombination-mediated removal and epigenetic silencing; however, increasing evidence suggests that silencing of TEs is often incomplete. The crow family experienced a recent radiation of LTR retrotransposons (LTRs), offering an opportunity to gain insight into the regulatory control of young, potentially still active TEs. We quantified the abundance of TE-derived transcripts across several tissues in 15 Eurasian crows (Corvus (corone) spp.) raised under common garden conditions and find evidence for ineffective TE suppression on the female-specific W Chromosome. Using RNA-seq data, we show that ~ 9.5% of all transcribed TEs had considerably greater (average: 16-fold) transcript abundance in female crows, and that more than 85% of these female-biased TEs originated on the W Chromosome. After accounting for differences in TE density among chromosomal classes, W-linked TEs were significantly more highly expressed than TEs residing on other chromosomes, consistent with ineffective silencing on the former. Together, our results suggest that the crow W Chromosome acts as a source of transcriptionally active TEs, with possible negative fitness consequences for female birds analogous to Drosophila (an X/Y system), where overexpression of Y-linked TEs is associated with male-specific aging and fitness loss ('toxic Y').

The diversity of endogenous viral elements in insects

We provide an overview of the currently known diversity of viral sequences integrated into insect genomes. Such endogenous viral elements (EVE) have so far been annotated in at least eight insect orders and can be assigned to at least three families of large double-stranded (ds) DNA viruses, at least 22 families of RNA viruses, and three families of single-stranded DNA viruses. The study of these EVE has already produced important insights into insect-virus interactions, including the discovery of a new form of adaptive antiviral immunity. Insect EVE diversity will continue to increase as new insect genomes and exogenous viruses are sequenced, which will continue to make paleovirology a vibrant research field in this group of animals in the years to come.

Host contributions to the force of Borrelia burgdorferi and Babesia microti transmission differ at edges of and within a small habitat patch

In the northeastern United States, the emergence of Lyme disease has been associated, in part, with the increase of small forest patches. Such disturbed habitat is exploited by generalist species, such as white-footed mice, which are considered the host with the greatest reservoir capacity for the agents of Lyme disease (Borrelia burgdorferi sensu stricto) and human babesiosis (Babesia microti). Spatial risk analyses have identified edge habitat as particularly risky. Using a retrotransposon-based quantitative PCR assay for host bloodmeal remnant identification, we directly measured whether the hosts upon which vector ticks fed differed at the edge or within the contiguous small habitat patch. Questing nymphal deer ticks, Ixodes dammini, the northern clade of Ixodes scapularis, were collected from either the edge or within a thicket on Nantucket Island over 3 transmission seasons and tested for evidence of infection as well as bloodmeal hosts. Tick bloodmeal hosts significantly differed by site as well as by year. Mice and deer were identified most often (49.9%), but shrews, rabbits and birds were also common. Ticks from the edge fed on a greater diversity of hosts than those from the thicket. Surprisingly, mice were not strongly associated with either infection at either sampling site (OR<2 for all). Although shrews were not the most common host utilized by ticks, they were highly associated with both infections at both sites (OR= 4.5 and 7.9 B. burgdorferi and 7.9 and 19.0 B. microti, edge and thicket). We conclude that reservoir hosts may differ in their contributions to infecting ticks between edge and contiguous vegetated patches. Importance Habitat fragmentation is thought to be a main factor in the emergence of Lyme disease and other of the deer tick-transmitted infections. The patchwork of forest and edges promotes altered biodiversity, favoring the abundance of generalist rodents such as white footed mice, heretofore considered a key tick and reservoir host in the northeastern U.S. We used tick bloodmeal analyses to directly identify the hosts from which nymphal deer ticks became infected. We demonstrate that there is considerable microfocality in host contributions to the cohort of infected ticks and that shrews, although they fed fewer ticks than mice, disproportionately influenced the force of pathogen transmission in our site. The venue of transmission of certain deer tick-transmitted agents may comprise a habitat scale of 10 meters or fewer and depend on alternative small mammal hosts such as shrews.

Paleovirology of the DNA viruses of eukaryotes

Paleovirology is the study of ancient viruses and how they have coevolved with their hosts. An increasingly detailed understanding of the diversity, origins, and evolution of the DNA viruses of eukaryotes has been obtained through the lens of paleovirology in recent years. Members of multiple viral families have been found integrated in the genomes of eukaryotes, providing a rich fossil record to study. These elements have extended our knowledge of exogenous viral diversity, host ranges, and the timing of viral evolution, and are revealing the existence of entire new families of eukaryotic integrating dsDNA viruses and transposons. Future work in paleovirology will continue to provide insights into antiviral immunity, viral diversity, and potential applications, and reveal other secrets of the viral world.

Multiple Infiltration and Cross-Species Transmission of Foamy Viruses across the Paleozoic to the Cenozoic Era

Foamy viruses (FVs) are complex retroviruses that can infect humans and other animals. In this study, by integrating transcriptomic and genomic data, we discovered 412 FVs from 6 lineages in amphibians, which significantly increased the known set of FVs in amphibians. Among these lineages, salamander FVs maintained a coevolutionary pattern with their hosts that could be dated back to the Paleozoic era, while in contrast, frog FVs were much more likely acquired from cross-species (class-level) transmission in the Cenozoic era. In addition, we found that three distinct FV lineages had integrated into the genome of a salamander. Unexpectedly, we identified a lineage of endogenous FVs in caecilians that expressed all complete major genes, demonstrating the potential existence of an exogenous form of FV outside of mammals. Our discovery of rare phenomena in amphibian FVs has significantly increased our understanding of the macroevolution of the complex retrovirus. IMPORTANCE Foamy viruses (FVs) represent, more so than other viruses, the best model of coevolution between a virus and a host. This study represents the largest investigation so far of amphibian FVs and reveals 412 FVs of 6 distinct lineages from three major orders of amphibians. Besides a coevolutionary pattern, cross-species and repeated infections were also observed during the evolution of amphibian FVs. Remarkably, expressed FVs including a potential exogenous form were discovered, suggesting that active FVs might be underestimated in nature. These findings revealed that the multiple origins and complex evolution of amphibian FVs started from the Paleozoic era.

100-My history of bornavirus infections hidden in vertebrate genomes

Although viruses have threatened our ancestors for millions of years, prehistoric epidemics of viruses are largely unknown. Endogenous bornavirus-like elements (EBLs) are ancient bornavirus sequences derived from the viral messenger RNAs that were reverse transcribed and inserted into animal genomes, most likely by retrotransposons. These elements can be used as molecular fossil records to trace past bornaviral infections. In this study, we systematically identified EBLs in vertebrate genomes and revealed the history of bornavirus infections over nearly 100 My. We confirmed that ancient bornaviral infections have occurred in diverse vertebrate lineages, especially in primate ancestors. Phylogenetic analyses indicated that primate ancestors were infected with various bornaviral lineages during evolution. EBLs in primate genomes formed clades according to their integration ages, suggesting that bornavirus lineages infected with primate ancestors had changed chronologically. However, some bornaviral lineages may have coexisted with primate ancestors and underwent repeated endogenizations for tens of millions of years. Moreover, a bornaviral lineage that coexisted with primate ancestors also endogenized in the genomes of some ancestral bats. The habitats of these bat ancestors have been reported to overlap with the migration route of primate ancestors. These results suggest that long-term virus-host coexistence expanded the geographic distributions of the bornaviral lineage along with primate migration and may have spread their infections to these bat ancestors. Our findings provide insight into the history of bornavirus infections over geological timescales that cannot be deduced from research using extant viruses alone, thus broadening our perspective on virus-host coevolution.

Ancient evolution of hepadnaviral paleoviruses and their impact on host genomes

Hepadnaviruses (family Hepadnaviviridae) are reverse-transcribing animal viruses that infect vertebrates. DNA sequences derived from ancient hepadnaviruses have been identified in the germline genome of numerous vertebrate species, and these ‘endogenous hepatitis B viruses’ (eHBVs) reveal aspects of the long-term coevolutionary relationship between hepadnaviruses and their vertebrate hosts. Here, we use a novel, data-oriented approach to recover and analyse the complete repertoire of eHBV elements in published animal genomes. We show that germline incorporation of hepadnaviruses is exclusive to a single vertebrate group (Sauria) and that the eHBVs contained in saurian genomes represent a far greater diversity of hepadnaviruses than previously recognized. Through in-depth characterization of eHBV elements, we establish the existence of four distinct subgroups within the genus Avihepadnavirus and trace their evolution through the Cenozoic Era. Furthermore, we provide a completely new perspective on hepadnavirus evolution by showing that the metahepadnaviruses (genus Metahepadnavirus) originated >300 million years ago in the Paleozoic Era and have historically infected a broad range of vertebrates. We also show that eHBVs have been intra-genomically amplified in some saurian lineages, and that eHBVs located at approximately equivalent genomic loci have been acquired in entirely distinct germline integration events. These findings indicate that selective forces have favoured the accumulation of hepadnaviral sequences at specific loci in the saurian germline. Our investigation provides a range of new insights into the long-term evolutionary history of reverse-transcribing DNA viruses and shows that germline incorporation of hepadnaviruses has played a role in shaping the evolution of saurian genomes.

Fibrillarin evolution through the Tree of Life: Comparative genomics and microsynteny network analyses provide new insights into the evolutionary history of Fibrillarin

Fibrillarin (FIB), a methyltransferase essential for life in the vast majority of eukaryotes, is involved in methylation of rRNA required for proper ribosome assembly, as well as methylation of histone H2A of promoter regions of rRNA genes. RNA viral progression that affects both plants and animals requires FIB proteins. Despite the importance and high conservation of fibrillarins, there little is known about the evolutionary dynamics of this small gene family. We applied a phylogenomic microsynteny-network approach to elucidate the evolutionary history of FIB proteins across the Tree of Life. We identified 1063 non-redundant FIB sequences across 1049 completely sequenced genomes from Viruses, Bacteria, Archaea, and Eukarya. FIB is a highly conserved single-copy gene through Archaea and Eukarya lineages, except for plants, which have a gene family expansion due to paleopolyploidy and tandem duplications. We found a high conservation of the FIB genomic context during plant evolution. Surprisingly, FIB in mammals duplicated after the Eutheria split (e.g., ruminants, felines, primates) from therian mammals (e.g., marsupials) to form two main groups of sequences, the FIB and FIB-like groups. The FIB-like group transposed to another genomic context and remained syntenic in all the eutherian mammals. This transposition correlates with differences in the expression patterns of FIB-like proteins and with elevated Ks values potentially due to reduced evolutionary constraints of the duplicated copy. Our results point to a unique evolutionary event in mammals, between FIB and FIB-like genes, that led to non-redundant roles of the vital processes in which this protein is involved.

Update on immunopathology of bornavirus infections in humans and animals

Knowledge on bornaviruses has expanded tremendously during the last decade through detection of novel bornaviruses and endogenous bornavirus-like elements in many eukaryote genomes, as well as by confirmation of insectivores as reservoir species for classical Borna disease virus 1 (BoDV-1). The most intriguing finding was the demonstration of the zoonotic potential of lethal human bornavirus infections caused by a novel bornavirus of different squirrel species (variegated squirrel 1 bornavirus, VSBV-1) and by BoDV-1 known as the causative agent for the classical Borna disease in horses and sheep. Whereas a T cell-mediated immunopathology has already been confirmed as key disease mechanism for infection with BoDV-1 by experimental studies in rodents, the underlying pathomechanisms remain less clear for human bornavirus infections, infection with other bornaviruses or infection of reservoir species. Thus, an overview of current knowledge on the pathogenesis of bornavirus infections focusing on immunopathology is given.

Identification and evolution of avian endogenous foamy viruses

A history of long-term co-divergence means that foamy viruses (family Retroviridae) provide an ideal framework to understanding virus-host evolution over extended time periods. Endogenous foamy viruses (EndFVs) are rare, and to date have only been described in a limited number of mammals, amphibians, reptiles and fish genomes. By screening 414 avian genomes we identified EndFVs in two bird species: the Maguari Stork (Ciconia maguari) and the Oriental Stork (Ciconia boyciana). Analyses of phylogenetic relationships, genome structures and flanking sequences revealed a single origin of EndFVs in Ciconia species. In addition, the marked incongruence between the virus and host phylogenies suggested that this integration event occurred independently in birds. In sum, by providing evidence that birds can be infected with foamy viruses, we fill the last major gap in the taxonomic distribution of foamy viruses and their animal hosts.

Genomic and transcriptional analyses of novel parvoviruses identified from dead peafowl

To identify potential pathogens responsible for a disease outbreak of cultured peafowls in China in 2013, metagenomic sequencing was conducted. The genomes of two closely related parvoviruses, namely peafowl parvovirus 1 (PePV1) and PePV2, were identified with size of 4428 bp and 4348 bp, respectively. Phylogenetic analysis revealed that both viruses are novel parvoviruses, belonging to the proposed genus Chapparvovirus of Parvoviridae. The transcriptional profile of PePV1 was analyzed by transfecting a nearly complete PePV1 genome into HEK-293T cells. Results revealed that PePV1 employs one promoter and two polyadenylation sites to start and terminate its transcriptions, with one donor site and two acceptor sites for pre-mRNA splicing. PePV1 DNA and structural protein were detected in several tissues of a dead peafowl, which appeared to have suffered enteritis, pneumonia and viremia. These results provide novel information of chapparvoviruses, and call for attention to the potential pathogens.

Time lapse: A glimpse into prehistoric genomics

For the purpose of this review, ‘time-lapse’ refers to the reconstruction of ancestral (in this case dinosaur) karyotypes using genome assemblies of extant species. Such reconstructions are only usually possible when genomes are assembled to ‘chromosome level’ i.e. a complete representation of all the sequences, correctly ordered contiguously on each of the chromosomes. Recent paleontological evidence is very clear that birds are living dinosaurs, the latest example of dinosaurs emerging from a catastrophic extinction event. Non-avian dinosaurs (ever present in the public imagination through art, and broadcast media) emerged some 240 million years ago and have displayed incredible phenotypic diversity. Here we report on our recent studies to infer the overall karyotype of the Theropod dinosaur lineage from extant avian chromosome level genome assemblies. Our work first focused on determining the likely karyotype of the avian ancestor (most likely a chicken-sized, two-legged, feathered, land dinosaur from the Jurassic period) finding karyotypic similarity to the chicken. We then took the work further to determine the likely karyotype of the bird-lizard ancestor and the chromosomal changes (chiefly translocations and inversions) that occurred between then and modern birds. A combination of bioinformatics and cross-species fluorescence in situ hybridization (zoo-FISH) uncovered a considerable number of translocations and fissions from a ‘lizard-like’ genome structure of 2n = 36–46 to one similar to that of soft-shelled turtles (2n = 66) from 275 to 255 million years ago (mya). Remarkable karyotypic similarities between some soft-shelled turtles and chicken suggests that there were few translocations from the bird-turtle ancestor (plus ∼7 fissions) through the dawn of the dinosaurs and pterosaurs, through the theropod linage and on to most to modern birds. In other words, an avian-like karyotype was in place about 240mya when the dinosaurs and pterosaurs first emerged. We mapped 49 chromosome inversions from then to the present day, uncovering some gene ontology enrichment in evolutionary breakpoint regions. This avian-like karyotype with its many (micro)chromosomes provides the basis for variation (the driver of natural selection) through increased random segregation and recombination. It may therefore contribute to the ability of dinosaurs to survive multiple extinction events, emerging each time as speciose and diverse.

A reptilian endogenous foamy virus sheds light on the early evolution of retroviruses

Endogenous retroviruses (ERVs) represent host genomic 'fossils' of ancient viruses. Foamy viruses, including those that form endogenous copies, provide strong evidence for virus-host co-divergence across the vertebrate phylogeny. Endogenous foamy viruses (EFVs) have previously been discovered in mammals, amphibians, and fish. Here we report a novel endogenous foamy virus, termed ERV-Spuma-Spu, in genome of the tuatara (Sphenodon punctatus), an endangered reptile species endemic to New Zealand. Phylogenetic analyses revealed that foamy viruses have likely co-diverged with their hosts over many millions of years. The discovery of ERV-Spuma-Spu fills a major gap in the fossil record of foamy viruses and provides important insights into the early evolution of retroviruses.

Discovery of endogenous retroviruses with mammalian envelopes in avian genomes uncovers long-term bird-mammal interaction

Endogenous retroviruses (ERVs) arise from the infection and integration of past retroviruses into animal hosts. We performed large-scale genomic mining of 101 avian genomes for discovery of ERVs having none-avian origin and investigated the cross-species transmission events. Phylogenetic analysis of the reverse transcriptase (RT) of polymerase gene (pol) and the transmembrane subunit (TM) of the envelope gene (env) supported that avian ERVs with a mammalian env gene existed in at least 15 avian species and can be divided into two major groups: Group-1 were of recombinant ERVs with an alpha-like pol gene and a gamma-like env gene, and Group-2 included ERVs with both gamma-like pol and env genes. Group-1 represented the avian alpharetroviral/mammalian gammaretroviral recombinant while Group-2 documented viral jump from mammals to birds. Molecular dating analysis suggested that Group-1 ERVs had integrated into avian genomes continuously, until recent past. We have expanded the knowledge of ERVs with cross-order transmission.

Immune genes are hotspots of shared positive selection across birds and mammals

Consistent patterns of positive selection in functionally similar genes can suggest a common selective pressure across a group of species. We use alignments of orthologous protein-coding genes from 39 species of birds to estimate parameters related to positive selection for 11,000 genes conserved across birds. We show that functional pathways related to the immune system, recombination, lipid metabolism, and phototransduction are enriched for positively selected genes. By comparing our results with mammalian data, we find a significant enrichment for positively selected genes shared between taxa, and that these shared selected genes are enriched for viral immune pathways. Using pathogen-challenge transcriptome data, we show that genes up-regulated in response to pathogens are also enriched for positively selected genes. Together, our results suggest that pathogens, particularly viruses, consistently target the same genes across divergent clades, and that these genes are hotspots of host-pathogen conflict over deep evolutionary time.

Functional Evolution of Avian RIG-I-Like Receptors

RIG-I-like receptors (retinoic acid-inducible gene-I-like receptors, or RLRs) are family of pattern-recognition receptors for RNA viruses, consisting of three members: retinoic acid-inducible gene I (RIG-I), melanoma differentiation-associated gene 5 (MDA5) and laboratory of genetics and physiology 2 (LGP2). To understand the role of RLRs in bird evolution, we performed molecular evolutionary analyses on the coding genes of avian RLRs using filtered predicted coding sequences from 62 bird species. Among the three RLRs, conservation score and dN/dS (ratio of nonsynonymous substitution rate over synonymous substitution rate) analyses indicate that avian MDA5 has the highest conservation level in the helicase domain but a lower level in the caspase recruitment domains (CARDs) region, which differs from mammals; LGP2, as a whole gene, has a lower conservation level than RIG-I or MDA5. We found evidence of positive selection across all bird lineages in RIG-I and MDA5 but only on the stem lineage of Galliformes in LGP2, which could be related to the loss of RIG-I in Galliformes. Analyses also suggest that selection relaxation may have occurred in LGP2 during the middle of bird evolution and the CARDs region of MDA5 contains many positively selected sites, which might explain its conservation level. Spearman’s correlation test indicates that species-to-ancestor dN/dS of RIG-I shows a negative correlation with endogenous retroviral abundance in bird genomes, suggesting the possibility of interaction between immunity and endogenous retroviruses during bird evolution.

Endogenous amdoparvovirus-related elements reveal insights into the biology and evolution of vertebrate parvoviruses

Amdoparvoviruses (family Parvoviridae: genus Amdoparvovirus) infect carnivores, and are a major cause of morbidity and mortality in farmed animals. In this study, we systematically screened animal genomes to identify endogenous parvoviral elements (EPVs) disclosing a high degree of similarity to amdoparvoviruses, and investigated their genomic, phylogenetic and protein structural features. We report the first examples of full-length, amdoparvovirus-derived EPVs in the genome of the Transcaucasian mole vole (Ellobius lutescens). We also identify four EPVs in mammal and reptile genomes that are intermediate between amdoparvoviruses and their sister genus (Protoparvovirus) in terms of their phylogenetic placement and genomic features. In particular, we identify a genome-length EPV in the genome of a pit viper (Protobothrops mucrosquamatus) that is more similar to a protoparvovirus than an amdoparvovirus in terms of its phylogenetic placement and the structural features of its capsid protein (as revealed by homology modeling), yet exhibits characteristically amdoparvovirus-like genome features including: (1) a putative middle ORF gene; (2) a capsid gene that lacks a phospholipase A2 domain; (3) a genome structure consistent with an amdoparvovirus-like mechanism of capsid gene expression. Our findings indicate that amdoparvovirus host range extends to rodents, and that parvovirus lineages possessing a mixture of proto- and amdoparvovirus-like characteristics have circulated in the past. In addition, we show that EPV sequences in the mole vole and pit viper encode intact, expressible replicase genes that have potentially been co-opted or exapted in these host species.

Reconstruction of the diapsid ancestral genome permits chromosome evolution tracing in avian and non-avian dinosaurs

Genomic organisation of extinct lineages can be inferred from extant chromosome-level genome assemblies. Here, we apply bioinformatic and molecular cytogenetic approaches to determine the genomic structure of the diapsid common ancestor. We then infer the events that likely occurred along this lineage from theropod dinosaurs through to modern birds. Our results suggest that most elements of a typical ‘avian-like’ karyotype (40 chromosome pairs, including 30 microchromosomes) were in place before the divergence of turtles from birds ~255 mya. This genome organisation therefore predates the emergence of early dinosaurs and pterosaurs and the evolution of flight. Remaining largely unchanged interchromosomally through the dinosaur–theropod route that led to modern birds, intrachromosomal changes nonetheless reveal evolutionary breakpoint regions enriched for genes with ontology terms related to chromatin organisation and transcription. This genomic structure therefore appears highly stable yet contributes to a large degree of phenotypic diversity, as well as underpinning adaptive responses to major environmental disruptions via intrachromosomal repatterning.

Ancient diapsids diverged into the lineages leading to turtles and birds over 250 million years ago. Here, the authors use genomic and molecular cytogenetic analyses of modern species to infer the genome structure of the diapsid common ancestor (DCA) and the changes occurring along the lineage to birds through theropod dinosaurs.

Paleovirology of bornaviruses: What can be learned from molecular fossils of bornaviruses

Endogenous viral elements (EVEs) are virus-derived sequences embedded in eukaryotic genomes formed by germline integration of viral sequences. As many EVEs were integrated into eukaryotic genomes millions of years ago, EVEs are considered molecular fossils of viruses. EVEs can be valuable informational sources about ancient viruses, including their time scale, geographical distribution, genetic information, and hosts. Although integration of viral sequences is not required for replications of viruses other than retroviruses, many non-retroviral EVEs have been reported to exist in eukaryotes. Investigation of these EVEs has expanded our knowledge regarding virus-host interactions, as well as provided information on ancient viruses. Among them, EVEs derived from bornaviruses, non-retroviral RNA viruses, have been relatively well studied. Bornavirus-derived EVEs are widely distributed in animal genomes, including the human genome, and the history of bornaviruses can be dated back to more than 65 million years. Although there are several reports focusing on the biological significance of bornavirus-derived sequences in mammals, paleovirology of bornaviruses has not yet been well described and summarized. In this paper, we describe what can be learned about bornaviruses from endogenous bornavirus-like elements from the view of paleovirology using published results and our novel data.

The evolution, distribution and diversity of endogenous circoviral elements in vertebrate genomes

•

Sequences derived from circoviruses occur sporadically in vertebrate genomes.

•

The vast majority of fixed circovirus sequences in vertebrate genomes are derived from rep genes.

•

Some ancient circoviruses were closely related to modern circoviruses.

•

Integrated circovirus sequences are highly duplicated in some carnivore lineages.

•

Sequences in the Mus spretus genome are closely related to modern canine circoviruses.

Circoviruses (family Circoviridae) are small, non-enveloped viruses that have short, single-stranded DNA genomes. Circovirus sequences are frequently recovered in metagenomic investigations, indicating that these viruses are widespread, yet they remain relatively poorly understood. Endogenous circoviral elements (CVe) are DNA sequences derived from circoviruses that occur in vertebrate genomes. CVe are a useful source of information about the biology and evolution of circoviruses. In this study, we screened 362 vertebrate genome assemblies in silico to generate a catalog of CVe loci. We identified a total of 179 CVe sequences, most of which have not been reported previously. We show that these CVe loci reflect at least 19 distinct germline integration events. We determine the structure of CVe loci, identifying some that show evidence of potential functionalization. We also identify orthologous copies of CVe in snakes, fish, birds, and mammals, allowing us to add new calibrations to the timeline of circovirus evolution. Finally, we observed that some ancient CVe group robustly with contemporary circoviruses in phylogenies, with all sequences within these groups being derived from the same host class or order, implying a hitherto underappreciated stability in circovirus-host relationships. The openly available dataset constructed in this investigation provides new insights into circovirus evolution, and can be used to facilitate further studies of circoviruses and CVe.

Divergent bornaviruses from Australian carpet pythons with neurological disease date the origin of extant Bornaviridae prior to the end-Cretaceous extinction

Tissue samples from Australian carpet pythons (Morelia spilota) with neurological disease were screened for viruses using next-generation sequencing. Coding complete genomes of two bornaviruses were identified with the gene order 3'-N-X-P-G-M-L, representing a transposition of the G and M genes compared to other bornaviruses and most mononegaviruses. Use of these viruses to search available vertebrate genomes enabled recognition of further endogenous bornavirus-like elements (EBLs) in diverse placental mammals, including humans. Codivergence patterns and shared integration sites revealed an ancestral laurasiatherian EBLG integration (77 million years ago [MYA]) and a previously identified afrotherian EBLG integration (83 MYA). The novel python bornaviruses clustered more closely with these EBLs than with other exogenous bornaviruses, suggesting that these viruses diverged from previously known bornaviruses prior to the end-Cretaceous (K-Pg) extinction, 66 MYA. It is possible that EBLs protected mammals from ancient bornaviral disease, providing a selective advantage in the recovery from the K-Pg extinction. A degenerate PCR primer set was developed to detect a highly conserved region of the bornaviral polymerase gene. It was used to detect 15 more genetically distinct bornaviruses from Australian pythons that represent a group that is likely to contain a number of novel species.

Extensive diversity and evolution of hepadnaviruses in bats in China

To better understand the evolution of hepadnaviruses, we sampled bats from Guizhou, Henan and Zhejiang provinces, China, and rodents from Zhejiang province. Genetically diverse hepadnaviruses were identified in a broad range of bat species, with an overall prevalence of 13.3%. In contrast, no rodent hepadnaviruses were identified. The newly discovered bat hepadnaviruses fell into two distinct phylogenetic groups. The viruses within the first group exhibited high diversity, with some closely related to viruses previously identified in Yunnan province. Strikingly, the newly discovered viruses sampled from Jiyuan city in the second phylogenetic group were most closely related to those found in bats from West Africa, suggestive of a long-term association between bats and hepadnaviruses. A co-phylogenetic analysis revealed frequent cross-species transmission among bats from different species, genera, and families. Overall, these data suggest that there are likely few barriers to the cross-species transmission of bat hepadnaviruses.

•

Diverse hepadnaviruses are identified in a broad range of bat species in China.

•

Some of them were closely related to those previously identified in China.

•

The viruses from Jiyuan were most closely related to Gabon bat hepadnaviruses.

•

Newly discovered viruses did not clustered by bat species or geographic location.

•

Frequent cross-species transmission among different bat species was observed.

New high copy tandem repeat in the content of the chicken W chromosome

The content of repetitive DNA in avian genomes is considerably less than in other investigated vertebrates. The first descriptions of tandem repeats were based on the results of routine biochemical and molecular biological experiments. Both satellite DNA and interspersed repetitive elements were annotated using library-based approach and de novo repeat identification in assembled genome. The development of deep-sequencing methods provides datasets of high quality without preassembly allowing one to annotate repetitive elements from unassembled part of genomes. In this work, we search the chicken assembly and annotate high copy number tandem repeats from unassembled short raw reads. Tandem repeat (GGAAA)_n has been identified and found to be the second after telomeric repeat (TTAGGG)n most abundant in the chicken genome. Furthermore, (GGAAA)_n repeat forms expanded arrays on the both arms of the chicken W chromosome. Our results highlight the complexity of repetitive sequences and update data about organization of sex W chromosome in chicken.

Low diversity, activity, and density of transposable elements in five avian genomes

In this study, we conducted the activity, diversity, and density analysis of transposable elements (TEs) across five avian genomes (budgerigar, chicken, turkey, medium ground finch, and zebra finch) to explore the potential reason of small genome sizes of birds. We found that these avian genomes exhibited low density of TEs by about 10% of genome coverages and low diversity of TEs with the TE landscapes dominated by CR1 and ERV elements, and contrasting proliferation dynamics both between TE types and between species were observed across the five avian genomes. Phylogenetic analysis revealed that CR1 clade was more diverse in the family structure compared with R2 clade in birds; avian ERVs were classified into four clades (alpha, beta, gamma, and ERV-L) and belonged to three classes of ERV with an uneven distributed in these lineages. The activities of DNA and SINE TEs were very low in the evolution history of avian genomes; most LINEs and LTRs were ancient copies with a substantial decrease of activity in recent, with only LTRs and LINEs in chicken and zebra finch exhibiting weak activity in very recent, and very few TEs were intact; however, the recent activity may be underestimated due to the sequencing/assembly technologies in some species. Overall, this study demonstrates low diversity, activity, and density of TEs in the five avian species; highlights the differences of TEs in these lineages; and suggests that the current and recent activity of TEs in avian genomes is very limited, which may be one of the reasons of small genome sizes in birds.

The biological significance of bornavirus-derived genes in mammals

The discoveries of sequences derived from non-retroviral viruses in eukaryotic genomes have significantly expanded our knowledge about virus evolution as well as the co-evolution between viruses and eukaryotes. However, the biological functions of such sequences in the host are largely unknown. Endogenous bornavirus-like elements (EBLs) have been relatively well studied by molecular biological methods, which have provided evidence that some EBLs have been co-opted by their hosts. This review highlights the current knowledge on the biological significance of EBLs, and discusses possible functions of EBLs. Further, we highlight the importance of extensive surveillance of exogenous viruses for a better understanding of endogenous viral sequences as well as the co-evolution of viruses and eukaryotes.

Evolution and Diversity of Transposable Elements in Vertebrate Genomes

Transposable elements (TEs) are selfish genetic elements that mobilize in genomes via transposition or retrotransposition and often make up large fractions of vertebrate genomes. Here, we review the current understanding of vertebrate TE diversity and evolution in the context of recent advances in genome sequencing and assembly techniques. TEs make up 4-60% of assembled vertebrate genomes, and deeply branching lineages such as ray-finned fishes and amphibians generally exhibit a higher TE diversity than the more recent radiations of birds and mammals. Furthermore, the list of taxa with exceptional TE landscapes is growing. We emphasize that the current bottleneck in genome analyses lies in the proper annotation of TEs and provide examples where superficial analyses led to misleading conclusions about genome evolution. Finally, recent advances in long-read sequencing will soon permit access to TE-rich genomic regions that previously resisted assembly including the gigantic, TE-rich genomes of salamanders and lungfishes.

Perspectives from the Avian Phylogenomics Project: Questions that Can Be Answered with Sequencing All Genomes of a Vertebrate Class

The rapid pace of advances in genome technology, with concomitant reductions in cost, makes it feasible that one day in our lifetime we will have available extant genomes of entire classes of species, including vertebrates. I recently helped cocoordinate the large-scale Avian Phylogenomics Project, which collected and sequenced genomes of 48 bird species representing most currently classified orders to address a range of questions in phylogenomics and comparative genomics. The consortium was able to answer questions not previously possible with just a few genomes. This success spurred on the creation of a project to sequence the genomes of at least one individual of all extant ∼10,500 bird species. The initiation of this project has led us to consider what questions now impossible to answer could be answered with all genomes, and could drive new questions now unimaginable. These include the generation of a highly resolved family tree of extant species, genome-wide association studies across species to identify genetic substrates of many complex traits, redefinition of species and the species concept, reconstruction of the genomes of common ancestors, and generation of new computational tools to address these questions. Here I present visions for the future by posing and answering questions regarding what scientists could potentially do with available genomes of an entire vertebrate class.

Complete mitochondrial genomes of living and extinct pigeons revise the timing of the columbiform radiation

BACKGROUND

Pigeons and doves (Columbiformes) are one of the oldest and most diverse extant lineages of birds. However, the nature and timing of the group's evolutionary radiation remains poorly resolved, despite recent advances in DNA sequencing and assembly and the growing database of pigeon mitochondrial genomes. One challenge has been to generate comparative data from the large number of extinct pigeon lineages, some of which are morphologically unique and therefore difficult to place in a phylogenetic context.

RESULTS

We used ancient DNA and next generation sequencing approaches to assemble complete mitochondrial genomes for eleven pigeons, including the extinct Ryukyu wood pigeon (Columba jouyi), the thick-billed ground dove (Alopecoenas salamonis), the spotted green pigeon (Caloenas maculata), the Rodrigues solitaire (Pezophaps solitaria), and the dodo (Raphus cucullatus). We used a Bayesian approach to infer the evolutionary relationships among 24 species of living and extinct pigeons and doves.

CONCLUSIONS

Our analyses indicate that the earliest radiation of the Columbidae crown group most likely occurred during the Oligocene, with continued divergence of major clades into the Miocene, suggesting that diversification within the Columbidae occurred more recently than has been reported previously.

Distinct Viral Lineages from Fish and Amphibians Reveal the Complex Evolutionary History of Hepadnaviruses

Hepadnaviruses (hepatitis B viruses [HBVs]) are the only animal viruses that replicate their DNA by reverse transcription of an RNA intermediate. Until recently, the known host range of hepadnaviruses was limited to mammals and birds. We obtained and analyzed the first amphibian HBV genome, as well as several prototype fish HBVs, which allow the first comprehensive comparative genomic analysis of hepadnaviruses from four classes of vertebrates. Bluegill hepadnavirus (BGHBV) was characterized from in-house viral metagenomic sequencing. The African cichlid hepadnavirus (ACHBV) and the Tibetan frog hepadnavirus (TFHBV) were discovered using in silico analyses of the whole-genome shotgun and transcriptome shotgun assembly databases. Residues in the hydrophobic base of the capsid (core) proteins, designated motifs I, II, and III, are highly conserved, suggesting that structural constraints for proper capsid folding are key to capsid protein evolution. Surface proteins in all vertebrate HBVs contain similar predicted membrane topologies, characterized by three transmembrane domains. Most striking was the fact that BGHBV, ACHBV, and the previously described white sucker hepadnavirus did not form a fish-specific monophyletic group in the phylogenetic analysis of all three hepadnaviral genes. Notably, BGHBV was more closely related to the mammalian hepadnaviruses, indicating that cross-species transmission events have played a major role in viral evolution. Evidence of cross-species transmission was also observed with TFHBV. Hence, these data indicate that the evolutionary history of the hepadnaviruses is more complex than previously realized and combines both virus-host codivergence over millions of years and host species jumping.

IMPORTANCE

Hepadnaviruses are responsible for significant disease in humans (hepatitis B virus) and have been reported from a diverse range of vertebrates as both exogenous and endogenous viruses. We report the full-length genome of a novel hepadnavirus from a fish and the first hepadnavirus genome from an amphibian. The novel fish hepadnavirus, sampled from bluegills, was more closely related to mammalian hepadnaviruses than to other fish viruses. This phylogenetic pattern reveals that, although hepadnaviruses have likely been associated with vertebrates for hundreds of millions of years, they have also been characterized by species jumping across wide phylogenetic distances.

Genomic Landscape of Long Terminal Repeat Retrotransposons (LTR-RTs) and Solo LTRs as Shaped by Ectopic Recombination in Chicken and Zebra Finch

Transposable elements (TEs) are nearly ubiquitous among eukaryotic genomes, but TE contents vary dramatically among phylogenetic lineages. Several mechanisms have been proposed as drivers of TE dynamics in genomes, including the fixation/loss of a particular TE insertion by selection or drift as well as structural changes in the genome due to mutation (e.g., recombination). In particular, recombination can have a significant and directional effect on the genomic TE landscape. For example, ectopic recombination removes internal regions of long terminal repeat retrotransposons (LTR-RTs) as well as one long terminal repeat (LTR), resulting in a solo LTR. In this study, we focus on the intra-species dynamics of LTR-RTs and solo LTRs in bird genomes. The distribution of LTR-RTs and solo LTRs in birds is intriguing because avian recombination rates vary widely within a given genome. We used published linkage maps and whole genome assemblies to study the relationship between recombination rates and LTR-removal events in the chicken and zebra finch. We hypothesized that regions with low recombination rates would harbor more full-length LTR-RTs (and fewer solo LTRs) than regions with high recombination rates. We tested this hypothesis by comparing the ratio of full-length LTR-RTs and solo LTRs across chromosomes, across non-overlapping megabase windows, and across physical features (i.e., centromeres and telomeres). The chicken data statistically supported the hypothesis that recombination rates are inversely correlated with the ratio of full-length to solo LTRs at both the chromosome level and in 1-Mb non-overlapping windows. We also found that the ratio of full-length to solo LTRs near chicken telomeres was significantly lower than those ratios near centromeres. Our results suggest a potential role of ectopic recombination in shaping the chicken LTR-RT genomic landscape.

Genomic Mining Reveals Deep Evolutionary Relationships between Bornaviruses and Bats

Bats globally harbor viruses in order Mononegavirales, such as lyssaviruses and henipaviruses; however, little is known about their relationships with bornaviruses. Previous studies showed that viral fossils of bornaviral origin are embedded in the genomes of several mammalian species such as primates, indicative of an ancient origin of exogenous bornaviruses. In this study, we mined the available 10 bat genomes and recreated a clear evolutionary relationship of endogenous bornaviral elements and bats. Comparative genomics showed that endogenization of bornaviral elements frequently occurred in vesper bats, harboring EBLLs (endogenous bornavirus-like L elements) in their genomes. Molecular dating uncovered a continuous bornavirus-bat interaction spanning 70 million years. We conclude that better understanding of modern exogenous bornaviral circulation in bat populations is warranted.

Characterization of a Novel Hepadnavirus in the White Sucker (Catostomus commersonii) from the Great Lakes Region of the United States

The white sucker Catostomus commersonii is a freshwater teleost often utilized as a resident sentinel. Here, we sequenced the full genome of a hepatitis B-like virus that infects white suckers from the Great Lakes Region of the United States. Dideoxy sequencing confirmed that the white sucker hepatitis B virus (WSHBV) has a circular genome (3,542 bp) with the prototypical codon organization of hepadnaviruses. Electron microscopy demonstrated that complete virions of approximately 40 nm were present in the plasma of infected fish. Compared to avi- and orthohepadnaviruses, sequence conservation of the core, polymerase, and surface proteins was low and ranged from 16 to 27% at the amino acid level. An X protein homologue common to the orthohepadnaviruses was not present. The WSHBV genome included an atypical, presumptively noncoding region absent in previously described hepadnaviruses. Phylogenetic analyses confirmed WSHBV as distinct from previously documented hepadnaviruses. The level of divergence in protein sequences between WSHBV and other hepadnaviruses and the identification of an HBV-like sequence in an African cichlid provide evidence that a novel genus of the family Hepadnaviridae may need to be established that includes these hepatitis B-like viruses in fishes. Viral transcription was observed in 9.5% (16 of 169) of white suckers evaluated. The prevalence of hepatic tumors in these fish was 4.9%, and only 2.4% of fish were positive for both virus and hepatic tumors. These results are not sufficient to draw inferences regarding the association of WSHBV and carcinogenesis in white sucker.

IMPORTANCE We report the first full-length genome of a hepadnavirus from fishes. Phylogenetic analysis of this genome indicates divergence from genomes of previously described hepadnaviruses from mammalian and avian hosts and supports the creation of a novel genus. The discovery of this novel virus may better our understanding of the evolutionary history of hepatitis B-like viruses of other hosts. In fishes, knowledge of this virus may provide insight regarding possible risk factors associated with hepatic neoplasia in the white sucker. This may also offer another model system for mechanistic research.