Search

Abstract The legume albumin-1 gene family, arising after nodulation, encodes linear a- and b-chain peptides for nutrient storage and defense. Intriguingly, in one prominent legume, Clitoria ternatea, the b-chains are replaced by domains producing ultra-stable cyclic peptides called cyclotides. The mechanism of this gene hijacking is until now unknown. Cyclotides require recruitment of ligase-type asparaginyl endopeptidases (AEPs) for maturation (cyclization), necessitating co-evolution of two gene families. Here we compare a chromosome-level C. ternatea genome with grain legumes to reveal an 8 to 40-fold expansion of the albumin-1 gene family, enabling the additional loci to undergo diversification. Iterative rounds of albumin-1 duplication and diversification create four albumin-1 enriched genomic islands encoding cyclotides, where they are physically grouped by similar pI and net charge values. We identify an ancestral hydrolytic AEP that exhibits neofunctionalization and multiple duplication events to yield two ligase-type AEPs. We propose cyclotides arise by convergence in C. ternatea where their presence enhances defense from biotic attack, thus increasing fitness compared to lineages with linear b-chains and ultimately driving the replacement of b-chains with cyclotides.

Abstract Wildlife harbour a diverse range of microorganisms that affect their health and development. Marsupials are born immunologically naïve and physiologically underdeveloped, with primary development occurring inside a pouch. Secretion of immunological compounds and antimicrobial peptides in the epithelial lining of the female’s pouch, pouch young skin, and through the milk, are thought to boost the neonate’s immune system and potentially alter the pouch skin microbiome. Here, using 16S rRNA amplicon sequencing, we characterised the Tasmanian devil pouch skin microbiome from 25 lactating and 30 non-lactating wild females to describe and compare across these reproductive stages. We found that the lactating pouch skin microbiome had significantly lower amplicon sequence variant richness and diversity than non-lactating pouches, however there was no overall dissimilarity in community structure between lactating and non-lactating pouches. The top five phyla were found to be consistent between both reproductive stages, with over 85% of the microbiome being comprised of Firmicutes, Proteobacteria, Fusobacteriota, Actinobacteriota, and Bacteroidota. The most abundant taxa remained consistent across all taxonomic ranks between lactating and non-lactating pouch types. This suggests that any potential immunological compounds or antimicrobial peptide secretions did not significantly influence the main community members. Of the more than 16,000 total identified amplicon sequence variants, 25 were recognised as differentially abundant between lactating and non-lactating pouches. It is proposed that the secretion of antimicrobial peptides in the pouch act to modulate these microbial communities. This study identifies candidate bacterial clades on which to test the activity of Tasmanian devil antimicrobial peptides and their role in pouch young protection, which in turn may lead to future therapeutic development for human diseases.

Abstract Conservation breeding programs are an effective approach to addressing biodiversity loss. Captive populations are managed to maintain genetic diversity, yet there remains an “implementation gap” in effectively translating molecular genetic data into management. Technological advancements are facilitating rapid generation of genetic data, increasing accessibility for breeding programs. In 2010, Frankham and colleagues proposed a six‐stage process for establishing successful conservation breeding and release programs. Here, we describe the conservation breeding program for the critically endangered Bellinger River turtle (Myuchelys georgesi) and characterize the value of genetic sampling for informing management actions. By generating a chromosome‐level genome and population genetic data, we investigated past and present diversity and assessed relatedness among captive founders. We present a framework modeled on Frankham and colleagues six stages to assist managers in implementing genetic data into actionable conservation strategies. This framework, and worked case study, for managers aims to better guide implementation of genetic approaches into conservation breeding programs.

Abstract* The swift parrot (Lathamus discolor) is a Critically Endangered migratory parrot that breeds in Tasmania and winters on the Australian mainland. Here we provide a reference genome assembly for the swift parrot. We sequence PacBio HiFi reads to create a high-quality reference assembly and identify a complete mitochondrial sequence. We also generate a reference transcriptome from five organs to inform genome annotation. The genome was 1.24 Gb in length and consisted of 847 contigs with a contig N50 of 18.97 Gb and L50 of 20 contigs. This study provides an annotated reference assembly and transcriptomic resources for the swift parrot to assist in future conservation genomic research.

Abstract The number of genome-level resources for non-model species continues to rapidly expand. However, frog species remain underrepresented, with up to 90% of frog genera having no genomic or transcriptomic data. Here, we assemble the first genomic and transcriptomic resources for the recently described southern stuttering frog (Mixophyes australis). The southern stuttering frog is ground-dwelling, inhabiting naturally vegetated riverbanks in south-eastern Australia. Using PacBio HiFi long-read sequencing and Hi-C scaffolding, we generated a high-quality genome assembly, with a scaffold N50 of 369.3 Mb and 95.1% of the genome contained in twelve scaffolds. Using this assembly, we identified the mitochondrial genome, and assembled six tissue-specific transcriptomes. We also bioinformatically characterised novel sequences of two families of antimicrobial peptides (AMPs) in the southern stuttering frog, the cathelicidins and β-defensins. While traditional peptidomic approaches to peptide discovery have typically identified one or two AMPs in a frog species from skin secretions, our bioinformatic approach discovered 12 cathelicidins and two β-defensins that were expressed in a range of tissues. We investigated the novelty of the peptides and found diverse predicted activities. Our bioinformatic approach highlights the benefits of multi-omics resources in peptide discovery and contributes valuable genomic resources in an under-represented taxon.

Abstract The Tasmanian devil (Sarcophilus harrisii) is endangered due to the spread of Devil Facial Tumour Disease (DFTD), a contagious cancer with no current treatment options. Here we test whether seven recently characterized Tasmanian devil cathelicidins are involved in cancer regulation. We measured DFTD cell viability in vitro following incubation with each of the seven peptides and describe the effect of each on gene expression in treated cells. Four cathelicidins (Saha-CATH3, 4, 5 and 6) were toxic to DFTD cells and caused general signs of cellular stress. The most toxic peptide (Saha-CATH5) also suppressed the ERBB and YAP1/TAZ signaling pathways, both of which have been identified as important drivers of cancer proliferation. Three cathelicidins induced inflammatory pathways in DFTD cells that may potentially recruit immune cells in vivo. This study suggests that devil cathelicidins have some anti-cancer and inflammatory functions and should be explored further to determine whether they have potential as treatment leads.

The Kroombit tinkerfrog (Taudactylus pleione) is a stream-dwelling amphibian of the Myobatrachidae family. It is listed as Critically Endangered and is at high risk of extinction due to chytridiomycosis. Here, we provide the first genome assembly of the evolutionarily distinct Taudactylus genus. We sequenced PacBio HiFi reads to assemble a high-quality long-read genome and identified the mitochondrial genome. We also generated a global transcriptome from a tadpole to improve gene annotation. The genome was 5.52 Gb in length and consisted of 4,196 contigs with a contig N50 of 8.853 Mb and an L50 of 153. This study provides the first genomic resources for the Kroombit tinkerfrog to assist in future phylogenetic, environmental DNA, conservation breeding, and disease susceptibility studies.

The Kroombit tinkerfrog (Taudactylus pleione) is a stream-dwelling amphibian of the Myobatrachidae family. It is listed as Critically Endangered and is at high risk of extinction due to chytridiomycosis. Here, we provide the first genome assembly of the evolutionarily distinct Taudactylus genus. We sequenced PacBio HiFi reads to assemble a high-quality long-read genome and identified the mitochondrial genome. We also generated a global transcriptome from a tadpole to improve gene annotation. The genome was 5.52 Gb in length and consisted of 4,196 contigs with a contig N50 of 8.853 Mb and an L50 of 153. This study provides the first genomic resources for the Kroombit tinkerfrog to assist in future phylogenetic, environmental DNA, conservation breeding, and disease susceptibility studies.

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.

Abstract Top carnivores are essential for maintaining ecosystem stability and biodiversity. Yet, carnivores are declining globally and current in situ threat mitigations cannot halt population declines. As such, translocations of carnivores to historic sites or those outside the species’ native range are becoming increasingly common. As carnivores are likely to impact herbivore and small predator populations, understanding how carnivores interact within an ecosystem following translocation is necessary to inform potential remedial management and future translocations. Dietary analyses provide a preliminary assessment of the direct influence of translocated carnivores on a recipient ecosystem. We used a metabarcoding approach to quantify the diet of Tasmanian devils introduced to Maria Island, Tasmania, a site outside the species’ native range. We extracted DNA from 96 scats and used a universal primer set targeting the vertebrate 12S rRNA gene to identify diet items. Tasmanian devils on Maria Island had an eclectic diet, with 63 consumed taxa identified. Cat DNA was detected in 14% of scats, providing the first instance of cats appearing as part of Tasmanian devil diets either via predation or scavenging. Short‐tail shearwaters and little penguins were commonly consumed, corresponding with previous surveys showing sharp population declines in these species since the introduction of Tasmanian devils. Our results indicate that the introduction of carnivores to novel ecosystems can be very successful for the focal species, but that commonly consumed species should be closely monitored to identify any vulnerable species in need of remedial management.

The numbat (Myrmecobius fasciatus) is an endangered Australian marsupial, and the last surviving member of the Myrmecobiidae family. The numbat regularly undergoes torpor and is unique amongst marsupials as it is the only diurnal and termitivorous species. Here we sequenced the first draft genome of the numbat using 10× Genomics Chromium linked-read technology, resulting in a 3.42 Gbp genome with a scaffold N50 of 223 kbp. A global transcriptome from liver, lung and tongue was also generated to aid genome annotation, identifying 21,465 protein-coding genes. To investigate adaptation to the numbat’s termitivorous diet and arid/semi-arid range, we interrogated the most highly expressed transcripts within the tongue and manually annotated taste, vomeronasal and aquaporin gene families. Antimicrobial proteins and proteins involved in digestion were highly expressed in the tongue, alongside umami taste receptors. However, sweet taste receptors were not expressed in this tissue, which combined with the putative contraction of the bitter taste receptor gene repertoire in the numbat genome, may indicate a potential evolutionary adaptation to their specialised termitivorous diet. Vomeronasal and aquaporin gene repertoires were similar to other marsupials. The draft numbat genome is a valuable tool for conservation and can be applied to population genetics/genomics studies and to investigate the unique biology of this interesting species.

Abstract Incorporating mate choice into conservation breeding programs can improve reproduction and the retention of natural behaviors. However, different types of genetic‐based mate choice can have varied consequences for genetic diversity management. As a result, it is important to examine mechanisms of mate choice in captivity to assess its costs and benefits. Most research in this area has focused on experimental pairing trials; however, this resource‐intensive approach is not always feasible in captive settings and can interfere with other management constraints. We used generalized linear mixed models and permutation approaches to investigate overall breeding success in group‐housed Tasmanian devils at three nonmutually exclusive mate choice hypotheses: (a) advantage of heterozygous individuals, (b) advantage of dissimilar mates, and (c) optimum genetic distance, using both 1,948 genome‐wide SNPs and 12 MHC‐linked microsatellites. The managed devil insurance population is the largest such breeding program in Australia and is known to have high variance in reproductive success. We found that nongenetic factors such as age were the best predictors of breeding success in a competitive breeding scenario, with younger females and older males being more successful. We found no evidence of mate choice under the hypotheses tested. Mate choice varies among species and across environments, so we advocate for more studies in realistic captive management contexts as experimental or wild studies may not apply. Conservation managers must weigh up the need to wait for adequate sample sizes to detect mate choice with the risk that genetic changes may occur during this time in captivity. Our study shows that examining and integrating mate choice into the captive management of species housed in realistic, semi‐natural group‐based contexts may be more difficult than previously considered.

Biodiversity is declining globally, and Australia has one of the worst extinction records for mammals. The development of sequencing technologies means that genomic approaches are now available as important tools for wildlife conservation and management. Despite this, genome sequences are available for only 5% of threatened Australian species. Here we report the first reference genome for the woylie (Bettongia penicillata ogilbyi), a critically endangered marsupial from Western Australia, and the first genome within the Potoroidae family. The woylie reference genome was generated using Pacific Biosciences HiFi long-reads, resulting in a 3.39 Gbp assembly with a scaffold N50 of 6.49 Mbp and 86.5% complete mammalian BUSCOs. Assembly of a global transcriptome from pouch skin, tongue, heart and blood RNA-seq reads was used to guide annotation with Fgenesh++, resulting in the annotation of 24,655 genes. The woylie reference genome is a valuable resource for conservation, management and investigations into disease-induced decline of this critically endangered marsupial.

The value of genetic management is increasingly recognised in natural resource management particularly in relation to threatened species conservation and translocations. Technical advancements in methodologies available to conservation managers for genetic management present a challenge for long-running programs that monitor genetic diversity over time. Given that microsatellite and single nucleotide polymorphism (SNP) data are not interchangeable, conservation managers are faced with two options: continue using older techniques for consistency (limits genetic knowledge); or discard previous data and move to more precise modern genetic markers (wastes resources and prevents temporal monitoring). To achieve a balance, we should consider whether conservation recommendations made are broadly similar between interpretations of different marker sets, even if estimates are not directly comparable. Here we generated 13,194 SNPs for the woylie (Bettongia penicillata), an Australian marsupial, to provide conservation recommendations across their current range (N = 13 indigenous, reintroduced and haven populations) and provide guidance to conservation managers on how these recommendations align with those previously made using 12 microsatellite markers. Observed heterozygosity was significantly less than expected heterozygosity reflecting the historical bottleneck for the species. Only 4 of the 13 populations showed evidence of inbreeding, and effective population sizes ranged from 9.9 to 268.9. Trends in observed heterozygosity were broadly similar between microsatellite and SNP marker sets. Herein we provide advice for conservation managers undertaking translocations on using genetic data, highlight the benefits of moving to newer SNP based approaches for monitoring conservation actions over time, and show that historical management actions are supported by new technologies.

Abstract Background Captivity presents extreme lifestyle changes relative to the wild, and evidence of microbiome dysbiosis in captive animals is growing. The gut microbiome plays a crucial role in host health. Whilst captive breeding and subsequent reintroduction to the wild is important for conservation, such efforts often have limited success. Post-release monitoring is essential for assessing translocation success, but changes to the microbiome of released individuals are poorly understood. The Tasmanian devil was previously shown to exhibit loss of microbiome diversity as a result of intense captive management. This current study examines changes in the devil gut microbiome in response to translocation and aims to determine if perturbations from captivity are permanent or reversible. Methods Using 16S rRNA amplicon sequencing, we conducted temporal monitoring of the gut microbiome of released devils during two translocation events, captive-to-wild and wild-to-wild. To investigate whether the microbiome of the released devils changed following translocation, we characterized their microbiome at multiple time points during the translocation process over the course of 6–12 months and compared them to the microbiome of wild incumbent devils (resident wild-born devils at the respective release sites). Results We showed that the pre-release microbiome was significantly different to the microbiome of wild incumbent animals, but that the microbiomes of animals post-release (as early as 3 to 4 weeks post-release) were similar to wild incumbents. The gut microbiome of released animals showed significant compositional shifts toward the wild incumbent microbiome of both translocation events. Conclusion Our results suggest that the devil gut microbiome is dynamic and that loss of microbiome diversity in captivity can be restored following release to the wild. We recommend the broader application of microbiome monitoring in wildlife translocation programs to assess the impacts of translocation on animal microbiomes.

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

Abstract For bottlenecked populations of threatened species, supplementation often leads to improved population metrics (genetic rescue), provided that guidelines can be followed to avoid negative outcomes. In cases where no “ideal” source populations exist, or there are other complicating factors such as prevailing disease, the benefit of supplementation becomes uncertain. Bringing multiple data and analysis types together to plan genetic management activities can help. Here, we consider three populations of Tasmanian devil, Sarcophilus harrisii, as candidates for genetic rescue. Since 1996, devil populations have been severely impacted by devil facial tumour disease (DFTD), causing significant population decline and fragmentation. Like many threatened species, the key threatening process for devils cannot currently be fully mitigated, so species management requires a multifaceted approach. We examined diversity of 31 putatively neutral and 11 MHC‐linked microsatellite loci of three remnant wild devil populations (one sampled at two time‐points), alongside computational diversity projections, parameterized by field data from DFTD‐present and DFTD‐absent sites. Results showed that populations had low diversity, connectivity was poor, and diversity has likely decreased over the last decade. Stochastic simulations projected further diversity losses. For a given population size, the effects of DFTD on population demography (including earlier age at death and increased female productivity) did not impact diversity retention, which was largely driven by final population size. Population sizes ≥500 (depending on the number of founders) were necessary for maintaining diversity in otherwise unmanaged populations, even if DFTD is present. Models indicated that smaller populations could maintain diversity with ongoing immigration. Taken together, our results illustrate how multiple analysis types can be combined to address complex population genetic challenges.

Devastating fires in Australia over 2019-20 decimated native fauna and flora, including koalas. The resulting population bottleneck, combined with significant loss of habitat, increases the vulnerability of remaining koala populations to threats which include disease. Chlamydia is one disease which causes significant morbidity and mortality in koalas. The predominant pathogenic species, Chlamydia pecorum, causes severe ocular, urogenital and reproductive tract disease. In marsupials, including the koala, gene expansions of an antimicrobial peptide family known as cathelicidins have enabled protection of immunologically naïve pouch young during early development. We propose that koala cathelicidins are active against Chlamydia and other bacteria and fungi. Here we describe ten koala cathelicidins, five of which contained full length coding sequences that were widely expressed in tissues throughout the body. Focusing on these five, we investigate their antimicrobial activity against two koala C. pecorum isolates from distinct serovars; MarsBar and IPTaLE, as well as other bacteria and fungi. One cathelicidin, PhciCath5, inactivated C. pecorum IPTaLE and MarsBar elementary bodies and significantly reduced the number of inclusions compared to the control (p

Coronavirus disease 2019 (COVID-19) and other epidemics (such as severe acute respiratory syndrome [SARS], Ebola, and H1N1) are stark reminders that knowledge of animal behavior and ecosystem health are key to controlling the spread of zoonotic diseases early in their onset. However, we have very limited information about the set of behavioral and ecological factors that promote viral spillover and the effects that has on ecosystem health and disease transmission. Thus, expanding our current knowledge of reservoir hosts and pandemics represents an urgent and critical tool in ecological epidemiology. We also propose to create an integrative database that ranks animal species in terms of their likelihood as hosts for specific infectious diseases. We call for a global and cooperative effort of field and laboratory scientists to create, maintain, and update this information in order to reduce the severity of future pandemics.

Antechinus are a genus of mouse-like marsupials that exhibit a rare reproductive strategy known as semelparity and also naturally develop age-related neuropathologies similar to those in humans. We provide the first annotated antechinus reference genome for the brown antechinus (Antechinus stuartii). The reference genome is 3.3 Gb in size with a scaffold N50 of 73Mb and 93.3% complete mammalian BUSCOs. Using bioinformatic methods we assign scaffolds to chromosomes and identify 0.78 Mb of Y-chromosome scaffolds. Comparative genomics revealed interesting expansions in the NMRK2 gene and the protocadherin gamma family, which have previously been associated with aging and age-related dementias respectively. Transcriptome data displayed expression of common Alzheimer’s related genes in the antechinus brain and highlight the potential of utilising the antechinus as a future disease model. The valuable genomic resources provided herein will enable future research to explore the genetic basis of semelparity and age-related processes in the antechinus.

Background Vulnerable species experiencing inbreeding depression are prone to localised extinctions because of their reduced fitness. For Tasmanian devils, the rapid spread of devil facial tumour disease (DFTD) has led to population declines and fragmentation across the species’ range. Here we show that one of the few remaining DFTD-free populations of Tasmanian devils is experiencing inbreeding depression. Moreover, this population has experienced a significant reduction in reproductive success over recent years. Methods We used 32 microsatellite loci to examine changes in genetic diversity and inbreeding in the wild population at Woolnorth, alongside field data on breeding success from females to test for inbreeding depression. Results Wefound that maternal internal relatedness has a negative impact on litter sizes. The results of this study imply that this population may be entering an extinction vortex and that to protect the population genetic rescue should be considered. This study provides conservation managers with useful information for managing wild devils and provides support for the “Wild Devil Recovery Program”, which is currently augmenting small, isolated populations.

The study of the gut microbiome in threatened wildlife species has enormous potential to improve conservation efforts and gain insights into host-microbe coevolution. Threatened species are often housed in captivity, and during this process undergo considerable changes to their gut microbiome. Studying the gut microbiome of captive animals therefore allows identification of dysbiosis and opportunities for improving management practices in captivity and for subsequent translocations. Manipulation of the gut microbiome through methods such as fecal transplant may offer an innovative means of restoring dysbiotic microbiomes in threatened species to provide health benefits. Finally, characterization of the gut microbiome (including the viral components, or virome) provides important baseline health information and may lead to discovery of significant microbial pathogens. Here we summarize our current understanding of microbiomes in Australian marsupial species.

Abstract Early pregnancy is a critical time for successful reproduction; up to half of human pregnancies fail before the development of the definitive chorioallantoic placenta. Unlike the situation in eutherian mammals, marsupial pregnancy is characterised by a long pre-implantation period prior to the development of the short-lived placenta, making them ideal models for study of the uterine environment promoting embryonic survival pre-implantation. Here we present a transcriptomic study of pre-implantation marsupial pregnancy, and identify differentially expressed genes in the Sminthopsis crassicaudata uterus involved in metabolism and biosynthesis, transport, immunity, tissue remodelling, and uterine receptivity. Interestingly, almost one quarter of the top 50 genes that are differentially upregulated in early pregnancy are putatively involved in histotrophy, highlighting the importance of nutrient transport to the conceptus prior to the development of the placenta. This work furthers our understanding of the mechanisms underlying survival of pre-implantation embryos in the earliest live bearing ancestors of mammals.

Abstract Background Retention of genetic diversity and demographic sustainability are the cornerstones of conservation breeding success. In theory, monogamous breeding with equal reproductive output will retain genetic diversity in insurance populations more effectively than group housing which allows mate choice or intrasexual competition. However, the ecological relevance of group housing to a species can outweigh the theoretical benefits of forced monogamy. Here we investigated the influence of different types of captive housing (group (mate choice) versus intensive (forced monogamy)) on reproductive success, litter size and genetic diversity in the endangered Tasmanian devil (Sarcophilus harrisii). Results For male Tasmanian devils, the proportion of individuals that failed to reproduce was significantly greater in group (maximum 10 individuals) than intensive housing. This suggests greater genetic diversity is retained when devils are bred in intensive housing. For male devils, body weight predicted reproductive success in group housing, suggesting certain individuals are dominant due to a larger body size, leading to unequal genetic contribution. We then used simulation models to predict rates of decline in genetic diversity and inbreeding accumulation over time comparing group and intensive housing. When managed independently, empirically observed reproductive outputs were predicted to result in large accumulations in inbreeding and loss of gene diversity in both housing types, although these effects were greater in group housing. Transferring individuals between the housing facilities decreased inbreeding accumulation and increased gene diversity in both housing types highlighting the importance of managing independent zoo populations collectively. Conclusions If conservation programs wish to provide mate choice opportunities through group housing, the impact intrasexual competition will have on dominance and sequential reproductive opportunities needs to be understood prior to commencement. Group housing is becoming increasingly topical as it provides potential ecological benefits, may decrease mate incompatibilities and increase offspring fitness, however it can also result in the loss of genetic diversity in already genetically depauperate species.

Abstract Inbreeding depression occurs when inbred individuals experience reduced fitness as a result of reduced genome-wide heterozygosity. The Tasmanian devil faces extinction due to a contagious cancer, devil facial tumour disease (DFTD). An insurance metapopulation was established in 2006 to ensure the survival of the species and to be used as a source population for re-wilding and genetic rescue. The emergence of DFTD and the rapid decline of wild devil populations have rendered the species at risk of inbreeding depression. We used 33 microsatellite loci to (1) reconstruct a pedigree for the insurance population and (2) estimate genome-wide heterozygosity for 200 individuals. Using heterozygosity-fitness correlations, we investigated the effect of heterozygosity on six diverse fitness measures (ulna length, asymmetry, weight-at-weaning, testes volume, reproductive success and survival). Despite statistically significant evidence of variation in individual inbreeding in this population, we found no associations between inbreeding and any of our six fitness measurements. We propose that the benign environment in captivity may decrease the intensity of inbreeding depression, relative to the stressful conditions in the wild. Future work will need to measure fitness of released animals to facilitate translation of this data to the broader conservation management of the species in its native range.

Abstract Captive breeding of threatened species, for release to the wild, is critical for conservation. This strategy, however, risks producing captive-raised animals with traits poorly suited to the wild. We describe the first study to characterise accumulated consequences of long-term captive breeding on behaviour, by following the release of Tasmanian devils to the wild. We test the impact of prolonged captive breeding on the probability that captive-raised animals are fatally struck by vehicles. Multiple generations of captive breeding increased the probability that individuals were fatally struck, a pattern that could not be explained by other confounding factors (e.g. age or release site). Our results imply that long-term captive breeding programs may produce animals that are naïve to the risks of the post-release environment. Our analyses have already induced changes in management policy of this endangered species, and serve as model of productive synergy between ecological monitoring and conservation strategy.

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

Extinction risk is increasing for a range of species due to a variety of threats, including disease. Emerging infectious diseases can cause severe declines in wild animal populations, increasing population fragmentation and reducing gene flow. Small, isolated, host populations may lose adaptive potential and become more susceptible to extinction due to other threats. Management of the genetic consequences of disease-induced population decline is often necessary. Whilst disease threats need to be addressed, they can be difficult to mitigate. Actions implemented to conserve the Tasmanian devil (Sarcophilus harrisii), which has suffered decline to the deadly devil facial tumour disease (DFTD), exemplify how genetic management can be used to reduce extinction risk in populations threatened by disease. Supplementation is an emerging conservation technique that may benefit populations threatened by disease by enabling gene flow and conserving their adaptive potential through genetic restoration. Other candidate species may benefit from genetic management via supplementation but concerns regarding outbreeding depression may prevent widespread incorporation of this technique into wildlife disease management. However, existing knowledge can be used to identify populations that would benefit from supplementation where risk of outbreeding depression is low. For populations threatened by disease and, in situations where disease eradication is not an option, wildlife managers should consider genetic management to buffer the host species against inbreeding and loss of genetic diversity.

Successful captive breeding programs are crucial to the long-term survival of many threatened species. However, pair incompatibility (breeding failure) limits sustainability of many captive populations. Understanding whether the drivers of this incompatibility are behavioral, genetic, or a combination of both, is crucial to improving breeding programs. We used 28 years of pairing data from the San Diego Zoo koala colony, plus genetic analyses using both major histocompatibility complex (MHC)-linked and non-MHC-linked microsatellite markers, to show that both genetic and non-genetic factors can influence mating success. Male age was reconfirmed to be a contributing factor to the likelihood of a koala pair copulating. This trend could also be related to a pair’s age difference, which was highly correlated with male age in our dataset. Familiarity was reconfirmed to increase the probability of a successful copulation. Our data provided evidence that females select mates based on MHC and genome-wide similarity. Male heterozygosity at MHC class II loci was associated with both pre- and post-copulatory female choice. Genome-wide similarity, and similarity at the MHC class II DAB locus, were also associated with female choice at the post-copulatory level. Finally, certain MHC-linked alleles were associated with either increased or decreased mating success. We predict that utilizing a variety of behavioral and MHC-dependent mate choice mechanisms improves female fitness through increased reproductive success. This study highlights the complexity of mate choice mechanisms in a species, and the importance of ascertaining mate choice mechanisms to improve the success of captive breeding programs.

Devil facial tumor disease (DFTD) is renowned for its successful evasion of the host immune system. Down regulation of the major histocompatabilty complex class I molecule (MHC-I) on the DFTD cells is a primary mechanism of immune escape. Immunization trials on captive Tasmanian devils have previously demonstrated that an immune response against DFTD can be induced, and that immune-mediated tumor regression can occur. However, these trials were limited by their small sample sizes. Here, we describe the results of two DFTD immunization trials on cohorts of devils prior to their wild release as part of the Tasmanian Government’s Wild Devil Recovery project. 95% of the devils developed anti-DFTD antibody responses. Given the relatively large sample sizes of the trials (N = 19 and N = 33), these responses are likely to reflect those of the general devil population. DFTD cells manipulated to express MHC-I were used as the antigenic basis of the immunizations in both trials. Although the adjuvant composition and number of immunizations differed between trials, similar anti-DFTD antibody levels were obtained. The first trial comprised DFTD cells and the adjuvant combination of ISCOMATRIX™, polyIC, and CpG with up to four immunizations given at monthly intervals. This compared to the second trial whereby two immunizations comprising DFTD cells and the adjuvant combination ISCOMATRIX™, polyICLC (Hiltonol®) and imiquimod were given a month apart, providing a shorter and, therefore, more practical protocol. Both trials incorporated a booster immunization given up to 5 months after the primary course. A key finding was that devils in the second trial responded more quickly and maintained their antibody levels for longer compared to devils in the first trial. The different adjuvant combination incorporating the RNAase resistant polyICLC and imiquimod used in the second trial is likely to be responsible. The seroconversion in the majority of devils in these anti-DFTD immunization trials was remarkable, especially as DFTD is hallmarked by its immune evasion mechanisms. Microsatellite analyzes of MHC revealed that some MHC-I microsatellites correlated to stronger immune responses. These trials signify the first step in the long-term objective of releasing devils with immunity to DFTD into the wild.

There are very few articles in the literature describing continuous models of bacterial infections that mimic disease pathogenesis in humans and animals without using separate cohorts of animals at each stage of disease. In this work, we developed bioluminescent mouse models of partial-thickness scald wound infection and sepsis that mimic disease pathogenesis in humans and animals using a recombinant luciferase-expressing Staphylococcus aureus strain (Xen29). Two days post-scald wound infection, mice were treated twice daily with a 2% topical mupirocin ointment for 7 days. For sepsis experiments, mice were treated intraperitoneally with 6 mg/kg daptomycin 2 h and 6 h post-infection and time to moribund monitored for 72 h. Consistent bacterial burden data were obtained from individual mice by regular photon intensity quantification on a Xenogen IVIS Lumina XRMS Series III biophotonic imaging system, with concomitant significant reduction in photon intensities in drug-treated mice. Post-mortem histopathological examination of wounds and bacterial counts in blood correlated closely with disease severity and total flux obtained from Xen29. The bioluminescent murine models provide a refinement to existing techniques of multiple bacterial enumeration during disease pathogenesis and promote animal usage reduction. The models also provide an efficient and information-rich platform for preclinical efficacy evaluation of new drug classes for treating acute and chronic human and animal bacterial infections.

The monotremes (platypuses and echidnas) represent one of only four extant venomous mammalian lineages. Until recently, monotreme venom was poorly understood. However, the availability of the platypus genome and increasingly sophisticated genomic tools has allowed us to characterize platypus toxins, and provides a means of reconstructing the evolutionary history of monotreme venom. Here we review the physiology of platypus and echidna crural (venom) systems as well as pharmacological and genomic studies of monotreme toxins. Further, we synthesize current ideas about the evolution of the venom system, which in the platypus is likely to have been retained from a venomous ancestor, whilst being lost in the echidnas. We also outline several research directions and outstanding questions that would be productive to address in future research. An improved characterization of mammalian venoms will not only yield new toxins with potential therapeutic uses, but will also aid in our understanding of the way that this unusual trait evolves.

Tasmanian devil (Sarcophilus harrisii) pouch young, like other marsupials, are born underdeveloped and immunologically naïve, and are unable to mount an adaptive immune response. The mother’s milk provides nutrients for growth and development as well as providing passive immunity. To better understand immune response in this endangered species, we set out to characterise the genes involved in passive immunity by sequencing and annotating the transcriptome of a devil milk sample collected during mid-lactation. At mid-lactation we expect the young to have heightened immune responses, as they have emerged from the pouch, encountering new pathogens. A total of 233,660 transcripts were identified, including approximately 17,827 unique protein-coding genes and 846 immune genes. The most highly expressed transcripts were dominated by milk protein genes such as those encoding early lactation protein, late lactation proteins, α-lactalbumin, α-casein and β-casein. There were numerous highly expressed immune genes including lysozyme, whey acidic protein, ferritin and major histocompatibility complex I and II. Genes encoding immunoglobulins, antimicrobial peptides, chemokines and immune cell receptors were also identified. The array of immune genes identified in this study reflects the importance of the milk in providing immune protection to Tasmanian devil young and provides the first insight into Tasmanian devil milk.

Toll-like receptors (TLRs) play a crucial role in the early defence against invading pathogens, yet our understanding of TLRs in marsupial immunity is limited. Here, we describe the characterisation of nine TLRs from a koala immune tissue transcriptome and one TLR from a draft sequence of the koala genome and the subsequent development of an assay to study genetic diversity in these genes. We surveyed genetic diversity in 20 koalas from New South Wales, Australia and showed that one gene, TLR10 is monomorphic, while the other nine TLR genes have between two and 12 alleles. 40 SNPs (16 non-synonymous) were identified across the ten TLR genes. These markers provide a springboard to future studies on innate immunity in the koala, a species under threat from two major infectious diseases.

* The swift parrot ( Lathamus discolor) is a Critically Endangered migratory parrot that breeds in Tasmania and winters on the Australian mainland. Here we provide a reference genome assembly for the swift parrot. We sequence PacBio HiFi reads to create a high-quality reference assembly and identify a complete mitochondrial sequence. We also generate a reference transcriptome from five organs to inform genome annotation. The genome was 1.24 Gb in length and consisted of 847 contigs with a contig N50 of 18.97 Gb and L50 of 20 contigs. This study provides an annotated reference assembly and transcriptomic resources for the swift parrot to assist in future conservation genomic research.

The Major Histocompatibility Complex (MHC) class I is a highly variable gene family that encodes cell-surface receptors vital for recognition of intracellular pathogens and initiation of immune responses. The MHC class I has yet to be characterised in bufonid toads (Order: Anura; Suborder: Neobatrachia; Family: Bufonidae), a large and diverse family of anurans. Here we describe the characterisation of a classical MHC class I gene in the Australian cane toad, Rhinella marina. From 25 individuals sampled from the Australian population, we found only 3 alleles at this classical class I locus. We also found large number of class I alpha 1 alleles, implying an expansion of class I loci in this species. The low classical class I genetic diversity is likely the result of repeated bottleneck events, which arose as a result of the cane toad's complex history of introductions as a biocontrol agent and its subsequent invasion across Australia.

Monotremes (echidna and platypus) are egg-laying mammals. One of their most unique characteristic is that males have venom/crural glands that are seasonally active. Male platypuses produce venom during the breeding season, delivered via spurs, to aid in competition against other males. Echidnas are not able to erect their spurs, but a milky secretion is produced by the gland during the breeding season. The function and molecular composition of echidna venom is as yet unknown. Hence, we compared the deeply sequenced transcriptome of an in-season echidna crural gland to that of a platypus and searched for putative venom genes to provide clues into the function of echidna venom and the evolutionary history of monotreme venom. We found that the echidna venom gland transcriptome was markedly different from the platypus with no correlation between the top 50 most highly expressed genes. Four peptides found in the venom of the platypus were detected in the echidna transcriptome. However, these genes were not highly expressed in echidna, suggesting that they are the remnants of the evolutionary history of the ancestral venom gland. Gene ontology terms associated with the top 100 most highly expressed genes in echidna, showed functional terms associated with steroidal and fatty acid production, suggesting that echidna "venom" may play a role in scent communication during the breeding season. The loss of the ability to erect the spur and other unknown evolutionary forces acting in the echidna lineage resulted in the gradual decay of venom components and the evolution of a new role for the crural gland.

Devil Facial Tumour Disease (DFTD) is a unique clonal cancer that threatens the world's largest carnivorous marsupial, the Tasmanian devil (Sarcophilus harrisii) with extinction. This transmissible cancer is passed between individual devils by cell implantation during social interactions. The tumour arose in a Schwann cell of a single devil over 15 years ago and since then has expanded clonally, without showing signs of replicative senescence; in stark contrast to a somatic cell that displays a finite capacity for replication, known as the "Hayflick limit".In the present study we investigate the role of telomere length, measured as Telomere Copy Number (TCN), and telomerase and shelterin gene expression, as well as telomerase activity in maintaining hyperproliferation of Devil Facial Tumour (DFT) cells. Our results show that DFT cells have short telomeres. DFTD TCN does not differ between geographic regions or between strains. However, TCN has increased over time. Unlimited cell proliferation is likely to have been achieved through the observed up-regulation of the catalytic subunit of telomerase (TERT) and concomitant activation of telomerase. Up-regulation of the central component of shelterin, the TRF1-intercating nuclear factor 2 (TINF2) provides DFT a mechanism for telomere length homeostasis. The higher expression of both TERT and TINF2 may also protect DFT cells from genomic instability and enhance tumour proliferation.DFT cells appear to monitor and regulate the length of individual telomeres: i.e. shorter telomeres are elongated by up-regulation of telomerase-related genes; longer telomeres are protected from further elongation by members of the shelterin complex, which may explain the lack of spatial and strain variation in DFT telomere copy number. The observed longitudinal increase in gene expression in DFT tissue samples and telomerase activity in DFT cell lines might indicate a selection for more stable tumours with higher proliferative potential.

Devil facial tumour disease (DFTD) is a fatal contagious cancer that has decimated Tasmanian devil populations. The tumour has spread without invoking immune responses, possibly due to low levels of Major Histocompatibility Complex (MHC) diversity in Tasmanian devils. Animals from a region in north-western Tasmania have lower infection rates than those in the east of the state. This area is a genetic transition zone between sub-populations, with individuals from north-western Tasmania displaying greater diversity than eastern devils at MHC genes, primarily through MHC class I gene copy number variation. Here we test the hypothesis that animals that remain healthy and tumour free show predictable differences at MHC loci compared to animals that develop the disease.We compared MHC class I sequences in 29 healthy and 22 diseased Tasmanian devils from West Pencil Pine, a population in north-western Tasmania exhibiting reduced disease impacts of DFTD. Amplified alleles were assigned to four loci, Saha-UA, Saha-UB, Saha-UC and Saha-UD based on recently obtained genomic sequence data. Copy number variation (caused by a deletion) at Saha-UA was confirmed using a PCR assay. No association between the frequency of this deletion and disease status was identified. All individuals had alleles at Saha-UD, disproving theories of disease susceptibility relating to copy number variation at this locus. Genetic variation between the two sub-groups (healthy and diseased) was also compared using eight MHC-linked microsatellite markers. No significant differences were identified in allele frequency, however differences were noted in the genotype frequencies of two microsatellites located near non-antigen presenting genes within the MHC.We did not find predictable differences in MHC class I copy number variation to account for differences in susceptibility to DFTD. Genotypic data was equivocal but indentified genomic areas for further study.

Devil facial tumour disease (DFTD) is a fatal, transmissible malignancy that threatens the world's largest marsupial carnivore, the Tasmanian devil, with extinction. First recognised in 1996, DFTD has had a catastrophic effect on wild devil numbers, and intense research efforts to understand and contain the disease have since demonstrated that the tumour is a clonal cell line transmitted by allograft. We used chromosome painting and gene mapping to deconstruct the DFTD karyotype and determine the chromosome and gene rearrangements involved in carcinogenesis. Chromosome painting on three different DFTD tumour strains determined the origins of marker chromosomes and provided a general overview of the rearrangement in DFTD karyotypes. Mapping of 105 BAC clones by fluorescence in situ hybridisation provided a finer level of resolution of genome rearrangements in DFTD strains. Our findings demonstrate that only limited regions of the genome, mainly chromosomes 1 and X, are rearranged in DFTD. Regions rearranged in DFTD are also highly rearranged between different marsupials. Differences between strains are limited, reflecting the unusually stable nature of DFTD. Finally, our detailed maps of both the devil and tumour karyotypes provide a physical framework for future genomic investigations into DFTD.

BackgroundTo overcome the increasing resistance of pathogens to existing antibiotics the 10×'20 Initiative declared the urgent need for a global commitment to develop 10 new antimicrobial drugs by the year 2020. Naturally occurring animal antibiotics are an obvious place to start. The recently sequenced genomes of mammals that are divergent from human and mouse, including the tammar wallaby and the platypus, provide an opportunity to discover novel antimicrobials. Marsupials and monotremes are ideal potential sources of new antimicrobials because they give birth to underdeveloped immunologically naïve young that develop outside the sterile confines of a uterus in harsh pathogen-laden environments. While their adaptive immune system develops innate immune factors produced either by the mother or by the young must play a key role in protecting the immune-compromised young. In this study we focus on the cathelicidins, a key family of antimicrobial peptide genes.Principal findingWe identified 14 cathelicidin genes in the tammar wallaby genome and 8 in the platypus genome. The tammar genes were expressed in the mammary gland during early lactation before the adaptive immune system of the young develops, as well as in the skin of the pouch young. Both platypus and tammar peptides were effective in killing a broad range of bacterial pathogens. One potent peptide, expressed in the early stages of tammar lactation, effectively killed multidrug-resistant clinical isolates of Pseudomonas aeruginosa, Klebsiella pneumoniae and Acinetobacter baumannii.Conclusions and significanceMarsupial and monotreme young are protected by antimicrobial peptides that are potent, broad spectrum and salt resistant. The genomes of our distant relatives may hold the key for the development of novel drugs to combat multidrug-resistant pathogens.

Tasmanian devils (Sarcophilus harrisii) are on the verge of extinction due to a transmissible cancer, devil facial tumour disease (DFTD). This tumour is an allograft that is transmitted between individuals without immune recognition of the tumour cells. The mechanism to explain this lack of immune recognition and acceptance is not well understood. It has been hypothesized that lack of genetic diversity at the Major Histocompatibility Complex (MHC) allowed the tumour cells to grow in genetically similar hosts without evoking an immune response to alloantigens. We conducted mixed lymphocyte reactions and skin grafts to measure functional MHC diversity in the Tasmanian devil population. The limited MHC diversity was sufficient to produce measurable mixed lymphocyte reactions. There was a wide range of responses, from low or no reaction to relatively strong responses. The highest responses occurred when lymphocytes from devils from the east of Tasmania were mixed with lymphocytes from devils from the west of Tasmania. All of the five successful skin allografts were rejected within 14 days after surgery, even though little or no MHC I and II mismatches were found. Extensive T-cell infiltration characterised the immune rejection. We conclude that Tasmanian devils are capable of allogeneic rejection. Consequently, a lack of functional allorecognition mechanisms in the devil population does not explain the transmission of a contagious cancer.

* The brown thornbill ( Acanthiza pusilla) is a songbird endemic to eastern Australia with five recognised subspecies within the brown thornbill. The most notable is the King Island brown thornbill ( Acanthiza pusilla magnirostris) of which there are less than 100 remaining and based on expert elicitation are the most likely Australian bird to become extinct in the next 20 years. We sequenced PacBio HiFi reads of the brown thornbill to generate a high-quality reference genome 1.25Gb in size and contig N50 of 20.1Mb. Additionally, we sequenced mRNA from three tissues to generate a global transcriptome to aid with genome annotation. The generation of a reference genome for the brown thornbill provides an important resource to align additional genomic data which will be produced in the near future.

* The swift parrot ( Lathamus discolor) is a Critically Endangered migratory parrot that breeds in Tasmania and winters on the Australian mainland. Here we provide a reference genome assembly for the swift parrot. We sequence PacBio HiFi reads to create a high-quality reference assembly and identify a complete mitochondrial sequence. We also generate a reference transcriptome from five organs to inform genome annotation. The genome was 1.24 Gb in length and consisted of 847 contigs with a contig N50 of 18.97 Gb and L50 of 20 contigs. This study provides an annotated reference assembly and transcriptomic resources for the swift parrot to assist in future conservation genomic research.

The first sequenced marsupial genome promises to reveal unparalleled insights into mammalian evolution. We have used the Monodelphis domestica (gray short-tailed opossum) sequence to construct the first map of a marsupial major histocompatibility complex (MHC). The MHC is the most gene-dense region of the mammalian genome and is critical to immunity and reproductive success. The marsupial MHC bridges the phylogenetic gap between the complex MHC of eutherian mammals and the minimal essential MHC of birds. Here we show that the opossum MHC is gene dense and complex, as in humans, but shares more organizational features with non-mammals. The Class I genes have amplified within the Class II region, resulting in a unique Class I/II region. We present a model of the organization of the MHC in ancestral mammals and its elaboration during mammalian evolution. The opossum genome, together with other extant genomes, reveals the existence of an ancestral "immune supercomplex" that contained genes of both types of natural killer receptors together with antigen processing genes and MHC genes.

The Kroombit tinkerfrog ( Taudactylus pleione) is a stream-dwelling amphibian of the Myobatrachidae family. It is listed as Critically Endangered and is at high risk of extinction due to chytridiomycosis. Here, we provide the first genome assembly of the evolutionarily distinct Taudactylus genus. We sequenced PacBio HiFi reads to assemble a high-quality long-read genome and identified the mitochondrial genome. We also generated a global transcriptome from a tadpole to improve gene annotation. The genome was 5.52 Gb in length and consisted of 4,196 contigs with a contig N50 of 8.853 Mb and an L50 of 153. This study provides the first genomic resources for the Kroombit tinkerfrog to assist in future phylogenetic, environmental DNA, conservation breeding, and disease susceptibility studies.

The Kroombit tinkerfrog ( Taudactylus pleione) is a stream-dwelling amphibian of the Myobatrachidae family. It is listed as Critically Endangered and is at high risk of extinction due to chytridiomycosis. Here, we provide the first genome assembly of the evolutionarily distinct Taudactylus genus. We sequenced PacBio HiFi reads to assemble a high-quality long-read genome and identified the mitochondrial genome. We also generated a global transcriptome from a tadpole to improve gene annotation. The genome was 5.52 Gb in length and consisted of 4,196 contigs with a contig N50 of 8.853 Mb and an L50 of 153. This study provides the first genomic resources for the Kroombit tinkerfrog to assist in future phylogenetic, environmental DNA, conservation breeding, and disease susceptibility studies.