- Research article
- Open Access
The effects of group versus intensive housing on the retention of genetic diversity in insurance populations
BMC Zoologyvolume 3, Article number: 2 (2018)
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).
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.
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.
Many threatened species require some form of population management as a conservation tool to prevent extinction. Topical reviews have found that conservation breeding or insurance populations have recently played a role in threat level reduction of endangered species [1, 2]. Combined with in-situ management, insurance populations prevent extinction until wild reintroductions are achievable as seen with California condors, Gymnogyps californianus , golden-lion tamarins, Leontopithecus rosalia  and black-footed ferrets, Mustela nigripes . The primary goal of an insurance population is to capture and maintain representative wild genetic diversity, as well as maintain demographic stability, for as long as possible [6,7,8]. Population genetic diversity has individual fitness benefits (e.g. reproductive success), in addition to population fitness benefits (e.g. potential to adapt to environmental changes), both being important for populations to be used as sources for release to the wild . Typically, insurance populations commence when a species is in dire straits and so accessing suitable numbers of unrelated founders is problematic as populations are usually fragmented and declining. Genetic diversity in captivity is further lost when founder representation is unequal . For decades the global zoo community has utilized a mean kinship strategy (the average relatedness of an individual to the population, calculated via pedigrees), pairing individuals who are least related to each other to maintain genetic diversity [11, 12]. Traditionally this meant that individual pairs were housed together for breeding purposes, a type of forced monogamy, regardless of the species life history and without the option of mate choice. This practice of forced monogamy is not realistic for all species (e.g. giant panda ), and often deviates from the natural ecology of the species’ wild habitat, so recent discussion has occurred around whether providing mate choice, via group housing of multiple individuals, to improve conservation breeding .
The literature provides numerous illustrations of difficulties in conservation breeding and the resulting impact on genetic diversity. Higher reproductive outputs from certain individuals results in higher genetic contributions to subsequent generations, and lower genetic contributions from less successful individuals . For example, in the clouded leopard (Neofelis nebulosa), incompatibilities among forced breeding pairs (aggressive behaviour resulting in death or serious injury) created unequal founder representation, as only 27% of males and 20% of females successfully reproduced . Furthermore, genetic diversity can also be lost as a result of genetic drift in small, finite populations . Genetic drift becomes increasingly influential with decreasing population sizes and insurance populations tend to be small in size . In the absence of selection from inadequate resources, disease and harsh weather, often encountered in a species' natural habitat, genetic drift leads to random fixation and loss of alleles in insurance populations . Simulations with population sizes similar to those considered ideal for insurance populations demonstrated that genetic drift can rapidly deplete genetic diversity; however, the immigration of only a few individuals annually had the potential to slow or halt the effects of genetic drift . The loss of genetic diversity arising from these two processes, mate incompatibilities and genetic drift, impedes conservation goals.
In conservation breeding, allowing mate choice may avoid reproductive failure that results from incompatible mate pairing, as breeders will have opportunity to select preferred mates, increasing overall productivity of the population. Mate choice can increase offspring fitness (quantity or quality) via direct behavioral benefits (parental investment) or indirect genetic benefits (parental compatibility) [19, 20]. In some species the reproductive rate of genetically desired individuals that fail to breed when force-paired, have increased with the addition of mate-choice management (e.g. cheetahs, Acinonyx jubatus  and zebra finch, Taeniopygia guttata ). A potential downside of allowing mate choice in an insurance population is increased reproductive skew among individuals. If certain individuals are consistently chosen as the preferred mate, it can result in the exclusion of non-preferred individuals  leading to unequal founder representation and loss of genetic diversity . Unequal founder representation in insurance populations occurs when certain founding individuals have low to absent genetic contributions to the subsequent generations in captivity, resulting in their genetic line being poorly represented relative to others, or completely extinguished. When founding lines become extinguished, the population has reduced adaptive potential, and so maintaining equal founder representation increases the likelihood of adaptation to novel environments, diseases and pollutants. Long-term survival for small populations can be potentially compromised if management is exclusively based on mate-choice breeding.
A manner in which to approach these issues in conservation breeding programs is to create a captive environment similar to a species’ wild environment which enables the inclusion of both mate-choice breeding preferences and some pedigree-based breeding. That is, housing low mean kinship, unrelated individuals together in group enclosures allows for individual mating preference to take place while still moderately controlling for equal genetic representation. Experimentally determining mate-choice preferences prior to breeding is logistically unachievable for large-scale insurance populations. As such, housing selected individuals in groups may be a way forward for conservation breeding that is efficient and financially achievable, minimises behavioural adaptation whilst boosting breeding success. In addition, reproductive success may improve as a result of mate preferences and intrasexual competition. A characteristic of male (and to a lesser degree, female) reproductive success is mate acquisition ability [23, 24]. Individuals in group enclosures therefore may obtain a disproportionately high reproductive output as a result of social hierarchies in addition to mate-preferences . Until recently, it has been difficult to explore these evolutionary processes in respect to conservation management, given the unknown complexity of mate choice across taxa and the logistic requirements within a captive-breeding setting.
In this study, we compare two conservation breeding approaches 1) forced monogamous pairings to equalize founder representation, and 2) opportunity for mate-choice. We use the Tasmanian devil (Sarcophilus harrisii) insurance metapopulation to assess how alternate management strategies, group housing (mate choice) versus intensive housing (forced monogamy), impact reproductive success and genetic representation. The Tasmanian devil has recently suffered a severe population decline (~85%) after the emergence of two transmissible cancer variants, Devil Facial Tumour 1 (DFT1) and Devil Facial Tumour 2 (DFT2), first observed in 1996 and 2015 respectively [25, 26]. An insurance population was established in 2006 under the management of the Save the Tasmanian Devil Program (STDP) in collaboration with the Zoo and Aquarium Association (ZAA), to provide insurance against extinction, while in situ conservation focused on disease control . Initial breeding in the insurance population was through intensive housing, with free-range enclosures (22 ha) commencing in 2009, group housing (2–3 ha) commencing in 2011 and an island population (120 ha) in 2012 .
We had three overall aims for this study. Firstly, we used pedigree data to determine the influence of intensive housing and group housing on reproductive success and litter size in the insurance population. Secondly, we determined the influence of body weight and age on reproductive success in male Tasmanian devils within group housing. Finally, we used stochastic population models to project population sustainability, loss of genetic diversity and accumulation of inbreeding in intensive housing and group housing over time and explored recommendations to resolve unequal genetic contributions.
The mating system of Tasmanian devil
Tasmanian devils live on average for 5–6 years [28, 29], and potentially up to 8 years in captivity . Both males and females reach sexual maturity between 1 and 2 years of age, and experience a reproductive life of 3–4 years [31, 32]. The annual breeding season occurs between February and June . Females are polygamous, polyestrous annual breeders, which can produce a maximum of four joeys each year, due to the maximum number of teats . Males are polygamous and have been reported to display mate guarding at dens after copulation .
Insurance population management of the Tasmanian devil
The Tasmanian devil insurance metapopulation is managed across a range of breeding facilities including: intensive zoos, group housing enclosures, free-range enclosures, an island population and a fenced peninsula . Breeding recommendations are issued to member institutions annually, indicating which devils should be paired together; not every Tasmanian devil is provided with a breeding opportunity every year. Zoo facilities are intensively managed, where breeding recommendations are most often carried out as single breeding pairs. Group housing (maximum 10 individuals) and free-range facilities (maximum 22 individuals) are less intensively managed. Breeding recommendations in group housing comprise a group of multiple males and females, allowing for mate choice to occur within a selected group of individuals . The island and fenced peninsula populations are monitored annually with population supplementation and removal undertaken as required for genetic and demographic sustainability, as well as source populations for wild release. In sites where breeding can be managed it is done using a traditional species management approach of minimising mean kinship  based on pedigree analysis, which is supplemented with molecular data.
In this study, we compared reproductive skew (a measure of unequal genetic contribution) and litter size (offspring survival to weaning, a measure of offspring fitness) between group and intensive housing. Group housing enclosures accommodate approximately 2.67 Tasmanian devils per hectare, with even sex ratios. The enclosures are designed to resemble the natural environment of the Tasmanian devil and provide multiple den sites, crucial for Tasmanian devil reproduction. The Tasmanian devils are housed together all year round, within recommended breeding groups, from December through to November. This allows devils to familiarize themselves with both the enclosure and the other Tasmanian devils prior to breeding season. Genetic parentage of offspring produced in group housing enclosures between 2011 and 2014 was previously determined using a panel of 33 polymorphic microsatellite markers . Parentage data for the 2015 cohort of joeys was determined herein following methods presented in ref  where all sexually mature Tasmanian devils residing in the group housing facility during 2015 and all joeys from the 2015 cohort were genotyped using 33 microsatellites and analysed in CERVUS  to determine parentage. Intensive facilities vary in enclosure size, with an average of 1–2 devils per enclosure. Females with breeding recommendations are housed independently throughout the year and potential mates are introduced to her enclosure once signs of estrous are evident. A female devil may have multiple males introduced throughout the year (in accordance with estrous cycles) if mating attempts are unsuccessful. The current dataset is limited to those females for which introductions were limited to one male at a time, but includes females that had males introduced to them sequentially. Parentage data for intensive facilities was determined from studbook records. The Zoo and Aquarium Association provide breeding recommendations to the facilities on an annual basis. All breeding recommendations are based on a minimisation of mean kinship strategy using studbook pedigree information.
Intensive (forced monogamy) vs. group housing (mate choice)
To investigate the effect of facility type (group housing or intensive housing) on reproductive skew and litter size we compared these outcomes between facilities (both defined below). Data from 2011 to 2015 were investigated for both housing types (Table 1). We predicted that group housing would have a higher reproductive skew (a smaller percentage of individuals achieving reproductive success) than intensive housing, as a result of mate choice and competition experienced within communal living. We predicted that females from group housing would have greater litter size than females from intensive housing as a result of the benefits of mate choice. We used generalized linear mixed models (GLMMs) to investigate the effect of facility type. All GLMMs were carried out in R  fitted with glmr using the R package lme4 . Male and female models were run separately as forced monogamy and mate-choice can affect sexes differently (for a review see ). In total, three models were run as described below.
Males and females were analysed separately. Individual reproductive success was a binomial response, where “achieved” was if they produced at least one offspring, and “failed” was if they failed to reproduce. The collective reproductive success and failure rate from each facility represents the reproductive skew of the facility. Along with facility (our predictor of interest), both models included age as a continuous predictor variable and year and individual ID as categorical random factors. We chose to include individual ID to account for repeated measures from individuals across the 5-year study period. In addition, the female reproductive skew model included body weight as a continuous predictable variable. This predictor was not included for the male reproductive skew model, as there were two years of the five where body weight data for both intensive and group housing males were not available.
To investigate the impact of facility type on litter size we quantified litter size using a binomial response variable, where the number of events (successes) was the number of offspring produced per litter, and the number of trials was the total number of offspring biologically possible (four due to pouch size ). Along with facility (our predictor of interest), age and body weight were included as continuous predictors and year and individual ID were included as categorical random factors.
Our male reproductive skew data indicated that males, within group housing, did not exhibit equal reproductive success among the males in a shared enclosure. To test whether reproductive success significantly differed from our null hypothesis (all males sire equal proportions of offspring in an enclosure) we performed a series of multinomial exact tests using the R-package XNomial . This approach allows us to examine each enclosure as an independent trial; necessary because the number of males and females varied among enclosures. Over the five-year study period we had 26 group housing enclosures (n-trials = 26). A substantial proportion of enclosures were significantly different from our null hypothesis (see Results).
To further understand what factors may influence reproductive variance we investigated whether standardized body weight determined male reproductive success in group housing enclosures using a GLMM. Our model represented reproductive success as a binomial response, where an individual achieved a success event if they sired at least one offspring, and achieved a fail event if they sired no offspring. The captive environment of group housing violates certain assumptions for male competition. Males exist in closed populations, where each individual male is in competition with only a subset of males, which can vary in body weight range across enclosures, and across years. We wanted to investigate whether male competition to be the superior male was a result of being the largest male per enclosure (each male is compared to the remaining males in said enclosure), or alternatively, if superior males had a body weight threshold (regardless of body weight ranking per enclosure, all individuals below a certain body weight are generally inferior males) (for a review on male body weight and reproductive success in invertebrates, reptiles, birds and mammals see ). To investigate this question our predictor of interest was standardised male body weight within enclosure (denoted “s.bodyweight”) for the three years body weight data was available, calculated as the difference between a males individual body weight and the mean body weight for all males within the enclosure (i.e. s.bodyweight = bodyweight – bodyweightenclosure mean). Age and s.bodyweight were included as our predictors of interest, along with year and individual ID as categorical random factors.
All global models, for all our GLMMs, were standardised to facilitate comparison of parameter estimates across models, using the arm package  in R. Model selection proceeded under information theory (following ref ). From each standardised global model we created a complete subset of models using the MuMIn package  in R. Each submodel was ranked using the Akaike’s information criterion (AICC), and conditional model averaging was used to take all models that fell within 2 AICC of the highest ranked model (see Additional file 1: Tables S2-S5). Inference was based on standardised effect sizes, their standard errors, and the relative importance (RI) scores for each predictor .
Population growth, genetic diversity and inbreeding modeling: Intensive housing (forced monogamy) vs. group housing (mate choice)
To investigate the long-term population differences between intensive housing and group housing we used stochastic population models to project population growth (or sustainability at carrying capacity), gene diversity and inbreeding accumulation over time using VORTEX v10 [45, 46]; all input parameters are shown in Additional file 1: Table S1. Models were run for 50 years, as this is the intended length of the Tasmanian devil insurance population , and a general milestone applied to insurance populations. The initial population size and long-term carrying capacity were equal between the intensive housing model and the group housing model (to allow comparisons over time). The carrying capacity was chosen to reflect a breeding age distribution similar to that maintained by conservation breeding programs broadly, in general terms more younger individuals than older animals and equal sex ratio. Female reproductive parameters were determined from our GLMM results, as opposed to directly applying our empirical data. We did this for two reasons; 1) female reproductive success was not influenced by facility type, and using our modeled results indicated for a single value for females across both management strategies, and 2) the modeling data summarises our dataset as a whole, providing a more accurate estimate of mean effects than data from any given subset (such as year). Mortality rates were set as equal between facilities, to enable strict comparison of reproductive processes alone, and were derived from direct census data from the Tasmanian devil insurance population. As age was found to be a significant predictor of reproductive success for female Tasmanian devils , the reproductive rate for females was entered as a function of age (Additional file 1: Table S1). For a baseline model, male reproductive rate had a 0% skew in both facilities (i.e. 100% of the males were contributing to the breeding pool), with all remaining parameters equal excluding mating opportunities (Additional file 1: Table S1), as this reflects the key difference between group housing and intensive housing. All models were run for 1000 iterations and average results reported.
To investigate how reproductive skews in intensive and group housing influences population growth (or sustainability at carrying capacity), gene diversity and inbreeding accumulation we ran a series of simulation models with gradually increasing reproductive skew. Eight models were run for each enclosure type, whereby reproductive skew increased at 10% increments from 0 to 70% (e.g. a model run with a 70% male reproductive skew operates with only 30% of the male population contributing to the breeding pool in a given year).
Finally, we investigated the impact of transfers between our two management types, in the face of average reproductive skews. All parameters were set to those listed in Additional file 1: Table S1, with the exception of male reproductive skews being set to the average observed percent for each enclosure type (58.12% for group housing and 43.90% for intensive housing). Eleven models were run; a baseline model with the aforementioned male reproductive skews for each facility and zero transfers; transferring 1, 3 or 5 males every year between the two facilities; transferring 1, 3, or 5 females every year between the two facilities; and transferring 3 or 5 males or females every 3 or 5 years between the two facilities. Male and female transfers were investigated separately because females were predicted to have a greater probability than males of contributing to the breeding pool in the presence of male reproductive skews. Our models were built so that an individual could only be transferred to its non-origin facility, and could not be transferred back in subsequent years. Raw data for the analyses presented in this paper are available in Additional file 2 (females) and Additional file 3 (males).
In total, empirical data from five mating seasons were analyzed (2011 to 2015) representing 69 breeding recommendations within intensive enclosures, of which 63.67% of females and 56.10% of males were successful (produced at least 1 offspring); and 105 breeding recommendations within group housing enclosures, of which 46.76% of females and 41.88% of males were successful. In group housing 19.17% of successful males were polygynous and 14.82% of successful females polyandrous (Table 1). Females did not have more than two mates per season, and only one male had three mates in a single breeding season, with the remaining polygynous males all having two mates. There were three occurrences of polygyny within the intensive facilities, when one male was given the opportunity to reproduce with two females on separate occasions (7.33% of mated males).
For males, reproductive skew ranged from 53.85% to 61.90% in group housing and 27.27% to 62.50% in intensive housing (Table 1). Our models provided poor support for an effect of group housing on male reproductive skew (Table 2). Female reproductive skew ranged from 43.48% to 63.64% in group housing and 21.43% to 53.33% in intensive housing (Table 2). Modeling female reproductive skew showed that age was the most accurate predictor of female reproductive success (Table 2), with a strong negative effect of age on reproductive output. Enclosure type (group or intensive) and body weight were present in models <2AIC of the top model, but were poorly supported as predictors of female reproductive skew (low RI and imprecise estimates with high standard error; Table 2).
There was no significant difference in the number of offspring produced in the group housing enclosures (2.60 ± 1.07 SD joeys/female) when compared to the intensive enclosures (2.64 ± 0.99 SD joeys/female) (Tables 1 and 2). Age appeared within 2 AIC of the top model for litter size, although evidence for this effect was poor (very low RI and high error; Table 2).
Of the 26 group housing enclosures, each enclosure had a minimum of one female reproduce. We found statistically significant deviations from equal siring success for males in 11 of the 26 group enclosures (Table 3). This effect was observed at a higher frequency when 4 joeys (or more) were produced per an enclosure, occurring in 11 out of 16 trials (i.e. a statistically significant deviation from uniform representation at α = 0.05 was seen in 68.75% of enclosures where at least 4 joeys were produced). In 10 enclosures where <4 joeys were produced, none deviated significantly from the null hypothesis of uniform representation (Table 3). Trials with low female reproductive success (small number of total joeys produced) likely limits our ability to detect male competition.
After model averaging, standardised body weight appeared in the top model for male Tasmanian devils in group housing, showing a strong positive effect of standardised weight on reproductive success (Table 2). Standard error was small (confidence intervals do not encompass zero) and the relative importance value was high (1.00; Table 2). Larger males within an enclosure were significantly more likely to achieve reproductive success than smaller males. Age appeared in the top model or a model within 2 AIC of the top model for male Tasmanian devil in group housing, but was poorly supported (Table 2).
Population growth, genetic diversity and inbreeding modeling: Intensive (forced monogamy) vs. group housing (mate choice)
Models, using Vortex, showed that after 50 years, in the absence of male reproductive skew (i.e. 100% of the males were contributing to the breeding population), population growth, gene diversity and inbreeding accumulation were very similar between group housing and intensive housing. The population size fluctuated mildly throughout the 50 year time period for both populations; group housing averaged 88.14 Tasmanian devils after 50 years, and intensive housing averaged 88.06 Tasmanian devils after 50 years. Gene diversity was similar between the two facilities, with group housing predicted to maintain 84.95% gene diversity and intensive housing predicted to maintain 85.07% gene diversity. Inbreeding coefficients were also similar, with group housing predicted to have a mean inbreeding coefficient of 0.13 after 50 years and intensive housing predicted to have a mean inbreeding coefficient of 0.14 after 50 years.
As male reproductive skew increased, we saw a decrease in genetic diversity and an increase in inbreeding accumulation, over 50 years, in both enclosure types (Fig. 1). Loss of gene diversity and inbreeding accumulation was most prominent for intensive enclosures when male reproductive skew was greater than 60% (i.e. less than 40% of the males were successfully reproducing). The results most reflective of the Tasmanian devil insurance population are displayed in Fig. 2. Group housing had a male reproductive skew of 58.12% and intensive housing had a male reproductive skew of 43.90%. Despite minimal differences existing between the two facilities under 0% male reproductive skew, observed levels of skew reveal a pattern whereby intensive housing accumulated less inbreeding and maintained a degree of higher gene diversity over the 50 year period than group housing (Fig. 2).
Simulated transfers of both male and female Tasmanian devils between group housing and intensive housing led to predicted increases in gene diversity and decreases in inbreeding accumulation. The change in gene diversity (increase) and inbreeding (decrease) was greatest when comparing the predicted outcome between zero transfers and one (male or female) annual transfer (Table 4). It should be noted that transferring females gave consistently greater increases in gene diversity and greater decreases in inbreeding accumulation relative to transferring males (Table 4).
For some critically endangered species, such as the giant panda  and cheetah , providing individuals with the opportunity for mate choice has led to successful reproduction in captive programs. However, here we have shown that an understanding of the impact of different housing on retention of genetic diversity is essential to the long-term management of any insurance population. As the global biodiversity crisis deepens and more species are under threat, conservation breeding programs and insurance populations are seen as viable solutions to species survival. Severe reproductive skews during the establishment of an insurance population will have long spanning consequences. These skews can arise as a result of forced monogamy, mate incompatibilities, intrasexual competition, or a combination of these. Teasing apart the mechanisms causing reproductive skews, and their consequences, in small populations is crucial for conservation management, particularly for genetically depauperate species such as Tasmanian devils . Ours is one of only a handful of studies that examines these effects empirically in a long-running, intensively managed population. We show that genetic diversity in the Tasmanian devil insurance metapopulation is retained more effectively in intensive housing (forced monogamy) when compared to group ( max 10 individuals) housing (mate choice). Furthermore, we show that purported benefits of mate choice, such as increased offspring survival, were not realized as both group and intensively housed animals showed similar litter sizes at weaning. This finding is in line with other work across the entire Tasmanian devil insurance population that showed intensive and group housing facilities produced on average 2.51 ± 0.23 joeys/female compared with larger free-range enclosures (3.05 ± 0.41 joeys/female) and Maria Island (3.10 ± 0.21 joeys/female) .
As part of the insurance metapopulation strategy Tasmanian devils have been housed in group scenarios to provide better mate choice and/or competition as well as maintain individuals in a more “wild-type” environment. The agile antechinus (Antechinus agilis), from the same Dasyurid family as the Tasmanian devil, has been shown to preferentially mate and produce young with genetically dissimilar males when provided with a range of potential mates [49, 50]. Chemosensory cues have been shown to play an influential role in this antechinus mate-choice process . Similar to the antechinus, the Tasmanian devil has been observed to regularly scent mark . The purpose of devil scent marking has not been investigated, but it is plausible that it plays a similar role in territory marking and may influence female mate choice. In theory, when provided with mate choice opportunities, female Tasmanian devils may favour dominant males (if they have engaged in more scent marking) or the most genetically diverse male (potentially determined via olfactory cues). However, the correlation between neutral diversity and male reproductive success has been previously investigated in group housing facilities with no significant influence being detected .
Providing Tasmanian devils with mate choice opportunities has resulted in a reproductive skew in the insurance population. Only 41.88% of all possible males have contributed genetically over the five-year period. Further, there was no fitness benefit in the productivity of offspring between group and intensive housing. The Tasmanian devil insurance metapopulation has been very successful to date with reproductive goals surpassed annually  and demographic sustainability achieved. However, allowing group housing for the purpose of mate choice in this species is not an optimal long-term genetic management strategy, as reproductive skew will hinder the maintenance of genetic diversity.
Due to the short reproductive life (1–3 years) of Tasmanian devils, the unknown parentage of offspring produced in group housing has meant annual breeding recommendations are issued before the mean kinship of any given year’s offspring can be determined through genetic analyses . This is because annual breeding recommendations are issued in September each year to permit the transfer of devils before the heat of the austral summer and breeding season the following March. Joeys wean from September through to December with parentage analysis only able to be completed by the following January/February. Whilst subtle departures from equal breeding success are expected to occur naturally in group enclosures (and in the wild), the results presented herein suggest male competition within enclosures creates significantly unequal reproductive success. Of the 26 enclosures examined, 11 had one male siring a significantly greater number of offspring (Table 3). In addition, five of these 11 enclosures had only one male reproducing and the remaining six enclosures had only two males reproducing. This means that each year as many as 50% - 80% of males housed in group enclosures and provided with the opportunity to breed produce no offspring. This has led to a reproductive skew in the insurance population potentially affecting the long-term retention of genetic diversity. Fortunately group housing has only been used in the insurance metapopulation since 2011 and so knowledge of this reproductive skew can now be accounted for in future breeding recommendations.
To better understand why only a small portion of available males were siring offspring we investigated the role of body weight on male reproductive success. Regardless of an individual’s age we found that large males are most likely to contribute the greatest proportion of joeys to the annual joey cohort. As the breeding yards are adjusted each year, we suggest that predetermined residency (or territory) may not be the greatest driver of reproductive success, but rather male body size. For many insurance populations experiencing difficulties with maintaining equal genetic representation, unavoidable overrepresentation of more manageable breeding pairs can play a role (e.g. mate aggression in the black rhinoceros (Diceros bicornis) inadvertently led to more frequent breeding opportunities between more timid or adapted individuals ). In other endangered species that experience reproductive skews and pair incompatibilities, such as pygmy rabbits (Brachylagus idahoensis) and clouded leopards (Neofelis nebulosa), long-term exposure to a genetically desirable breeding partner, prior to breeding season, increased reproductive success in both species [16, 54]. Reproductive dominance in males has been observed in the endangered gopher tortoise (Gopherus polyphemus) where larger males achieve the greatest reproductive output annually in a heavily monitored rewilding site . In the Tasmanian devil, the annual breeding recommendations for group housing enclosures provide for a staggered age structure, as housing individuals of the same age together for breeding may lead to over-aggressive behaviours (DPIPWE/ZAA 2013; STDP pers. comm.). Given body weight, rather than age, is a greater determinant of reproductive success in group housed male Tasmanian devils, different management strategies will be required to balance genetic representation in the future. For male Tasmanian devils within intensive facilities, both age and body weight had negligible effect on individual reproductive success. Further analysis into the potential reproductive skew in free-range enclosures (1 devil/ha) and at the island site (0.5–0.8 devils/ha) will provide a better understanding of the interplay of weight and density on the reproductive success of male Tasmanian devils.
The overarching aim of the Tasmanian devil insurance population is to maintain at least 95% of wild-sourced genetic diversity for 50 years . By modelling varying reproductive success of male Tasmanian devils in both group and intensive housing we were able to assess how these differing approaches would sustain population size, retain gene diversity and accumulate inbreeding over the 50 years. There was a gradual increase of inbreeding accumulation in group housing (relative to intensive) and gradual decrease in gene diversity in group housing (relative to intensive) over 50 years (Fig. 2). However, as the Tasmanian devil insurance population operates as a metapopulation we modelled the impact transfers would have on the gene retention and inbreeding accumulation under differing male reproductive skews. In the face of “average” male reproductive skews within both facilities, the transfer of a single male or female Tasmanian devil annually restored a significant degree of genetic diversity to the population and greatly decreased the accumulation of inbreeding within the population (Table 4). This demonstrates the value and necessity of cooperative management for insurance populations.
The concept of group housing in conservation breeding programs is becoming increasingly topical, providing many behavioural and welfare advantages, in addition for mate choice opportunities. However, management teams should be cautious of using solely a group housing approach, as this can potentially obstruct the long-term genetic maintenance goal for the species. Captive breeding programs traditionally aim to maintain diversity by identifying the breeding combinations that would maximise genetic diversity of the population, and then housing those breeders together, a type of forced monogamy, regardless of species life history. Our study is one of the few studies to have tested the immediate and long-term consequences of providing mate choice opportunities in an insurance program. We found that the purported benefits of these mate choice opportunities were not realised in Tasmanian devils. Larger litter sizes did not occur in group housing, and there was a high degree of reproductive skew detected in group enclosures raising serious concerns about maintenance of genetic diversity and accumulation of inbreeding. Larger males were significantly more likely to achieve reproductive success in group enclosures compared to smaller males, suggesting intrasexual competition can influence group housing outcomes. Using computational simulations, our population projections reinforced these conclusions, and showed that negative genetic effects can be mitigated through frequent transfers among sites. Group housing during the breeding season is appropriate for many species of conservation concern for welfare and other reasons; we show that acute awareness of reproductive skew, and what is the potential causative for this skew, is essential to the effective use of group housing in conservation breeding programs.
Akaike’s Information Criterion
Devil Facial Tumour 1
Devil Facial Tumour 2
Generalised Linear Mixed Model
Save the Tasmanian Devil Program
Zoo and Aquarium Association
Conde DA, Flesness N, Colchero F, Jones OR, Scheuerlein A. An emerging role of zoos to conserve biodiversity. Science. 2011;331:1390–1.
Hoffmann M, Hilton-Taylor C, Angulo A, Böhm M, Brooks TM, Butchart SHM, Carpenter KE, Chanson J, Collen B, Cox NA, et al. The impact of conservation on the status of the world’s vertebrates. Science. 2010;330:1503–9.
Toone WD, Wallace MP. The extinction in the wild and reintroduction of the California condor (Gymnogyps californianus). In: Olney PJS, Mace GM, Feistner ATC editors. Creative Conservation. Springer Dordrecht; 1994. p. 411–419
Kierulff M, Ruiz-Miranda C, Oliveira PP, Beck B, Martins A, Dietz J, Rambaldi D, Baker A. The golden lion tamarin Leontopithecus rosalia: a conservation success story. Int Zoo Yearb. 2012;46(1):36–45.
Santymire RM, Livieri TM, Branvold-Faber H, Marinari PE. The Black-Footed Ferret: On the Brink of Recovery? In: Holt WV, Brown JL, Comizzoli P, editors. Reproductive Sciences in Animal Conservation: Progress and Prospects. New York, NY: Springer New York; 2014. p. 119–34.
Ebenhard T. Conservation breeding as a tool for saving animal species from extinction. Trends Ecol Evol. 1995;10(11):438–43.
Ralls K, Ballou J. Captive breeding programs for populations with a small number of founders. Trends Ecol Evol. 1986;1(1):19–22.
Fa JE, Funk SM. O'Connell D. Zoo Conservation Biology. UK: Cambridge University Press; 2011.
Frankham R, Ballou JD, Briscoe DA. Introduction to conservation genetics (2nd edition), 2nd edition edn. Cambridge: UK Cambridge University Press; 2010.
Lacy RC. Analysis of founder representation in pedigrees: founder equivalents and founder genome equivalents. Zoo Biol. 1989;8:111–23.
Ballou JD, Lacy RC. Identifying genetically important individuals for management of genetic diversity in pedigreed populations. In: Ballou JD, Gilpin M, Foose TJ, editors. Population management for survival and recovery. New York, NY: Columbia University Press; 1995. p. 76–111.
Ballou JD, Lees CM, Faust LJ, Long S, Lynch C, Bingaman-Lackey L, Foose TJ. Demographic and genetic management of captive populations. Wild mammals in captivity. Chicago, London: University of Chicago; 2010.
Martin-Wintle MS, Shepherdson D, Zhang G, Zhang H, Li D, Zhou X, Li R, Swaisgood RR. Free mate choice enhances conservation breeding in the endangered giant panda. Nat Commun. 2015;6:10125.
Asa CS, Traylor-Holzer K, Lacy RC. Can conservation-breeding programmes be improved by incorporating mate choice? Int Zoo Yearb. 2011;45(1):203–12.
Schulte-Hostedde AI, Mastromonaco GF. Integrating evolution in the management of captive zoo populations. Evol Appl. 2015;8(5):413–22.
Yamada JK, Durrant BS. Reproductive parameters of clouded leopards (Neofelis nebulosa). Zoo Biol. 1989;8:223–31.
Fuerst P, Maruyama T. Considerations on the conservation of alleles and of genic heterozygosity in small managed populations. Zoo Biology. 1986;5(2):171–9.
Lacy R. Loss of genetic diversity from managed populations: interacting effects of drift, mutation, immigration, selection, and population subdivision. Conserv Biol. 1987;1(2):143–58.
Neff B, Pitcher T. Genetic quality and sexual selection: an integrated framework for good genes and compatible genes. Mol Ecol. 2005;14:19–38.
Ihle M, Kempenaers B, Forstmeier W. Fitness benefits of mate choice for compatibility in a socially monogamous species. PLoS Biol. 2015;13(9):e1002248.
Batkay D. Mate selection preferences of captive female cheetahs (Acinonyx jubatus): York University; 2005.
Miller KA, Nelson NJ, Smith HG, Moore JA. How do reproductive skew and founder group size affect genetic diversity in reintroduced populations? Mol Ecol. 2009;18(18):3792–802.
Kokko H, Rankin DJ. Lonely hearts or sex in the city? Density-dependent effects in mating systems. Philos Trans Royal Soc B. 2006;361(1466):319–34.
Hunt J, Breuker CJ, Sadowski JA, Moore AJ. Male-male competition, female mate choice and their interaction: determining total sexual selection. J Evol Biol. 2009;22(1):13–26.
Loh R, Bergfeld J, Hayes D, O’Hara A, Pyecroft S, Raidal S. The pathology of devil facial tumor disease (DFTD) in Tasmanian devils (Sarcophilus harrisii). Vet Pathol. 2006;43(6):890–5.
Pye RJ, Pemberton D, Tovar C, Tubio JMC, Dunn KA, Fox S, Darby J, Hayes D, Knowles GW, Kreiss A, et al. A second transmissible cancer in Tasmanian devils. Proc Natl Acad Sci U S A. 2016;113(2):374–9.
Hogg CJ, Lee AV, Srb C, Hibbard C. Metapopulation management of an endangered species with limited genetic diversity in the presence of disease: the Tasmanian devil Sarcophilus harrisii. Int Zoo Yearb. 2016;51:1–17.
Guiler ER. Observations on the Tasmanian devil, Sarcophilus harrisii (Marsupialia: Dasyuridae). I. Numbers, home range, movements, and food in two populations. Aust J Zool. 1970;18:49–62.
Guiler ER. Observations on the Tasmanian devil, Sarcophilus harrisii (Dasyuridae: Marsupialia) at Granville harbour, 1966-75. Pap Proc R Soc Tasmania. 1978;112:161–88.
Srb C: Tasmanian devil studbook, Healesville sanctuary on behalf of the Zoo and Aquarium Association. 2016.
Guiler E. Observations on the Tasmanian devil, Sarcophilus harrisii (Marsupialia: Dasyuridae) II. Reproduction, breeding and growth of pouch young. Aust J Zool. 1970;18(1):63–70.
Keeley T, O'Brien JK, Fanson BG, Masters K, McGreevy PD. The reproductive cycle of the Tasmanian devil (Sarcophilus harrisii) and factors associated with reproductive success in captivity. Gen Comp Endocrinol. 2012;176(2):182–91.
Keeley T, Russell T, Carmody K, Kirk G, Eastley T, Britt-Lewis A, Post M, Burridge M, Eccleston S, Faulkner T, et al. Seasonality and breeding success of captive and wild Tasmanian devils (Sarcophilus harrisii). Theriogenology. 2017;95:33–41.
Hughes RL. Reproduction in the Tasmanian devil Sarcophilus harrisii (Dasyuridae, Marsupialia). In: Archer M, editor. Carnivorous Marsupials. Sydney: Surrey Beatty and Sons; 1982. p. 49–63.
Pemberton D. Social organisation and behaviour of the Tasmanian devil, Sarcophilus harrisii: University of Tasmania; 1990.
Gooley R, Hogg C, Belov K, Grueber C. No evidence of inbreeding depression in a Tasmanian devil insurance population despite significant variation in inbreeding. Sci Rep. 2017;7(1):1830
Kalinowski S, Taper M, Marshall T. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Mol Ecol. 2007;16(5):1099–106.
R Core Team: R: A language and environment for statistical computing. Vienna. Available: http://www.r-project.org/. 2015.
Bates D, Maechler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67(1):1–48.
Edward DA. The description of mate choice. Behav Ecol. 2014;26(2):301–10.
Edwards K, Walker S, Dunham A, Pilgrim M, Okita-Ouma B, Shultz S. Low birth rates and reproductive skew limit the viability of Europe’s captive eastern black rhinoceros, Diceros bicornis michaeli. Biodivers Conserv. 2015;24:2831–52.
Gelman A, Su Y-S, Yajima M, Hill J, Pittau MG, Kerman J, Zheng T, Dorie V: arm: Data analysis using regression and multilevel/hierarchical models. http://CRAN.R-project.org/package=arm 2015.
Bartoń K: MuMIn: Multi-model inference. https://CRAN.R-project.org/package=MuMIn 2016.
Grueber CE, Nakagawa S, Laws RJ, Jamieson IG. Multimodel inference in ecology and evolution: challenges and solutions. J Evol Biol. 2011;24(4):699–711.
Lacy RC. VORTEX: a computer simulation model for population viability analysis. Wildl Res. 1993;20:45–65.
Lacy RC. Structure of the VORTEX simulation model for population viability analysis. Ecol Bull. 2000;48:191–203.
CBSG. Tasmanian devil PHVA final report. In: IUCN/SSC Conservation Breeding Specialist Group: Apple Valley, MN; 2008.
Miller W, Hayes VM, Ratan A, Petersen DC, Wittekindt NE, Miller J, Walenz B, Knight J, Qi J, Zhao FQ, et al. Genetic diversity and population structure of the endangered marsupial Sarcophilus harrisii (Tasmanian devil). Proc Natl Acad Sci U S A. 2012;109(45):18625–5.
Parrott ML, Ward SJ, Temple-Smith PD. Olfactory cues, genetic relatedness and female mate choice in the agile antechinus (Antechinus agilis). Behav Ecol Sociobiol. 2007;61(7):1075–9.
Parrott ML, Ward SJ, Temple-Smith PD, Selwood L. Females choose mates based on genetic relatedness in a small Dasyurid marsupial, the agile antechinus (Antechinus agilis). PLoS One. 2015;10(4):e0122381.
Lazenby BT, Dickman CR. Patterns of detection and capture are associated with cohabiting predators and prey. PLoS One. 2013;8(4):e59846.
Hogg CJ, Ivy JA, Srb C, Hockley J, Lees C. Influence of genetic provenance and birth origin on productivity of the Tasmanian devil insurance population. Conserv Genet. 2015;16:1465–73.
Edwards KL, Walker SL, Dunham AE, Pilgrim M, Okita-Ouma B, Shultz S. Low birth rates and reproductive skew limit the viability of Europe’s captive eastern black rhinoceros, Diceros bicornis michaeli. Biodivers Conserv. 2015;24(11):2831–52.
Martin MS, Shepherdson DJ. Role of familiarity and preference in reproductive success in ex situ breeding programs. Conserv Biol. 2012;26(4):649–56.
Tuberville TD, Norton TM, Waffa BJ, Hagen C, Glenn TC. Mating system in a gopher tortoise population established through multiple translocations: apparent advantage of prior residence. Biol Conserv. 2011;144(1):175–83.
We thank John Weigel, Tim Faulkner, Liz Vella, Dean Reid, Carla Srb, Kathy Starr, Androo Kelly, Paul Ralley, Anthea Guisenburg and all the other devil keepers at the 37 zoos involved in the Tasmanian devil insurance population. This work would not have been possible without their hard work and dedication to the insurance program.
This research was supported by funding from the Australian Research Council, grant numbers LP140100508 (KB, CEG, CJH), and DP170101253 (CEG, KB), the Save the Tasmanian Devil Program, the Zoo and Aquarium Association Australasia and San Diego Zoo Global.
Availability of data and materials
All data generated or analysed during this study are included in this published article [and its Additional files].
Sampling for this study occurred under approval from University of Sydney Ethics Committee #2014/550.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Vortex parameters for the Tasmanian devil insurance population across two housing types; group housing and intensive housing. Table S2. Top model set (top 2AICc) of generalised linear mixed models for male reproductive skew. Table S3. Top model set (top 2AICc) of generalised linear mixed models for female reproductive skew. Table S4. Top model set (top 2AICc) of generalised linear mixed models for litter size. Table S5. Top model set (top 2AICc) of generalised linear mixed models for group housing male reproductive success. (DOCX 16 kb)
Raw Data Female Tasmanian Devils. (CSV 4 kb)
Raw Data Male Tasmanian Devils. (CSV 5 kb)