# MHC diversity and female age underpin reproductive success in an Australian icon; the Tasmanian Devil

Tracey Russell, Simeon Lisovski, Mats Olsson, Gregory Brown, Rebecca Spindler, Amanda Lane, Tamara Keeley, Chris Hibbard, Carolyn J Hogg, Frédéric Thomas, et al.

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Tracey Russell, Simeon Lisovski, Mats Olsson, Gregory Brown, Rebecca Spindler, et al. MHC diversity and female age underpin reproductive success in an Australian icon; the Tasmanian Devil. Scientific Reports, 2018, 8, pp.4175. 10.1038/s41598-018-20934-9. hal-02504096

HAL Id: hal-02504096

https://hal.umontpellier.fr/hal-02504096

Submitted on 10 Mar 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

OPEN

# MHC diversity and female age underpin reproductive success in an Australian icon; the Tasmanian Devil

Received: 7 June 2017

Accepted: 24 January 2018

Published online: 08 March 2018

Tracey Russell<fnref n="1" />, Simeon Lisovski<fnref n="1" />,², Mats Olsson³, Gregory Brown<fnref n="1" />, Rebecca Spindler⁴, Amanda Lane<fnref n="1" />, Tamara Keeley⁵, Chris Hibbard⁶, Carolyn J. Hogg<fnref n="1" />, Frédéric Thomas⁷, Katherine Belov<fnref n="1" />, Beata Ujvari<fnref n="1" />,⁸ & Thomas Madsen<fnref n="1" />,⁸,⁹

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.

The major histocompatibility complex (MHC) has been shown to be one of the key molecular determinants of mate choice in numerous vertebrates including humans<fnref n="1" />⁻⁷. 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 system⁸,⁹. 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 neutrality¹⁰.

Pathogen-driven selection has been documented to result in two, not mutually exclusive, MHC responses: (i) selection for specific MHC alleles¹¹⁻¹⁴ and/or (ii) selection for enhanced MHC polymorphism¹³⁻¹⁶. 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 combinations¹⁷,¹⁸. Consequently, diversifying MHC selection might be counteracted by the deleterious effects caused by autoimmunity¹⁹. Indeed, theoretical models as well as empirical studies have shown that optimal pathogen resistance often occurs at an optimal intermediate level of MHC polymorphism¹³,²⁰⁻²⁷.

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

<footnote n="1">¹School of Life and Environmental Sciences, University of Sydney, Sydney, 2006, NSW, Australia. ²University of California Davis, Department of Neurobiology Physiology and Behavior, Davis, CA, 95616, USA. ³Dept of Biological and Environmental Sciences, University of Gothenburg, Box 463, SE 405 30, Gothenburg, Sweden. ⁴Taronga Conservation Society Australia, Mosman, 2088, NSW, Australia. ⁵School of Agriculture and Food Science, University of Queensland, Gatton, 4343, QLD, Australia. ⁶Zoo and Aquarium Association Australasia, Mosman, 2088, NSW, Australia. ⁷CREEC/MIVEGEC, UMR IRD/CNRS/UM 5290, 911 Avenue Agropolis, BP 64501, 34394, Montpellier Cedex 5, France. ⁸Centre for Integrative Ecology, School of Life and Environmental Science, Deakin University, Geelong, 3216, VIC, Australia. ⁹School of Biological Sciences, University of Wollongong, Wollongong, 2522, NSW, Australia. Correspondence and requests for materials should be addressed to T.M. (email: madsen@uow.edu.au)</footnote>
lifetime reproductive success regardless of an individual's additional genotype²⁸, and/or (ii) acquiring optimal genetic compatibility²⁹, and/or (iii) achieving enhanced genetic diversity³⁰.

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 depression³¹. 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*)³². 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 species³³,³⁴. This devastating disease is spread among devils via biting during social interactions³⁵. 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 expression³⁶,³⁷. Metabolic failure, tumour related cachexia and metastases result in devil death within 6 to 9 months of the emergence of the first lesions³⁸. In order to ensure that the Tasmanian devil will not face extinction a large scale captive breeding program commenced in 2005³⁹, but overall breeding success is still not optimal⁴⁰,⁴¹. 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).

## 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 years³³. 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 MHC⁴²⁻⁴⁵. The low level of genetic polymorphism has been suggested to be caused by population bottlenecks during the Pleistocene and Holocene⁴⁶,⁴⁷ and by pathogens such as a canine-distemper-like disease in the early twentieth century⁴⁴. The genetic diversity of the captive devils used in the present study was similar to that recorded in devils captured in the wild⁴⁸. 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*)⁴⁹, savannah sparrows (*Passerculus sandwichensis*)⁵⁰, southern dunlins (*Calidris alpina schinzii*)⁵¹, great tits (*Parus major*)⁵² and the cynomolgus macaque (*Macaca fascicularis*)⁵³ 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 function¹³,²⁰⁻²⁷. 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 habitats⁵⁴, 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 humans¹⁻⁷,⁵⁵. MHC-based post-copulatory female cryptic choice of male sperm has also been documented in several vertebrates such as fish⁵⁶, lizards⁴⁹, birds⁵⁷ and mammals⁵⁸. As female devils in both the successful and unsuccessful pair combinations were observed mating⁵⁹⁻⁶⁵, 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 genotype⁶⁶,⁶⁷. However, if a similar mechanism may influence male fertilization success in internal fertilizers, such as mammals, is to our knowledge unknown.
<table><thead><tr><th>number of female heterozygous loci</th><th>female age (months)</th><th>female ID</th><th>number of male heterozygous loci</th><th>male age (months)</th><th>male ID</th><th>pair combination</th></tr></thead><tbody><tr><td>4</td><td>48</td><td>FD1</td><td>1</td><td>48</td><td>MD1</td><td>successful</td></tr><tr><td>4</td><td>36</td><td>FD1</td><td>2</td><td>36</td><td>MD6</td><td>unsuccessful</td></tr><tr><td>4</td><td>36</td><td>FD1</td><td>3</td><td>36</td><td>MD21</td><td>unsuccessful</td></tr><tr><td>3</td><td>24</td><td>FD10</td><td>2</td><td>36</td><td>MD40</td><td>unsuccessful</td></tr><tr><td>4</td><td>36</td><td>FD11</td><td>1</td><td>60</td><td>MD8</td><td>successful</td></tr><tr><td>4</td><td>24</td><td>FD11</td><td>3</td><td>36</td><td>MD22</td><td>successful</td></tr><tr><td>6</td><td>48</td><td>FD12</td><td>1</td><td>60</td><td>MD8</td><td>successful</td></tr><tr><td>6</td><td>24</td><td>FD12</td><td>3</td><td>36</td><td>MD10</td><td>successful</td></tr><tr><td>3</td><td>24</td><td>FD13</td><td>4</td><td>24</td><td>MD15</td><td>unsuccessful</td></tr><tr><td>3</td><td>24</td><td>FD13</td><td>3</td><td>36</td><td>MD27</td><td>unsuccessful</td></tr><tr><td>4</td><td>36</td><td>FD14</td><td>2</td><td>24</td><td>MD30</td><td>unsuccessful</td></tr><tr><td>4</td><td>36</td><td>FD14</td><td>2</td><td>36</td><td>MD31</td><td>unsuccessful</td></tr><tr><td>4</td><td>36</td><td>FD14</td><td>3</td><td>24</td><td>MD37</td><td>unsuccessful</td></tr><tr><td>4</td><td>36</td><td>FD14</td><td>3</td><td>48</td><td>MD41</td><td>unsuccessful</td></tr><tr><td>4</td><td>24</td><td>FD15</td><td>1</td><td>24</td><td>MD32</td><td>successful</td></tr><tr><td>6</td><td>24</td><td>FD16</td><td>2</td><td>24</td><td>MD11</td><td>successful</td></tr><tr><td>6</td><td>48</td><td>FD16</td><td>3</td><td>48</td><td>MD43</td><td>successful</td></tr><tr><td>5</td><td>24</td><td>FD17</td><td>1</td><td>36</td><td>MD18</td><td>successful</td></tr><tr><td>5</td><td>36</td><td>FD17</td><td>2</td><td>36</td><td>MD31</td><td>successful</td></tr><tr><td>3</td><td>36</td><td>FD18</td><td>1</td><td>24</td><td>MD8</td><td>unsuccessful</td></tr><tr><td>3</td><td>24</td><td>FD18</td><td>4</td><td>48</td><td>MD15</td><td>successful</td></tr><tr><td>1</td><td>24</td><td>FD19</td><td>4</td><td>36</td><td>MD5</td><td>successful</td></tr><tr><td>1</td><td>36</td><td>FD19</td><td>3</td><td>36</td><td>MD33</td><td>successful</td></tr><tr><td>5</td><td>24</td><td>FD2</td><td>2</td><td>36</td><td>MD29</td><td>successful</td></tr><tr><td>5</td><td>24</td><td>FD20</td><td>1</td><td>36</td><td>MD8</td><td>successful</td></tr><tr><td>5</td><td>48</td><td>FD20</td><td>1</td><td>36</td><td>MD18</td><td>unsuccessful</td></tr><tr><td>5</td><td>48</td><td>FD20</td><td>2</td><td>36</td><td>MD19</td><td>unsuccessful</td></tr><tr><td>2</td><td>24</td><td>FD21</td><td>3</td><td>24</td><td>MD34</td><td>successful</td></tr><tr><td>4</td><td>48</td><td>FD22</td><td>3</td><td>48</td><td>MD43</td><td>successful</td></tr><tr><td>5</td><td>24</td><td>FD23</td><td>0</td><td>24</td><td>MD36</td><td>successful</td></tr><tr><td>3</td><td>24</td><td>FD24</td><td>4</td><td>24</td><td>MD12</td><td>successful</td></tr><tr><td>1</td><td>36</td><td>FD25</td><td>1</td><td>60</td><td>MD1</td><td>unsuccessful</td></tr><tr><td>1</td><td>24</td><td>FD25</td><td>4</td><td>24</td><td>MD38</td><td>successful</td></tr><tr><td>1</td><td>36</td><td>FD25</td><td>3</td><td>36</td><td>MD39</td><td>unsuccessful</td></tr><tr><td>2</td><td>24</td><td>FD26</td><td>4</td><td>24</td><td>MD12</td><td>successful</td></tr><tr><td>5</td><td>36</td><td>FD27</td><td>4</td><td>48</td><td>MD3</td><td>unsuccessful</td></tr><tr><td>5</td><td>24</td><td>FD27</td><td>1</td><td>24</td><td>MD32</td><td>unsuccessful</td></tr><tr><td>1</td><td>24</td><td>FD28</td><td>4</td><td>36</td><td>MD17</td><td>successful</td></tr><tr><td>2</td><td>36</td><td>FD29</td><td>4</td><td>48</td><td>MD26</td><td>unsuccessful</td></tr><tr><td>2</td><td>48</td><td>FD29</td><td>2</td><td>48</td><td>MD31</td><td>unsuccessful</td></tr><tr><td>3</td><td>24</td><td>FD3</td><td>4</td><td>36</td><td>MD5</td><td>successful</td></tr><tr><td>5</td><td>36</td><td>FD30</td><td>2</td><td>24</td><td>MD40</td><td>successful</td></tr><tr><td>5</td><td>60</td><td>FD31</td><td>2</td><td>60</td><td>MD23</td><td>unsuccessful</td></tr><tr><td>3</td><td>36</td><td>FD32</td><td>6</td><td>48</td><td>MD16</td><td>unsuccessful</td></tr><tr><td>3</td><td>36</td><td>FD32</td><td>4</td><td>48</td><td>MD26</td><td>unsuccessful</td></tr><tr><td>3</td><td>36</td><td>FD33</td><td>2</td><td>36</td><td>MD4</td><td>successful</td></tr><tr><td>3</td><td>48</td><td>FD33</td><td>4</td><td>48</td><td>MD25</td><td>successful</td></tr><tr><td>3</td><td>24</td><td>FD34</td><td>4</td><td>36</td><td>MD3</td><td>successful</td></tr><tr><td>3</td><td>36</td><td>FD34</td><td>1</td><td>36</td><td>MD32</td><td>unsuccessful</td></tr><tr><td>1</td><td>24</td><td>FD35</td><td>4</td><td>24</td><td>MD12</td><td>successful</td></tr><tr><td>2</td><td>24</td><td>FD36</td><td>4</td><td>36</td><td>MD7</td><td>successful</td></tr><tr><td>3</td><td>24</td><td>FD37</td><td>1</td><td>24</td><td>MD8</td><td>unsuccessful</td></tr><tr><td>1</td><td>36</td><td>FD4</td><td>4</td><td>36</td><td>MD42</td><td>successful</td></tr><tr><td>4</td><td>48</td><td>FD5</td><td>3</td><td>36</td><td>MD13</td><td>successful</td></tr><tr><td>4</td><td>36</td><td>FD5</td><td>3</td><td>24</td><td>MD33</td><td>successful</td></tr><tr><td>3</td><td>48</td><td>FD6</td><td>4</td><td>48</td><td>MD2</td><td>unsuccessful</td></tr></tbody></table>

Continued
<table><caption>**Table 1.** Female – male pair combinations.</caption><thead><tr><th>number of female heterozygous loci</th><th>female age (months)</th><th>female ID</th><th>number of male heterozygous loci</th><th>male age (months)</th><th>male ID</th><th>pair combination</th></tr></thead><tbody><tr><td>3</td><td>48</td><td>FD6</td><td>5</td><td>48</td><td>MD14</td><td>unsuccessful</td></tr><tr><td>3</td><td>48</td><td>FD6</td><td>5</td><td>48</td><td>MD35</td><td>unsuccessful</td></tr><tr><td>2</td><td>36</td><td>FD7</td><td>2</td><td>48</td><td>MD9</td><td>unsuccessful</td></tr><tr><td>2</td><td>24</td><td>FD7</td><td>2</td><td>36</td><td>MD20</td><td>successful</td></tr><tr><td>2</td><td>36</td><td>FD7</td><td>2</td><td>36</td><td>MD24</td><td>unsuccessful</td></tr><tr><td>5</td><td>24</td><td>FD8</td><td>4</td><td>36</td><td>MD38</td><td>successful</td></tr><tr><td>5</td><td>36</td><td>FD8</td><td>3</td><td>48</td><td>MD39</td><td>successful</td></tr><tr><td>3</td><td>60</td><td>FD9</td><td>2</td><td>60</td><td>MD23</td><td>unsuccessful</td></tr><tr><td>3</td><td>36</td><td>FD9</td><td>2</td><td>36</td><td>MD28</td><td>unsuccessful</td></tr></tbody></table>

<table><caption>**Table 2.** Final results of logistic mixed model.</caption><thead><tr><th>Effect</th><th>Estimate</th><th>SE</th><th>DF</th><th>t</th><th>Pr &gt; |t|</th></tr></thead><tbody><tr><td colspan="6"><strong>Parameter estimate</strong></td></tr><tr><td>Intercept</td><td>1.9734</td><td>1.1183</td><td>62</td><td>1.76</td><td>0.0826</td></tr><tr><td>Difference in number of heterozygous loci</td><td>0.6204</td><td>0.2539</td><td>62</td><td>2.44</td><td>0.0174</td></tr><tr><td>Female age</td><td>-0.0871</td><td>0.0316</td><td>62</td><td>-2.75</td><td>0.0077</td></tr><tr><td colspan="6"><strong>Type III Tests of Fixed Effects</strong></td></tr><tr><td>Effect</td><td>Num DF</td><td>Den DF</td><td>F</td><td>Pr &gt; F</td><td></td></tr><tr><td>Difference in number of heterozygous loci</td><td>1</td><td>62</td><td>5.97</td><td>0.017</td><td></td></tr><tr><td>Female age</td><td>1</td><td>62</td><td>7.58</td><td>0.008</td><td></td></tr></tbody></table>

<figure><img src="image_3.png" /><figcaption>**Figure 1.** 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).</figcaption></figure>

## 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.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** (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.</figcaption></figure>

**Material and Methods**

**Study animals.** Data on pairings were obtained from the Tasmanian devil studbook⁶⁸ and the insurance population annual reports⁵⁹⁻⁶⁵. 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 appetite⁶⁹. 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 & Belov⁷⁰; Sh-MHCI101, Sh-MHCI102, Sh-MHCI105, Sh-MHCI106, Sh-MHCI107, and Sh-MHCI108. For further details on size range and primer sequences see Cheng & Belov⁷⁰. Polymerase chain reactions (PCR) were carried out using the protocols of Cheng & Belov⁷⁰ 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-Checker⁷¹. The number of alleles per locus together with observed and expected heterozygosity were estimated using the software program ARLEQUIN 3.1⁷². 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.1⁷² (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 GLIMMIX⁷³,⁷⁴ 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**

1. 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).

2. 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).

3. 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).

4. Olsson, M. et al. Major histocompatibility complex and mate choice in sand lizards. *Proc. R. Soc. B* **270**, S254–S256 (2003).

5. Sin, Y. W. et al. MHC class II-assortative mate choice in European badgers *Meles meles*. *Mol. Ecol.* **24**, 3138–3150 (2015).

6. 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).

7. Yamazaki, K. & Beauchamp, G. K. (2007). Genetic basis for MHC dependent mate choice. *Adv. Genet.* **59**, 129–145 (2007).

8. Klein, J. Natural history of the Major Histocompatibility Complex. (Wiley, New York, 1986).

9. Zinkernagel, R. M. Associations between major histocompatibility antigens and susceptibility to disease. *Ann. Rev. Microbiol.* **33**, 201–213 (1979).

10. Edwards, S. V. & Hedrick, P. W. Evolution and ecology of MHC molecules, from genomics to sexual selection. *Trends Ecol. Evol.* **13**, 305–311 (1998).

11. 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).

12. Lohm, J. et al. Experimental evidence for majorhistocompatibility complex-allele-specific resistance to a bacterial infection. *Proc. R. Soc. B* **269**, 2029–2033 (2002).

13. Madsen, T. & Ujvari, B. MHC Class I variation associates with parasite resistance and longevity in tropical pythons. *J. Evol. Biol.* **19**, 1973–1978 (2006).

14. Westerdahl, H. et al. Association between malaria and MHC genes in a migratory songbird. *Proc. R. Soc. B* **272**, 1511–1518 (2005).

15. 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).

16. 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).

17. 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).

18. Vidovic, D. & Matzinger, P. Unresponsiveness to a foreign antigen can be caused by self-tolerance. *Nature* **336**, 222–225 (1988).

19. Fernando, M. M. A. et al. Defining the role of the MHC in autoimmunity, a review and pooled analysis. *PLoS Genet.* **4**, e1000024 (2008).

20. 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).

21. De Boer, R. J. & Perelson, A. S. How diverse should the immune system be? *Proc. R. Soc. B* **252**, 171–175 (1993).

22. Kalbe, M. et al. Lifetime reproductive success is maximized with optimal MHC diversity. *Proc. R. Soc. B* **276**, 925–934 (2009).

23. Milinski, M. The major histocompatibility complex, sexual selection and mate choice. *Ann. Rev. Ecol. Syst.* **37**, 159–186 (2006).

24. 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).

25. 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).

26. Wegner, K. M., Kalbe, M., Kurtz, J., Reusch, T. R. B. & Milinski, M. Parasite selection for immunogenic optimality. *Science* **301**, 1343 (2003).

27. 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).
28. Olsson, M. et al. MHC, health, color, and reproductive success in sand lizards. *Behav. Ecol. Sociobiol.* **58**, 289–294 (2005).

29. Kempenaers, B. Mate choice and genetic quality: a review of the heterozygosity theory. *Adv. Study Behav.* **37**, 189–278 (2007).

30. 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).

31. Madsen, T., Shine, R., Olsson, M. & Witzell, H. Restoration of an inbred adder population. *Nature* **402**, 34–35 (1999).

32. Belov, K. The role of the major histocompatibility complex in the spread of contagious cancers. *Mammal. Genome* **22**, 83–90 (2011).

33. Jones, M. et al. Life-history change in disease-ravaged Tasmanian devil populations. *Proc. Natl. Acad. Sci. USA* **105**, 10023–10027 (2008).

34. 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).

35. 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).

36. 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).

37. Karu, N. et al. Discovery of biomarkers for Tasmanian devil cancer (DFTD) by metabolic profiling of serum. *J. Proteome Res.* **15**, 3827–3840 (2016).

38. Pearse, A.-M. & Swift, K. Allograft theory transmission of devil facial-tumour disease. *Nature* **439**, 549 (2006).

39. 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).

40. Hockley, J. DPIPWE/ZAA Husbandry guidelines for Tasmanian devil, *Sarcophilus harrisii*. Zoo and Aquarium Association, Australia pp 46 (2011).

41. 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).

42. 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).

43. Epstein, B. et al. Rapid evolutionary response to a transmissible cancer in Tasmanian devils. *Nature Commun.* **7**, 12684 (2016).

44. 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).

45. 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).

46. Brown, O. Tasmanian devil *Sarcophilus harrisii* extinction on the Australian mainland in the Mid-Holocene, multicausality and ENSO intensification. *Alcheringa* **30**, 49–57 (2006).

47. 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).

48. Lane, A. et al. New Insights into the role of MHC diversity in devil facial tumour disease. *PLoS One* **7**, e36955 (2012).

49. Olsson, M., Shine, R., Madsen, T., Gullberg, A. & Tegelström, H. Sperm selection by females. *Nature* **383**, 585 (1996).

50. 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).

51. 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).

52. 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).

53. 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).

54. 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).

55. 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).

56. 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).

57. 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).

58. 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).

59. Lees, C. & Srb, C. Annual report for the DPIW-ARAZPA Tasmanian devil insurance population. Zoo and Aquarium Association. pp 27 (2007).

60. Lees, C. & Srb, C. Annual report for the DPIW-ARAZPA Tasmanian devil insurance population. Zoo and Aquarium Association. pp 49 (2008).

61. Hibbard, C. J., Ford, C. & Srb, C. Annual report for the DPIW-ZAA Tasmanian Devil insurance population. Zoo and Aquarium Association. pp 57 (2009).

62. 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).

63. 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).

64. Hogg, C. J., Srb, C. & Hockley, J. Annual report for the DPIPWE-ZAA Tasmanian devil insurance population. Zoo and Aquarium Association. pp 43 (2012).

65. Hogg, C. J., Srb, C. & Hockley, J. Annual Report for the DPIPWE-ZAA Tasmanian devil insurance population. Zoo and Aquarium Association. pp 42 (2013).

66. 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).

67. 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).

68. Srb, C. Australasian regional Tasmanian devil studbook. (Healesville, Victoria, 2014) (2014).

69. Brice, S. Zoo aquarium association husbandry guidelines for Tasmanian devil, *Sarcophilus harrisii*. Zoo and Aquarium Association, Australia pp 36 (2006).

70. 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).

71. 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).
## 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 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.

## Additional Information

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_5.png" /></figure>

**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/.

© The Author(s) 2018

<footnote n="72">72. 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).</footnote>

<footnote n="73">73. SAS/STAT 9.22. User's Guide The GLIMMIX Procedure. (SAS Institute Inc., Cary, NC 2010).</footnote>

<footnote n="74">74. Wang, J., Xie, H. & Fisher, J.H. Multilevel models – Applications using SAS. De Gruyter Higher Education Press, Berlin (2012).</footnote>