Genome-wide SNP loci reveal novel insights into koala (Phascolarctos cinereus) population variability across its range

Using bioinformatics to investigate functional diversity: a case study of MHC diversity in koalas

Conservation genomics can greatly improve conservation outcomes of threatened populations, including those impacted by disease. Understanding diversity within immune gene families, including the major histocompatibility complex (MHC) and toll-like receptors (TLR), is important due to the role they play in disease resilience and susceptibility. With recent advancements in sequencing technologies and bioinformatic tools, the cost of generating high-quality sequence data has significantly decreased and made it possible to investigate diversity across entire gene families in large numbers of individuals compared to investigating only a few genes or a few populations previously. Here, we use the koala as a case study for investigating functional diversity across populations. We utilised previous target enrichment data and 438 whole genomes to firstly, determine the level of sequencing depth required to investigate MHC diversity and, secondly, determine the current level of diversity in MHC genes in koala populations. We determined for low complexity, conserved genes such as TLR genes 10 × sequencing depth is sufficient to reliably genotype more than 90% of variants, whereas for complex genes such as the MHC greater than 20 × and preferably 30 × sequencing depth is required. We used whole genome data to identify 270 biallelic SNPs across 24 MHC genes as well as copy number variation (CNV) within class I and class II genes and conduct supertype analysis. Overall, we have provided a bioinformatic workflow for investigating variation in a complex immune gene family from whole genome sequencing data and determined current levels of diversity within koala MHC genes.

The Utility of the Koala Scat: A Scoping Review

The use of samples or scats to provide important ecological, genetic, disease and physiology details on free-range populations is gaining popularity as an alternative non-invasive methodology. Koala populations in SE Queensland and NSW have recently been listed as endangered and continue to face anthropomorphic and stochastic environmental impacts that could potentially lead to their extinction. This scoping review examines the current and potential utility of the koala scat to contribute data relevant to the assessment of koala conservation status and decision making. Although we demonstrate that there is great potential for this methodology in providing details for both individual wild animal and population biology (distribution, abundance, sex ratio, immigration/emigration, genetic diversity, evolutionary significant unit, disease epidemiology, nutrition, reproductive status and stress physiology), the calibre of this information is likely to be a function of the quality of the scat that is sampled.

Small-Mammal Genomics Highlights Viaducts as Potential Dispersal Conduits for Fragmented Populations

Wildlife crossings are implemented in many countries to facilitate the dispersal of animals among habitats fragmented by roads. However, the efficacy of different types of habitat corridors remains poorly understood. We used a comprehensive sampling regime in two lowland dipterocarp forest areas in peninsular Malaysia to sample pairs of small mammal individuals in three treatment types: (1) viaduct sites, at which sampling locations were separated by a highway but connected by a vegetated viaduct; (2) non-viaduct sites, at which sampling locations were separated by a highway and not connected by a viaduct; and (3) control sites, at which there was no highway fragmenting the forest. For four small mammal species, the common tree shrew Tupaia glis, Rajah's spiny rat Maxomys rajah, Whitehead's spiny rat Maxomys whiteheadi and dark-tailed tree rat Niviventer cremoriventer, we used genome-wide markers to assess genetic diversity, gene flow and genetic structure. The differences in genetic distance across sampling settings among the four species indicate that they respond differently to the presence of highways and viaducts. Viaducts connecting forests separated by highways appear to maintain higher population connectivity than forest fragments without viaducts, at least in M. whiteheadi, but apparently not in the other species.

Contrasting whole-genome and reduced representation sequencing for population demographic and adaptive inference: an alpine mammal case study

Genomes capture the adaptive and demographic history of a species, but the choice of sequencing strategy and sample size can impact such inferences. We compared whole genome and reduced representation sequencing approaches to study the population demographic and adaptive signals of the North American mountain goat (Oreamnos americanus). We applied the restriction site-associated DNA sequencing (RADseq) approach to 254 individuals and whole genome resequencing (WGS) approach to 35 individuals across the species range at mid-level coverage (9X) and to 5 individuals at high coverage (30X). We used ANGSD to estimate the genotype likelihoods and estimated the effective population size (Ne), population structure, and explicitly modelled the demographic history with δaδi and MSMC2. The data sets were overall concordant in supporting a glacial induced vicariance and extremely low Ne in mountain goats. We evaluated a set of climatic variables and geographic location as predictors of genetic diversity using redundancy analysis. A moderate proportion of total variance (36% for WGS and 21% for RADseq data sets) was explained by geography and climate variables; both data sets support a large impact of drift and some degree of local adaptation. The empirical similarities of WGS and RADseq presented herein reassuringly suggest that both approaches will recover large demographic and adaptive signals in a population; however, WGS offers several advantages over RADseq, such as inferring adaptive processes and calculating runs-of-homozygosity estimates. Considering the predicted climate-induced changes in alpine environments and the genetically depauperate mountain goat, the long-term adaptive capabilities of this enigmatic species are questionable.

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.

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.

Genome-wide SNP analysis of three moose subspecies at the southern range limit in the contiguous United States

Genome-wide evaluations of genetic diversity and population structure are important for informing management and conservation of trailing-edge populations. North American moose (Alces alces) are declining along portions of the southern edge of their range due to disease, species interactions, and marginal habitat, all of which may be exacerbated by climate change. We employed a genotyping by sequencing (GBS) approach in an effort to collect baseline information on the genetic variation of moose inhabiting the species’ southern range periphery in the contiguous United States. We identified 1920 single nucleotide polymorphisms (SNPs) from 155 moose representing three subspecies from five states: A. a. americana (New Hampshire), A. a. andersoni (Minnesota), and A. a. shirasi (Idaho, Montana, and Wyoming). Molecular analyses supported three geographically isolated clusters, congruent with currently recognized subspecies. Additionally, while moderately low genetic diversity was observed, there was little evidence of inbreeding. Results also indicated > 20% shared ancestry proportions between A. a. shirasi samples from northern Montana and A. a. andersoni samples from Minnesota, indicating a putative hybrid zone warranting further investigation. GBS has proven to be a simple and effective method for genome-wide SNP discovery in moose and provides robust data for informing herd management and conservation priorities. With increasing disease, predation, and climate related pressure on range edge moose populations in the United States, the use of SNP data to identify gene flow between subspecies may prove a powerful tool for moose management and recovery, particularly if hybrid moose are more able to adapt.

Estimating koala density from incidental koala sightings in South-East Queensland, Australia (1997-2013), using a self-exciting spatio-temporal point process model

The koala, Phascolarctos cinereus, is an iconic Australian wildlife species facing a rapid decline in South-East Queensland (SEQLD). For conservation planning, the ability to estimate the size of koala populations is crucial. Systematic surveys are the most common approach to estimate koala populations but because of their cost they are often restricted to small geographic areas and are conducted infrequently. Public interest and participation in the collection of koala sighting data is increasing in popularity, but such data are generally not used for population estimation. We modeled monthly sightings of koalas reported by members of the public from 1997 to 2013 in SEQLD by developing a self-exciting spatio-temporal point process model. This allowed us to account for characteristics that are associated with koala presence (which vary over both space and time) while accounting for detection bias in the koala sighting process and addressing spatial clustering of observations. The density of koalas varied spatially due to the heterogeneous nature of koala habitat in SEQLD, with a mean density of 0.0019 koalas per km2 over the study period. The percentage of land areas with very low densities (0-0.0005 koalas per km2) remained similar throughout the study period representing, on average, 66% of the total study area. The approach described in this paper provides a useful starting point to allow greater use to be made of incidental koala sighting data. We propose that the model presented here could be used to combine systematic koala survey data (which is spatially restricted, but more precise) with koala sighting data (which is incidental and often biased by nature, but often collected over large geographical areas). Our approach could also be adopted for modeling the density of other wildlife species where data is collected in the same manner.

Koala Retrovirus in Northern Australia Shows a Mixture of Stable Endogenization and Exogenous Lineage Diversification within Fragmented Koala Populations

The koala population in northern Australia has become increasingly fragmented due to natural and man-made barriers and interventions. This situation has created a unique opportunity to study both endogenous and exogenous koala retrovirus (KoRV). To determine the impact that population isolation has had on KoRV diversity in Queensland, 272 koalas from six fragmented koala populations were profiled for their KoRV provirus across two natural biogeographical barriers (the St Lawrence Gap and the Brisbane Valley Barrier), one man-made geographical barrier (the city of Brisbane) and two translocation events (the single movement of koalas to an island and the repeated movement of koalas into a koala sanctuary). Analysis revealed that all koalas tested were KoRV-A positive, with 90 - 96% of the detected KoRV provirus from each koala representing a single, likely endogenous, KoRV-A strain. The next most abundant proviral sequence was a defective variant of the dominant KoRV-A strain, accounting for 3 - 10% of detected provirus. The remaining KoRV provirus represented expected exogenous strains of KoRV and included geographically localized patterns of KoRV-B, -C, -D, -F, -G, and -I. These results indicate that lineage diversification of exogenous KoRV is actively ongoing. In addition, comparison of KoRV provirus within known dam-sire-joey family groups from the koala sanctuary revealed that joeys consistently had KoRV proviral patterns more similar to their dams than their sires in KoRV-B, -C and -D provirus composition. Collectively, this study highlights both the consistency of endogenous KoRV and the diversity of exogenous KoRV across the fragmented koala populations in northern Australia.IMPORTANCE KoRV infection has become a permanent part of koalas in northern Australia. With KoRV presence and abundance linked to more severe chlamydial disease and neoplasia in these koalas, understanding how KoRV exists throughout an increasingly fragmented koala population is a key first step in designing conservation and management strategies. This survey of KoRV provirus in Queensland koalas indicates that endogenous KoRV provirus is ubiquitous and consistent throughout the state while exogenous KoRV provirus is diverse and distinct in fragmented koala populations. Understanding the prevalence and impact of both endogenous and exogenous KoRV will be needed to ensure a future for all koala populations.

Double-digest RAD-sequencing: do pre- and post-sequencing protocol parameters impact biological results?

Next-generation sequencing technologies have opened a new era of research in population genetics. Following these new sequencing opportunities, the use of restriction enzyme-based genotyping techniques, such as restriction site-associated DNA sequencing (RAD-seq) or double-digest RAD-sequencing (ddRAD-seq), has dramatically increased in the last decade. From DNA sampling to SNP calling, the laboratory and bioinformatic parameters of enzyme-based techniques have been investigated in the literature. However, the impact of those parameters on downstream analyses and biological results remains less documented. In this study, we investigated the effects of sevral pre- and post-sequencing settings on ddRAD-seq results for two biological systems: a complex of butterfly species (Coenonympha sp.) and several populations of common beech (Fagus sylvatica). Our results suggest that pre-sequencing parameters (i.e., DNA quantity, number of PCR cycles during library preparation) have a significant impact on the number of recovered reads and SNPs, on the number of unique alleles and on individual heterozygosity. In the same way, we found that post-sequencing settings (i.e., clustering and minimum coverage thresholds) influenced loci reconstruction (e.g., number of loci, mean coverage) and SNP calling (e.g., number of SNPs; heterozygosity) but had only a marginal impact on downstream analyses (e.g., measure of genetic differentiation, estimation of individual admixture, and demographic inferences). In addition, replication analyses confirmed the reproducibility of the ddRAD-seq procedure. Overall, this study assesses the degree of sensitivity of ddRAD-seq data to pre- and post-sequencing protocols, and illustrates its robustness when studying population genetics.

Fine-scale barriers to connectivity across a fragmented South-East Asian landscape in six songbird species

Habitat fragmentation is a major extinction driver. Despite dramatically increasing fragmentation across the globe, its specific impacts on population connectivity across species with differing life histories remain difficult to characterize, let alone quantify. Here, we investigate patterns of population connectivity in six songbird species from Singapore, a highly fragmented tropical rainforest island. Using massive panels of genome-wide single nucleotide polymorphisms across dozens of samples per species, we examined population genetic diversity, inbreeding, gene flow and connectivity among species along a spectrum of ecological specificities. We found a higher resilience to habitat fragmentation in edge-tolerant and forest-canopy species as compared to forest-dependent understorey insectivores. The latter exhibited levels of genetic diversity up to three times lower in Singapore than in populations from contiguous forest elsewhere. Using dense genomic and geographic sampling, we identified individual barriers such as reservoirs that effectively minimize gene flow in sensitive understorey birds, revealing that terrestrial forest species may exhibit levels of sensitivity to fragmentation far greater than previously expected. This study provides a blueprint for conservation genomics at small scales with a view to identifying preferred locations for habitat corridors, flagging candidate populations for restocking with translocated individuals and improving the design of future reserves.

Pathological Findings in Koala Retrovirus-positive Koalas (Phascolarctos cinereus) from Northern and Southern Australia

Koala retrovirus (KoRV) infection shows differences in prevalence and load between northern and southern Australian koala populations; however, the effect of this on diseases such as lymphoma and chlamydial disease is unclear. This study compared clinicopathological findings, haematology and splenic lymphoid area of KoRV-positive koalas from northern (Queensland [Qld], n = 67) and southern (South Australia [SA], n = 92) populations in order to provide further insight into KoRV pathogenesis. Blood was collected for routine haematology and for measurement of KoRV proviral load by quantitative polymerase chain reaction (qPCR). Plasma samples were assessed for KoRV viral load by reverse transcriptase qPCR and conjunctival and cloacal swabs were collected for measurement of the load of Chlamydia pecorum (qPCR). During necropsy examination, spleen was collected for lymphoid area analysis. Lymphoma was morphologically similar between the populations and occurred in koalas with the highest KoRV proviral and viral loads. Severe ocular chlamydial disease was observed in both populations, but urinary tract disease was more severe in Qld, despite similar C. pecorum loads. No associations between KoRV and chlamydial disease severity or load were observed, except in SA where viral load correlated positively with chlamydial disease severity. In both populations, proviral and viral loads correlated positively with lymphocyte and metarubricyte counts and correlated negatively with erythrocyte and neutrophil counts. Splenic lymphoid area was correlated positively with viral load. This study has shown further evidence for KoRV-induced oncogenesis and highlighted that lymphocytes and splenic lymphoid tissue may be key sites for KoRV replication. However, KoRV infection appears to be highly complex and continued investigation is required to fully understand its pathogenesis.

Population differentiation and historical demography of the threatened snowy plover Charadrius nivosus (Cassin, 1858)

Delineating conservation units is a complex and often controversial process that is particularly challenging for highly vagile species. Here, we reassess population genetic structure and identify those populations of highest conservation value in the threatened snowy plover (Charadrius nivosus, Cassin, 1858), a partial migrant shorebird endemic to the Americas. We use four categories of genetic data—mitochondrial DNA (mtDNA), microsatellites, Z-linked and autosomal single nucleotide polymorphisms (SNPs)—to: (1) assess subspecies delineation and examine population structure (2) compare the sensitivity of the different types of genetic data to detect spatial genetic patterns, and (3) reconstruct demographic history of the populations analysed. Delineation of two traditionally recognised subspecies was broadly supported by all data. In addition, microsatellite and SNPs but not mtDNA supported the recognition of Caribbean snowy plovers (C. n. tenuirostris) and Floridian populations (eastern C. n. nivosus) as distinct genetic lineage and deme, respectively. Low migration rates estimated from autosomal SNPs (m < 0.03) reflect a general paucity of exchange between genetic lineages. In contrast, we detected strong unidirectional migration (m = 0.26) from the western into the eastern nivosus deme. Within western nivosus, we found no genetic differentiation between coastal Pacific and inland populations. The correlation between geographic and genetic distances was weak but significant for all genetic data sets. All demes showed signatures of bottlenecks occurring during the past 1000 years. We conclude that at least four snowy plover conservation units are warranted: in addition to subspecies nivosus and occidentalis, a third unit comprises the Caribbean tenuirostris lineage and a fourth unit the distinct eastern nivosus deme.

Haematological reference intervals of wild southern Australian koalas (Phascolarctos cinereus)

OBJECTIVE

Current haematology reference intervals (RIs) for koalas were developed in northern Australian koalas, using low numbers and/or individuals of unknown Chlamydia pecorum and koala retrovirus (KoRV) status. This study developed haematological RIs for wild, clinically healthy southern Australian koalas of known C. pecorum and KoRV infection status and investigated the effects of population, age and sex.

METHODS

Haematological RIs were determined for 138 clinically healthy South Australian koalas (Mount Lofty Ranges [MLR], n = 68; Kangaroo Island, n = 70) examined in April 2016 and February 2017, respectively. C. pecorum and KoRV status were determined by PCR.

RESULTS

RIs for southern koala haematological parameters were established for all koalas based on the finding that there were limited differences in haematological values in koalas with subclinical C. pecorum or KoRV infections (P > 0.05), except KoRV-infected koalas had a lower haematocrit than noninfected koalas. MLR koalas had significantly lower erythrocyte mass and leucocyte counts than Kangaroo Island koalas. Young koalas had significantly lower haemoglobin, haematocrit and higher mean cellular haemoglobin concentration and lymphocyte counts than adult koalas. MLR male koalas had elevated erythrocyte, leucocyte and neutrophil counts compared with MLR females.

CONCLUSION

The haematological RIs developed in this study are based on a large number of clinically healthy koalas, where subclinical C. pecorum and KoRV infections had no effect on haematological values and will be a valuable tool during clinical examination and disease investigation by veterinarians and researchers Australia-wide.

Koala retrovirus viral load and disease burden in distinct northern and southern koala populations

Koala retrovirus (KoRV) displays features of both an endogenous and exogenous virus and is linked to neoplasia and immunosuppression in koalas. This study explores the apparent differences in the nature and impact of KoRV infection between geographically and genetically separated "northern" and "southern" koala populations, by investigating the disease status, completeness of the KoRV genome and the proviral (DNA) and viral (RNA) loads of 71 northern and 97 southern koalas. All northern animals were positive for all KoRV genes (gag, pro-pol and env) in both DNA and RNA forms, whereas many southern animals were missing one or more KoRV genes. There was a significant relationship between the completeness of the KoRV genome and clinical status in this population. The proviral and viral loads of the northern population were significantly higher than those of the southern population (P < 0.0001), and many provirus-positive southern animals failed to express any detectable KoRV RNA. Across both populations there was a positive association between proviral load and neoplasia (P = 0.009). Potential reasons for the differences in the nature of KoRV infection between the two populations are discussed.

Population Genomic Signatures of Genetic Structure and Environmental Selection in the Catadromous Roughskin Sculpin Trachidermus fasciatus

Understanding the patterns of genetic diversity and adaptation across species' range is crucial to assess its long-term persistence and determine appropriate conservation measures. The impacts of human activities on the genetic diversity and genetic adaptation to heterogeneous environments remain poorly understood in the marine realm. The roughskin sculpin (Trachidermus fasciatus) is a small catadromous fish, and has been listed as a second-class state protected aquatic animal since 1988 in China. To elucidate the underlying mechanism of population genetic structuring and genetic adaptations to local environments, RAD tags were sequenced for 202 individuals in nine populations across the range of T. fasciatus in China. The pairwise FST values over 9,271 filtered SNPs were significant except that between Dongying and Weifang. All the genetic clustering analysis revealed significant population structure with high support for eight distinct genetic clusters. Both the minor allele frequency spectra and Ne estimations suggested extremely small Ne in some populations (e.g., Qinhuangdao, Rongcheng, Wendeng, and Qingdao), which might result from recent population bottleneck. The strong genetic structure can be partly attributed to genetic drift and habitat fragmentation, likely due to the anthropogenic activities. Annotations of candidate adaptive loci suggested that genes involved in metabolism, development, and osmoregulation were critical for adaptation to spatially heterogenous environment of local populations. In the context of anthropogenic activities and environmental change, results of the present population genomic work provided important contributions to the understanding of genetic differentiation and adaptation to changing environments.

Intraspecific DNA contamination distorts subtle population structure in a marine fish: Decontamination of herring samples before restriction-site associated sequencing and its effects on population genetic statistics

Wild specimens are often collected in challenging field conditions, where samples may be contaminated with the DNA of conspecific individuals. This contamination can result in false genotype calls, which are difficult to detect, but may also cause inaccurate estimates of heterozygosity, allele frequencies and genetic differentiation. Marine broadcast spawners are especially problematic, because population genetic differentiation is low and samples are often collected in bulk and sometimes from active spawning aggregations. Here, we used contaminated and clean Pacific herring (Clupea pallasi) samples to test (a) the efficacy of bleach decontamination, (b) the effect of decontamination on RAD genotypes and (c) the consequences of contaminated samples on population genetic analyses. We collected fin tissue samples from actively spawning (and thus contaminated) wild herring and nonspawning (uncontaminated) herring. Samples were soaked for 10 min in bleach or left untreated, and extracted DNA was used to prepare DNA libraries using a restriction site-associated DNA (RAD) approach. Our results demonstrate that intraspecific DNA contamination affects patterns of individual and population variability, causes an excess of heterozygotes and biases estimates of population structure. Bleach decontamination was effective at removing intraspecific DNA contamination and compatible with RAD sequencing, producing high-quality sequences, reproducible genotypes and low levels of missing data. Although sperm contamination may be specific to broadcast spawners, intraspecific contamination of samples may be common and difficult to detect from high-throughput sequencing data and can impact downstream analyses.

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

The Australian koala is an iconic marsupial with highly specific dietary requirements distributed across heterogeneous environments, over a large geographic range. The distribution and genetic structure of koala populations has been heavily influenced by human actions, specifically habitat modification, hunting and translocation of koalas. There is currently limited information on population diversity and gene flow at a species-wide scale, or with consideration to the potential impacts of local adaptation. Using species-wide sampling across heterogeneous environments, and high-density genome-wide markers (SNPs and PAVs), we show that most koala populations display levels of diversity comparable to other outbred species, except for those populations impacted by population reductions. Genetic clustering analysis and phylogenetic reconstruction reveals a lack of support for current taxonomic classification of three koala subspecies, with only a single evolutionary significant unit supported. Furthermore, ~70% of genetic variance is accounted for at the individual level. The Sydney Basin region is highlighted as a unique reservoir of genetic diversity, having higher diversity levels (i.e., Blue Mountains region; AvHecorr=0.20, PL% = 68.6). Broad-scale population differentiation is primarily driven by an isolation by distance genetic structure model (49% of genetic variance), with clinal local adaptation corresponding to habitat bioregions. Signatures of selection were detected between bioregions, with no single region returning evidence of strong selection. The results of this study show that although the koala is widely considered to be a dietary-specialist species, this apparent specialisation has not limited the koala's ability to maintain gene flow and adapt across divergent environments as long as the required food source is available.

Genome-wide SNP analyses reveal high gene flow and signatures of local adaptation among the scalloped spiny lobster (Panulirus homarus) along the Omani coastline

BACKGROUND

The scalloped spiny lobster (Panulirus homarus) is a popular seafood commodity worldwide and an important export item from Oman. Annual catches in commercial fisheries are in serious decline, which has resulted in calls for the development of an integrated stock management approach. In Oman, the scalloped spiny lobster is currently treated as a single management unit (MU) or stock and there is an absence of information on the genetic population structure of the species that can inform management decisions, particularly at a fine-scale level. This work is the first to identify genome-wide single nucleotide polymorphisms (SNPs) for P. homarus using Diversity Arrays Technology sequencing (DArT-seq) and to elucidate any stock structure in the species.

RESULTS

After stringent filtering, 7988 high utility SNPs were discovered and used to assess the genetic diversity, connectivity and structure of P. homarus populations from Al Ashkharah, Masirah Island, Duqm, Ras Madrakah, Haitam, Ashuwaymiyah, Mirbat and Dhalkut landing sites. Pairwise FST estimates revealed low differentiation among populations (pairwise FST range = - 0.0008 - 0.0021). Analysis of genetic variation using putatively directional FST outliers (504 SNPs) revealed higher and significant pairwise differentiation (p < 0.01) for all locations, with Ashuwaymiyah being the most diverged population (Ashuwaymiyah pairwise FST range = 0.0288-0.0736). Analysis of population structure using Discriminant Analysis of Principal Components (DAPC) revealed a broad admixture among P. homarus, however, Ashuwaymiyah stock appeared to be potentially under local adaptive pressures. Fine scale analysis using Netview R provided further support for the general admixture of P. homarus.

CONCLUSIONS

Findings here suggested that stocks of P. homarus along the Omani coastline are admixed. Yet, fishery managers need to treat the lobster stock from Ashuwaymiyah with caution as it might be subject to local adaptive pressures. We emphasize further study with larger number of samples to confirm the genetic status of the Ashuwaymiyah stock. The approach utilised in this study has high transferability in conservation and management of other marine stocks with similar biological and ecological attributes.

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.

Polygyny does not explain the superior competitive ability of dominant ant associates in the African ant-plant, Acacia (Vachellia) drepanolobium

The Acacia drepanolobium (also known as Vachellia drepanolobium) ant-plant symbiosis is considered a classic case of species coexistence, in which four species of tree-defending ants compete for nesting space in a single host tree species. Coexistence in this system has been explained by trade-offs in the ability of the ant associates to compete with each other for occupied trees versus the ability to colonize unoccupied trees. We seek to understand the proximal reasons for how and why the ant species vary in competitive or colonizing abilities, which are largely unknown. In this study, we use RADseq-derived SNPs to identify relatedness of workers in colonies to test the hypothesis that competitively dominant ants reach large colony sizes due to polygyny, that is, the presence of multiple egg-laying queens in a single colony. We find that variation in polygyny is not associated with competitive ability; in fact, the most dominant species, unexpectedly, showed little evidence of polygyny. We also use these markers to investigate variation in mating behavior among the ant species and find that different species vary in the number of males fathering the offspring of each colony. Finally, we show that the nature of polygyny varies between the two commonly polygynous species, Crematogaster mimosae and Tetraponera penzigi: in C. mimosae, queens in the same colony are often related, while this is not the case for T. penzigi. These results shed light on factors influencing the evolution of species coexistence in an ant-plant mutualism, as well as demonstrating the effectiveness of RADseq-derived SNPs for parentage analysis.

Single-nucleotide polymorphism discovery and panel characterization in the African forest elephant

The continuing decline in forest elephant (Loxodonta cyclotis) numbers due to poaching and habitat reduction is driving the search for new tools to inform management and conservation. For dense rainforest species, basic ecological data on populations and threats can be challenging and expensive to collect, impeding conservation action in the field. As such, genetic monitoring is being increasingly implemented to complement or replace more burdensome field techniques. Single-nucleotide polymorphisms (SNPs) are particularly cost-effective and informative markers that can be used for a range of practical applications, including population census, assessment of human impact on social and genetic structure, and investigation of the illegal wildlife trade. SNP resources for elephants are scarce, but next-generation sequencing provides the opportunity for rapid, inexpensive generation of SNP markers in nonmodel species. Here, we sourced forest elephant DNA from 23 samples collected from 10 locations within Gabon, Central Africa, and applied double-digest restriction-site-associated DNA (ddRAD) sequencing to discover 31,851 tags containing SNPs that were reduced to a set of 1,365 high-quality candidate SNP markers. A subset of 115 candidate SNPs was then selected for assay design and validation using 56 additional samples. Genotyping resulted in a high conversion rate (93%) and a low per allele error rate (0.07%). This study provides the first panel of 107 validated SNP markers for forest elephants. This resource presents great potential for new genetic tools to produce reliable data and underpin a step-change in conservation policies for this elusive species.

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.

Koalas (Phascolarctos cinereus) From Queensland Are Genetically Distinct From 2 Populations in Victoria

The koala (Phascolarctos cinereus) suffered population declines and local extirpation due to hunting in the early 20th century, especially in southern Australia. Koalas were subsequently reintroduced to the Brisbane Ranges (BR) and Stony Rises (SR) by translocating individuals from a population on French Island descended from a small number of founders. To examine genetic diversity and north-south differentiation, we genotyped 13 microsatellite markers in 46 wild koalas from the BR and SR, and 27 Queensland koalas kept at the US zoos. The Queensland koalas displayed much higher heterozygosity (H _O = 0.73) than the 2 southern Australian koala populations examined: H _O = 0.49 in the BR, whereas H _O = 0.41 in the SR. This is consistent with the historical accounts of bottlenecks and founder events affecting the southern populations and contrasts with reports of high genetic diversity in some southern populations. The 2 southern Australian koala populations were genetically similar (F _ST = 0.018, P = 0.052). By contrast, northern and southern Australian koalas were highly differentiated (F _ST = 0.27, P < 0.001), thereby suggesting that geographic structuring should be considered in the conservation management of koalas. Sequencing of 648bp of the mtDNA control region in Queensland koalas found 8 distinct haplotypes, one of which had not been previously detected among koalas. Queensland koalas displayed high mitochondrial haplotype diversity (H = 0.753) and nucleotide diversity (π = 0.0072), indicating along with the microsatellite data that North American zoos have maintained high levels of genetic diversity among their Queensland koalas.

Population genetics of the koala (Phascolarctos cinereus) in north-eastern New South Wales and south-eastern Queensland

Habitat loss and fragmentation are key threats to local koala (Phascolarctos cinereus) populations. Broad-scale management is suboptimal for koalas because distribution models are not easily generalised across regions. Therefore, it is imperative that data relevant to local management bodies are available. Genetic data provides important information on gene flow and potential habitat barriers, including anthropogenic disturbances. Little genetic data are available for nationally significant koala populations in north-eastern New South Wales, despite reported declines due to urbanisation and habitat loss. In this study, we develop 14 novel microsatellite loci to investigate koala populations in north-eastern New South Wales (Port Macquarie, Coffs Harbour, Tyagarah, Ballina) and south-eastern Queensland (Coomera). All locations were significantly differentiated (FST = 0.096–0.213; FʹST = 0.282–0.582), and this pattern was not consistent with isolation by distance (R2 = 0.228, P = 0.058). Population assignment clustered the more northern populations (Ballina, Tyagarah and Coomera), suggesting contemporary gene flow among these sites. For all locations, low molecular variation among (16%) rather than within (84%) sites suggests historical connectivity. These results suggest that koala populations in north-eastern New South Wales and south-eastern Queensland are experiencing contemporary impediments to gene flow, and highlight the importance of maintaining habitat connectivity across this region.