The human zoo: endogenous retroviruses in the human genome

A recent gibbon ape leukemia virus germline integration in a rodent from New Guinea

Germline colonization by retroviruses results in the formation of endogenous retroviruses (ERVs). Most colonization's occurred millions of years ago. However, in the Australo-Papuan region (Australia and New Guinea), several recent germline colonization events have been discovered. The Wallace Line separates much of Southeast Asia from the Australo-Papuan region restricting faunal and pathogen dispersion. West of the Wallace Line, gibbon ape leukemia viruses (GALVs) have been isolated from captive gibbons. Two microbat species from China appear to have been infected naturally. East of Wallace's Line, the woolly monkey virus (a GALV) and the closely related koala retrovirus (KoRV) have been detected in eutherians and marsupials in the Australo-Papuan region, often vertically transmitted. The detected vertically transmitted GALV-like viruses in Australo-Papuan fauna compared to sporadic horizontal transmission in Southeast Asia and China suggest the GALV-KoRV clade originates in the former region and further models of early-stage genome colonization may be found. We screened 278 samples, seven bat and one rodent family endemic to the Australo-Papuan region and bat and rodent species found on both sides of the Wallace Line. We identified two rodents (Melomys) from Australia and Papua New Guinea and no bat species harboring GALV-like retroviruses. Melomys leucogaster from New Guinea harbored a genomically complete replication-competent retrovirus with a shared integration site among individuals. The integration was only present in some individuals of the species indicating this retrovirus is at the earliest stages of germline colonization of the Melomys genome, providing a new small wild mammal model of early-stage genome colonization.

A novel Australian flying-fox retrovirus shares an evolutionary ancestor with Koala, Gibbon and Melomys gamma-retroviruses

A novel gamma-retroviral sequence (7912 bp), inclusive of both partial 5' and 3' long terminal repeat regions, was identified from the brain of a black flying-fox (Pteropus alecto), Queensland, Australia. The sequence was distinct from other retroviral sequences identified in bats and showed greater identity to Koala, Gibbon ape leukaemia, Melomys burtoni and Woolly monkey retroviruses, forming their own phylogenetic clade. This finding suggests that these retroviruses may have an unknown common ancestor and that further investigation into the diversity of gamma-retroviruses in Australian Pteropus species may elucidate their evolutionary origins.

Degradation and remobilization of endogenous retroviruses by recombination during the earliest stages of a germ-line invasion

Endogenous retroviruses (ERVs) are proviral sequences that result from host germ-line invasion by exogenous retroviruses. The majority of ERVs are degraded. Using the koala retrovirus (KoRV) as a model system, we demonstrate that recombination with an ancient koala retroelement disables KoRV, and that recombination occurs frequently and early in the invasion process. Recombinant KoRVs (recKoRVs) are then able to proliferate in the koala germ line. This may in part explain the generally degraded nature of ERVs in vertebrate genomes and suggests that degradation via recombination is one of the earliest processes shaping retroviral genomic invasions.

Endogenous retroviruses (ERVs) are proviral sequences that result from colonization of the host germ line by exogenous retroviruses. The majority of ERVs represent defective retroviral copies. However, for most ERVs, endogenization occurred millions of years ago, obscuring the stages by which ERVs become defective and the changes in both virus and host important to the process. The koala retrovirus, KoRV, only recently began invading the germ line of the koala (Phascolarctos cinereus), permitting analysis of retroviral endogenization on a prospective basis. Here, we report that recombination with host genomic elements disrupts retroviruses during the earliest stages of germ-line invasion. One type of recombinant, designated recKoRV1, was formed by recombination of KoRV with an older degraded retroelement. Many genomic copies of recKoRV1 were detected across koalas. The prevalence of recKoRV1 was higher in northern than in southern Australian koalas, as is the case for KoRV, with differences in recKoRV1 prevalence, but not KoRV prevalence, between inland and coastal New South Wales. At least 15 additional different recombination events between KoRV and the older endogenous retroelement generated distinct recKoRVs with different geographic distributions. All of the identified recombinant viruses appear to have arisen independently and have highly disrupted ORFs, which suggests that recombination with existing degraded endogenous retroelements may be a means by which replication-competent ERVs that enter the germ line are degraded.

Koala immunology and infectious diseases: How much can the koala bear?

Infectious diseases are contributing to the decline of the iconic Australian marsupial, the koala (Phascolarctos cinereus). Infections with the obligate intracellular bacteria, Chlamydia pecorum, cause debilitating ocular and urogenital-tract disease while the koala-retrovirus (KoRV) has been implicated in host immunosuppression and exacerbation of chlamydial pathogenesis. Although histological studies have provided insight into the basic architecture of koala immune tissues, our understanding of the koala immune response to infectious disease has been limited, until recently, by a lack of species-specific immune reagents. Recent advances in the characterisation of key immune genes have focused on advancing our understanding of the immune response to Chlamydia infection, revealing commonalities in disease pathologies and immunity between koalas and other hosts and paving the way for the development of a koala Chlamydia vaccine. This review summarises these recent findings and highlights key aspects of the koala immune system requiring further attention with particular regard to their most prominent infectious diseases.

Transmission, Evolution, and Endogenization: Lessons Learned from Recent Retroviral Invasions

Viruses of the subfamily Orthoretrovirinae are defined by the ability to reverse transcribe an RNA genome into DNA that integrates into the host cell genome during the intracellular virus life cycle. Exogenous retroviruses (XRVs) are horizontally transmitted between host individuals, with disease outcome depending on interactions between the retrovirus and the host organism. When retroviruses infect germ line cells of the host, they may become endogenous retroviruses (ERVs), which are permanent elements in the host germ line that are subject to vertical transmission. These ERVs sometimes remain infectious and can themselves give rise to XRVs. This review integrates recent developments in the phylogenetic classification of retroviruses and the identification of retroviral receptors to elucidate the origins and evolution of XRVs and ERVs. We consider whether ERVs may recurrently pressure XRVs to shift receptor usage to sidestep ERV interference. We discuss how related retroviruses undergo alternative fates in different host lineages after endogenization, with koala retrovirus (KoRV) receiving notable interest as a recent invader of its host germ line. KoRV is heritable but also infectious, which provides insights into the early stages of germ line invasions as well as XRV generation from ERVs. The relationship of KoRV to primate and other retroviruses is placed in the context of host biogeography and the potential role of bats and rodents as vectors for interspecies viral transmission. Combining studies of extant XRVs and "fossil" endogenous retroviruses in koalas and other Australasian species has broadened our understanding of the evolution of retroviruses and host-retrovirus interactions.

Koala Retroviruses: Evolution and Disease Dynamics

A retroviral etiology for malignant neoplasias in koalas has long been suspected. Evidence for retroviral involvement was bolstered in 2000 by the isolation of a koala retrovirus (KoRV), now termed KoRV-A. KoRV-A is an endogenous retrovirus-a retrovirus that infects germ cells-a feature that makes it a permanent resident of the koala genome. KoRV-A lacks the genetic diversity of an exogenous retrovirus, a quality associated with the ability of a retrovirus to cause neoplasias. In 2013, a second KoRV isolate, KoRV-B, was obtained from koalas with lymphomas in the Los Angeles Zoo. Unlike KoRV-A, which is present in the genomes of all koalas in the United States, KoRV-B is restricted in its distribution and is associated with host pathology (neoplastic disease). Here, our current understanding of the evolution of endogenous and exogenous KoRVs, and the relationship between them, is reviewed to build a perspective on the future impact of these viruses on koala sustainability.

A new look at the LTR retrotransposon content of the chicken genome

Background

LTR retrotransposons contribute approximately 10 % of the mammalian genome, but it has been previously reported that there is a deficit of these elements in the chicken relative to both mammals and other birds. A novel LTR retrotransposon classification pipeline, LocaTR, was developed and subsequently utilised to re-examine the chicken LTR retrotransposon annotation, and determine if the proposed chicken deficit is biologically accurate or simply a technical artefact.

Results

Using LocaTR 3.01 % of the chicken galGal4 genome assembly was annotated as LTR retrotransposon-derived elements (nearly double the previous annotation), including 1,073 that were structurally intact. Element distribution is significantly correlated with chromosome size and is non-random within each chromosome. Elements are significantly depleted within coding regions and enriched in gene sparse areas of the genome. Over 40 % of intact elements are found in clusters, unrelated by age or genera, generally in poorly recombining regions. The transcription of most LTR retrotransposons were suppressed or incomplete, but individual domain and full length retroviral transcripts were produced in some cases, although mostly with regularly interspersed stop codons in all reading frames. Furthermore, RNAseq data from 23 diverse tissues enabled greater characterisation of the co-opted endogenous retrovirus Ovex1. This gene was shown to be expressed ubiquitously but at variable levels across different tissues. LTR retrotransposon content was found to be very variable across the avian lineage and did not correlate with either genome size or phylogenetic position. However, the extent of previous, species-specific LTR retrotransposon annotation appears to be a confounding factor.

Conclusions

Use of the novel LocaTR pipeline has nearly doubled the annotated LTR retrotransposon content of the chicken genome compared to previous estimates. Further analysis has described element distribution, clustering patterns and degree of expression in a variety of adult tissues, as well as in three embryonic stages. This study also enabled better characterisation of the co-opted gamma retroviral envelope gene Ovex1. Additionally, this work suggests that there is no deficit of LTR retrotransposons within the Galliformes relative to other birds, or to mammalian genomes when scaled for the three-fold difference in genome size.

Electronic supplementary material

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

Comprehensive profiling of retroviral integration sites using target enrichment methods from historical koala samples without an assembled reference genome

Background. Retroviral integration into the host germline results in permanent viral colonization of vertebrate genomes. The koala retrovirus (KoRV) is currently invading the germline of the koala (Phascolarctos cinereus) and provides a unique opportunity for studying retroviral endogenization. Previous analysis of KoRV integration patterns in modern koalas demonstrate that they share integration sites primarily if they are related, indicating that the process is currently driven by vertical transmission rather than infection. However, due to methodological challenges, KoRV integrations have not been comprehensively characterized.

Results. To overcome these challenges, we applied and compared three target enrichment techniques coupled with next generation sequencing (NGS) and a newly customized sequence-clustering based computational pipeline to determine the integration sites for 10 museum Queensland and New South Wales (NSW) koala samples collected between the 1870s and late 1980s. A secondary aim of this study sought to identify common integration sites across modern and historical specimens by comparing our dataset to previously published studies. Several million sequences were processed, and the KoRV integration sites in each koala were characterized.

Conclusions. Although the three enrichment methods each exhibited bias in integration site retrieval, a combination of two methods, Primer Extension Capture and hybridization capture is recommended for future studies on historical samples. Moreover, identification of integration sites shows that the proportion of integration sites shared between any two koalas is quite small.

Genetic diversity of koala retroviral envelopes

Genetic diversity, attributable to the low fidelity of reverse transcription, recombination and mutation, is an important feature of infectious retroviruses. Under selective pressure, such as that imposed by superinfection interference, gammaretroviruses commonly adapt their envelope proteins to use alternative receptors to overcome this entry block. The first characterized koala retroviruses KoRV subgroup A (KoRV-A) were remarkable in their absence of envelope genetic variability. Once it was determined that KoRV-A was present in all koalas in US zoos, regardless of their disease status, we sought to isolate a KoRV variant whose presence correlated with neoplastic malignancies. More than a decade after the identification of KoRV-A, we isolated a second subgroup of KoRV, KoRV-B from koalas with lymphomas. The envelope proteins of KoRV-A and KoRV-B are sufficiently divergent to confer the ability to bind and employ distinct receptors for infection. We have now obtained a number of additional KoRV envelope variants. In the present studies we report these variants, and show that they differ from KoRV-A and KoRV-B envelopes in their host range and superinfection interference properties. Thus, there appears to be considerable variation among KoRVs envelope genes suggesting genetic diversity is a factor following the KoRV-A infection process.

Sequence variation of koala retrovirus transmembrane protein p15E among koalas from different geographic regions

The koala retrovirus (KoRV), which is transitioning from an exogenous to an endogenous form, has been associated with high mortality in koalas. For other retroviruses, the envelope protein p15E has been considered a candidate for vaccine development. We therefore examined proviral sequence variation of KoRV p15E in a captive Queensland and three wild southern Australian koalas. We generated 163 sequences with intact open reading frames, which grouped into 39 distinct haplotypes. Sixteen distinct haplotypes comprising 139 of the sequences (85%) coded for the same polypeptide. Among the remaining 23 haplotypes, 22 were detected only once among the sequences, and each had 1 or 2 non-synonymous differences from the majority sequence. Several analyses suggested that p15E was under purifying selection. Important epitopes and domains were highly conserved across the p15E sequences and in previously reported exogenous KoRVs. Overall, these results support the potential use of p15E for KoRV vaccine development.

Proliferation of endogenous retroviruses in the early stages of a host germ line invasion

Endogenous retroviruses (ERVs) comprise 8% of the human genome and are common in all vertebrate genomes. The only retrovirus known to be currently transitioning from exogenous to endogenous form is the koala retrovirus (KoRV), making koalas (Phascolarctos cinereus) ideal for examining the early stages of retroviral endogenization. To distinguish endogenous from exogenous KoRV proviruses, we isolated koala genomic regions flanking KoRV integration sites. In three wild southern Australian koalas, there were fewer KoRV loci than in three captive Queensland koalas, consistent with reports that southern Australian koalas carry fewer KoRVs. Of 39 distinct KoRV proviral loci examined in a sire-dam-progeny triad, all proved to be vertically transmitted and endogenous; none was exogenous. Of the 39 endogenous KoRVs (enKoRVs), only one was present in the genomes of both the sire and the dam, suggesting that, at this early stage in the retroviral invasion of a host germ line, very large numbers of ERVs have proliferated at very low frequencies in the koala population. Sequence divergence between the 5'- and 3'-long terminal repeats (LTRs) of a provirus can be used as a molecular clock. Within each of ten enKoRVs, the 5'-LTR sequence was identical to the 3'-LTR sequence, suggesting a maximum age for enKoRV invasion of the koala germ line of approximately 22,200-49,900 years ago, although a much younger age is possible. Across the ten proviruses, seven LTR haplotypes were detected, indicating that at least seven different retroviral sequences had entered the koala germ line.

Discovery of a novel retrovirus sequence in an Australian native rodent (Melomys burtoni): a putative link between gibbon ape leukemia virus and koala retrovirus

Gibbon ape leukaemia virus (GALV) and koala retrovirus (KoRV) share a remarkably close sequence identity despite the fact that they occur in distantly related mammals on different continents. It has previously been suggested that infection of their respective hosts may have occurred as a result of a species jump from another, as yet unidentified vertebrate host. To investigate possible sources of these retroviruses in the Australian context, DNA samples were obtained from 42 vertebrate species and screened using PCR in order to detect proviral sequences closely related to KoRV and GALV. Four proviral partial sequences totalling 2880 bases which share a strong similarity with KoRV and GALV were detected in DNA from a native Australian rodent, the grassland melomys, Melomys burtoni. We have designated this novel gammaretrovirus Melomys burtoni retrovirus (MbRV). The concatenated nucleotide sequence of MbRV shares 93% identity with the corresponding sequence from GALV-SEATO and 83% identity with KoRV. The geographic ranges of the grassland melomys and of the koala partially overlap. Thus a species jump by MbRV from melomys to koalas is conceivable. However the genus Melomys does not occur in mainland South East Asia and so it appears most likely that another as yet unidentified host was the source of GALV.

Prevalence of koala retrovirus in geographically diverse populations in Australia

OBJECTIVE

To determine the prevalence of koala retrovirus (KoRV) in selected koala populations and to estimate proviral copy number in a subset of koalas.

METHODS

Blood or tissue samples from 708 koalas in Queensland, New South Wales, Victoria and South Australia were tested for KoRV pol provirus gene using standard polymerase chain reaction (PCR), nested PCR and real-time PCR (qPCR).

RESULTS

Prevalence of KoRV provirus-positive koalas was 100% in four regions of Queensland and New South Wales, 72.2% in mainland Victoria, 26.6% on four Victorian islands and 14.8% on Kangaroo Island, South Australia. Estimated proviral copy number per cell in four groups of koalas from Queensland and Victoria showed marked variation, ranging from a mean of 165 copies per cell in the Queensland group to 1.29 × 10(-4) copies per cell in one group of Victorian koalas.

CONCLUSIONS

The higher prevalence of KoRV-positive koalas in the north of Australia and high proviral loads in Queensland koalas may indicate KoRV entered and became endogenous in the north and is spreading southwards. It is also possible there are genetic differences between koalas in northern and southern Australia that affect susceptibility to KoRV infection or endogenisation, or that environmental factors affecting transmission in northern states are absent or uncommon in southern regions. Although further studies are required, the finding of proviral copy numbers orders of magnitude lower than what would be expected for the presence of a single copy in every cell for many Victorian animals suggests that KoRV is not endogenous in these animals and likely reflects ongoing exogenous infection.

Viruses are essential agents within the roots and stem of the tree of life

In contrast with former definitions of life limited to membrane-bound cellular life forms which feed, grow, metabolise and replicate (i) a role of viruses as genetic symbionts, (ii) along with peripheral phenomena such as cryptobiosis and (iii) the horizontal nature of genetic information acquisition and processing broaden our view of the tree of life. Some researchers insist on the traditional textbook conviction of what is part of the community of life. In a recent review [Moreira, D., Lopez-Garcia, P., 2009. Ten reasons to exclude viruses from the tree of life. Nat. Rev. Microbiol. 7, 306-311.] they assemble four main arguments which should exclude viruses from the tree of life because of their inability to self-sustain and self-replicate, their polyphyly, the cellular origin of their cell-like genes and the volatility of their genomes. In this article we will show that these features are not coherent with current knowledge about viruses but that viral agents play key roles within the roots and stem of the tree of life.

Phylogeny of Banana Streak Virus reveals recent and repetitive endogenization in the genome of its banana host (Musa 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.

Retroviral promoters in the human genome

MOTIVATION

Endogenous retrovirus (ERV) elements have been shown to contribute promoter sequences that can initiate transcription of adjacent human genes. However, the extent to which retroviral sequences initiate transcription within the human genome is currently unknown. We analyzed genome sequence and high-throughput expression data to systematically evaluate the presence of retroviral promoters in the human genome.

RESULTS

We report the existence of 51,197 ERV-derived promoter sequences that initiate transcription within the human genome, including 1743 cases where transcription is initiated from ERV sequences that are located in gene proximal promoter or 5' untranslated regions (UTRs). A total of 114 of the ERV-derived transcription start sites can be demonstrated to drive transcription of 97 human genes, producing chimeric transcripts that are initiated within ERV long terminal repeat (LTR) sequences and read-through into known gene sequences. ERV promoters drive tissue-specific and lineage-specific patterns of gene expression and contribute to expression divergence between paralogs. These data illustrate the potential of retroviral sequences to regulate human transcription on a large scale consistent with a substantial effect of ERVs on the function and evolution of the human genome.

The heterochromatic copies of the LTR retrotransposons as a record of the genomic events that have shaped the Drosophila melanogaster genome

Transposable elements, which are major components of most genomes, are known to accumulate in heterochromatic regions in which they have progressively diverged in sequence by mutations and internal deletions and insertions (indels) during the course of evolution. They therefore provide a record of the genomic events that have shaped the genomes, some of which could correspond to speciation events. Using the sequence divergence between the long terminal repeats (LTRs), we estimated the date of the insertion events of the LTR retrotransposon copies embedded within the heterochromatin regions of the Drosophila melanogaster genome. We did not detect traces of any specific waves of mobilization of retrotransposons within heterochromatin, apart from a very recent wave, which corresponds to the numerous LTR retrotransposon copies found in euchromatin.

Retroviral invasion of the koala genome

Endogenous retroviruses are a common ancestral feature of mammalian genomes with most having been inactivated over time through mutation and deletion. A group of more intact endogenous retroviruses are considered to have entered the genomes of some species more recently, through infection by exogenous viruses, but this event has never been directly proved. We have previously reported koala retrovirus (KoRV) to be a functional virus that is associated with neoplasia. Here we show that KoRV also shows features of a recently inserted endogenous retrovirus that is vertically transmitted. The finding that some isolated koala populations have not yet incorporated KoRV into their genomes, combined with its high level of activity and variability in individual koalas, suggests that KoRV is a virus in transition between an exogenous and endogenous element. This ongoing dynamic interaction with a wild species provides an exciting opportunity to study the process and consequences of retroviral endogenization in action, and is an attractive model for studying the evolutionary event in which a retrovirus invades a mammalian genome.

In vitro characterization of a koala retrovirus

Recently, a new endogenous koala gammaretrovirus, designated KoRV, was isolated from koalas. The KoRV genome shares 78% nucleotide identity with another gammaretrovirus, gibbon ape leukemia virus (GALV). KoRV is endogenous in koalas, while GALV is exogenous, suggesting that KoRV predates GALV and that gibbons and koalas acquired the virus at different times from a common source. We have determined that subtle adaptive differences between the KoRV and GALV envelope genes account for differences in their receptor utilization properties. KoRV represents a unique example of a gammaretrovirus whose envelope has evolved to allow for its expanded host range and zoonotic potential.

Impact of transposable elements on the evolution of mammalian gene regulation

Transposable elements (TEs) are present in all organisms and nearly half of the human and mouse genome is derived from ancient transpositions. This fact alone suggests that TEs have played a major role in genome organization and evolution. Studies undertaken over the last two decades or so clearly show that TEs of various kinds have played an important role in organism evolution. Here we review the impact TEs have on the evolution of gene regulation and gene function with an emphasis on humans. Understanding the mechanisms resulting in genomic change is central to our understanding of gene regulation, genetic disease and genome evolution. Full comprehension of these biological processes is not possible without an in depth knowledge of how TEs impact upon the genome.

Comprehensive analysis of human endogenous retrovirus transcriptional activity in human tissues with a retrovirus-specific microarray

Retrovirus-like sequences account for 8 to 9% of the human genome. Among these sequences, about 8,000 pol-containing proviral elements have been identified to date. As part of our ongoing search for active and possibly disease-relevant human endogenous retroviruses (HERVs), we have recently developed an oligonucleotide-based microarray. The assay allows for both the detection and the identification of most known retroviral reverse transcriptase (RT)-related nucleic acids in biological samples. In the present study, we have investigated the transcriptional activity of representative members of 20 HERV families in 19 different normal human tissues. Qualitative evaluation of chip hybridization signals and quantitative analysis by real-time RT-PCR revealed distinct HERV activity in the human tissues under investigation, suggesting that HERV elements are active in human cells in a tissue-specific manner. Most active members of HERV families were found in mRNA prepared from skin, thyroid gland, placenta, and tissues of reproductive organs. In contrast, only few active HERVs were detectable in muscle cells. Human tissues that lack HERV transcription could not be found, confirming that human endogenous retroviruses are permanent components of the human transcriptome. Distinct activity patterns may reflect the characteristics of the regulatory machinery in these cells, e.g., cell type-dependent occurrence of transcriptional regulatory factors.

The modern molecular clock

The discovery of the molecular clock--a relatively constant rate of molecular evolution--provided an insight into the mechanisms of molecular evolution, and created one of the most useful new tools in biology. The unexpected constancy of rate was explained by assuming that most changes to genes are effectively neutral. Theory predicts several sources of variation in the rate of molecular evolution. However, even an approximate clock allows time estimates of events in evolutionary history, which provides a method for testing a wide range of biological hypotheses ranging from the origins of the animal kingdom to the emergence of new viral epidemics.