The association between height at adulthood and lifetime reproductive success in humans has attracted long-time scientific interest. It has been suggested that women show a preference for taller than average men when seeking sexual partners (Courtiol et al., Reference Courtiol, Raymond, Godelle and Ferdy2010), and taller men also have a higher probability of getting married (Nettle, Reference Nettle2002a). Since tallness is associated with better health and socioeconomic position (Batty et al., Reference Batty, Shipley, Gunnell, Huxley, Kivimäki, Woodward and Smith2009), favoring tall men may offer long-term benefits for the offspring of women, thus increasing their evolutionary fitness. In females, tall stature may reflect better physical health and thus increased capability of having children, especially in the presence of environmental stress (Pollet & Nettle, Reference Pollet and Nettle2008). There is evidence that tall stature is associated with better reproductive health in women, such as reduced need for emergency cesarean section (Stulp et al., Reference Stulp, Verhulst, Pollet, Nettle and Buunk2012b). However, in addition to the tentative advantage of taller stature on reproductive success, it may also have disadvantages. Growth is an energy-consuming process, and thus reaching tall stature consumes energy and time that could otherwise be used for reproduction. Thus, allocating energy to growth versus early reproduction may reflect different adaptive strategies (Walker et al., Reference Walker, Gurven, Hill, Migliano, Chagnon, De Souza and Yamauchi2006). These strategies probably offer different costs and benefits in different environments and may contribute to the large variation in height both between and within human populations (Eveleth & Tanner, Reference Eveleth and Tanner2003).
Previous studies in Western countries have shown that tall stature is associated with later timing of menarche (McIntyre & Kacerosky, Reference McIntyre and Kacerosky2011; Onland-Moret et al., Reference Onland-Moret, Peeters, van Gils, Clavel-Chapelon, Key, Tjonneland and Riboli2005), and a previous Finnish study also reported that tall stature is associated with older age at first birth (AFB) in women (Helle, Reference Helle2008). This may be, however, a special feature of developed societies as in agrarian societies tall stature has been found to predict earlier rather than later menarche (McIntyre & Kacerosky, Reference McIntyre and Kacerosky2011). When studying lifetime reproductive success in males, tallness was associated with a higher number of children in a study of US cadets (Mueller & Mazur, Reference Mueller and Mazur2001) and in the Polish population (Pawlowski et al., Reference Pawlowski, Dunbar and Lipowicz2000), but no association was found in a study of US college students (Scott & Bajema, Reference Scott and Bajema1982), in Mexicans living in the United States (Goldstein & Kobyliansky, Reference Goldstein and Kobyliansky1984), in the United Kingdom (Nettle, Reference Nettle2002a), in China (Fielding et al., Reference Fielding, Schooling, Adab, Cheng, Lao, Jiang and Lam2008), or in Gambia (Sear, Reference Sear2006). In another cohort of US college students (Vetta, Reference Vetta1975) and high-school graduates (Stulp et al., Reference Stulp, Pollet, Verhults and Buunk2012a), average stature was associated with the highest number of children. When studying childlessness, the results are more consistent: short stature in men was associated with higher probability of childlessness in the United Kingdom (Nettle, Reference Nettle2002a) and Poland (Pawlowski et al., Reference Pawlowski, Dunbar and Lipowicz2000). In females, a positive association between stature and lifetime reproductive success was found in Gambia (Sear, Reference Sear2006) and Guatemala (Pollet & Nettle, Reference Pollet and Nettle2008), whereas in the United Kingdom the association was curvilinear and women who were shorter than average had higher reproductive success than the shortest and taller than average women (Nettle, Reference Nettle2002b); no association was found in US college students (Scott & Bajema, Reference Scott and Bajema1982), in Mexicans living in the United States (Goldstein & Kobyliansky, Reference Goldstein and Kobyliansky1984), in China (Fielding et al., Reference Fielding, Schooling, Adab, Cheng, Lao, Jiang and Lam2008), or in Finland (Helle, Reference Helle2008).
A shortcoming in these previous studies on the association between height and reproduction is, however, that they have not properly adjusted for socioeconomic or other environmental covariates which could contribute to this association. Childhood socioeconomic background is especially important to consider as growth takes place during childhood and adolescence, and childhood environment can affect adult life in multiple ways. As height is associated with childhood socioeconomic background (Silventoinen, Reference Silventoinen2003) and high social position is found to be associated with higher fertility in males and lower fertility in females (Kravdal & Rindfuss, Reference Kravdal and Rindfuss2008; Nettle & Pollet, Reference Nettle and Pollet2008), these factors rather than height itself may explain this association. Adjusting the results for measured socioeconomic factors offers only a partial solution, as no study can fully observe all relevant social factors. In this study, we analyze this question utilizing a twin study design. Using twin pairs raised in the same family who differ in the number of children they have as adults allows an analysis of the role of height on reproductive success when family background can be regarded as standardized. Furthermore, the twin study design offers the opportunity to quantify the role of genetic and environmental factors behind the variation of different traits and covariation of their mutual associations. We used a large and nationally representative Finnish twin data set that allows participants to be followed up to 51–59 years of age, and thus we have been able to analyze the lifetime number of children also in men.
Data and Methods
The data were derived from the older cohort of the Finnish Twin Cohort Study (Kaprio & Koskenvuo, Reference Kaprio and Koskenvuo2002). Questionnaires were mailed in 1975 and 1981 to all same-sex Finnish twin pairs born before 1958 and with both co-twins alive in 1975. The response rates to these questionnaires were 89% and 84%, respectively. Height, marital history, and education, classified into eight classes described elsewhere (Silventoinen et al., Reference Silventoinen, Kaprio and Lahelma2000), were sent identical questionnaires at both times. Those respondents still studying when reporting their level of education were assumed to have reached the next possible category given in the questionnaire. For height and education, we used the questionnaire responses in 1981, and if these were not available then the data from the questionnaire in 1975 were used. Marital history was classified as currently living or previously lived with a spouse (married, co-habiting, divorced, and widowed) in 1981 or in 1975. Self-reported height showed high correlations with measured height in a subsample of 100 men and 127 women in this same cohort (Pearson r = .98 in men and .96 in women), showing good reliability of self-reported height (Silventoinen et al., Reference Silventoinen, Kaprio, Lahelma and Koskenvuo2000).
Information on live births is available for participants born after 1950 in the Finnish population register. This register was updated until June 2009 and linked to the cohorts born between 1950 and 1957 in our study sample, using unique personal identification numbers assigned to all Finnish citizens. The participants alive in 2009 were 51–59 years of age at the time of the end of the follow-up and therefore effectively post-reproductive. From 9,006 twin individuals who answered the questionnaires in 1975 or 1981, we removed 1,168 twin individuals who were living with their co-twin in 1981. This is because a large proportion of these twins remained childless and therefore this group may not be representative to the general population in terms of reproduction. However, we replicated the analyses in the whole data set to confirm whether this restriction affected the results. Furthermore, individuals from triplet or quadruplet births were excluded from the data (N = 8). In the final study sample, we had 3,595 male and 4,235 female twin individuals, including 1,076 monozygotic (MZ) and 2,445 dizygotic (DZ) complete twin pairs, and 787 twin individuals without information on their co-twin. The study was approved by the ethical committee of the Department of Public Health, University of Helsinki. Informed consent was obtained from all respondents.
We analyzed the association between height and lifetime number of children among twin individuals using both Poisson regression and multinomial logistic regression models producing odds ratios (OR). Poisson regression was used to make the results comparable with previous studies. However, we also conducted multinominal logistic regression models because the association between height and number of children may not be similar at different levels of fertility. For example, it is possible that different factors may affect the transition to first, second, third, and further parities; for example, mate choice is probably most crucial in having any children. For this reason, height was classified into sex-specific quintiles and the number of children into three categories (no children, 1–2 children, and 3 children or more) in the multinominal logistic regression models. We then adjusted the results for education used as a categorized covariate because of possible non-linearity and after that additionally for marital history. As the twins were sampled as pairs, there is a lack of statistical independence between twins in a pair. Therefore, within-pair correlations were taken into account in the estimation of the standard errors using a survey cluster option (Williams, Reference Williams2000). The analyses were continued by studying pairwise associations in twin pairs discordant for height, that is, belonging to a different quintile of height, and having any children using a conditional logistic regression model. Since twins reared together share their childhood family environment, and MZ twins also share their gene sequence, any within-pair associations would indicate that the association between height and reproduction is not because of familial factors. After these analyses, we conducted similar analyses for marital history to get more insight into the role of this factor in the formation of the association between height and reproduction. These analyses were carried out using the Stata/SE statistical software, version 11.2 (StataCorp, College Station, Texas, USA).
We continued the analysis by studying the association of height with lifetime number of children and AFB by applying genetic twin models if a linear association in regression models was found between them. This was done in order to disentangle the effects of genetic and environmental factors. The classical twin analysis is based on the comparison of similarity of MZ twins, who are genetically identical at gene sequence level, to DZ twins who share, on average, 50% of their segregating genes. The variation within traits and covariation between two traits were decomposed to additive genetic factors (A), which is the sum of the effects of all alleles affecting the trait; common environmental factors (C), including the effect of environment shared by co-twins; and specific environmental factors (E), including environmental effects unique to each twin individual as well as all measurement errors. This was done by using bivariate Cholesky decomposition. The assumptions of twin modeling, that is, equal means and variances for MZ and DZ twins, were tested by comparing the genetic models with saturated models, which do not make these assumptions. The comparisons in model fit between the nested models were based on the change of χ2-goodness-of-fit statistics and degrees of freedom (Δχ2degrees of freedom). The genetic modeling was carried out with the Mx statistical package, version 1.7.03, using the raw data analysis option (Neale, Reference Neale2003). As the distribution of AFB was somewhat skewed, we normalized it by using logarithmic transformation and used this in the statistical modeling.
Results
Table 1 presents descriptive statistics for height, AFB, and number of children. In childlessness and other categories of number of children, no linearity was seen over the height quintiles in men or women. An increasing trend in AFB over height quintiles was seen in women but not in men. In men, the regression coefficient of height was not statistically significant (β for 10 cm = 0.004; 95% CI = -0.008, 0.016) when predicting log AFB. In women, a clear negative association was found (β for 10 cm = -0.699; 95% CI = -1.229, -0.169). However, there was also some leveling off in this association in the tallest women, and a weak but statistically significant positive regression coefficient for the square of height was seen (β for 10 cm = 0.022; 95% CI = 0.006, 0.038).
Figure 1 presents the mean number of children by height. Some non-linearity was seen in both men and women, but the confidence intervals were wide. When we modeled the number of children using Poisson regression, we found a weak positive but statistically non-significant trend in men (β for 10 cm = 0.006; 95% CI = -0.002, 0.010) whereas in women the association virtually did not exist (β for 10 cm = -0.01; 95% CI = -0.05, 0.03).
We then conducted a multinominal logistic regression model for height quintiles when predicting the transitions in the categories of the number of children born (Table 2). The third quintile and having no children were used as the reference categories in the models (OR = 1.00). In males, the highest probability of having one to two children was in the third quintile.
Model 1 = adjusted for birth year and zygosity; model 2 = additionally adjusted for education; model 3 = additionally adjusted for marital history.
When adjusted for education and marital history, the difference in other categories decreased slightly. When analyzing the probability of having three or more children, the tallest and third quintiles showed the highest probabilities, but these associations were not statistically significant when education and marital history were adjusted for. When we replicated the results in the whole data set, the differences were slightly smaller, but the similar pattern was seen (Table S1 of the Supplementary Material). In women, no statistically significant associations were found between height and number of children.
We then continued the analyses by studying the probability of having any children within twin pairs discordant for height quintile using the conditional logistic regression model (Table 3). Number of children was used in these analyses as a dichotomized variable as only discordant pairs are informative, and thus this analysis has less power than when using all twins as individuals. The probability of having any children was highest in the third quintile in males, and this association decreased slightly when education and marital history were adjusted for. In the whole data set, the differences were again somewhat smaller, but still the similar pattern was seen (Table S2). In women, no association was found.
Model 1 = adjusted for birth year and zygosity; model 2 = additionally adjusted for education; model 3 = additionally adjusted for marital history.
We also conducted similar analyses for marital history (Table 4). In males, the highest probability of having a spouse was in the third category and lowest in the shortest category. The pattern was similar in the analyses using twins as individuals and for discordant twin pairs, but the differences were larger in the latter analysis. When using the whole data set, the results were largely similar (Table S3). In females, no differences were found.
Model 1 = adjusted for birth year and zygosity; model 2 = additionally adjusted for education.
Finally, we analyzed the genetics of height, log AFB, number of children, and their mutual correlations by using genetic twin modeling. Compared with the saturated models, the fit of univariate models was good in men (Δχ26 = 4.6, p = .59 for height and Δχ26 = 8.7, p = .19 for log AFB) and women (Δχ26 = 11.9, p = .06 for height, Δχ26 = 5.5, p = .48 for log AFB and Δχ26 = 2.6, p = .86 for number of children), suggesting that the assumptions of twin modeling were not violated. Only for the number of children in men the fit was non-optimal (Δχ26 = 21.5, p = .002), but was still close to the Bonferroni-corrected p value of six tests (p = .008). To find the best-fitting model, we removed common environmental effect for log AFB (Δχ21 = 2.6, p = .11 in men and Δχ21 = 1.8, p = .18 in women) and for the number of children (Δχ21 = 0, p = 1.00 and Δχ21 = 0.79, p = .37, respectively). For height, this effect was statistically significant (Δχ21 = 5.3, p = .02 and Δχ21 = 6.7, p = .01, respectively), showing that the common environmental component is needed in the model. For females, we also conducted Cholesky decomposition to analyze the correlation between height and log AFB. We were able to drop a specific environmental correlation between log AFB and height (Δχ21 = 0.2, p = .63). In men, no association between height and log AFB was found, and thus we were not able to make similar analysis. Furthermore, because in men there was evidence for non-linear association between height and the number of children and in women no association was found, it was not possible to conduct similar genetic modeling for this association.
Table 5 shows the results for the best-fitting genetic models. Genetic factors were important for height and they explained 78% of the variation in men and 80% of the variation in women. For AFB and the number of children, the specific environmental factors were more important, but still a substantial proportion of the variation of AFB (30% in men and 42% in women) and number of children (42% and 48%, respectively) was explained by genetic factors. In women, the genetic correlation between height and AFB was 0.18, suggesting that 3% of the genetic variation of these two traits is shared.
aLogarithmic transformation was used.
Discussion
Our results do not support the hypothesis that tallness is linearly associated with reproductive success in males. In contrast, lifetime reproductive success was highest in men with average stature. This result remained when adjusting for education as a measure of adulthood socioeconomic status and marital history. Importantly, we also found a similar association when analyzing twin pairs discordant for height. Since twin pairs are the same age and were raised in the same family, this is strong evidence that this association is not affected by childhood family background. Many studies conducted in Western developed countries (Goldstein & Kobyliansky, Reference Goldstein and Kobyliansky1984; Nettle, Reference Nettle2002a; Scott & Bajema, Reference Scott and Bajema1982), in China (Fielding et al., Reference Fielding, Schooling, Adab, Cheng, Lao, Jiang and Lam2008), and in a traditional Gambian population (Sear, Reference Sear2006) have not found a linear association between adult height and lifetime reproductive success in males. Thus, our findings are consistent with previous empirical evidence suggesting that the association between tallness and men's reproductive success is at best weak. The long-lasting hypothesis of the better reproductive success of tall men is thus probably based more on theoretical considerations than on solid empirical evidence. Sexual selection for physical traits is very common in the animal kingdom (Andersson, Reference Andersson1994), and a correlation between reproductive success and large body size is found in many higher primates (Kappeler & van Schaik, Reference Kappeler, van Schaik, Kappeler and van Schaik2004). Thus, it is tempting to seek a similar association also in humans. However, while sexual selection on body size has led to large sexual dimorphism in many higher primates (Plavcan, Reference Plavcan2001), in humans the difference in height between men and women is generally relatively small (Gray & Wolfe, Reference Gray and Wolfe1980). Therefore, it may not be surprising that tallness in males is not strongly associated with reproductive success in humans.
We found in this study that males with average stature had the best reproductive success. The same result has been found for the number of children in an early study of US males who were college students between 1880 and 1912 (Vetta, Reference Vetta1975) and US high-school students who graduated in 1957 (Stulp et al., Reference Stulp, Pollet, Verhults and Buunk2012a). Furthermore, Stulp and colleagues found that in eight of the ten previous studies reporting no association between height and reproductive success in males, there was evidence for the curvilinear association (Stulp et al., Reference Stulp, Pollet, Verhults and Buunk2012a). These studies and our results together suggest that sexual selection may balance average height in the population rather than favor tall or short stature. A possible explanation for this result may be assortative mating for height, which is well known in humans (Spuhler, Reference Spuhler1982). Thus, females are not seeking a partner who is very tall, but modify their preferences according to their own height. A Finnish twin study showed that approximately half of the spousal correlation in height is because of same socioeconomic background — that is, social homogamy — and half because of phenotypic assortment (Silventoinen et al., Reference Silventoinen, Kaprio, Lahelma, Viken and Rose2003). Even when this tendency is not fixed and especially short and tall persons can modify their mate preferences to maximize the number of potential partners (Pawlowski, Reference Pawlowski2003), very short and tall men may still have problems with finding a suitable partner. We found some evidence for this, because the probability of having a spouse was also highest in the men with average stature. However, differences in the reproductive success remained, even when marital history was adjusted for.
In common with a previous Finnish study (Helle, Reference Helle2008), we found that short stature in females was associated with earlier AFB. This finding can be associated with the timing of puberty as short stature has been reported to be associated with earlier age at menarche in several European populations (Onland-Moret et al., Reference Onland-Moret, Peeters, van Gils, Clavel-Chapelon, Key, Tjonneland and Riboli2005). However, the reason for this association is not clear. It has been suggested that the association between height and AFB could be because of a trade-off between allocating energy to continuing growth or to beginning reproduction (Allal et al., Reference Allal, Sear, Prentice and Mace2004). However, in our study cohort, only 1% of women had the first birth under 18 years of age, and even in the shortest quartile of women AFB was 25 years. Since growth virtually ceases in girls before 17 years of age (Thodberg et al., Reference Thodberg, Juul, Lomholt, Martin, Jenni, Caflisch, Kreiborg and Preedy2012), it is likely that only a few women had their first birth while still growing. Furthermore, a recent review suggested that the positive association between female stature and age at puberty is specific to highly developed societies, whereas in agrarian societies the association is reversed (McIntyre & Kacerosky, Reference McIntyre and Kacerosky2011); these results do not support the hypothesis that the energy trade-off would be contributing to this association. In our study cohorts, contraception was also generally used, and at least the large majority of the pregnancies were planned. This may suggest that this association is rather because short stature is associated with a decision to begin reproduction early than because of any causal association. Utilizing genetic twin modeling, we found in this study that the association between height and AFB was because of common genetic factors. The trait correlation between these two traits was, however, only modest, and thus only 3% of the genetic variation was shared between height and AFB. Because all causal associations between these traits would be modeled as part of specific environment, the hypothesis of the energy trade-off between growth and reproduction is not supported by our study (Walker et al., Reference Walker, Gurven, Hill, Migliano, Chagnon, De Souza and Yamauchi2006). One explanation would be that the same hormonal factors, which are strongly genetically regulated, contribute to early puberty, leading to shorter adult stature, and also contribute to stronger childbearing motivation at an early age.
It is noteworthy that not all persons wishing to have children can do so because of biological reasons. The 12-month prevalence of infertility was estimated to be 9% in a large international study, with little difference between developed and less developed countries (Boivin et al., Reference Boivin, Bunting, Collins and Nygren2007). Thus, the proportion of persons with no children is a mix of those desiring them but unable to have them because of social or biological reasons, and those deciding not to have children. Given that better-off persons have more access to health care and are generally taller, this may be a partial confounder despite adjustment for education in our models. Infertility may be due to acquired factors, such as infections of the reproductive organs, or idiopathic biological factors. At the extremes of biological factors are chromosomal abnormalities such as Klinefelters syndrome, who are tall and infertile men, and Turner syndrome, who are short and infertile women (Frühmesser & Kotzot, Reference Frühmesser and Kotzot2011). The influence of hormonal factors can be inherited, such as congenital adrenal hyperplasia, or can be external as the result of hormonal treatments to induce growth or reduce growth in childhood (Chiavaroli et al., Reference Chiavaroli, Liberati, D'Antonio, Masuccio, Capanna, Verrotti and Mohn2010; Hendriks et al., Reference Hendriks, Laven, Valkenburg, Fong, Fauser, de Ridder and Drop2011).
While these specific conditions are rare, genetic variation in the genes underlying the chromosomal abnormalities is common in the population. As some 200 genes are now known to be associated with height in the populations (Lango Allen et al., Reference Lango Allen, Estrada, Lettre, Berndt, Weedon, Rivadeneira and Hirschhorn2010) and copy-number variants particularly associated with short stature have also been identified (Dauber et al., Reference Dauber, Yu, Turchin, Chiang, Meng, Demerath and Hirschhorn2010), it is possible that some of these also affect fertility, given the complex relationships between growth, puberty, and sexual maturation. However, in our data, stature in males was associated in a similar way with having any children than with further parities, and in females no association was found. Potentially different biological and social pathways between height and reproduction emphasize that it should not be expected a priori that height is associated with all parities in a similar way. Alternatively, the number of children could be considered as a categorical variable, allowing that the association between height and probability of having any further children is not independent of pre-existing parity, such as we have done in this study. We found some evidence for this because the association between height and reproduction was more visible when using the number of children as categorized variable than when using the total number of children. This suggests that in higher parities other factors than height become more dominant and thus can confound the association.
Our study has strengths as well as limitations. A major strength of our study is that we had information on a large number of twins that allowed us to optimally take into account the family background and also to analyze the role of genetic and environmental factors behind the association between height and AFB. Furthermore, the study participants were 51–59 years of age at the end of the follow-up, and thus we can safely say that we have information on virtually all participants for the lifetime number of children also in males. Also, the information on children is register based and thus is not prone to reporting bias. A limitation is that twins may not be fully representative to the general population. We found some evidence for this, because a part of twins were living with their co-twin, and in these twins childlessness was more common than in the general population. We removed these pairs from the analyses as they may not be representative of the whole Finnish population. After this exclusion, the fertility rate in our data is very close to the general Finnish population in these cohorts, suggesting that our data of twins are representative (Miettinen, Reference Miettinen2010). When we replicated the analysis in the whole data set, we found that the associations were slightly weaker than in the restricted data. This suggests that other mechanisms may affect the number of children in those twins living with their co-twin.
To conclude, our results do not support the hypothesis that tall stature is associated with higher reproductive success. Instead, we showed that among men average stature may give some advantage for reproduction. Among women, short stature was associated with an early age at reproduction, and this association was because of shared genetic factors. Height is thus associated with reproduction, but these associations and their background may be more complex than previously assumed.
Acknowledgments
This study is supported by the Kone foundation (KS and JN), the Turku Collegium for Science and Medicine (SH), and the Academy of Finland Centre of Excellence in Complex Disease Genetics (grant numbers 213506 and 129680 to JK).
Supplementary Material
To view supplementary material for this article, please visit http://dx.doi.org/10.1017/thg.2012.150.