- Research article
- Open Access
Rainfall-driven sex-ratio genes in African buffalo suggested by correlations between Y-chromosomal haplotype frequencies and foetal sex ratio
BMC Evolutionary Biology volume 10, Article number: 106 (2010)
The Y-chromosomal diversity in the African buffalo (Syncerus caffer) population of Kruger National Park (KNP) is characterized by rainfall-driven haplotype frequency shifts between year cohorts. Stable Y-chromosomal polymorphism is difficult to reconcile with haplotype frequency variations without assuming frequency-dependent selection or specific interactions in the population dynamics of X- and Y-chromosomal genes, since otherwise the fittest haplotype would inevitably sweep to fixation. Stable Y-chromosomal polymorphism due one of these factors only seems possible when there are Y-chromosomal distorters of an equal sex ratio, which act by negatively affecting X-gametes, or Y-chromosomal suppressors of a female-biased sex ratio. These sex-ratio (SR) genes modify (suppress) gamete transmission in their own favour at a fitness cost, allowing for stable polymorphism.
Here we show temporal correlations between Y-chromosomal haplotype frequencies and foetal sex ratios in the KNP buffalo population, suggesting SR genes. Frequencies varied by a factor of five; too high to be alternatively explained by Y-chromosomal effects on pregnancy loss. Sex ratios were male-biased during wet and female-biased during dry periods (male proportion: 0.47-0.53), seasonally and annually. Both wet and dry periods were associated with a specific haplotype indicating a SR distorter and SR suppressor, respectively.
The distinctive properties suggested for explaining Y-chromosomal polymorphism in African buffalo may not be restricted to this species alone. SR genes may play a broader and largely overlooked role in mammalian sex-ratio variation.
Sex-ratio (SR) distorters are genes that cause an unequal representation of the sex chromosomes among the offspring through differences in number, quality or function of X- and Y-bearing spermatozoa. These SR distorters are typically located on one of the sex chromosomes. They act during spermatogenesis by distorting the meiosis in their own favour (meiotic drive) or by impairing the function or increasing the mortality of the opposite-sex gametes, which gives the SR distorters a selective advantage [1–3]. Because of this selective advantage, newly mutated SR distorters can quickly invade a population. In principal, SR distorters may exert their effect not only during spermatogenesis but also post-copulation in the female genital tract, by influencing sperm transport, sperm survival or fertilization capability. SR distorters are assumed to be selfish and non-adaptive for their carriers [1, 4, 5].
SR distorters exert strong selection pressure for suppressors of the distortion, which are genes occurring on the opposite sex chromosome or on one of the autosomal chromosomes, that neutralize the action of the distorter and often lead to a genetic arms race between the two types of SR genes . A stable gene polymorphism is possible when SR-distorter carriers have a lower fitness, thereby preventing the driving distorter alleles from going to fixation . (A SR suppressor does not prevent SR-distorter fixation since the SR distorter still has a selective advantage compared to the wild-type allele, given that it is only being neutralized when co-occurring with a SR suppressor in the same individual.) Most examples of SR distorters are found in Drosophila spp. . Their discovery required sophisticated crossing schemes, often between individuals from different populations or incipient species . Although sex chromosomes with SR distorting elements are known in Drosophila spp. since 1925, the identity of causal genes remained unknown until the late 1990s because the SR distorting elements are usually associated with complex chromosome inversions, impeding their genetic mapping [7–12]. Nevertheless, SR distorters are thought to be relatively common, even in mammals, but are rarely observed, probably because of rapid gene fixation and SR suppressors that mask their activities [4, 5]. (When a SR distorter and SR suppressor are co-occurring in a population both will go to fixation, unless they have a lowered fitness.)
Stable presence of multiple alleles (per locus) or haplotypes (combinations of alleles at different loci along a chromosome inherited en block) on the haploid Y-chromosome in the absence of population subdivision or metapopulation structure only seems possible with negative frequency-dependent selection, since otherwise the fittest allele or haplotype would inevitably sweep to fixation [1, 13–15]. SR distorters generally reduce the fitness of their bearer by a reduction of the sperm count making negative frequency-dependent selection possible through associations between mating rate and polygamy on one hand and sperm depletion and sperm competition on the other [1, 5, 7, 16, 17]. In this scenario increased SR distorter frequencies result in an increasingly biased population sex ratio and in a decrease of the average semen quality, which can cause an increase in the female or male mating rate. An increased mating rate benefits wild-type alleles at the expense of SR distorters as their carriers perform better in sperm competition (with increased female mating rate) and show less decrease in fertility after repeated matings due to sperm depletion (with increased male mating rate). Also SR suppressors are hypothesised to have a fitness cost considering that they are polymorphic in Drosophila spp. . A fitness cost is necessary to prevent gene fixation in view of their selective advantage against wild-type alleles. (Wild-type alleles have a low transmission rate when co-occurring with a SR distorter.) Alternatively, stable gene polymorphism at the Y chromosome may result from interactions in the population dynamics of SR distorters and SR suppressors as indicated in simulation studies, albeit in a small fraction of parameter space [13, 19–21].
The African buffalo (Syncerus caffer) population of the Kruger National Park (KNP, South Africa, 22°-25° S, 31°-32° E, 20,000 km2) consists of some 30,000 individuals [22, 23]. The population was completely fenced from neighbouring populations in Zimbabwe and Mozambique between 1976 and 2004. The KNP buffalo are polygamous and live in both small bachelor herds and mixed-sex herds of a few dozen to over 1,000 individuals [24, 25]. Reproduction is strongly related to the availability and seasonality of resources with around 75% of all births occurring late in the wet season, from January to March . The foetal sex ratio for the population as a whole does not differ significantly from equality across decades (1969-1998: Pexact = 0.52, sex ratio (male proportion) = 0.506, 95% CI = 0.487-0.526, n = 2626) in agreement with "Fisher's Principle", which predicts the evolution of an equal birth sex ratio, assuming equal cost of daughters and sons [27, 28].
High Y-chromosomal microsatellite polymorphism, characterized by rainfall-correlated haplotype-frequency shifts between year cohorts, has been observed not only in the KNP population but also in the Hluhluwe-iMfolozi Park (HiP) population located some 300 km farther south (gene diversity H, probability of randomly sampling two different alleles or haplotypes; H = 0.48 and 0.74 in HiP and KNP respectively; haplotype number = 5 and 15 in HiP and KNP, respectively) . The gene diversity in the KNP population is similar to that recently observed in Eurasian cattle (all breeds combined: H = 0.82; mean of the individual breeds: H = 0.63), although the estimate for cattle is positively biased due to a larger number of microsatellites (five vs. three in African buffalo) . When looking at the individual microsatellites, gene diversity and allelic diversity are higher in the single African buffalo population of KNP than in all Eurasian cattle breeds combined (African buffalo: H = 0.60, 0.66, and 0.67 with 5, 7 and 9 alleles for microsatellites UM1113, UM0304 and INRA189, respectively; cattle: H = 0.05, 0.04, 0.36, 0.47 and 0.52 with 2, 3, 4, 4, and 8 alleles for microsatellites INRA124, BM861, BYM1, DYZ1 and INRA189, respectively; supplementary info in ).
The observation of high genetic diversity despite yearly haplotype frequency fluctuations was surprising as the latter is expected to quickly result in gene fixation . The presence of multiple haplotypes cannot be attributed to population subdivision. There is no significant genetic differentiation between Gonarezhou NP (Zimbabwe), northern KNP and southern KNP indicating extensive (historical) gene flow (autosomal microsatellites: 95% CI FST = -0.017 - 0.031) . This is supported by radio-tracking data of over 130 individual KNP buffalo that show all sex and age groups move between mixed-sex herds, with monthly emigration rates of 9-26% for adult males and 2% for females and subadults . The high haplotype diversity also cannot be explained by recent mutations as three of the five major haplotypes in KNP (112, 557 and 436; total frequency = 74%) differ by at least seven mutational events (Figure 1 in ). The presence of multiple haplotypes therefore indicates some kind of balancing selection. Balancing selection is expected to result in an increased gene diversity relative to haplotype (allele) diversity because of a relatively even haplotype (allele) frequency distribution. This is supported by the observation that Y-chromosomal gene diversity is higher in African buffalo than in Eurasian cattle when comparing microsatellites with a similar number of alleles (H = 0.66-0.67 with 5-7 alleles in the KNP buffalo vs. H = 0.52 with 8 alleles in cattle; H = 0.46 with 2 and 3 alleles for microsatellites UM1113 and UM0304, respectively, in the HiP buffalo vs. H = 0.04-0.05 with 2 alleles in cattle).
We hypothesise that the frequency shifts in Y-chromosomal haplotypes between year cohorts in the African buffalo population of KNP are indicative of SR genes. Accordingly, we expected to find temporal correlations between Y-chromosomal haplotype frequencies and sex ratios of foetuses and calves (even though the foetal sex ratio for the population as a whole does not differ significantly from equality).
A temporal genetic analysis of genetic differentiation in calves using Y-chromosomal microsatellite data demonstrated highly significant (Prandomisation-based = 0.00089, n = 74) differences between male calves conceived during a dry season (age cohorts from 1998: 0-6 months, 1-1.5 years) and those conceived during a wet season (age cohorts from 1998: 6-12 months, 1.5-2 years). This was mainly due to an opposite (negative) correlation between haplotypes 112 and 557 (Figure 1). Haplotype 557 was observed only among males conceived during the observed dry seasons, whilst haplotype 112 was observed only among males conceived during the observed wet seasons. A negative (112 and 557 in different cohorts) or positive (112 and 557 in the same cohort) correlation of this magnitude by chance for these two haplotypes, either between seasons (wet vs. dry) or years (0-year vs. 1-year old calves), is highly unlikely (Pexact = 0.00022, sum of all four possibilities, n = 74).
If haplotypes 112 and 557 are linked to SR distorters or SR suppressors then one would expect an alternating seasonal pattern with sex ratios (male proportion) as well, which was indeed observed (Figure 2). In the same age cohorts as were used in the genetic analyses the sex ratio was higher among individuals conceived during the wet seasons than among those conceived during the dry seasons (Pexact = 0.096, n = 291; calves culled in 1998 plus foetuses collected in 1996). There was an identical pattern in earlier age cohorts. Among 0-4 months old calves culled between June 1982 and February 1983, the sex ratio was only 0.27 (8/30) among those conceived in the dry season of 1981 (culled in June-November) but as high as 0.69 (9/13) among those conceived in the subsequent wet season (culled in January-February) (Pexact = 0.016) . Most importantly, a significant seasonal difference in sex ratio, with an effect size of 0.045 (absolute difference between two sex ratios), was observed among the foetuses collected between 1978 and 1998 (Figure 3 and Table 1; Pexact = 0.025; data pooled across years). The sex ratio was female-biased during the dry season (PWilsonscore = 0.049, one-sided) and male-biased during the wet season (PWilsonscore = 0.056, one-sided). The largest difference was observed at the change of seasons between March and April. The seasonal difference in sex ratio was already present among foetuses ≤ 4 months old (Table 1; Pexact = 0.035; data pooled across years), ruling out late-term foetal loss as an explanation.
It was previously noted that the frequency of haplotype 557 per year cohort is correlated with the mean annual rainfall in the three years before birth, showing high frequencies after dry years (Plogisticregression < 0.0001, years running from Sept. till Aug., n = 201) . Additional analyses show that the same is true for haplotype 112, showing high frequencies after wet years (Plogisticregression = 0.044, years running from Sept. till Aug., n = 201). Thus haplotypes 557 and 112 show an opposite correlation not only across seasons, but also across years (Figure 4). Again, this pattern is associated with sex-ratio variation. There was a positive correlation between the annual foetal sex ratio (years running from Nov. till Oct.) and the mean annual rainfall in the previous two or three years, including or excluding part of the wet season of the focal year (previous two years: Pspearman = 0.0065-0.059, n = 14; for each of the four 2-year periods starting at Nov., Dec., Jan. or Feb.; previous three years: Pspearman = 0.037-0.088, n = 14; for each of the six 3-year periods starting at Nov., Dec., Jan., Feb. Mar. or Apr.). With a rainfall period of three years including the complete wet season of the focal year (years running from Apr. till Mar.), the effect size of the sex ratio variation between dry and wet periods was 0.047 (Figure 4; Pspearman = 0.041, n = 14), which is similar to that between seasons (rainfall > 520 mm: sex ratio = 0.531, 95% CI = 0.503-0.560, n = 1178; rainfall ≤ 520 mm: sex ratio = 0.484, 95% CI = 0.461-0.508, n = 1767). The foetal sex ratio was significantly female-biased in 1990 (PWilsonscore = 0.0005, n = 332, one-sided) but significantly male-biased in 1998 (PWilsonscore = 0.046, n = 114, one-sided). The average foetal sex ratio across years did not deviate from equality (Pexact = 0.65, sex ratio = 0.504, 95% CI = 0.486-0.522, n = 2953; data pooled across years) as has also been observed in an earlier study .
The frequency of haplotype 557 was ≤ 0.08 (Pexact = 0.90, n = 28, one-sided) among males conceived during the wet seasons of 1996 and 1997, but it was ≥ 0.40 (PWilsonscore = 0.90, n = 38, one-sided) among those conceived between 1992 and 1994 following the severe drought of 1992 (rainfall: 243 mm vs. long-term average of 542 mm in the period 1961-2004, years running from Sept. till Aug.). This constitutes a difference of at least a factor of five in half a generation (assuming a generation time of 7.5 years ).
The conception peak occurred between February and April, corresponding to 67.6% (95% CI = 65.9-69.1%, n = 3298) of all conceptions. It differs by one month compared to the birth peak as the gestation period is around 11 months. No more than 33.7% (95% CI = 32.1-35.3%, n = 3298) of the conceptions occurred during the dry season (Apr.-Oct.). There was a negative correlation between the annual pregnancy rate among ≥ 2-year old females (at the moment of sampling) and the annual rainfall (Pspearman ≤ 0.021, n = 15; for each of five 1-year periods spanning previous and focal year and starting at May, Jun., Jul, Aug. or Sept.; Figure 5 with 1-year period Sept.-Aug., Pspearman = 0.011). The mean pregnancy rate was 51.6% (95% CI = 50.4-52.9%, n = 6439), with the yearly rate varying between 19% (1992, the driest year in the period 1961-2004) and 68% (1986, 1995, 1996). The mean adult pregnancy rate (≥ 6-year old) was as high as 69.8% (95% CI = 68.3-71.2%, n = 3821).
A plausible explanation for the associations between Y-chromosomal haplotype frequencies and foetal sex ratio is that haplotypes 112 and 557 are linked to SR genes. Haplotype 112 would be linked to a Y-chromosomal SR distorter, resulting in high population sex ratios during wet periods, and haplotype 557 to a Y-chromosomal SR suppressor, causing this haplotype to increase its frequency during dry periods when population sex ratios are low (schematic overview of the proposed mechanism in Table 2). The low sex ratios cannot be explained by mere inactivation of the SR distorter considering that they were significantly female-biased in 1990 and during the dry seasons. The SR-suppressor carrying males produce equal numbers of sons and daughters during dry periods (in dry seasons and after dry years), considering the high frequency of haplotype 557 while the sex ratio is low, but very few offspring during wet periods (in wet seasons and after wet years), as the then high sex ratio precludes a low frequency due to female-biased conception. On the other hand, SR-distorter carrying males produce substantially more sons than daughters during wet periods considering the then high frequency of haplotype 112. In turn, we argue that they produce very few offspring during dry periods considering the low frequency of haplotype 112 at those times. Although the frequency of haplotype 112 is low (0.11 after wet years and 0.18 during wet seasons), it is high enough to explain sex ratios as high as 0.53-0.55 during wet periods, assuming a frequency of 0.08-0.10 among fathers and 95% male progeny (a frequency of 0.08-0.10 among fathers results in a frequency of 0.11-0.18 among sons).
In the scenario described above, the fitness cost of temporarily reduced fertility prevents the SR distorter and SR suppressor from going to fixation. Our failure to observe haplotype 112 and 557 in dry and wet season cohorts, respectively, may even indicate a temporal reduction in fertility that nears sterility. A strong temporal reduction in fertility co-occurring with SR distortion is not implausible considering that reductions of more than 50% and even sterility have been observed with SR distorters in Drosophila spp., with the latter being caused by a single gene [5, 11].
Female-biased sex ratios and the presence of a Y-chromosomal SR suppressor during dry periods indicate that an X-chromosomal SR distorter may also be present in KNP buffalo (Table 2). The two SR distorters and the SR suppressor may have a common origin, as it has recently been shown that a SR suppressor (Nmy) in Drosophila simulans has originated as a retro-transposed duplication of a SR distorter (Dox) [7, 8, 12]. However, we cannot exclude the possibility that the female-biased sex ratios are not due to an X-chromosomal SR distorter but instead are the result of an inherent property of Y-bearing spermatozoa making them relatively vulnerable to stress and aging [35–38]. The Y-chromosomal SR suppressor in such a case would act by making the X-bearing spermatozoa equally vulnerable to stress and aging.
The SR distorter and SR suppressor are likely to exert their effect during spermatogenesis, which is supported by two observations. Firstly, many Y-chromosomal genes are known to be involved in spermatogenesis [39–43]. They may act by influencing the transcription of various genes through epigenetic and epistatic mechanisms [7, 39, 44]. It is also possible that the SR distorter and SR suppressor exert their main effect post-copulation in the female genital tract by influencing sperm transport, sperm survival or fertilization capability. Secondly, in our study population it was shown that ejaculate volume, sperm motility and proportion of morphologically normal spermatozoa decrease significantly during the dry season (April-May vs. October-November) . Studies on cattle (Bos Taurus), red deer (Cervus elaphus), guinea pig (Cavia porcellus), mice (Mus musculus) and humans have shown correlations between semen quality parameters on one hand and spermatozoal sex ratio, foetal sex ratio, offspring sex ratio and sex-differential fertilizing capability of spermatozoa on the other [37, 38, 46–50].
In the aforementioned mammals a low semen quality is associated with a low sex ratio, which also is the case for the African buffalo (i.e. low foetal sex ratio in dry season when semen quality is low). These observations may point towards a general mechanism in mammals whereby semen-quality related sex-ratio variation is driven by SR genes (distorters and suppressors). A possible explanation is that a positive association evolves more readily as males with high-quality semen may benefit from producing sons, who can inherit their father's semen quality, while males with low-quality semen benefit from producing daughters . This would increase the transmission rate of newly mutated SR genes, which can be important for their invasion success as stochastic events increase the risk of establishment failure .
In African buffalo the association between SR genes and semen quality may have set in motion the co-evolution of Y-chromosomal and X-chromosomal SR distorters and their suppressors whose temporal variation of activity (distortion/suppression in one season, low fertility in the other) is indirectly triggered by the amount of rainfall. Here we have shown that Y-chromosomal haplotype frequencies, foetal sex ratios and pregnancy rates (indicator of female body condition) are all associated with rainfall. The decrease in semen quality in the African buffalo during the dry season is likely due to decreasing availability and quality of food resources, which are directly related to the amount of rainfall . Rainfall is one of the most critical environmental factors affecting resource availability in savannah ecosystems [26, 53]. In other African bovids sperm production is reduced by 30-90% during the non-breeding season .
The characteristics of the SR genes in African buffalo are different from those described in other species to date. The SR distorter and SR suppressor occur on a single sex chromosome and show a temporal and opposite variation in gene activity. Furthermore, there are indications that both an X-chromosomal and Y-chromosomal SR distorter are present in the population. Another unique feature is that female buffaloes in KNP are also able to influence the foetal sex in relation to body condition (in males semen quality, which likely correlates with male body condition). The effect is opposite to that of males, with females in above average body condition (although population average may be low) investing in daughters and those in below average body condition investing in sons .
Constant fitness models show that the co-occurrence of an X- and Y-chromosomal SR distorter can result in a stable gene polymorphism at the Y chromosome . Although this only seems possible in 2% of the parameter space, extra stability may be provided by additional factors not included in the models. Such factors may be presence of a Y-chromosomal suppressor, opposite effects of the two parents on foetal sex ratio and frequency-dependent selection through an association between mating rate and sperm quality (depletion and competition). Furthermore, the constant fitness models show cycling of allele (haplotype) frequencies, which in our study population is observed in relation to rainfall.
We consider it unlikely that the postulated Y-chromosomal SR-correlated genes are associated with male mating success or semen quality (fertilization success) as such since these cannot explain sex-ratio variation. That the postulated genes exert their main effect post-conception by influencing pregnancy loss (embryonic and foetal mortality) seems equally unlikely. To prevent gene fixation, pregnancy losses with wild-type haplotypes should be similar to those with haplotypes 112 and 557, which would have to be considerable (at least 50%) to explain the large haplotype frequency variations (up to a factor of five). Moreover, additional losses of female embryos and foetuses must be assumed to explain the male-biased sex ratios. However, high pregnancy losses are difficult to reconcile with adult pregnancy rates of 70%. Furthermore, this scenario would imply that haplotype 557 is associated with below average pregnancy loss during dry periods while the opposite is true for the remaining haplotypes, as indicated by the then low pregnancy rates and conception numbers for the population as a whole. It is difficult to envisage a physiological mechanism capable of making haplotype 557 so successful under stressful conditions (dry periods), yet so unsuccessful under beneficial conditions (wet periods). Cryptic female choice at least seems unlikely since this would imply a positive adaptive value of choosing male offspring on the basis of just a single Y-chromosomal haplotype, unlinked to the autosomal genome. Furthermore, haplotype-specific pregnancy loss is difficult to relate to frequency-dependent selection. This would require that the physiological mechanism by which the haplotypes are associated with pregnancy loss in individual females is modulated, directly or indirectly, by their frequency in the population at large. A final alternative explanation would be a co-occurrence of two independent mechanisms: rainfall-dependent and haplotype-specific mating or fertilization success in the fathers on one hand and rainfall-dependent and sex-specific pregnancy loss in the mothers on the other. Again, this scenario is difficult to relate to frequency-dependent selection while it also requires a mechanism that causes haplotype 557 to be most successful under stressful conditions but least successful under beneficial conditions.
Although we consider the alternative explanations for the correlations between Y-chromosomal haplotype frequencies and foetal sex ratio unlikely, and propose the presence of SR genes as the most plausible explanation, this will need direct confirmation in a future studies. This may be done by analyzing associations between semen quality and spermatozoal sex ratio in individual bulls with known Y-chromosomal genotype.
Polymorphic active SR genes in natural populations of mammals have rarely been observed but may be more common than previously thought. We suggest that in many mammals semen-quality related primary sex-ratio variation is influenced by SR genes. Stable gene polymorphism may be maintained with relative ease when SR genes fulfil one or more of the specific characteristics observed in African buffalo.
From 1969 onwards, culling of African buffalo has been used in KNP as a population management and disease monitoring tool [27, 54]. The data are from 3184 foetuses and their mothers collected during culling events in 1978 and between 1984 and 1996 (data from SANParks Scientific Services Department), mostly between May and November (99.4%), from 218 0-1 year-old calves and from 114 foetuses and their mothers culled from 32 herds between September and November 1998, and from published data from 43 0-4 month-old calves culled between June 1982 and February 1983 [27, 33, 54].
The foetuses from 1998 are the same as those reported previously . Foetus sex was determined using external genitalia characteristics. Sexed foetuses ≤ 18 g (age ± 70 days, n = 150) were excluded from the data analyses because they appeared to be outliers, considering that they were predominantly diagnosed as male, irrespective of season or year (76% males, 114/150; 95% CI: 69-82%). This is most likely to be attributed to sexing errors and ascertainment bias. In the closely related Asian buffalo (Bubalus bubalis) reliable foetal sex differentiation is possible no sooner than the 56th day of pregnancy with macroscopic diagnosis and the 70th day of pregnancy with ultrasonographic diagnosis, with external genitalia occurring sooner in males than in females [55, 56]. Furthermore, 195 foetuses from 1978-1996 > 18 g were of unknown sex, leaving 2839 foetuses from 1978-1996 and 114 foetuses from 1998 for sex-ratio analyses. The sex ratio of these foetuses did not significantly deviate from equality (see Results section) indicating that there was no systematic bias in foetus sexing. Conception dates were assigned using a weight-at-age regression (8.57*weight (g)1/3 + 48.5) .
The age of the animals was estimated in years on the basis of dental wear patterns, number of erupted incisor teeth, body size, and horn development . The age of the calves from 1998 was estimated in 3-month and half-year cohorts (age categories: 0-3 months, n = 5; 3-6 months, n = 63; 6-9 months, n = 43; 9-12 months, n = 5; 1-1.5 year, n = 73; 1.5-2 year, n = 29). Conception dates of the calves were estimated on the basis of the midpoint age per age cohort (i.e. 1.5, 4.5, 7.5, 10.5, 15 and 21 months for the 1998 data and 2 months for the 1982-1983 data), assuming a gestation period of 340-343 days [25, 58].
Y-chromosomal microsatellite marker data were available from the males culled in 1998, consisting of 74 0-1 year-old calves and 127 subadults and adults (≥ 2 years old). In an earlier study, which gives detailed information on the molecular analyses, these data were used to assess the genetic diversity per year cohort . Three microsatellites were analyzed whose combination of alleles constitute a haplotype, since the Y chromosome is a single, haploid non-recombining unit (except for a small pseudoautosomal region). A 3-number code was assigned to each haplotype, with each number corresponding to a microsatellite allele. Fifteen haplotypes were observed in total. No foetus samples were available for DNA analyses since they were not stored after analysis in the field. Here we used the microsatellite data to test for genetic differentiation among season cohorts in calves (Figure 1). Randomisation-based tests for genetic differentiation with 100,000 repetitions and no pooling of haplotypes were performed using FSTAT 18.104.22.168 [59, 60]. The program generates contingency tables of haplotypes by samples (age cohorts) and the statistic that is used to classify them is the log-likelihood statistic G.
Effect sizes were estimated as the absolute difference between sex ratios, which is defined as the proportion of males. Confidence intervals of the sex ratios were estimated using the Wilson score formula . The significance level of differences between two sex ratios or haplotype frequencies was estimated with the Fisher's exact test. Haplotype frequencies were regressed against rainfall in forward and backward stepwise (conditional) logistic regression models using SPSS 22.214.171.124. Correlations were analysed using the Spearman rank test. All significance tests in this study were performed two-sided, unless otherwise indicated.
Monthly rainfall data for September 1960 to December 2004 were averaged across 13 South African Weather Service (SAWS) rainfall stations in KNP. The body condition of adult female buffalo in KNP coincides with monthly Normalized Difference Vegetation Index (NDVI) values, which lags the rainfall by one month and which is a proxy for vegetation productivity . NDVI values tend to increase sharply in November and reach their peak in February-March . The wet season was therefore defined as the period from November to March. During this period 78% of the rainfall occurs, which is five times higher on a monthly basis than in the dry season.
The raw data used in this study are available in an additional file (Additional file 1).
Jaenike J: Sex chromosome meiotic drive. Annu Rev Ecol Syst. 2001, 32: 25-49. 10.1146/annurev.ecolsys.32.081501.113958.
Joseph SB, Kirkpatrick M: Haploid selection in animals. Trends Ecol Evol. 2004, 19 (11): 592-597. 10.1016/j.tree.2004.08.004.
Lyon MF: Transmission ratio distortion in mice. Annu Rev Genet. 2003, 37: 393-408. 10.1146/annurev.genet.37.110801.143030.
Kleene KC: Sexual selection, genetic conflict, selfish genes, and the atypical patterns of gene expression in spermatogenic cells. Dev Biol. 2005, 277 (1): 16-26. 10.1016/j.ydbio.2004.09.031.
Price TAR, Wedell N: Selfish genetic elements and sexual selection: their impact on male fertility. Genetica. 2008, 132 (3): 295-307. 10.1007/s10709-007-9173-2.
Jaenike J: Suppression of sex-ratio meiotic drive and the maintenance of Y-chromosome polymorphism in Drosophila . Evolution. 1999, 53 (1): 164-174. 10.2307/2640929.
Tao Y, Masly JP, Araripe L, Ke Y, Hartl DL: A sex-ratio meiotic drive system in Drosophila simulans. I: an autosomal suppressor. PLoS Biol. 2007, 5 (11): 2560-2575.
Ferree PM, Barbash DA: Distorted sex ratios: a window into RNAi-mediated silencing. PLoS Biol. 2007, 5 (11): 2453-2457. 10.1371/journal.pbio.0050303.
Merrill C, Bayraktaroglu L, Kusano A, Ganetzky B: Truncated RanGAP encoded by the segregation distorter locus of Drosophila . Science. 1999, 283 (5408): 1742-1745. 10.1126/science.283.5408.1742.
Montchamp-Moreau C: Sex-ratio meiotic drive in Drosophila simulans : cellular mechanism, candidate genes and evolution. Biochem Soc Trans. 2006, 34: 562-565. 10.1042/BST0340562.
Phadnis N, Orr HA: A single gene causes both male sterility and segregation distortion in Drosophila hybrids. Science. 2009, 323 (5912): 376-379. 10.1126/science.1163934.
Tao Y, Araripe L, Kingan SB, Ke Y, Xiao H, Hartl DL: A sex-ratio meiotic drive system in Drosophila simulans. II: an X-linked distorter. PLoS Biol. 2007, 5 (11): 2576-2588.
Clark AG: Natural selection and Y-linked polymorphism. Genetics. 1987, 115 (3): 569-577.
Clark AG: Two tests of Y chromosomal variation in male fertility of Drosophila melanogaster . Genetics. 1990, 125 (3): 527-534.
Dean AM: Protecting haploid polymorphisms in temporally variable environments. Genetics. 2005, 169 (2): 1147-1156. 10.1534/genetics.104.036053.
Bernasconi G, Ashman TL, Birkhead TR, Bishop JDD, Grossniklaus U, Kubli E, Marshall DL, Schmid B, Skogsmyr I, Snook RR, Taylor D, Till-Bottraud I, Ward PI, Zeh DW, Hellriegel B: Evolutionary ecology of the prezygotic stage. Science. 2004, 303 (5660): 971-975. 10.1126/science.1092180.
Price TAR, Bretman AJ, Avent TD, Snook RR, Hurst GDD, Wedell N: Sex ratio distorter reduces sperm competitive ability in an insect. Evolution. 2008, 62 (7): 1644-1652. 10.1111/j.1558-5646.2008.00386.x.
Montchamp-Moreau C, Ginhoux V, Atlan A: The Y chromosomes of Drosophila simulans are highly polymorphic for their ability to suppress sex-ratio drive. Evolution. 2001, 55 (4): 728-737. 10.1554/0014-3820(2001)055[0728:TYCODS]2.0.CO;2.
Carvalho AB, Vaz SC, Klaczko LB: Polymorphism for Y-linked suppressors of sex-ratio in two natural populations of Drosophila mediopunctata . Genetics. 1997, 146 (3): 891-902.
Taylor JE, Jaenike J: Sperm competition and the dynamics of X chromosome drive: Stability and extinction. Genetics. 2002, 160 (4): 1721-1731.
Hall DW: Meiotic drive and sex chromosome cycling. Evolution. 2004, 58 (5): 925-931.
Rodwell TC, Whyte IJ, Boyce WM: Evaluation of population effects of bovine tuberculosis in free-ranging African buffalo (Syncerus caffer). J Mammal. 2001, 82 (1): 231-238. 10.1644/1545-1542(2001)082<0231:EOPEOB>2.0.CO;2.
Whyte I: Kruger NP buffalo census. Gnusletter. 2005, 24 (2): 14-16.
Cross PC, Lloyd-Smith JO, Getz WM: Disentangling association patterns in fission-fusion societies using African buffalo as an example. Anim Behav. 2005, 69: 499-506. 10.1016/j.anbehav.2004.08.006.
Prins HHT: Ecology and behaviour of the African buffalo: social inequality and decision making. 1996, London: Chapman and Hall
Ryan SJ, Knechtel CU, Getz WM: Ecological cues, gestation length, and birth timing in African buffalo (Syncerus caffer). Behav Ecol. 2007, 18 (4): 635-644. 10.1093/beheco/arm028.
Visscher DR, van Aarde RJ, Whyte I: Environmental and maternal correlates of foetal sex ratios in the African buffalo (Syncerus caffer) and savanna elephant (Loxodonta africana). J Zool. 2004, 264: 111-116. 10.1017/S0952836904005436.
Fisher RA: The genetical theory of natural selection. 1930, Oxford: Claredon Press
van Hooft P, Greyling BJ, Prins HHT, Getz WM, Jolles AE, Bastos ADS: Selection at the Y chromosome of the African buffalo driven by rainfall. PLoS ONE. 2007, 10.1371/journal.pone.0001086(10)
Kantanen J, Edwards CJ, Bradley DG, Viinalass H, Thessler S, Ivanova Z, Kiselyova T, Cinkulov M, Popov R, Stojanovic S, Ammosov I, Vilkki J: Maternal and paternal genealogy of Eurasian taurine cattle (Bos taurus). Heredity. 2009, 103 (5): 404-415. 10.1038/hdy.2009.68.
van Hooft WF, Groen AF, Prins HHT: Microsatellite analysis of genetic diversity in African buffalo (Syncerus caffer) populations throughout Africa. Mol Ecol. 2000, 9 (12): 2017-2025. 10.1046/j.1365-294X.2000.01101.x.
Cross PC, Getz WM: Assessing vaccination as a control strategy in an ongoing epidemic: Bovine tuberculosis in African buffalo. Ecol Model. 2006, 196 (3-4): 494-504. 10.1016/j.ecolmodel.2006.02.009.
Mason DR: Sex and age data from cropping of buffalo Syncerus caffer in the Kruger National Park. Koedoe. 1984, 27: 73-78.
O'Ryan C, Harley EH, Bruford MW, Beaumont M, Wayne RK, Cherry MI: Microsatellite analysis of genetic diversity in fragmented South African buffalo populations. Animal Conservation. 1998, 1 (2): 85-94. 10.1111/j.1469-1795.1998.tb00015.x.
James WH: Evidence that mammalian sex ratios at birth are partially controlled by parental hormone levels around the time of conception. J Endocrinol. 2008, 198 (1): 3-15. 10.1677/JOE-07-0446.
James WH: Some comments on the paper of Grant (2007). J Theor Biol. 2008, 253 (2): 401-404. 10.1016/j.jtbi.2008.03.014.
Hendricks KEM, Martins L, Hansen PJ: Consequences for the bovine embryo of being derived from a spermatozoon subjected to post-ejaculatory aging and heat shock: development to the blastocyst stage and sex ratio. J Reprod Dev. 2009, 55 (1): 69-74. 10.1262/jrd.20097.
Perez-Crespo M, Pintado B, Gutierrez-Adan A: Scrotal heat stress effects on sperm viability, sperm DNA integrity, and the offspring sex ratio in mice. Mol Reprod Dev. 2008, 75 (1): 40-47. 10.1002/mrd.20759.
Chippindale AK, Rice WR: Y chromosome polymorphism is a strong determinant of male fitness in Drosophila melanogaster . Proc Natl Acad Sci USA. 2001, 98 (10): 5677-5682. 10.1073/pnas.101456898.
Cummins J: Mitochondrial DNA and the Y chromosome: parallels and paradoxes. Reprod Fertil Dev. 2001, 13 (7-8): 533-542. 10.1071/RD01064.
Quintana-Murci L, Krausz C, McElreavey K: The human Y chromosome: function, evolution and disease. Forensic SciInt. 2001, 118 (2-3): 169-181. 10.1016/S0379-0738(01)00387-5.
Roldan ERS, Gomendio M: The Y chromosome as a battle ground for sexual selection. Trends Ecol Evol. 1999, 14 (2): 58-62. 10.1016/S0169-5347(98)01567-5.
Skaletsky H, Kuroda-Kawaguchi T, Minx PJ, Cordum HS, Hillier L, Brown LG, Repping S, Pyntikova T, Ali J, Bieri T, Chinwalla A, Delehaunty A, Delehaunty K, Du H, Fewell G, Fulton L, Fulton R, Graves T, Hou SF, Latrielle P, Leonard S, Mardis E, Maupin R, McPherson J, Miner T, Nash W, Nguyen C, Ozersky P, Pepin K, Rock S, Rohlfing T, Scott K, Schultz B, Strong C, Tin-Wollam A, Yang SP, Waterston RH, Wilson RK, Rozen S, Page DC: The male-specific region of the human Y chromosome is a mosaic of discrete sequence classes. Nature. 2003, 423 (6942): 825-U822. 10.1038/nature01722.
Lemos B, Araripe LO, Hartl DL: Polymorphic Y chromosomes harbor cryptic variation with manifold functional consequences. Science. 2008, 319 (5859): 91-93. 10.1126/science.1148861.
Brown JL, Wildt DE, Raath JR, Devos V, Howard JG, Janssen DL, Citino SB, Bush M: Impact of season on seminal characteristics and endocrine status of adult free-ranging African buffalo (Syncerus caffer). J Reprod Fertil. 1991, 92 (1): 47-57.
Chandler JE, Canal AM, Paul JB, Moser EB: Collection frequency affects percent Y-chromosome bearing sperm, sperm head area and quality of bovine ejaculates. Theriogenology. 2002, 57 (4): 1327-1346. 10.1016/S0093-691X(01)00721-X.
Gomendio M, Malo AF, Soler AJ, Fernandez-Santos MR, Esteso MC, Garcia AJ, Roldan ERS, Garde J: Male fertility and sex ratio at birth in red deer. Science. 2006, 314 (5804): 1445-1447. 10.1126/science.1133064.
Hwang SY, Kim WJ, Wee JJ, Choi JS, Kim SK: Panax ginseng improves survival and sperm quality in guinea pigs exposed to 2,3,7,8-tetrachlorodibenzo-p-dioxin. BJU International. 2004, 94 (4): 663-668. 10.1111/j.1464-410X.2004.05019.x.
Jimenez A, Fernandez R, Madrid-Bury N, Moreira PN, Borque C, Pintado B, Gutierrez-Adan A: Experimental demonstration that pre- and post-conceptional mechanisms influence sex ratio in mouse embryos. Mol Reprod Dev. 2003, 66 (2): 162-165. 10.1002/mrd.10345.
Zorn B, Sucur V, Stare J, Meden-Vrtovec H: Decline in sex ratio at birth after 10-day war in Slovenia. Hum Reprod. 2002, 17 (12): 3173-3177. 10.1093/humrep/17.12.3173.
Gomendio M, Malo AF, Garde J, Roldan ERS: Sperm traits and male fertility in natural populations. Reproduction. 2007, 134 (1): 19-29. 10.1530/REP-07-0143.
Brown BW: A review of nutritional influences on reproduction in boars, bulls and rams. Reprod Nutr Dev. 1994, 34 (2): 89-114. 10.1051/rnd:19940201.
Owen-Smith N, Mason DR, Ogutu JO: Correlates of survival rates for 10 African ungulate populations: density, rainfall and predation. J Anim Ecol. 2005, 74 (4): 774-788. 10.1111/j.1365-2656.2005.00974.x.
Rodwell TC, Kriek NP, Bengis RG, Whyte IJ, Viljoen PC, de Vos V, Boyce WM: Prevalence of bovine tuberculosis in African buffalo at Kruger National Park. J Wildl Dis. 2001, 37 (2): 258-264.
Ranjbar R, Rashidi SH, Alboghobeish N, Sadrkhanloo RA, Mazaheri Y: Study of fetal sex determination based on external genitalia and gonadal differentiation in water buffalos. Internet J Vet Med. 2007, 3:
Ali A, Fahmy S: Ultrasonographic fetometry and determination of fetal sex in buffaloes (Bubalus bubalis). Anim Reprod Sci. 2008, 106 (1-2): 90-99. 10.1016/j.anireprosci.2007.04.010.
Reeves JT, Gentry M, Daoud FS: Growth of fetal calf and its arterial pressure, blood gases and hematologic data. J Appl Physiol. 1972, 32 (2): 240-244.
Schmidt S, Gerber D, Soley JT, Aire TA, Boos A: Histo-morphology of the uterus and early placenta of the African buffalo (Syncerus caffer) and comparative placentome morphology of the African buffalo and cattle (Bos taurus). Placenta. 2006, 27 (8): 899-911. 10.1016/j.placenta.2005.09.008.
Goudet J: FSTAT, a program to estimate and test gene diversities and fixation indices (version 2.9.3). Updated from Goudet (1995). 2001
Goudet J: FSTAT (Version 1.2): A computer program to calculate F-statistics. J Hered. 1995, 86 (6): 485-486.
Brown LD, Cai TT, DasGupta A, Agresti A, Coull BA, Casella G, Corcoran C, Mehta C, Ghosh M, Santner TJ, Brown LD, Cai TT, DasGupta A: Interval estimation for a binomial proportion - Comment - Rejoinder. Stat Sci. 2001, 16 (2): 101-133.
We thank Dr. Roy Bengis, Dr. Dave Cooper, Dr. Ian Whyte and SANParks Scientific Services Department for the field data and assistance with sample collection, and Prof. Tim Clutton-Brock, Dr. Arjan de Visser and Dr. Paul Funston for providing helpful comments. The molecular DNA analyses were supported by US NIH/NSF Ecology of Infectious Disease Grant GM83863 awarded to WMG and performed by PvH on a post-doctoral research fellowship held at the University of Pretoria.
PvH designed and performed the molecular DNA analyses, analyzed the data and wrote the paper. PDvH and ADSB contributed samples. ADSB contributed laboratory reagents and facilities. All authors discussed and commented on the paper.