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 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.

Correction: Phylogeography of the iconic Australian red-tailed black-cockatoo (Calyptorhynchus banksii) and implications for its conservation.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Phylogeography of the iconic Australian red-tailed black-cockatoo (Calyptorhynchus banksii) and implications for its conservation.

Advances in sequencing technologies have revolutionized wildlife conservation genetics. Analysis of genomic data sets can provide high-resolution estimates of genetic structure, genetic diversity, gene flow, and evolutionary history. These data can be used to characterize conservation units and to effectively manage the genetic health of species in a broad evolutionary context. Here we utilize thousands of genome-wide single-nucleotide polymorphisms (SNPs) and mitochondrial DNA to provide the first genetic assessment of the Australian red-tailed black-cockatoo (Calyptorhynchus banksii), a widespread bird species comprising populations of varying conservation concern. We identified five evolutionarily significant units, which are estimated to have diverged during the Pleistocene. These units are only partially congruent with the existing morphology-based subspecies taxonomy. Genetic clusters inferred from mitochondrial DNA differed from those based on SNPs and were less resolved. Our study has a range of conservation and taxonomic implications for this species. In particular, we provide advice on the potential genetic rescue of the Endangered and restricted-range subspecies C. b. graptogyne, and propose that the western C. b. samueli population is diagnosable as a separate subspecies. The results of our study highlight the utility of considering the phylogeographic relationships inferred from genome-wide SNPs when characterizing conservation units and management priorities, which is particularly relevant as genomic data sets become increasingly accessible.

Correction to "Insights into the Effects of Violating Statistical Assumptions for Dimensionality Reduction for Chemical "-omics" Data with Multiple Explanatory Variables".

[This corrects the article DOI: 10.1021/acsomega.3c01613.].

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.

A forensically validated genetic toolkit for the species and lineage identification of the highly trafficked shingleback lizard (Tiliqua rugosa).

Shingleback lizards (Tiliqua rugosa) are among the most trafficked native fauna from Australia in the illegal pet trade. There are four morphologically recognised subspecies of shinglebacks, all with differing overseas market values. Shinglebacks from different geographic locales are often trafficked and housed together, which may complicate identifying the State jurisdiction where the poaching event occurred. Additionally, shinglebacks can be housed and trafficked with other species within the same genus, which may complicate DNA analysis, especially in scenarios where indirect evidence (e.g. swabs, faeces) is taken for analysis. In this study, a forensic genetic toolkit was designed and validated to target shingleback DNA for species identification and geographic origin. To do this, field sampling across Australia was conducted to expand the phylogeographic sampling of shinglebacks across their species range and include populations suspected to be poaching hotspots. A commonly used universal reptile primer set (ND4/LEU) was then validated for use in forensic casework related to the genus Tiliqua. Two additional ND4 primer sets were designed and validated. The first primer set was designed and demonstrated to preferentially amplify an ∼510 bp region of the genus Tiliqua over other reptiles and builds on existing data to expand the available phylogeographic database. The second primer set was designed and demonstrated to solely amplify an ∼220 bp region of T. rugosa ND4 over any other reptile species. Through the validation process, all primers were demonstrated to amplify T. rugosa DNA from a variety of sample types (e.g. degraded, low quality and mixed). Two of the primer sets were able to distinguish the genetic lineage of T. rugosa from the phylogeographic database. This work provides the first forensically validated toolkit and phylogeographic genetic database for Squatmate lizards.

Insights into the Effects of Violating Statistical Assumptions for Dimensionality Reduction for Chemical "-omics" Data with Multiple Explanatory Variables.

Biological volatilome analysis is inherently complex due to the considerable number of compounds (i.e., dimensions) and differences in peak areas by orders of magnitude, between and within compounds found within datasets. Traditional volatilome analysis relies on dimensionality reduction techniques which aid in the selection of compounds that are considered relevant to respective research questions prior to further analysis. Currently, compounds of interest are identified using either supervised or unsupervised statistical methods which assume the data residuals are normally distributed and exhibit linearity. However, biological data often violate the statistical assumptions of these models related to normality and the presence of multiple explanatory variables which are innate to biological samples. In an attempt to address deviations from normality, volatilome data can be log transformed. However, whether the effects of each assessed variable are additive or multiplicative should be considered prior to transformation, as this will impact the effect of each variable on the data. If assumptions of normality and variable effects are not investigated prior to dimensionality reduction, ineffective or erroneous compound dimensionality reduction can impact downstream analyses. It is the aim of this manuscript to assess the impact of single and multivariable statistical models with and without the log transformation to volatilome dimensionality reduction prior to any supervised or unsupervised classification analysis. As a proof of concept, Shingleback lizard (Tiliqua rugosa) volatilomes were collected across their species distribution and from captivity and were assessed. Shingleback volatilomes are suspected to be influenced by multiple explanatory variables related to habitat (Bioregion), sex, parasite presence, total body volume, and captive status. This work determined that the exclusion of relevant multiple explanatory variables from analysis overestimates the effect of Bioregion and the identification of significant compounds. The log transformation increased the number of compounds that were identified as significant, as did analyses that assumed that residuals were normally distributed. Among the methods considered in this work, the most conservative form of dimensionality reduction was achieved through analyzing untransformed data using Monte Carlo tests with multiple explanatory variables.

Novel insights into the phylogenetic relationships of the endangered marsupial genus Potorous.

The three extant potoroo species of the marsupial genus Potorous -Potorous tridactylus, P. longipes and P. gilbertii - are all of conservation concern due to introduced predators and habitat loss associated with the European settlement of Australia. Robust phylogenies can be useful to inform conservation management, but past phylogenetic studies on potoroos have been unable to fully resolve relationships within the genus. Here, a multi-locus approach was employed, using three mitochondrial DNA (mtDNA): NADH dehydrogenase subunit 2, cytochrome c oxidase subunit 1 and 12S rRNA and four nuclear DNA (nuDNA) gene regions: breast and ovarian cancer susceptibility gene, recombination activating gene-1, apolipoprotein B and omega globin. This was coupled with widespread geographic sampling of the broadly distributed P. tridactylus, to investigate the phylogenetic relationships within this genus. Analyses of the mtDNA identified five distinct and highly divergent lineages including, P. longipes, P. gilbertii and three distinct lineages within P. tridactylus (northern mainland, southern mainland and Tasmanian). P. tridactylus was paraphyletic with the P. gilbertii lineage, suggesting that cryptic taxa may exist within P. tridactylus. NuDNA sequences lacked the resolution of mtDNA. Although they resolved the three currently recognised species, they were unable to differentiate lineages within P. tridactylus. Current management of P. tridactylus as two sub-species (mainland and Tasmania) does not recognise the full scope of genetic diversity within this species, especially that of the mainland populations. Until data from more informative nuDNA markers are available, we recommend this species be managed as the following three subspecies: Potorous tridactylus tridactylus (southern Queensland and northern New South Wales); Potorous tridactylus trisulcatus (southern New South Wales and Victoria) Potorous tridactylus apicalis (Tasmania). Molecular dating estimated that divergences within Potorous occurred in the late Miocene through to the early Pliocene.

Museum specimens provide reliable SNP data for population genomic analysis of a widely distributed but threatened cockatoo species.

Natural history museums harbour a plethora of biological specimens which are of potential use in population and conservation genetic studies. Although technical advancements in museum genomics have enabled genome-wide markers to be generated from aged museum specimens, the suitability of these data for robust biological inference is not well characterized. The aim of this study was to test the utility of museum specimens in population and conservation genomics by assessing the biological and technical validity of single nucleotide polymorphism (SNP) data derived from such samples. To achieve this, we generated thousands of SNPs from 47 red-tailed black cockatoo (Calyptorhychus banksii) traditional museum samples (i.e. samples that were not collected with the primary intent of DNA analysis) and 113 fresh tissue samples (cryopreserved liver/muscle) using a restriction site-associated DNA marker approach (DArTseq™ ). Thousands of SNPs were successfully generated from most of the traditional museum samples (with a mean age of 44 years, ranging from 5 to 123 years), although 38% did not provide useful data. These SNPs exhibited higher error rates and contained significantly more missing data compared with SNPs from fresh tissue samples, likely due to considerable DNA fragmentation. However, based on simulation results, the level of genotyping error had a negligible effect on inference of population structure in this species. We did identify a bias towards low diversity SNPs in older samples that appears to compromise temporal inferences of genetic diversity. This study demonstrates the utility of a RADseq-based method to produce reliable genome-wide SNP data from traditional museum specimens.

Unmasking the complexity of species identification in Australasian flying-foxes.

Pteropus (flying-foxes) are a speciose group of non-echolocating large bats, with five extant Australian species and 24 additional species distributed amongst the Pacific Islands. In 2015, an injured flying-fox with unusual facial markings was found in Sydney, Australia, following severe and widespread storms. Based on an initial assessment, the individual belonged to Pteropus but could not be readily identified to species. As a consequence, four hypotheses for its identification/origin were posited: the specimen represented (1) an undescribed Australian species; or (2) a morphological variant of a recognised Australian species; or (3) a hybrid individual; or (4) a vagrant from the nearby Southwest Pacific Islands. We used a combination of morphological and both mitochondrial- and nuclear DNA-based identification methods to assess these hypotheses. Based on the results, we propose that this morphologically unique Pteropus most likely represents an unusual P. alecto (black flying-fox) potentially resulting from introgression from another Pteropus species. Unexpectedly, this individual, and the addition of reference sequence data from newly vouchered specimens, revealed a previously unreported P. alecto mitochondrial DNA lineage. This lineage was distinct from currently available haplotypes. It also suggests long-term hybridisation commonly occurs between P. alecto and P. conspicillatus (spectacled flying-fox). This highlights the importance of extensive reference data, and the inclusion of multiple vouchered specimens for each species to encompass both intraspecific and interspecific variation to provide accurate and robust species identification. Moreover, our additional reference data further demonstrates the complexity of Pteropus species relationships, including hybridisation, and potential intraspecific biogeographical structure that may impact on their management and conservation.