Reed bunting females increase fitness through extra-pair mating with genetically dissimilar males

Stefan M Suter, Martin Keiser, Raoul Feignoux, Dietrich R Meyer

Abstract

Females of many socially monogamous species accept or even actively seek copulations outside the social pair bond. As females cannot increase the number of offspring with promiscuous behaviour, the question arises why they engage in extra-pair mating. We used microsatellite data to determine paternity, heterozygosity and genetic relatedness in the reed bunting (Emberiza schoeniclus), a species with high levels of extra-pair paternity (EPP). We found that extra-pair young (EPY) were more heterozygous than within-pair young (WPY). The high heterozygosity of the EPY resulted from a low genetic similarity between females and their extra-pair mates. EPY were heavier and larger when compared with their maternal half-siblings shortly before they left the nest. Recapture data indicated a higher fledgling survival of EPY compared with WPY. Our data suggest that reed bunting females increase the viability of their offspring and thus fitness through extra-pair mating with genetically dissimilar males.

Keywords:

1. Introduction

Promiscuity, whereby males and females mate with other than their social partner, is common in socially monogamous bird species (Griffith et al. 2002). Males may initiate extra-pair copulations (Westneat & Stewart 2003) but females have pre- and post-copulatory control over fertilization (Andersson 1994; Birkhead & Pizzari 2002). While males can increase the number of offspring via extra-pair fertilizations, it is less obvious how females could gain fitness due to promiscuous behaviour. It has been proposed that females increase offspring heterozygosity and, thus, the viability of their offspring through extra-pair fertilizations (Brown 1997). Evidence has been found in several species (Foerster et al. 2003; Masters et al. 2003; Bishop et al. 2007). However, other studies do not support these findings (Schmoll et al. 2005; Smith et al. 2005). A high heterozygosity reduces the risk that recessive deleterious alleles are expressed and prevents other negative effects of inbreeding, which would lead to reduced fitness. However, females may not only engage in extra-pair copulations to avoid inbreeding but also benefit from mating with more dissimilar mates (Amos et al. 2001). In that case, high heterozygosity in the offspring will be advantageous owing to potentially useful gene products. A higher heterozygosity at the major histocompatibility complex (MHC) does, for example, lead to an increased resistance against diseases (McClelland et al. 2003). In bluethroats (Luscinia s. svecica), Johnsen et al. (2000) found that extra-pair young (EPY) have a higher immunocompetence when compared with within-pair young (WPY). Females can maximize the heterozygosity of their offspring if they mate to males that are either highly heterozygous or genetically different from themselves (Masters et al. 2003).

Females that engage in extra-pair copulations may envisage different costs such as the loss of the social mate's investment in the offspring (Dixon et al. 1994), enhanced risk of sexually transmitted diseases, injury and predation. Therefore, females should engage in extra-pair mating only if it leads to substantial benefits. In a meta-analysis, Arnqvist & Kirkpatrick (2005) suggest that infidelity may not be adaptive for females because costs caused by depressed paternal investment outweigh indirect genetic benefits, but see Griffith (2007).

In this study we show that EPY were more heterozygous, heavier and larger when compared with WPY in a population of reed buntings (Emberiza schoeniclus) in Switzerland. Additionally we found that extra-pair fledglings survive at a higher rate when compared with within-pair fledglings. The observed extra-pair fertilizations resulted from matings between genetically dissimilar mates. Our findings suggest that reed bunting females increase the viability of their offspring through extra-pair mating.

2. Material and methods

(a) Field methods

The study was conducted on a reed bunting population in the Grande Cariçaie at the southeast shore of the lake of Neuchâtel in Switzerland for 6 years (during 2000–2002 and 2004–2006). Adult birds were caught with mist nets either at the nest site or (for some males) with additional help of playback within the territory. For individual identification, the birds were ringed with three coloured plastic rings and a numbered aluminium ring from the Swiss Ornithological Institute. From all birds we took 5–50μl blood from the brachial or leg vein for paternity analyses. Blood samples were collected in a 70 μl capillary tube and put on ice in the field and stored in a freezer at −18°C on the same day. Nestlings were bled from day 2 onwards and measured and ringed between days 6 and 11 after hatching (day 0=hatching day). Weight was measured with a Pesola balance (to the nearest 0.1 g) and tarsus length with a calliper (to the nearest 0.1 mm). Dead nestlings and eggs that failed to hatch were collected. Nests were monitored daily to determine hatch dates. However, some nests were found after hatching. We confirmed the social parents of a brood while they were feeding the young either by direct observation or with a colour camera installed near the nest.

Some of the birds that had been ringed as nestlings were recaptured later in the season from July to October. In the analysis, we used only data from birds that were recaptured at least 30 days after they had been ringed at the nest, thus, they were independent from their parents at the time of recapture. Birds were recaptured by our own research group at the study site at the end of the breeding season or by different other research groups during the autumn migration. All birds were recaptured along the southeast shore of the lake of Neuchâtel (0–25 km) except two of which one was caught at La Touvière (Geneva; 73 km) and the other at Romans-sur-Isère (France; 250 km).

(b) Molecular methods

DNA was extracted from blood and tissue samples using a Promega protocol during 2000–2002 and Peqgold blood DNA isolation kit (Peqlab) during 2004–2006.

For paternity analyses, we used fluorescently labelled primers for six different variable microsatellite loci: Escμ1; Escμ3; Escμ4; Escμ6 (Hanotte et al. 1994); Pdoμ5 (Griffith et al. 1999); and Ppi2 (Martinez et al. 1999). They were amplified by PCR and products were run on an ABI PRISM 310 Genetic Analyzer (Applied Biosystems Instrument). The alleles were determined using DNA fragment analysis software (Genescan v. 3.1). The microsatellite loci did not deviate from Hardy–Weinberg equilibrium except Escμ3 indicating null alleles at this locus. The combined exclusion probability of the six microsatellite was higher than 0.995 for the first parent and 0.999 for the second parent using the program Cervus (v. 3.0; Kalinowski et al. 2007).

To sex the young we used the primers P2 (reverse) and P8 (forward) that anneal to conserved regions and amplify across introns that differ in length between the CHD-W and the CHD-Z gene (Griffiths et al. 1998).

(c) Statistical analyses

We used a general linear mixed model to test the effect of paternity on tarsus length and weight of nestlings from broods with mixed paternity. Sex was included as a fixed factor in the model because males are heavier and larger than females. The age of the brood was included as a linear covariate because broods were measured between days 6 and 8 after hatching. Brood identity was included as a random factor in the model to take into account the seasonal differences (e.g. food abundance) and the quality of the social parents.

Standardized individual heterozygosity (SIH) was calculated as the proportion of heterozygous loci divided by the mean heterozygosity of the typed loci in the population (Coltman et al. 1999). Genetic similarity between individuals was calculated using the program ML-Relate that takes into account the presence of null alleles (Kalinowski et al. 2006). A permutation test was conducted to test whether females mate randomly to their social mate with regard to the genetic relatedness. The observed number of social pairs was randomly formed 1000 times using all adult birds present at the site. The observed mean relatedness of the social pairs was then compared to the distribution of the means that had been generated under the assumption of random mating. To test whether the extra-pair fertilizations resulted from matings between genetically dissimilar mates, we compared the genetic relatedness between the female and the actual extra-pair sire or sires to the mean genetic relatedness between the female and the four nearest residents that did not sire EPY in the focal female's nest. Since almost all EPY were sired by males that have their territories in the neighbourhood of the female's nest (S. Suter 2007, unpublished data), this procedure allows a comparison between actual and potential extra-pair mates (Masters et al. 2003).

If there was more than one extra-pair father within the same brood, then the mean of all extra-pair sires was used in the analysis. Pairs where social partners changed between years were treated as independent. Statistical analysis was performed using JMP v. 5.0.1, R v. 2.4.1 and Excel. Non-parametric tests were used in case of non-normal distribution of the data. All tests are two tailed with a significance level of p<0.05.

3. Results

(a) Frequency and distribution of extra-pair paternity

We genotyped 915 offspring, of which 65 were fledged young, 786 nestlings and 64 embryos from eggs that did not hatch. The genetic father was determined for 835 (91%) of the young and the social father for 801 (88%). Overall there were 486 (61%) WPY and 315 (39%) EPY. The proportion of EPY in different years varied between 33 and 46% and the proportion of broods that contained EPY varied between 54 and 72%. No significant differences between years were found (G=8.23 and 4.09, d.f.=5, p=0.14 and 0.54, respectively). We could determine the sex of 792 out of the 801 young where the social father was known. Overall, EPY were not more likely to be males: there were 47% (224/481) WPY males and 51% EPY males (159/311; Fisher's exact test, p=0.22).

The EPY were not equally distributed among broods; there were 36 entire extra-pair broods, 98 mixed broods and 79 broods containing only WPY. In each year, the distribution of EPY among broods differed from what is expected under a binomial distribution (G-tests, all p<0.01), except in the year 2001 (G=4.92, d.f.=5, p=0.42). No egg dumping was detected as the social mother of a brood always corresponded to the genetic mother. A summary of the microsatellite loci used for paternity analysis can be found in table 3 of electronic supplementary material.

(b) Paternity and offspring heterozygosity

SIH of the 915 offspring varied between 0.28 and 1.18 (mean SIH±s.d., 1.00±0.17). EPY were more heterozygous when compared with WPY (mean SIH±s.e., EPY n=315, 1.02±0.01; WPY n=486, 0.98±0.01, Wilcoxon rank–sum test, Z=3.55, p=0.0004).

However, if we look only at mixed broods, the difference between maternal half-siblings is not significant (Wilcoxon signed-rank test, Z=0.55, n=57, p=0.58). The high overall SIH of EPY resulted from high SIH of the young that came from entire extra-pair broods and the low overall SIH of WPY resulted from low SIH of the young that came from entire within-pair broods (figure 1). In a pairwise comparison, the mean SIH of extra-pair offspring showed a tendency to be higher compared with that of within-pair offspring a female produced (Wilcoxon signed-rank test, Z=1.67, n=67, p=0.09). A paired test of the mean SIH between offspring of entire within- (mean SIH ± s.e.: 0.99±0.03) and entire extra-pair broods (1.05±0.03) of the same female revealed a significant difference (Wilcoxon signed-rank test, Z=2.03, n=13, p=0.04). Allele frequencies of different microsatellite loci are presented in table 4 of the electronic supplementary material.

Figure 1

Mean standardized individual heterozygosity (SIH) of female's offspring over different brood categories and paternities (mean values at the top, number of females at the bottom of the bars). Different brood categories are: a, entire within-pair broods; b, mixed broods; c, entire extra-pair broods. Mixed broods are separated by paternity: within-pair young (WPY) and extra-pair young (EPY). Offspring of a female may occur in more than one brood category, in total 117 different mothers are included.

(c) Paternity and offspring size parameters

Paternity had a significant effect on weight and tarsus length, controlling for sex, age and brood identity in a general linear mixed model (table 1). EPY were on average 0.5 g (3.3%) heavier when compared with WPY, and male nestlings were on average 0.8 g (5.8%) heavier compared with female nestlings. The effect of paternity on tarsus length was greater than the effect of sex. The tarsi of EPY were on average 0.24 mm (1.3%) longer than that of WPY, and the tarsi of male nestlings were on average 0.15 mm (0.8%) longer than that of female nestlings. Only mixed broods where the size parameters had been measured between days 6 and 8 after hatching were included in the analyses.

View this table:
Table 1

General linear mixed model to test the effect of paternity on weight and tarsus length in nestlings from 83 broods with mixed paternity. (Extra-pair young (n=163) were heavier and they had longer tarsi compared with within-pair young (n=165). The different factors in the model are: paternity; sex; and age. Brood identity had a significant effect and was included as a random factor in the model. Interactions between factors were not significant.)

(d) Paternity, offspring survival and recapture data

Overall there was no difference in hatching success between WPY (450/486) and EPY (299/315; Fisher's exact test p=0.24) or in fledging success WPY (385/486), EPY (264/315; Fisher's exact test p=0.12). A pairwise comparison in mixed broods revealed no differences in mortality from egg to fledging between maternal half-siblings (Wilcoxon signed-rank test, Z=0.54, n=18, p=0.58).

However, a larger percentage of EPY compared with WPY was recaptured later in the season (figure 2). This pattern is consistent over the 6 years. A pairwise comparison in mixed broods revealed a significantly higher recapture rate for EPY compared with WPY (Wilcoxon signed-rank test, Z=2.31, n=13, p=0.02). Also a pairwise comparison between within- and extra-pair offspring of the same female shows that more EPY were recaptured (Wilcoxon signed-rank test, Z=2.03, n=22, p=0.05).

Figure 2

The percentage of fledglings recaptured in different years. In all years a larger percentage of EPY (grey bars) compared with WPY (white bars) was recaptured. Pooled data over all years show that more EPY were recaptured (WPY 3.4% (13/385), EPY 7.6% (20/264); Fisher's exact test p=0.02). The numbers above bars are the number of ringed nestlings.

Birds that fledged early in the season were not more likely to be recaptured (logistic regression: Χ12=0.92, p=0.34). EPY did not hatch earlier in the season compared with WPY (mean relative hatching date±s.e., WPY n=450, 22.25±0.86; EPY n=299, 21.74±1.05, Wilcoxon rank–sum test, Z=−0.09, p=0.93). The relative hatching date of a brood is the difference in days between the hatching date of the focused brood and the hatching date of the first brood that hatched in a given year. No difference in recapture rate was recorded between male and female fledglings (recaptured by sex: males 5.03% (16/318), females 5.14% (17/331); Fisher's exact test p=1).

(e) Adult heterozygosity and extra-pair paternity

All tests within single year led to the same conclusions. Therefore, only results of pooled data are shown here. Males that sired EPY were not more heterozygous compared with males that did not (mean SIH±s.e. of males that sired: no EPY n=67, 1.01±0.02, EPY n=58, 1.02±0.02; Wilcoxon rank–sum test, Z=0.29, p=0.78). A paired test revealed no difference in heterozygosity between social males and their cuckolders (Wilcoxon signed-rank test, Z=−0.02, n=96, p=0.98) nor were extra-pair males that gained entire paternity in a brood more heterozygous compared with the corresponding social males (Wilcoxon signed-rank test, Z=0.25, n=31, p=0.80). There was no difference in heterozygosity between females that produced EPY and females that produced only WPY (mean SIH±s.e. of females with: no EPY n=35, 1.02±0.03, EPY n=81, 1.02±0.02; Wilcoxon rank–sum test, Z=−0.40, p=0.69).

(f) Relatedness and extra-pair paternity

A higher genetic relatedness between parents led to a lower heterozygosity in the offspring (rs=−0.20, n=258, p=0.002). Females did not choose genetically dissimilar males as breeding partners as the mean relatedness of females to their social mates was not different from that of randomly drawn pairs (table 2a). The relatedness of a female to her social mate did not predict the likelihood of having extra-pair offspring within a breeding season. For this analysis, we pooled the data from all years and included each pair only once (binomial logistic regression: Χ12=1.52, p=0.22, n=145). Successful extra-pair mates were genetically more dissimilar to the females than nearby potential unsuccessful extra-pair mates (table 2b). A test with pooled data over all years where each female was only included once revealed a significant difference (Wilcoxon signed-rank test, Z=3.84, n=73, p=0.0003).

View this table:
Table 2

Mean genetic relatedness±s.e. (n) between females and males. ((a) Comparison between female relatedness to all theoretically possible mates and female relatedness to their social mates (p values from permutation tests). (b) Comparison between female relatedness to potential extra-pair (ep) mates and female relatedness to actual extra-pair mates. (p values from Wilcoxon signed-rank test).)

In a pairwise comparison, extra-pair mates showed a tendency to be less related to females than the social mates (Wilcoxon signed-rank test, Z=1.77, n=83, p=0.07). The relatedness between female and social mate was not related to the percentage of EPP in a brood (rs=0.07, n=145, p=0.39). To test whether the percentage of EPP a brood contains depends on the relatedness between female and the extra-pair mate, we used one randomly chosen brood per female where a single male sired all EPY. The percentage of EPP a brood contained did not depend on the genetic relatedness between female and extra-pair mate (rs=0.16, n=73, p=0.16). The difference of genetic relatedness between female to the social mate and female to the extra-pair mate did not correlate with the amount of EPP (rs=0.05, n=73, p=0.68).

4. Discussion

Female fitness gains of extra-pair fertilizations could be assessed from different survival rates between EPY and WPY. Yet, dispersal makes it hard to get any long-term data on the survival of young birds in order to investigate differences based on paternity, but see Hasselquist et al. (1996), Lubjuhn et al. (1999) and Schmoll et al. (2003). A higher survival of EPY from hatching to fledging has been found in blue tits (Parus caeruleus; Kempenaers et al. 1997; Charmantier & Perret 2004). However, differences in nestling mortality may be less important compared with that in fledgling mortality. The time when a young bird leaves the nest until it becomes independent from its parents is known to be a critical phase in the life of a bird, and high mortality during this period has been reported (Ringsby et al. 1998; Naef-Daenzer et al. 2001). In our study, we assume that a within-pair fledgling that survived was as likely to be recaptured as an extra-pair fledgling. From the higher recapture rate of extra-pair fledglings, we suggest that they are more viable and therefore survive better compared with within-pair fledglings. Alternatively it might be argued that different dispersal patterns between EPY and WPY explain the higher recapture rate of EPY. Different migration patterns between males and females exist in the reed bunting (Schmitz & Steiner 2006). However, in our study we recaptured similar percentages of young males and females and the sex ratio did not differ between EPY and WPY.

EPY were heavier and had longer tarsi compared with their maternal half-siblings shortly before they left the nest. Weight has been shown to be important for fledgling survival (Magrath 1991; Naef-Daenzer et al. 2001). Longer tarsi of EPY compared with WPY have also been found in a population of reed buntings in The Netherlands (Bouwman et al. 2007) as well as in savannah sparrows (Passerculus sandwichensis; Freeman-Gallant et al. 2006). The tarsus was not full grown when we took measurements from the nestlings. It is therefore not a measure of the final size but an indicator of the development stage of a young. Fast development and thus the ability to leave the nest early are important in ground breeding birds owing to the high predation risk at nest sites. From the differences in the size parameters between EPY and WPY, we suggest that EPY develop faster and are therefore more viable. However, the size differences between extra- and within-pair nestlings could be based on other than genetic factors. The first eggs laid in a clutch might be more likely to contain EPY and therefore hatch earlier (Krist et al. 2005), but see Westneat et al. (1995), Whittingham et al. (2003) and Barber & Robertson (2007). Another alternative explanation would be that non-genetic maternal effects like increased female investment in the nutrient or hormonal content of the eggs that contain EPY cause these size differences (Gil et al. 1999; Cunningham & Russell 2000).

Females did choose their social partners randomly with regard to genetic relatedness. The choice of a social partner might be based on resources like territory quality and potential paternal care (Andersson 1994). Furthermore, female's choice for a social partner is restricted because males can already be paired to other females. To become a secondary female is disadvantageous owing to the reduced paternal assistance at nest (Dixon et al. 1994). A better strategy would be to pair with an unpaired male to assure the help at the nest and then go for extra-pair copulations to correct for low genetic quality or low compatibility of the social partner. In our study the genetic relatedness between breeding partners did not predict the occurrence of EPP, but see Blomqvist et al. (2002) and Freeman-Gallant et al. (2006). For the females even the choice of an extra-pair mate can be restricted. The social male has an interest to assure his paternity and could achieve it through intense mate guarding and frequent copulations.

We do not know whether reed bunting females are able to estimate genetic relatedness or compatibility of a male. Our finding that successful extra-pair mates are less related to females than potential but unsuccessful extra-pair mates is based on the observed extra-pair fertilizations. As we do not have any observational data on extra-pair copulations, we do not know whether females engage selectively in extra-pair copulations with genetically dissimilar males. In humans, mice and fishes, females were able to estimate genetic compatibility through odours in studies investigating mating preferences for genes of the MHC (Wedekind et al. 1995; Penn & Potts 1998; Reusch et al. 2001). Whether odour is important for mate choice in birds has not been tested yet although it is known that birds have the necessary olfactory ability (Roper 1999; Bonadonna & Nevitt 2004; Mennerat et al. 2005). It has been shown that song plays a role in kin recognition in birds (Komdeur & Hatchwell 1999; Sharp et al. 2005), therefore song could be a potential trait on which females estimate the genetic relatedness of a male.

Cryptic female choice could play an important role in the fertilization process (Birkhead & Pizzari 2002). Females could simply eject sperm of less desirable mates (Davies et al. 1996; Pizzari & Birkhead 2000) or the selection of sperm could take place within the female's reproductive tract (Wilson et al. 1997). Artificial insemination experiments have to be done to find out more about possible sperm selection within the female's reproductive tract.

Studies on reed bunting populations in Norway and The Netherlands did not reveal differences in heterozygosity between EPY and WPY (Kleven & Lifjeld 2005; Bouwman et al. 2007). However, the authors were focusing on mixed broods. In our population, there was likewise no significant difference in heterozygosity between maternal half-siblings of mixed broods. The higher heterozygosity of EPY in our study is based on the high heterozygosity of EPY that came from entire extra-pair broods. This raises the question whether the percentage of EPP in a brood depends on the genetic relatedness between the female and the extra-pair mate. We could not establish such a correlation. It is possible that an existing relationship is hidden by confounding factors which influence female ability to engage in extra-pair copulations, such as weather conditions (Bouwman & Komdeur 2006) or paternity assurance behaviour (Marthinsen et al. 2005).

The two main explanations for indirect genetic benefits of female infidelity are the good genes hypothesis and the compatible genes hypothesis, reviewed by different authors (Jennions & Petrie 2000; Griffith et al. 2002; Mays & Hill 2004; Neff & Pitcher 2005). In the reed bunting, it has been shown that mainly old males gain extra-pair paternity (Bouwman et al. 2007). The fact that old males have shown their ability to survive might be an indicator of their good genes (Kokko & Lindstrom 1996; Brooks & Kemp 2001). Our findings that EPY are more heterozygous and that this is based on a low relatedness between females and their extra-pair mates support the genetic compatibility hypothesis. However, whether dissimilarity equals compatibility can be questioned (Mays & Hill 2004; Puurtinen et al. 2005). Mate choice based on dissimilarity maintains or even increases genetic diversity within populations contrary to mate choice based on good genes. Both good and compatible genes would lead to higher viability in the offspring. Therefore, female fitness gains of extra-pair fertilizations can be based on both mechanisms because they are not mutually exclusive (Colegrave et al. 2002; Neff & Pitcher 2005).

In conclusion, our findings show that females can increase the quality of their offspring and thus fitness through extra-pair mating. Our study provides evidence against the meta-analysis of Arnqvist & Kirkpatrick (2005). The positive indirect selection may outweigh the cost of extra-pair copulation behaviour and therefore infidelity could be adaptive for females in the reed bunting population we studied.

Acknowledgements

Access to nature reserves was granted by the cantons of Fribourg and Vaud. The birds were caught under licence of the Federal Office for the Environment and blood samples were taken under licence of the cantonal ethics committee on animal experiments.

We thank Joanna Bielanska, Lukas Escher, David Ermacora, Vanessa Heierli, Manuel Lingg, Stephanie Michler, Nadia Rieille, Sabine Röthlin-Spillman, Ludivine Strambini, Elham Tarbush and Karin Safi-Widmer for their help in the field and in the laboratory. We thank Tadeusz Kawecki and three anonymous referees for their comments on the manuscript. We also thank all the people who caught birds and the Swiss Ornithological Station for information about recaptured birds. Fieldwork was coordinated with Michel Antoniazza from the Groupe d'étude et de gestion de la Grande Cariçaie (GEG).

Footnotes

References

View Abstract