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.

Genomic insights into the critically endangered King Island scrubtit.

Small, fragmented, or isolated populations are at risk of population decline due to fitness costs associated with inbreeding and genetic drift. The King Island scrubtit Acanthornis magna greeniana is a critically endangered subspecies of the nominate Tasmanian scrubtit A. m. magna, with an estimated population of <100 individuals persisting in three patches of swamp forest. The Tasmanian scrubtit is widespread in wet forests on mainland Tasmania. We sequenced the scrubtit genome using PacBio HiFi and undertook a population genomic study of the King Island and Tasmanian scrubtits using a double-digest restriction site-associated DNA (ddRAD) dataset of 5,239 SNP loci. The genome was 1.48 Gb long, comprising 1,518 contigs with an N50 of 7.715 Mb. King Island scrubtits formed one of four overall genetic clusters, but separated into three distinct subpopulations when analyzed independently of the Tasmanian scrubtit. Pairwise FST values were greater among the King Island scrubtit subpopulations than among most Tasmanian scrubtit subpopulations. Genetic diversity was lower and inbreeding coefficients were higher in the King Island scrubtit than all except one of the Tasmanian scrubtit subpopulations. We observed crown baldness in 8/15 King Island scrubtits, but 0/55 Tasmanian scrubtits. Six loci were significantly associated with baldness, including one within the DOCK11 gene which is linked to early feather development. Contemporary gene flow between King Island scrubtit subpopulations is unlikely, with further field monitoring required to quantify the fitness consequences of its small population size, low genetic diversity, and high inbreeding. Evidence-based conservation actions can then be implemented before the taxon goes extinct.

Phylogenomic and morphological data reveal hidden patterns of diversity in the national tree of Brazil, Paubrasilia echinata.

PREMISE: Paubrasilia echinata (common names, pau brasil, brazilwood) is the national tree of Brazil and an endangered species endemic to the Brazilian Atlantic Forest. Over its wide distribution of 2000 km, its leaflets morphology exhibits extensive plasticity. Three morphotypes are commonly identified based on leaf size, but it is unclear if they represent distinct taxa or a single polymorphic species. This study aims to clarify the taxonomic position of the three morphotypes to inform conservation decisions. METHODS: A morphometric study of leaf characters of herbarium specimens was coupled with genetic analyses using genotype-by-sequencing data. We used maximum-likelihood and coalescent methods to evaluate the phylogenetic and population structure of the species. We compared these with a morphological dendrogram built from hierarchical clustering. RESULTS: Two of the three morphotypes formed separately evolving lineages, the third morphotype formed two geographically separate lineages, and northern trees with intermediate leaf morphology formed a separate fifth lineage. Leaflet size varied by over 35-fold, and although morphological clustering generally matched the genetic patterns, there were some overlaps, highlighting the cryptic diversity within this group. CONCLUSIONS: Our genetic and morphological results provide some evidence that cultivated trees from different states in Brazil seem to have a limited genetic origin and do not reflect the broader genetic and geographical diversity of the species. As a result, more care is likely needed to preserve the overall genomic diversity of this endangered and iconic species.

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.

Distinct Genetic Structure Reflects Ploidy Level Differentiation in Newly Discovered, Extremely Small Populations of Xanthocyparis vietnamensis from Southwestern China.

Population genetic assessment is crucial for the conservation and management of threatened species. Xanthocyparis vietnamensis is an endangered species that is currently restricted to karst mountains in southwestern China and Vietnam. This rare conifer was first recorded in 2002 from northern Vietnam and then in 2013 from Guangxi, China, yet nothing is known about its genetic diversity nor ploidy level variation, although previous cytological study suggest that Vietnamese populations are tetraploids. There have been about 45 individuals found to date in Guangxi, China. Here, we genotyped 33 X. vietnamensis individuals using 20 newly developed, polymorphic microsatellite loci, to assess the genetic variability of its extremely small populations. The genetic diversity of X. vietnamensis (H E = 0.511) was lower than that of two other heliophile species, Calocedrus macrolepis and Fokienia hodginsii, which have similar distribution ranges. This is consistent with the signature of a genetic bottleneck detected in X. vietnamensis. Although the population genetic differentiation coefficient across loci is moderate (F ST = 0.125), STRUCTURE analysis revealed two distinct genetic clusters, namely the northern and southern population groups; DAPC analysis grouped the southern populations together in one cluster separate from the northern populations; AMOVA analysis detected a significant genetic differentiation between the two population groups (F RT = 0.089, p < 0.05), and BARRIER analysis detected a genetic barrier between them. Moreover, we detected differentiation in ploidy level between northern and southern populations, sampled individuals from the former and the later are all diploid and tetraploid cytotypes with mean genome sizes of 26.08 and 48.02 pg/2C, respectively. We deduced that heterogeneous geomorphology and historical events (e.g., human deforestation, Quaternary climate oscillations) may have contributed to population fragmentation and small population size in X. vietnamensis. Considering both genetic and ploidy level differentiation, we propose that two different management units (northern and southern) should be considered and a combination of in situ and ex situ conservation measures should be employed to preserve populations of this endangered species in southwestern China in the light of our findings.

DNA barcoding identifies cryptic animal tool materials.

Some animals fashion tools or constructions out of plant materials to aid foraging, reproduction, self-maintenance, or protection. Their choice of raw materials can affect the structure and properties of the resulting artifacts, with considerable fitness consequences. Documenting animals' material preferences is challenging, however, as manufacture behavior is often difficult to observe directly, and materials may be processed so heavily that they lack identifying features. Here, we use DNA barcoding to identify, from just a few recovered tool specimens, the plant species New Caledonian crows (Corvus moneduloides) use for crafting elaborate hooked stick tools in one of our long-term study populations. The method succeeded where extensive fieldwork using an array of conventional approaches-including targeted observations, camera traps, radio-tracking, bird-mounted video cameras, and behavioral experiments with wild and temporarily captive subjects-had failed. We believe that DNA barcoding will prove useful for investigating many other tool and construction behaviors, helping to unlock significant research potential across a wide range of study systems.

Understanding Historical Demographic Processes to Inform Contemporary Conservation of an Arid zone Specialist: The Yellow-Footed Rock-Wallaby.

Little genetic research has been undertaken on mammals across the vast expanse of the arid biome in Australia, despite continuing species decline and need for conservation management. Here, we evaluate the contemporary and historical genetic connectivity of the yellow-footed rock-wallaby, Petrogalexanthopusxanthopus, a threatened macropodid which inhabits rocky outcrops across the disconnected mountain range systems of the southern arid biome. We use 17 microsatellite loci together with mitochondrial control region data to determine the genetic diversity of populations and the evolutionary processes shaping contemporary population dynamics on which to base conservation recommendations. Our results indicate the highly fragmented populations have reduced diversity and limited contemporary gene flow, with most populations having been through population bottlenecks. Despite limited contemporary gene flow, the phylogeographic relationships of the mitochondrial control region indicate a lack of structure and suggests greater historical connectivity. This is an emerging outcome for mammals across this arid region. On the basis of our results, we recommend augmentation of populations of P. x.xanthopus, mixing populations from disjunct mountain range systems to reduce the chance of continued diversity loss and inbreeding depression, and therefore maximize the potential for populations to adapt and survive into the future.

Could do better! A high school market survey of fish labelling in Sydney, Australia, using DNA barcodes.

BACKGROUND: Processed seafood products are not readily identifiable based on physical characteristics, which leaves the industry vulnerable to high levels of product mislabelling (globally estimated at 5-30% mislabelled). This is both a food safety issue and a consumer protection issue as cheaper species could be substituted for more expensive species. DNA barcoding is proving to be a valuable tool for authentication of fish products. We worked with high school students to perform a market survey and subsequent species assessment via DNA barcoding to investigate the accuracy of fish product names used by retailers in Sydney, Australia. METHODS: Sixty-eight fish samples, sold under 50 different common names, were purchased anonymously from two retailers in Sydney. Each product name was recorded and reconciled with the Australian Fish Names Standard (AFNS). Samples were DNA barcoded and resulting sequences were deposited in the online Barcode of Life Data system using the simplified Student Data Portal interface. RESULTS: Forty percent of the fish names did not comply with the AFNS, however, half of these were either spelling errors or vendors supplied more information than the standard requires. The other half of the non-compliant samples were given common names not listed on the AFNS. Despite this lack of standardization, DNA barcode data confirmed the retailers' identifications for 93% of samples and 90% of species sampled. DISCUSSION: The level of mislabelling we report for Sydney retailers (7% of samples or 10% of species) compares favorably with the global rates of 5-30%, but unfavorably with the only previous DNA barcode fish authentication study for Australia, which found no confirmed mislabelling in Hobart. Our study sampled mostly Australian produce, only two retailers and no restaurants. Results of our limited sample suggest that although many Sydney fish retailers attempt to implement the voluntary fish name standards, the standards are inadequate. As Australia imports 75% of its seafood, and in other countries restaurants generally show lower levels of compliance than retailers, broader surveys are needed before generalizing these results. DNA barcoding is a powerful yet simple method supported by accessible online analytical tools. Incorporation of fish barcoding into high school science classes provided students with valuable firsthand experience in scientific research and drew together different strands of the NSW curriculum relating to genetics and sustainability. Given the techniques, equipment, and reagents are now readily accessible, we expect to see greater uptake of DNA barcoding technology by high schools, citizen scientists and consumer groups in Australia in future. However, there remains much scope for further development of DNA barcode diagnostics (both data and analytical methods) for commercial fish species.

Adaptation and conservation insights from the koala genome.

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.

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.

Analysis of Phylogenomic Tree Space Resolves Relationships Among Marsupial Families.

A fundamental challenge in resolving evolutionary relationships across the tree of life is to account for heterogeneity in the evolutionary signal across loci. Studies of marsupial mammals have demonstrated that this heterogeneity can be substantial, leaving considerable uncertainty in the evolutionary timescale and relationships within the group. Using simulations and a new phylogenomic data set comprising nucleotide sequences of 1550 loci from 18 of the 22 extant marsupial families, we demonstrate the power of a method for identifying clusters of loci that support different phylogenetic trees. We find two distinct clusters of loci, each providing an estimate of the species tree that matches previously proposed resolutions of the marsupial phylogeny. We also identify a well-supported placement for the enigmatic marsupial moles (Notoryctes) that contradicts previous molecular estimates but is consistent with morphological evidence. The pattern of gene-tree variation across tree-space is characterized by changes in information content, GC content, substitution-model adequacy, and signatures of purifying selection in the data. In a simulation study, we show that incomplete lineage sorting can explain the division of loci into the two tree-topology clusters, as found in our phylogenomic analysis of marsupials. We also demonstrate the potential benefits of minimizing uncertainty from phylogenetic conflict for molecular dating. Our analyses reveal that Australasian marsupials appeared in the early Paleocene, whereas the diversification of present-day families occurred primarily during the late Eocene and early Oligocene. Our methods provide an intuitive framework for improving the accuracy and precision of phylogenetic inference and molecular dating using genome-scale data.

Extensive genetic differentiation detected within a model marsupial, the tammar wallaby (Notamacropus eugenii).

The tammar wallaby (Notamacropus eugenii) is one of the most intensively studied of all macropodids and was the first Australasian marsupial to have its genome sequenced. However, comparatively little is known about genetic diversity and differentiation amongst the morphologically distinct allopatric populations of tammar wallabies found in Western (WA) and South Australia (SA). Here we compare autosomal and Y-linked microsatellite genotypes, as well as sequence data (~600 bp) from the mitochondrial DNA (mtDNA) control region (CR) in tammar wallabies from across its distribution. Levels of diversity at autosomal microsatellite loci were typically high in the WA mainland and Kangaroo Island (SA) populations (A = 8.9-10.6; He = 0.77-0.78) but significantly reduced in other endemic island populations (A = 3.8-4.1; He = 0.41-0.48). Autosomal and Y-linked microsatellite loci revealed a pattern of significant differentiation amongst populations, especially between SA and WA. The Kangaroo Island and introduced New Zealand population showed limited differentiation. Multiple divergent mtDNA CR haplotypes were identified within both SA and WA populations. The CR haplotypes of tammar wallabies from SA and WA show reciprocal monophyly and are highly divergent (14.5%), with levels of sequence divergence more typical of different species. Within WA tammar wallabies, island populations each have unique clusters of highly related CR haplotypes and each is most closely related to different WA mainland haplotypes. Y-linked microsatellite haplotypes show a similar pattern of divergence although levels of diversity are lower. In light of these differences, we suggest that two subspecies of tammar wallaby be recognized; Notamacropus eugenii eugenii in SA and N. eugenii derbianus in WA. The extensive neutral genetic diversity and inter-population differentiation identified within tammar wallabies should further increase the species value and usefulness as a model organism.

Understanding and monitoring the consequences of human impacts on intraspecific variation.

Intraspecific variation is a major component of biodiversity, yet it has received relatively little attention from governmental and nongovernmental organizations, especially with regard to conservation plans and the management of wild species. This omission is ill-advised because phenotypic and genetic variations within and among populations can have dramatic effects on ecological and evolutionary processes, including responses to environmental change, the maintenance of species diversity, and ecological stability and resilience. At the same time, environmental changes associated with many human activities, such as land use and climate change, have dramatic and often negative impacts on intraspecific variation. We argue for the need for local, regional, and global programs to monitor intraspecific genetic variation. We suggest that such monitoring should include two main strategies: (i) intensive monitoring of multiple types of genetic variation in selected species and (ii) broad-brush modeling for representative species for predicting changes in variation as a function of changes in population size and range extent. Overall, we call for collaborative efforts to initiate the urgently needed monitoring of intraspecific variation.

Phylogeography of the Koala, (Phascolarctos cinereus), and Harmonising Data to Inform Conservation.

The Australian continent exhibits complex biogeographic patterns but studies of the impacts of Pleistocene climatic oscillation on the mesic environments of the Southern Hemisphere are limited. The koala (Phascolarctos cinereus), one of Australia's most iconic species, was historically widely distributed throughout much of eastern Australia but currently represents a complex conservation challenge. To better understand the challenges to koala genetic health, we assessed the phylogeographic history of the koala. Variation in the maternally inherited mitochondrial DNA (mtDNA) Control Region (CR) was examined in 662 koalas sampled throughout their distribution. In addition, koala CR haplotypes accessioned to Genbank were evaluated and consolidated. A total of 53 unique CR haplotypes have been isolated from koalas to date (including 15 haplotypes novel to this study). The relationships among koala CR haplotypes were indicative of a single Evolutionary Significant Unit and do not support the recognition of subspecies, but were separated into four weakly differentiated lineages which correspond to three geographic clusters: a central lineage, a southern lineage and two northern lineages co-occurring north of Brisbane. The three geographic clusters were separated by known Pleistocene biogeographic barriers: the Brisbane River Valley and Clarence River Valley, although there was evidence of mixing amongst clusters. While there is evidence for historical connectivity, current koala populations exhibit greater structure, suggesting habitat fragmentation may have restricted female-mediated gene flow. Since mtDNA data informs conservation planning, we provide a summary of existing CR haplotypes, standardise nomenclature and make recommendations for future studies to harmonise existing datasets. This holistic approach is critical to ensuring management is effective and small scale local population studies can be integrated into a wider species context.

Genetic and ecological outcomes of Inga vera subsp. affinis (Leguminosae) tree plantations in a fragmented tropical landscape.

Planting of native trees for habitat restoration is a widespread practice, but the consequences for the retention and transmission of genetic diversity in planted and natural populations are unclear. Using Inga vera subsp. affinis as a model species, we genotyped five natural and five planted populations in the Atlantic forest of northeastern Brazil at polymorphic microsatellite loci. We studied the breeding system and population structure to test how much genetic diversity is retained in planted relative to natural populations. We then genotyped seedlings from these populations to test whether genetic diversity in planted populations is restored by outcrossing to natural populations of I. vera. The breeding system of natural I. vera populations was confirmed to be highly outcrossing (t = 0.92; FIS = -0.061, P = 0.04), with populations showing weak population substructure (FST = 0.028). Genetic diversity in planted populations was 50% less than that of natural populations (planted: AR = 14.9, HO = 0.865 and natural: AR = 30.8, HO = 0.655). However, seedlings from planted populations showed a 30% higher allelic richness relative to their parents (seedlings AR = 10.5, parents AR = 7.6). Understanding the processes and interactions that shape this system are necessary to provide ecologically sensible goals and successfully restore hyper-fragmented habitats. Future restoration plans for I. vera must consider the genetic diversity of planted populations and the potential for gene flow between natural populations in the landscape, in order to preserve ecological interactions (i.e. pollination), and promote opportunities for outcrossing.

Landscape discontinuities influence gene flow and genetic structure in a large, vagile Australian mammal, Macropus fuliginosus.

Large vagile mammals typically exhibit little genetic structuring across their range, particularly when their habitat is essentially continuous. We investigated the population genetic structure of a large vagile Australian macropodid, Macropus fuliginosus, which is continuously distributed across most of southern Australia, using nine highly polymorphic nuclear microsatellite loci. Five distinct genetic units were identified across the range, four on the mainland and one on Kangaroo Island. In addition to the predicted historic Nullarbor Plain Barrier, two unexpected mainland barriers to gene flow were identified. Both were associated with landscape discontinuities (Swan River, Flinders Ranges), which appear within the dispersal capabilities of M. fuliginosus. Typical of large vagile mammals, M. fuliginosus displays high genetic diversity (with the exception of an insular population) and weak genetic structuring (within genetic units). However, the expansion of M. fuliginosus from southwestern Australia during the Pleistocene has resulted in significantly reduced genetic diversity in eastern populations. No significant sex-biased dispersal was detected, although differences in habitat, densities and climatic conditions between the eastern and western regions of the range appear to influence dispersal with the effects of isolation by distance only evident in the west. These results suggest that the biogeography of southern Australia is more complex than previously thought and reveal that seemingly minor landscape features can significantly impact genetic structuring in large vagile mammals.