Characterization of a Xenopus tropicalis endogenous retrovirus with developmental and stress-dependent expression

The Stress Response of the Holothurian Central Nervous System: A Transcriptomic Analysis

Injury to the central nervous system (CNS) results in permanent damage and lack of function in most vertebrate animals, due to their limited regenerative capacities. In contrast, echinoderms can fully regenerate their radial nerve cord (RNC) following transection, with little to no scarring. Investigators have associated the regenerative capacity of some organisms to the stress response and inflammation produced by the injury. Here, we explore the gene activation profile of the stressed holothurian CNS. To do this, we performed RNA sequencing on isolated RNC explants submitted to the stress of transection and enzyme dissection and compared them with explants kept in culture for 3 days following dissection. We describe stress-associated genes, including members of heat-shock families, ubiquitin-related pathways, transposons, and apoptosis that were differentially expressed. Surprisingly, the stress response does not induce apoptosis in this system. Other genes associated with stress in other animal models, such as hero proteins and those associated with the integrated stress response, were not found to be differentially expressed either. Our results provide a new viewpoint on the stress response in the nervous system of an organism with amazing regenerative capacities. This is the first step in deciphering the molecular processes that allow echinoderms to undergo fully functional CNS regeneration, and also provides a comparative view of the stress response in other organisms.

Revealing the uncharacterised diversity of amphibian and reptile viruses

Amphibians and non-avian reptiles represent a significant proportion of terrestrial vertebrates, however knowledge of their viruses is not proportional to their abundance. Many amphibians and reptiles have strict habitual environments and localised populations and are vulnerable to viral outbreaks and potential elimination as a result. We sought to identify viruses that were hidden in amphibian and reptile metatranscriptomic data by screening 235 RNA-sequencing datasets from a 122 species covering 25 countries. We identified 26 novel viruses and eight previously characterised viruses from fifteen different viral families. Twenty-five viruses had RNA genomes with identity to Arteriviridae, Tobaniviridae, Hantaviridae, Rhabdoviridae, Astroviridae, Arenaviridae, Hepeviridae, Picornaviridae, Orthomyxoviridae, Reoviridae, Flaviviridae and Caliciviridae. In addition to RNA viruses, we also screened datasets for DNA viral transcripts, which are commonly excluded from transcriptomic analysis. We identified ten DNA viruses with identity to Papillomaviridae, Parvoviridae, Circoviridae and Adomaviridae. With the addition of these viruses, we expand the global amphibian and reptile virome and identify new potentially pathogenic viruses that could challenge populations. We speculate that amphibian viruses often have simpler genomes than those in amniotes, as in the case of the Secondpapillomavirinae and Orthomyxoviridae viruses identified in this study. In addition, we find evidence of inter-family recombination in RNA viruses, and we also identify new members of the recombinant Adomaviridae family. Overall, we provide insights into the uncharacterised diversity of amphibian and reptile viruses with the aim of improving population management, treatment and conservation into the future.

A novel class III endogenous retrovirus with a class I envelope gene in African frogs with an intact genome and developmentally regulated transcripts in Xenopus tropicalis

BACKGROUND

Retroviruses exist as exogenous infectious agents and as endogenous retroviruses (ERVs) integrated into host chromosomes. Such endogenous retroviruses (ERVs) are grouped into three classes roughly corresponding to the seven genera of infectious retroviruses: class I (gamma-, epsilonretroviruses), class II (alpha-, beta-, delta-, lentiretroviruses) and class III (spumaretroviruses). Some ERVs have counterparts among the known infectious retroviruses, while others represent paleovirological relics of extinct or undiscovered retroviruses.

RESULTS

Here we identify an intact ERV in the Anuran amphibian, Xenopus tropicalis. XtERV-S has open reading frames (ORFs) for gag, pol (polymerase) and env (envelope) genes, with a small additional ORF in pol and a serine tRNA primer binding site. It has unusual features and domain relationships to known retroviruses. Analyses based on phylogeny and functional motifs establish that XtERV-S gag and pol genes are related to the ancient env-less class III ERV-L family but the surface subunit of env is unrelated to known retroviruses while its transmembrane subunit is class I-like. LTR constructs show transcriptional activity, and XtERV-S transcripts are detected in embryos after the maternal to zygotic mid-blastula transition and before the late tailbud stage. Tagged Gag protein shows typical subcellular localization. The presence of ORFs in all three protein-coding regions along with identical 5' and 3' LTRs (long terminal repeats) indicate this is a very recent germline acquisition. There are older, full-length, nonorthologous, defective copies in Xenopus laevis and the distantly related African bullfrog, Pyxicephalus adspersus. Additional older, internally deleted copies in X. tropicalis carry a 300 bp LTR substitution.

CONCLUSIONS

XtERV-S represents a genera-spanning member of the largely env-less class III ERV that has ancient and modern copies in Anurans. This provirus has an env ORF with a surface subunit unrelated to known retroviruses and a transmembrane subunit related to class I gammaretroviruses in sequence and organization, and is expressed in early embryogenesis. Additional XtERV-S-related but defective copies are present in X. tropicalis and other African frog taxa. XtERV-S is an unusual class III ERV variant, and it may represent an important transitional retroviral form that has been spreading in African frogs for tens of millions of years.

Viral Discovery in the Invasive Australian Cane Toad (Rhinella marina) Using Metatranscriptomic and Genomic Approaches

Cane toads are poisonous amphibians that were introduced to Australia in 1935 for insect control. Since then, their population has increased dramatically, and they now threaten many native Australian species. One potential method to control the population is to release a cane toad virus with high mortality rates, yet few cane toad viruses have been characterized. This study samples cane toads from different Australian locations and uses an RNA sequencing and computational approach to find new viruses. We report novel complete picornavirus and retrovirus sequences that were genetically similar to viruses infecting frogs, reptiles, and fish. Using data generated in other studies, we show that these viral sequences are present in cane toads from distinct Australian locations. Three sequences related to circoviruses were also found in the toad genome. The identification of new viral sequences will aid future studies that investigate their prevalence and potential as agents for biocontrol.

Cane toads are a notorious invasive species, inhabiting over 1.2 million km² of Australia and threatening native biodiversity. The release of pathogenic cane toad viruses is one possible biocontrol strategy yet is currently hindered by the poorly described cane toad virome. Metatranscriptomic analysis of 16 cane toad livers revealed the presence of a novel and full-length picornavirus, Rhimavirus A (RhiV-A), a member of a reptile- and amphibian-specific cluster of the Picornaviridae basal to the Kobuvirus-like group. In the combined liver transcriptome, we also identified a complete genome sequence of a distinct epsilonretrovirus, Rhinella marina endogenous retrovirus (RMERV). The recently sequenced cane toad genome contains 8 complete RMERV proviruses as well as 21 additional truncated insertions. The oldest full-length RMERV provirus was estimated to have inserted 1.9 million years ago (MYA). To screen for these viral sequences in additional toads, we analyzed publicly available transcriptomes from six diverse Australian locations. RhiV-A transcripts were identified in toads sampled from three locations across 1,000 km of Australia, stretching to the current Western Australia (WA) invasion front, while RMERV transcripts were observed at all six sites. Finally, we scanned the cane toad genome for nonretroviral endogenous viral elements, finding three sequences related to small DNA viruses in the family Circoviridae. This shows ancestral circoviral infection with subsequent genomic integration. The identification of these current and past viral infections enriches our knowledge of the cane toad virome, an understanding of which will facilitate future work on infection and disease in this important invasive species.

IMPORTANCE Cane toads are poisonous amphibians that were introduced to Australia in 1935 for insect control. Since then, their population has increased dramatically, and they now threaten many native Australian species. One potential method to control the population is to release a cane toad virus with high mortality rates, yet few cane toad viruses have been characterized. This study samples cane toads from different Australian locations and uses an RNA sequencing and computational approach to find new viruses. We report novel complete picornavirus and retrovirus sequences that were genetically similar to viruses infecting frogs, reptiles, and fish. Using data generated in other studies, we show that these viral sequences are present in cane toads from distinct Australian locations. Three sequences related to circoviruses were also found in the toad genome. The identification of new viral sequences will aid future studies that investigate their prevalence and potential as agents for biocontrol.

Co-option of endogenous viral sequences for host cell function

Eukaryotic genomes are littered with sequences of diverse viral origins, termed endogenous viral elements (EVEs). Here we used examples primarily drawn from mammalian endogenous retroviruses to document how the influx of EVEs has provided a source of prefabricated coding and regulatory sequences that were formerly utilized for viral infection and replication, but have been occasionally repurposed for cellular function. While EVE co-option has benefited a variety of host biological functions, there appears to be a disproportionate contribution to immunity and antiviral defense. The mammalian embryo and placenta offer opportunistic routes of viral transmission to the next host generation and as such they represent hotbeds for EVE cooption. Based on these observations, we propose that EVE cooption is initially driven as a mean to mitigate conflicts between host and viruses, which in turn acts as a stepping-stone toward the evolution of cellular innovations serving host physiology and development.

Origin of the retroviruses: when, where, and how?

•

Retroviruses are a virus family of considerable medical and veterinary importance.

•

Until recently, very little was known about deep retroviral origins.

•

New research supports a marine origin of retroviruses, ∼460–550 million years ago.

•

The evolutionary events leading to the origin of retroviruses remain obscure.

•

Improved understanding of Metaviridae diversity and evolution are required for this.

Retroviruses are a virus family of considerable medical and veterinary importance. Additionally, it is now clear that endogenous retroviruses (ERVs) comprise significant portions of vertebrate genomes. Until recently, very little was known about the deep evolutionary origins of retroviruses. However, advances in genomics and bioinformatics have opened the way for great strides in understanding. Recent research employing a wide variety of bioinformatic approaches has demonstrated that retroviruses evolved during the early Palaeozoic Era, between 460 and 550 million years ago, providing the oldest inferred date estimate for any virus group. This finding presents an important framework to investigate the evolutionary transitions that led to the emergence of the retroviruses, offering potential insights into the infectious origins of a major group of pathogenic viruses.

The contribution of transposable elements to size variations between four teleost genomes

Background

Teleosts are unique among vertebrates, with a wide range of haploid genome sizes in very close lineages, varying from less than 400 mega base pairs (Mb) for pufferfish to over 3000 Mb for salmon. The cause of the difference in genome size remains largely unexplained.

Results

In this study, we reveal that the differential success of transposable elements (TEs) correlates with the variation of genome size across four representative teleost species (zebrafish, medaka, stickleback, and tetraodon). The larger genomes represent a higher diversity within each clade (superfamily) and family and a greater abundance of TEs compared with the smaller genomes; zebrafish, representing the largest genome, shows the highest diversity and abundance of TEs in its genome, followed by medaka and stickleback; while the tetraodon, representing the most compact genome, displays the lowest diversity and density of TEs in its genome. Both of Class I (retrotransposons) and Class II TEs (DNA transposons) contribute to the difference of TE accumulation of teleost genomes, however, Class II TEs are the major component of the larger teleost genomes analyzed and the most important contributors to genome size variation across teleost lineages. The hAT and Tc1/Mariner superfamilies are the major DNA transposons of all four investigated teleosts. Divergence distribution revealed contrasting proliferation dynamics both between clades of retrotransposons and between species. The TEs within the larger genomes of the zebrafish and medaka represent relatively stronger activity with an extended time period during the evolution history, in contrast with the very young activity in the smaller stickleback genome, or the very low level of activity in the tetraodon genome.

Conclusion

Overall, our data shows that teleosts represent contrasting profiles of mobilomes with a differential density, diversity and activity of TEs. The differences in TE accumulation, dominated by DNA transposons, explain the main size variations of genomes across the investigated teleost species, and the species differences in both diversity and activity of TEs contributed to the variations of TE accumulations across the four teleost species. TEs play major roles in teleost genome evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s13100-016-0059-7) contains supplementary material, which is available to authorized users.

Endogenous viruses: Connecting recent and ancient viral evolution

The rapid rates of viral evolution allow us to reconstruct the recent history of viruses in great detail. This feature, however, also results in rapid erosion of evolutionary signal within viral molecular data, impeding studies of their deep history. Thus, the further back in time, the less accurate the inference becomes. Furthermore, reconstructing complex histories of transmission can be challenging, especially where extinct viral lineages are concerned. This problem has been partially solved by the discovery of viruses embedded in host genomes, known as endogenous viral elements (EVEs). Some of these endogenous viruses are derived from ancient relatives of extant viruses, allowing us to better examine ancient viral host range, geographical distribution and transmission routes. Moreover, our knowledge of viral evolutionary timescales and rate dynamics has also been greatly improved by their discovery, thereby bridging the gap between recent and ancient viral evolution.

Insights on genome size evolution from a miniature inverted repeat transposon driving a satellite DNA

The genome size in eukaryotes does not correlate well with the number of genes they contain. We can observe this so-called C-value paradox in amphibian species. By analyzing an amphibian genome we asked how repetitive DNA can impact genome size and architecture. We describe here our discovery of a Tc1/mariner miniature inverted-repeat transposon family present in Xenopus frogs. These transposons named miDNA4 are unique since they contain a satellite DNA motif. We found that miDNA4 measured 331 bp, contained 25 bp long inverted terminal repeat sequences and a sequence motif of 119 bp present as a unique copy or as an array of 2-47 copies. We characterized the structure, dynamics, impact and evolution of the miDNA4 family and its satellite DNA in Xenopus frog genomes. This led us to propose a model for the evolution of these two repeated sequences and how they can synergize to increase genome size.

Multiple groups of endogenous epsilon-like retroviruses conserved across primates

Several types of cancer in fish are caused by retroviruses, including those responsible for major outbreaks of disease, such as walleye dermal sarcoma virus and salmon swim bladder sarcoma virus. These viruses form a phylogenetic group often described as the epsilonretrovirus genus. Epsilon-like retroviruses have become endogenous retroviruses (ERVs) on several occasions, integrating into germ line cells to become part of the host genome, and sections of fish and amphibian genomes are derived from epsilon-like retroviruses. However, epsilon-like ERVs have been identified in very few mammals. We have developed a pipeline to screen full genomes for ERVs, and using this pipeline, we have located over 800 endogenous epsilon-like ERV fragments in primate genomes. Genomes from 32 species of mammals and birds were screened, and epsilon-like ERV fragments were found in all primate and tree shrew genomes but no others. These viruses appear to have entered the genome of a common ancestor of Old and New World monkeys between 42 million and 65 million years ago. Based on these results, there is an ancient evolutionary relationship between epsilon-like retroviruses and primates. Clearly, these viruses had the potential to infect the ancestors of primates and were at some point a common pathogen in these hosts. Therefore, this result raises questions about the potential of epsilonretroviruses to infect humans and other primates and about the evolutionary history of these retroviruses.

IMPORTANCE

Epsilonretroviruses are a group of retroviruses that cause several important diseases in fish. Retroviruses have the ability to become a permanent part of the DNA of their host by entering the germ line as endogenous retroviruses (ERVs), where they lose their infectivity over time but can be recognized as retroviruses for millions of years. Very few mammals are known to have epsilon-like ERVs; however, we have identified over 800 fragments of endogenous epsilon-like ERVs in the genomes of all major groups of primates, including humans. These viruses seem to have circulated and infected primate ancestors 42 to 65 million years ago. We are now interested in how these viruses have evolved and whether they have the potential to infect modern humans or other primates.

LTR retroelements are intrinsic components of transcriptional networks in frogs

Background

LTR retroelements (LTR REs) constitute a major group of transposable elements widely distributed in eukaryotic genomes. Through their own mechanism of retrotranscription LTR REs enrich the genomic landscape by providing genetic variability, thus contributing to genome structure and organization. Nonetheless, transcriptomic activity of LTR REs still remains an obscure domain within cell, developmental, and organism biology.

Results

Here we present a first comparative analysis of LTR REs for anuran amphibians based on a full depth coverage transcriptome of the European pool frog, Pelophylax lessonae, the genome of the African clawed frog, Silurana tropicalis (release v7.1), and additional transcriptomes of S. tropicalis and Cyclorana alboguttata. We identified over 1000 copies of LTR REs from all four families (Bel/Pao, Ty1/Copia, Ty3/Gypsy, Retroviridae) in the genome of S. tropicalis and discovered transcripts of several of these elements in all RNA-seq datasets analyzed. Elements of the Ty3/Gypsy family were most active, especially Amn-san elements, which accounted for approximately 0.27% of the genome in Silurana. Some elements exhibited tissue specific expression patterns, for example Hydra1.1 and MuERV-like elements in Pelophylax. In S. tropicalis considerable transcription of LTR REs was observed during embryogenesis as soon as the embryonic genome became activated, i.e. at midblastula transition. In the course of embryonic development the spectrum of transcribed LTR REs changed; during gastrulation and neurulation MuERV-like and SnRV like retroviruses were abundantly transcribed while during organogenesis transcripts of the XEN1 retroviruses became much more active.

Conclusions

The differential expression of LTR REs during embryogenesis in concert with their tissue-specificity and the protein domains they encode are evidence for the functional roles these elements play as integrative parts of complex regulatory networks. Our results support the meanwhile widely accepted concept that retroelements are not simple “junk DNA” or “harmful genomic parasites” but essential components of the transcriptomic machinery in vertebrates.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-626) contains supplementary material, which is available to authorized users.

Characterisation of retroviruses in the horse genome and their transcriptional activity via transcriptome sequencing

The recently released draft horse genome is incompletely characterised in terms of its repetitive element profile. This paper presents characterisation of the endogenous retrovirus (ERVs) of the horse genome based on a data-mining strategy using murine leukaemia virus proteins as queries. 978 ERV gene sequences were identified. Sequences were identified from the gamma, epsilon and betaretrovirus genera. At least one full length gammaretroviral locus was identified, though the gammaretroviral sequences are very degenerate. Using these data the RNA expression of these ERVs were derived from RNA transcriptome data from a variety of equine tissues. Unlike the well studied human and murine ERVs there do not appear to be particular phylogenetic groups of equine ERVs that are more transcriptionally active. Using this novel approach provided a more technically feasible method to characterise ERV expression than previous studies.