Reversing the decline of threatened koala (Phascolarctos cinereus) populations in New South Wales: Using genomics to enhance conservation outcomes

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.

Koala Genome Survey: An Open Data Resource to Improve Conservation Planning

Genome sequencing is a powerful tool that can inform the management of threatened species. Koalas (Phascolarctos cinereus) are a globally recognized species that captured the hearts and minds of the world during the 2019/2020 Australian megafires. In 2022, koalas were listed as ‘Endangered’ in Queensland, New South Wales, and the Australian Capital Territory. Populations have declined because of various threats such as land clearing, habitat fragmentation, and disease, all of which are exacerbated by climate change. Here, we present the Koala Genome Survey, an open data resource that was developed after the Australian megafires. A systematic review conducted in 2020 demonstrated that our understanding of genomic diversity within koala populations was scant, with only a handful of SNP studies conducted. Interrogating data showed that only 6 of 49 New South Wales areas of regional koala significance had meaningful genome-wide data, with only 7 locations in Queensland with SNP data and 4 locations in Victoria. In 2021, we launched the Koala Genome Survey to generate resequenced genomes across the Australian east coast. We have publicly released 430 koala genomes (average coverage: 32.25X, range: 11.3–66.8X) on the Amazon Web Services Open Data platform to accelerate research that can inform current and future conservation planning.

Future-proofing the koala: synergising genomic and environmental data for effective species management

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-300 kya), a major climatic reorganization 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 threatened 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.

How microbiomes can help inform conservation: landscape characterisation of gut microbiota helps shed light on additional population structure in a specialist folivore

BACKGROUND

The koala (Phascolarctos cinereus), an iconic yet endangered specialised folivore experiencing widespread decline across Australia, is the focus of many conservation programs. Whilst animal translocation and progressive conservation strategies such as faecal inoculations may be required to bring this species back from the brink of extinction, insight into the variation of host-associated gut microbiota and the factors that shape this variation are fundamental for their success. Despite this, very little is known about the landscape variability and factors affecting koala gut microbial community dynamics. We used large scale field surveys to evaluate the variation and diversity of koala gut microbiotas and compared these diversity patterns to those detected using a population genetics approach. Scat samples were collected from five locations across South East Queensland with microbiota analysed using 16S rRNA gene amplicon sequencing.

RESULTS

Across the landscape koala gut microbial profiles showed large variability, with location having a large effect on bacterial community composition and bacterial diversity. Certain bacteria were found to be significantly differentially abundant amongst locations; koalas from Noosa showed a depletion in two bacterial orders (Gastranaerophilales and Bacteroidales) which have been shown to provide beneficial properties to their host. Koala gut microbial patterns were also not found to mirror population genetic patterns, a molecular tool often used to design conservation initiatives.

CONCLUSIONS

Our data shows that koala gut microbiotas are extremely variable across the landscape, displaying complex micro- and macro- spatial variation. By detecting locations which lack certain bacteria we identified koala populations that may be under threat from future microbial imbalance or dysbiosis. Additionally, the mismatching of gut microbiota and host population genetic patterns exposed important population structure that has previously gone undetected across South East Queensland. Overall, this baseline data highlights the importance of integrating microbiota research into conservation biology in order to guide successful conservation programs such as species translocation and the implementation of faecal inoculations.

Physiological stress levels in wild koala sub-populations facing anthropogenic induced environmental trauma and disease

Australian small mammals such as koalas must cope with immense pressure from anthropogenic induced stressors or trauma such as bushfires, vehicle collision impacts and habitat disturbance and land clearance. In addition, they must cope with diseases such as chlamydia. To date, there is no published literature on physiological stress levels in wild koala populations compared with identified environmental stressors. This study investigated physiological stress levels within sub-populations of wild koalas encountering environmental trauma and disease from New South Wales (NSW), Queensland (QLD) and South Australia (SA). Physiological stress was determined using a faecal glucocorticoid (or cortisol) metabolites (FGMs) enzyme-immunoassay (EIA) from 291 fresh faecal samples collected from wild koalas at the point of rescue. A healthy breeding sub-population from a forest reserve in QLD acted as a control group. Clearance of prime Eucalyptus habitat had the largest impact on FGMs, followed by bushfire related factors (e.g. flat demeanour, dehydration and burns injury). Koalas with other sources of physical injury (dog-attacks and vehicle collisions) and disease (chlamydia) also had higher FGMs compared to healthy wild koalas. Healthy wild koalas expressed the lowest median levels of FGMs. Overall, the results highlight that anthropogenic-induced stressors tend to increase physiological stress in wild koalas. Thus, the ultimate stressors such as habitat clearance and bush fire events could increase the incidences of proximate stressors such as dog attacks and vehicle collisions, and increase risks of foliage shortage, diseases and mortality. Therefore, there is need for ecological monitoring, conservation management actions and policy changes to curb the koala population crisis, especially within on-going and future land and road development programs.

Disentangling the mechanisms of mate choice in a captive koala population

Successful captive breeding programs are crucial to the long-term survival of many threatened species. However, pair incompatibility (breeding failure) limits sustainability of many captive populations. Understanding whether the drivers of this incompatibility are behavioral, genetic, or a combination of both, is crucial to improving breeding programs. We used 28 years of pairing data from the San Diego Zoo koala colony, plus genetic analyses using both major histocompatibility complex (MHC)-linked and non-MHC-linked microsatellite markers, to show that both genetic and non-genetic factors can influence mating success. Male age was reconfirmed to be a contributing factor to the likelihood of a koala pair copulating. This trend could also be related to a pair's age difference, which was highly correlated with male age in our dataset. Familiarity was reconfirmed to increase the probability of a successful copulation. Our data provided evidence that females select mates based on MHC and genome-wide similarity. Male heterozygosity at MHC class II loci was associated with both pre- and post-copulatory female choice. Genome-wide similarity, and similarity at the MHC class II DAB locus, were also associated with female choice at the post-copulatory level. Finally, certain MHC-linked alleles were associated with either increased or decreased mating success. We predict that utilizing a variety of behavioral and MHC-dependent mate choice mechanisms improves female fitness through increased reproductive success. This study highlights the complexity of mate choice mechanisms in a species, and the importance of ascertaining mate choice mechanisms to improve the success of captive breeding programs.

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.

Fresh is best: Accurate SNP genotyping from koala scats

Maintaining genetic diversity is a crucial component in conserving threatened species. For the iconic Australian koala, there is little genetic information on wild populations that is not either skewed by biased sampling methods (e.g., sampling effort skewed toward urban areas) or of limited usefulness due to low numbers of microsatellites used. The ability to genotype DNA extracted from koala scats using next‐generation sequencing technology will not only help resolve location sample bias but also improve the accuracy and scope of genetic analyses (e.g., neutral vs. adaptive genetic diversity, inbreeding, and effective population size). Here, we present the successful SNP genotyping (1272 SNP loci) of koala DNA extracted from scat, using a proprietary DArTseq^™ protocol. We compare genotype results from two‐day‐old scat DNA and 14‐day‐old scat DNA to a blood DNA template, to test accuracy of scat genotyping. We find that DNA from fresher scat results in fewer loci with missing information than DNA from older scat; however, 14‐day‐old scat can still provide useful genetic information, depending on the research question. We also find that a subset of 209 conserved loci can accurately identify individual koalas, even from older scat samples. In addition, we find that DNA sequences identified from scat samples through the DArTseq^™ process can provide genetic identification of koala diet species, bacterial and viral pathogens, and parasitic organisms.