No evidence of inbreeding depression in a Tasmanian devil insurance population despite significant variation in inbreeding

Population genomics of the muskox' resilience in the near absence of genetic variation

Genomic studies of species threatened by extinction are providing crucial information about evolutionary mechanisms and genetic consequences of population declines and bottlenecks. However, to understand how species avoid the extinction vortex, insights can be drawn by studying species that thrive despite past declines. Here, we studied the population genomics of the muskox (Ovibos moschatus), an Ice Age relict that was at the brink of extinction for thousands of years at the end of the Pleistocene yet appears to be thriving today. We analysed 108 whole genomes, including present-day individuals representing the current native range of both muskox subspecies, the white-faced and the barren-ground muskox (O. moschatus wardi and O. moschatus moschatus) and a ~21,000-year-old ancient individual from Siberia. We found that the muskox' demographic history was profoundly shaped by past climate changes and post-glacial re-colonizations. In particular, the white-faced muskox has the lowest genome-wide heterozygosity recorded in an ungulate. Yet, there is no evidence of inbreeding depression in native muskox populations. We hypothesize that this can be explained by the effect of long-term gradual population declines that allowed for purging of strongly deleterious mutations. This study provides insights into how species with a history of population bottlenecks, small population sizes and low genetic diversity survive against all odds.

Genomics reveals introgression and purging of deleterious mutations in the Arabian leopard (Panthera pardus nimr)

In endangered species, low-genetic variation and inbreeding result from recent population declines. Genetic screenings in endangered populations help to assess their vulnerability to extinction and to create informed management actions toward their conservation efforts. The leopard, Panthera pardus, is a highly generalist predator with currently eight different subspecies. Yet, genomic data are still lacking for the Critically Endangered Arabian leopard (P. p. nimr). Here, we sequenced the whole genome of two Arabian leopards and assembled the most complete genomic dataset for leopards to date. Our phylogenomic analyses show that leopards are divided into two deeply divergent clades: the African and the Asian. Conservation genomic analyses indicate a prolonged population decline, which has led to an increase in inbreeding and runs of homozygosity, with consequent purging of deleterious mutations in both Arabian individuals. Our study represents the first attempt to genetically inform captive breeding programmes for this Critically Endangered subspecies.

The relevance of pedigrees in the conservation genomics era

Over the past 50 years conservation genetics has developed a substantive toolbox to inform species management. One of the most long‐standing tools available to manage genetics—the pedigree—has been widely used to characterize diversity and maximize evolutionary potential in threatened populations. Now, with the ability to use high throughput sequencing to estimate relatedness, inbreeding, and genome‐wide functional diversity, some have asked whether it is warranted for conservation biologists to continue collecting and collating pedigrees for species management. In this perspective, we argue that pedigrees remain a relevant tool, and when combined with genomic data, create an invaluable resource for conservation genomic management. Genomic data can address pedigree pitfalls (e.g., founder relatedness, missing data, uncertainty), and in return robust pedigrees allow for more nuanced research design, including well‐informed sampling strategies and quantitative analyses (e.g., heritability, linkage) to better inform genomic inquiry. We further contend that building and maintaining pedigrees provides an opportunity to strengthen trusted relationships among conservation researchers, practitioners, Indigenous Peoples, and Local Communities.

First evidence of deviation from Mendelian proportions in a conservation programme

Classic Mendelian inheritance is the bedrock of population genetics and underpins pedigree-based management of animal populations. However, assumptions of Mendelian inheritance might not be upheld in conservation breeding programmes if early viability selection occurs, even when efforts are made to equalise genetic contributions of breeders. To test this possibility, we investigated deviations from Mendelian proportions in a captive metapopulation of the endangered Tasmanian devil. This marsupial population is ideal for addressing evolutionary questions in conservation due to its large size, range of enclosure types (varying in environmental conditions), good genomic resources (which aid interpretation), and the species' biology. Devil mothers give birth to more offspring than they can nurse in the pouch, providing the potential for intense viability selection amongst embryos. We used data from 140 known sire-dam-offspring triads to isolate within-family selection from population-level mechanisms (such as mate choice or inbreeding), and compared observed offspring genotypes at 123 targeted SNPs to neutral (i.e., Mendelian) expectations. We found lower offspring heterozygosity than expected, and subtle patterns that varied across a gradient of management intensity from zoo-like enclosures to semi-wild environments for some loci. Meiotic drive or maternal-foetal incompatibilities are consistent with our results, although we cannot statistically confirm these mechanisms. We found some evidence that maternal genotype affects annual litter size, suggesting that family-level patterns are driven by differential offspring mortality before birth or during early development. Our results show that deviations from Mendelian inheritance can occur in conservation programmes, despite best-practice management to prevent selection.

A Tasmanian devil breeding program to support wild recovery

Tasmanian devils are threatened in the wild by devil facial tumour disease: a transmissible cancer with a high fatality rate. In response, the Save the Tasmanian Devil Program (STDP) established an 'insurance population' to enable the preservation of genetic diversity and natural behaviours of devils. This breeding program includes a range of institutions and facilities, from zoo-based intensive enclosures to larger, more natural environments, and a strategic approach has been required to capture and maintain genetic diversity, natural behaviours and to ensure reproductive success. Laboratory-based research, particularly genetics, in tandem with adaptive management has helped the STDP reach its goals, and has directly contributed to the conservation of the species in the wild. Here we review this work and show that the Tasmanian devil breeding program is a powerful example of how genetic research can be used to understand and improve reproductive success in a threatened species.

Comparison of genomic diversity and structure of sable antelope (Hippotragus niger) in zoos, conservation centers, and private ranches in North America

As we enter the sixth mass extinction, many species that are no longer self-sustaining in their natural habitat will require ex situ management. Zoos have finite resources for ex situ management, and there is a need for holistic conservation programs between the public and private sector. Ex situ populations of sable antelope, Hippotragus niger, have existed in zoos and privately owned ranches in North America since the 1910s. Unknown founder representation and relatedness has made the genetic management of this species challenging within zoos, while populations on privately owned ranches are managed independently and retain minimal-to-no pedigree history. Consequences of such challenges include an increased risk of inbreeding and a loss of genetic diversity. Here, we developed and applied a customized targeted sequence capture panel based on 5,000 genomewide single-nucleotide polymorphisms to investigate the genomic diversity present in these uniquely managed populations. We genotyped 111 sable antelope: 23 from zoos, 43 from a single conservation center, and 45 from ranches. We found significantly higher genetic diversity and significantly lower inbreeding in herds housed in zoos and conservation centers, when compared to those in privately owned ranches, likely due to genetic-based breeding recommendations implemented in the former populations. Genetic clustering was strong among all three populations, possibly as a result of genetic drift. We propose that the North American ex situ population of sable antelope would benefit from a metapopulation management system, to halt genetic drift, reduce the occurrence of inbreeding, and enable sustainable population sizes to be managed ex situ.

Genetic monitoring of the greater stick-nest rat meta-population for strategic supplementation planning

Translocation is an increasingly common component of species conservation efforts. However, translocated populations often suffer from loss of genetic diversity and increased inbreeding, and thus may require active management to establish gene flow across isolated populations. Assisted gene flow can be laborious and costly, so recipient and source populations should be carefully chosen to maximise genetic diversity outcomes. The greater stick-nest rat (GSNR, Leporillus conditor), a threatened Australian rodent, has been the focus of a translocation program since 1985, resulting in five extant translocated populations (St Peter Island, Reevesby Island, Arid Recovery, Salutation Island and Mt Gibson), all derived from a remnant wild population on the East and West Franklin Islands. We evaluated the genetic diversity in all extant GSNR populations using a large single nucleotide polymorphism dataset with the explicit purpose of informing future translocation planning. Our results show varying levels of genetic divergence, inbreeding and loss of genetic diversity in all translocated populations relative to the remnant source on the Franklin Islands. All translocated populations would benefit from supplementation to increase genetic diversity, but two—Salutation Island and Mt Gibson—are of highest priority. We recommend a targeted admixture approach, in which animals for supplementation are sourced from populations that have low relatedness to the recipient population. Subject to assessment of contemporary genetic diversity, St Peter Island and Arid Recovery are the most appropriate source populations for genetic supplementation. Our study demonstrates an effective use of genetic surveys for data-driven management of threatened species.

Inbreeding depression in one of the last DFTD-free wild populations of Tasmanian devils

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.

A comparison of pedigree, genetic and genomic estimates of relatedness for informing pairing decisions in two critically endangered birds: Implications for conservation breeding programmes worldwide

Conservation management strategies for many highly threatened species include conservation breeding to prevent extinction and enhance recovery. Pairing decisions for these conservation breeding programmes can be informed by pedigree data to minimize relatedness between individuals in an effort to avoid inbreeding, maximize diversity and maintain evolutionary potential. However, conservation breeding programmes struggle to use this approach when pedigrees are shallow or incomplete. While genetic data (i.e., microsatellites) can be used to estimate relatedness to inform pairing decisions, emerging evidence indicates this approach may lack precision in genetically depauperate species, and more effective estimates will likely be obtained from genomic data (i.e., thousands of genome-wide single nucleotide polymorphisms, or SNPs). Here, we compare relatedness estimates and subsequent pairing decisions using pedigrees, microsatellites and SNPs from whole-genome resequencing approaches in two critically endangered birds endemic to New Zealand: kakī/black stilt (Himantopus novaezelandiae) and kākāriki karaka/orange-fronted parakeet (Cyanoramphus malherbi). Our findings indicate that SNPs provide more precise estimates of relatedness than microsatellites when assessing empirical parent-offspring and full sibling relationships. Further, our results show that relatedness estimates and subsequent pairing recommendations using PMx are most similar between pedigree- and SNP-based approaches. These combined results indicate that in lieu of robust pedigrees, SNPs are an effective tool for informing pairing decisions, which has important implications for many poorly pedigreed conservation breeding programmes worldwide.

The Value of Reference Genomes in the Conservation of Threatened Species

Conservation initiatives are now more crucial than ever-over a million plant and animal species are at risk of extinction over the coming decades. The genetic management of threatened species held in insurance programs is recommended; however, few are taking advantage of the full range of genomic technologies available today. Less than 1% of the 13505 species currently listed as threated by the International Union for Conservation of Nature (IUCN) have a published genome. While there has been much discussion in the literature about the importance of genomics for conservation, there are limited examples of how having a reference genome has changed conservation management practice. The Tasmanian devil (Sarcophilus harrisii), is an endangered Australian marsupial, threatened by an infectious clonal cancer devil facial tumor disease (DFTD). Populations have declined by 80% since the disease was first recorded in 1996. A reference genome for this species was published in 2012 and has been crucial for understanding DFTD and the management of the species in the wild. Here we use the Tasmanian devil as an example of how a reference genome has influenced management actions in the conservation of a species.

Genomic pedigree reconstruction identifies predictors of mating and reproductive success in an invasive vertebrate

The persistence of an invasive species is influenced by its reproductive ecology, and a successful control program must operate on this premise. However, the reproductive ecology of invasive species may be enigmatic due to factors that also limit their management, such as cryptic coloration and behavior. We explored the mating and reproductive ecology of the invasive Brown Treesnake (BTS: Boiga irregularis) by reconstructing a multigenerational genomic pedigree based on 654 single nucleotide polymorphisms for a geographically closed population established in 2004 on Guam (N = 426). The pedigree allowed annual estimates of individual mating and reproductive success to be inferred for snakes in the study population over a 14-year period. We then employed generalized linear mixed models to gauge how well phenotypic and genomic data could predict sex-specific annual mating and reproductive success. Average snout-vent length (SVL), average body condition index (BCI), and trappability were significantly related to annual mating success for males, with average SVL also related to annual mating success for females. Male and female annual reproductive success was positively affected by SVL, BCI, and trappability. Surprisingly, the degree to which individuals were inbred had no effect on annual mating or reproductive success. When juxtaposed with current control methods, these results indicate that baited traps, a common interdiction tool, may target fecund BTS in some regards but not others. Our study emphasizes the importance of reproductive ecology as a focus for improving BTS control and promotes genomic pedigree reconstruction for such an endeavor in this invasive species and others.

MHC-associated mate choice under competitive conditions in captive versus wild Tasmanian devils

Mate choice contributes to driving evolutionary processes when animals choose breeding partners that confer genetic advantages to offspring, such as increased immunocompetence. The major histocompatibility complex (MHC) is an important group of immunological molecules, as MHC antigens bind and present foreign peptides to T-cells. Recent studies suggest that mates may be selected based on their MHC profile, leading to an association between an individual’s MHC diversity and their breeding success. In conservation, it may be important to consider mate choice in captive breeding programs, as this mechanism may improve reproductive rates. We investigated the reproductive success of Tasmanian devils in a group housing facility to determine whether increased MHC-based heterozygosity led individuals to secure more mating partners and produce more offspring. We also compared the breeding success of captive females to a wild devil population. MHC diversity was quantified using 12 MHC-linked microsatellite markers, including 11 previously characterized markers and one newly identified marker. Our analyses revealed that there was no relationship between MHC-linked heterozygosity and reproductive success either in captivity or the wild. The results of this study suggest that, for Tasmanian devils, MHC-based heterozygosity does not produce greater breeding success and that no specific changes to current captive management strategies are required with respect to preserving MHC diversity.

Multiple paternity and precocial breeding in wild Tasmanian devils, Sarcophilus harrisii (Marsupialia: Dasyuridae)

Polyandry, a common reproductive strategy in various animal species, has potential female benefits, which include enhanced offspring fitness. Benefits can be direct, such as reduced risk of male infanticide of offspring, or indirect, such as increased genetic diversity of offspring and the acquisition of ‘good genes’. Multiple paternity of litters has been recorded in numerous marsupial species but has not been reported in Tasmanian devils, Sarcophilus harrisii (Boitard). We investigated whether multiple paternity occurred in litters within a wild population of Tasmanian devils. Using major histocompatibility complex-linked and neutral microsatellite markers, the paternity of nine litters was analysed. We found multiple paternity in four out of nine litters and that yearling (> 1, < 2 years old) male devils were siring offspring. This is the first record of multiple paternity and of male precocial breeding in wild Tasmanian devils. To date, there are no data relating to the subsequent survival of devils from single- vs. multiple-sired litters; therefore, we do not know whether multiple paternity increases offspring survival in the wild. These results have implications for the Tasmanian devil captive insurance programme, because group housing can lead to multiple-sired litters, making the maintenance of genetic diversity over time difficult to manage.

Complex problems need detailed solutions: Harnessing multiple data types to inform genetic management in the wild

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.

Conserving adaptive potential: lessons from Tasmanian devils and their transmissible cancer

Maintenance of adaptive genetic variation has long been a goal of management of natural populations, but only recently have genomic tools allowed identification of specific loci associated with fitness-related traits in species of conservation concern. This raises the possibility of managing for genetic variation directly relevant to specific threats, such as those due to climate change or emerging infectious disease. Tasmanian devils (Sarcophilus harrisii) face the threat of a transmissible cancer, devil facial tumor disease (DFTD), that has decimated wild populations and led to intensive management efforts. Recent discoveries from genomic and modeling studies reveal how natural devil populations are responding to DFTD, and can inform management of both captive and wild devil populations. Notably, recent studies have documented genetic variation for disease-related traits and rapid evolution in response to DFTD, as well as potential mechanisms for disease resistance such as immune response and tumor regression in wild devils. Recent models predict dynamic persistence of devils with or without DFTD under a variety of modeling scenarios, although at much lower population densities than before DFTD emerged, contrary to previous predictions of extinction. As a result, current management that focuses on captive breeding and release for maintaining genome-wide genetic diversity or demographic supplementation of populations could have negative consequences. Translocations of captive devils into wild populations evolving with DFTD can cause outbreeding depression and/or increases in the force of infection and thereby the severity of the epidemic, and we argue that these risks outweigh any benefits of demographic supplementation in wild populations. We also argue that genetic variation at loci associated with DFTD should be monitored in both captive and wild populations, and that as our understanding of DFTD-related genetic variation improves, considering genetic management approaches to target this variation is warranted in developing conservation strategies for Tasmanian devils.

Thylacine and Tasmanian devil: between hope and reality – a lesson to be learnt from Google Trends search data

The two iconic Tasmanian species, the Tasmanian devil (Sarcophilus harrisii) and the thylacine (Thylacinus cynocephalus), are of great interest to the general public and the media. The most likely extinct Tasmanian wolf or tiger, the thylacine, symbolises human responsibility for nature and species conservation and inspired the ‘National Threatened Species Day’, which commemorates the death of the last thylacine at Beaumaris Zoo in Hobart on 7 September 1936 to raise awareness of endangered plants and animals. Since the spread of the Devil Facial Tumour Disease critically endangered the survival of the largest remaining native carnivore (S. harrisii) today, this has generated both scientific interest and the interest of the general public. Google Trends has already been used as a tool for documenting and investigating the information needs and concerns of the population, as has been shown using the example of diseases. In this study, Google Trends data were used to examine the seasonality of the search term ‘thylacine sightings’ and the development of the frequency of different search terms in the period between 2004 and 2020. As a result, relative search intensities for ‘thylacine cloning’ and ‘cloning extinct species’ have shown a decrease over time. While Google Trends cannot clearly determine search motivation, search terms can be selected for the examinations that document more hope or a rational need for information or concern.

The International Mouse Phenotyping Consortium (IMPC): a functional catalogue of the mammalian genome that informs conservation

The International Mouse Phenotyping Consortium (IMPC) is building a catalogue of mammalian gene function by producing and phenotyping a knockout mouse line for every protein-coding gene. To date, the IMPC has generated and characterised 5186 mutant lines. One-third of the lines have been found to be non-viable and over 300 new mouse models of human disease have been identified thus far. While current bioinformatics efforts are focused on translating results to better understand human disease processes, IMPC data also aids understanding genetic function and processes in other species. Here we show, using gorilla genomic data, how genes essential to development in mice can be used to help assess the potentially deleterious impact of gene variants in other species. This type of analyses could be used to select optimal breeders in endangered species to maintain or increase fitness and avoid variants associated to impaired-health phenotypes or loss-of-function mutations in genes of critical importance. We also show, using selected examples from various mammal species, how IMPC data can aid in the identification of candidate genes for studying a condition of interest, deliver information about the mechanisms involved, or support predictions for the function of genes that may play a role in adaptation. With genotyping costs decreasing and the continued improvements of bioinformatics tools, the analyses we demonstrate can be routinely applied.

Population monitoring of an urban gliding mammal in eastern Australia

Long-term monitoring is an important element of species conservation. This study describes changes in the size of a squirrel glider (Petaurus norfolcensis) population over a 10-year period. The population occupied a 45-ha forest remnant within the urban area of Brisbane. Gliders were tagged from 25 nights of trapping during 2006–08 and from 16 nights of trapping in 2015. Population modelling was used to estimate adult population size. This suggested the adult population comprised 30–40 individuals at the beginning and end of the 10-year period. It reached a peak of 70 individuals in mid-2007. These data suggest that the study area contains a small population that is prone to interannual variation but there was no evidence of it being in decline. Survival estimates during 2006–08 were equivalent to those estimated for a larger population in Victoria. Population monitoring should be continued to determine how resilient this population is to population decline and to investigate factors that may cause decline. This study provides an example of an approach that could be used to monitor threatened populations of the squirrel glider.

Evidence of opposing fitness effects of parental heterozygosity and relatedness in a critically endangered marine turtle?

How individual genetic variability relates to fitness is important in understanding evolution and the processes affecting populations of conservation concern. Heterozygosity-fitness correlations (HFCs) have been widely used to study this link in wild populations, where key parameters that affect both variability and fitness, such as inbreeding, can be difficult to measure. We used estimates of parental heterozygosity and genetic similarity ('relatedness') derived from 32 microsatellite markers to explore the relationship between genetic variability and fitness in a population of the critically endangered hawksbill turtle, Eretmochelys imbricata. We found no effect of maternal MLH (multilocus heterozygosity) on clutch size or egg success rate, and no single-locus effects. However, we found effects of paternal MLH and parental relatedness on egg success rate that interacted in a way that may result in both positive and negative effects of genetic variability. Multicollinearity in these tests was within safe limits, and null simulations suggested that the effect was not an artefact of using paternal genotypes reconstructed from large samples of offspring. Our results could imply a tension between inbreeding and outbreeding depression in this system, which is biologically feasible in turtles: female-biased natal philopatry may elevate inbreeding risk and local adaptation, and both processes may be disrupted by male-biased dispersal. Although this conclusion should be treated with caution due to a lack of significant identity disequilibrium, our study shows the importance of considering both positive and negative effects when assessing how variation in genetic variability affects fitness in wild systems.