Associations between High Protein Intake, Linear Growth, and Stunting in Children and Adolescents: A Cross-Sectional Study

Background and aims: Childhood and adolescence are critical periods for linear growth and preventing stunting. Current evidence indicates that dietary protein intake in children and adolescents is often two to three times higher than the recommendations in many regions worldwide. However, few studies have focused on the association between high protein intake and linear growth and stunting in this population. We aim to investigate this association in children and adolescents aged 6 to 18 years in a population with relatively high protein consumption. Methods: We conducted a large cross-sectional study involving 3299 participants from Shenzhen, a modern metropolis of China. Protein intake, including total protein, animal protein, and plant protein, was evaluated by a food-frequency questionnaire and expressed as grams per kilogram of body weight per day (g·kg−1·d−1) and as a percentage of total energy intake (%E). The primary outcomes were body height and height-for-age Z score (HAZ). Generalized linear models and logistic regression analyses were employed to examine the associations between protein intake and outcomes. We also conducted stratified analyses across different genders and pubertal stages in the aforementioned associations. Results: The mean protein intake was 1.81 g·kg−1·d−1 (17% E). After adjusting for serum calcium, zinc, vitamin D3, vitamin A levels, birth outcomes, lifestyle, and parental characteristics, each standard deviation increase of 1 in protein intake (0.64 kg−1·d−1) is found to be associated with a −5.78 cm change in body height (95% CI: −6.12, −5.45) and a −0.79 change in HAZ (95% CI: −0.84, −0.74). Consistent results were observed when protein intake was expressed as %E or specifically as animal or plant protein. Moreover, the relationship between protein intake and linear growth remained consistent across genders in different pubertal stages, similar to that of the overall participants. Conclusions: Our findings highlight the potential hazards of high protein intake on linear growth in children and adolescents. Caution should be exercised when promoting increased protein consumption in children and adolescents who already have a high intake of protein.


Introduction
Childhood and adolescence are critical periods for linear growth, which refers to the increase in body height.According to the Child Growth Standards established by the World Health Organization (WHO), individuals whose height falls below two standard deviations of the expected height for their gender and age are considered to have stunting [1].This condition affects approximately one-quarter of children worldwide [2].Stunted children and adolescents face a higher risk of mortality, morbidity, and suboptimal cognitive and motor development [3].Consequently, stunting is recognized as one of the six global nutrition targets endorsed by the WHO to improve child nutrition by 2025 [4].Therefore, it is crucial to maximize the potential for linear growth and reduce the occurrence of stunting in children and adolescents.
It is widely recognized that nutritional status directly impacts the growth and development of children and adolescents.Protein, which is an essential macronutrient, plays a crucial role in linear growth by meeting the metabolic demand for amino acids needed for tissue growth.Additionally, it elevates hormone levels such as insulin and insulinlike growth factor 1 (IGF-1), which promotes linear growth by stimulating endochondral ossification [5].Protein deficiency has been proven to lead to significant decreases in IGF-1, both in animal [6] and human [7] studies, and this growth-limiting effect can be neutralized by increasing dietary protein intake [6].Therefore, there is clear evidence that protein deficiency hampers linear growth.Interestingly, current evidence suggests that the protein intake levels in most developed countries exceed their dietary protein recommendations [8,9].This trend is even more pronounced among children and adolescents, with their protein intakes being 2-to 3-fold higher than the recommended levels in many developed countries [10].In China, for instance, 52.4% of children and adolescents have daily protein intakes that meet or exceed the dietary protein recommendations.Moreover, the ratios are significantly higher in urban areas compared to rural areas [11].All of the evidence indicate that a considerable number of children and adolescents consume protein beyond the recommended levels worldwide.However, there is limited research focusing on the association between high protein intake and linear growth within this population, particularly considering the prevailing trend of high dietary protein consumption in many developed countries and many cities of China.
In this context, we aimed to explore the association between protein intake and linear growth in a population with relatively high protein intake in a modern metropolis on the southeast coast of China.Considering the critical impact of gender and puberty stage on the correlation between protein intake and linear growth, we also examined this association in different pubertal stages across genders within the same population.

Study Participants
The study participants were from the Evaluation and Monitoring on School-based Nutrition and Growth in Shenzhen (EMSNGS).The EMSNGS study is an ongoing cohort study conducted in Shenzhen, a modern metropolis on the southeast coast of China.Its objective is to assess and monitor the nutritional status, growth and development, and other health outcomes in children and adolescents.The baseline survey for this study was completed in December 2021.
To select the participants, a multistage stratified cluster random sampling method was employed.Primary and secondary schools in ten administrative districts of the city were sampled.The sampling process consisted of two stages.In the first stage, one primary school, one junior high school, and one senior high school were chosen from each district, resulting in a total of 27 schools (including three nine-year schools).In the second stage, one class from each grade (excluding classes with fewer than 40 students) was randomly selected from each school.Overall, 112 classes were included, with a final participant count of 5153 school students.
All study participants and their guardians were informed about the study's protocol and provided informed consent before enrollment.The study received approval from the ethics review committee of the Shenzhen Center for Chronic Disease Control (SZCCC-2021-037-01-PJ) and was registered at the China Clinical Trials Registry (ChiCTR2100051722; www.chictr.org.cn,accessed on 1 January 2022).
For the analysis, certain participants were excluded based on specific criteria (see Supplementary Figure S1).Initially, 35 participants without basic information, 129 participants with diseases that could affect growth, 44 participants without available disease data, and 15 participants taking medications that might impact growth were excluded.Additionally, we excluded 1361 individuals without dietary data (food frequency questionnaires were not conducted for students in grades 1 through 3, totaling 1359) and those with implausible total energy intake (total daily energy intake less than 0.5 times or greater than 2 times the estimated energy requirements of Chinese residents of different ages, totaling 173), non-fasting participants (totaling 84), participants without height data (totaling 2), and participants aged below 6 or above 18 years old (totaling 13).Ultimately, the analysis included 3299 students.

Dietary Assessment
Dietary intake was collected using a validated Food Frequency Questionnaire (FFQ) consisting of 61 food items, carefully considering the eating habits and food composition specific to Chinese individuals [12].Experienced nutritionists provided guidance to participants, who estimated their food intake using food models representing standard portion sizes and food pictures depicting different portion sizes for all food items.Participants were asked to report their usual frequency of consuming each food item and the average serving size over the past month.The calculation of total energy intake and nutrient values relied primarily on the China Food Composition Database (6th edition), which served as the foundation for the food composition database [13].To account for individual variations in energy intake, protein intake and other dietary nutrients were adjusted using the nutrient residual method [14].Adjusting for total energy intake has been demonstrated to mitigate the impact of measurement errors inherent in FFQs when assessing protein intake.Furthermore, in this study, dietary protein intake was calculated per kilogram of body weight by incorporating energy-adjusted protein intake and body weights.This methodology has been previously validated in other studies [15].As a result, protein intake was expressed as grams per kilogram of body weight per day (g•kg −1 •d −1 ), as well as a percentage of total energy intake (%E).

Outcome Assessment
The primary outcomes included body height, height-for-age Z score (HAZ), and the prevalence of stunting.Qualified doctors and nurses who underwent rigorous training conducted the measurements of each participant's height and weight.Participants were required to wear loose clothing, be barefoot, feel comfortable, and extend both legs during the measurements.Height was recorded to the nearest 0.1 cm, and weight to the nearest 0.1 kg.HAZ was calculated using the gender-and age-specific growth standards provided by WHO in 2006 [16].According to the WHO Child Growth Standards, stunting is defined as a HAZ value below −2 [1].

Covariates
Information regarding age, gender, medical history, medication usage, lifestyle, family income, parental height, and education level was collected through structured questionnaires.Central obesity was defined as a waist circumference equal to or greater than the 90th percentile specific to age and gender [17].A pediatrician from the endocrinology department assessed pubertal development stages based on Tanner stages (Tanner stage I-V), categorizing them as pre-puberty (Tanner stage I), mid-puberty (Tanner stage II-III), and post-puberty (Tanner stage IV or higher) [18][19][20].Parents' educational levels were classified as ≤9, 10-12, 13-15, and ≥16 years of completed schooling.Household income levels were categorized as <120,000, 120,000-249,999, and ≥250,000 Chinese yuan (CNY) (1 CNY ≈ 0.13 EUR; 1 CNY ≈ $ 0.14 USD), based on the average yearly household income.Nighttime sleep duration was determined by calculating the time spent sleeping at night until waking up in the morning.Moderate to vigorous physical activity (MVPA) was classified as <1.0, 1.0-2.9, and ≥3 h per day.Smoking status was categorized as never or yes, including <1 cigarette per month, ≥1 cigarette per month, and e-cigarette use.Drinking status was categorized as never, <1 standard drink per month, and ≥1 standard drink per month.

Statistical Analysis
All covariates were assessed for normality and analyzed using nonparametric methods if necessary.Descriptive statistics at baseline included means ± standard deviations (SDs) for continuous variables and percentages for categorical variables.Associations between participant characteristics, protein intake, and outcome variables were examined using the χ 2 test for categorical variables and analysis of variance, and the Kruskal-Wallis rank test for continuous variables.Spline smoothing [21] was used to visually depict the dose-response relationship between protein intake and linear growth, as well as the rate of stunting.
Generalized linear models were employed to investigate the association between protein intake and height or height-for-age Z score (HAZ).Results were presented as coefficients (β) with 95% confidence intervals (CIs).Logistic regression analyses were conducted to calculate odds ratios (ORs) and 95% CIs, assessing the relationship between the prevalence of stunting and protein intake.Protein intake was evaluated as either a continuous variable (per 1 SD increase) or a categorical variable (quartiles of protein intake).Quartile 1 (Q1) was chosen as the reference group in each model when examining protein intake by quartiles.The linear trend across quartiles of protein intake was tested by assigning the median value for each quartile and analyzing it as a continuous variable in multivariate models.Protein intake was expressed in grams per kilogram weight per day (g•kg −1 •d −1 ), as well as a percentage of total energy intake (%E).
Covariates associated with protein intake and linear growth, as reported in previous studies, were included as potential confounders in the analyses.In the association of characteristics of study participants with protein intake and outcomes, confounders such as serum vitamin D 3 , serum Zn, serum vitamin A, dietary fat intake, dietary carbohydrate intake, and mothers' age were grouped according to tertiles.The following regression models were used: crude model (unadjusted); Model I, which was adjusted for gender, ethnicity, age (as a continuous variable), pubertal stages, premature, birth length, maternal age (as a continuous variable), parental education levels (as a continuous variable), parental height (as a continuous variable), household income levels, smoking status, drinking status, moderate to vigorous physical activity (MVPA), central obesity, length of nighttime sleep (as a continuous variable), dietary fat intake (as a continuous variable), and dietary carbohydrate intake (as a continuous variable); and Model II with additional adjustment for serum calcium, zinc, vitamin D 3 , and vitamin A levels, all as continuous variables.The premature status was missing for some of the participants (2/3299).We included the two participants with missing information on premature status in a separate missing data group in our main analysis.Similar procedures were applied to other missing variables, including household income, MVPA, length of nighttime sleep, fathers' height, fathers' education, and mothers' education.
To assess potential modifiers of linear growth, such as gender (male or female) and pubertal stages (pre-puberty, mid-puberty, post-puberty), stratified analyses were conducted, and p values for interaction were estimated using the log-likelihood ratio test.Analyses were also performed separately for animal and plant protein intake in relation to body height and HAZ.
All statistical analyses were conducted using Empower (R) software based on R (version 3.4.3)from X&Y Solutions, Inc., Boston, MA, USA.Significance was determined at two-tailed p values of < 0.05 for all analyses.

Characteristics of the Participants
A total of 3299 participants were included in the study, and Table 1 presents their characteristics.The mean protein intake was 1.91 ± 0.64 g•kg −1 •d −1 (17% ± 6%).The prevalence of stunting was 1.58% (52/3299).Protein intake, assessed in grams per kilogram weight per day or as a percentage of total energy intake, was significantly associated with gender, age, drinking status, MVPA, central obesity, length of nighttime sleep, serum zinc, vitamin D 3 , and vitamin A levels, as well as fat and carbohydrate intake (all p < 0.05).When assessing protein intake in grams per kilogram weight per day, additional associations were observed with pubertal stage, serum calcium level, parental height, and parental educational levels (all p < 0.05).Boys had higher height (p < 0.001) and HAZ (p < 0.001) than girls, resulting in a higher prevalence of stunting among females.Taller participants were older and had lower protein intake, more time spent on MVPA, greater incidence of central obesity, longer length of nighttime sleep, higher serum calcium, zinc, vitamin D 3 , and vitamin A levels, taller parents, and higher parental educational levels.HAZ was associated with protein intake, age, household income levels, MVPA time, central obesity, length of nighttime sleep, vitamin D 3 level, serum calcium level, parental height, and parental educational levels.The prevalence of stunting was associated with gender, age, ethnicity, paternal height, MVPA time, central obesity, protein intake, and parental educational levels.More details are provided in Table 1.

Association of Protein Intake with Linear Growth in Children and Adolescents
Figure 1 illustrates the dose-response relationship between protein intake and linear growth of the participants.In the fully adjusted Model II, a significantly negative association was observed between protein intake (g•kg −1 •d −1 ) and body height (Figure 1A) (p < 0.001), as well as HAZ (Figure 1B) (p < 0.001).A similar relationship was found between protein intake (%E) and body height (Figure 1C) (p < 0.001) and HAZ (Figure 1D) (p < 0.001).Results from the generalized linear model (Table 2) indicated that each unit increase in protein intake (g•kg −1 •d −1 ) was associated with a change in body height by −9.45 cm (95% CI: −9.72, −9.18) and a change in HAZ by −0.21 (95% CI: −0.25, −0.18) in the crude model.After adjusting for participant demographics, lifestyle factors, and paternal factors in Model I, each unit increase in protein intake was found to be associated with a change in height by −5.79 cm (95% CI: −6.12, −5.45) and a change in HAZ by −0.78 (95% CI: −0.83, −0.73).Further adjustment for serum calcium, zinc, vitamin D 3 , and vitamin A levels in Model I yielded consistent results, with each unit increase in protein intake (g•kg −1 •d −1 ) linked to a change in linear growth by −5.78 cm (95% CI: −6.12, −5.45) and a change in HAZ by −0.79 (95% CI: −0.84, −0.74).Similar results were obtained when assessing protein intake as a percentage of total energy intake (%E) (Table 2).We also separately examined animal and plant protein intake, both of which showed negative associations with linear growth, similar to the overall negative relationship with total protein intake (Supplementary Figure S2 and Supplementary Table S1).
Considering the possible impact of gender and puberty stage on the correlation between protein intake and linear growth, we conducted further analysis on the association between dietary protein intake and linear growth in different pubertal stages across different genders and found a similar significant negative correlation.More details are provided in Table 3.
In addition, to assess potential modifiers of linear growth, such as gender (boy or girl) and pubertal stages (pre-puberty, mid-puberty, post-puberty), stratified analyses were conducted.The relationship between protein intake, linear growth, and HAZ remained consistent in stratified analyses by gender and pubertal stages.However, we observed that the association between dietary protein intake and height was more pronounced in males (p for interaction < 0.001) and those in the post-puberty stage (p for interaction < 0.001) (Supplementary Figures S3 and S4 and Supplementary Table S2).Considering the possible impact of gender and puberty stage on the correlation tween protein intake and linear growth, we conducted further analysis on the associa between dietary protein intake and linear growth in different pubertal stages ac The association between protein intake and linear growth.Spline smoothing, based on a generalized additive model, was adjusted for gender ethnicity, age (as a continuous variable), pubertal stages, premature, birth length (as a continuous variable), maternal age (as a continuous variable), parental education levels (as a continuous variable), parental height (as a continuous variable), household income levels, smoking status, drinking status, moderate to vigorous physical activity, central obesity, length of nighttime sleep (as a continuous variable), dietary fat intake (as a continuous variable), dietary carbohydrate intake (as a continuous variable), and serum calcium, zinc, vitamin D 3 , and vitamin A levels (as continuous variables).(A) The association between protein intake (g•kg −1 •d −1 ) with body height; (B) The association between protein intake (g•kg −1 •d −1 ) with HAZ; (C) The association between protein intake (E%) with body height; (D) The association between protein intake (E%) with body height.

Discussion
In a large cross-sectional study using multistage stratified cluster random sampling, we observed for the first time a consistent negative association between protein intake and linear growth in children and adolescents after adjusting for potential confounders.Moreover, we found that higher protein intake was associated with a higher prevalence of stunting.
While previous studies have primarily focused on body weight, body composition, and related outcomes such as body mass index or fat mass index [10], our study specifically examined linear growth and its association with protein intake.Our findings contradict two earlier cross-sectional studies that reported a positive association between dietary protein intake and linear growth [21,22].Some prospective studies also reported a positive association between protein intake with height [23][24][25][26].For instance, Catherine S Berkey et al. contributed to research on the impact of milk on growth in a cohort of 5101 premenarchal girls from throughout the United States completed annual surveys (1996)(1997)(1998)(1999)(2000)(2001)2003).They found, that of the foods/nutrients studied, dairy protein had the strongest association with height growth.But the study finally suggested that the factor in the non-lipid phase of milk, but not the protein itself, has growth-promoting action in girls [23].Aliya Alimujiang et al. designed their study for a different purpose, but it also reported a positive association between animal protein intake and height at age 13, as well as adult height [24].It is noteworthy that the 229 participants were Caucasians of predominantly northern European descent from low-to middle-class families in Boston in the 1930s.The protein intake levels were not high.This was fundamentally different from the protein intake levels in our study.In contrast to our findings, Yifan Hua et al. also observed a direct relationship between protein intake and linear length in 189 girls [25].However, their average protein intake, assessed as a percentage of total energy intake (%E), was between 12.5% and 13.3%.In our study, the participants were from Shenzhen, a modern metropolis on the southeast coast of China, where the average protein intake was as high as 17.0%, 1.81 g•kg −1 •d −1 , significantly higher than the aforementioned study.According to the Recommended Dietary Allowance in the US, the recommended protein intake for age groups of 4-13 years and 14-17 years is 0.95 g•kg −1 •d −1 and 0.85 g•kg −1 •d −1 , respectively.Hence, the protein intake in our study significantly exceeded the recommended dietary protein intake in the US and closely resembled the current protein intake levels in children and adolescents in developed countries in Europe and the US [8].Therefore, we speculate that the previously observed discrepancies in associations between protein intake and linear growth might be attributed to different average protein intake levels.However, Kim Ve Braun et al. conducted a study where protein intake, comparable to our study at 1.9 g/kg BW per day, still showed a significant association.They found that a 10-g higher total protein intake per day at 1 year was significantly associated with a greater height of 0.03 SD [26].However, it was an early programming longitudinal study.As we know, individuals at a younger age have a greater demand for protein; therefore, a higher protein intake may not necessarily exert negative effects on early linear growth.However, two other studies did not find any significant association between protein intake and linear growth [23,27].
Therefore, the discrepancies may be attributed to variations in protein intake assessment methods (grams per day, grams per kilogram of body weight, or as a percentage of total energy intake), types of proteins studied (total proteins, animal proteins, plant proteins, etc.), and average protein intake levels.In this study, we primarily calculated protein intake in grams per kilogram of body weight to minimize any confounding effects of body weight status.Additionally, we assessed protein intake as a percentage of total energy intake (%E) to investigate whether different ways of measuring protein intake could explain previously observed discrepancies in associations with linear growth.We found a similar but weaker association between protein intake and linear growth when assessing protein intake as a percentage of total energy intake instead of grams per kilogram of body weight.Furthermore, we examined protein intake from both animal and plant sources separately, and both showed negative associations with body height and HAZ, similar to the associations observed with total protein intake.Therefore, neither the different methods of assessing protein intake nor the specific types of proteins accounted for the discrepancies between our study and previous studies.
In addition, a large population-based study in the US, which had relatively high protein intake, found that elevated protein levels during pregnancy were associated with shorter offspring birth length and slower linear growth into mid-childhood [15], and these findings align with our study results.
Although the precise molecular mechanisms underlying the relationship are unclear, our results are biologically plausible.On the one hand, a recent study demonstrated that a high maternal protein diet can impair offspring's bone mass through miR-24-1-5pmediated targeting of recombinant SMAD family member 5 (SMAD5) [28].SMAD5 is a well-characterized downstream component of the bone morphogenetic protein signaling pathway in osteoblasts, and its impairment leads to reduced osteoblast maturation and diminished bone mass [29].This finding suggests that high protein intake may negatively affect osteoblast maturation and bone mass.Therefore, this study provides a potential explanation for our results, pending confirmation in children and adolescents in future research.On the other hand, excessive protein intake can lead to an abundance of amino acids, which may accelerate amino acid oxidation and result in toxic levels of ammonia in the plasma.Ultimately, this could potentially hinder linear growth [30].
The strength of our study lies in being the first to observe a negative association between high protein intake and linear growth in children and adolescents, drawing more attention to this issue.Moreover, we assessed protein intake both as grams per kilogram of body weight and as a percentage of total energy intake (%E), which helps eliminate the influence of different methods of protein intake evaluation on the results.Furthermore, we adjusted for numerous confounders such as pubertal development stages, serum vitamin D 3 levels, serum vitamin A levels, serum calcium and zinc levels, physical exercise intensity, birth outcomes, parental height, etc., thus minimizing the impact of confounding factors on our results.Lastly, the large sample size of our study allowed us to conduct stratified analyses.
However, there are several limitations in our study.It is a cross-sectional study, making it difficult to establish causal relationships.Additionally, dietary assessments were performed using the FFQ, and although extensively validated in various populations [12], there may be some measurement error in estimating absolute nutrient intake.Nevertheless, we employed multiple approaches to address this error.Firstly, consistent results were obtained when ranking participants into quartiles of protein intake and when analyzing protein intake as a continuous variable.This method is considered valid for analyzing data obtained from FFQs [14].Furthermore, adjusting for total energy intake is a common technique for improving the validity of FFQs [14], and our nutrient intake estimates were all energy-adjusted using the nutrient residual method.

Conclusions
In conclusion, this large cross-sectional study in children and adolescents revealed that higher protein intake was associated with shorter linear growth and a higher prevalence of stunting.These findings suggest that increased dietary protein intake among wellnourished children and adolescents may not only fail to promote linear growth, but also impair it.Therefore, considering the potential risks of excessive protein intake is important.More studies are needed to confirm the causal relationship between them in the future.

Supplementary Materials:
The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nu15224821/s1, Figure S1: Flow of participants in current study.; Figure S2: The association of different sources of protein intake from with body height and HAZ.; Figure S3: The gender-specific association between dietary protein intake and linear growth.; Figure S4: The puberty stage-specific association between dietary protein intake and linear growth.;Table S1: Association of protein intake from different sources with body height and HAZ.; Table S2.Stratified analyses of potential modification effect for the association of dietary protein intake on linear growth.

Institutional Review Board Statement:
The study was conducted in accordance with the Declaration of Helsinki, and received approval from the ethics review committee of the Shenzhen Cen-ter for Chronic Disease Control (SZCCC-2021-037-01-PJ) and was registered at the China Clinical Trials Registry (ChiCTR2100051722; www.chictr.org.cn,accessed on 1 January 2022).
Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

Figure 1 .Table 2 .
Figure1.The association between protein intake and linear growth.Spline smoothing, based generalized additive model, was adjusted for gender ethnicity, age (as a continuous variable), bertal stages, premature, birth length (as a continuous variable), maternal age (as a continuous iable), parental education levels (as a continuous variable), parental height (as a continuous v ble), household income levels, smoking status, drinking status, moderate to vigorous physica tivity, central obesity, length of nigh ime sleep (as a continuous variable), dietary fat intake ( continuous variable), dietary carbohydrate intake (as a continuous variable), and serum calc zinc, vitamin D3, and vitamin A levels (as continuous variables).(A) The association between pro intake (g•kg −1 •d −1 ) with body height; (B) The association between protein intake (g•kg −1 •d −1 ) HAZ; (C) The association between protein intake (E%) with body height; (D) The association tween protein intake (E%) with body height.Table 2. Association of protein intake (per SD) with body height and HAZ.

Figure 1 .
Figure1.The association between protein intake and linear growth.Spline smoothing, based on a generalized additive model, was adjusted for gender ethnicity, age (as a continuous variable), pubertal stages, premature, birth length (as a continuous variable), maternal age (as a continuous variable), parental education levels (as a continuous variable), parental height (as a continuous variable), household income levels, smoking status, drinking status, moderate to vigorous physical activity, central obesity, length of nighttime sleep (as a continuous variable), dietary fat intake (as a continuous variable), dietary carbohydrate intake (as a continuous variable), and serum calcium, zinc, vitamin D 3 , and vitamin A levels (as continuous variables).(A) The association between protein intake (g•kg −1 •d −1 ) with body height; (B) The association between protein intake (g•kg −1 •d −1 ) with HAZ; (C) The association between protein intake (E%) with body height; (D) The association between protein intake (E%) with body height.
X.: visualization and validation.Y.L. (Yan Li), B.K. and Z.Z.: investigation and data curation.L.C., Y.T. and P.Y.: writing, reviewing, editing, and supervision.All authors have read and agreed to the published version of the manuscript.Funding: Ting Xiong is supported by the Young Scientists Fund of the National Natural Science Foundation of China (grant/award no.81903313), Guangzhou Science and Technology Plan Project (202201011543), and Plan on enhancing scientific research in GMU.Yuanjue Wu is supported by the Young Scientists Fund of the National Natural Science Foundation of China (grant/award no.82103835) and Plan on enhancing scientific research in GMU.Yanyan Li was supported by Shenzhen Fundamental Research Program (KCXFZ20201221173400002) and Shenzhen Medical Key Discipline Construction Fund (SZXK065).

Table 1 .
Association of characteristics of study participants with protein intake and outcomes (N = 3299).

Table 2 .
Association of protein intake (per SD) with body height and HAZ.

Table 3 .
The association of protein intake (per SD) with body height and HAZ in different pubertal stages across different genders.HAZ, height-for-age z score; SD, standard deviation.A generalized linear model was used.Protein intake was set to 1 SD increase of grams per kilogram weight per day (g•kg −1 •d−1 ). Adjusted for ethnicity, age (as a continuous variable), pubertal stages, premature, birth length (as a continuous variable), maternal age (as a continuous variable), parental education levels (as a continuous variable), parental height (as a continuous variable), household income levels, smoking status, drinking status, moderate to vigorous physical activity (MVPA), central obesity, length of nighttime sleep (as a continuous variable), dietary fat intake (as a continuous variable), and dietary carbohydrate intake (as a continuous variable), and serum calcium, zinc, vitamin D 3 , and vitamin A levels as continuous variables.