Retroviruses that colonize the host germline can be passed on as inherited genetic variants. The koala (Phascolarctos cinereus) is currently experiencing germline colonization by two retroviruses, the koala retrovirus (KoRV) and phaCin-β. We analyze the integration site segregation and diversity of endogenous KoRV, phaCin-β, and the related phaCin-β-like in 111 pedigreed koalas from the San Diego Zoo Wildlife Alliance and seven European Zoos. The use of multigenerational pedigrees and the inclusion of health information for each individual koala reveal elimination of retroviruses from proto-oncogenes and the generation and spread of new germline integrations. Seven-hundred-and-fourteen integrations do not persist in the living population. For the 55 triads examined, 21 unique integrations identified in individual koalas are absent in their parents. Retroviral integrations associated with leukemia, fertility, and longevity are used to estimate genetic risk scores and develop a longevity breeding index to minimize neoplasia risk in the captive koala population.

Poncelet’s giant rat (Solomys ponceleti Troughton 1935), is a rare, large murine rodent that is endemic to the Solomon Islands Archipelago in the southwest Pacific Ocean. The species is only known from the adjacent islands of Bougainville and Choiseul and is the largest member of the genus. Here, we describe the complete mitochondrial genome of S. ponceleti and compare it to other Rodentia. The S. ponceleti circular mitogenome was 16,246 bp and contained 13 protein-coding genes, two rRNA genes, 22 tRNAs, and a control region (D-loop). Phylogenetic analysis of selected, published mitogenomes reveals a close relationship to other Australo-Papuan murine rodents.

Genetic management is a critical component of threatened species conservation. Understanding spatial patterns of genetic diversity is essential for evaluating the resilience of fragmented populations to accelerating anthropogenic threats. Nowhere is this more relevant than on the Australian continent, which is experiencing an ongoing loss of biodiversity that exceeds any other developed nation. Using a proprietary genome complexity reduction-based method (DArTSeq), we generated a data set of 3239 high quality Single Nucleotide Polymorphisms (SNPs) to investigate spatial patterns and indices of genetic diversity in the koala (Phascolarctos cinereus), a highly specialised folivorous marsupial that is experiencing rapid and widespread population declines across much of its former range. Our findings demonstrate that current management divisions across the state of New South Wales (NSW) do not fully represent the distribution of genetic diversity among extant koala populations, and that care must be taken to ensure that translocation paradigms based on these frameworks do not inadvertently restrict gene flow between populations and regions that were historically interconnected. We also recommend that koala populations should be prioritised for conservation action based on the scale and severity of the threatening processes that they are currently faced with, rather than placing too much emphasis on their perceived value (e.g., as reservoirs of potentially adaptive alleles), as our data indicate that existing genetic variation in koalas is primarily partitioned among individual animals. As such, the extirpation of koalas from any part of their range represents a potentially critical reduction of genetic diversity for this iconic Australian species.

Climatic and evolutionary processes are inextricably linked to conservation. Avoiding extinction in rapidly changing environments often depends upon a species' capacity to adapt in the face of extreme selective pressures. Here, we employed exon capture and high-throughput next-generation sequencing to investigate the mechanisms underlying population structure and adaptive genetic variation in the koala (Phascolarctos cinereus), an iconic Australian marsupial that represents a unique conservation challenge because it is not uniformly threatened across its range. An examination of 250 specimens representing 91 wild source locations revealed that five major genetic clusters currently exist on a continental scale. The initial divergence of these clusters appears to have been concordant with the Mid-Brunhes Transition (~430 to 300 kya), a major climatic reorganisation that increased the amplitude of Pleistocene glacial-interglacial cycles. While signatures of polygenic selection and environmental adaptation were detected, strong evidence for repeated, climate-associated range contractions and demographic bottleneck events suggests that geographically isolated refugia may have played a more significant role in the survival of the koala through the Pleistocene glaciation than in situ adaptation. Consequently, the conservation of genome-wide genetic variation must be aligned with the protection of core koala habitat to increase the resilience of vulnerable populations to accelerating anthropogenic threats. Finally, we propose that the five major genetic clusters identified in this study should be accounted for in future koala conservation efforts (e.g., guiding translocations), as existing management divisions in the states of Queensland and New South Wales do not reflect historic or contemporary population structure.

Maintenance of genetic diversity within species is a key objective of biodiversity conservation, and small, isolated populations are particularly vulnerable to genetic erosion. Conservation management actions such as predator removal, captive breeding and reintroduction can facilitate numerical recovery of a population, but species often remain at risk from depleted genetic diversity and inbreeding. We investigated dispersal, genetic bottlenecks and genetic population structuring in the island‐dwelling Lord Howe woodhen, a species that came perilously close to extinction in the 1970s. Analyses of mark‐resighting records and variable genetic markers (single‐nucleotide polymorphisms) collected from the contemporary population and 100‐year‐old museum specimens found strong evidence of restricted dispersal at fine spatial scales, with both the contemporary and historic populations exhibiting strong population structuring between mountain and lowland/slopes sites. Additionally, genetic comparison of the contemporary population and historic specimens demonstrated a decline in genetic diversity over the past century. Specifically for the Lord Howe woodhen, we recommend ongoing genetic monitoring and translocations to increase genetic diversity within the re‐established lowland subpopulation. More generally, our results demonstrate how pronounced genetic erosion can arise in species subject to human persecution and predation by introduced predators, and how genetic fragmentation of natural populations can be present at fine geographical scales (less than hundreds of metres). Moreover, without prior information about genetic structure and subsequent genetic monitoring, conservation management can have unexpected negative consequences for the genetic health of populations. Therefore, genetic monitoring and management early in the recovery of populations is desirable to maximize their adaptive potential.

As the global biodiversity crisis deepens, with increasing habitat fragmentation and a changing climate, innovative options for conserving species are being explored. One such conservation action is genetic rescue: introduction of new alleles to promote population fitness. However, for critically endangered species where only one viable population remains, options for introducing new alleles are limited. Interspecies hybridization offers a potential solution but requires resolution of evolutionary relationships, a sound understanding of species biology, social license, and permissive legislative frameworks. Here, we show how phylogenetics and species biology can inform genetic rescue options for the orange‐bellied parrot (OBP; Neophema chrysogaster ), a critically endangered Australian bird with one small remaining wild population. Our phylogenetic analysis of mitochondrial genomes and nuclear loci for all congeneric species provided strong support for OBPs being the sister species to a group comprising elegant, rock, and blue‐winged parrots. Accounting for species distribution, behavior, and ecology, a captive trial of interspecific hybridization with the blue‐winged parrot is recommended, including assessment of the fitness of hybrid individuals. Introduction of new alleles into the OBP genome would achieve the conservation goal of improving genetic diversity in a critically endangered species. Concurrently, legislative issues will need to be resolved.

There are few available tools to comprehensively and economically identify uncharacterized flanking regions that are not extremely labour intensive and which exploit the advantages of emerging long‐read sequencing platforms. We describe SIP; a sonication‐based inverse PCR high‐throughput sequencing strategy to investigate uncharacterized flanking region sequences, including those flanking mobile DNA. SIP combines unbiased fragmentation by sonication and target enrichment by coupling outward facing PCR priming with long‐read sequencing technologies. We demonstrate the effectiveness of SIP by determining retroviral integrations which are high copy and challenging to characterize. We further describe SIP's workflow, examine retroviral (proviral) enrichment and characterize viral structural variants identified. When SIP was coupled with long‐read sequencing using the PacBio RS II platform, proviral integration was extensively characterized at high sequence depth per integration. By interrogating the sequence data, we were also able to test several intrinsic factors including SIP's propensity to form chimeric sequences and adapter ligation efficiencies. SIP is an adaption of a traditional molecular biology technique that can be used to characterize any unknown genomic flanking sequence or to extend any sequence for which only minimal sequence information is available. SIP can be applied broadly to study complex biological systems such as mobile genetic elements with high throughput.

X chromosome inactivation (XCI) mediated by differential DNA methylation between sexes is an iconic example of epigenetic regulation. Although XCI is shared between eutherians and marsupials, the role of DNA methylation in marsupial XCI remains contested. Here, we examine genome-wide signatures of DNA methylation across fives tissues from a male and female koala ( Phascolarctos cinereus ), and present the first whole-genome, multi-tissue marsupial ‘methylome atlas’. Using these novel data, we elucidate divergent versus common features of representative marsupial and eutherian DNA methylation. First, tissue-specific differential DNA methylation in koalas primarily occurs in gene bodies. Second, females show significant global reduction (hypomethylation) of X chromosome DNA methylation compared to males. We show that this pattern is also observed in eutherians. Third, on average, promoter DNA methylation shows little difference between male and female koala X chromosomes, a pattern distinct from that of eutherians. Fourth, the sex-specific DNA methylation landscape upstream of Rsx , the primary lnc RNA associated with marsupial XCI, is consistent with the epigenetic regulation of female-specific (and presumably inactive X chromosome-specific) expression. Finally, we use the prominent female X chromosome hypomethylation and classify 98 previously unplaced scaffolds as X-linked, contributing an additional 14.6 Mb (21.5%) to genomic data annotated as the koala X chromosome. Our work demonstrates evolutionarily divergent pathways leading to functionally conserved patterns of XCI in two deep branches of mammals.

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.

X chromosome inactivation (XCI) mediated by differential DNA methylation between sexes is well characterized in eutherian mammals. Although XCI is shared between eutherians and marsupials, the role of DNA methylation in marsupial XCI remains contested. Here we examine genome-wide signatures of DNA methylation from methylation maps across fives tissues from a male and female koala ( Phascolarctos cinereus ) and present the first whole genome, multi-tissue marsupial “methylome atlas.” Using these novel data, we elucidate divergent versus common features of marsupial and eutherian DNA methylation. First, tissue-specific differential DNA methylation in marsupials primarily occurs in gene bodies. Second, females show significant global reduction (hypomethylation) of X chromosome DNA methylation compared to males. We show that this pattern is also observed in eutherians. Third, on average, promoter DNA methylation shows little difference between male and female koala X chromosomes, a pattern distinct from that of eutherians. Fourth, the sex-specific DNA methylation landscape upstream of Rsx , the primary lnc RNA associated with marsupial XCI, is consistent with the epigenetic regulation of female-(and presumably inactive X chromosome-) specific expression. Finally, we utilize the prominent female X chromosome hypomethylation and classify 98 previously unplaced scaffolds as X-linked, contributing an additional 14.6 Mb (21.5 %) to genomic data annotated as the koala X chromosome. Our work demonstrates evolutionarily divergent pathways leading to functionally conserved patterns of XCI in two deep branches of mammals.### Competing Interest StatementThe authors have declared no competing interest.

OBJECTIVES:Sample preparation for high-throughput sequencing (HTS) includes treatment with various laboratory components, potentially carrying viral nucleic acids, the extent of which has not been thoroughly investigated. Our aim was to systematically examine a diverse repertoire of laboratory components used to prepare samples for HTS in order to identify contaminating viral sequences. METHODS:A total of 322 samples of mainly human origin were analysed using eight protocols, applying a wide variety of laboratory components. Several samples (60% of human specimens) were processed using different protocols. In total, 712 sequencing libraries were investigated for viral sequence contamination. RESULTS:Among sequences showing similarity to viruses, 493 were significantly associated with the use of laboratory components. Each of these viral sequences had sporadic appearance, only being identified in a subset of the samples treated with the linked laboratory component, and some were not identified in the non-template control samples. Remarkably, more than 65% of all viral sequences identified were within viral clusters linked to the use of laboratory components. CONCLUSIONS:We show that high prevalence of contaminating viral sequences can be expected in HTS-based virome data and provide an extensive list of novel contaminating viral sequences that can be used for evaluation of viral findings in future virome and metagenome studies. Moreover, we show that detection can be problematic due to stochastic appearance and limited non-template controls. Although the exact origin of these viral sequences requires further research, our results support laboratory-component-linked viral sequence contamination of both biological and synthetic origin.

BACKGROUND:Viruses and other infectious agents cause more than 15% of human cancer cases. High-throughput sequencing-based studies of virus-cancer associations have mainly focused on cancer transcriptome data. METHODS:In this study, we applied a diverse selection of presequencing enrichment methods targeting all major viral groups, to characterize the viruses present in 197 samples from 18 sample types of cancerous origin. Using high-throughput sequencing, we generated 710 datasets constituting 57 billion sequencing reads. RESULTS:Detailed in silico investigation of the viral content, including exclusion of viral artefacts, from de novo assembled contigs and individual sequencing reads yielded a map of the viruses detected. Our data reveal a virome dominated by papillomaviruses, anelloviruses, herpesviruses, and parvoviruses. More than half of the included samples contained 1 or more viruses; however, no link between specific viruses and cancer types were found. CONCLUSIONS:Our study sheds light on viral presence in cancers and provides highly relevant virome data for future reference.

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.

In recent years, the availability of reduced representation library (RRL) methods has catalysed an expansion of genome-scale studies to characterize both model and non-model organisms. Most of these methods rely on the use of restriction enzymes to obtain DNA sequences at a genome-wide level. These approaches have been widely used to sequence thousands of markers across individuals for many organisms at a reasonable cost, revolutionizing the field of population genomics. However, there are still some limitations associated with these methods, in particular, the high molecular weight DNA required as starting material, the reduced number of common loci among investigated samples, and the short length of the sequenced site-associated DNA. Here, we present MobiSeq, a RRL protocol exploiting simple laboratory techniques, that generates genomic data based on PCR targeted-enrichment of transposable elements and the sequencing of the associated flanking region. We validate its performance across 103 DNA extracts derived from three mammalian species: grey wolf ( Canis lupus ), red deer complex ( Cervus sp .), and brown rat ( Rattus norvegicus ). MobiSeq enables the sequencing of hundreds of thousands loci across the genome, and performs SNP discovery with relatively low rates of clonality. Given the ease and flexibility of MobiSeq protocol, the method has the potential to be implemented for marker discovery and population genomics across a wide range of organisms – enabling the exploration of diverse evolutionary and conservation questions.

The koala, the only extant species of the marsupial family Phascolarctidae, is classified as 'vulnerable' due to habitat loss and widespread disease. We sequenced the koala genome, producing a complete and contiguous marsupial reference genome, including centromeres. We reveal that the koala's ability to detoxify eucalypt foliage may be due to expansions within a cytochrome P450 gene family, and its ability to smell, taste and moderate ingestion of plant secondary metabolites may be due to expansions in the vomeronasal and taste receptors. We characterized novel lactation proteins that protect young in the pouch and annotated immune genes important for response to chlamydial disease. Historical demography showed a substantial population crash coincident with the decline of Australian megafauna, while contemporary populations had biogeographic boundaries and increased inbreeding in populations affected by historic translocations. We identified genetically diverse populations that require habitat corridors and instituting of translocation programs to aid the koala's survival in the wild.

Here we present a new method for the identification of recurrent genomic entities that play a causative role in the onset of disease. Our approach is particularly amenable for the analyses high-throughput sequencing data. Existing approaches often follow a bottom-up approach where taxonomic determination necessarily takes place before associations to disease can be determined; naturally failing to establish the causality of novel pathogens not present in reference databases. To overcome this intrinsic limitation, we have developed a species-agnostic top-bottom approach that clusters sequences and identifies co-occurrence in multiple patients, associates recurrent sequences to disease and, finally, determines the taxonomic content where existing knowledge permits. We analysed 686 sequencing libraries from 252 cancer specimens and 56 controls. Recurrent sequences were statistically associated to biological, methodological and technical features to identify novel pathogens and contaminants stemming from laboratory reagents. We provide examples of identified inhabitants of the healthy tissue flora, known experimental contaminants and uncharacterised sequences that co-occur with high statistical significance with disease. The latter represent a category that can only be addressed by a species-independent approach. Thus, our method helps to chart the unknown sequence-space where novel pathogens can be identified.

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.

Outbreaks of zoonotic diseases in humans and livestock are not uncommon, and an important component in containment of such emerging viral diseases is rapid and reliable diagnostics. Such methods are often PCR-based and hence require the availability of sequence data from the pathogen. Rattus norvegicus (R. norvegicus) is a known reservoir for important zoonotic pathogens. Transmission may be direct via contact with the animal, for example, through exposure to its faecal matter, or indirectly mediated by arthropod vectors. Here we investigated the viral content in rat faecal matter (n= 29) collected from two continents by analyzing 2.2 billion next-generation sequencing reads derived from both DNA and RNA. Among other virus families, we found sequences from members of the Picornaviridae to be abundant in the microbiome of all the samples. Here we describe the diversity of the picornavirus-like contigs including near-full-length genomes closely related to the Boone cardiovirus and Theiler's encephalomyelitis virus. From this study, we conclude that picornaviruses within R. norvegicus are more diverse than previously recognized. The virome of R. norvegicus should be investigated further to assess the full potential for zoonotic virus transmission.

Virus discovery from high throughput sequencing data often follows a bottom-up approach where taxonomic annotation takes place prior to association to disease. Albeit effective in some cases, the approach fails to detect novel pathogens and remote variants not present in reference databases. We have developed a species independent pipeline that utilises sequence clustering for the identification of nucleotide sequences that co-occur across multiple sequencing data instances. We applied the workflow to 686 sequencing libraries from 252 cancer samples of different cancer and tissue types, 32 non-template controls, and 24 test samples. Recurrent sequences were statistically associated to biological, methodological or technical features with the aim to identify novel pathogens or plausible contaminants that may associate to a particular kit or method. We provide examples of identified inhabitants of the healthy tissue flora as well as experimental contaminants. Unmapped sequences that co-occur with high statistical significance potentially represent the unknown sequence space where novel pathogens can be identified.

Herbarium genomics is proving promising as next-generation sequencing approaches are well suited to deal with the usually fragmented nature of archival DNA. We show that routine assembly of partial plastome sequences from herbarium specimens is feasible, from total DNA extracts and with specimens up to 146 years old. We use genome skimming and an automated assembly pipeline, Iterative Organelle Genome Assembly, that assembles paired-end reads into a series of candidate assemblies, the best one of which is selected based on likelihood estimation. We used 93 specimens from 12 different Angiosperm families, 73 of which were from herbarium material with ages up to 146 years old. For 84 specimens, a sufficient number of paired-end reads were generated (in total 9.4 × 1012 nucleotides), yielding successful plastome assemblies for 74 specimens. Those derived from herbarium specimens have lower fractions of plastome-derived reads compared with those from fresh and silica-gel-dried specimens, but total herbarium assembly lengths are only slightly shorter. Specimens from wet-tropical conditions appear to have a higher number of contigs per assembly and lower N50 values. We find no significant correlation between plastome coverage and nuclear genome size (C value) in our samples, but the range of C values included is limited. Finally, we conclude that routine plastome sequencing from herbarium specimens is feasible and cost-effective (compared with Sanger sequencing or plastome-enrichment approaches), and can be performed with limited sample destruction.