Genomic comparisons reveal biogeographic and anthropogenic impacts in the koala (Phascolarctos cinereus): a dietary-specialist species distributed across heterogeneous environments

Health Assessments of Koalas after Wildfire: A Temporal Comparison of Rehabilitated and Non-Rescued Resident Individuals

Many koalas (Phascolarctos cinereus) required rehabilitation after the 2019/20 Australian megafires. Little is known about how the post-release health of rehabilitated koalas compares to non-rescued resident koalas. We evaluated health parameters in rehabilitated koalas and resident koalas in burnt and unburnt habitat in southern New South Wales, Australia. Health checks were undertaken within six weeks of fire (rehabilitated group), 5-9 months post-fire and 12-16 months post-fire. Body condition improved significantly over time in rehabilitated koalas, with similar condition between all groups at 12-16 months. Rehabilitated koalas therefore gained body condition at similar rates to koalas who remained and survived in the wild. The prevalence of Chlamydia pecorum was also similar between groups and timepoints, suggesting wildfire and rehabilitation did not exacerbate disease in this population. While there was some variation in measured serum biochemistry and haematology parameters between groups and timepoints, most were within normal reference ranges. Our findings show that koalas were generally healthy at the time of release and when recaptured nine months later. Landscapes in the Monaro region exhibiting a mosaic of burn severity can support koalas, and rehabilitated koala health is not compromised by returning them to burnt habitats 4-6 months post-fire.

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.

Evaluating the genetic consequences of population subdivision as it unfolds and how to best mitigate them: A rare story about koalas

The genetic consequences of the subdivision of populations are regarded as significant to long-term evolution, and research has shown that the scale and speed at which this is now occurring is critically reducing the adaptive potential of most species which inhabit human-impacted landscapes. Here, we provide a rare and, to our knowledge, the first analysis of this process while it is happening and demonstrate a method of evaluating the effect of mitigation measures such as fauna crossings. We did this by using an extensive genetic data set collected from a koala population which was intensely monitored during the construction of linear transport infrastructure which resulted in the subdivision of their population. First, we found that both allelic richness and effective population size decreased through the process of population subdivision. Second, we predicted the extent to which genetic drift could impact genetic diversity over time and showed that after only 10 generations the resulting two subdivided populations could experience between 12% and 69% loss in genetic diversity. Lastly, using forward simulations we estimated that a minimum of eight koalas would need to disperse from each side of the subdivision per generation to maintain genetic connectivity close to zero but that 16 koalas would ensure that both genetic connectivity and diversity remained unchanged. These results have important consequences for the genetic management of species in human-impacted landscapes by showing which genetic metrics are best to identify immediate loss in genetic diversity and how to evaluate the effectiveness of any mitigation measures.

Susceptibility to a sexually transmitted disease in a wild koala population shows heritable genetic variance but no inbreeding depression

The koala, one of the most iconic Australian wildlife species, is facing several concomitant threats that are driving population declines. Some threats are well known and have clear methods of prevention (e.g., habitat loss can be reduced with stronger land-clearing control), whereas others are less easily addressed. One of the major current threats to koalas is chlamydial disease, which can have major impacts on individual survival and reproduction rates and can translate into population declines. Effective management strategies for the disease in the wild are currently lacking, and, to date, we know little about the determinants of individual susceptibility to disease. Here, we investigated the genetic basis of variation in susceptibility to chlamydia using one of the most intensively studied wild koala populations. We combined data from veterinary examinations, chlamydia testing, genetic sampling and movement monitoring. Out of our sample of 342 wild koalas, 60 were found to have chlamydia. Using genotype information on 5007 SNPs to investigate the role of genetic variation in determining disease status, we found no evidence of inbreeding depression, but a heritability of 0.11 (95% CI: 0.06-0.23) for the probability that koalas had chlamydia. Heritability of susceptibility to chlamydia could be relevant for future disease management, as it suggests adaptive potential for the population.

Geographic patterns of koala retrovirus genetic diversity, endogenization, and subtype distributions

Koala retrovirus (KoRV) subtype A (KoRV-A) is currently in transition from exogenous virus to endogenous viral element, providing an ideal system to elucidate retroviral-host coevolution. We characterized KoRV geography using fecal DNA from 192 samples across 20 populations throughout the koala's range. We reveal an abrupt change in KoRV genetics and incidence at the Victoria/New South Wales state border. In northern koalas, pol gene copies were ubiquitously present at above five per cell, consistent with endogenous KoRV. In southern koalas, pol copies were detected in only 25.8% of koalas and always at copy numbers below one, while the env gene was detected in all animals and in a majority at copy numbers above one per cell. These results suggest that southern koalas carry partial endogenous KoRV-like sequences. Deep sequencing of the env hypervariable region revealed three putatively endogenous KoRV-A sequences in northern koalas and a single, distinct sequence present in all southern koalas. Among northern populations, env sequence diversity decreased with distance from the equator, suggesting infectious KoRV-A invaded the koala genome in northern Australia and then spread south. The exogenous KoRV subtypes (B to K), two novel subtypes, and intermediate subtypes were detected in all northern koala populations but were strikingly absent from all southern animals tested. Apart from KoRV subtype D, these exogenous subtypes were generally locally prevalent but geographically restricted, producing KoRV genetic differentiation among northern populations. This suggests that sporadic evolution and local transmission of the exogenous subtypes have occurred within northern Australia, but this has not extended into animals within southern Australia.

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.

Genome-wide SNPs detect fine-scale genetic structure in threatened populations of squirrel glider Petaurus norfolcensis

Australian arboreal mammals are experiencing significant population declines, particularly due to land clearing and resulting habitat fragmentation. The squirrel glider, Petaurus norfolcensis, is a threatened species in New South Wales, with a stronghold population in the Lake Macquarie Local Government Area (LGA) where fragmentation due to urbanization is an ongoing problem for the species conservation. Here we report on the use of squirrel glider mitochondrial (385 bp cytochrome b gene, 70 individuals) and nuclear DNA (6,834 SNPs, 87 individuals) markers to assess their population genetic structure and connectivity across 14 locations sampled in the Lake Macquarie LGA. The mitochondrial DNA sequences detected evidence of a historical genetic bottleneck, while the genome-wide SNPs detected significant population structure in the Lake Macquarie squirrel glider populations at scales as fine as one kilometer. There was no evidence of inbreeding within patches, however there were clear effects of habitat fragmentation and biogeographical barriers on gene flow. A least cost path analysis identified thin linear corridors that have high priority for conservation. These areas should be protected to avoid further isolation of squirrel glider populations and the loss of genetic diversity through genetic drift.

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.

Testing the effectiveness of genetic monitoring using genetic non-invasive sampling

Effective conservation requires accurate data on population genetic diversity, inbreeding, and genetic structure. Increasingly, scientists are adopting genetic non-invasive sampling (gNIS) as a cost-effective population-wide genetic monitoring approach. gNIS has, however, known limitations which may impact the accuracy of downstream genetic analyses. Here, using high-quality single nucleotide polymorphism (SNP) data from blood/tissue sampling of a free-ranging koala population (n = 430), we investigated how the reduced SNP panel size and call rate typical of genetic non-invasive samples (derived from experimental and field trials) impacts the accuracy of genetic measures, and also the effect of sampling intensity on these measures. We found that gNIS at small sample sizes (14% of population) can provide accurate population diversity measures, but slightly underestimated population inbreeding coefficients. Accurate measures of internal relatedness required at least 33% of the population to be sampled. Accurate geographic and genetic spatial autocorrelation analysis requires between 28% and 51% of the population to be sampled. We show that gNIS at low sample sizes can provide a powerful tool to aid conservation decision-making and provide recommendations for researchers looking to apply these techniques to free-ranging systems.

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.

Koala tree selection in a mixed-tenure landscape and post-fire implications

Abstract ContextThe more frequent and intense bushfires predicted under climate change are likely to heavily impact koalas living inside protected areas and intact forests, which makes koala habitat in asset protection zones more important. Understanding how koalas use different habitats in a mixed-tenure landscape can inform effective conservation management. AimsThe aims of this study were to determine (1) the species and size of trees used by koalas, (2) whether choice of daytime trees was affected by available proportions of tree species and sizes in koala home ranges, and (3) whether developed land was used as frequently as native forest. MethodsTree use by koalas was determined by radio-tracking 10 koalas for an average of 12 months and recording the species and diameter of trees they occupied. To compare the proportions of tree species used by koalas with those available, tree availability was measured using random quadrats in the home ranges of five koalas that used forested areas extensively. The habitats used by 10 koalas were classified as native forest or developed land to investigate the importance of human-modified areas. Key resultsEucalyptus punctata was a preferred tree species, but each individual used four to nine species. Tree species were not chosen solely on the basis of available proportions; some species were selected preferentially. Half of the koalas used more developed land than native forest. Koalas preferentially used trees with larger diameters than the mean of available trees, and selected larger trees in developed areas than in native forest and when re-using trees. ConclusionsDespite the higher availability of trees in protected native forest at the site, 70% of the koalas used developed land in and around asset protection zones. Koalas whose home ranges were geographically close, but had different soil types and vegetation communities, were able to utilise different tree species. ImplicationsKoalas would benefit from protection of remnant native forests containing preferred trees on shale cap soil and conservation of native vegetation corridors along fence lines and in paddocks in developed areas because they are valuable resources and connect patchy landscapes. Protecting koalas on developed land improves their likelihood of surviving bushfires, allowing recolonisation of surrounding protected areas.

Comparing Genomic Signatures of Selection Between the Abbassa Strain and Eight Wild Populations of Nile Tilapia (Oreochromis niloticus) in Egypt

Domestication to captive rearing conditions, along with targeted selective breeding have genetic consequences that vary from those in wild environments. Nile tilapia (Oreochromis niloticus) is one of the most translocated and farmed aquaculture species globally, farmed throughout Asia, North and South America, and its African native range. In Egypt, a breeding program established the Abbassa Strain of Nile tilapia (AS) in 2002 based on local broodstock sourced from the Nile River. The AS has been intensively selected for growth and has gone through genetic bottlenecks which have likely shifted levels and composition of genetic diversity within the strain. Consequently, there are questions on the possible genetic impact AS escapees may have on endemic populations of Nile tilapia. However, to date there have been no genetic studies comparing genetic changes in the domesticated AS to local wild populations. This study used 9,827 genome-wide SNPs to investigate population genetic structure and signatures of selection in the AS (generations 9–11) and eight wild Nile tilapia populations from Egypt. SNP analyses identified two major genetic clusters (captive and wild populations), with wild populations showing evidence of isolation-by-distance among the Nile Delta and upstream riverine populations. Between genetic clusters, approximately 6.9% of SNPs were identified as outliers with outliers identified on all 22 O. niloticus chromosomes. A lack of localized outlier clustering on the genome suggests that no genes of major effect were presently detected. The AS has retained high levels of genetic diversity (H_{o_All} = 0.21 ± 0.01; H_{e_All} = 0.23 ± 0.01) when compared to wild populations (H_{o_All} = 0.18 ± 0.01; H_{e_All} = 0.17 ± 0.01) after 11 years of domestication and selective breeding. Additionally, 565 SNPs were unique within the AS line. While these private SNPs may be due to domestication signals or founder effects, it is suspected that introgression with blue tilapia (Oreochromis aureus) has occurred. This study highlights the importance of understanding the effects of domestication in addition to wild population structure to inform future management and dissemination decisions. Furthermore, by conducting a baseline genetic study of wild populations prior to the dissemination of a domestic line, the effects of aquaculture on these populations can be monitored over time.

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.

Histological survey for oxalate nephrosis in Victorian koalas (Phascolarctos cinereus)

The Mount Lofty Ranges koala (Phascolarctos cinereus) population in South Australia has a high prevalence of the renal disease oxalate nephrosis, for which an underlying genetic cause is suspected. South Australian koalas primarily originate from those in French Island, Victoria; however, oxalate nephrosis has not previously been reported in Victorian koalas. Examination of kidney tissue sections from 63 koalas across Victoria found that nine koalas were affected by oxalate nephrosis (14.3%). These included 2/5 koalas from French Island (40%), 4/14 koalas from the western regions (29%), 2/11 Raymond Island koalas (18%), and 1/13 Cape Otway koalas (8%). There were no cases of oxalate nephrosis identified in the Strzelecki koalas (n = 12). These findings suggest that oxalate nephrosis occurs in koalas from French Island and populations that have received significant influx of koalas from French Island, but not in the Strzelecki region, which has little to no French Island input. This lends support to the theory that an inherited abnormality of oxalate metabolism could underlie the high prevalence of oxalate nephrosis in the Mount Lofty Ranges koala population, and molecular investigations are currently underway to investigate a genetic cause.

Inbreeding and disease avoidance in a free-ranging koala population

Habitat destruction and fragmentation are increasing globally, forcing surviving species into small, isolated populations. Isolated populations typically experience heightened inbreeding risk and associated inbreeding depression and population decline; although individuals in these populations may mitigate these risks through inbreeding avoidance strategies. For koalas, as dietary specialists already under threat in the northern parts of their range, increased habitat fragmentation and associated inbreeding costs are of great conservation concern. Koalas are known to display passive inbreeding avoidance through sex-biased dispersal, although population isolation will reduce dispersal pathways. We tested whether free-ranging koalas display active inbreeding avoidance behaviours. We used VHF tracking data, parentage reconstruction, and veterinary examination results to test whether free-ranging female koalas avoid mating with (a) more closely related males; and (b) males infected with sexually transmitted Chlamydia pecorum. We found no evidence that female koalas avoid mating with relatively more related available mates. In fact, as the relatedness of potential mates increases, so did inbreeding events. We also found no evidence that female koalas can avoid mating with males infected with C. pecorum. The absence of active inbreeding avoidance mechanisms in koalas is concerning from a conservation perspective, as small, isolated populations may be at even higher risk of inbreeding depression than expected. At risk koala populations may require urgent conservation interventions to augment gene flow and reduce inbreeding risks. Similarly, if koalas are not avoiding mating with individuals with chlamydial disease, populations may be at higher risk from disease than anticipated, further impacting population viability.

Infection dynamics, dispersal, and adaptation: understanding the lack of recovery in a remnant frog population following a disease outbreak

Emerging infectious diseases can cause dramatic declines in wildlife populations. Sometimes, these declines are followed by recovery, but many populations do not recover. Studying differential recovery patterns may yield important information for managing disease-afflicted populations and facilitating population recoveries. In the late 1980s, a chytridiomycosis outbreak caused multiple frog species in Australia's Wet Tropics to decline. Populations of some species (e.g., Litoria nannotis) subsequently recovered, while others (e.g., Litoria dayi) did not. We examined the population genetics and current infection status of L. dayi, to test several hypotheses regarding the failure of its populations to recover: (1) a lack of individual dispersal abilities has prevented recolonization of previously occupied locations, (2) a loss of genetic variation has resulted in limited adaptive potential, and (3) L. dayi is currently adapting to chytridiomycosis. We found moderate-to-high levels of gene flow and diversity (Fst range: <0.01-0.15; minor allele frequency (MAF): 0.192-0.245), which were similar to previously published levels for recovered L. nannotis populations. This suggests that dispersal ability and genetic diversity do not limit the ability of L. dayi to recolonize upland sites. Further, infection intensity and prevalence increased with elevation, suggesting that chytridiomycosis is still limiting the elevational range of L. dayi. Outlier tests comparing infected and uninfected individuals consistently identified 18 markers as putatively under selection, and several of those markers matched genes that were previously implicated in infection. This suggests that L. dayi has genetic variation for genes that affect infection dynamics and may be undergoing adaptation.