Mate choice for major histocompatibility complex (MHC) complementarity in the Yellow-rumped Flycatcher (Ficedula zanthopygia)

Genes of the major histocompatibility complex (MHC) are an important component of the vertebrate immune system and play a significant role in mate choice in animal populations. However, the MHC genetic targets of female mate choice have not been clearly identified, and whether female mate choice is based on neutral genetic characteristics remains an open question. Here, we focus on the effects of morphological traits and genetic similarity among individuals in MHC class IIB (MHC IIB) exon 2 on mating in a sexually dimorphic songbird that exhibits social monogamy with extra-pair paternity (EPP). We sequenced 64 parent–offspring triads sampled over a 3-year period using two MHC class II loci to detect disassortative mating in the Yellow-rumped Flycatcher (Ficedula zanthopygia). We found that MHC similarity in social pairs was lower than that in random pairs. Extra-pair mate choice according to MHC IIB was observed, in which females’ extra-pair mates had fewer MHC alleles than their within-pair mates, but there was no significant band-sharing between extra-pair sires and potential extra-pair mates. However, the interaction between the MHC diversity of females and that of the social males affected the occurrence of EPP. Our results support the “optimality hypothesis” of MHC-based social and extra-pair choice. Female choice probably maintains a certain level of MHC diversity in offspring in the Yellow-rumped Flycatcher.


Background
The prevalence and significance of precopulatory mate choice remains a keenly debated topic in sexual selection because of its potentially important influence on fitness, which is generally considered an adaptive evolutionary mechanism (Forstmeier et al. 2014).
Mate choice is controversial mainly due to difficulties in quantifying the evolutionary costs and benefits of being "choosy" (Prokop et al. 2012;Forstmeier et al. 2014). Choosy females may obtain benefits derived from obtaining direct benefits (e.g., resources, territory or paternal care) (Griffin et al. 2013;Brouwer and Griffith 2019) and indirect benefits (e.g., genetic material) (Sheldon et al. 1997). Indirect benefits may not immediately contribute to the female's fitness but may benefit their progeny (genetic quality or compatibility). Paternity loss caused by extra-pair copulation

Open Access
Avian Research *Correspondence: wanght402@nenu.edu.cn 2 Jilin Engineering Laboratory for Avian Ecology and Conservation Genetics, School of Life Sciences, Northeast Normal University, Changchun 130024, China Full list of author information is available at the end of the article (EPC) is common among monogamous birds. EPCs are considered a way by which females could modify their initial mate choice (Brouwer and Griffith 2019). Several hypotheses have been proposed to explain the influence of male genetic quality on extra-pair mating, where quality is defined as the contribution of alleles or a genotype to individual fitness (Mays et al. 2008). According to genetic benefit hypotheses, females paired with high-quality males may be less prone to promiscuity, whereas those paired with poor-quality males may actively pursue EPCs to increase their offspring's fitness (Whittingham and Dunn 2010;Ferretti et al. 2011;Carlos and Botero 2012).
The major histocompatibility complex (MHC) contributes to the adaptive immune system and is a potential genetic marker of self-recognition versus non-self-recognition, pathogen-host interactions, and local adaptation (Eizaguirre et al. 2009;Spurgin et al. 2011). The MHC is regarded as an important candidate for mating preferences due to the potential genetic benefits that may be conferred to offspring (Penn and Potts 1998;Bos et al. 2009;Natsch et al. 2010). MHC-based mate-choice decisions have been reported across taxa, including mammals (Wang et al. 2010;Huchard et al. 2013;Ferrandiz-Rovira et al. 2016), birds (Juola and Dearborn 2012;Rymesova et al. 2017), reptiles (Olsson et al. 2003;Jaeger et al. 2016), amphibians (Bos et al. 2009) and fishes (Milinski et al. 2005;Agbali et al. 2010).
Theory suggests that females may gain indirect benefits, such as acquisition of "good" paternal genes (Hamilton and Zuk 1982) and enhanced survival of offspring, by choosing males with special or compatible genes (Brown 1997;Mays and Hill 2004). With respect to the MHC, females may choose among males to confer good paternal genes to their offspring, either by preferring "good" MHC haplotypes associated with enhanced resistance to pathogens or by avoiding "bad" MHC haplotypes associated with deleterious effects (Abplanalp et al. 1992;Wedekind et al. 1996). Alternatively, but not exclusively, females may seek the compatibility of maternal and paternal MHC genes to improve offspring MHC characteristics, either by choosing males with dissimilar MHC genes to maximize offspring MHC diversity (Brown 1997;E et al. 2017) or by choosing males with which they share an intermediate number of MHC alleles to optimize offspring MHC diversity at an intermediate level (Kloch et al. 2010;Hawley and Fleischer 2012). The optimality hypothesis is regarded as an extension of the genetic compatibility hypothesis (Han et al. 2019), which states that the optimal level of MHC diversity in offspring may be that conferring the most resistance against parasites and pathogens (Milinski et al. 2005;Cockburn et al. 2013).
In many socially monogamous species, individuals have been found to engage in extra-pair copulations (EPCs), which results in extra-pair paternity (EPP) (Westneat and Stewart 2003). EPP has been recorded in approximately 90% of investigated avian species (Griffith et al. 2002) and in other animal taxa. Social mate choice might be constrained by the need to acquire a dominant position before reproduction (Vaclav and Hoi 2007) and extrapair mate choice may be a way for females to adjust their mate choice and copulate with males possessing specific genetic characteristics (Gasparini et al. 2015). Extra-pair mating presents an opportunity to examine the processes governing mate choice in the absence of any potential direct benefits since the extra-pair male does not contribute to parental care (Dearborn et al. 2001). It is often predicted that females will choose extra-pair mates that are more ornamented than their social mates to gain superior genes for their offspring and, thus, indirect genetic benefits (Tschirren et al. 2012). Several studies have shown that MHC genes are important for mate choice, including extra-pair mate choice, in several passerine species, with females benefiting directly by mating with "healthy" males (e.g., being less likely to be infected with pathogens from their mates), but very few studies have addressed this question to date (Richardson et al. 2005;Santos et al. 2017;Rekdal et al. 2019). MHC genes thus provide an ideal study system to identify the specific functional loci effects of female mate choice (Penn and Potts 1998;Wedekind and Evanno 2010).
The Yellow-rumped Flycatcher (Ficedula zanthopygia) is a migratory, insectivorous, solitary-nesting species with biparental care (Alcaide et al. 2007;E et al. 2017), a high degree of sexual dimorphism in socially monogamous mating and moderate rates of EPP [with almost 50% of nests contain extra-pair young (EPY)] (E et al. 2017). According to microsatellite genotyping studies (E et al. 2017), female flycatchers are likely to choose males that are larger and relatively more heterozygous than their social mates as extra-pair mates. To assess MHC genotypes, we sequenced the second exon of the MHC class II B locus, which aids in the recognition of extracellular parasites and codes for the most polymorphic segment of the peptide-binding region in vertebrate genomes (Sallaberry et al. 2016). We predicted that the highly restrictive nature of social mate choice would prevent MHC-based social mate choice but would favour the evolution of MHC-based extra-pair mate choice. In this study, we compared the MHC genotypes of females with those of their social and extra-pair mates to determine whether within-pair or extra-pair matings were related to the overall MHC diversity of males or the compatibility of female and male genotypes. We specifically tested whether females mated preferentially with (1) more MHC-dissimilar mates, (2) more MHC-diverse males, or (3) males that provided an intermediate (optimal) level of MHC diversity for offspring. Finally, we investigated whether links were evident between the MHC II B genotype and phenotypic traits or body condition in the males preferred by females.

Study population
We studied a population of Yellow-rumped Flycatchers nesting at Zuojia Nature Reserve (126°1′-127°2′ N, 44°6′-45°5′ E) in Jilin Province, China. A detailed description of the study site and field procedures is given in E et al. (2017). The dataset included samples of 453 individuals (representing 64 nests and comprising 64 females, 64 males and 325 nestlings) collected during the breeding season in 2011-2013. The geographical position of each nest was recorded using a GPS device (N400PLUS, BHCnav, Beijing, China), and these data were used to calculate the linear distances between all breeding pairs in the study plot. To reduce the likelihood of nest abandonment as a result of early interference, a blood sample (approximately 20 μL) was collected from all adults and 6-7-dayold chicks and stored in absolute ethanol at −20 °C until DNA extraction.

MHC class II genotyping
Adult and nestling DNA was extracted from blood samples using a UNIQ-10 column animal genomic DNA isolation kit (SK1206, Sangon, Shanghai). To assess paternity, we genotyped all individuals at ten polymorphic microsatellite loci. The amplification conditions in this study are described in Leder et al. (2008); for more details on parentage analysis see E et al. (2017). We used the 326 and 325 degenerate primers (Ekblom et al. 2003) to perform partial amplification of exon 2. Using the resulting sequence, we designed a new primer pair (13-MHCNF: CTAA CGGC ACGG AGAA GG and 13-MHCNR: TCTG CGAT CCAC GCTG AA) that amplified a 240 bp sequence comprising exon 2, the region that contributes the most polymorphism to the peptide-binding site. For genotyping, we carried out denaturing gradient gel electrophoresis (DGGE) (Bonneaud et al. 2004). Each DGGE band was considered to correspond to one allele (Milinski et al. 2005). We grouped heterozygotes based on similarity in DGGE profile and then cloned one to three samples for each group using the PGEM-T Easy Vector System (Promega, Madison, WI, USA). All sequences were corrected, aligned and compared with those previously described for the same species (Wittzell et al.1999; Baratti et al. 2010).

Statistical analyses
We used CHROMAS version 2.01 (Technelesyum Pty Ltd, Brisbane, QLD, Australia) to correct the chromatograms. We translated the sequences, which had been previously aligned with CLUSTAL X version 1.83, and calculated the proportion of amino acid differences between mates for each allele (p-distance) in MEGA v. 5.05 (Tamura et al. 2013). Then, we averaged those values for all the pairwise comparisons of alleles between individuals. Our estimates of MHC diversity included the total number of alleles in an individual. Allele sharing was calculated to estimate MHC similarity between males and females forming pair bonds. Allele sharing was estimated as twice the number of shared alleles divided by the number of different alleles of each individual, i.e., D = 2F ab /(F a + F b ) (Wetton et al. 1987). Allele sharing does not take into account the functional similarity of individuals, so we also obtained a second estimate based on amino acids. Body condition was evaluated using the residuals from a linear regression of body mass against tarsus length. Generalized linear mixed models (GLMMs) were used to assess the effect of allelic diversity (number of identified MHC class I alleles) of social fathers and mothers on the occurrence of EPP. Monitoring year and nest box ID were included as random factors to avoid pseudo replication by incorporating repeated measures, EPP was included as the dependent variable, and MHC band-sharing between the two social mates, clutch size, and the date of onset of laying were included as independent variables. The GLMMs were fitted using the "lme4" package of R.
We examined the MHC diversity of these males and its similarity to that of females in four sets of comparisons. First, we used paired t tests to compare the social mate with one or more nearby "unpaired" males that a female could have paired with but that were not chosen and apparently were less preferred (n = 44) (the 44 nests included in the analysis represented those cases where a nearest neighbour existed). These "unpaired" males did not have a mate at the time when the focal female settled, but most eventually gained a social mate. In those cases where multiple unpaired males were present (up to three), the values were averaged. Previous paternity analyses indicated that 95% of EPY are not sired by immediate neighbours (E et al. 2017). Second, we used independent-sample t tests to compare MHC diversity and similarity between nests without EPY (n = 29, representing the males that did not lose paternity) and without EPY (n = 35, containing at least one EPY), assuming that males with some EPY were less preferred than males without. Third, we used paired t tests to examine MHC diversity and similarity for within-pair males and extrapair sires for the same nest. Last, we used paired t tests to compare the MHC diversity and similarity of known extra-pair mates with those of potential extra-pair mates (neighbouring males).
To test whether partners were more dissimilar at the MHC locus than expected on the basis of random pairs, differences in MHC amino acid distance, functional distance and allele sharing between true pairs and randomly chosen pairs were analysed. Bootstrap confidence tests were performed with 10,000 permutations. We tested for correlation between phenotypic traits and the number of MHC alleles of the males, the occurrence of EPP and MHC band sharing between social mates using Spearman's correlation test. All statistical analyses were performed in R 3.5.3 (R Core Development Team, http:// cran.r-proje ct. org/).

Social mate choice
Compared with the social fathers, randomly assigned males did not have a greater number of MHC class II alleles that differed from those of the female. We found no significant difference in the mean or variance of amino acid distance of social mates from that of random dyads (t = 1.809, df = 63, p = 0.075; random: mean = 0.168, social: mean = 0.165) (Fig. 2). However, we found that social mates had higher band-sharing than random dyads (t = − 3.135, df = 63, p = 0.003; random: mean = 0.568, social: mean = 0.638) (Fig. 3). In addition, within-pair mates and unpaired males did not differ in the total number of class II alleles. However, with respect to relative similarity between females and their mates, compared with unpaired males, within-pair males shared more MHC alleles with the female (paired t test: t 44 = 3.111, p = 0.003) and exhibited lower amino acid distances from the female (paired t test: t 44 = 2.076, p = 0.044) ( Table 1).

EPP
Compared with extra-pair mates, within-pair mates shared more MHC alleles with the females (paired t test: t 34 = 2.644, p = 0.012) and had significantly fewer amino acid differences from the females in their MHC alleles (paired t test: t 34 = 2.337, p = 0.026) ( Table 2). However, there was no significant difference in MHC diversity between within-pair males and extra-pair males (paired t test: t 34 = 1.108, p = 0.276). Neither the female's MHC diversity (GLMM: χ 2 = 0.2381, df = 1, p = 0.626) nor the male's MHC diversity (GLMM: χ 2 = 0.354, df = 1, p = 0.552) significantly influenced the Fig. 1 Percentage of EPY in a brood in relation to the MHC class II amino MHC band-sharing between the female and her within-pair mate Fig. 2 Distribution of the 10,000 bootstrapped MHC distance values for pairs assembled at random from the whole male population (i.e., mated and non-mated males) occurrence of EPP. The interaction between the female's MHC diversity and the social male's MHC diversity had a significant effect on the occurrence of EPP (GLMM: χ 2 = − 0.337, df = 1, p = 0.027) ( Table 3; Fig. 4).

Correlates of MHC diversity
Phenotypes that may be linked to male quality were then tested for a relationship with MHC diversity to determine whether they could be the direct cause of the link between extra-pair mate choice and MHC diversity. There was no correlation between MHC diversity and tarsal (Spearman's rho: r = 0.264, n = 64, p = 0.086), bill (Spearman's rho: r = 0.223, n = 64, p = 0.076), tail (Spearman's rho: r = − 0.036, n = 64, p = 0.777), wing (Spearman's rho: r = − 0.002, n = 64, p = 0.377) or body (Spearman's rho: r = − 0.010, n = 64, p = 0.851) lengths traits. The body condition index derived from the regression of weight against tarsal length was not related to the MHC diversity of the male (Spearman's rho: r = 0.037, n = 64, p = 0.771).

Discussion
In our study, we found evidence that female Yellowrumped Flycatchers tended to mate with males showing intermediate MHC diversity. It has been suggested that both natural and sexual selection may be driving the evolution of intermediate levels of diversity as a result of trade-offs between inbreeding and outbreeding depression (Bonneaud et al. 2006;Hansson and Westerberg 2008;Gagnon et al. 2020). In small populations, genetic diversity will be depleted if the majority of successful pairings have high pairwise similarities (leading to inbreeding), and will increase with dissimilar pairings Fig. 3 Distribution of the 10,000 bootstrapped MHC band-sharing values for pairs assembled at random from the whole male population (i.e., mated and non-mated males)

Table 1 Paired comparisons of MHC variation between within-pair males and nearest unpaired mates
Comparisons were made between within-pair mates and nearest unpaired mates (male diversity) and between the similarity estimates for the female and each type of male (similarity to the female). *Significant with p < 0.05 as per the advantage asserted by the dissimilar mates' hypothesis (although outbreeding depression may be a potential risk; Chargé et al. 2014). We also compared the occurrence of EPP and MHC band sharing between social mates. We found a negative relationship of the incidence of EPP and the relatedness between social mates, which indicates the potential occurrence of outbreeding or some level of local mate competition (LMC), in this mating system (Aoki et al. 2007).
Several studies have found that it is possible for a population to maintain high levels of MHC diversity even in the face of extremely severe bottlenecks (i.e., in genetically depauperate populations) and/or under high inbreeding (Richardson et al. 2005;Ellison et al. 2011).  (Wegner et al. 2003). After longdistance migration, Yellow-rumped Flycatchers may be in generally poorer condition than they were before migration and because of susceptibility to novel parasites in the new habitat, more infected and less fit (Stutz and Bolnick 2017). Another possibility is that females choose similar mates to avoid disruption of locally coadapted genes (Bonneaud et al. 2006;Eizaguirre et al. 2009), although this possibility awaits detailed study. However, we cannot exclude the possibility that male-male competition and cryptic female choice bias fertilization patterns (Andersson 1994; Jennions and Petrie 2000; Kamiya et al. 2014) and contributed to the observed pattern in our study. Some studies have revealed correlations among MHC variation, external traits and mate choice. Females may assess the MHC genotypes of males by using phenotypic traits as cues. For example, in the endangered Crested Ibis (Nipponia nippon), mates are selected using UV visual capability ). In the Common Pheasant (Phasianus colchicus), relationships between male phenotypic traits and certain MHC genotypes have been reported (Ekblom et al. 2003;Baratti et al. 2010;Hale et al. 2014;Yu et al. 2018). Thus, our study has some limitations. Although we found evidence that female Yellow-rumped Flycatchers prefer males with intermediate MHC II diversity, the phenotypic traits mediating the direct recognition of partners remain unknown and are difficult to observe. It is possible that such correlations could arise because both MHC and external traits are correlated with individual body condition, which is determined by other factors. In studies of mammal or fish taxa with well-developed olfaction, individuals are able to recognize both their own genotype and the genotype of potential partners, as observed in humans (Thornhill et al. 2003), mice (Firman and Simmons 2015), fish (Milinski et al. 2005;Gahr et al. 2018), and lizards (Olsson et al. 2003). However, most bird species are considered microsmatic or even anosmic (Roper 1999). Nevertheless, recent experimental work revealed MHC-dependent mate choice suggestive of an olfactory capability in Blue Petrels (Halobaena caerulea) (Strandh et al. 2012). Furthermore, Leclaire et al. (2017) showed that similarity in preen secretion chemicals is positively correlated with MHC relatedness in Black-legged Kittiwakes (Rissa tridactyla). The associations between the MHC, odours, and mating preferences in birds are a promising avenue for further exploration.
Associations between MHC genotype and social mate choice or extra-pair mate choice have been reported in a few species of vertebrates, and suggest a preference for the most heterozygous individuals or the use of MHC-based cues to evaluate the degree of relatedness of a potential partner to avoid inbreeding or outbreeding. However, the optimal level of genomic divergence for maximizing reproductive success be an intermediate level (Neff 2005;Slade et al. 2017). In this outbred population of Yellow-rumped Flycatchers, female choice probably maintains a certain level of diversity in the offspring while avoiding the costs of the disruption of local adaptations. We are aware that our sample size is modest in order to discovery of (dis-)assortative mating, although given the placement of our empirical mean in relation to the distribution generated by randomisations. Nevertheless, precise estimates (obtainable through the use of large sample sizes and/or very informative genetic markers) are required in further studies.

Conclusion
Our results support the optimality hypothesis and contribute to our knowledge of MHC-based social and extrapair choice in passerine birds. It appears that intermediate levels of MHC dissimilarity allow females to increase their inclusive fitness (Bichet et al. 2014) in the absence of significant inbreeding depression (Lehtonen and Kokko 2015;Sepil et al. 2015). However, our study gives rise to additional unsolved questions regarding MHC-mediated mate choice. For example, olfactory, colour or morphological traits are often found as reliable signals of the MHC profile, but our analysis shows that such traits, which might reflect male genetic quality were not correlated with male MHC diversity in the Yellow-rumped Flycatcher. Thus, MHC-mediated mate choice in this species may be determined by other factors, a possibility that warrants further study.