Abstract
Devil Facial Tumour Disease (DFTD), a highly contagious cancer, has decimated Tasmanian devil (Sarcophilus harrisii) numbers in the wild. To ensure its long-term survival, a captive breeding program was implemented but has not been as successful as envisaged at its launch in 2005. We therefore investigated the reproductive success of 65 captive devil pair combinations, of which 35 produced offspring (successful pairs) whereas the remaining 30 pairs, despite being observed mating, produced no offspring (unsuccessful pairs). The devils were screened at six MHC Class I-linked microsatellite loci. Our analyses revealed that younger females had a higher probability of being successful than older females. In the successful pairs we also observed a higher difference in total number of heterozygous loci, i.e. when one devil had a high total number of heterozygous loci, its partner had low numbers. Our results therefore suggest that devil reproductive success is subject to disruptive MHC selection, which to our knowledge has never been recorded in any vertebrate. In order to enhance the success of the captive breeding program the results from the present study show the importance of using young (2-year old) females as well as subjecting the devils to MHC genotyping.
Similar content being viewed by others
Introduction
The major histocompatibility complex (MHC) has been shown to be one of the key molecular determinants of mate choice in numerous vertebrates including humans1,2,3,4,5,6,7. Class I and II MHC molecules are responsible for the processing and presentation of intra- and extra-cellular peptides derived from invading pathogens to cytotoxic T cells and helper T cells and, hence, constitute a crucial part of the vertebrate immune system8,9. Due to the ability to recognize and present peptides from a wide array of rapidly evolving pathogens, the MHC encompasses the most variable set of genes with heterozygosity values exceeding those predicted by neutrality10.
Pathogen-driven selection has been documented to result in two, not mutually exclusive, MHC responses: (i) selection for specific MHC alleles11,12,13,14 and/or (ii) selection for enhanced MHC polymorphism13,14,15,16. Although increased levels of MHC polymorphism enable wider recognition of pathogens, it might also lead to an inability to eliminate T cells reacting with self-peptide-MHC combinations17,18. Consequently, diversifying MHC selection might be counteracted by the deleterious effects caused by autoimmunity19. Indeed, theoretical models as well as empirical studies have shown that optimal pathogen resistance often occurs at an optimal intermediate level of MHC polymorphism13,20,21,22,23,24,25,26,27.
MHC has also been shown to be involved in individual mate choice and providing offspring with indirect genetic benefits in at least three ways via (i) acquisition of “good genes” i.e. genetic elements that contribute to lifetime reproductive success regardless of an individual’s additional genotype28, and/or (ii) acquiring optimal genetic compatibility29, and/or (iii) achieving enhanced genetic diversity30.
Importantly, significant reduction in MHC diversity may not only affect resistance to pathogens but has also been shown to result in an increased risk of extinction due to inbreeding depression31. Moreover, reduced MHC polymorphism has been suggested to have contributed to the emergence of a clonally transmissible cancer, Devil Facial Tumour Disease (DFTD) in the world’s largest living marsupial carnivore: the Tasmanian devil (Sarcophilus harrisii)32. DFTD was first reported in north-eastern Tasmania in 1996 but since then has caused severe reductions in devil numbers (>70%) questioning the long-term survival of this iconic species33,34. This devastating disease is spread among devils via biting during social interactions35. The devil’s immune system is unable to mount an effective immune response to the tumour as DFTD cells are able to avoid immune recognition by down-regulating MHC expression36,37. Metabolic failure, tumour related cachexia and metastases result in devil death within 6 to 9 months of the emergence of the first lesions38. In order to ensure that the Tasmanian devil will not face extinction a large scale captive breeding program commenced in 200539, but overall breeding success is still not optimal40,41. In the present study we therefore explore how traits such as devil age and MHC polymorphism affected reproductive success in the captive devil population.
Results
Of a total of 65 captive devil pair combinations, 35 produced offspring whereas the remaining 30 pair combinations did not produce any offspring (see Table 1 for a detailed description of the pair combinations). The genetic diversity analyses did not reveal any allele dropouts or null alleles. The full logistic mixed model revealed that male age, identification number, number of female heterozygous loci, number of male heterozygous loci and number of similar alleles did not affect devil reproductive success (P > 0.20), and those predictors were therefore backwards-eliminated. The final model revealed that younger females had a higher probability of being successful than older females as indicated by a negative regression parameter estimate for female age = −0.1068 ± 0.04, SE, p = 0.012 (Table 2). In fact, the final analysis predicted a decrease in the chance of producing offspring from 0.68 (0.82–0.51, 95% confidence interval) for females at the age of 24 months to 0.19 (0.5–0.05) for females older than 60 months (Fig. 1). As mentioned above, male age did not affect reproductive success. Interestingly, in successful pairs we observed a higher difference in total number of heterozygous loci, i.e. when one devil had a high total number of heterozygous loci, its partner had a low number, as evident by a positive regression parameter estimate (0.8202 ± 0.3329, p = 0.0165). Thus, pairs with opposing total numbers of heterozygous loci were found to have a higher probability of being successful reproducers than pairs with similar numbers (Fig. 2).
Relationship between female age and the success of reproduction in Tasmanian devil pairs (a). Points with error bars depict reproductive success (±95% CI) of pairs across female age. Grey line shows the binomial model prediction of this relationship with the 95% confidence interval (grey polygon).
(a) Mode prediction of the interaction between the number of male and female heterozygous loci. The color scale (within the data range) and contour lines (across the entire range of female and male heterozygous loci) indicate the probability of being successful according to the combination of the number of male and female heterozygous loci. (b) Diagonal section of the model prediction shown in (a) across the optimal pair combinations with 95% confidence interval.
Discussion
Prior to the emergence of DFTD female devils in the wild commenced reproduction at an age of two years and thereafter produced offspring for the next three years, becoming senescent at an age of five to six years33. The results from the present study, however, show that a decline in captive female reproductive success may already occur at an age of three years suggesting that some captive females might be affected by reproductive senescence at an earlier age than that recorded in the wild.
Previous studies have shown that devils are subjected to low genetic diversity at both neutral (microsatellite) and coding genomic regions at the MHC42,43,44,45. The low level of genetic polymorphism has been suggested to be caused by population bottlenecks during the Pleistocene and Holocene46,47 and by pathogens such as a canine-distemper-like disease in the early twentieth century44. The genetic diversity of the captive devils used in the present study was similar to that recorded in devils captured in the wild48. As we did not observe any overall significant difference in genetic diversity among the successful and the unsuccessful pair combinations we find it unlikely that the difference in reproductive success was caused by a concomitant discrepancy in genetic diversity.
In vertebrates such as sand lizards (Lacerta agilis)49, savannah sparrows (Passerculus sandwichensis)50, southern dunlins (Calidris alpina schinzii)51, great tits (Parus major)52 and the cynomolgus macaque (Macaca fascicularis)53 reproductive success has been negatively linked to male - female genetic similarity. However, the results from the present study suggest that within-pair genetic similarity did not affect devil reproductive success.
As mentioned above, both theoretical and empirical studies have shown that selection for an intermediate and optimal number of MHC alleles may result in both increased reproductive success and immune function13,20,21,22,23,24,25,26,27. However, our results show that pairs with higher difference in total number of heterozygous loci had an increased probability of being successful than pairs with similar numbers of heterozygous loci.
If devils have been under selection for an intermediate, optimal number of MHC alleles (e.g., 6 alleles), we would expect reproductive success to be the same in pairs in which partners both have 3 heterozygous loci as in pairs in which one partner has 1 and the other 5 heterozygous loci. Our results therefore suggest that devil MHC may be subject to disruptive selection. Although similar MHC processes have been suggested to drive MHC evolution in allopatric taxa in different habitats54, it has to our knowledge, never been recorded to affect reproductive success.
We propose two, not mutually exclusive processes, to underpin the significant effect of selection for MHC-driven devil reproductive success: pre-copulatory and/or post-copulatory (cryptic) female choice of male sperm. As mentioned above, MHC-based pre-copulatory mate choice, sometimes based on olfactory clues, has been recorded in wide range of vertebrates including humans1,2,3,4,5,6,7,55. MHC-based post-copulatory female cryptic choice of male sperm has also been documented in several vertebrates such as fish56, lizards49, birds57 and mammals58. As female devils in both the successful and unsuccessful pair combinations were observed mating59,60,61,62,63,64,65, we suggest that the significant difference in reproductive success between the two groups might be caused by post-copulatory cryptic female choice. In externally fertilizing fishes, the ovarian fluid released by the female may bias fertilization success towards males with a particular genotype66,67. However, if a similar mechanism may influence male fertilization success in internal fertilizers, such as mammals, is to our knowledge unknown.
Conclusion
Our results suggest that the combination of female age (i.e. younger females having higher reproductive success compared to older females) and MHC diversity constitutes a significant determinant of reproductive success in two groups of devils with different MHC traits and may provide an explanation for the relatively low reproductive success recorded in captive breeding programs. In order to enhance the success of this iconic species, our results advocate (i) the use of young (2 year old) female devils, and (ii) that both male and female devils are subjected to MHC genotyping with pair combinations maximizing the total numbers of heterozygous loci at opposite ends of the heterozygosity scale in order to maximize breeding success.
Material and Methods
Study animals
Data on pairings were obtained from the Tasmanian devil studbook68 and the insurance population annual reports59,60,61,62,63,64,65. Single female devils were paired with a single male when females showed signs of estrus such as a fluid filled neck roll, fat stores in the tail, reduced aggression and loss of appetite69. Consequently, the females were not given a choice of partner. The 37 females and 43 males used in the present study were kept across 11 institutions in Australia and the pairings were conducted from 2007 to 2012. All ear biopsies were collected by veterinarians and trained staff during the devil’s annual exams at the zoological institutions where they were housed and were hence carried out in accordance with relevant guidelines and regulations. The research was conducted under University of Sydney animal ethics approval number N00/8-2011/1/5584.
Genetic analyses
DNA was extracted from ear biopsies using a Qiagen DNeasy® extraction kit and subsequently genotyped at six MHC Class I linked microsatellite loci described by Cheng & Belov70; Sh-MHCI101, Sh-MHCI102, Sh-MHCI105, Sh-MHCI106, Sh-MHCI107, and Sh-MHCI108. For further details on size range and primer sequences see Cheng & Belov70. Polymerase chain reactions (PCR) were carried out using the protocols of Cheng & Belov70 and the PCR products were quality tested on a 2% agarose gel. The samples were analyzed at the Ramaciotti Centre (University of New South Wales, Australia) and microsatellite profiles of the individual devils were subsequently scored using Peak Scanner 2.0 (Applied Biosystems 2012).
Estimates of allele frequency, heterozygosity and linkage disequilibrium
Analyses of dropout and the presence of null alleles were conducted using the software Micro-Checker71. The number of alleles per locus together with observed and expected heterozygosity were estimated using the software program arlequin 3.172. Tests for deviations from Hardy–Weinberg equilibrium and linkage disequilibrium, as well as the analysis of genetic structure and differentiation of the two groups were carried out using arlequin 3.172 (HWE parameters included: number of steps in Markov chain = 1,000,000, number of dememorisation steps = 100,000, number of permutations = 10,000; AMOVA parameters included 999 permutations and 3000 Markov steps). Sequential Bonferroni corrections were subsequently conducted on the Hardy-Weinberg equilibrium and the linkage disequilibrium tests. Within-pair genetic similarity was conducted by recording the number of shared/identical alleles in each of the six loci.
Statistical analyses
Logistic mixed models analyses were performed in SAS 9.4 using Proc GLIMMIX73,74 with institution and female ID number as random effects as 22 of the 36 females were paired on more than one occasion. Proc GLIMMIX uses restricted pseudolikelihoods and the full model included male age and identification number, number of female heterozygous loci, number of male heterozygous loci and number of similar alleles. The final model included institution as random effect and absolute number of different heterozygous loci and female age as fixed effects.
References
Bonneaud, C., Chastel, O., Federici, P., Westerdahl, H. & Sorci, G. Complex MHC-based mate choice in a wild passerinje. Proc. R. Soc. B 273, 1111–1116 (2006).
Griggio, M., Biard, C., Penn, D. J. & Hoi, H. Female house sparrows “count on” male genes: experimental evidence for MHC-dependent mate preference in birds. BMC Evol. Biol. 11, 44 (2011).
Milinski, M. et al. Mate choice decisions of stickleback females predictably modified by MHC peptide ligands. Proc. Natl. Acad. Sci. USA 102, 4414–4418 (2005).
Olsson, M. et al. Major histocompatibility complex and mate choice in sand lizards. Proc. R. Soc. B 270, S254–S256 (2003).
Sin, Y. W. et al. MHC class II-assortative mate choice in European badgers Meles meles. Mol. Ecol. 24, 3138–3150 (2015).
Wedekind, C. & Füri, S. Body odour preferences in men and women: do they aim for specific MHC combinations or simply heterozygosity? Proc. R. Soc. B 264, 1471–1479 (1997).
Yamazaki, K. & Beauchamp, G. K. (2007). Genetic basis for MHC dependent mate choice. Adv. Genet. 59, 129–145 (2007).
Klein, J. Natural history of the Major Histocompatibility Complex. (Wiley, New York, 1986).
Zinkernagel, R. M. Associations between major histocompatibility antigens and susceptibility to disease. Ann. Rev. Microbiol. 33, 201–213 (1979).
Edwards, S. V. & Hedrick, P. W. Evolution and ecology of MHC molecules, from genomics to sexual selection. Trends Ecol. Evol. 13, 305–311 (1998).
Harf, R. & Sommer, S. Association between major histocompatibility complex class II DRB alleles and parasite load in the hairy-footed gerbil, Gerbillus paeba, in southern Kalahari. Mol. Ecol. 14, 85–91 (2005).
Lohm, J. et al. Experimental evidence for majorhistocompatibility complex-allele-specific resistance to a bacterial infection. Proc. R. Soc. B 269, 2029–2033 (2002).
Madsen, T. & Ujvari, B. MHC Class I variation associates with parasite resistance and longevity in tropical pythons. J. Evol. Biol. 19, 1973–1978 (2006).
Westerdahl, H. et al. Association between malaria and MHC genes in a migratory songbird. Proc. R. Soc. B 272, 1511–1518 (2005).
Penn, D. J., Damjanovich, K. & Potts, W. MHC heterozygosity confers a selective advantage against multiple-strain infections. Proc. Natl. Acad. Sci. USA 99, 11260–11264 (2002).
Thursz, M. R., Thomas, H. C., Greenwood, B. M. & Hill, A. V. Heterozygote advantage for HLA class-II type in hepatitis B virus infection. Nature Genet. 17, 11–12 (1997).
Borghans, J. A., Noest, A. J. & De Boer, R. J. Thymic selection does not limit the individual MHC diversity. Eur. J. Immunol. 33, 3353–3358 (2003).
Vidovic, D. & Matzinger, P. Unresponsiveness to a foreign antigen can be caused by self-tolerance. Nature 336, 222–225 (1988).
Fernando, M. M. A. et al. Defining the role of the MHC in autoimmunity, a review and pooled analysis. PLoS Genet. 4, e1000024 (2008).
Aeschlimann, P. B., Häberli, M. A., Reusch, T. B. H., Boehm, T. & Milinski, M. Female sticklebacks Gasterosteus aculeatus use self-reference to optimize MHC allele number during mate selection. Behav. Ecol. Sociobiol. 54, 119–126 (2003).
De Boer, R. J. & Perelson, A. S. How diverse should the immune system be? Proc. R. Soc. B 252, 171–175 (1993).
Kalbe, M. et al. Lifetime reproductive success is maximized with optimal MHC diversity. Proc. R. Soc. B 276, 925–934 (2009).
Milinski, M. The major histocompatibility complex, sexual selection and mate choice. Ann. Rev. Ecol. Syst. 37, 159–186 (2006).
Nowak, M. A., Tarcy-Hornoch, K. & Austyn, J. The optimal number of major histocompatibility complex molecules in an individual. Proc. Natl. Acad. Sci. USA 89, 10896–10899 (1992).
Reusch, T. B. H., Häberli, M. A., Aeschlimann, P. B. & Milinski, M. Female sticklebacks count alleles in a strategy of sexual selection explaining MHC polymorphism. Nature 414, 300–302 (2001).
Wegner, K. M., Kalbe, M., Kurtz, J., Reusch, T. R. B. & Milinski, M. Parasite selection for immunogenic optimality. Science 301, 1343 (2003).
Woelfing, B., Traulsen, A., Milinski, M. & Boehm, T. Does intra-individual major histocompatibility complex diversity keep a golden mean? Phil. Trans. R. Soc. B 364, 117–128 (2009).
Olsson, M. et al. MHC, health, color, and reproductive success in sand lizards. Behav. Ecol. Sociobiol. 58, 289–294 (2005).
Kempenaers, B. Mate choice and genetic quality, a review of the heterozygosity theory. Adv. Study Behav. 37, 189–278 (2007).
Ejsmond, M. J., Radwan, J. & Wilson, A. B. Sexual selection and the evolutionary dynamics of the major histocompatibility complex. Proc. R. Soc. B 281, 20141662 (2014).
Madsen, T., Shine, R., Olsson, M. & Witzell, H. Restoration of an inbred adder population. Nature 402, 34–35 (1999).
Belov, K. The role of the major histocompatibility complex in the spread of contagious cancers. Mammal. Genome 22, 83–90 (2011).
Jones, M. et al. Life-history change in disease-ravaged Tasmanian devil populations. Proc. Natl. Acad. Sci. USA 105, 10023–10027 (2008).
Lachish, S., Miller, K. J., Storfer, A., Goldizen, A. W. & Jones, M. E. Evidence that disease-induced population decline changes genetic structure and alters dispersal patterns in the Tasmanian devil. Heredity 106, 172–182 (2011).
Hamede, R. K., McCallum, H. & Jones, M. Seasonal, demographic and density-related patterns of contact between Tasmanian devils (Sarcophilus harrisii). Implications for transmission of devil facial tumour disease. Austral Ecol. 33, 614–622 (2008).
Siddle, H. V. et al. Reversible epigenetic down-regulation of MHC molecules by devil facial tumour disease illustrates immune escape by a contagious cancer. Proc. Natl. Acad. Sci. USA 110, 5103–5108 (2013).
Karu, N. et al. Discovery of biomarkers for Tasmanian devil cancer (DFTD) by metabolic profiling of serum. J. Proteome Res. 15, 3827–3840 (2016).
Pearse, A.-M. & Swift, K. Allograft theory transmission of devil facial-tumour disease. Nature 439, 549 (2006).
Hogg, C. J. et al. Influence of genetic provenance and birth origin on productivity of the Tasmanian devil insurance population. Conserv. Genet. 16, 1465–1473 (2015).
Hockley, J. DPIPWE/ZAA Husbandry guidelines for Tasmanian devil, Sarcophilus harrisii. Zoo and Aquarium Association, Australia pp 46. (2011).
Keeley, T., O’Brien, J. K., Fanson, B. G., Masters, K. & McGreevy, P. D. The reproductive cycle of the Tasmanian devil Sarcophilus harrisii and factors associated with reproductive success in captivity. Gen. Comp. Endocrinol. 176, 182–191 (2012).
Jones, M. E., Paetkau, D., Geffen, E. & Moritz, C. Genetic diversity and population structure of Tasmanian devils, the largest marsupial carnivore. Mol. Ecol. 13, 2197–2209 (2004).
Epstein, B. et al. Rapid evolutionary response to a transmissible cancer in Tasmanian devils. Nature Commun. 7, 12684 (2016).
Morris, K., Austin, J. J. & Belov, K. Low major histocompatibility complex diversity in the Tasmanian devil predates European settlement and may explain susceptibility to disease epidemics. Biol. Lett. 9, 20120900 (2012).
Morris, K., Wright, B., Grueber, C. E., Hogg, C. & Belov, K. Lack of genetic diversity across diverse immune genes in an endangered mammal, the Tasmanian devil Sarcophilus harrisii. Mol. Ecol. 24, 3860–3872 (2015).
Brown, O. Tasmanian devil Sarcophilus harrisii extinction on the Australian mainland in the Mid- Holocene, multicausality and ENSO intensification. Alcheringa 30, 49–57 (2006).
Bruniche-Olsen, A., Jones, M. E., Austin, J. J., Burridge, C. P. & Holland, B. R. Extensive population decline in the Tasmanian devil predates European settlement and devil facial tumour disease. Biol. Lett. 10, 20140619 (2014).
Lane, A. et al. New Insights into the role of MHC diversity in devil facial tumour disease. PLoS One 7, e36955 (2012).
Olsson, M., Shine, R., Madsen, T., Gullberg, A. & Tegelström, H. Sperm selection by females. Nature 383, 585 (1996).
Freeman-Gallant, C. R., Meguerdichian, M., Wheelwright, N. T. & Sollecito, S. V. Social pairing and female mating fidelity predicted by restriction fragment length polymorphism similarity at the major histocompatibility complex in a songbird. Mol. Ecol. 12, 3077–3083 (2003).
Blomqvist, D., Angela Pauliny, A., Mikael Larsson, M. & Flodin, L.-Å. Trapped in the extinction vortex? Strong genetic effects in a declining vertebrate population. BMC Evol. Biol. 10, 33 (2010).
García-Navas, V., Cáliz-Campal, C., Ferrer, E. S., Sanz, J. J. & Ortego, J. Heterozygosity at a single locus explains a large proportion of variation in two fitness-related traits in great tits: a general or a local effect? J. Evol. Biol. 27, 2807–2819 (2014).
Aarnink, A. et al. Deleterious impact of feto-maternal MHC compatibility on the success of pregnancy in a macaque model. Immunogenetics 66, 105–113 (2014).
Matthews, B., Harmon, L. J., M’Gonigle, L., Marchinko, K. B. & Schaschl, H. Sympatric and allopatric divergence of MHC genes in threespine stickleback. PLoS ONE 5, e10948 (2010).
Kamiya, T., O’Dwyer, K., Westerdahl, H., Senior, A. & Nakagawa, S. A quantitative review of MHC-based mating preference: the role of diversity and dissimilarity. Mol. Ecol. 23, 5151–5163 (2014).
Yeates, S. E. et al. Atlantic salmon eggs favour sperm in competition that have similar major histocompatibility alleles. Proc. R. Soc. B 276, 559–566 (2009).
Løvlie, H., Gillingham, M. A. F., Worley, K., Pizzari, T. & Richardson, D. S. Cryptic female choice favours sperm from major histocompatibility complex-dissimilar males. Proc. R. Soc. B 280, 0131296 (2013).
Schwensow, N., Eberle, M. & Sommer, S. Compatibility counts, MHC-associated mate choice in a wild promiscuous primate. Proc. R. Soc. B 275, 555–564 (2008).
Lees, C. & Srb, C. Annual report for the DPIW-ARAZPA Tasmanian devil insurance population. Zoo and Aquarium Association. pp 27 (2007).
Lees, C. & Srb, C. Annual report for the DPIW-ARAZPA Tasmanian devil insurance population. Zoo and Aquarium Association. pp 49 (2008).
Hibbard, C. J., Ford. C. & Srb, C. Annual report for the DPIW-ZAA Tasmanian Devil insurance population. Zoo and Aquarium Association. pp 57 (2009).
Hibbard, C. J., Hogg, C. J., Srb, C. & Hockley, J. Annual report for the DPIPWE-ZAA Tasmanian devil insurance population. Zoo and Aquarium Association. pp 35 (2010).
Hibbard, C. J., Hogg, C. J., Srb, C., & Hockley, J. Annual report for the DPIPWE-ZAA Tasmanian devil insurance population. Zoo and Aquarium Association. pp 46 (2011).
Hogg, C. J., Srb, C. & Hockley, J. Annual report for the DPIPWE-ZAA Tasmanian devil insurance population. Zoo and Aquarium Association. pp 43 (2012).
Hogg, C. J., Srb, C. & Hockley, J. Annual Report for the DPIPWE-ZAA Tasmanian devil insurance population. Zoo and Aquarium Association. pp 42 (2013).
Alonzo, S. H., Stiver, K. A. & Marsh-Rollo, S. E. (2016). Ovarian fluid allows directional cryptic female choice despite external fertilization. Nature Commun. 7, 12452 (2012).
Gasparini, C. & Pilastro, A. Cryptic female preference for genetically unrelated males is mediated by ovarian fluid in the guppy. Proc. R. Soc. B 278, 2495–2501 (2011).
Srb, C. Australasian regional Tasmanian devil studbook. (Healesville, Victoria, 2014) (2014).
Brice, S. Zoo aquarium association husbandry guidelines for Tasmanian devil, Sarcophilus harrisii. Zoo and Aquarium Association, Australia pp 36 (2006).
Cheng, Y. & Belov, K. Isolation and characterisation of 11 MHC-linked microsatellite loci in the Tasmanian devil Sarcophilus harrisii. Conserv. Genet. Res. 4, 463–465 (2012).
van Oosterhout, C., Hutchinson, W. F., Wills, D. P. M. & Shipley, P. Micro-checker: software for identifying and correcting genotyping errors in microsatellite data. Mol. Ecol. Notes 4, 535–538 (2004).
Excoffier, L., Laval, G. & Schneider, S. ARLEQUIN Version 31, An Integrated Software Package for Population Genetics Data Analysis Computational and Molecular Population Genetics Lab CMPG, Institute of Zoology, University of Berne, Switzerland CMPG website, http//cmpgunibech/software/arlequin3 (2006).
SAS/STAT 9.22. User’s Guide The GLIMMIX Procedure. (SAS Institute Inc., Cary, NC 2010).
Wang, J., Xie, H. & Fisher, J.H. Multilevel models – Applications using SAS. De Gruyter Higher Education Press, Berlin (2012).
Acknowledgements
This work was supported by a Morris Animal Foundation grant to R.S. and K.B. ARC Discovery and Linkage grant (DP140103260, DP110102731, LP0989727) to K.B. and the ANR EVOCAN, André Hoffmann and la Fondation MAVA to F.T. We thank the Zoo and Aquarium Association and the Save the Tasmanian Devil Program and all participating zoos and wildlife parks for access to samples and stud book records.
Author information
Authors and Affiliations
Contributions
K.B. and R.S. designed the experiment. C.J.H. and C.H. collected the samples. T.R. conducted the molecular analyses. S.L., G.B. and M.O. conducted the statistical analyses. T.R., T.M. and B.U. wrote the manuscript. S.L, G.B., M.O., R.S., A.L., T.K., C.H., C.J.H., F.T. and K.B. contributed to the scientific discussion and revision of the manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Russell, T., Lisovski, S., Olsson, M. et al. MHC diversity and female age underpin reproductive success in an Australian icon; the Tasmanian Devil. Sci Rep 8, 4175 (2018). https://doi.org/10.1038/s41598-018-20934-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-018-20934-9
This article is cited by
-
Transcriptome annotation reveals minimal immunogenetic diversity among Wyoming toads, Anaxyrus baxteri
Conservation Genetics (2022)
-
Meta-analytic evidence that sexual selection improves population fitness
Nature Communications (2019)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
