Which of the following best describes the relationship between genetic and environmental influences on adolescent behavior?

Height variation is known to be determined by both genetic and environmental factors, but a systematic description of how their influences differ by sex, age and global regions is lacking. We conducted an individual-based pooled analysis of 45 twin cohorts from 20 countries, including 180,520 paired measurements at ages 1–19 years. The proportion of height variation explained by shared environmental factors was greatest in early childhood, but these effects remained present until early adulthood. Accordingly, the relative genetic contribution increased with age and was greatest in adolescence (up to 0.83 in boys and 0.76 in girls). Comparing geographic-cultural regions (Europe, North-America and Australia and East-Asia), genetic variance was greatest in North-America and Australia and lowest in East-Asia, but the relative proportion of genetic variation was roughly similar across these regions. Our findings provide further insights into height variation during childhood and adolescence in populations representing different ethnicities and exposed to different environments.

Human height is a classic anthropometric quantitative trait for its ease of measurement, approximately normal distribution and relative stability in adulthood and thus has been the target of extensive research across many fields of science. The study of height has a long standing tradition in genetics; in fact, the field of quantitative genetics was born out of studies of human height in the late 19th and early 20th centuries. Galton1 published data as early as 1886 on the relationship between parent and offspring height and inferred that “when dealing with the transmission of stature from parents to children, the average height of the two parents is, … all we need care to know about them”. Later on, Pearson and Lee2 presented correlations of height between relatives, also providing evidence for the inheritance of height. In 1918 Fisher3 calculated the first heritability estimate of height, i.e. the proportion of total variation explained by genetic variation; in this seminal paper presenting the statistical principles of quantitative genetics, he demonstrated that continuous characters are caused by a combination of many genetic loci with small effects (polygenic inheritance), replacing the blending inheritance hypothesis proposed by Galton. Since then many lines of evidence such as twin, adoption and family studies have estimated the role of genetic factors in the determination of height, showing that it is one of the most heritable human quantitative phenotypes4. Interest in the genetic influences on height was renewed when genetic linkage studies enabled research into genetic effects over the whole genome5 and genome-wide association (GWA) studies allowed identification of loci consistently associated with height in populations of different ancestry6,7,8,9,10.

Beside the genetic factors, a multitude of environmental factors can affect height. They can operate during the whole growth period, but infancy is probably the most sensitive phase regarding external influences11,12. In the presence of adverse environmental conditions, the physical growth of children can decline and even adult height be affected12,13,14. Nutrition and especially lack of dietary protein is universally the most important environmental factor influencing height, but also childhood diseases, in particular infections, can affect growth11. These and other proximate biological determinants are further associated with social and economic conditions manifesting as socio-economic differences in height both within and between populations12.

Although the heritable nature of height has been recognized for more than one hundred years, only a few studies have explored in detail the genetic variation of height during childhood and adolescence. Twin studies have consistently estimated that the heritability of height is lowest (0.2–0.5) in infancy15,16,17, rapidly increases in childhood with varying values15,16,18 and reaches estimates ranging from 0.70 to 0.90 in adolescence and adulthood15,17,19,20,21. However, these studies leave unclear whether environmental factors shared by co-twins, which are generally important in infancy and childhood, persist in adolescence or after the cessation of growth15,16,17,18,19,20,21,22,23. A study in four countries with over 12,000 twin pairs from birth to 19 years of age showed that the effect of shared environment remained up through 12 years and was present again at 16 years15. Somewhat different results were observed in a longitudinal study of two Finnish twin cohorts, which found that common environmental factors affected height at different ages in adolescence and early adulthood20.

Height is also a classic example of a sexually dimorphic trait; on average, men are taller than women in all human populations13. However, much less is known about sex-differences in genetic and environmental contributions to height variation. Greater heritability estimates for males than for females in childhood15 and adulthood21 have been reported. Also sex-specific genetic effects have been found for height, but the results are inconsistent across studies18,19,20,21,24.

Further, a greater mean height has been consistently observed in Western populations as compared with East-Asian populations13, but most studies on the genetic and environmental factors influencing height variation to date are based on Western populations. A multinational study on adolescent twins from eight countries showed that even when the total variation of height was higher in Western populations, the heritability estimates were largely similar between Western and East-Asian populations19. However, these studies did not address the possible differences in the genetic variation pattern between ethnic-cultural groups in childhood and late adolescence/early adulthood.

Using height measures obtained from 45 twin cohorts in 20 countries participating in the COllaborative project of Development of Anthropometrical measures in Twins (CODATwins), we conducted an individual-based analysis of pooled twin cohorts (i) to analyze the genetic and environmental contribution to variation of height from 1 to 19 years of age; (ii) to explore sex-differences in these contributions over each year of age; and (iii) to assess whether this age pattern varies by geographic-cultural region (Europe, North-America and Australia and East-Asia).

Descriptive statistics of height by age and sex for the pooled data (all cohorts together) and by geographic-cultural region are presented in Table 1. Mean height expectedly increased with age in both sexes with the exception of the slight decrease observed at 18/19 years of age, which reflects differences in the distribution of different cohorts within each age group. Mean height was greater in boys than in girls; only at the age of 11 and 12 years were girls slightly taller than boys, reflecting the earlier onset of pubertal growth in girls. The difference in mean height between consecutive age groups was very similar in boys and girls during childhood; these mean height differences started to decrease considerably from 12 years in girls and 14 years in boys. The variation of height increased with age and reached the peak at 12 years in girls and 13 in boys and then decreased slightly. When comparing geographic-cultural regions, mean height was tallest in Europe, somewhat shorter in North-America and Australia and shortest in East-Asia at all ages in boys and girls. The variation of height showed a less clear pattern but was generally greatest in North-America and Australia and lowest in East-Asia.

Table 1 Descriptive statistics of height by age, sex and geographic-cultural region.

Figure 1 presents the proportions of height variation explained by additive genetic, common (shared) environmental and unique environmental factors from 1 to 19 years of age in the pooled data (estimates with 95% confidence intervals (CIs) are available in Supplementary Table S1). The proportion of environmental variation shared by co-twins was greatest at age 1 (0.48 in boys and 0.49 in girls), decreased over childhood and stabilized in adolescence, remaining considerable until 19 years (except at ages 14 and 16 in boys) with values generally lower than 0.2. Accordingly, heritability was lowest at age 1 (0.40 in boys and 0. 38 in girls) increased with age in early and middle childhood (~2–5 and 6–8 years of age, respectively) and was generally greater than 0.7 in late childhood (~9–11 years of age) and adolescence; the greatest heritability estimates were found for boys at ages 14 and 16 (0.83 and 0.82, respectively). The proportion of height variation explained by environmental factors unique to each twin individual, which also includes measurement error, did not show any clear age pattern and was largely similar at all ages (0.05–0.14). In spite of the observed sex differences in the relative variance components at most of ages (See Supplementary Table S2), the age pattern was generally similar in boys and girls; the biggest sex-differences were found in late adolescence when the heritability estimates were slightly greater in boys. The point estimates for the genetic correlations within opposite-sex DZ pairs were generally lower than 0.5 suggesting sex-specific genetic effects. There was a trend for these correlations to be lowest in adolescence, although the largest 95% confidence intervals were estimated at 14 and 16 years of age (Fig. 2).

Figure 1

Proportion of height variation explained by additive genetic (a,b), shared environmental (c,d) and specific environmental (e,f) factors with 95% confidence intervals by age and sex for all cohorts together.

Figure 2

Additive genetic correlations with 95% confidence intervals within opposite-sex DZ pairs by age for all cohorts together.

Univariate models for height were then conducted separately in the three geographic-cultural regions. Only the estimates of additive genetic factors are presented in Fig. 3, but all estimates with 95% CIs are available in Supplementary Table S1. The three geographic-cultural regions showed the general trend of increasing proportion of additive genetic factors with age during childhood. Explained by its largest sample size, the pattern in Europe was practically the same to that observed for all cohorts together, but with slightly greater heritability estimates at most ages. In North-America and Australia and East-Asia, heritability estimates in childhood were generally somewhat lower than in Europe. In East-Asia the pattern in adolescence was not so clear because of the smaller sample size leading to wider 95% CIs. In spite of the roughly similar age patterns, the proportions of height variation explained by genetic and environmental factors were different between the geographic-cultural regions (See Supplementary Table S2). The Chinese National Twin Registry was excluded from these analyses because the heritability estimates in that cohort were substantially lower than in other East-Asian cohorts. When data from this cohort was included in the analyses for East-Asia, the proportion of genetic factors decreased and common environmental factors increased considerably; the change in heritability estimates was from 0.1 to 0.3 units depending on the age group (data available on request).

Figure 3

Proportion of height variation with 95% confidence intervals explained by additive genetic factors separately by age and sex in Europe (a,b), North America and Australia (c,d) and East Asia (e,f).

Finally, we studied how age modifies the genetic and environmental variances of height by using gene-age interaction analysis, with data pooled across all age groups. Figure 4 shows the change in the predicted raw genetic and environmental variances in height as a function of age (parameter estimates with 95% CIs are available in Supplementary Table S3). Additive genetic variation increased steadily from age 1, reached its peak at 14 years in boys and 13 years in girls and then decreased again in the pooled data; however, common and unique environmental variation were largely similar across ages. When stratified by geographic-cultural region, genetic variation was largest in North-America and Australia, somewhat lower in Europe and lowest in East-Asia, particularly for boys. The pattern of genetic variance increasing to a maximum and thereafter decreasing was consistent across the regions. Also common environmental variation was greatest in North-America and Australia, reaching the peak at 10 years in boys and 7 years in girls, whereas in Europe and East-Asia it was similar across ages. Unique environmental variation showed a similar pattern and magnitude in the three geographic-cultural regions. The differences between the regions were highly statistically significant in boys [difference in −2 log-likelihood values (Δ − 2LL) = 1257, difference in degrees of freedom (Δd.f.) = 30, p-value < 0.0001)] and girls (Δ − 2LL = 1364, Δd.f. = 30, p-value < 0.0001). When comparing sexes, in Europe and North-America and Australia there was a trend toward a greater genetic variation for boys than for girls, which increased with age. In East-Asia, however, genetic variation was slightly greater for girls until 14 years of age and for boys in late adolescence.

Figure 4

Change of additive genetic (dash line), common environmental (solid line) and unique environmental (dot line) variance with increasing age in quadratic gene-environment interaction model in Europe, North America and Australia and East Asia.

The present study of 180,520 paired measurements from 86,037 complete twin pairs in 20 countries revealed that environmental factors shared by co-twins contribute to the inter-individual variation in height from infancy to early adulthood. The relative proportion of common environmental factors was greatest during the first years of life, representing almost half of the variation at age 1 and decreased over childhood and adolescence. The interpretation of these results, however, deserves some caution. It has been questioned whether twin studies are suitable for estimating heritability of height in infancy, since early growth patterns in twins differ considerably from singleton growth patterns25. Prenatal environmental factors can act very differently on MZ twins leading to differences in body size within pairs (the most extreme case is the twin-to-twin transfusion syndrome). This is an important issue because in the classical twin design heritability is estimated by comparing the resemblance of MZ and DZ twin pairs and thus body size differences in MZ pairs will result in lower heritability estimates. Since children may take several years to fully catch-up after birth, the high proportion of height variation explained by the shared environment in infancy may still reflect these prenatal environmental factors. Among other possible explanations, it might be that the shared environment represents the effects of gestational age or the effects of the higher measurement error (correlated in twins) at earlier ages.

The influence of the shared environment on height variation up to 19 years, which is consistent with previous studies in adolescents19 and adults21 with enough statistical power to detect this component, suggests that adult height variation reflects childhood living conditions. Studies have shown that the secular trend in adult height occurs during the first two years of life mainly due to increases in leg length26. A plausible explanation is that the period of most rapid growth, when the effect of an adverse environment is strongest, coincides with the period when most growth takes place in the long bones of the legs26. Multinational studies analyzing the genetic and environmental influences on body length segments, particularly leg length, are thus needed to disentangle the aetiology of total height variation. The small but considerable effect of unique environment on height variation, very similar across ages, may partly be due to measurement error, which is modelled as part of unique environmental factors. However, it is likely that it also reflects real environmental factors, for example, different exposure to childhood diseases.

A recent and large meta-analysis of twin correlations and variance components for 17,804 traits carried out separately in four age groups (0–11, 12–17, 18–64 and 65+ years) showed that the heritability estimate of height at 12–17 years was considerably greater than at 0–11 years27. Given the rapid growth that occurs in infancy, childhood and adolescence, in this individual-based pooled analysis we analyzed the heritability of height in one year age groups. We found that genetic contributions increase over childhood with heritability estimates in the range of previous studies in children and adults15,16,18,20,21. GWA studies have identified many common genetic variants for adult height. The most recent GWA meta-analysis in 253,288 individuals of European ancestry identified 697 genome-wide significant SNPs in 423 loci that together explained one-fifth of the heritability for adult height10. Further, in a study using whole-genome sequencing data from 44,126 unrelated individuals, all imputed variants explained 56% of variance for height suggesting that missing heritability is negligible for human height28. However, much less is known on the genetics of height in children. Van der Valk et al.29 found that polygenic scores based on 180 SNPs previously associated with adult height explained 2.95% of the variance of infant length and that of 180 known adult height loci, only 11 were genome-wide significantly associated with infant length.

The pattern of total height variation across ages was largely driven by genetic variance. The most consistent result is the increasing genetic variance with age, reaching its peak at around 13 years in girls and 14 years in boys. After that point, even if mean height continued to increase, genetic variance started to decrease in such a way that in late adolescence the magnitude was similar to that before pubertal events start. Adolescence is characterized by the onset of puberty and the occurrence of growth spurts. Although a secular and population-dependent decline has been observed in the age at onset of pubertal growth spurt and peak height velocity since the mid 1900s30,31, the pubertal height spurt generally begins at age 10–11 years in girls and 11–13 years in boys and reaches peak height velocity at about 12 years and 14 years, respectively13,30. In this study, twins within age groups are at various stages of puberty. In addition to the substantial heritability reported for pubertal timing32, a genome-wide genetic correlation (0.13) between age at menarche and adult height has also been found33. In fact, a genome-wide association meta-analysis showed that five loci associated with pubertal timing impacted multiple aspects of growth, both before and during puberty34. Therefore, it is possible that some of the genetic variance in height at these ages is confounded with genetic variance in pubertal events.

In spite of the largely similar age patterns observed in boys and girls, boys showed somewhat greater heritability estimates and genetic variation, especially in late adolescence. Greater heritability estimates in boys than in girls have previously been reported from birth through 19 years15 and in adulthood21. Moreover, some studies have shown a sex-specific genetic effect on height variation in adolescents19 and adults24. It is clear that both of the sex chromosomes are implicated in determining mean height. Short stature has been demonstrated in females with Turner syndrome who have only one X chromosome35 and taller stature seen in XYY men compared with XY men36. However, sex chromosomes have also been associated with height variation; for example, Gudbjartsson et al.37 identified 27 regions of the genome including a locus on X chromosome that together explained around 3.7% of the population variation in height. In our multinational data, the lowest genetic correlations within opposite-sex DZ pairs were found at 14–16 years of age and again at 18 years, suggesting that sex-specific genes have a role in the genetic variation of height not only during puberty, but also in late adolescence.

Comparison between geographic-cultural regions showed that mean height was greatest in Europe, somewhat shorter in North-America and Australia and shortest in East-Asia, but total variance was largest in North-America and Australia. Accordingly, genetic variation was also greatest in North-America and Australia and lowest in East-Asia. However, the relative proportions of additive and environmental variations were more similar in the different geographic-cultural regions. These results are consistent with a previous comparative twin study which found that the mean and variance of height were larger in Caucasian than in East-Asian populations in adolescence, but the heritability estimates were still at the same level19. An important proportion of the differences in total variances between geographic-cultural regions were attributable to genetic differences. It may be that allelic frequencies and effects of the genes involved in height vary between Europeans, North-Americans and Australians and East-Asians, leading to differences in genetic variation between the three population groups. A recent study across 14 European countries found that many independent loci contribute to population genetic differences in height and estimated that these differences account for 24% of the captured additive genetic variance38. However, a major part of the differences in genetic variation may also be because of gene-environment interactions modelled as part of the additive genetic component in our model. That is, the higher genetic variation observed in Caucasians could arise because there is a set of genes expressed more strongly in Western environments. For example, a study of adults of Japanese descent living in the United States and native Japanese found that Japanese men and women were shorter than Japanese-Americans, suggesting that environmental factors play a role in physical growth39. Analyzing this question in detail would require collection of twins or GWA studies in unrelated individuals with East-Asian origin living in a Western environment.

The study in Caucasian and East-Asian populations showed that approximately 91% of the differences in the total variance between these two population groups was attributable to genetic variances19. However, our study found that shared environmental variance also differed between geographic-cultural regions. The lower shared environmental variance observed in East Asian girls and greater in North-America and Australia during childhood may reflect cultural differences in terms of nutrition and other environmental resources. It is also important to note that we limited our East-Asian cohorts to affluent East-Asian populations including the Shandong and Guangdong provinces but excluding poorer areas of China. As reported previously, the heritability estimates of height were considerably lower and common environmental estimates higher in the poorer areas40, which may indicate larger differences between families in nutrition and infection history in these areas of China. This emphasizes the need to collect data on twins living under different environmental exposures.

The main strength of the present study is the very large sample size of our multinational database of twin cohorts, with height data from 1 through 19 years of age, allowing a more detailed investigation of the genetic and environmental contributions to individual differences in height during childhood and adolescence than in the previous studies. Twin participants are from 20 different countries, thereby making it possible to stratify the analyses by regions representing different ethnicities and environments. Important advantages of individual-based data are better opportunities for statistical modelling and lack of publication bias. However, our study also has limitations. The equal-environment assumption, upon which twin methodology is based, assumes that MZ and DZ twins are equally exposed to environmental factors relevant to the outcome. If equal-environment assumption is violated, it should be seen as differences in variances between MZ and DZ twins, but we did not find such evidence. In the classical twin design phenotypic assortment increases DZ correlations and thus inflates the common environmental component when not accounted for in the modelling. Assortative mating is well recognized for height and when the potential underestimation of heritability estimates was corrected using a sample of twins and their parents41, these authors showed that doing so increased the heritability estimates from 0.75 to 0.85. In our database we do not have information on parental height and thus could not take into account assortative mating, which may thus explain part of the shared environmental variation. A recent study showed that increased homozygosity, which is influenced by inbreeding, was associated with decreased height and that the effect sizes were similar across different continental groups and populations with different degrees of genome-wide homozygosity42. These authors thus suggested that homozygosity, rather than confounding as a result of environmental or additive genetic effects, directly contributes to phenotypic variance42. Further, most of the height measures were self-reported43, which are prone to error and can bias our analyses toward lower heritability estimates and higher estimates of unique environmental effects. Finally, countries and/or ethnic-cultural regions are not equally represented and the database is heavily weighted towards populations following the Westernized lifestyle; even when the large majority of the twin cohorts in the world participated in this project, our data still had limited power for East-Asia especially in adolescence. An even bigger problem is that there are few data available from South-Asia, Middle-East and Africa and no data from South-America. This demonstrates the need for new data collections in these regions.

Our findings provide further insights into height variation during childhood and adolescence in populations representing different ethnicities and exposed to different environments. Worthwhile objectives for future research are to study whether the same genetic and environmental factors contributing to height variation operate throughout time or new genes or new environmental factors start to operate at different ages and to analyze the heritability of growth in height. Further, a major challenge in future studies with more information on birth and pregnancy related variables is to explore the reasons for the low heritability of height at young ages.

In conclusion, environmental factors shared by co-twins exert their strongest influence on height variation in childhood, but these effects remain until the onset of adulthood. Genetic variation in height increased steadily during childhood and reached its peak at around 13 years in girls and 14 years in boys, which may be confounded with genetic variation in pubertal events. Especially in adolescence, there was a trend toward somewhat greater genetic variation in boys than in girls and part of the genetic variation of height was sex-specific. Genetic variation of height was larger in North-America and Australia and Europe compared with East-Asia, but the relative proportions of genetic and environmental variations between these three geographic-cultural regions were roughly similar. These findings suggest that, in spite of different ethnicities and environmental exposures, genetic factors play a major role on height variation in adolescence and early adulthood, but environmental factors shared by co-twins are also important.

All participants were volunteers and gave their informed consent when participating in the study. No experimental data were asked and thus we did not ask ethical approval. Only a limited set of observational variables and anonymized data were delivered to the data management centre at University of Helsinki. The pooled analysis was approved by the ethical committee of Department of Public Health, University of Helsinki and the methods were carried out in accordance with the approved guidelines.

Sample

This study is based on the data from the CODATwins project described elsewhere43. Briefly, the CODATwins project was intended to recruit all twin projects in the world with information on height and weight measurements. For the present analyses, we selected height measurements at ages from 0.5 to 19.5 years (n = 420,707). Age was classified to single-year age groups (e.g., age 1 refers to 0.5–1.5 years range). Impossible values and outliers were checked by visual inspection of histograms for each age and sex group and were removed to obtain an approximately normal distribution (0.3% of the measurements). Since individuals in longitudinal studies have more than one measurement over time, analyses were restricted to one observation per individual in each age group. Analyses were additionally restricted to having at least 50 measurements per cohort. Finally we had data from 45 cohorts in 20 countries: one cohort from Africa (Guinea Bissau Twin Study), two cohorts from Australia (Peri/Postnatal Epigenetic Twins Study and Queensland Twin Register), seven cohorts from East-Asia (Chinese National Twin Registry, Guangzhou Twin Eye Study, Japanese Twin Cohort, Mongolian Twin Registry, Qingdao Twin Registry of Children, South Korea Twin Registry and West Japan Twins and Higher Order Multiple Births Registry), 20 cohorts from Europe (Adult Netherlands Twin Registry, Berlin Twin Register, Bielefeld Longitudinal Study of Adult Twins, Danish Twin Cohort, East Flanders Prospective Twin Survey, Finnish Older Twin Cohort, FinnTwin12, FinnTwin16, Gemini Study, Genesis 12–19 Study, Hungarian Twin Registry, Italian Twin Registry, Portugal Twin Cohort, Swedish Twin Cohorts, Swedish Young Male Twins Study of Adults, Swedish Young Male Twins Study of Children, TCHAD-study, Twins Early Developmental Study, TwinsUK and Young Netherlands Twin Registry), two cohorts from Middle-East (Longitudinal Israeli Study of Twins and Turkish Twin Study) and 13 cohorts from North-America (Boston University Twin Project, California Twin Program, Colorado Twin Registry, Michigan Twins Study, Mid Atlantic Twin Registry, Minnesota Twin Family Study, NAS-NRC Twin Registry, Quebec Newborn Twin Study, SRI-international, Texas Twin Project, University of British Columbia Twin Project, University of Southern California Twin Study and University of Washington Twin Registry). A more detailed description of the participating twin cohorts was presented previously43.

During the course of the study, we found that the heritability estimates of height were substantially lower in the Chinese National Twin Registry than in other East-Asian cohorts, as also reported previously40. Because of this heterogeneity, we did not include the data from the Chinese National Twin Registry in the reported analyses but tested how it would change the results in East-Asia. The final database comprised information on 86,037 different complete twin pairs with a total of 180,520 paired measurements (39% monozygotic (MZ), 34% same- sex dizygotic (SSDZ) and 27% opposite-sex dizygotic (OSDZ) twin pairs); that is, since some twin pairs have measurements at different ages, our database is based on measurement pairs. In order to analyze possible ethnic-cultural differences in the genetic and environmental contribution on height, cohorts were pooled in three groups according to their geographical and cultural characteristics: Europe (20 cohorts), North-America and Australia (15 cohorts) and East-Asia (6 cohorts) with 131,856, 29,856 and 17,924 paired measurements, respectively. In the additional analyses including the Chinese National Twin Registry, the number of pairs in East-Asia was 27,067. The cohort from Africa and the two from Middle-East were not included in these sub-analyses by geographic-cultural region because the data is too small to study these two areas separately.

Statistical analyses

To analyze genetic and environmental influences on the variation of height, we used classic twin modelling based on linear structural equations44. Briefly, the analysis is based on the fact that MZ twins share the same gene sequence, whereas DZ twins share, on average, 50% of their genes identical-by-descent. On this basis, it is possible to divide the total variation of height into variance due to additive genetic effects (A: correlated 1.0 for MZ and 0.5 for DZ pairs), dominance genetic effects (D: 1.0 for MZ and 0.25 for DZ pairs), common (shared) environmental effects (C: by definition, correlated 1.0 for MZ and DZ pairs) and unique (non-shared) environmental effects (E: by definition, uncorrelated in MZ and DZ pairs). However, since our data included only twins reared together, we cannot simultaneously estimate common environmental and dominance genetic effects. All genetic models were fitted by the OpenMx package (version 2.0.1) in the R statistical platform45 using the maximum likelihood method.

Prior to conducting the modelling, height values were adjusted using linear regressions and the resulting residuals were used as input phenotypes for the following analyses. Adjustment was carried out for birth year, exact age at the time of the measurement and study cohort within one-year and sex groups in univariate analyses and for birth year and cohort in gene-environment interaction analyses.

The ACE sex-limitation model was selected as a starting point of the univariate modelling based on the following criteria: (i) MZ within-pair correlations were clearly higher than DZ correlations consistent with the influence of genetic effects, (ii) the magnitude of the difference between MZ and DZ correlations (rDZ > 1/2 rMZ) indicated the presence of common environmental effects and (iii) the lower within-pair correlations for OSDZ than for SSDZ twins observed at most ages suggested the presence of sex-specific genetic effects (See Supplementary Table S4). Previous findings from this international database showed that both male and female DZ twins have greater height than MZ twins in these age groups46 and thus different means for MZ and DZ twins were allowed. The fit of the univariate models for height at each one-year age group is shown in Supplementary Table S2. In the present study, the equal-environment assumption was tested by comparing the ACE model to the saturated model. The fit of the model after Bonferroni correction of multiple testing did not worsen at most ages, which suggested that the assumption of equality of variances between MZ and DZ twins was not violated. When fixing A, C and E parameters to be the same in boys and girls, the fit of the model was poorer at many ages, particularly from the beginning of adolescence, suggesting that these variance components differ between sexes. We additionally fitted a scale model allowing for different sizes of variance components but fixing the relative size of these components to be equal. Since this model also showed differences, we decided to present the results separately for boys and girls. Sex-specific genetic effects were considerable at some ages and thus all modelling results are presented in sex-limited form for consistency. Finally, comparative model fitting revealed that the C parameter could be not excluded from the model without a significant (P < 0.001) deterioration in fit for all ages.

In order to study how age modifies the genetic and environmental variances of height, we additionally conducted gene–environment interaction modelling using age as an environmental modification factor47. The advantage of the gene-environment interaction model, as compared to fitting a series of univariate models, is that it estimates a fewer number of parameters and thus has more statistical power to analyze the age patterns in genetic and environmental variances. The parameter estimates differed in boys and girls (Δ − 2LL = 407, Δd.f. = 15, p-value < 0.0001) and thus sex-specific models were conducted. Because the size of the sex-specific genetic effect varied according to age, only same-sex pairs were included in these age-moderation analyses. In addition to linear effects of age, quadratic age effects were also included on the variance components since they were highly statistically significant in boys (Δ − 2LL = 3298, Δd.f. = 3, p-value < 0.0001) and girls (Δ − 2LL = 4091, Δd.f. = 3, p-value < 0.0001).

How to cite this article: Jelenkovic, A. et al. Genetic and environmental influences on height from infancy to early adulthood: An individual-based pooled analysis of 45 twin cohorts. Sci. Rep. 6, 28496; doi: 10.1038/srep28496 (2016).

• Galton, F. Regression towards mediocrity in heriditary stature. Journal of the Anthropological Institute 15, 246–262 (1886).

  Google Scholar 

• Pearson, K. & Lee, A. On the laws on inheritance in man. Biometrika 2, 356–462 (1903).

  Google Scholar 

• Fisher, R. A. The correlation between relatives on the supposition of Mendelian inheritance. Transactions of the Royal Society of Edinburgh 52, 399–433 (1918).

  Google Scholar 

• Silventoinen, K. Determinants of variation in adult body height. J. Biosoc. Sci. 35, 263–285 (2003).

  PubMed  Google Scholar 

• Perola, M. et al. Combined genome scans for body stature in 6,602 European twins: evidence for common Caucasian loci. PLoS Genet. 3, e97 (2007).

  PubMed  PubMed Central  Google Scholar 

• Cho, Y. S. et al. A large-scale genome-wide association study of Asian populations uncovers genetic factors influencing eight quantitative traits. Nat. Genet. 41, 527–534 (2009).

  CAS  ADS  PubMed  Google Scholar 

• Hao, Y. et al. Genome-wide association study in Han Chinese identifies three novel loci for human height. Hum. Genet. 132, 681–689 (2013).

  PubMed  Google Scholar 

• Lango Allen, H. et al. Hundreds of variants clustered in genomic loci and biological pathways affect human height. Nature 467, 832–838 (2010).

  CAS  ADS  PubMed  PubMed Central  Google Scholar 

• N’Diaye, A. et al. Identification, replication and fine-mapping of Loci associated with adult height in individuals of african ancestry. PLoS Genet. 7, e1002298 (2011).

  PubMed  PubMed Central  Google Scholar 

• Wood, A. R. et al. Defining the role of common variation in the genomic and biological architecture of adult human height. Nat. Genet. 46, 1173–1186 (2014).

  CAS  PubMed  PubMed Central  Google Scholar 

• Bozzoli, C., Deaton, A. & Quintana-Domeque, C. Adult height and childhood disease. Demography 46, 647–669 (2009).

  PubMed  PubMed Central  Google Scholar 

• Steckel, R. H. Heights and human welfare: Recent developments and new directions. Explorations in Economic History 46, 1–23 (2009).

  Google Scholar 

• Eveleth, P. B. & Tanner, J. M. Worldwide variation in human growth (Cambridge University Press, Cambridge, 1990).

• Bogin, B. The growth of humanity (Wiley-Liss, New York, 2001).

• Dubois, L. et al. Genetic and environmental contributions to weight, height and BMI from birth to 19 years of age: an international study of over 12,000 twin pairs. PLoS One 7, e30153 (2012).

  CAS  ADS  PubMed  PubMed Central  Google Scholar 

• Mook-Kanamori, D. O. et al. Heritability estimates of body size in fetal life and early childhood. PLoS One 7, e39901 (2012).

  CAS  ADS  PubMed  PubMed Central  Google Scholar 

• Silventoinen, K. et al. Genetic regulation of growth from birth to 18 years of age: the Swedish young male twins study. Am. J. Hum. Biol. 20, 292–298 (2008).

  PubMed  Google Scholar 

• Silventoinen, K. et al. Genetic regulation of growth in height and weight from 3 to 12 years of age: a longitudinal study of Dutch twin children. Twin Res. Hum. Genet. 10, 354–363 (2007).

  PubMed  Google Scholar 

• Hur, Y. M. et al. Genetic influences on the difference in variability of height, weight and body mass index between Caucasian and East Asian adolescent twins. Int. J. Obes. (Lond) 32, 1455–1467 (2008).

  Google Scholar 

• Jelenkovic, A. et al. Genetic and environmental influences on growth from late childhood to adulthood: a longitudinal study of two Finnish twin cohorts. Am. J. Hum. Biol. 23, 764–773 (2011).

  PubMed  Google Scholar 

• Silventoinen, K. et al. Heritability of adult body height: a comparative study of twin cohorts in eight countries. Twin Res. 6, 399–408 (2003).

  PubMed  Google Scholar 

• Dubois, L. et al. Genetic and environmental influences on body size in early childhood: a twin birth-cohort study. Twin Res. Hum. Genet. 10, 479–485 (2007).

  PubMed  Google Scholar 

• Estourgie-van Burk, G. F., Bartels, M., van Beijsterveldt, T. C., Delemarre-van de Waal, H. A. & Boomsma, D. I. Body size in five-year-old twins: heritability and comparison to singleton standards. Twin Res. Hum. Genet. 9, 646–655 (2006).

  PubMed  Google Scholar 

• Silventoinen, K., Kaprio, J., Lahelma, E., Viken, R. J. & Rose, R. J. Sex differences in genetic and environmental factors contributing to body-height. Twin Res. 4, 25–29 (2001).

  CAS  PubMed  Google Scholar 

• Phillips, D. I. Twin studies in medical research: can they tell us whether diseases are genetically determined? Lancet 341, 1008–1009 (1993).

  CAS  PubMed  Google Scholar 

• Cole, T. J. The secular trend in human physical growth: a biological view. Econ. Hum. Biol. 1, 161–168 (2003).

  CAS  PubMed  Google Scholar 

• Polderman, T. J. et al. Meta-analysis of the heritability of human traits based on fifty years of twin studies. Nat. Genet. 47, 702–709 (2015).

  CAS  PubMed  Google Scholar 

• Yang, J. et al. Genetic variance estimation with imputed variants finds negligible missing heritability for human height and body mass index. Nat. Genet. 47, 1114–1120 (2015).

  CAS  PubMed  PubMed Central  Google Scholar 

• van der Valk, R. J. et al. A novel common variant in DCST2 is associated with length in early life and height in adulthood. Hum. Mol. Genet. 24, 1155–1168 (2015).

  CAS  PubMed  Google Scholar 

• Aksglaede, L., Olsen, L. W., Sørensen, T. I. & Juul, A. Forty years trends in timing of pubertal growth spurt in 157,000 Danish school children. PLoS One 3, e2728 (2008).

  ADS  PubMed  PubMed Central  Google Scholar 

• Hermanussen, M. Studying human growth and development (ed Hermanussen, M. ) (Schweizerbart Science Publishers, Stuttgart, 2013).

• Wehkalampi, K. et al. Genetic and environmental influences on pubertal timing assessed by height growth. Am. J. Hum. Biol. 20, 417–423 (2008).

  PubMed  PubMed Central  Google Scholar 

• Bulik-Sullivan, B. et al. An atlas of genetic correlations across human diseases and traits. Nat. Genet. 47, 1236–1241 (2015).

  CAS  PubMed  PubMed Central  Google Scholar 

• Cousminer, D. L. et al. Genome-wide association and longitudinal analyses reveal genetic loci linking pubertal height growth, pubertal timing and childhood adiposity. Hum. Mol. Genet. 22, 2735–2747 (2013).

  CAS  PubMed  PubMed Central  Google Scholar 

• Rao, E. et al. Pseudoautosomal deletions encompassing a novel homeobox gene cause growth failure in idiopathic short stature and Turner syndrome. Nat. Genet. 16, 54–63 (1997).

  CAS  PubMed  Google Scholar 

• Ratcliffe, S. G., Pan, H. & McKie, M. Growth during puberty in the XYY boy. Ann. Hum. Biol. 19, 579–587 (1992).

  CAS  PubMed  Google Scholar 

• Gudbjartsson, D. F. et al. Many sequence variants affecting diversity of adult human height. Nat. Genet. 40, 609–615 (2008).

  CAS  PubMed  Google Scholar 

• Robinson, M. R. et al. Population genetic differentiation of height and body mass index across Europe. Nat. Genet. 47, 1357–1362 (2015).

  CAS  PubMed  PubMed Central  Google Scholar 

• Fujiwara, S. et al. Differences in health characteristics between native Japanese and Japanese-Americans. J. Cross Cult. Gerontol. 14, 273–287 (1999).

  CAS  PubMed  Google Scholar 

• Liu, Q. et al. Genetic and Environmental Effects on Weight, Height and BMI Under 18 Years in a Chinese Population-Based Twin Sample. Twin Res. Hum. Genet. 18, 571–580 (2015).

  PubMed  Google Scholar 

• Zietsch, B. P., Verweij, K. J., Heath, A. C. & Martin, N. G. Variation in human mate choice: simultaneously investigating heritability, parental influence, sexual imprinting and assortative mating. Am. Nat. 177, 605–616 (2011).

  PubMed  PubMed Central  Google Scholar 

• Joshi, P. K. et al. Directional dominance on stature and cognition in diverse human populations. Nature 523, 459–462 (2015).

  PubMed  PubMed Central  Google Scholar 

• Silventoinen, K. et al. The CODATwins Project: The Cohort Description of Collaborative Project of Development of Anthropometrical Measures in Twins to Study Macro-Environmental Variation in Genetic and Environmental Effects on Anthropometric Traits. Twin Res. Hum. Genet. 18, 348–360 (2015).

  PubMed  PubMed Central  Google Scholar 

• Neale, M. C. & Cardon, L. R. Methodology for genetic studies of twins and families (Kluwer Academic Publishers, Dordrecht, 1992).

• Boker, S. et al. OpenMx: An Open Source Extended Structural Equation Modeling Framework. Psychometrika 76, 306–317 (2011).

  MathSciNet  PubMed  PubMed Central  MATH  Google Scholar 

• Jelenkovic, A. et al. Zygosity Differences in Height and Body Mass Index of Twins From Infancy to Old Age: A Study of the CODATwins Project. Twin Res. Hum. Genet. 18, 557–570 (2015).

  PubMed  PubMed Central  Google Scholar 

• Purcell, S. Variance components models for gene-environment interaction in twin analysis. Twin Res. 5, 554–571 (2002).

  PubMed  Google Scholar 

Download references

This study was conducted within the CODATwins project (Academy of Finland #266592). Support for participating twin projects: The Boston University Twin Project is funded by grants (#R01 HD068435 #R01 MH062375) from the National Institutes of Health to K. Saudino. California Twin Program was supported by The California Tobacco-Related Disease Research Program (7RT-0134H, 8RT-0107H, 6RT-0354H) and the National Institutes of Health (1R01ESO15150-01). Chinese National Twin Registry is funded by Special Fund for Health Scientific Research in the Public Welfare (Project No: 201502006), China. Colorado Twin Registry is funded by NIDA funded center grant DA011015, & Longitudinal Twin Study HD10333; Author Huibregtse is supported by to 5T32DA017637-11. Danish Twin Registry is supported by the National Program for Research Infrastructure 2007 from the Danish Agency for Science, Technology and Innovation, The Research Council for Health and Disease, the Velux Foundation and the US National Institute of Health (P01 AG08761). Since its origin the East Flanders Prospective Survey has been partly supported by grants from the Fund of Scientific Research, Flanders and Twins, a non-profit Association for Scientific Research in Multiple Births (Belgium). Data collection and analyses in Finnish twin cohorts have been supported by ENGAGE–European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145 and AA-09203 to RJ Rose, the Academy of Finland Center of Excellence in Complex Disease Genetics (grant numbers: 213506, 129680) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 265240, 263278 and 264146 to J. Kaprio). K. Silventoinen is supported by Osaka University’s International Joint Research Promotion Program. Gemini was supported by a grant from Cancer Research UK (C1418/A7974). Waves 1–3 of Genesis 12–19 were funded by the WT Grant Foundation, the University of London Central Research fund and a Medical Research Council Training Fellowship (G81/343) and Career Development Award (G120/635) to Thalia C. Eley. Wave 4 was supported by grants from the Economic and Social Research Council (RES-000-22-2206) and the Institute of Social Psychiatry (06/07–11) to Alice M. Gregory who was also supported at that time by a Leverhulme Research Fellowship (RF/2/RFG/2008/0145). Wave 5 was supported by funding to Alice M. Gregory from Goldsmiths, University of London. Guangzhou Twin Eye Study is supported by National Natural Science Foundation of China (grant #81125007). Anthropometric measurements of the Hungarian twins were supported by Medexpert Ltd., Budapest, Hungary. Longitudinal Israeli Study of Twins was funded by the Starting Grant no. 240994 from the European Research Council (ERC) to Ariel Knafo. The Michigan State University Twin Registry has been supported by Michigan State University, as well as grants R01-MH081813, R01-MH0820-54, R01-MH092377-02, R21-MH070542-01, R03-MH63851-01 from the National Institute of Mental Health (NIMH), R01-HD066040 from the Eunice Kennedy Shriver National Institute for Child Health and Human Development (NICHD) and 11-SPG-2518 from the MSU Foundation. The content of this manuscript is solely the responsibility of the authors and does not necessarily represent the official views of the NIMH, the NICHD, or the National Institutes of Health. The NAS-NRC Twin Registry acknowledges financial support from the National Institutes of Health grant number R21 AG039572. Netherlands Twin Register acknowledges the Netherlands Organization for Scientific Research (NWO) and MagW/ZonMW grants 904-61-090, 985-10-002, 912-10-020, 904-61-193,480-04-004, 463-06-001, 451-04-034, 400-05-717, Addiction-31160008, Middelgroot-911-09-032, Spinozapremie 56-464-14192; VU University’s Institute for Health and Care Research (EMGO+); the European Research Council (ERC-230374), the Avera Institute, Sioux Falls, South Dakota (USA). PETS was supported by grants from the Australian National Health and Medical Research Council (grant numbers 437015 and 607358 to J.C. and R.S.), the Bonnie Babes Foundation (grant number BBF20704 to J.M.C.), the Financial Markets Foundation for Children (grant no. 032-2007 to J.M.C.) and by the Victorian Government’s Operational Infrastructure Support Program. Madeira data comes from the following project: Genetic and environmental influences on physical activity, fitness and health: the Madeira family study Project reference: POCI/DES/56834/2004 Founded by the Portuguese agency for research (The Foundation for Science and Technology [FCT]). The Quebec Newborn Twin Study acknowledges financial support from the Fonds Québécois de la Recherche sur la Société et la Culture, the Fonds de la Recherche en Santé du Québec, the Social Science and Humanities Research Council of Canada, the National Health Research Development Program, the Canadian Institutes for Health Research, Sainte-Justine Hospital’s Research Center and the Canada Research Chair Program (Michel Boivin). South Korea Twin Registry is supported by National Research Foundation of Korea (NRF-371-2011-1 B00047). The Texas Twin Project is currently funded by grants AA023322 and HD081437 from the National Institutes of Health. The Twins Early Development Study (TEDS) is supported by a program grant (G0901245) from the UK Medical Research Council and the work on obesity in TEDS is supported in part by a grant from the UK Biotechnology and Biological Sciences Research Council (31/D19086). S.Y. Öncel and F. Aliev are supported by Kırıkkale University Research Grant: KKU, 2009/43 and TUBITAK grant 114C117. TwinsUK was funded by the Wellcome Trust; European Community’s Seventh Framework Programme (FP7/2007-2013). The study also receives support from the National Institute for Health Research (NIHR) BioResource Clinical Research Facility and Biomedical Research Centre based at Guy’s and St Thomas’ NHS Foundation Trust and King’s College London. The University of Southern California Twin Study is funded by a grant from the National Institute of Mental Health (R01 MH58354). Washington State Twin Registry (formerly the University of Washington Twin Registry) was supported in part by grant NIH RC2 HL103416 (D. Buchwald, PI). The West Japan Twins and Higher Order Multiple Births Registry was supported by Grant-in-Aid for Scientific Research (B) (grant number 15H05105) from the Japan Society for the Promotion of Science.

Author notes

1. Department of Social Research, University of Helsinki, Helsinki, Finland

   Aline Jelenkovic, Reijo Sund, Sari Aaltonen & Karri Silventoinen

2. Department of Genetics, Physical Anthropology and Animal Physiology, University of the Basque Country UPV/EHU, Leioa, Spain

   Aline Jelenkovic & Esther Rebato

3. Department of Education, Mokpo National University, Jeonnam, South Korea

   Yoon-Mi Hur

4. Department of Public Health Nursing, Osaka City University, Osaka, Japan

   Yoshie Yokoyama

5. Department of Public Health, The Danish Twin Registry, Epidemiology, Biostatistics & Biodemography, University of Southern Denmark, Odense, Denmark

   Jacob v. B. Hjelmborg, Sören Möller, Axel Skytthe & Kaare Christensen

6. Osaka University Graduate School of Medicine, Osaka University, Osaka, Japan

   Chika Honda & Karri Silventoinen

7. Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Stockholm, Sweden

   Patrik K. E. Magnusson, Nancy L. Pedersen, Margaret Gatz & Paul Lichtenstein

8. Department of Health Science, Ishikawa Prefectural Nursing University, Kahoku, Ishikawa, Japan

   Syuichi Ooki

9. Department of Public Health, University of Helsinki, Helsinki, Finland

   Sari Aaltonen, Kauko Heikkilä & Jaakko Kaprio

10. Istituto Superiore di Sanità - National Center for Epidemiology, Surveillance and Health Promotion, Rome, Italy

    Maria A. Stazi, Corrado Fagnani & Cristina D’Ippolito

11. Department of Physical Education and Sport, University of Madeira, Funchal, Portugal

    Duarte L. Freitas

12. CIFI2D, Faculty of Sport, Porto, University of Porto, Portugal

    José Antonio Maia

13. Department of Noncommunicable Diseases Prevention, Qingdao Centers for Disease Control and Prevention, Qingdao, China

    Fuling Ji, Feng Ning & Zengchang Pang

14. HealthTwiSt GmbH, Berlin, Germany

    Andreas Busjahn

15. Department of Psychology, Bielefeld University, Bielefeld, Germany

    Christian Kandler

16. Department of Psychological and Brain Sciencies, Boston University, Boston, MA, USA

    Kimberly J. Saudino

17. Department of Psychiatry, University of British Columbia, Vancouver, BC, Canada

    Kerry L. Jang

18. Department of Preventive Medicine, Keck School of Medicine of USC, University of Southern California, Los Angeles, California, USA

    Wendy Cozen, Amie E. Hwang & Thomas M. Mack

19. USC Norris Comprehensive Cancer Center, Los Angeles, California, USA

    Wendy Cozen & Thomas M. Mack

20. Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China

    Wenjing Gao, Canqing Yu & Liming Li

21. Institute for Behavioral Genetics, University of Colorado, Boulder, Colorado, USA

    Robin P. Corley & Brooke M. Huibregtse

22. Centre of Human Genetics, University Hospitals Leuven, Leuven, Belgium

    Catherine A. Derom & Robert F. Vlietinck

23. Department of Obstetrics and Gynaecology, Ghent University Hospitals, Ghent, Belgium

    Catherine A. Derom

24. The Charles Bronfman Institute for Personalized Medicine, The Mindich Child Health and Development Institute, Icahn School of Medicine at Mount Sinai, New York, NY, USA

    Ruth J. F. Loos

25. Department of Epidemiology and Public Health, Health Behaviour Research Centre, Institute of Epidemiology and Health Care, University College London, London, UK

    Jane Wardle, Clare H. Llewellyn & Abigail Fisher

26. King’s College London, MRC Social, Genetic & Developmental Psychiatry Centre, Institute of Psychiatry, Psychology & Neuroscience, London, UK

    Tom A. McAdams, Thalia C. Eley & Robert Plomin

27. Department of Psychology, Goldsmiths, University of London, London, UK

    Alice M. Gregory

28. State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou, China

    Mingguang He & Xiaohu Ding

29. Centre for Eye Research Australia, University of Melbourne, Melbourne, Australia

    Mingguang He

30. Bandim Health Project, INDEPTH Network, Bissau, Guinea-Bissau

    Morten Bjerregaard-Andersen

31. Research Center for Vitamins and Vaccines, Statens Serum Institute, Copenhagen, Denmark

    Morten Bjerregaard-Andersen

32. Department of Endocrinology, Odense University Hospital, Odense, Denmark

    Morten Bjerregaard-Andersen & Henning Beck-Nielsen

33. Department of Infectious Diseases, Odense University Hospital, Odense, Denmark

    Morten Sodemann

34. Department of Radiology and Oncotherapy, Semmelweis University, Budapest, Hungary

    Adam D. Tarnoki & David L. Tarnoki

35. Hungarian Twin Registry, Budapest, Hungary

    Adam D. Tarnoki & David L. Tarnoki

36. The Hebrew University of Jerusalem, Jerusalem, Israel

    Ariel Knafo-Noam & Lior Abramson

37. Hadassah Hospital Obstetrics and Gynecology Department, Hebrew University Medical School, Jerusalem, Israel

    David Mankuta

38. Michigan State University, East Lansing, Michigan, USA

    S. Alexandra Burt & Kelly L. Klump

39. Department of Human and Molecular Genetics, Virginia Institute for Psychiatric and Behavioral Genetics, Virginia Commonwealth University, Richmond, Virginia, USA

    Judy L. Silberg & Lindon J. Eaves

40. Department of Human and Molecular Genetics, Psychiatry & Massey Cancer Center, Virginia Commonwealth University, Richmond, Virginia, USA

    Hermine H. Maes

41. Department of Psychology, University of Minnesota, Minneapolis, MN, USA

    Robert F. Krueger, Matt McGue & Shandell Pahlen

42. Department of Psychology, University of Southern California, Los Angeles, CA, USA

    Margaret Gatz, Laura A. Baker & Catherine Tuvblad

43. Institute of Medicine, National Academy of Sciences, Washington, DC, USA

    David A. Butler

44. Department of Biological Psychology, VU University Amsterdam, Amsterdam, The Netherlands

    Meike Bartels, Toos C. E. M. van Beijsterveldt, Gonneke Willemsen & Dorret I. Boomsma

45. Murdoch Childrens Research Institute, Royal Children’s Hospital, Parkville, Victoria, Australia

    Jeffrey M. Craig & Richard Saffery

46. Department of Paediatrics, University of Melbourne, Parkville, Victoria, Australia

    Jeffrey M. Craig & Richard Saffery

47. School of Epidemiology, Public Health and Preventive Medicine, University of Ottawa, Ottawa, Ontario, Canada

    Lise Dubois

48. École de psychologie, Université Laval, Québec, Canada

    Michel Boivin & Ginette Dionne

49. Institute of Genetic, Neurobiological and Social Foundations of Child Development, Tomsk State University, Russian Federation

    Michel Boivin

50. Département de psychologie, Université du Québec à Montréal, Montréal, Québec, Canada

    Mara Brendgen

51. École de psychoéducation, Université de Montréal, Montréal, Québec, Canada

    Frank Vitaro

52. Genetic Epidemiology Department, QIMR Berghofer Medical Research Institute, Brisbane, Australia

    Nicholas G. Martin & Sarah E. Medland

53. Molecular Epidemiology Department, QIMR Berghofer Medical Research Institute, Brisbane, Australia

    Grant W. Montgomery

54. Department of Medicine, Stanford Prevention Research Center, Stanford University School of Medicine, Stanford, CA, USA

    Gary E. Swan

55. Center for Health Sciences, SRI International, Menlo Park, CA, USA

    Ruth Krasnow

56. Department of Public Health Sciences, Karolinska Institutet, Stockholm, Sweden

    Per Tynelius & Finn Rasmussen

57. MRC Integrative Epidemiology Unit, University of Bristol, Bristol, UK

    Claire M. A. Haworth

58. Healthy Twin Association of Mongolia, Ulaanbaatar, Mongolia

    Gombojav Bayasgalan & Danshiitsoodol Narandalai

59. Graduate School of Biomedical and Health Sciences, Hiroshima University, Hiroshima, Japan

    Danshiitsoodol Narandalai

60. Department of Psychology, University of Texas at Austin, Austin, TX, USA

    K. Paige Harden & Elliot M. Tucker-Drob

61. Department of Twin Research and Genetic Epidemiology, King’s College, London, UK

    Timothy Spector, Massimo Mangino & Genevieve Lachance

62. School of Law, Psychology and Social Work, Örebro University, Sweden

    Catherine Tuvblad

63. Washington State Twin Registry, Washington State University - Health Sciences Spokane, Spokane, WA, USA

    Glen E. Duncan

64. Washington State Twin Registry, Washington State University, Seattle, WA, USA

    Dedra Buchwald

65. Department of Clinical Research, University of Southern Denmark, Odense, Denmark

    Kirsten O. Kyvik

66. Odense Patient data Explorative Network (OPEN), Odense University Hospital, Odense, Denmark

    Kirsten O. Kyvik

67. Department of Clinical Biochemistry and Pharmacology and Department of Clinical Genetics, Odense University Hospital, Odense, Denmark

    Kaare Christensen

68. Department of Statistics, Faculty of Arts and Sciences, Kırıkkale University, Kırıkkale, Turkey

    Sevgi Y. Öncel

69. Departments of Psychiatry, Psychology and Human and Molecular Genetics, Virginia Institute for Psychiatric and Behavioral Genetics, Virginia Commonwealth University, Richmond, USA

    Fazil Aliev

70. Department of Epidemiology, School of Public Health, University of Washington, Seattle, WA, USA

    Jack H. Goldberg

71. Novo Nordisk Foundation Centre for Basic Metabolic Research (Section on Metabolic Genetics) and Department of Public Health, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark

    Thorkild I. A. Sørensen

72. Institute of Preventive Medicine, Bispebjerg and Frederiksberg Hospitals, The Capital Region, Copenhagen, Denmark

    Thorkild I. A. Sørensen

73. National Institute for Health and Welfare, Helsinki, Finland

    Jaakko Kaprio

74. Institute for Molecular Medicine FIMM, Helsinki, Finland

    Jaakko Kaprio

K.S., J.K., T.I.A.S., D.I.B., F.R., K.O.K., Y.-M.H. and Y.Y. planned the study design of the CODATwins project. A.B., C.K., K.J.S., K.L.J., W.C., A.E.H., T.M.M., W.G., C.Y., L.L., R.P.C., B.M.H., K.C., A.S., K.O.K., T.I.A.S., C.A.D., R.F.V., R.J.F.L., J.K., K.H., J.W., C.H.L., A.F., T.A.M., T.C.E., A.M.G., M.H., X.D., M.B.-A., H.B.-N., M.S., A.D.T., D.L.T., M.A.S., C.F., C.D., S.O., A.K.-N., D.M., L.A., S.A.B., K.L.K., H.H.M., L.J.E., J.L.S., R.F.K., M.M., S.P., M.G., D.A.B., D.I.B., G.W., T.C.E.M.v.B., M.B., C.H., J.M.C., R.S., D.L.F., J.A.M., F.J., F.N., Z.P., L.D., M.B., M.B., G.D., F.V., N.G.M., S.E.M., G.W.M., Y.-M.H., J.H.G., S.M., J.v.B.H., S.A., R.S., G.E.S., R.K., P.K.E.M., N.L.P., F.R., P.T., P.L., C.M.A.H., R.P., K.P.H., E.M.T.-D., S.Y.O., F.A., G.B., D.N., T.S., M.M., G.L., L.A.B., C.T., G.E.D., D.B. and Y.Y. collected the data used in this study. K.S. and A.J. were in charge of data management. A.J. conducted the analyses, wrote the first draft of the manuscript and has primary responsibility of the final content. All authors have commented the manuscript and read and approved the final version of the manuscript.

The authors declare no competing financial interests.

• BMC Medicine (2022)

• Nutrition Journal (2022)

• Scientific Reports (2022)

• Cliometrica (2022)

• Scientific Reports (2021)

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.