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.

Devil facial tumour disease (DFTD) has decimated wild populations of Tasmanian devils (Sarcophilus harrisii) due to its ability to avoid immune detection and pass from host to host by biting. A small number of devils have been observed to spontaneously recover from the disease which is otherwise fatal. We have sequenced the genomes of these rare cases and compared them to the genomes of devils who succumbed to the disease. Genome-wide association, based on this limited sampling, highlighted two key genomic regions potentially associated with ability to survive DFTD. Following targeted genotyping in additional samples, both of these loci remain significantly different between cases and controls, with the PAX3 locus retaining significance at the 0.001 level, though genome-wide significance was not achieved. We propose that PAX3 may be involved in a regulatory pathway that influences the slowing of tumour growth and may allow more time for an immune response to be mounted in animals with regressed tumours. This provides an intriguing hypothesis for further research and could provide a novel route of treatment for this devastating disease.

A heated debate has recently emerged between tiger farmers and conservationists about the potential consequences of lifting the ban on trade in farmed tiger products in China. This debate has caused unfounded speculation about the extent of the potential market for tiger products. To fill this knowledge gap, we surveyed 1880 residents from a total of six Chinese cities to understand Urban Chinese tiger consumption behavior, knowledge of trade issues and attitudes towards tiger conservation. We found that 43% of respondents had consumed some product alleged to contain tiger parts. Within this user-group, 71% said that they preferred wild products over farmed ones. The two predominant products used were tiger bone plasters (38%) and tiger bone wine (6.4%). 88% of respondents knew that it was illegal to buy or sell tiger products, and 93% agreed that a ban in trade of tiger parts was necessary to conserve wild tigers. These results indicate that while Urban Chinese people are generally supportive of tiger conservation, there is a huge residual demand for tiger products that could resurge if the ban on trade in tiger parts is lifted in China. We suspect that the current supply of the market is predominantly met by fakes or substitutes branded as tiger medicines, but not listing tiger as an ingredient. We suggest that the Traditional Chinese Medicine community should consider re-branding these products as bone-healing medicines in order to reduce the residual demand for real tiger parts over the long-term. The lifting of the current ban on trade in farmed tiger parts may cause a surge in demand for wild tiger parts that consumers say are better. Because of the low input costs associated with poaching, wild-sourced parts would consistently undercut the prices of farmed tigers that could easily be laundered on a legal market. We therefore recommend that the Chinese authorities maintain the ban on trade in tiger parts, and work to improve the enforcement of the existing ban.

Conservation breeding programmes have become widespread as natural habitats shrink, and have been historically managed using pedigree data and an assumption that population founders are unrelated. Molecular genotyping is able to determine founder relatedness, but is rarely used. To empirically test the impact of assuming founders to be unrelated, we utilized data from 203 founding individuals and 11 subsequent years of breeding records for the Tasmanian devil (Sarcophilus harrisii) insurance population. We integrated molecular data (N = 119 founders) and detailed trapping information (N = 203 founders) to test how founder relationship assumptions impacted the genetic characteristics of the population over time (N = 942 unique individuals). We developed a method to combine molecular kinship (using the TrioML estimator), year of birth and trapping location, and integrated the resulting empirical kinship estimates into population management software. We tested the effect of using pedigree data only, versus our integrated approach, on population outcomes. Inbreeding coefficients evaluated using the integrated approach were significantly higher than pedigree-only inbreeding coefficients in the first few years following population establishment. A geographic distanceonly approach showed an association between kinship and probability of successful breeding. Our results show the value in using field and/or molecular data combined with pedigree data in conservation breeding programmes to provide new information for managing crucial populations and improving their success. We caution population managers against commencing expensive conservation breeding programmes in the absence of understanding founder relatedness, especially when wild augmentation is a goal. Long-term costs of assuming founders are unrelated include financial costs, reduced productivity and release of potentially highly inbred individuals.

As species extinction rates increase, genomics provides a powerful tool to support intensive management of threatened species. We use the Tasmanian devil (Sarcophilus harrisii) to demonstrate how conservation genomics can be implemented in threatened species management. We conducted whole genome sequencing (WGS) of 25 individuals from the captive breeding programme and reduced-representation sequencing (RRS) of 98 founders of the same programme. A subset of the WGS samples was also sequenced by RRS, allowing us to directly compare genome-wide heterozygosity with estimates from RRS data. We found good congruence in interindividual variation and gene-ontology classifications between the two data sets, indicating that our RRS data reflect the genome well. We also attempted genome-wide association studies with both data sets (regarding breeding success), but the genomic data suffered from small sample size, while the RRS data suffered from lack of precision, highlighting a key trade-off in the design of conservation genomic research. Nevertheless, we identified a number of candidate genes that may be associated with variation in breeding success. Individual heterozygosity, as measured by WGS or RRS, was not associated with breeding success in captivity but was negatively associated with litter sizes of breeding females in the RRS data set. Our findings enable conservation managers to have confidence in RRS data while understanding its limitations, and provide avenues for further investigation into which processes underlie variation in breeding success in captive Tasmanian devils. We caution, however, that deep functional insights using RRS may be impaired by a lack of precision, especially when marker density is low.

The Tasmanian devil (Sarcophilus harrisii) is threatened with extinction due to the spread of devil facial tumour disease. Polymorphisms in immune genes can provide adaptive potential to resist diseases. Previous studies in diversity at immune loci in wild species have almost exclusively focused on genes of the major histocompatibility complex (MHC); however, these genes only account for a fraction of immune gene diversity. Devils lack diversity at functionally important immunity loci, including MHC and Toll-like receptor genes. Whether there are polymorphisms at devil immune genes outside these two families is unknown. Here, we identify polymorphisms in a wide range of key immune genes, and develop assays to type single nucleotide polymorphisms (SNPs) within a subset of these genes. A total of 167 immune genes were examined, including cytokines, chemokines and natural killer cell receptors. Using genome-level data from ten devils, SNPs within coding regions, introns and 10 kb flanking genes of interest were identified. We found low polymorphism across 167 immune genes examined bioinformatically using whole-genome data. From this data, we developed long amplicon assays to target nine genes. These amplicons were sequenced in 29-220 devils and found to contain 78 SNPs, including eight SNPS within exons. Despite the extreme paucity of genetic diversity within these genes, signatures of balancing selection were exhibited by one chemokine gene, suggesting that remaining diversity may hold adaptive potential. The low functional diversity may leave devils highly vulnerable to infectious disease, and therefore, monitoring and preserving remaining diversity will be critical for the long-term management of this species. Examining genetic variation in diverse immune genes should be a priority for threatened wildlife species. This study can act as a model for broad-scale immunogenetic diversity analysis in threatened species.

Background Recent advances in genomics have greatly increased research opportunities for non-model species. For wildlife, a growing availability of reference genomes means that population genetics is no longer restricted to a small set of anonymous loci. When used in conjunction with a reference genome, reduced-representation sequencing (RRS) provides a cost-effective method for obtaining reliable diversity information for population genetics. Many software tools have been developed to process RRS data, though few studies of non-model species incorporate genome alignment in calling loci. A commonly-used RRS analysis pipeline, Stacks, has this capacity and so it is timely to compare its utility with existing software originally designed for alignment and analysis of whole genome sequencing data. Here we examine population genetic inferences from two species for which reference-aligned reduced-representation data have been collected. Our two study species are a threatened Australian marsupial (Tasmanian devil Sarcophilus harrisii ; declining population) and an Arctic-circle migrant bird (pink-footed goose Anser brachyrhynchus ; expanding population). Analyses of these data are compared using Stacks versus two widely-used genomics packages, SAMtools and GATK. We also introduce a custom R script to improve the reliability of single nucleotide polymorphism (SNP) calls in all pipelines and conduct population genetic inferences for non-model species with reference genomes. Results Although we identified orders of magnitude fewer SNPs in our devil dataset than for goose, we found remarkable symmetry between the two species in our assessment of software performance. For both datasets, all three methods were able to delineate population structure, even with varying numbers of loci. For both species, population structure inferences were influenced by the percent of missing data. Conclusions For studies of non-model species with a reference genome, we recommend combining Stacks output with further filtering (as included in our R pipeline) for population genetic studies, paying particular attention to potential impact of missing data thresholds. We recognise SAMtools as a viable alternative for researchers more familiar with this software. We caution against the use of GATK in studies with limited computational resources or time. Electronic supplementary material The online version of this article (10.1186/s12864-019-5806-y) contains supplementary material, which is available to authorized users.

Farewell Spit gannetry in New Zealand. Our results showed intra-and inter-specific variation in the protein, lipid and water composition of prey captured by our sample of 111 Australasian gannets. In addition, we observed significant differences in the Australasian gannets' nutritional niche between seasons. We provide evidence of sex-specific macronutrient foraging strategies in a successful marine predator in the wild. We have shown that in spite of fluctuations in the nutritional composition of foods available to Australasian gannets, males consistently capture prey with higher protein-to-lipid ratios and lower lipid-to-water ratios than females. These results aid to better understand the evolutionary relationship between macronutrient selection and sex-specific traits in wild animals. They also suggest an incentive for these predators to combine individually imbalanced but nutritionally complementary foods to achieve dietary balance, further highlighting the likelihood that prey selection is guided by the balance of macronutrients, rather than energy alone.