Contrasting segregation patterns among endogenous retroviruses across the koala population.

Koalas (Phascolarctos cinereus) have experienced a history of retroviral epidemics leaving their trace as heritable endogenous retroviruses (ERVs) in their genomes. A recently identified ERV lineage, named phaCin-ß, shows a pattern of recent, possibly current, activity with high insertional polymorphism in the population. Here, we investigate geographic patterns of three focal ERV lineages of increasing estimated ages, from the koala retrovirus (KoRV) to phaCin-ß and to phaCin-ß-like, using the whole-genome sequencing of 430 koalas from the Koala Genome Survey. Thousands of ERV loci were found across the population, with contrasting patterns of polymorphism. Northern individuals had thousands of KoRV integrations and hundreds of phaCin-ß ERVs. In contrast, southern individuals had higher phaCin-ß frequencies, possibly reflecting more recent activity and a founder effect. Overall, our findings suggest high ERV burden in koalas, reflecting historic retrovirus-host interactions. Importantly, the ERV catalogue supplies improved markers for conservation genetics in this endangered species.

Genome Reference Assembly for Bottlenecked Southern Australian Koalas.

Koala populations show marked differences in inbreeding levels and in the presence or absence of the endogenous Koala retrovirus (KoRV). These genetic differences among populations may lead to severe disease impacts threatening koala population viability. In addition, the recent colonization of the koala genome by KoRV provides a unique opportunity to study the process of retroviral adaptation to vertebrate genomes and the impact this has on speciation, genome structure, and function. The genome build described here is from an animal from the bottlenecked Southern population free of endogenous and exogenous KoRV. It provides a more contiguous genome build than the previous koala reference derived from an animal from a more outbred Northern population and is the first koala genome from a KoRV polymerase-free animal.

Expansion of a retrovirus lineage in the koala genome.

Retroviruses have left their legacy in host genomes over millions of years as endogenous retroviruses (ERVs), and their structure, diversity, and prevalence provide insights into the historical dynamics of retrovirus-host interactions. In bioinformatic analyses of koala (Phascolarctos cinereus) whole-genome sequences, we identify a recently expanded ERV lineage (phaCin-ß) that is related to the New World squirrel monkey retrovirus. This ERV expansion shares many parallels with the ongoing koala retrovirus (KoRV) invasion of the koala genome, including highly similar and mostly intact sequences, and polymorphic ERV loci in the sampled koala population. The recent phaCin-ß ERV colonization of the koala genome appears to predate the current KoRV invasion, but polymorphic ERVs and divergence comparisons between these two lineages predict a currently uncharacterized, possibly still extant, phaCin-ß retrovirus. The genomics approach to ERV-guided discovery of novel retroviruses in host species provides a strong incentive to search for phaCin-ß retroviruses in the Australasian fauna.

A targeted approach to investigating immune genes of an iconic Australian marsupial.

Disease is a contributing factor to the decline of wildlife populations across the globe. Koalas, iconic yet declining Australian marsupials, are predominantly impacted by two pathogens, Chlamydia and koala retrovirus. Chlamydia is an obligate intracellular bacterium and one of the most widespread sexually transmitted infections in humans worldwide. In koalas, Chlamydia infections can present as asymptomatic or can cause a range of ocular and urogenital disease signs, such as conjunctivitis, cystitis and infertility. In this study, we looked at differences in response to Chlamydia in two northern populations of koalas using a targeted gene sequencing of 1209 immune genes in addition to genome-wide reduced representation data. We identified two MHC Class I genes associated with Chlamydia disease progression as well as 25 single nucleotide polymorphisms across 17 genes that were associated with resolution of Chlamydia infection. These genes are involved in the innate immune response (TLR5) and defence (TLR5, IFNγ, SERPINE1, STAT2 and STX4). This study deepens our understanding of the role that genetics plays in disease progression in koalas and leads into future work that will use whole genome resequencing of a larger sample set to investigate in greater detail regions identified in this study. Elucidation of the role of host genetics in disease progression and resolution in koalas will directly contribute to better design of Chlamydia vaccines and management of koala populations which have recently been listed as "endangered."

Subtype distribution and expression of the koala retrovirus in the Japanese zoo koala population.

We investigated the proviral copies and RNA expression in koala retrovirus (KoRV)-infected koalas. To ascertain any variation in viral load by institution, age, sex, or body condition score, we quantified KoRV proviral DNA and RNA loads in captive koalas (n = 37) reared in Japanese zoos. All koalas were positive for KoRV genes (pol, LTRs, and env of KoRV-A) in genomic DNA (gDNA), and 91.89% were positive for the pol gene in RNA. In contrast, the distribution rates of KoRV-B, KoRV-C, KoRV-D, and KoRV-F env genes in gDNA were 94.59%, 27.03%, 67.57%, and 54.05%, respectively. A wide inter-individual variation and/or a significant inter-institutional difference in proviral DNA (p < 0.0001) and RNA (p < 0.001) amounts (copies/103 koala ß-actin copies) were observed in Awaji Farm England Hill Zoo koalas, which were obtained from southern koala populations, suggesting exogenous incorporation of KoRV in these koalas. Significant (p < 0.05) age differences were noted in KoRV RNA load (p < 0.05) and median total RNA load (p < 0.001), with loads higher in younger koalas (joeys and juveniles). Thus, the current study provides the distribution of KoRV subtypes in Japanese zoo koala populations and identifies several additional risk factors (sex, age, and body condition) associated with KoRV expression.

Retroviral integrations contribute to elevated host cancer rates during germline invasion.

Repeated retroviral infections of vertebrate germlines have made endogenous retroviruses ubiquitous features of mammalian genomes. However, millions of years of evolution obscure many of the immediate repercussions of retroviral endogenisation on host health. Here we examine retroviral endogenisation during its earliest stages in the koala (Phascolarctos cinereus), a species undergoing germline invasion by koala retrovirus (KoRV) and affected by high cancer prevalence. We characterise KoRV integration sites (IS) in tumour and healthy tissues from 10 koalas, detecting 1002 unique IS, with hotspots of integration occurring in the vicinity of known cancer genes. We find that tumours accumulate novel IS, with proximate genes over-represented for cancer associations. We detect dysregulation of genes containing IS and identify a highly-expressed transduced oncogene. Our data provide insights into the tremendous mutational load suffered by the host during active retroviral germline invasion, a process repeatedly experienced and overcome during the evolution of vertebrate lineages.

High-throughput immunogenetic typing of koalas suggests possible link between MHC alleles and cancers.

Characterizing the allelic diversity within major histocompatibility complex (MHC) genes is an important way of determining the potential genetic resilience of a population to infectious and ecological pressures. For the koala (Phascolarctos cinereus), endemic diseases, anthropogenic factors and climate change are all placing increased pressure on this vulnerable marsupial. To increase the ability of researchers to study MHC genetics in koalas, this study developed and tested a high-throughput immunogenetic profiling methodology for targeting MHC class I UA and UC genes and MHC class II DAB, DBB, DCB and DMB genes in a population of 82 captive koalas. This approach was validated by comparing the determined allelic profiles from 36 koala family units (18 dam-sire-joey units and 18 parent-joey pairs), finding 96% overall congruence within family profiles. Cancers are a significant cause of morbidity in koalas and the risk factors remain undetermined. Our analysis of this captive population revealed several novel MHC alleles, including a potential link between the DBB*03 allele and a risk of developing cancer. This method offers a reliable, high-throughput protocol for expanded study into koala immunogenetics.

Koala immunogenetics and chlamydial strain type are more directly involved in chlamydial disease progression in koalas from two south east Queensland koala populations than koala retrovirus subtypes.

Chlamydial disease control is increasingly utilised as a management tool to stabilise declining koala populations, and yet we have a limited understanding of the factors that contribute to disease progression. To examine the impact of host and pathogen genetics, we selected two geographically separated south east Queensland koala populations, differentially affected by chlamydial disease, and analysed koala major histocompatibility complex (MHC) genes, circulating strains of Chlamydia pecorum and koala retrovirus (KoRV) subtypes in longitudinally sampled, well-defined clinical groups. We found that koala immunogenetics and chlamydial genotypes differed between the populations. Disease progression was associated with specific MHC alleles, and we identified two putative susceptibility (DCb 03, DBb 04) and protective (DAb 10, UC 01:01) variants. Chlamydial genotypes belonging to both Multi-Locus Sequence Typing sequence type (ST) 69 and ompA genotype F were associated with disease progression, whereas ST 281 was associated with the absence of disease. We also detected different ompA genotypes, but not different STs, when long-term infections were monitored over time. By comparison, KoRV profiles were not significantly associated with disease progression. These findings suggest that chlamydial genotypes vary in pathogenicity and that koala immunogenetics and chlamydial strains are more directly involved in disease progression than KoRV subtypes.

The piRNA Response to Retroviral Invasion of the Koala Genome.

Antisense Piwi-interacting RNAs (piRNAs) guide silencing of established transposons during germline development, and sense piRNAs drive ping-pong amplification of the antisense pool, but how the germline responds to genome invasion is not understood. The KoRV-A gammaretrovirus infects the soma and germline and is sweeping through wild koalas by a combination of horizontal and vertical transfer, allowing direct analysis of retroviral invasion of the germline genome. Gammaretroviruses produce spliced Env mRNAs and unspliced transcripts encoding Gag, Pol, and the viral genome, but KoRV-A piRNAs are almost exclusively derived from unspliced genomic transcripts and are strongly sense-strand biased. Significantly, selective piRNA processing of unspliced proviral transcripts is conserved from insects to placental mammals. We speculate that bypassed splicing generates a conserved molecular pattern that directs proviral genomic transcripts to the piRNA biogenesis machinery and that this "innate" piRNA response suppresses transposition until antisense piRNAs are produced, establishing sequence-specific adaptive immunity.

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.

Long-read genome sequence assembly provides insight into ongoing retroviral invasion of the koala germline.

The koala retrovirus (KoRV) is implicated in several diseases affecting the koala (Phascolarctos cinereus). KoRV provirus can be present in the genome of koalas as an endogenous retrovirus (present in all cells via germline integration) or as exogenous retrovirus responsible for somatic integrations of proviral KoRV (present in a limited number of cells). This ongoing invasion of the koala germline by KoRV provides a powerful opportunity to assess the viral strategies used by KoRV in an individual. Analysis of a high-quality genome sequence of a single koala revealed 133 KoRV integration sites. Most integrations contain full-length, endogenous provirus; KoRV-A subtype. The second most frequent integrations contain an endogenous recombinant element (recKoRV) in which most of the KoRV protein-coding region has been replaced with an ancient, endogenous retroelement. A third set of integrations, with very low sequence coverage, may represent somatic cell integrations of KoRV-A, KoRV-B and two recently designated additional subgroups, KoRV-D and KoRV-E. KoRV-D and KoRV-E are missing several genes required for viral processing, suggesting they have been transmitted as defective viruses. Our results represent the first comprehensive analyses of KoRV integration and variation in a single animal and provide further insights into the process of retroviral-host species interactions.

Infection with koala retrovirus subgroup B (KoRV-B), but not KoRV-A, is associated with chlamydial disease in free-ranging koalas (Phascolarctos cinereus).

The virulence of chlamydial infection in wild koalas is highly variable between individuals. Some koalas can be infected (PCR positive) with Chlamydia for long periods but remain asymptomatic, whereas others develop clinical disease. Chlamydia in the koala has traditionally been studied without regard to coinfection with other pathogens, although koalas are usually subject to infection with koala retrovirus (KoRV). Retroviruses can be immunosuppressive, and there is evidence of an immunosuppressive effect of KoRV in vitro. Originally thought to be a single endogenous strain, a new, potentially more virulent exogenous variant (KoRV-B) was recently reported. We hypothesized that KoRV-B might significantly alter chlamydial disease outcomes in koalas, presumably via immunosuppression. By studying sub-groups of Chlamydia and KoRV infected koalas in the wild, we found that neither total KoRV load (either viraemia or proviral copies per genome), nor chlamydial infection level or strain type, was significantly associated with chlamydial disease risk. However, PCR positivity with KoRV-B was significantly associated with chlamydial disease in koalas (p = 0.02961). This represents an example of a recently evolved virus variant that may be predisposing its host (the koala) to overt clinical disease when co-infected with an otherwise asymptomatic bacterial pathogen (Chlamydia).

Characterization of the antimicrobial peptide family defensins in the Tasmanian devil (Sarcophilus harrisii), koala (Phascolarctos cinereus), and tammar wallaby (Macropus eugenii).

Defensins comprise a family of cysteine-rich antimicrobial peptides with important roles in innate and adaptive immune defense in vertebrates. We characterized alpha and beta defensin genes in three Australian marsupials: the Tasmanian devil (Sarcophilus harrisii), koala (Phascolarctos cinereus), and tammar wallaby (Macropus eugenii) and identified 48, 34, and 39 defensins, respectively. One hundred and twelve have the classical antimicrobial peptides characteristics required for pathogen membrane targeting, including cationic charge (between 1+ and 15+) and a high proportion of hydrophobic residues (>30%). Phylogenetic analysis shows that gene duplication has driven unique and species-specific expansions of devil, koala, and tammar wallaby beta defensins and devil alpha defensins. Defensin genes are arranged in three genomic clusters in marsupials, whereas further duplications and translocations have occurred in eutherians resulting in four and five gene clusters in mice and humans, respectively. Marsupial defensins are generally under purifying selection, particularly residues essential for defensin structural stability. Certain hydrophobic or positively charged sites, predominantly found in the defensin loop, are positively selected, which may have functional significance in defensin-target interaction and membrane insertion.

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.

One hundred twenty years of koala retrovirus evolution determined from museum skins.

Although endogenous retroviruses are common across vertebrate genomes, the koala retrovirus (KoRV) is the only retrovirus known to be currently invading the germ line of its host. KoRV is believed to have first infected koalas in northern Australia less than two centuries ago. We examined KoRV in 28 koala museum skins collected in the late 19th and 20th centuries and deep sequenced the complete proviral envelope region from five northern Australian specimens. Strikingly, KoRV env sequences were conserved among koalas collected over the span of a century, and two functional motifs that affect viral infectivity were fixed across the museum koala specimens. We detected only 20 env polymorphisms among the koalas, likely representing derived mutations subject to purifying selection. Among northern Australian koalas, KoRV was already ubiquitous by the late 19th century, suggesting that KoRV evolved and spread among koala populations more slowly than previously believed. Given that museum and modern koalas share nearly identical KoRV sequences, it is likely that koala populations, for more than a century, have experienced increased susceptibility to diseases caused by viral pathogenesis.

Biology and evolution of the endogenous koala retrovirus.

Although endogenous retroviruses are ubiquitous features of all mammalian genomes, the process of initial germ line invasion and subsequent inactivation from a pathogenic element has not yet been observed in a wild species. Koala retrovirus (KoRV) provides a unique opportunity to study this process of endogenisation in action as it still appears to be spreading through the koala population. Ongoing expression of the endogenous sequence and consequent high levels of viraemia have been linked to neoplasia and immunosuppression in koalas. This apparently recent invader of the koala genome shares a remarkably close sequence relationship with the pathogenic exogenous Gibbon ape leukaemia virus (GALV), and comparative analyses of KoRV and GALVare helping to shed light on how retroviruses in general adapt to a relatively benign or at least less pathogenic existence within a new host genome. (Part of a multi-author review).

Koala retrovirus: a genome invasion in real time.

Koalas are currently undergoing a wave of germline infections by the retrovirus KoRV. Study of this phenomenon not only provides an opportunity for understanding the processes regulating retrovirus endogenization but may also be essential to preventing the extinction of the species.

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.