# Increasing gap in human height between rich and poor countries associated to their different intakes of N and P

Josep Penuelas, Ivan Janssens, Philippe Ciais, Michael Obersteiner, Tamás

Krisztin, Shilong Piao, Jordi Sardans

► To cite this version:

Josep Penuelas, Ivan Janssens, Philippe Ciais, Michael Obersteiner, Tamás Krisztin, et al. Increasing gap in human height between rich and poor countries associated to their different intakes of N and P. Scientific Reports, 2017, 7 (1), pp.17671. 10.1038/s41598-017-17880-3. hal-02903975

HAL Id: hal-02903975

https://hal.science/hal-02903975v1

Submitted on 17 Sep 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
SCIENTIFIC REPORTS

OPEN

Received: 2 May 2017

Accepted: 2 December 2017

Published online: 15 December 2017

# Increasing gap in human height between rich and poor countries associated to their different intakes of N and P

Josep Peñuelas<fnref n="1" />,², Ivan A. Janssens³, Philippe Ciais⁴, Michael Obersteiner⁵, Tamás Krisztin<fnref n="1" />,⁵, Shilong Piao⁶ & Jordi Sardans<fnref n="1" />,²

We analyzed mean height of men born in the 1960s, 1970s and 1980s in 80 countries. Both height and the change in height during the last decades were correlated with N and P intake, as well as the N:P intake ratio. Rich countries had higher per capita N and P intake than poor countries (on average $19.5 \pm 0.3$ versus $9.66 \pm 0.18$ kg N y⁻<fnref n="1" /> and $2.17 \pm 0.04$ versus $1.35 \pm 0.02$ kg P y⁻<fnref n="1" />), and also larger increases in per capita N intake ($12.1 \pm 2.0\%$ vs. $7.0 \pm 2.1\%$) and P intake ($7.6 \pm 1.0\%$ vs $6.01 \pm 0.7\%$), during the period 1961–2009. The increasing gap in height trends between rich and poor countries is associated with an increasing gap in nutrition, so a more varied diet with higher N, P, and N:P intake is a key factor to improve food intake quality in poor countries and thus shorten the gap with rich countries. More N and P are needed with the consequent requirements for a better management of the socioeconomic and environmental associated problems.

Height is a fundamental human trait. Within countries, taller people tend to have, on average, more education, higher earnings and possibly even higher social position<fnref n="1" />⁻³. They live longer and have a lower risk of cardiovascular and respiratory diseases, but a higher risk of some cancers (references in NCD-RISC, 2016⁴; Paajanen et al.⁵; World Cancer Research Foundation/American Institute for Cancer Research 2007⁶; Davies et al.⁷; Zhang et al.⁸).

Human height is a quantitative trait controlled by multiple genes⁹,¹⁰, but it also depends on environmental factors¹⁰⁻¹³. Among the latter, nutrition, particularly the intake of the main bio-elements, N and P, during childhood and adolescence, should be particularly important determinants of height during adulthood. Several studies have observed the link between infant growth and N-metabolism imbalances¹⁴,¹⁵ and with low P intake¹⁶,¹⁷. Moreover, the weight of infants at birth has been positively related to the P concentrations in blood of mothers during gestation¹⁸ and the growth rates of new-born infants have been linked to P concentration in mother's milk¹⁹. Improvements in growth rates in infants have been associated to higher capacity to retain N in their body²⁰.

The main source of N for humans comes from proteins and amino acids and at lesser extent from nucleic acids in food intakes²¹,²². Whereas in affluent countries the easy supply of N-fertilizers helps to produce an excess of food in general, and animal food in particular, in most underdeveloped countries the use of N-fertilizers makes the difference between adequate diet and malnutrition²². The main food sources of P are milk, chocolate and their derivatives, cereal grains, mollusks, beans and nuts, carbonate and cola beverages, fish and meat, and specially organ meats²³,²⁴. Most absorbed phosphorus comes from organic-P bound *in vivo* to animal and plant proteins although most plant phosphorus is mostly associated with phytates and is less absorbed in gut²³. Several industrial processed food products include food additives rich in P²⁵, which can even drive to an excess of P dietary intake that promotes higher serum parathyroid hormone (PTH) intake and lower serum Ca concentrations with negative impacts on bone health²⁵,²⁶.

<footnote n="1">¹CSIC, Global Ecology Unit CREAF-CEAB-UAB, Cerdanyola del Valles, 08193, Catalonia, Spain. ²CREAF, Cerdanyola del Valles, 08193, Catalonia, Spain. ³Research Group of Plant and Vegetation Ecology (PLECO), Department of Biology, University of Antwerp, B-2610, Wilrijk, Belgium. ⁴Laboratoire des Sciences du Climat et de l'Environnement, IPSL, 91191, Gif-sur-Yvette, France. ⁵International Institute for Applied Systems Analysis (IIASA), Ecosystems Services and Management, Schlossplatz 1, A-2361, Laxenburg, Austria. ⁶Sino-French Institute for Earth System Science, College of Urban and Environmental Sciences, Peking University, Beijing, 100871, China. Correspondence and requests for materials should be addressed to J.P. (email: josep.penuelas@uab.cat)</footnote>
In addition to the differences in total amount of food intake (calories), we hypothesized a greater gap in the per capita N and P food intake resulting from the higher per capita food intake from animal sources, richer in N, and the larger use of food additives, rich in P, in developed than in developing and undeveloped countries. We also hypothesized that this gap in the global N, P and the N:P intake ratio between developed and less developed countries is increasing and that it would go associated with an increasing gap in growth and adult height.

Here, we thus aimed to test the hypotheses that, across countries, i) there is a gap in the human height and the N, P and N:P intake ratio between developed and underdeveloped countries, ii) this gap is increasing, and iii) the variability in human height and in recent trends therein is associated with the N and P contents (and their ratio N:P) in food intake and trends therein, while controlling for other factors which might influence height such as calories intake, total GDP, % of low weight infants at birth (less than 2500 g), the increase in urban population which might have an influence on average human height through an increase in pollution, and the Human Development Index (HDI) as indicators of economic, health and development conditions. HDI is a summary measure of average achievement in key dimensions of human development: a long and healthy life, being knowledgeable and have a decent standard of living. The HDI is the geometric mean of normalized indices for each of the three dimensions. The health dimension is assessed by life expectancy at birth, the education dimension by mean of years of schooling and the standard of living dimension by gross national income per capita.

To test these hypotheses, we used human height and food intake databases from FAO²⁷, DTU²⁸, USDA²⁹, INFOODS³⁰ and Universität Tübingen³¹, and analyzed the relationships of adult male height with per capita intake of N and P, the N:P ratio, food origin (animal or vegetal) and their changes in recent decades. We also considered the possible role of GDP, % of low weight infants at birth³², percentage of urban population, and per capita calories intake using databases from World Bank and FAO. HDI was obtained from Human Development Programme³³.

We defined two groups of countries to separate those with more animal-derived N and P intake from those with more plant-derived N and P intake, mostly coinciding with those with more N and P intake and those with less N and P intake respectively (Fig. 1): a) countries with (animal N/plant N) + (animal P/plant P) > 2 and those with this ratio <2. Figure 1 shows the relationship between (animal P/plant P) per capita intake versus (animal N/plant N) per capita intake in different countries and human cohorts (1960s, 1970s and 1980s).

## Results

Averaged over entire populations, men born in richer, higher-income countries had higher intake of animal than of vegetal products, had higher N, P and N:P intake, and were taller than men born in poorer countries (Figs 2–5). Across the 80 countries included in our analysis, the tallest adult men were found in Denmark and the Netherlands, ca. 183 cm, and the shortest were found in Guatemala and Vietnam, ca. 160 cm (Fig. 2).

Mean adult men height across different countries and human cohorts was correlated with per capita N intake, per capita P intake, N:P intake and Calories intake (R² = 0.48, R² = 0.35, R² = 0.37 and R² = 0.45, respectively) (Fig. 3, Figure S1a). The relationships of adult men height with GDP, HDI and % of low weight infants at birth were not linear. Men height increased strongly with increasing GDP at low values of GDP of poor countries and tends to progressively increase less when GDP reach values above 20000 $ per capita in developed richer countries (Figure S1b). The relationship of adult men height with HDI was positive only after HDI values around 0.4 mostly in developed countries (Figure S1c). Mean height strongly decreased with % of low weight infants at birth at the low values of this variable characteristic of developed countries but the relationship disappeared at the higher values characteristic of poorer countries (Figure S1d).

When analyzed by multiple regressions, the relationships between the absolute human male height differences among countries of 1980s cohort and the differences in N intake and GDP (Fig. 4a), the model explained 55% of total human male height differences, and total N intake had a stronger significance in the model (F = 521) than GDP (F = 82). Similarly, when considering the relationships between the absolute human male height differences among countries of 1980s cohort and the corresponding differences in per capita N intake and calories intake (Fig. 4b) the model explained a 54% of height differences, and total per N intake differences had a stronger significance in the model (F = 794) than the differences in calories intake (F = 82). The best SEM model explaining the relationships of height differences with differences in all the studied variables explained a 65% of height differences among countries (Figure S2). The SEM analysis showed once more that the relationships of height differences among countries with per capita N intake differences were stronger than with socioeconomic and sanitary variables such as GDP, % of low weight infants at birth and HDI differences. Moreover, the positive relationships of these socioeconomic variables with height differences among countries were mainly due to indirect effects through the effects on per capita N, P and N:P intake differences (Figure S2).

Bayesian analyses of these data showed that the N:P ratio had the highest significant impact on human male height in poor and especially rich countries (Fig. 5). Rich countries also displayed a significant positive effect of GDP, whereas poor countries did not. The results of the Bayesian analyses considering also HDI and % of low weight infants at birth also reinforce that the most important factors for human height is the ratio of the per capita N and P intake (posterior inclusion probability of 97), as well as the per capita N intake (posterior inclusion probability of 81.6) (Table S5, Fig. 6).

Relative to men born in the 1960s, the largest increases in adult height in men born in the 1980s were found in Japan, ca. 6 cm, and Denmark and Portugal, ca. 3 cm (Fig. 7). In contrast, adult height decreased in some countries, mainly from the African continent, with the largest decreases in Togo, ca. -3 cm, followed by Madagascar, ca. -2 cm. Height changes differed greatly between rich and poor countries, with an average increase in height in rich countries of 1.5 cm in the cohort born in the 1980s relative to the cohort born in the 1960s, and with no increase in average height across the poor countries during the same period (Fig. 7).

Although less sensitive (lower slopes in the regression lines) than the absolute height, also the change in height between men born in the 1960s and in the 1980s was correlated with the change in N and P intake during the
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Relationships between (animal P/plant P) and (animal N/plant N) ratios per capita intake in different countries and human cohorts (1960s, 1970s and 1980s). The country abbreviations are: Armenia (Ar), Australia (A), Austria (Au), Azerbaijan (Az), Bangladesh (Ba), Belgium (Be), Benin (B), Bolivia (Bo), Botswana (Bt), Brazil (Bz), Burkina Faso (Bu), Cabo Verde (CV), Cambodia (Cm), Cameroon (Ca), Canada (CN), Central Africa Republic (Ct), Chad (CH), Chile (CI), China (C), Colombia (Co), Congo (Cg), Costa Rica (CR), Czechoslovakia (Cz), Denmark (DN), Djibouti (D), Dominican Republic (Do), Egypt (E), United Arab Emirates (Em), Estonia (Es), Ethiopia (Et), Finland (FN), France (FR), Gabon (Ga), Germany (GE), Ghana (Gh), Greece (GR), Guatemala (Gu), Guinea (Gi), Guyana (Gy), Haiti (Ha), Honduras (Ho), Hungary (H), India (In), Indonesia (Id), Iran (Ia), Ireland (IR), Italy (I), Jamaica (JA), Japan (Jp), Jordan (Jo), Kazahkstan (Ka), Kenya (Ke), North Korea (K), Kyrgyzstan (Ky), Lesotho (Le), Liberia (Li), Madagascar (Ma), Malawi (Mw), Mali (Ml), Mexico (Me), Moldavia (Mo), Morocco (Mr), Mozambique (Mz), Namibia (Na), Nepal (Ne), Netherlands (NT), Nicaragua (Ni), Niger (Ng), Nigeria (Nr), Norway (No), Panama (P), Peru (Pe), Phillipines (PH), Poland (PL), Portugal (Po), Russia (R), Rwanda (Rw), Saudi Arabia (As), Senegal (Se), Sierra Leone (Si), Slovenia (SL), South Africa (Sa), Spain (Sp), Sri Lanka (Sk), Swaziland (Sz), Sweden (SW), Taiwan (Tw), Tanzania (Ta), Togo (To), Trinidad Tobago (Tr), Tunisia (Tu), Turkey (Tk), Uganda (Ug) United Kingdom (UK), United States of America (US), Uzbekistan (Uz), Vietnam (VN), Yemen (Y), Zambia (Z), Zimbabwe (Zi). The countries with higher N and P intake from animal than vegetal products are indicated in red, and the countries with higher N and P intake from vegetal products are indicated in blue.</figcaption></figure>

same period (Fig. 7). Rich countries increased the intake of N (on average $2.3 \pm 0.4 \text{ kg y}^{-1}$, i.e. $12.1 \pm 2.0\%$) more than the intake of P (on average $0.17 \pm 0.03 \text{ kg y}^{-1}$, i.e. $7.6 \pm 1.0\%$), so the intake N:P ratio or imbalance between N and P has increased ($0.3 \pm 0.08$, i.e. $3.53 \pm 0.73\%$). Poor countries increased significantly less their intake of N ($0.7 \pm 0.2 \text{ kg y}^{-1}$, i.e. $7.0 \pm 2.1\%$) or P ($0.08 \pm 0.2 \text{ kg y}^{-1}$, i.e. $6.01 \pm 0.7\%$), and their N:P ratio did not change statistically significantly. These minor increases in nutrient intakes occurred in parallel with a lack of an overall increase in height, with many individual countries, mostly African, even exhibiting a decreasing trend in height (Fig. 7).

Differences in N:P, P and N intake had the largest posterior impact on differences in human height between cohorts (Figure S4). Based on the densities, we can reject the hypothesis, that – when controlling for N and P intake – the differences in calories or urban population have an impact that is significantly different from zero on the differences of human height. The results of the Bayesian analyses including also HDI and % of low weight infants at birth also reinforced the main role of N:P and N and P intake (Figure S4).

## Discussion

Average male height was highly variable among countries, with a difference of ca. 23 cm between the tallest and shortest country means. This difference was associated with differences in the per capita N, P and N:P intake among countries, more strongly than with other possible explanatory factors such as daily calories intake associated to energy for living but not directly to building blocks for growth.

In recent decades, adult height has changed substantially in many countries. The change has been substantial but less sensitive to the recent change in N and P intake than expected from the general relationships between height and N and P intake (lower slopes in the regression lines). This slower response to changes in N and P intake at short term indicates the existence of other scaling factors integrated at longer time scales in the relationship height-N and P intake. Genetic and epigenetic changes are involved at these longer term scales. Genetic variation
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Mean heights of the men born in the 1960s (blue), 1970s (red) and 1980s (green) by country. The countries with higher N and P intake from animal than vegetal products are indicated in red, and the countries with higher N and P intake from vegetal products are indicated in blue. See the caption for Fig. 1 for the country abbreviations.</figcaption></figure>

explains the reduced slope of the relationship between increased intake of nutrition factors and differences in height between time intervals: there is a limit to the phenotypic response to improved environmental conditions that a fixed genetic set up can produce. The increasing differences in immigration rates among countries, though, will generate different levels of gene flow. This gene flow can introduce new variation that can change the phenotypic response in the new cohorts.

Although height trends varied strongly among countries, the difference in adult male height between the rich and poor countries increased, on average, by ca. 1.5 cm in the last three decades. This increasing height difference was the result of men in richer countries becoming even taller (on average), whereas the men in poor countries did not increase in height (on average). This divergent height growth pattern was associated with a different trend in N and P intake between rich and poor countries, and specially with increasing N:P ratio of food intake, indicative of increasingly larger relative intake of food of animal origin in rich countries, and indicative of absence of saturation in height and in N and P intake. Food of animal origin has, on average, higher N and P contents and N:P ratios than food from plants, and the percentage of animal-derived food of total intake is higher in developed than less developed countries (Fig. 1). In addition to N and P intake and nutrition, other factors may impact on the increase of height as highlighted by the very high increase in Japan, much above the explained value by the overall relationships with N, P and N:P uptake (Figs 3 and 7).

Our results, both of current mean height and its changes in recent decades, are consistent with available data and trends previously recorded in many countries, usually treated as indicators of social and environmental conditions that affect nutrition, health and economic prosperity in each country and generation. For example, previous studies have documented an increase in stunting in men in some poor countries, such as those of sub-Saharan Africa⁴,³⁴. Rising population, coupled with worsening economic status, may underlie the decrease in height in these African countries³⁴⁻³⁷.

Height is not a neutral trait. In addition to behavioral and social links, higher adult height of men within the population of a same country or geopolitical region is known to be associated with lower occurrence of cardiovascular and respiratory diseases, but also with higher occurrence of colorectal and possibly pancreatic and prostate cancers, independently of its inverse correlation with body mass index⁶,⁷,³⁸. If these associations would be causal³⁹, then the ca. 23 cm height range globally could be associated with a 17% lower risk of cardiovascular mortality and 20-40% higher risk of various site-specific cancers in countries with tall people relative to countries with short people⁴⁰. Increases in mean population height in successive cohorts are associated with lower mortality in middle and older ages in countries with reliable mortality data⁴⁰, consistent with individual-level evidence of the association between taller height and lower mortality from all causes in adults⁴¹. Caution must, however, clearly drive the application of these associations when comparing among countries and geopolitical regions instead of individuals of the same population group.

The average height of young adult men in different countries is also an indicator of sustainable human development⁴⁰. Adult height reflects not only fetal and early childhood nutrition, which was included in the Millennium Development Goals, but also that of adolescents³⁴. Adult height is also associated with these early-life experiences and with longevity, education and earnings. It can easily be measured in health surveys and can be used to investigate differences across countries and trends over time, as we have, and to investigate within-country inequalities.

The association between height and nutrition, particularly with N and P intakes and the corresponding N:P ratio, highlights the importance of these two bio-elements. The increases in N and P intake in rich countries favor
<figure><img src="image_4.png" /></figure>

**Figure 3.** Relationships of mean male height with the average annual total N (a), P (b) and N:P ratio (c) intake per person (kg y⁻¹) in the corresponding country and cohort. See the caption for Fig. 1 for the country abbreviations. Men born in the 1960s, 1970s and 1980s are identified by 6, 7 and 8 in the symbols, respectively. The countries with higher N and P intake from animal than vegetal products are indicated in red (also red arrow for the mean), and the countries with higher N and P intake from vegetal products are indicated in blue (also the arrow for the mean). Different letters on the axes indicate significant differences (*P* < 0.05). For detailed numerical values for each country and cohort see Table S2 (Supporting information).

<figure><img src="image_5.png" /></figure>

**Figure 4.** Relationships of mean male height differences (men cohort of 1980s) between all pairwise comparisons of the 80 countries studied with the corresponding differences in annual N intake and GDP (a), and with the corresponding differences in annual N intake and daily calories intake (b).

rapid protein synthesis and maximizes the genetic capacity for attaining a maximal height. Several other factors in which rich and poor countries differ, such as the quality and accessibility of sanitary services, can indirectly influence the differences in height, but it is food that provides the main sources, the building blocks, to grow while other factors can directly influence more other variables such as life expectancy and health. More N and P intake is a direct consequence of the intake of more proteins, which are the molecules underlying growth processes at cellular scale. The relationship between socioeconomic variables and the gap of human height between rich and poor countries is mostly due to their indirect effects through total per capita N and P intake and N:P intake ratio.

The increase in height is thus associated to larger N and P intakes. For N-rich crops N fertilization may not be a problem, but most soils globally contain only low concentrations of available forms of P, so P inputs in the form of fertilizers are essential for food production around the globe. The majority of agricultural P is mined from finite and therefore exhaustible⁴²⁻⁴⁵ sedimentary deposits of phosphorite. The growing human and livestock populations requiring food and feed, and the growing demand for biofuels may exhaust these reserves in the next 40 to 400 years⁴²⁻⁴⁴. P is thus expected to become economically inaccessible to low-income and food-deficient countries, which will likely further widen the gap in nutrition and therefore in height and health between rich and
<figure><img src="image_6.png" /><figcaption>**Figure 5.** Posterior densities of the coefficient impacts for absolute human male height in the 1960s, 1970s and 1980s cohorts. The countries with higher N and P intake from animal than vegetal products are indicated in red (i), and the countries with higher N and P intake from vegetal products are indicated in blue (ii). The figures are ordered by the sum of the posterior inclusion probability of the displayed coefficients. A relative, high density mass at zero (signified by the vertical dashed line) corresponds to a high probability of the coefficient being excluded from the model.</figcaption></figure>

poor countries. The implications of this yield gap to food security add a key additional requirement of fundamental changes in agricultural and food management to those already posed by the need of mitigating climate change.

Our results show that increases in N and P intakes are more important than the increases in calories intake to increase height in a given country. Thus, a more varied diet with higher N and P intake is a key factor to improve food quality in poor countries and thus shorten the gap with rich countries. The range of calories intake between countries (from 1879 to 3721 daily Cal intake) is much narrower than the differences in N (from 3.3 to 23.7 per capita kg y⁻¹) and P intake (from 0.41 to 2.76 per capita kg y⁻¹), meaning that not only calories intake must be increased but mainly N and P intake and food variability (more proteins and vitamins). This requires increasing crops with more N and P content, and/or animal derived food production, which by itself has severe sustainable drawbacks⁴⁶. Consequently more sources of N and P will be needed to improve the food quality and variability of poor countries with the consequent socioeconomic and environmental problems to be solved.
<figure><img src="image_7.png" /><figcaption>**Figure 6.** Posterior densities of the coefficient impacts for the differences between countries in the 1980s cohort, ordered by posterior inclusion probability. A high density mass at zero (signified by the vertical dashed line) corresponds to a high probability of the coefficient being excluded from the model. The bold, continuous colored lines denote the median posterior impact of the coefficient, conditional on its posterior inclusion probability. The dashed colored lines depict +/− two posterior standard deviations, conditional on the coefficient being included in the model.</figcaption></figure>

## Materials and Methods

We estimated the mean height of young adults born in the 1960s, 1970s and 1980s (i.e. people who had reached their 21st birthday between 1981 and 2010) in 80 countries. We used the data from University of Tubingen-World Height databases (www.uni-tuebingen.de) that provided the average adult male heights for each of these three decades. This database contains all data available from official sources including WHO/UNESCO, U.S. Department of Health and Eurobarometer, among others. We organized countries based on per capita N and P intakes of animal and plant origin (Supplementary Table 1).

We calculated the annual food intake of N and P per person for the three cohorts with data of adult height for each country as:

$$ (\Sigma \text{ annual intake of each food group[1]} \times \text{mean N or P concentration of each food group[2]})/\text{number of years} $$

* Data from FAOSTAT (2014)

* Data from INFOODS_FOOD_DATABASE FOR BIODIVERSITY (2016), USDA (2016) and DTU Fodeva-reinstituttet (2016).
<figure><img src="image_8.png" /><figcaption>**Figure 7.** Relationships of the change in male height (in cm) from the 1960s cohort to the 1980s cohort in the countries with the increase in annual N (a), P (b) and N:P ratio (c) intake during the growth period (1960s-1980s for the 1960s cohort, and 1980s-2000s for the 1980s cohort). Red arrows represent the countries with higher N and P intake from animal products, and blue arrows represent the countries with higher N and P intake from vegetal products. Different letters on the axes indicate significant differences (P < 0.05). See the caption for Fig. 1 for the country abbreviations.</figcaption></figure>

To estimate the N and P concentrations of each food group of FAO databases we used the data bases [2]. In these data bases there are the N and P concentrations of the different food items. We grouped these different food items in different sets corresponding to the distinct FAO food groups, and made the corresponding average of each of these groups. We used the average as a final value when data for N and/or P concentrations were provided by more than one database for the same food group.

To obtain the mean annual intake of N and P per person in each country during the period of growth of the three human cohorts with data for male height, we used the average from 1961 to 1989 for the cohort born in the 1960s, from 1971 to 1999 for the cohort born in the 1970s and from 1981 to 2009N for the cohort born in the 1980s. The N:P ratios were estimated on a mass basis.

The increases in annual P and N intake and adult male height for each country in men born in the 1980s relative to those born in the 1960s were estimated for all countries and cohorts for which information was available.

We obtained the data of GDP and percentage of urban population from World Bank⁴⁷. To obtain data of daily per capita calories intake we used data from FAO (http://www.fao.org/economics/ess/food-security-statistics/food-security-statistics-metadata/en).

We used reduced major axis analysis to assess the relationships of height and height change with N, P and N:P ratio intake for all countries and territories. We used general linear models to analyze the relationship of men height differences between all the pairwise comparisons among the studied countries with the corresponding differences in per capita calories intake, GDP and per capita N and P annual intake. We also conducted general linear models to analyze the relationships between the country men height with the corresponding country GDP, and per capita daily calories, N, P and N:P intakes. We similarly analyzed the relationships of the differences in men adult height between the cohort of the 1980s and the cohort of the 1960s with the changes in GDP, per capita calories, N, P and N:P intake during their respective growth periods.

Due to the limited number of observations and the lack of a strictly theoretical framework, we also used a Bayesian model averaging⁴⁸ approach to explore a large number of potential models and draw inference over the relative importance of individual covariates. Moreover, the resulting parameter estimates stem from averages over the models with the highest posterior probability, thus providing robust estimates.

We modeled human height in absolute levels and across the cohorts with a Bayesian model averaging approach. For this purpose we had 208 observations, ranging over 86 countries. The average male human height in a country *i* in cohort *t* was modeled as:

$$ height_{it} = \sum_{k=1}^{2} [\beta_1 + \alpha_t + \beta_2 N_{it} + \beta_3 P_{it} + \beta_4 GDP_{it} + \beta_5 GDP_{g_{it}} \\ + \beta_6 Cal_{it} + \beta_7 N_{it} \\ + \beta_8 P_{it} + \beta_9 \frac{N}{P_{it}} \\ + \beta_{10} P_{animal} / P_{vegetables_{it}} + \beta_{11} N_{animal} / N_{vegetables_{it}} \\ + \beta_{12} Urban_{it}] I_{(\omega_{itk}=1)} + \varepsilon_{it} \quad (1) $$
where height<sub>it</sub> denotes the height of the t-th cohort in country i. As in the previous equation, β<sub>1</sub> is the intercept, α<sub>1</sub> a period specific dummy, N<sub>it</sub> and P<sub>it</sub> are the levels of N and P intake per capita in cohort t GDP<sub>it</sub>, GDP<sub>it</sub>, Cal<sub>it</sub>, and Urban<sub>it</sub> denote the level of GDP, GDP growth, calories and percentage of urban population, respectively. The ratios of N and P intake, as well as N and P intake from animals in proportion to vegetables are denoted by P<sub>animal</sub>/P<sub>vegetables</sub>, and N<sub>animal</sub>/N<sub>vegetables</sub>, respectively. β<sub>2</sub> to β<sub>12</sub> are the corresponding coefficients, and ε<sub>it</sub> is the i.i.d. normally distributed error term with zero mean and σ<sup>2</sup> variance. I<sub>(c)</sub> denotes an indicator function, which takes on the value of one if the condition ω<sub>itk</sub> = 1 is fulfilled and is zero otherwise. ω<sub>itl</sub> takes on the value of one if the sum of the N intake from animals per the N intake from vegetables plus the P intake from animals per the P intake from vegetables exceeds two and is zero otherwise. In this fashion we can model the impacts of both groups in one joint model, with a joint variance function.

We had complete data for the three cohorts for 54 countries in total. Taking the differences between the cohorts, results in n = 108 observations. We analysed these data with a similar Bayesian model where i denotes a specific country (with i = 1, ..., n), and t the specific cohort:

$$
\begin{align}
\Delta(\text{height})_{it} ={}& \beta_1 + \alpha_t + \beta_2 N_{it} + \beta_3 P_{it} + \beta_4 \Delta(\text{GDP})_{it} \nonumber \\
& + \beta_5 \Delta(\text{GDPg})_{it} + \beta_6 \Delta(\text{Cal})_{it} \nonumber \\
& + \beta_7 \Delta(N)_{it} + \beta_8 \Delta(P)_{it} + \beta_9 \Delta\left(\frac{N}{P}\right)_{it} \nonumber \\
& + \beta_{10} \Delta\left(\frac{P_{\text{animal}}}{P_{\text{vegetables}}}\right)_{it} + \beta_{11} \Delta\left(\frac{N_{\text{animal}}}{N_{\text{vegetables}}}\right)_{it} \nonumber \\
& + \beta_{12} \Delta(\text{Urban})_{it} + \varepsilon_{it} \tag{2}
\end{align}
$$

where the Δ(·) denotes the differences between cohorts t and t + 1. Δ (height)$_{it}$ denotes the difference in height between cohorts. β₁ is the intercept, α₁ a period specific dummy, Nᵢₜ and Pᵢₜ are the levels of N and P intake per capita in cohort t Δ(GDP)ᵢₜ, Δ(GDPg)ᵢₜ, Δ(Cal)ᵢₜ, Δ(N)ᵢₜ, Δ(P)ᵢₜ and Δ(Urban)ᵢₜ denote the difference in GDP, GDP growth, calories, N intake, P intake and percentage of urban population, respectively. The differences in ratios of N and P intake, as well as N and P intake from animals in proportion to vegetables are denoted by Δ(P<sub>animal</sub>/P<sub>vegetables</sub>)ᵢₜ, and Δ(N<sub>animal</sub>/N<sub>vegetables</sub>)ᵢₜ, respectively. β₂ to β₁₂ are the corresponding coefficients, and εᵢₜ is the i.i.d. normally distributed error term with zero mean and σ² variance.

We finally also conducted these same Bayesian analyses for those countries of the data base for which we had these same variables and additional variables HDI and % of low weight infants at birth.

References

1. Adair, L. S. et al. Associations of linear growth and relative weight gain during early life with adult health and human capital in countries of low and middle income: Findings from five birth cohort studies. *Lancet* **382**, 525–534 (2013).

2. Stulp, G., Buunk, A. P., Verhulst, S. & Pollet, T. V. Human height Is positively related to interpersonal dominance in dyadic interactions. *PLoS One* **10**, 1–18 (2015).

3. Case, A. & Paxson, C. Stature and status: height, ability, and labor market outcomes. *J. Polit. Econ.* **116**, 499–532 (2008).

4. NCD-Risk. Trends in adult body-mass index in 200 countries from 1975 to 2014: A pooled analysis of 1698 population-based measurement studies with 19.2 million participants. *Lancet* **387**, 1377–1396 (2016).

5. Paajanen, T. A., Oksala, N. K. J., Kuukasjärvi, P. & Karhunen, P. J. Short stature is associated with coronary heart disease: A systematic review of the literature and a meta-analysis. *Eur. Heart J* **31**, 1802–1809 (2010).

6. World Cancer Research Fund/American Institute for Cancer Research. *Food, Nutrition, Physical Activity, and the Prevention of Cancer: a Global Perspective*. (WorldCancer Research Fund/American Institute for Cancer Research,) at www.aicr.org/assets/docs/pdf/reports/Second_Expert_Report.pdf (2007)

7. Davies, N. M. et al. The effects of height and BMI on prostate cancer incidence and mortality: a Mendelian randomization study in 20,848 cases and 20,214 controls from the PRACTICAL consortium. *Cancer Causes Control* **26**, 1603–1616 (2015).

8. Zhang, B. et al. Height and Breast Cancer Risk: Evidence From Prospective Studies and Mendelian Randomization. *J Natl Cancer Inst* **107**, 1–17 (2015).

9. Fisher, B. A. The correlation between relatives on the supposition of Mendelian inheritance. *Trans. R. Soc. Edinb.* **399**–433 (1918).

10. Lettre, G. Recent progress in the study of the genetics of height. *Hum. Genet.* **129**, 465–472 (2011).

11. Deaton, A. Height, health, and development. *Proc. Natl. Acad. Sci.* **104**, 13232–13237 (2007).

12. Silventoinen, K., Kaprio, J., Lahelma, E. & Koskenvuo, M. Relative Effect of Genetic and Environmental Factors on Body Height: Differences Across Birth Cohorts Among Finnish Men and Women. *Am. J. Public Health* **90**, 627–630 (2000).

13. Haeffner, L. S., Barbieri, M. A., Rona, R. J., Bettiol, H. & Silva, A. A. The relative strength of weight and length at birth in contrast to social factors as determinants of height at 18 years in Brazil. *Ann. Hum. Biol.* **29**, 627–640 (2002).

14. Wu, G. et al. Arginine metabolism and nutrition in growth, health and disease. *Amino Acids* **37**, 153–168 (2009).

15. Rönnholm, K. A. R., Perheentupa, J. & Siimes, M. A. Supplementation with human milk protein improves growth of small premature infants fed human milk. *Pediatrics* **77**, 654–658 (1986).

16. Rowe, J. & Carey, D. Phosphorus deficiency syndrome in very low birth weight infants. *Pediatr. Clin. North Am.* **34**, 997–1017 (1987).

17. Loui, A., Raab, A., Obladen, M. & Brätter, P. Calcium, phosphorus and magnesium balance: FM 85 fortification of human milk does not meet mineral needs of extremely low birthweight infants. *Eur. J. Clin. Nutr.* **56**, 228–235 (2002).

18. Mataloun, M. M. & Leone, C. R. Human milk mineral intake and serum concentrations of calcium and phosphorus in newborn term infants: influence of intrauterine growth restriction. *Acta Paediatr.* **89**, 1093–7 (2000).

19. Khorana, M. & Jiamsajjamongkhon, C. Pilot study on growth parameters and nutritional biochemical markers in very low birth weight preterm infants fed human milk fortified with either human milk fortifier or post discharge formula. *J. Med. Assoc. Thail.* **97**, S164–S175 (2014).

20. Dempsher, D. P. et al. Whole body nitrogen kinetics and their relationship to growth in short children treated with recombinant human growth hormone. *Pediatr. Res.* **28**, 394–400 (1990).

21. Tomé, D. & Bos, C. Dietary Protein and Nitrogen Utilization. *J. Nutr.* **130**, 1868–1873 (2000).

22. Smil, V. Nitrogen and Food Production: Proteins for Human Diets. *AMBIO A J. Hum. Environ.* **31**, 126–131 (2002).

23. Noori, N. et al. Organic and inorganic dietary phosphorus and its management in chronic kidney disease. *Iran. J. Kidney Dis.* **4**, 89–100 (2010).
24. D'Alessandro, C., Piccoli, G. B. & Cupisti, A. The 'phosphorus pyramid': a visual tool for dietary phosphate management in dialysis and CKD patients. *BMC Nephrol.* **16**, 1–6 (2015).

25. Gutiérrez, O. Sodium and phosphorus-based food additives: persistent but surmountable hurdles in the management of nutrition in chronic kidney disease. *Adv. Chronic Kidney Dis.* **20**, 150–156 (2013).

26. Takeda, E., Yamamoto, H., Yamanaka-Okumura, H. & Taketani, Y. Increasing Dietary Phosphorus Intake from Food Additives: Potential for Negative Impact. *Adv. Nutr.* - *Am. Soc. Nutr.* **5**, 92–97 (2014).

27. Food and Agriculture organization of the United Nations. Food and Agriculture Organization of the United Nations. At http://faostat.fao.org/default.aspx (2014).

28. DTU - Technical University of Denmark. Fodevareinstitutet (http://www.food.dtu.dk/english). 21-22 (2016).

29. USDA-US Department of Agriculture. Foods List (https://ndb.nal.usda.gov/ndb/search/list). 1-12 (2016).

30. INFOODS. FAO/INFOODS Food Composition Database for Biodiversity. (2013).

31. Tübingen, U. (http://www.uni-tuebingen.de/en/faculties/economics-and-social-sciences/subjects/school-of-business-and-economics/faculty/economics/wirtschaftsgeschichte/data-hub-height.html). 1-4 (2016).

32. UNICEF. UNICEF Data: Monitoring the Situation of Children and Women TOPIC + COUNTRY. Unicef At, https://data.unicef.org/topic/nutrition/low-birthweight/ (2014).

33. Human Development Programme. Human Development Report 2016 Human Development for Everyone (2016).

34. Stevens, G. A. et al. Trends in mild, moderate, and severe stunting and underweight, and progress towards MDG 1 in 141 developing countries: A systematic analysis of population representative data. *Lancet* **380**, 824–834 (2012).

35. Weil, D. E. C., Alic busan, A. P., Wilson, J. F., Reich, M. R. & Bradley, D. J. *The impact of development policies on health A review of the literature*. (1990).

36. Pongou, R., Salomon, J. A. & Ezzati, M. Health impacts of macroeconomic crises and policies: Determinants of variation in childhood malnutrition trends in Cameroon. *Int. J. Epidemiol.* **35**, 648–656 (2006).

37. Sundberg, S. Agriculture, poverty and growth in Africa. *CAB Rev. Perspect. Agric. Vet. Sci. Nutr. Nat. Resour* **4**, 11–12 (2009).

38. Green, J. et al. Height and cancer incidence in the Million Women Study: prospective cohort, and meta-analysis of prospective studies of height and total cancer risk. *Lancet Oncol.* **12**, 785–794 (2011).

39. Nüesch, E. et al. Adult height, coronary heart disease and stroke: a multi-locus Mendelian randomization meta-analysis. *Int. J. Epidemiol.* dyv074, https://doi.org/10.1093/ije/dyv074 (2015).

40. NCD-RISC. A century of trends in adult human height. *Elife* **5**, 1–29 (2016).

41. Emerging Risk Factors. Adult height and the risk of cause-specific death and vascular morbidity in 1 million people: Individual participant meta-analysis. *Int. J. Epidemiol.* **41**, 1419–1433 (2012).

42. Cordell, D. & White, S. Peak phosphorus: Clarifying the key issues of a vigorous debate about long-term phosphorus security. *Sustainability* **3**, 2027–2049 (2011).

43. Elser, J. & Bennett, E. A broken biogeochemical cycle. *Nature* **478**, 29–31 (2011).

44. Obersteiner, M., Peñuelas, J., Ciais, P., van der Velde, M. & Janssens, I. A. The Phosphorus Trilemma. *Nat. Geosci.* **6**, 897–898 (2013).

45. Peñuelas, J. et al. Human-induced nitrogen-phosphorus imbalances alter natural and managed ecosystems across the globe. *Nat. Commun.* **4**, 2934 (2013).

46. Tilman, D. & Clark, M. Global diets link environmental sustainability and human health. *Nature* **515**, 518–522 (2014).

47. World Development Indicators Database. Urban population data. At http://data.worldbank.org/indicator/NY.GDP.PCAP.CD (2016).

48. Raftery, A. E., Madigan, D. & Hoeting, J. Bayesian model averaging for linear regression models. *J. Am. Stat. Assoc.* **92**, 179–191 (1997).

## Acknowledgements

The authors would like to acknowledge the financial support from the European Research Council Synergy grant ERC-SyG-2013-610028 IMBALANCE-P, the Spanish Government grant CGL2016-79835-P and the Catalan Government grant SGR 2014-274.

## Author Contributions

J.P. and J.S. designed and drafted the article and were in charge of the analyses together with T.K. and M.O. All the other authors (I.A.J., P.C., M.O., T.K., S.P.) discussed and improved the research and the article.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-017-17880-3.

**Competing Interests:** The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_9.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2017