Endogenous hepadnaviruses,bornaviruses and circoviruses in snakes

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

Low frequency of paleoviral infiltration across the avian phylogeny

Background

Mammalian genomes commonly harbor endogenous viral elements. Due to a lack of comparable genome-scale sequence data, far less is known about endogenous viral elements in avian species, even though their small genomes may enable important insights into the patterns and processes of endogenous viral element evolution.

Results

Through a systematic screening of the genomes of 48 species sampled across the avian phylogeny we reveal that birds harbor a limited number of endogenous viral elements compared to mammals, with only five viral families observed: Retroviridae, Hepadnaviridae, Bornaviridae, Circoviridae, and Parvoviridae. All nonretroviral endogenous viral elements are present at low copy numbers and in few species, with only endogenous hepadnaviruses widely distributed, although these have been purged in some cases. We also provide the first evidence for endogenous bornaviruses and circoviruses in avian genomes, although at very low copy numbers. A comparative analysis of vertebrate genomes revealed a simple linear relationship between endogenous viral element abundance and host genome size, such that the occurrence of endogenous viral elements in bird genomes is 6- to 13-fold less frequent than in mammals.

Conclusions

These results reveal that avian genomes harbor relatively small numbers of endogenous viruses, particularly those derived from RNA viruses, and hence are either less susceptible to viral invasions or purge them more effectively.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0539-3) contains supplementary material, which is available to authorized users.     

Early Mesozoic Coexistence of Amniotes and Hepadnaviridae

Hepadnaviridae are double-stranded DNA viruses that infect some species of birds and mammals. This includes humans, where hepatitis B viruses (HBVs) are prevalent pathogens in considerable parts of the global population. Recently, endogenized sequences of HBVs (eHBVs) have been discovered in bird genomes where they constitute direct evidence for the coexistence of these viruses and their hosts from the late Mesozoic until present. Nevertheless, virtually nothing is known about the ancient host range of this virus family in other animals. Here we report the first eHBVs from crocodilian, snake, and turtle genomes, including a turtle eHBV that endogenized >207 million years ago. This genomic “fossil” is >125 million years older than the oldest avian eHBV and provides the first direct evidence that Hepadnaviridae already existed during the Early Mesozoic. This implies that the Mesozoic fossil record of HBV infection spans three of the five major groups of land vertebrates, namely birds, crocodilians, and turtles. We show that the deep phylogenetic relationships of HBVs are largely congruent with the deep phylogeny of their amniote hosts, which suggests an ancient amniote–HBV coexistence and codivergence, at least since the Early Mesozoic. Notably, the organization of overlapping genes as well as the structure of elements involved in viral replication has remained highly conserved among HBVs along that time span, except for the presence of the X gene. We provide multiple lines of evidence that the tumor-promoting X protein of mammalian HBVs lacks a homolog in all other hepadnaviruses and propose a novel scenario for the emergence of X via segmental duplication and overprinting of pre-existing reading frames in the ancestor of mammalian HBVs. Our study reveals an unforeseen host range of prehistoric HBVs and provides novel insights into the genome evolution of hepadnaviruses throughout their long-lasting association with amniote hosts.     

Endogenous hepadnaviruses in the genome of the budgerigar (Melopsittacus undulatus) and the evolution of avian hepadnaviruses

Endogenous hepadnaviruses (hepatitis B viruses [HBVs]) were recently discovered in the genomes of passerine birds. We mined six additional avian genomes and discovered multiple copies of endogenous HBVs in the budgerigar (order Psittaciformes), designated eBHBV. A phylogenetic analysis reveals that the endogenous hepadnaviruses are more diverse than their exogenous counterparts and that the endogenous and exogenous hepadnaviruses form distinct lineages even when sampled from the same avian order, indicative of multiple genomic integration events.     

Distinct families of cis-acting RNA replication elements epsilon from hepatitis B viruses

The hepadnavirus encapsidation signal, epsilon (ε), is an RNA structure located at the 5' end of the viral pregenomic RNA. It is essential for viral replication and functions in polymerase protein binding and priming. This structure could also have potential regulatory roles in controlling the expression of viral replicative proteins. In addition to its structure, the primary sequence of this RNA element has crucial functional roles in the viral lifecycle. Although the ε elements in hepadnaviruses share common critical functions, there are some significant differences in mammalian and avian hepadnaviruses, which include both sequence and structural variations. Here we present several covariance models for ε elements from the Hepadnaviridae. The model building included experimentally determined data from previous studies using chemical probing and NMR analysis. These models have sufficient similarity to comprise a clan. The clan has in common a highly conserved overall structure consisting of a lower-stem, bulge, upper-stem and apical-loop. The models differ in functionally critical regions-notably the two types of avian ε elements have a tetra-loop (UGUU) including a non-canonical UU base pair, while the hepatitis B virus (HBV) epsilon has a tri-loop (UGU). The avian epsilon elements have a less stable dynamic structure in the upper stem. Comparisons between these models and all other Rfam models, and searches of genomes, showed these structures are specific to the Hepadnaviridae. Two family models and the clan are available from the Rfam database.     

Isolation and characterization of a hepatitis B virus endemic in herons.

A new hepadnavirus (designated heron hepatitis B virus [HHBV]) has been isolated; this virus is endemic in grey herons (Ardea cinerea) in Germany and closely related to duck hepatitis B virus (DHBV) by morphology of viral particles and size of the genome and of the major viral envelope and core proteins. Despite its striking similarities to DHBV, HHBV cannot be transmitted to ducks by infection or by transfection with cloned viral DNA. After the viral genome was cloned and sequenced, a comparative sequence analysis revealed an identical genome organization of HHBV and DHBV (pre-C/C-, pre-S/S-, and pol-ORFs). An open reading frame, designated X in mammalian hepadnaviruses, is not present in DHBV. DHBV and HHBV differ by 21.6% base exchanges, and thus they are less closely related than the two known rodent hepatitis B viruses (16.4%). The nucleocapsid protein and the 17-kilodalton envelope protein sequences of DHBV and HHBV are well conserved. In contrast, the pre-S part of the 34-kilodalton envelope protein which is believed to mediate virus attachment to the cell is highly divergent (less than 50% homology). The availability of two closely related avian hepadnaviruses will now allow us to test recombinant viruses in vivo and in vitro for host specificity-determining sequences.     

Endogenous bornavirus elements in mammalian genome

Approximately 8% of our genome is made up of endogenous retroviral elements. Endogenous retrovirus is a fossil record of ancient retrovirus infection and, therefore, gives important insights into the evolutional relationship between retroviruses and their hosts. On the other hand, until recently, it has been believed that no endogenous non-retroviral viruses exist in animal genomes. We lately discovered endogenous elements homologous to the nucleoprotein of bornaviruses, a negative-strand RNA virus, in the genomes of many mammalian species, including humans. We also demonstrated that mRNA of extant mammalian bornavirus, Borna disease virus, is reverse-transcribed and integrated into the host genome DNA. These findings provided novel insights not only into the interaction between RNA viruses and their hosts, but also into the mechanism underlying the gain of novelty in mammalian genomes. In this review, we will briefly summarize our recent knowledge about endogenous bornavirus elements and also introduce some recent discoveries regarding endogenous elements of non-retroviral viruses in vertebrate genomes.     

Phylogeny of <Emphasis Type="Italic">Banana Streak Virus</Emphasis> Reveals Recent and Repetitive Endogenization in the Genome of Its Banana Host (<Emphasis Type="Italic">Musa</Emphasis> sp.)

Banana streak virus (BSV) is a plant dsDNA pararetrovirus (family Caulimoviridae, genus badnavirus). Although integration is not an essential step in the BSV replication cycle, the nuclear genome of banana (Musa sp.) contains BSV endogenous pararetrovirus sequences (BSV EPRVs). Some BSV EPRVs are infectious by reconstituting a functional viral genome. Recent studies revealed a large molecular diversity of episomal BSV viruses (i.e., nonintegrated) while others focused on BSV EPRV sequences only. In this study, the evolutionary history of badnavirus integration in banana was inferred from phylogenetic relationships between BSV and BSV EPRVs. The relative evolution rates and selective pressures (d_N/d_S ratio) were also compared between endogenous and episomal viral sequences. At least 27 recent independent integration events occurred after the divergence of three banana species, indicating that viral integration is a recent and frequent phenomenon. Relaxation of selective pressure on badnaviral sequences that experienced neutral evolution after integration in the plant genome was recorded. Additionally, a significant decrease (35%) in the EPRV evolution rate was observed compared to BSV, reflecting the difference in the evolution rate between episomal dsDNA viruses and plant genome. The comparison of our results with the evolution rate of the Musa genome and other reverse-transcribing viruses suggests that EPRVs play an active role in episomal BSV diversity and evolution.     

Extant‐only comparative methods fail to recover the disparity preserved in the bird fossil record

Most extant species are in clades with poor fossil records, and recent studies of comparative methods show they have low power to infer even highly simplified models of trait evolution without fossil data. Birds are a well‐studied radiation, yet their early evolutionary patterns are still contentious. The fossil record suggests that birds underwent a rapid ecological radiation after the end‐Cretaceous mass extinction, and several smaller, subsequent radiations. This hypothesized series of repeated radiations from fossil data is difficult to test using extant data alone. By uniting morphological and phylogenetic data on 604 extant genera of birds with morphological data on 58 species of extinct birds from 50 million years ago, the “halfway point” of avian evolution, I have been able to test how well extant‐only methods predict the diversity of fossil forms. All extant‐only methods underestimate the disparity, although the ratio of within‐ to between‐clade disparity does suggest high early rates. The failure of standard models to predict high early disparity suggests that recent radiations are obscuring deep time patterns in the evolution of birds. Metrics from different models can be used in conjunction to provide more valuable insights than simply finding the model with the highest relative fit.     

An Endogenous Foamy-like Viral Element in the Coelacanth Genome

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.     

Unexpected Inheritance: Multiple Integrations of Ancient Bornavirus and Ebolavirus/Marburgvirus Sequences in Vertebrate Genomes

Vertebrate genomes contain numerous copies of retroviral sequences, acquired over the course of evolution. Until recently they were thought to be the only type of RNA viruses to be so represented, because integration of a DNA copy of their genome is required for their replication. In this study, an extensive sequence comparison was conducted in which 5,666 viral genes from all known non-retroviral families with single-stranded RNA genomes were matched against the germline genomes of 48 vertebrate species, to determine if such viruses could also contribute to the vertebrate genetic heritage. In 19 of the tested vertebrate species, we discovered as many as 80 high-confidence examples of genomic DNA sequences that appear to be derived, as long ago as 40 million years, from ancestral members of 4 currently circulating virus families with single strand RNA genomes. Surprisingly, almost all of the sequences are related to only two families in the Order Mononegavirales: the Bornaviruses and the Filoviruses, which cause lethal neurological disease and hemorrhagic fevers, respectively. Based on signature landmarks some, and perhaps all, of the endogenous virus-like DNA sequences appear to be LINE element-facilitated integrations derived from viral mRNAs. The integrations represent genes that encode viral nucleocapsid, RNA-dependent-RNA-polymerase, matrix and, possibly, glycoproteins. Integrations are generally limited to one or very few copies of a related viral gene per species, suggesting that once the initial germline integration was obtained (or selected), later integrations failed or provided little advantage to the host. The conservation of relatively long open reading frames for several of the endogenous sequences, the virus-like protein regions represented, and a potential correlation between their presence and a species'' resistance to the diseases caused by these pathogens, are consistent with the notion that their products provide some important biological advantage to the species. In addition, the viruses could also benefit, as some resistant species (e.g. bats) may serve as natural reservoirs for their persistence and transmission. Given the stringent limitations imposed in this informatics search, the examples described here should be considered a low estimate of the number of such integration events that have persisted over evolutionary time scales. Clearly, the sources of genetic information in vertebrate genomes are much more diverse than previously suspected.     

Phylogenetic analysis of pelecaniformes (aves) based on osteological data: implications for waterbird phylogeny and fossil calibration studies

Background

Debate regarding the monophyly and relationships of the avian order Pelecaniformes represents a classic example of discord between morphological and molecular estimates of phylogeny. This lack of consensus hampers interpretation of the group''s fossil record, which has major implications for understanding patterns of character evolution (e.g., the evolution of wing-propelled diving) and temporal diversification (e.g., the origins of modern families). Relationships of the Pelecaniformes were inferred through parsimony analyses of an osteological dataset encompassing 59 taxa and 464 characters. The relationships of the Plotopteridae, an extinct family of wing-propelled divers, and several other fossil pelecaniforms (Limnofregata, Prophaethon, Lithoptila, ?Borvocarbo stoeffelensis) were also assessed. The antiquity of these taxa and their purported status as stem members of extant families makes them valuable for studies of higher-level avian diversification.

Methodology/Principal Findings

Pelecaniform monophyly is not recovered, with Phaethontidae recovered as distantly related to all other pelecaniforms, which are supported as a monophyletic Steganopodes. Some anatomical partitions of the dataset possess different phylogenetic signals, and partitioned analyses reveal that these discrepancies are localized outside of Steganopodes, and primarily due to a few labile taxa. The Plotopteridae are recovered as the sister taxon to Phalacrocoracoidea, and the relationships of other fossil pelecaniforms representing key calibration points are well supported, including Limnofregata (sister taxon to Fregatidae), Prophaethon and Lithoptila (successive sister taxa to Phaethontidae), and ?Borvocarbo stoeffelensis (sister taxon to Phalacrocoracidae). These relationships are invariant when ‘backbone’ constraints based on recent avian phylogenies are imposed.

Conclusions/Significance

Relationships of extant pelecaniforms inferred from morphology are more congruent with molecular phylogenies than previously assumed, though notable conflicts remain. The phylogenetic position of the Plotopteridae implies that wing-propelled diving evolved independently in plotopterids and penguins, representing a remarkable case of convergent evolution. Despite robust support for the placement of fossil taxa representing key calibration points, the successive outgroup relationships of several “stem fossil + crown family” clades are variable and poorly supported across recent studies of avian phylogeny. Thus, the impact these fossils have on inferred patterns of temporal diversification depends heavily on the resolution of deep nodes in avian phylogeny.     

Fossils of parasites: what can the fossil record tell us about the evolution of parasitism?

Parasites are common in many ecosystems, yet because of their nature, they do not fossilise readily and are very rare in the geological record. This makes it challenging to study the evolutionary transition that led to the evolution of parasitism in different taxa. Most studies on the evolution of parasites are based on phylogenies of extant species that were constructed based on morphological and molecular data, but they give us an incomplete picture and offer little information on many important details of parasite–host interactions. The lack of fossil parasites also means we know very little about the roles that parasites played in ecosystems of the past even though it is known that parasites have significant influences on many ecosystems. The goal of this review is to bring attention to known fossils of parasites and parasitism, and provide a conceptual framework for how research on fossil parasites can develop in the future. Despite their rarity, there are some fossil parasites which have been described from different geological eras. These fossils include the free‐living stage of parasites, parasites which became fossilised with their hosts, parasite eggs and propagules in coprolites, and traces of pathology inflicted by parasites on the host's body. Judging from the fossil record, while there were some parasite–host relationships which no longer exist in the present day, many parasite taxa which are known from the fossil record seem to have remained relatively unchanged in their general morphology and their patterns of host association over tens or even hundreds of millions of years. It also appears that major evolutionary and ecological transitions throughout the history of life on Earth coincided with the appearance of certain parasite taxa, as the appearance of new host groups also provided new niches for potential parasites. As such, fossil parasites can provide additional data regarding the ecology of their extinct hosts, since many parasites have specific life cycles and transmission modes which reflect certain aspects of the host's ecology. The study of fossil parasites can be conducted using existing techniques in palaeontology and palaeoecology, and microscopic examination of potential material such as coprolites may uncover more fossil evidence of parasitism. However, I also urge caution when interpreting fossils as examples of parasites or parasitism‐induced traces. I point out a number of cases where parasitism has been spuriously attributed to some fossil specimens which, upon re‐examination, display traits which are just as (if not more) likely to be found in free‐living taxa. The study of parasite fossils can provide a more complete picture of the ecosystems and evolution of life throughout Earth's history.     

Surface antigenic determinants of mammalian "hepadnaviruses" defined by group- and class-specific monoclonal antibodies

The hepatitis B-like viruses (human hepatitis B virus, woodchuck hepatitis virus, ground squirrel hepatitis virus, and duck hepatitis B virus) are hepatotropic DNA viruses which have been referred to collectively as "hepadnaviruses." Using a murine monoclonal antibody (101-2) to the surface antigen of woodchuck hepatitis virus, we have shown that the surface antigens of mammalian hepadnaviruses (HBsAg, WHsAg, and GSHsAg) are antigenically related via a common determinant (HV/101). Furthermore, analysis with other monoclonal antibodies to WHsAg revealed that WHsAg and GHsAg are antigenically distinct, although the antigens had more determinants in common with each other than with HBsAg. The hepadnavirus group-specific antibody (101-2) reacted with HBsAg subtypic variants in a group-specific rather than subtype-specific manner. In conjunction with observations with an HBsAg-specific, group-reactive monoclonal antibody (BX259), the present data suggest that there are at least two group-reactive epitopes of HBsAg: one which is virus specific (HBV/259) and one which is common to two other mammalian hepadnaviruses (HV/101).     

Why coelacanths are not ‘living fossils’

A series of recent studies on extant coelacanths has emphasised the slow rate of molecular and morphological evolution in these species. These studies were based on the assumption that a coelacanth is a ‘living fossil’ that has shown little morphological change since the Devonian, and they proposed a causal link between low molecular evolutionary rate and morphological stasis. Here, we have examined the available molecular and morphological data and show that: (i) low intra‐specific molecular diversity does not imply low mutation rate, (ii) studies not showing low substitution rates in coelacanth are often neglected, (iii) the morphological stability of coelacanths is not supported by paleontological evidence. We recall that intra‐species levels of molecular diversity, inter‐species genome divergence rates and morphological divergence rates are under different constraints and they are not necessarily correlated. Finally, we emphasise that concepts such as ‘living fossil’, ‘basal lineage’, or ‘primitive extant species’ do not make sense from a tree‐thinking perspective. Editor's suggested further reading in BioEssays Tree thinking for all biology: the problem with reading phylogenies as ladders of progress Abstract     

Morphological Integration and the Interpretation of Fossil Hominin Diversity

The fossil record of primate and human evolution cannot provide accurate estimates of within species variation and integration. This means that we cannot directly observe how patterns of integration have evolved over time in this lineage. And yet, our interpretations of fossil diversity are awash with assumptions about variation patterning in precisely these fossil taxa. Most commonly, researchers rely on extant models of variation for interpreting past diversity, by assuming equality of variation (and occasionally covariation) among extant and fossil populations. Yet one of the things we know from studies of integration in primates is that patterns of morphological covariation can differ among even closely related taxa, indicating that they have diverged over evolutionary time, either in response to selection or as the result of neutral evolution. At the same time, overall patterns of integration remain remarkably similar, meaning that in many respects they are highly conserved evolutionarily. Taken together, these seemingly contradictory observations offer an important conceptual framework for interpreting patterns that we observe in the fossil past. This framework dictates that while we can use patterns of covariation in extant taxa as proxies for extinct diversity, and indeed their conserved nature makes them superior to approaches that rely on variation alone, we also need to account for the fact that such patterns change over time, and incorporate that into our models. Here I provide examples using covariation patterns estimated from modern humans and African great apes to demonstrate the extent to which divergence in covariance structure might affect our interpretations of hominin diversity.

Protein X of hepatitis B virus: origin and structure similarity with the central domain of DNA glycosylase

Orthohepadnavirus (mammalian hosts) and avihepadnavirus (avian hosts) constitute the family of Hepadnaviridae and differ by their capability and inability for expression of protein X, respectively. Origin and functions of X are unclear. The evolutionary analysis at issue of X indicates that present strains of orthohepadnavirus started to diverge about 25,000 years ago, simultaneously with the onset of avihepadnavirus diversification. These evolutionary events were preceded by a much longer period during which orthohepadnavirus developed a functional protein X while avihepadnavirus evolved without X. An in silico generated 3D-model of orthohepadnaviral X protein displayed considerable similarity to the tertiary structure of DNA glycosylases (key enzymes of base excision DNA repair pathways). Similarity is confined to the central domain of MUG proteins with the typical DNA-binding facilities but without the capability of DNA glycosylase enzymatic activity. The hypothetical translation product of a vestigial X reading frame in the genome of duck hepadnavirus could also been folded into a DNA glycosylase-like 3D-structure. In conclusion, the most recent common ancestor of ortho- and avihepadnavirus carried an X sequence with orthology to the central domain of DNA glycosylase.