Endogenous type D retrovirus in a marsupial, the common brushtail possum (Trichosurus vulpecula)

Long-term host-pathogen evolution of endogenous beta- and gammaretroviruses in mouse lemurs with little evidence of recent retroviral introgression

Madagascar's flora and fauna have evolved in relative isolation since the island split from the African and Indian continents. When the last common ancestors of lemurs left Africa between 40 and 70 million years ago, they carried a subset of the viral diversity of the mainland population within them, which continued to evolve throughout the lemur radiation. Relative to other primate radiations, we know very little about the past or present viral diversity of lemurs, particularly mouse lemurs. Using high-throughput sequencing, we identified two gammaretroviruses and three betaretroviruses in the genomes of four species of wild mouse lemurs. The two gammaretroviruses and two betaretroviruses have not previously been described. One betaretrovirus was previously identified. All identified viruses are present in both Lorisiformes and Lemuriformes but absent from haplorrhine primates. The estimated ages of these viruses are consistent with the estimated divergence dates of the host lineages, suggesting they colonized the lemur genome after the Haplorrhine-Strepsirrhine split, but before the Lorisiformes-Lemuriformes split and before the colonization of Madagascar. The viral phylogenies connect multiple lineages of retroviruses from non-lemur and non-Madagascar-native species, suggesting substantial cross-species transmission occurred deep in the primate clade prior to its geographic dispersal. These phylogenies provide novel insights into known retroviral clades. They suggest that the origin of gammaretroviruses in rodents or bats may be premature and that the Jaagsiekte sheep virus clade may be older and more broadly distributed among mammals than previously thought.

Identification of a novel polyomavirus from a marsupial host

We report the identification and analysis of a full sequence of a novel polyomavirus from a brushtail possum (Trichosurus vulpecula ) termed possum polyomavirus (PPyV). The sequence was obtained from the next-generation sequencing assembly during an investigation into the aetiological agent for a neurological disease of possums termed wobbly possum disease (WPD), but the virus was not aetiologically involved in WPD. The PPyV genome was 5,224 nt long with the organisation typical for polyomaviruses, including early (large and small T antigens) and late (Viral Protein 1 (VP1), VP2, and VP3) coding regions separated by the non-coding control region of 465 nt. PPyV clustered with betapolyomaviruses in the WUKI clade but showed less than 60 per cent identity to any of the members of this clade. We propose that PPyV is classified within a new species in the genus Betapolyomavirus . These data add to our limited knowledge of marsupial viruses and their evolution.

Detection and Characterisation of an Endogenous Betaretrovirus in Australian Wild Deer

Endogenous retroviruses (ERVs) are the remnants of past retroviral infections that once invaded the host’s germline and were vertically transmitted. ERV sequences have been reported in mammals, but their distribution and diversity in cervids are unclear. Using next-generation sequencing, we identified a nearly complete genome of an endogenous betaretrovirus in fallow deer (Dama dama). Further genomic analysis showed that this provirus, tentatively named cervid endogenous betaretrovirus 1 (CERV β1), has typical betaretroviral genome features (gag-pro-pol-env) and the betaretrovirus-specific dUTPase domain. In addition, CERV β1 pol sequences were detected by PCR in the six non-native deer species with wild populations in Australia. Phylogenetic analyses demonstrated that CERV β1 sequences from subfamily Cervinae clustered as sister taxa to ERV-like sequences in species of subfamily Muntiacinae. These findings, therefore, suggest that CERV β1 endogenisation occurred after the split of these two subfamilies (between 3.3 and 5 million years ago). Our results provide important insights into the evolution of betaretroviruses in cervids.

An endogenous retrovirus presumed to have been endogenized or relocated recently in a marsupial, the red-necked wallaby

An albino infant wallaby was born to a mother with the wild-type body color. PCR and sequencing analyses of TYR (encoding tyrosinase, which is essential for melanin biosynthesis) of this albino wallaby revealed a 7.1-kb-long DNA fragment inserted in the first exon. Because the fragment carried long terminal repeats, we assumed it to be a copy of an endogenous retrovirus, which we named walb. We cloned other walb copies residing in the genomes of this species and another wallaby species. The copies exhibited length variation, and the longest copy (>8.0 kb) contained open reading frames whose deduced amino acid sequences were well aligned with those of gag, pol, and env of retroviruses. It is not known through which of the following likely processes the walb copy was inserted into TYR: endogenization (infection of a germline cell by an exogenous virus), reinfection (infection by a virus produced from a previously endogenized provirus), or retrotransposition (intracellular relocation of a provirus). In any case, the insertion into TYR is considered to have been a recent event on an evolutionary timescale because albino mutant alleles generally do not persist for long because of their deleterious effects in wild circumstances. .

Transcriptomic and genomic variants between koala populations reveals underlying genetic components to disorders in a bottlenecked population

Historical hunting pressures on koalas in the southern part of their range in Australia have led to a marked genetic bottleneck when compared with their northern counterparts. There are a range of suspected genetic disorders such as testicular abnormalities, oxalate nephrosis and microcephaly reported at higher prevalence in these genetically restricted southern animals. This paper reports analysis of differential expression of genes from RNAseq of lymph nodes, SNPs present in genes and the fixation index (population differentiation due to genetic structure) of these SNPs from two populations, one in south east Queensland, representative of the northern genotype and one in the Mount Lofty Ranges South Australia, representative of the southern genotype. SNPs that differ between these two populations were significantly enriched in genes associated with brain diseases. Genes which were differentially expressed between the two populations included many associated with brain development or disease, and in addition a number associated with testicular development, including the androgen receptor. Finally, one of the 8 genes both differentially expressed and with a statistical difference in SNP frequency between populations was SLC26A6 (solute carrier family 26 member 6), an anion transporter that was upregulated in SA koalas and is associated with oxalate transport and calcium oxalate uroliths in humans. Together the differences in SNPs and gene expression described in this paper suggest an underlying genetic basis for several disorders commonly seen in southern Australian koalas, supporting the need for further research into the genetic basis of these conditions, and highlighting that genetic selection in managed populations may need to be considered in the future.

Retroviruses of the RDR superinfection interference group: ancient origins and broad host distribution of a promiscuous Env gene

Due to recombination, different regions of a retrovirus genome can have distinct phylogenetic histories. The RD114-and-D-type-retrovirus (RDR) interference group provides an extreme example: the RDR group comprises a variety of taxonomically distinct retroviruses, isolated from diverse mammalian and avian hosts, that share a homologous env gene and use the same cell-surface entry receptor. RDR env homologs are also found among ancient endogenous retrovirus (ERV) sequences, including the syncytin genes of humans and rabbits, indicating that RDR Env glycoproteins have likely mediated endogenization on multiple occasions in diverse vertebrate lineages. The distribution of RDR env among exogenous and endogenous retroviruses indicates that it has been swapped between viruses many times, and that it likely facilitated multiple cross-species transmission events spanning millions of years of vertebrate evolution.

Genome-Wide Screening of Retroviral Envelope Genes in the Nine-Banded Armadillo (Dasypus novemcinctus, Xenarthra) Reveals an Unfixed Chimeric Endogenous Betaretrovirus Using the ASCT2 Receptor

Retroviruses enter host cells through the interaction of their envelope (Env) protein with a cell surface receptor, which triggers the fusion of viral and cellular membranes. The sodium-dependent neutral amino acid transporter ASCT2 is the common receptor of the large RD114 retrovirus interference group, whose members display frequent env recombination events. Germ line retrovirus infections have led to numerous inherited endogenous retroviruses (ERVs) in vertebrate genomes, which provide useful insights into the coevolutionary history of retroviruses and their hosts. Rare ERV-derived genes display conserved viral functions, as illustrated by the fusogenic syncytin env genes involved in placentation. Here, we searched for functional env genes in the nine-banded armadillo (Dasypus novemcinctus) genome and identified dasy-env1.1, which clusters with RD114 interference group env genes and with two syncytin genes sharing ASCT2 receptor usage. Using ex vivo pseudotyping and cell-cell fusion assays, we demonstrated that the Dasy-Env1.1 protein is fusogenic and can use both human and armadillo ASCT2s as receptors. This gammaretroviral env gene belongs to a provirus with betaretrovirus-like features, suggesting acquisition through recombination. Provirus insertion was found in several Dasypus species, where it has not reached fixation, whereas related family members integrated before diversification of the genus Dasypus >12 million years ago (Mya). This newly described ERV lineage is potentially useful as a population genetic marker. Our results extend the usage of ASCT2 as a retrovirus receptor to the mammalian clade Xenarthra and suggest that the acquisition of an ASCT2-interacting env gene is a major selective force driving the emergence of numerous chimeric viruses in vertebrates.

IMPORTANCE

Retroviral infection is initiated by the binding of the viral envelope glycoprotein to a host cell receptor(s), triggering membrane fusion. Ancient germ line infections have generated numerous endogenous retroviruses (ERVs) in nearly all vertebrate genomes. Here, we report a previously uncharacterized ERV lineage from the genome of a xenarthran species, the nine-banded armadillo (Dasypus novemcinctus). It entered the Dasypus genus >12 Mya, with one element being inserted more recently in some Dasypus species, where it could serve as a useful marker for population genetics. This element exhibits an env gene, acquired by recombination events, with conserved viral fusogenic properties through binding to ASCT2, a receptor used by a wide range of recombinant retroviruses infecting other vertebrate orders. This specifies the ASCT2 transporter as a successful receptor for ERV endogenization and suggests that ASCT2-binding env acquisition events have favored the emergence of numerous chimeric viruses in a wide range of species.

Characterization of a novel full-length bovine endogenous retrovirus, BERV-β1

Recent studies have suggested that certain classes of endogenous retroviruses (ERVs) may be present in cattle. The aim of this study was increase the scope of knowledge regarding bovine ERVs. The ovine ERV β1 pro/pol sequence was used to design a primer set for polymerase chain reaction (PCR) amplification of a similar sequence in the bovine genome. Through phylogenetic and bioinformatic analysis of the PCR product sequence together with its flanking region, a sequence 8107 bp in length was characterized. This sequence had a typical 5'-LTR-gag-pro-pol-env-LTR-3' organization, and phylogenetic investigation defined it as a bovine ERV β1. Thus, we were able to identify a novel bovine endogenous retrovirus element.

Suppression subtractive hybridization method for the identification of a new strain of murine hepatitis virus from xenografted SCID mice

During attempts to clone retroviral determinants associated with a mouse model of Langerhans cell histiocytosis (LCH), suppression subtractive hybridization (SSH) was used to identify unique viruses in the liver of severe combined immunodeficiency (SCID) mice transplanted with LCH tissues. A partial genomic sequence of a murine coronavirus was identified, and the whole genome (31428 bp) of the coronavirus was subsequently sequenced using PCR cloning techniques. Nucleotide sequence comparisons revealed that the genome sequence of the new virus was 91-93 % identical to those of known murine hepatitis viruses (MHVs). The predicted open reading frame from the nucleotide sequence encoded all known proteins of MHVs. Analysis at the protein level showed that the virus was closely related to the highly virulent MHV-JHM strain. The virus strain was named MHV-MI. No type D retroviruses were found. Degenerate PCR targeting of type D retrovirus and 5′-RACE targeting of other types of retroviruses confirmed the absence of any retroviral association with the LCH xenografted SCID mice.

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.

A novel endogenous betaretrovirus in the common vampire bat (Desmodus rotundus) suggests multiple independent infection and cross-species transmission events

The Desmodus rotundus endogenous betaretrovirus (DrERV) is fixed in the vampire bat D. rotundus population and in other phyllostomid bats but is not present in all species from this family. DrERV is not phylogenetically related to Old World bat betaretroviruses but to betaretroviruses from rodents and New World primates, suggesting recent cross-species transmission. A recent integration age estimation of the provirus in some taxa indicates that an exogenous counterpart might have been in recent circulation.

A transcriptome resource for the koala (Phascolarctos cinereus): insights into koala retrovirus transcription and sequence diversity

Background

The koala, Phascolarctos cinereus, is a biologically unique and evolutionarily distinct Australian arboreal marsupial. The goal of this study was to sequence the transcriptome from several tissues of two geographically separate koalas, and to create the first comprehensive catalog of annotated transcripts for this species, enabling detailed analysis of the unique attributes of this threatened native marsupial, including infection by the koala retrovirus.

Results

RNA-Seq data was generated from a range of tissues from one male and one female koala and assembled de novo into transcripts using Velvet-Oases. Transcript abundance in each tissue was estimated. Transcripts were searched for likely protein-coding regions and a non-redundant set of 117,563 putative protein sequences was produced. In similarity searches there were 84,907 (72%) sequences that aligned to at least one sequence in the NCBI nr protein database. The best alignments were to sequences from other marsupials. After applying a reciprocal best hit requirement of koala sequences to those from tammar wallaby, Tasmanian devil and the gray short-tailed opossum, we estimate that our transcriptome dataset represents approximately 15,000 koala genes. The marsupial alignment information was used to look for potential gene duplications and we report evidence for copy number expansion of the alpha amylase gene, and of an aldehyde reductase gene.

Koala retrovirus (KoRV) transcripts were detected in the transcriptomes. These were analysed in detail and the structure of the spliced envelope gene transcript was determined. There was appreciable sequence diversity within KoRV, with 233 sites in the KoRV genome showing small insertions/deletions or single nucleotide polymorphisms. Both koalas had sequences from the KoRV-A subtype, but the male koala transcriptome has, in addition, sequences more closely related to the KoRV-B subtype. This is the first report of a KoRV-B-like sequence in a wild population.

Conclusions

This transcriptomic dataset is a useful resource for molecular genetic studies of the koala, for evolutionary genetic studies of marsupials, for validation and annotation of the koala genome sequence, and for investigation of koala retrovirus. Annotated transcripts can be browsed and queried at http://koalagenome.org.

Electronic supplementary material

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

Identification of diverse full-length endogenous betaretroviruses in megabats and microbats

BACKGROUND

Betaretroviruses infect a wide range of species including primates, rodents, ruminants, and marsupials. They exist in both endogenous and exogenous forms and are implicated in animal diseases such as lung cancer in sheep, and in human disease, with members of the human endogenous retrovirus-K (HERV-K) group of endogenous betaretroviruses (βERVs) associated with human cancers and autoimmune diseases. To improve our understanding of betaretroviruses in an evolutionarily distinct host species, we characterized βERVs present in the genomes and transcriptomes of mega- and microbats, which are an important reservoir of emerging viruses.

RESULTS

A diverse range of full-length βERVs were discovered in mega- and microbat genomes and transcriptomes including the first identified intact endogenous retrovirus in a bat. Our analysis revealed that the genus Betaretrovirus can be divided into eight distinct sub-groups with evidence of cross-species transmission. Betaretroviruses are revealed to be a complex retrovirus group, within which one sub-group has evolved from complex to simple genomic organization through the acquisition of an env gene from the genus Gammaretrovirus. Molecular dating suggests that bats have contended with betaretroviral infections for over 30 million years.

CONCLUSIONS

Our study reveals that a diverse range of betaretroviruses have circulated in bats for most of their evolutionary history, and cluster with extant betaretroviruses of divergent mammalian lineages suggesting that their distribution may be largely unrestricted by host species barriers. The presence of βERVs with the ability to transcribe active viral elements in a major animal reservoir for viral pathogens has potential implications for public health.

Identification and classification of feline endogenous retroviruses in the cat genome using degenerate PCR and in silico data analysis

The purpose of this study was to identify and classify endogenous retroviruses (ERVs) in the cat genome. Pooled DNA from five domestic cats was subjected to degenerate PCR with primers specific to the conserved retroviral pro/pol region. The 59 amplified retroviral sequences were used for in silico analysis of the cat genome (Felis_catus-6.2). We identified 219 ERV γ and β elements from cat genome contigs, which were classified into 42 ERV γ and 4 β families and further analysed. Among them, 99 γ and 5 β ERV elements contained the complete retroviral structure. Furthermore, we identified 757 spuma-like ERV elements based on the sequence homology to murine (Mu)ERV-L and human (H)ERV-L. To the best of our knowledge, this is the first detailed genome-scale analysis examining Felis catus endogenous retroviruses (FcERV) and providing advanced insights into their structural characteristics, localization in the genome, and diversity.

Characterisation of two enteroviruses isolated from Australian brushtail possums (Trichosurus vulpecula) in New Zealand

Two enteroviruses, designated W1 and W6, were isolated from intestinal contents of Australian brushtail possums (Trichosurus vulpecula) in New Zealand. The genomic sequences of W1 and W6 were 7390 and 7391 nucleotides (nt), respectively. Genetically, possum isolates W1 and W6 were related to bovine enterovirus serotype 2 (BEV-2) strains, especially to the strain PS87/Belfast, based on the capsid protein sequence. However, W1 and W6 formed a clade that was distinct from PS87Belfast based on nucleotide sequences of the 3' and 5'-non-translated region and in the amino acid sequences of 2A, 3C and 3D. Possum isolates W1 and W6 grew more readily in possum kidney cells than in Madin-Darby bovine kidney (MDBK) cells, suggesting that co-evolution of W1 and W6 with possums has made them more adapted to possum cells.

Application of pharmaceutical drug delivery for biological control of the common brushtail possum in New Zealand: a review

The common brushtail possum (Trichosurus vulpecula) is the most significant vertebrate pest in New Zealand, being a major ecological threat to the indigenous biodiversity and an economic threat as a vector for bovine tuberculosis. Novel and effective strategies to reduce the population of T. vulpecula are needed urgently. Several biocontrol agents are currently being assessed and from research to date it is likely that the biocontrol agents will be peptide or protein molecules. It is not possible to administer such biocontrol agents alone because they would be degraded rapidly in the animal, especially if delivered orally. Technologies used in the pharmaceutical industry to design efficacious drug-delivery systems for humans and animals can be applied to the design of delivery systems for biocontrol agents used in wildlife management, although there are some unique challenges that must be overcome.

The pp24 phosphoprotein of Mason-Pfizer monkey virus contributes to viral genome packaging

Background

The Gag protein of Mason-Pfizer monkey virus, a betaretrovirus, contains a phosphoprotein that is cleaved into the Np24 protein and the phosphoprotein pp16/18 during virus maturation. Previous studies by Yasuda and Hunter (J. Virology. 1998. 72:4095–4103) have demonstrated that pp16/18 contains a viral late domain required for budding and that the Np24 protein plays a role during the virus life cycle since deletion of this N-terminal domain blocked virus replication. The function of the Np24 domain, however, is not known.

Results

Here we identify a region of basic residues (KKPKR) within the Np24 domain that is highly conserved among the phosphoproteins of various betaretroviruses. We show that this KKPKR motif is required for virus replication yet dispensable for procapsid assembly, membrane targeting, budding and release, particle maturation, or viral glycoprotein packaging. Additional experiments indicated that deletion of this motif reduced viral RNA packaging 6–8 fold and affected the transient association of Gag with nuclear pores.

Conclusion

These results demonstrate that the Np24 domain plays an important role in RNA packaging and is in agreement with evidence that suggests that correct intracellular targeting of Gag to the nuclear compartment is an fundamental step in the retroviral life cycle.

Evolution and distribution of class II-related endogenous retroviruses

Endogenous retroviruses (ERVs) are widespread in vertebrate genomes and have been loosely grouped into "classes" on the basis of their phylogenetic relatedness to the established genera of exogenous retroviruses. Four of these genera-the lentiviruses, alpharetroviruses, betaretroviruses, and deltaretroviruses-form a well-supported clade in retroviral phylogenies, and ERVs that group with these genera have been termed class II ERVs. We used PCR amplification and sequencing of retroviral fragments from more than 130 vertebrate taxa to investigate the evolution of the class II retroviruses in detail. We confirm that class II retroviruses are largely confined to mammalian and avian hosts and provide evidence for a major novel group of avian retroviruses, and we identify additional members of both the alpha- and the betaretrovirus genera. Phylogenetic analyses demonstrated that the avian and mammalian viruses form distinct monophyletic groups, implying that interclass transmission has occurred only rarely during the evolution of the class II retroviruses. In contrast to previous reports, the lentiviruses clustered as sister taxa to several endogenous retroviruses derived from rodents and insectivores. This topology was further supported by the shared loss of both the class II PR-Pol frameshift site and the class II retrovirus G-patch domain.

Real-time reverse transcriptase PCR for the endogenous koala retrovirus reveals an association between plasma viral load and neoplastic disease in koalas

Koala retrovirus (KoRV) is a newly described endogenous retrovirus and is unusual in that inserts comprise a full-length replication competent genome. As koalas are known to suffer from an extremely high incidence of leukaemia/lymphoma, the association between this retrovirus and disease in koalas was examined. Using quantitative real-time reverse transcriptase PCR it was demonstrated that KoRV RNA levels in plasma are significantly increased in animals suffering from leukaemia or lymphoma when compared with healthy animals. Increased levels of KoRV were also seen for animals with clinical chlamydiosis. A significant positive association between viral RNA levels and age was also demonstrated. Real-time PCR demonstrated as much as 5 log variation in KoRV proviral DNA levels in genomic DNA extracted from whole blood from different animals. Taken together these data indicate that KoRV is an active endogenous retrovirus and suggests that it may be causally linked to neoplastic disease in koalas.

Multiple groups of endogenous betaretroviruses in mice, rats, and other mammals

Betaretroviruses exist in endogenous and exogenous forms in hosts that are widely distributed and evolutionarily distantly related. Here we report the discovery and characterization of several previously unknown betaretrovirus groups in the genomes of Mus musculus and Rattus norvegicus. Each group contains both mouse and rat elements, and several of the groups are more closely related to previously known betaretroviruses from nonmurine hosts. Some of the groups also include members from hosts which were not previously known to harbor betaretroviruses, such as the gray mouse lemur (Microcebus murinus) and Seba's short-tailed bat (Carollia perspicillata). Some of the mouse and rat elements possess intact open reading frames for gag, pro, pol, and/or env genes and display characteristics of having retrotransposed recently. We propose a model whereby betaretroviruses have been evolving within the genomes of murid rodents for at least the last 20 million years and, subsequent to (or concomitant with) the global spread of their murid hosts, have occasionally been transmitted to other species.

The evolution, distribution and diversity of endogenous retroviruses

The retroviral capacity for integration into the host genome can give rise to endogenous retroviruses (ERVs): retroviral sequences that are transmitted vertically as part of the host germ line, within which they may continue to replicate and evolve. ERVs represent both a unique archive of ancient viral sequence information and a dynamic component of host genomes. As such they hold great potential as informative markers for studies of both virus evolution and host genome evolution. Numerous novel ERVs have been described in recent years, particularly as genome sequencing projects have advanced. This review discusses the evolution of ERV lineages, considering the processes by which ERV distribution and diversity is generated. The diversity of ERVs isolated so far is summarised in terms of both their distribution across host taxa, and their relationships to recognised retroviral genera. Finally the relevance of ERVs to studies of genome evolution, host disease and viral ecology is considered, and recent findings discussed.

Molecular confirmation of an adenovirus in brushtail possums (Trichosurus vulpecula)

Partial genome characterisation of a non-cultivable marsupial adenovirus is described. Adenovirus-like particles were found by electron microscopy (EM) in the intestinal contents of brushtail possums (Trichosurus vulpecula) in New Zealand. Using degenerate PCR primers complementary to the most conserved genome regions of adenoviruses, the complete nucleotide sequence of the penton base gene, and partial nucleotide sequences of the DNA polymerase, hexon, and pVII genes were obtained. Phylogenetic analysis of the penton base gene strongly suggested that the brushtail possum adenovirus (candidate PoAdV-1) belongs to the recently proposed genus Atadenovirus. Sequence analysis of the PCR products amplified from the intestinal contents of brushtail possums originating from different geographical regions of New Zealand identified a single genotype. This is the first report of molecular confirmation of an adenovirus in a marsupial.