Skip to main content

Evidence that pairing with genetically similar mates is maladaptive in a monogamous bird



Evidence of multiple genetic criteria of mate choice is accumulating in numerous taxa. In many species, females have been shown to pair with genetically dissimilar mates or with extra-pair partners that are more genetically compatible than their social mates, thereby increasing their offsprings' heterozygosity which often correlates with offspring fitness. While most studies have focused on genetically promiscuous species, few studies have addressed genetically monogamous species, in which mate choice tends to be mutual.


Here, we used microsatellite markers to assess individual global heterozygosity and genetic similarity of pairs in a socially and genetically monogamous seabird, the black-legged kittiwake Rissa tridactyla. We found that pairs were more genetically dissimilar than expected by chance. We also identified fitness costs of breeding with genetically similar partners: (i) genetic similarity of pairs was negatively correlated with the number of chicks hatched, and (ii) offspring heterozygosity was positively correlated with growth rate and survival.


These findings provide evidence that breeders in a genetically monogamous species may avoid the fitness costs of reproducing with a genetically similar mate. In such species that lack the opportunity to obtain extra-pair fertilizations, mate choice may therefore be under high selective pressure.


Numerous traits influence mate choice that may produce non-random mating patterns in many species. Although most studies have focused on morphological and behavioral traits [14]), there is rapidly growing evidence of multiple genetic criteria of mate choice (reviewed in [5, 6]). For example, females might choose the most heterozygous males [7], which may increase the resistance of offspring to parasites [8, 9]. Alternatively, females might choose males carrying alleles that are compatible with their own genotypes. The main driving forces of mate choice would then be to maintain equilibrium between co-adapted genes [10] or alternatively, to enhance the genetic variability of offspring [1116].

In species with biparental care, both sexes should be choosy in selecting a mate [17, 18]. Blomqvist et al. [14] reported that in socially monogamous shorebirds, both sexes obtained extra-pair fertilizations when mates were genetically similar. In blue tits (Parus caeruleus), females acquired extra-pair fertilizations that enhanced the heterozygosity and fitness of their offspring [19]. In superb starlings (Lamprotornis superbus), the benefits of extra-pair fertilizations may differ according to the genetic similarity of the extra-pair partner [20]. However, few studies have focused on the mating patterns in genetically monogamous species which lack extra-pair fertilizations to diminish the costs of pairing with a suboptimal mate.

The black-legged kittiwake (Rissa tridactyla) is a long-lived, monogamous seabird with no extra-pair fertilizations [21] and high between-year repairing rates [2225]. In this species, mate choice may thus profoundly affect reproductive success throughout an individual's lifetime. To examine whether mating patterns in kittiwakes are influenced by genetic criteria, we assessed individual heterozygosity and genetic similarity of mates with microsatellite markers.

Our first aim was to examine three hypotheses of mating patterns driven by individual genotypes. Breeders may be paired with: (1) heterozygous mates ("preference of heterozygous mates" hypothesis; [15]), (2) genetically dissimilar mates, in order to increase the genetic variability of offspring ("genetic similarity avoidance" hypothesis; [14, 26, 27]), or (3) genetically similar mates, in order to preserve the link between locally co-adapted genes ("genetic similarity preference" hypothesis; [10, 28]). According to Hypothesis 1, the most heterozygous mates are of better quality because they may provide direct benefits (e.g., better parental care) and/or indirect benefits in term of more heterozygous offspring. The hypothesis predicts that paired individuals are more heterozygous than unpaired ones, and that there a positive correlation between male and female heterozygosity [7]. Hypothesis 2 predicts that the observed mean genetic similarity between pair members is lower than expected through random matings. Hypothesis 3 is the reverse of Hypothesis 2 and predicts that mates share more alleles than expected by chance, meaning that the observed mean genetic similarity between pair members will be higher than expected through random matings.

Our second aim was to examine the fitness consequences of breeding with genetically similar or dissimilar individuals. We searched for possible relationships between genetic characteristics of pairs and fitness components such as clutch size and hatching success [29, 30]. Because genetically similar pairs are more likely to produce homozygous offspring than dissimilar pairs, we also examined predicted relationships between offspring heterozygosity and their growth and survival.

Microsatellites are generally assumed to be neutral genetic markers and have been widely used to estimate relatedness, individual heterozygosity and population level of inbreeding [3134]. If heterozygosity at certain selected loci enhances fitness [8, 15], heterozygosity at microsatellite loci may be a good surrogate of the overall genetic quality of an individual, especially in wild species where little is known about genes under selection. However, Lieutenant-Gosselin & Bernatchez [35] and Tiira et al. [36] have shown that global heterozygosity might be highly influenced by heterozygosity at certain specific loci. Such loci may be physically linked to fitness loci, and one should therefore distinguish effects of global heterozygosity from effects of heterozygosity at each microsatellite marker. Thus we performed our analyses both with the global heterozygosity and with the heterozygosity at each locus in order to distinguish the effects of global heterozygosity from that of specific loci that might be physically linked to fitness loci.


We monitored 348 genotyped adults in 2003 and 2004: 241 were seen alive in 2003 and 289 in 2004. Adults formed 74 pairs in 2003 and 72 in 2004; the remaining adults for each year corresponded to unpaired adults or adults paired with non-genotyped mates. All these adults were included in the bootstrap analyses because they were alive in the considered year and thus potentially available for pairing.

For chicks with genotyped parents, we found that Phm xy of the parents was closely related to all indices of offspring heterozygosity (82 chicks, p < 0.007), confirming that Phm xy might be a reliable estimate of the probability of a given pair of producing homozygous offspring. Furthermore, Phm xy of pairs used in our study exhibited a wide, six-fold range of variation (from 0.06 to 0.39, mean: 0.18 +/- 0.06).

Mating pattern and genetics

To correct for linkages between loci, OHW loci and K31 were excluded from calculations of H, SH and IR, but results did not differ with all loci. The "preference of heterozygous mates" hypothesis (hypothesis 1) predicted assortative mating by heterozygosity. However, male and female H were not correlated for any year (r ranging from 0.021 to 0.10, p ranging from 0.38 to 0.86). Results were the same for SH (0.02 ≤ r ≤ 0.13, 0.26 ≤ p ≤ 0.87) and IR (0.05 ≤ r ≤ 0.10, 0.38 ≤ p ≤ 0.67). Heterozygosity of paired and unpaired individuals did not differ (for all indices: in 2003, t2,241 < 1.93, p > 0.06; in 2004, t2,241 < 1.25, p > 0.21).

The "genetic similarity avoidance" hypothesis (hypothesis 2) predicted that observed pairs comprise less genetically similar individuals than expected by chance. To test this, for each year we ran 10,000 bootstraps using the observed individuals. For each run, we used either a calculation across all loci, or excluding OHW loci (Table 1). We found a significant difference in terms of genetic similarity between observed and simulated pairs in 2003 and 2004 (Figure 1), indicating that pairs were formed of more genetically dissimilar individuals than expected by chance.

Figure 1
figure 1

Distribution of genetic similarity (calculated with the Phm index) in pairs observed in 2003 and/or 2004. Black bars represent the percentage of observed pairs (n = 92) and white bars represent a random distribution of the Phm of pairs obtained by simulating 10,000 pairings. Phm was calculated with all ten loci. Kittiwake pairs were less genetically similar than expected by chance (p = 0.016).

Table 1 Differences between observed and simulated means of genetic similarity of mates

When pooling pairs seen in 2003 and 2004 (one observation for each pair, column "all pairs" in Table 1), we still found a significant difference in genetic dissimilarity between observed and simulated pairs. Results were similar when taking only pairs of kittiwakes breeding together both years (considered as "stable pairs", Table 1). Because of a low divorce rate, we lacked a sufficient sample to compare the genetic similarity of divorced versus reunited pairs.

The "genetic similarity preference" hypothesis (hypothesis 3) predicts that mates share more alleles than expected by chance. This means that observed pairs should comprise more genetically similar individuals than expected by chance. We would then expect to have a p-value for Phm xy higher than 0.95, which we did not find (Table 1).

Reproductive success and genetic similarity

We used only 2003 data to determine the number of eggs laid and the number and proportion of eggs hatched, because the reproductive success of some of the pairs in 2004 and subsequent years may have been affected by other experiments after pair formation. Genetic similarity indices were computed without OHW loci and K31. The number of eggs laid was not correlated with pair genetic similarity (χ272 = 0.01, p = 0.91, n = 74). However, for pairs that laid eggs, the number of hatched chicks was correlated with Phm xy 269 = 4.0, p = 0.045, n = 71), with the mean number of chicks hatched being lower in more genetically similar pairs (Figure 2). Similarly, hatching rate was also negatively correlated with Phm xy 269 = 4.0, p = 0.045, n = 71).

Figure 2
figure 2

Relationship between number of hatched chicks and Phm in 2003. Here, the probability of producing a homozygous offspring (Phm) was calculated excluding out of Hardy-Weinberg (OHW) loci, but results were similarly significant with all ten loci. The number of hatched chicks decreased significantly with genetic similarity of pairs (see the text for statistical tests).

Offspring growth and survival in relation to offspring heterozygosity

Chicks were not genotyped for K67 (because of low variability). In 2005, we assessed the survival and growth in body weight, tarsus and wing length of 82 chicks until age 25 days. These life history parameters were highly correlated to hatching rank (χ2 > 12, p < 0.001 in all analyses). The hatching rank*chick heterozygosity interactions were non-significant for all heterozygosity indices (p > 0.6) and were thus discarded from the models. Chick survival was positively correlated with chick heterozygosity for all indices (H: χ2 = 5.2, p = 0.022, see Figure 3; SH: χ2 = 4.7, p < 0.03; IR: χ2 = 4.5, p = 0.03), but became non-significant when removing OHW loci and K31 from the computation of chick heterozygosity (H: χ2 = 3.6, p = 0.06; SH: χ2 = 2.5, p = 0.11; IR: χ2 = 3.3, p = 0.07).

Figure 3
figure 3

Relationship between heterozygosity and chick survival in 2005. We show the correlation between chick survival up to 25 days of age and chick heterozygosity estimated by the H index calculated over all loci. Chicks that survived 25 days were significantly more heterozygous than chicks that died (see the text for statistical tests).

When analyzing chick growth in body weight, wing and tarsus length, we found a significant interaction of Age*Chick heterozygosity*Hatching rank (p < 0.017, when chick heterozygosity was calculated over all loci, p < 0.089 when chick heterozygosity was calculated without OHW loci and K31; Table 2 and Additional file 1). This effect of chick heterozygosity on chick growth was however only evident for B-chicks (Table 2, see Figure 4 for a representation of this effect on chick growth in wing length).

Figure 4
figure 4

Relationship between chick heterozygosity, age and growth rate in 2005. We plotted the predicted value of chick growth in wing length (in mm) according to the model A+H+R+AH+AR+HR+ARH+A2 where A is the age of the chick (measured in days), R its hatching rank (binary factor equaling A for first hatched, or B for second hatched chicks) and H its heterozygosity (according to the H index). For clarity, the random parameter (chick's identity) was removed from the model. Significance of the statistical tests is given in Table 2.

Table 2 Relationships between chick age, hatching rank, heterozygosity and growth in body weight, wing and tarsus length

Local effects and global effects of microsatellite loci

For mating patterns in 2003, we found that pairs were formed of genetically dissimilar individuals according to three of the 10 loci (Phm index for K32, K6 and K71). Hatching success was negatively correlated to genetic similarity of pairs in 2003 for two loci (Phm index for K32 and RBG20). Chick heterozygosity and survival until 25 days (in 2005) were positively correlated only for RBG29 (H index, p = 0.027; p > 0.1 for all other loci). For chick growth, we found that the parameters associated with the interactions Age*Heterozygosity and Age*Heterozygosity*Rank were significant in explaining chick growth in body weight and wing length when chick heterozygosity was estimated through three of the nine loci (H index for K16, K31, K32). Similarly, the parameters associated with the interactions Age*Heterozygosity and Age*Heterozygosity*Rank were significant in explaining chick growth in tarsus length when chick heterozygosity was calculated with two loci (H index for K31 and K32).


Our main objective was to examine predicted relationships between genetic variables and mating patterns in a genetically monogamous species. Selective pressures on mate choice may be stronger in such species, which do not mitigate the costs of homozygous offspring with extra-pair fertilizations. To our knowledge, there is only one other study of this kind of genetically monogamous species, the New Zealand robin (Petroica australis) and saddlebacks (Philesturnus carunculatus), in which pairings were random in relation to genetic similarity of mates [37]. Our findings in kittiwakes therefore appear to be the first evidence in a strictly monogamous species for the genetic similarity avoidance hypothesis.

Indeed, we found that kittiwake breeders were not paired randomly, but with mates that were less genetically similar than expected by chance. Observed pairs had a lower probability of producing homozygous offspring than expected, a finding that fits the predictions of the "genetic similarity avoidance" hypothesis. In contrast, the "genetic similarity preference" hypothesis predicted opposite results, and the "preference of heterozygous mates" hypothesis predicted assortative mating according to heterozygosity, a pattern we did not find.

Our findings show that genetically monogamous species may avoid pairing with genetically similar mates. Non-random mating with respect to genetic similarity has been reported in three shorebird species [14], ruffs (Philomachus pugnax, [11]), sand lizards (Lacerta agilis, [12]), mice [38] and humans [39]. In contrast, some other species behave consistently with the "good-genes as heterozygosity" hypothesis (e.g. house sparrows Passer domesticus, [7], and Seychelles warblers Acrocephalus sechellensis, [40]), or with the "genetic similarity preference" hypothesis (e.g. pied flycatcher Ficedula hypoleuca, [41], and great frigate birds Fregata minor, [10]), while other species display no significant relationships between mating patterns and genetic dissimilarity indices (e.g. great reed warbler Acrocephalus arundinaceus, [42]).

As expected, we also found deleterious effects of genetic similarity and homozygosity of offspring. This set of results (offspring hatchability, growth and survival) is statistically independent of results about mating patterns, but are biologically linked and consistent. Genetically similar pairs hatched fewer offspring than more dissimilar pairs (Figure 2), which is similar to studies that found significant effects of inbreeding on egg hatchability [29, 30]. Likewise in these studies, hatching success was found to correlate negatively with pair genetic similarity. This pattern could be due to at least two non-exclusive effects. First, because both sexes incubate, cooperation between mates may be maximized when the genetic quality of the pair is high. Second, the overall genetic quality of heterozygous offspring may be higher and thus increase the success of early development. Given the relatively weak relationships between measures of reproductive success and genetic similarity, further evaluations of these hypotheses are warranted. Another apparent cost of homozygosity is that more homozygous offspring were less likely to reach 25 days of age than less homozygous ones. Again this pattern may be explained by different effects. For example, more heterozygous individuals may cope better with pathogens [8]. Furthermore, homozygous offspring are more likely to be produced by genetically similar pairs that may invest less in chick rearing, thus reinforcing the deleterious effects of homozygosity on chick growth and survival. We also found that chick growth slowed as homozygosity increased, which may also explain the observed differences in survival: slower growing chicks may be in worse condition than faster growing chicks and are thus less efficient at defending themselves against pathogens and environmental stress. Interestingly, this effect was only evident in second hatched chicks (Figure 4), which might be explained by the fact these chicks suffer from more severe selective pressures due to sibling competition and brood reduction (unpublished data, see also [43]).

The use of multi-locus estimates of heterozygosity and genetic similarity has been criticized by authors arguing that global heterozygosity-fitness correlations may be driven by certain loci that are physically linked to fitness-affecting loci [35, 36, 44]. However, such local effects are expected to be weak in this population, since we found no locus that correlated with all components of fitness. We are also aware that the number of microsatellites we used (7–10) is low given that the number of loci needed to achieve an accurate estimation of individual global heterozygosity may be much higher [45]. However, a small number of microsatellite loci should diminish our capacity to detect any relationship between heterozygosity and fitness components, making our analyses conservative. The fact that we found negative effects of genetic similarity and homozygosity on different components of fitness (mate choice, hatching success, offspring growth and survival) despite our relatively small number of loci suggests that genetic similarity is costly, and therefore selected against, in this population.

Little is known about mate choice in kittiwakes, and the pattern we describe here may result from different mechanisms. It could be due to a passive process with, for example, genetically similar pairs being more likely to fail in their breeding attempt, which in turn would make them more prone to divorcing than successful pairs. In our population, this seems unlikely since divorce after reproductive failure is not systematic [25]. Furthermore, the correlation between genetic similarity of mates and reproductive success may not be strong enough to make divorces significantly correlated to genetic similarity. Alternatively, kittiwakes may choose to pair with individuals from different areas within a structured population, but we think this is unlikely. Genetic structuring is low at the scale of the whole North Atlantic kittiwake population [46, 47]). At a local scale, movements due to re-nesting of previously failed breeders into new areas are likely to diminish preexisting genetic structuring in this population. A third mechanism could be that prospective breeders may actively search for genetically dissimilar mates. The evolution of this strategy could be driven by the fitness costs of genetic similarity. Active mate choice would allow individuals to encounter a more adaptive mate faster than a passive process. However, experiments are needed to demonstrate that active mate choice has produced the observed pattern. Long-term monitoring and genotyping may also allow us to compare pairs that divorced with pairs that stayed together in relation to the genetic similarity of the first and second mate.

Active choice of genetically dissimilar mates has been demonstrated in mammals (e.g., humans, [48], mice [49] or fur seals, [50]) where genetic similarity is detecTable 2ecause of correlations between MHC alleles and body odors (humans, [51]; lemurs, [52]; mice [53]). Although the use of smell remains poorly known in birds, Antarctic prions (Pachiptila desolata) have been shown to recognize their mates through odors [54]. Genetically driven odors have not yet been shown in birds, but various aspects of social interactions in kittiwakes may allow them to recognize and choose their mates according to odor. Experiments are needed to determine the potential mechanisms for birds to estimate their relatedness to potential mates.


Black-legged kittiwakes pair with individuals that are genetically more dissimilar than expected by chance, a pattern that is consistent with our observation of the existence of potential costs of producing homozygous progeny, such as decreases in hatching success and survival. In such genetically monogamous species, mate choice therefore seems to be under selection. Thus, the genetic criteria of mate and extra-pair mate choice in various genetically polygamous species may also be a general feature in genetically monogamous species. Indeed, breeders of strictly monogamous species may experience the highest selective pressure to choose genetically dissimilar mates.


Study species and population

We conducted our study on Middleton Island (Gulf of Alaska, 58°25' N, 146°19' W, May-July 2003–2006). This island supports a large declining population of black-legged kittiwakes (from 166,000 birds in 1981 to fewer than 25,000 in 1999; [55]). We studied kittiwakes nesting on the ledges of an abandoned U.S. Air Force radar tower that has been modified to enable close observations and easy capture. The study plot is characterized by vertical walls and uniform nest spacing, with breeders nesting on wooden ledges built specifically for cliff-nesting seabirds [55].

Nest sites were observed twice daily from mid-May to late-August to assess individual attendance and reproductive success. We arrived at the field site too late to assess the arrival dates of adults. During each visit, we recorded the color-band combination of attending adults, the stage of nest building, the presence of eggs and chicks and the behavior of adults in the form of incubation, nest building and copulation, variables that we used to confirm pair identity. Each chick was marked at hatching and banded at 25 days old. Adults received five bands (one coded U.S. metal and four color-bands). We could thus precisely determine laying and hatching dates of focal pairs, as well as the hatching rank of chicks (recorded as A for first hatched and B for second hatched; brood size rarely exceeds two chicks). Copulation behavior was also monitored in order to determine the sex of every paired adult; for unpaired adults, the sex was identified by morphological values such as bill length, tarsus length, bill width and head+bill length [56]. We sampled blood from breeding and non-breeding individuals in 2003–2004 in order to analyze mating patterns. Results for mating patterns in 2005 were not analyzed because of other experiments performed in 2004 that affected mating patterns and divorce rates in subsequent years. Therefore, in 2005, we used our dataset to analyze the fitness costs of homozygosity in chicks hatched from pairs that were not manipulated in 2004 or 2005. In order to analyze the fitness costs of homozygosity, chicks of such pairs were blood-sampled at hatching and were weighed and measured (for wing and tarsus length) every 5 days from hatching to 25 days old. We did not remain at the field station long enough to record fledging success, which occurs between 35 and 40 days old. However, mortality is very low between 25 days old and fledging in this species, and thus survival at 25 days old is a reliable estimation of survival at fledging (unpublished data). Blood was kept in a preservation buffer solution (Longmire buffer, [57]) that allows the storage of samples without refrigeration.

Genetic analysis

DNA extraction

Whole genomic DNA was extracted from each blood sample using a "salting out" protocol described in [58]), modified by substituting Pronase E for Protease K and incubating at 37°C in the lysis phase, and substituting 0.7 volumes of 2-propanol in place of 2 volumes of ETOH in the DNA precipitation phase. Genomic DNA extractions were quantified using fluorometry and diluted to 50 ng/μL working solutions.

PCR and electrophoresis

Samples were genotyped at 10 microsatellite loci. Loci K6, K16, K31, K32, K67 and K71 were first described in the black-legged kittiwake [59], see Table 3); RBG20, RBG27, RBG29 and RBG39 were developed from the red-billed gull, Larus novaehollandiae scopulinus [60], see Table 3). Three additional loci (K56, [59], RBG18 and RBG13, [60]) were also tested but not used because the first two were found to be unreliable (see the quality controls below), and the third was monomorphic in our population.

Table 3 Summary of the ten microsatellite loci.

The forward primer in each primer pair was synthesized with a modified 19- to 20-bp universal tail (M13F, M13R or SP6) added to the 5' end of the oligonucleotide [61]. We used a complementary fluorescently labeled (IRD700 or IRD800) primer, identical to the specific tail used to modify the forward primer, to detect alleles at each loci. We carried out amplifications in a final volume of 10 μL that contained 50 ng DNA extract, 0.2 mM dNTPs, 0.1 mg BSA, 1× PCR buffer (Perkin Elmer Cetus I; PE Biosystems, Forest City, California), 10.0 pmoles forward and reverse unlabeled primers, 1.0 pmole fluorescently labeled primer, and 0.2 units of Taq polymerase. PCR reactions began at 94°C for 90 seconds, and continued with 40 cycles each of 94 C for 15 s, 50°C for 15 s and 74°C for 30 s.

Amplification products were separated on a 48-well 25-cm 6% polyacrylamide gel on a LI-COR 4200 LR automated sequencer, using Base ImagIR™ (LI-COR, Inc., Lincoln, Nebraska). Allele sizes for specific samples at each locus were determined relative to the M13 phage single nucleotide ladder. These samples were later used as internal size standards to score new genotypes, using Gene ImagIR™ 4.05 software (Scanalytics, Inc., Fairfax, Virginia). For quality control purposes, we reprocessed a minimum of 10% of the samples for all markers, and only kept markers giving reliable scores. Two markers (K56, [59], and RBG18 [60]) that gave inconsistent scores at each PCR were thus discarded.

Genetic Diversity and Tests of Equilibrium

Mean number of alleles (A) and observed and expected heterozygosities (HO and HE) were calculated in GENEPOP Version 3.1 [62]. We also used this program to test linkage disequilibria and deviation from Hardy-Weinberg equilibrium (Markov chain parameters: 10,000 dememorization steps, 100 batches, and 5,000 iterations per batch). These different tests guided final marker selection.

A total of 645 adults were genotyped in the Middleton population from 2003 to 2006. The number of alleles per locus varied from 6 to 35. After correcting for multiple tests, two loci appeared to be out of Hardy-Weinberg equilibrium: K16 (p < 0.0001) and RBG29 (p = 0.005). Because of its lower heterozygosity, K67 was not genotyped in all samples. Thus, for each analysis, we present the results obtained using all loci, or excluding K16, K67 and RBG29, which are designated as the "OHW" loci, for "out of Hardy-Weinberg" equilibrium. Only K31 and K32 were genetically linked (p < 0.05 after correcting for multiple tests). Since our genetic similarity analyses are based on bootstrapping on the multi-locus genotypes (such that linkages between alleles are conserved for every bootstrap) and not on alleles, both loci were kept. We also conducted every allelic-based analysis with all ten loci and without "OHW" loci and K31 (7 remaining loci).

Genetic similarity indices

The genetic similarity between two given pair members (x and y) was estimated by the probability for this pair of producing homozygous offspring (Phm). For each locus (l), this probability is equal to:

following Lynch & Ritland's [63]) notation, where s ij is a Boolean factor equal to 1 if alleles i and j are similar (i and j standing for any other letter), and 0 otherwise. An index based on such probabilities was first proposed by Mathieu et al. [64], but has seldom been used. Belkhir et al. [65] showed that when the number of alleles per loci is low, such indices have a low variance in their estimations of relatedness, a property relevant to our study in which some loci were not highly variable (see Table 3).

To obtain a value for Phm across all loci, a weighted average of the Phm xy (l) is calculated using the formula proposed by Mathieu et al. [64]:

where p l is the probability of an individual being homozygous by chance at locus l. Therefore Phm xy will be close to 1 if the pair has a high probability of producing a homozygous offspring (pair members are closely genetically similar), and close to 0 otherwise.

Individual heterozygosity

We used three different indices to estimate individual global heterozygosity: the direct heterozygosity H (proportion of heterozygous loci in a given individual), the standardized heterozygosity SH and the internal relatedness IR, all described in [66]). All these indices were calculated using the same program computed in Delphi (Borland Delphi 5.0, ©1983, 1989 Inprise Corporation; Additional file 2 and 3).

Monte-Carlo analyses

Genetic indices were performed using a program computed (Additional file 2 and 3) in Delphi (Borland Delphi 5.0, ©1983, 1989 Inprise Corporation) that was previously tested for no significant differences with Queller & Goodnight estimations obtained from the IDENTIX software [65]. This program allowed us to calculate Phm xy and, given the configuration of our datasets, it makes our Monte-Carlo analyses much easier than previously published programs.

For each given sample of individuals, we re-mated adults randomly according to their sex 10,000 times (i.e. 10,000 bootstraps). These bootstraps gave a set of values that we used as the distribution under the assumption of random mating (in further analyses, we only give the mean of the genetic indices of all these bootstraps). The observed values of the genetic indices were then tested for significant departure from the expected distribution obtained by the bootstraps. When bootstrapping, we used only adults seen in the considered year, and in each bootstrap we randomly made the same number of pairs that had been genotyped in the population using all live adults, as they were all considered as potentially available for pairing. Allelic frequencies were calculated for this sub-sample of individuals. For Hypothesis 2, that mates are genetically less similar than expected by chance, statistical significance was assessed by the proportion of the bootstraps having a mean genetic index lower than the observed mean. We statistically tested Hypothesis 3 that mates share more alleles than expected by chance, by assessing the proportion of the bootstrap values with a mean genetic similarity index higher than the observed mean (i.e. 1 minus the p-value calculated under Hypothesis 2). By bootstrapping on observed individuals and not on randomly created individuals generated via allelic frequencies, we might also expect to correct for biases due to linkage disequilibria or genotyping errors. Individuals with incomplete genotypes had thus the same probability of being included in observed pairs as in bootstrapped pairs. Thus, biases will be similar for observed and bootstrapped pairs.

General statistics

General statistical analyses were made using SAS® package (©SAS Institute Inc.1999, Cary, NC, USA). Correlations between the two indices of genetic similarity, or between male and female heterozygosity, were analyzed through standard linear models. To analyze the links between reproductive success (number of eggs laid and the number and proportion of eggs hatched) and genetic variables, we used a generalized linear model (allowing fitting standard linear models to discrete variables) since the dependent variables never showed more than three levels. For each of these three levels, the genetic variables and the residuals from the model were normally distributed. Similarly, we analyzed the correlations between chick survival and heterozygosity using generalized linear models, since both chick survival (alive or dead) and hatching rank (A or B) are binary effects. For each rank and each survival value, the residuals from the model procedure were normally distributed.

For the analyses of growth variables (body weight, tarsus or wing length), we have repeated measures of individuals, having measured each chick every 5 days. We therefore conducted mixed linear models with chick identity as a random parameter in order to account for non-independence. We made mixed linear models with tarsus length, wing length or body weight as dependent variables, chick identity as a random parameter, and chick age (A), age*age (A2), heterozygosity (H) and hatching rank (R) as explaining variables. We added A2 because growth of morphological traits is non-linear. For each variable and each heterozygosity index, we ran several models and kept only the one with the lowest AIC value. Given that the three growth variables are correlated, and given that the different indices of heterozygosity are also correlated, the structure of the model with the lowest AIC was always the same, in the form of Growth variable = A + H + R + A*H + A*R + H*R + A*H*R + A2 + Chick identity as a random parameter. We will only give the results obtained with the latter model, since the difference in AIC between this model and other models was always higher than 2. The significance of the interaction A*H was then tested in each rank separately.


  1. Burley NT, Foster V: Variation in female choice of mates: condition influences selectivity. Anim Behav. 2006, 72: 713-719. 10.1016/j.anbehav.2006.01.017.

    Article  Google Scholar 

  2. Johnsen A, Delhey K, Schlicht E, Peters A, Kempenaers B: Male sexual attractiveness and parental effort in Blue tits: a test of the differential allocation hypothesis. Anim Behav. 2005, 70: 877-888. 10.1016/j.anbehav.2005.01.005.

    Article  Google Scholar 

  3. Poesel A, Kunc HP, Foerster K, Johnsen A, Kempenaers B: Early birds are sexy: male age, dawn song and extrapair paternity in Blue tits, Cyanistes (formerly Parus) caeruleus. Anim Behav. 2006, 72: 531-538. 10.1016/j.anbehav.2005.10.022.

    Article  Google Scholar 

  4. Rantala MJ, Eriksson CJP, Vainikka A, Kortet R: Male steroid hormones and female preference for male body odor. Evolution and Human Behavior. 2006, 27: 259-269. 10.1016/j.evolhumbehav.2005.11.002.

    Article  Google Scholar 

  5. Neff BD, Pitcher TE: Genetic quality and sexual selection: an integrated framework for good genes and compatible genes. Molec Ecol. 2005, 14: 19-38. 10.1111/j.1365-294X.2004.02395.x.

    Article  CAS  Google Scholar 

  6. Mays HL, Hill GE: Choosing mates: good genes versus genes that are a good fit. TREE. 2004, 19: 554-559.

    PubMed  Google Scholar 

  7. Bonneaud C, Chastel O, Federici P, Westerdahl H, Sorci G: Complex MHC-based mate choice in a wild passerine. Proc R Soc B. 2006, 273: 1111-1116. 10.1098/rspb.2005.3325.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Penn DJ, Damjanovich K, Potts WK: MHC heterozygosity confers a selective advantage against multiple-strain infections. PNAS. 2002, 99: 11260-11264. 10.1073/pnas.162006499.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Westerdahl H, Waldenström J, Hansson B, Hasselquist D, Von Schantz T, Bensch S: Associations between malaria and MHC genes in a migratory songbird. Proc R Soc B. 2005, 272: 1511-1518. 10.1098/rspb.2005.3113.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Cohen LB, Dearborn DC: Great frigatebirds, Fregata minor, choose mates that are genetically similar. Anim Behav. 2004, 68: 1229-1236. 10.1016/j.anbehav.2003.12.021.

    Article  Google Scholar 

  11. Thuman KA, Griffith SC: Genetic similarity and the nonrandom distribution of paternity in a genetically highly polyandrous shorebird. Anim Behav. 2005, 69: 765-770. 10.1016/j.anbehav.2004.10.003.

    Article  Google Scholar 

  12. Olsson M, Madsen T, Nordby J, Wapstra E, Ujvari B, Wittsell H: Major histocompatibility complex and mate choice in sand lizards. Proc R Soc B. 2003, S1-S3. Supplement

  13. Freeman-Gallant CR, Meguerdichian M, Wheelwright NT, Sollecito SV: Social pairing and female mating fidelity predicted by restriction fragment length polymorphism similarity at the major histocompatibility complex in a songbird. Molec Ecol. 2003, 12: 3077-3083. 10.1046/j.1365-294X.2003.01968.x.

    Article  Google Scholar 

  14. Blomqvist D, Andersson M, Küpper C, Cuthill IC, Kis J, Lanctot RB, Sandercock BK, Székely T, Wallander J, Kempenaers B: Genetic simimarity between mates and extra-pair parentage in three species of shorebirds. Nature. 2002, 419: 613-615. 10.1038/nature01104.

    Article  CAS  PubMed  Google Scholar 

  15. Landry C, Garant D, Duchesne P, Bernatchez L: 'Good genes as heterozygosity': the major histocompatibility complex and mate choice in Atlantic salmon (Salmo salar). Proc R Soc B. 2001, 268: 1279-1285. 10.1098/rspb.2001.1659.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Tregenza T, Wedell N: Genetic compatibility, mate choice and patterns of parentage: Invited review. Molec Ecol. 2000, 9: 1013-1027. 10.1046/j.1365-294x.2000.00964.x.

    Article  CAS  Google Scholar 

  17. Trivers RL: Parental investment and sexual selection. Sexual selection and the descent of man, 1871–1971. Edited by: Campbell B. 1972, Chicago, IL: Aldine Publishing Company, 136-179.

    Google Scholar 

  18. Jones IL, Hunter FM: Mutual sexual selection in a monogamous seabird. Nature. 1993, 362: 238-239. 10.1038/362238a0.

    Article  Google Scholar 

  19. Foerster K, Delhey K, Johnsen A, Lifjeld JT, Kempenaers B: Females increase offspring heterozygosity and fitness through extra-pair matings. Nature. 2003, 425: 714-717. 10.1038/nature01969.

    Article  CAS  PubMed  Google Scholar 

  20. Rubenstein DR: Female extrapair mate choice in a cooperative breeder: trading sex for help and increasing offspring heterozygosity. Proc R Soc B. 2007, 274: 1895-1903. 10.1098/rspb.2007.0424.

    Article  PubMed Central  PubMed  Google Scholar 

  21. Helfenstein F, Tirard C, Danchin E, Wagner RH: Low frequency of extra-pair paternity and high frequency of adoption in Black-legged kittiwakes. Condor. 2004, 106: 149-155. 10.1650/7337.

    Article  Google Scholar 

  22. Coulson JC, Thomas CS: Mate choice in the kittiwake gull. Mate choice. Edited by: Bateson P. 1983, Cambridge: Cambridge University Press, 361-376.

    Google Scholar 

  23. Fairweather JA, Coulson JC: Mate retention in the kittiwake, Rissa tridactyla, and the significance of nest site tenacity. Anim Behav. 1995, 50: 455-464. 10.1006/anbe.1995.0259.

    Article  Google Scholar 

  24. Hatch SA, Roberts BD, Fadely BS: Adult survival of Black-legged kittiwakes Rissa tridactyla in a Pacific colony. Ibis. 1993, 135: 247-254. 10.1111/j.1474-919X.1993.tb02841.x.

    Article  Google Scholar 

  25. Naves LC, Monnat J-Y, Cam E: Breeding performance, mate fidelity, and nest site fidelity in a long-lived seabird: behaving against the current?. OIKOS. 2006, 115: 263-276. 10.1111/j.2006.0030-1299.14883.x.

    Article  Google Scholar 

  26. Bensch S, Hasselquist D, Von Schantz T: Genetic similarity between parents predicts hatching failure nonincestuous inbreeding in the Great reed warbler?. Evolution. 1994, 48: 317-326. 10.2307/2410095.

    Article  Google Scholar 

  27. Brown JL: A theory of mate choice based on heterozygosity. Behav Ecol. 1996, 8: 60-65. 10.1093/beheco/8.1.60.

    Article  Google Scholar 

  28. Kokko H, Ots I: When not to avoid inbreeding. Evolution. 2006, 60: 467-475.

    Article  PubMed  Google Scholar 

  29. Kempenaers B, Adriaensen F, Van Noordwijk AJ, Dhondt AA: Genetic similarity, inbreeding and hatching failure in blue tits: are unhatched eggs infertile?. Proc R Soc B. 1996, 263: 179-185. 10.1098/rspb.1996.0029.

    Article  Google Scholar 

  30. Casteele Van de T, Galbusera P, Schenck T, Matthysen E: Seasonal and lifetime reproductive consequences of inbreeding in the great tit Parus major . Behav Ecol. 2003, 14: 165-174. 10.1093/beheco/14.2.165.

    Article  Google Scholar 

  31. Coulson TN, Albon SD, Slate J, Pemberton JM: Microsatellite loci reveal sex-dependent responses to inbreeding and outbreeding in red deer calves. Evolution. 1999, 53: 1951-1960. 10.2307/2640453.

    Article  Google Scholar 

  32. Coltman DW, Slate J: Microsatellite measures of inbreeding: a meta-analysis. Evolution. 2003, 57: 971-983.

    Article  CAS  PubMed  Google Scholar 

  33. Hansson B: Marker-based relatedness predicts egg-hatching failure in great reed warblers. Conserv Gen. 2004, 5: 339-348. 10.1023/B:COGE.0000031143.51063.cb.

    Article  CAS  Google Scholar 

  34. Höglund J, Piertney SB, Alatalo RV, Lindell J, Lundberg A, Rintamäki PT: Inbreeding depression and male fitness in black grouse. Proc R Soc B. 2002, 269: 711-715. 10.1098/rspb.2001.1937.

    Article  PubMed Central  PubMed  Google Scholar 

  35. Lieutenant-Gosselin M, Bernatchez L: Local heterozygosity-fitness correlations with global positive effects of fitness in threespine stickleback. Evolution. 2006, 60: 1658-1668.

    Article  CAS  PubMed  Google Scholar 

  36. Tiira K, Laurila A, Enberg K, Piironen J, Aikio S, Ranta E, Primmer CR: Do dominants have higher heterozygosity? Social status and genetic variation in brown trout, Salmo trutta. Behav Ecol Sociobiol. 2005, 59: 657-665. 10.1007/s00265-005-0094-8.

    Article  Google Scholar 

  37. Jamieson IG, Taylor SS, Tracy LN, Kokko H, Armstrong DP: Why some species of birds do not avoid inbreeding: insights from New Zealand robins and saddlebacks. Behav Ecol. 2009, 20: 575-584. 10.1093/beheco/arp034.

    Article  Google Scholar 

  38. Roberts SC, Gosling LM: Genetic similarity and quality interact in mate choice decisions by female mice. Nat Gen. 2003, 35: 103-106. 10.1038/ng1231.

    Article  CAS  Google Scholar 

  39. Milinski M, Wedekind C: Evidence for MHC-correlated perfume preferences in humans. Behav Ecol. 2001, 12: 140-149. 10.1093/beheco/12.2.140.

    Article  Google Scholar 

  40. Richardson DS, Komdeur J, Burke T, Von Schantz T: MHC-based patterns of social and extra-pair mate choice in the Seychelles warbler. Proc R Soc B. 2005, 272: 759-767. 10.1098/rspb.2004.3028.

    Article  PubMed Central  PubMed  Google Scholar 

  41. Rätti O, Hovi M, Lundberg A, Tegelström H, Alatalo RV: Extra-pair paternity and male characteristics in the pied flycatcher. Behav Ecol Sociobiol. 1995, 37: 419-425. 10.1007/BF00170590.

    Article  Google Scholar 

  42. Westerdahl H: No evidence of an MHC-based female mating preference in great reed warblers. Molec Ecol. 2004, 13: 2465-2470. 10.1111/j.1365-294X.2004.02238.x.

    Article  CAS  Google Scholar 

  43. Braun BM, Hunt GLJ: Brood reduction in black-legged kittiwakes. Auk. 1983, 100: 469-476.

    Google Scholar 

  44. Hansson B, Westerberg L: On the correlation between heterozygosity and fitness in natural populations. Molec Ecol. 2002, 11: 2467-2474. 10.1046/j.1365-294X.2002.01644.x.

    Article  Google Scholar 

  45. Cooper G, Amos W, Bellamy R, Ruby Siddiqui M, Frodsham A, Hill AVS, Rubinsztein DC: An empirical exploration of the (δμ)2 genetic distance for 213 human microsatellite markers. Am J Hum Genet. 1999, 65: 1125-1133. 10.1086/302574.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  46. McCoy KD, Boulinier T, Tirard C: Comparative host-parasite population structures: disentangling prospecting and dispersal in the Black-legged kittiwake Rissa tridactyla . Molec Ecol. 2005, 14: 2825-2838. 10.1111/j.1365-294X.2005.02631.x.

    Article  CAS  Google Scholar 

  47. Friesen VL, Burg TM, McCoy KD: Mechanisms of population differenciation in seabirds. Molec Ecol. 2007, 16: 1765-1785. 10.1111/j.1365-294X.2006.03197.x.

    Article  CAS  Google Scholar 

  48. Wedekind C, Seebeck T, Bettens F, Paepke AJ: MHC-dependent mate preferences in humans. Proc R Soc B. 1995, 260: 245-249. 10.1098/rspb.1995.0087.

    Article  CAS  PubMed  Google Scholar 

  49. Penn DJ, Potts WK: Untrained mice discriminate MHC-determined odors. Physiol & Behav. 1998, 63: 235-243. 10.1016/S0031-9384(98)00052-3.

    Article  Google Scholar 

  50. Hoffman JI, Fordaca J, Trathan PN, Amos W: Female fur seals show active choice for males that are heterozygous and unrelated. Nature. 2007, 445: 912-914. 10.1038/nature05558.

    Article  CAS  PubMed  Google Scholar 

  51. McClintock M, Schumm P, Jacob S, Zelano B, Ober C: (in reply of) The MHC and body odors: arbitrary effects caused by shifts of mean pleasantness. Nature. 2002, 31: 237-238.

    CAS  Google Scholar 

  52. Knapp LA, Robson J, Waterhouse JS: Olfactory signals and the MHC: a review and a case study in Lemur catta . American Journal of Primatology. 2006, 68: 568-584. 10.1002/ajp.20253.

    Article  CAS  PubMed  Google Scholar 

  53. Yamazaki K, Beauchamp GK, Imai Y, Bard J, Phelan SP, Thomas L, Boyse EA: Odortypes determined by the major histocompatibility complex in germfree mice. PNAS. 1990, 87: 8413-8416. 10.1073/pnas.87.21.8413.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  54. Bonadonna F, Nevitt GA: Partner-specific odor recognition in an Antarctic seabird. Science. 2004, 306: 835-10.1126/science.1103001.

    Article  CAS  PubMed  Google Scholar 

  55. Gill VA, Hatch SA: Components of productivity in black-legged kittiwakes Rissa tridactyla : response to supplemental feeding. J Av Biol. 2002, 33: 113-126. 10.1034/j.1600-048X.2002.330201.x.

    Article  Google Scholar 

  56. Jodice PGR, Lanctot RB, Gill VA, Roby DD, Hatch S: Sexing adult black-legged kittiwakes by DNA, behavior, and morphology. Waterbirds. 2000, 23: 405-415.

    Article  Google Scholar 

  57. Longmire JL, Lewis AK, Brown NC, Buckingham JM, Clark LM, Jones MD, Meincke LJ, Meyne J, Ratliff RL, Ray FA, et al: Isolation and molecular characterization of a highly polymorphic centromeric tandem repeat in the family Falconidae. Genomics. 1988, 2: 14-24. 10.1016/0888-7543(88)90104-8.

    Article  CAS  PubMed  Google Scholar 

  58. Medrano JF, Aasen E, Sharrow L: DNA extraction from nucleated red blood cells. Biotechniques. 1990, 8: 43-

    CAS  PubMed  Google Scholar 

  59. Tirard C, Helfenstein F, Danchin E: Polymorphic microsatellites in the Black-legged kittiwake Rissa tridactyla . Molec Ecol Notes. 2002, 2: 431-433. 10.1046/j.1471-8286.2002.00258.x.

    Article  CAS  Google Scholar 

  60. Given AD, Mills A, Baker J: Isolation of polymorphic microsatellite loci from the red-billed gull (Larus novaehollandiae scopulinus) and amplification in related species. Molec Ecol Notes. 2002, 2: 416-418. 10.1046/j.1471-8286.2002.00261.x.

    Article  CAS  Google Scholar 

  61. Oetting WS, Lee HK, Flanders DJ, Wiesner GL, Sellers TA, King RA: Linkage analysis with multiplexed short tandem repeat polymorphisms using infrared fluorescence and M13 taied primers. Genomics. 1995, 30: 450-458. 10.1006/geno.1995.1264.

    Article  CAS  PubMed  Google Scholar 

  62. Raymond M, Rousset F: GENEPOP (ver. 1.2): a population genetics software for exact test and ecumenicism. Heredity. 1995, 86: 248-249.

    Google Scholar 

  63. Lynch M, Ritland K: Estimation of pairwise relatedness with molecular markers. Genetics. 1999, 152: 1753-1766.

    PubMed Central  CAS  PubMed  Google Scholar 

  64. Mathieu E, Autem M, Roux M, Bonhomme F: Epreuves de validation dans l'analyse de structures génétiques multivariées: comment tester l'équilibre panmictique?. Revue de Statistique Appliquée. 1990, 38: 47-66.

    Google Scholar 

  65. Belkhir K, Castric V, Bonhomme F: IDENTIX, a software to test for relatedness in a population using permutation methods. Molec Ecol Notes. 2002, 2: 611-614. 10.1046/j.1471-8286.2002.00273.x.

    Article  Google Scholar 

  66. Amos W, Worthington Wilmer J, Fullard K, Burg TM, Croxall JP, Bloch D, Coulson T: The influence of parental relatedness on reproductive success. Proc R Soc B. 2001, 268: 2021-2027. 10.1098/rspb.2001.1751.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references


We thank Kevin Sage and Judy Gust for their help in the lab, as well as several assistants for their help in the field. Karen McCoy, Donald Blomqvist and Bart Kempenaers gave valuable comments on earlier drafts. Experiments were carried out in accordance with American regulations and under permits from the U.S. Fish and Wildlife Service (MB789758-2) and the U.S. Department of Agriculture (43860). This study was financed by a four year grant from the French Polar Institute Paul-Emile Victor ("Programme Arctique 429" 2004–2007), a France-Austria PICS (Programme International de Coopération Scientifique), and by the U.S. Geological Survey, Alaska Science Center. This project was part of the French GDR2155 "Écologie comportementale".

Author information

Authors and Affiliations


Corresponding author

Correspondence to Hervé Mulard.

Additional information

Authors' contributions

HM carried out the genetic analyses, monitored the birds on Middleton Island and drafted the manuscript with RHW. ED and RHW conceived the study and participated in its design and coordination. JFW and FH helped in monitoring the birds and in the conception of the study in the field. SLT and AMR participated in the microsatellite genotyping and scoring. SAH is responsible of the long-term monitoring on Middleton Island and gave access to the study area. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Table 2with effect size of heterozygosity. Extended version of the Table 2, including mean effect size of heterozygosity for all analyses. (DOC 91 KB)


Additional file 2: Code of the program Nausicaa. This file contains the complete code of the program used for calculation of genetic indices and Monte-Carlo analyses. (TXT 172 KB)


Additional file 3: GNU-GPL Licence. GNU-GPL Licence for the program Nausicaa contained in additional file 2. (TXT 32 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Mulard, H., Danchin, E., Talbot, S.L. et al. Evidence that pairing with genetically similar mates is maladaptive in a monogamous bird. BMC Evol Biol 9, 147 (2009).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: