Introduction

The prevalence of obesity has been increasing worldwide and is considered to represent a pandemic situation requiring urgent action1,2,3. In 2016, more than 1.9 billion adults aged 18 years or older (corresponding to 39% of adults) were overweight and more than 650 million (corresponding to 13% of adults) were obese4. The risk of all-cause mortality increases even in overweight adults: every 5 unit increase in body mass index (BMI) above 25 kg/m2 is associated with an approximately 31% higher risk of mortality5. Thus, interventions are urgently needed to reduce the prevalence of overweight and obesity. The most important intervention for obesity is prevention (especially during childhood) rather than treatment6,7,8. Simmons et al. showed that about 55% of obese children remained obese during adolescence and about 80% of obese adolescents remained obese in adulthood. Therefore, interventions to reduce and prevent obesity during childhood and adolescence are needed. Understanding BMI trajectories during development can provide useful information for prevention efforts.

The BMI trajectory during development has been evaluated in many previous studies. However, most studies focused on BMI trends in children during segmented periods such as preschool, school age, or preadolescence. Only a few large cohort studies have evaluated BMI trends longitudinally from birth to adolescence9. In addition, although BMI is defined as weight (in kilograms) divided by height (in meters squared), few studies have considered the role of height changes in defining BMI trajectories10,11. Because puberty has been reported to occur earlier in obese children, accelerated height changes during puberty should be taken into consideration to understand BMI trends during that period12.

Moreover, there have been considerable racial differences observed in obesity studies based on BMI13. Thus, BMI trajectories in various racial groups must be delineated based on large longitudinal birth cohort studies. Obese Asian individuals have been found to have higher risks of hypertension and cardiovascular disease compared with obese white Europeans as well as higher risks of early death from cardiovascular disease or any cause14,15. However, few studies have characterized the BMI trajectories of Asian children and their relationships with obesity during adolescence16,17. The Lancet World Report 2007 already highlighted the growing epidemic of obesity in Japan: obesity is a critical concern for understanding future national patterns of disease. One in four men aged 20–69 years in Japan was obese in 2000, whereas this figure had risen to one in three men by 200718. Investigating the long-term BMI trends during childhood that lead to adolescent obesity is of major public health significance in Asia, especially in Japan.

The Longitudinal Survey of Newborns in the 21st Century is a national birth cohort study that has low susceptibility to cohort effects because it has been conducted for all births in Japan during specific weeks of 200119. In the present study, we investigated BMI trajectories from childhood to adolescence by BMI status at adolescence using data from this large nationwide birth cohort in Japan. We examined associations between the timing of puberty (age at peak height velocity, APHV) and obesity status in adolescence.

Methods

Participants

The Ministry of Health, Labour, and Welfare of Japan has been conducting The Longitudinal Survey of Newborns in the 21st Century since 2001 to establish strategies to counter the declining birthrate in Japan. The survey targeted all babies born in Japan between January 10 and 17 or between July 10 and 17 of 2001. Baseline questionnaires were sent to a total of 53,575 families when eligible babies reached the age of 6 months and 47,015 families initially completed the baseline questionnaire (88% response rate). These respondents were mailed follow-up questionnaires to investigate medical conditions and behaviors when children reached the ages of 1.5, 2.5, 3.5, 4.5, 5.5, 7, 8, 9, 10, 11, 12, 13, 14, and 15 years20,21,22,23. Birth record data from Vital Statistics of Japan are also linked for each child participating in the study. The current study included data for children/families who responded both to the baseline questionnaire and the fifteenth questionnaire at age 15 years.

The baseline survey at age 6 months included questions regarding children’s perinatal status as well as household and socioeconomic factors such as parental academic attainment, parental smoking status, and daycare attendance. The subsequent annual surveys starting at age 1.5 years included questions regarding each child’s height, weight and health status. We excluded 2382 children born before 37 weeks of pregnancy and one child with responses only for the baseline survey and the survey at age 15 years. A total of 26,778 children (315,581 data points) were included in the final analysis. A total of 11,141 children (41.61%) had responses to all 15 questionnaires between the ages of 6 months and 15 years, and responses to more than 12 questionnaires were available for the majority (91.94%) of children (Fig. 1, Table S1).

Figure 1
figure 1

Flowchart of study participants.

Measures

We calculated BMI based on each participant’s reported annual height and weight. Each participant’s annual BMI was converted to a BMI Z-score using smoothed L, M, and S values for BMI standards from a representative population of Japanese children24. Briefly, the LMS (lambda–mu–sigma) method is a method proposed by Cole et al. to monitor changes in the skewness of the distribution during childhood as a way of constructing normalized growth standards25. Participants were then classified into four BMI categories based on the World Health Organization (WHO) criteria26: underweight (BMI standard deviation [SD] score of − 5 or more but less than − 2), normal weight (BMI SD score of − 2 or more but less than 1), overweight (BMI SD score of 1 or more but less than 2), and obese (BMI SD score of 2 or more but less than 5). The definitions of overweight and obesity were different for children under 5 years of age: a BMI Z-score of 2 SD or more was categorized as overweight and a BMI Z-score of 3 SD or more was categorized as obese. BMI category at age 15 years was the main outcome of interest in the current study.

We also calculated annual height growth for each participant by subtracting the height reported at the previous survey from that reported in the current survey. For annual height growth between 5.5 and 7 years of age, this value was multiplied by 2/3 because of the 1.5-year interval between surveys.

Statistical analyses

We first compared baseline characteristics among the four BMI categories (underweight, normal weight, overweight and obese) at age 15 years. To evaluate potential selection bias resulting from losses to follow-up, we also compared the baseline characteristics of children included in the analysis and those of children lost to follow-up through to the fifteenth survey (at age 15 years).

We retrospectively examined annual aggregate categorical changes in individuals of the four BMI categories (groups) at age 15 years. For each group, the proportion of each BMI category at each survey between the ages of 1.5 and 14 years was calculated. In addition, we prospectively calculated the proportion of children in each BMI category at each survey between the ages of 1.5 and 14 years who eventually became underweight, normal weight, overweight, or obese at age 15 years. Note that these analyses were based on aggregate data and do not describe individual BMI changes and were performed using only the data obtained without imputation of missing values.

Under the assumption that missing data were missing at random, mixed effect models with natural cubic regression splines were applied to calculate the trajectories of BMI Z-scores and annual BMI Z-score changes through age 15 years for participants of each BMI category at age 15 years. Knots at seven locations were placed in percentiles of age to yield a sufficient number of measurements between each consecutive knot (age 1.5, 3.5, 5.5, 8.5, 11, 13 and 15 years), as recommended by Harrell27. The mixed effect model is useful for describing population average growth trajectories and individual growth trajectories even when data are not available for all children at all ages28,29,30,31. Briefly, the population average growth trajectory was modeled with fixed effects, while the individual variability is represented as random effects.

After fitting individual BMI trajectories using a mixed-effects model with natural cubic spline function, we estimated individual adiposity rebound timing as the age where the first derivative of the trajectory reached its minimum and the second derivative was positive32. Children were then classified into five categories (1.5–2.5 years, 3.5–4.5 years, 5.5–7 years, 8–10 years, and 11 years or older) for analysis of adiposity rebound timing33,34. The distribution of adiposity rebound timing was calculated for individuals of each BMI status at age 15 years overall and by gender.

Finally, we modelled annual height change and its associations with BMI status at age 15 years separately for each gender using mixed-effects models with natural cubic regression splines.

All statistical analyses were performed using Stata version 16 (StataCorp LLC, College Station, TX, USA). This study was approved by the Institutional Review Board at Okayama University Graduate School of Medicine, Dentistry, and Pharmaceutical Sciences (No.1506-073) and was conducted in accordance with the 1964 Helsinki Declaration and Ethical Guidelines for Medical and Health Research Involving Human Subjects. Informed consent was obtained by the opt-out method on the university's website.

Results

Demographic characteristics

Participants’ demographic characteristics according to BMI status at age 15 years are shown in Table 1. Obese adolescents tended to be boys, to be large for gestational age at birth, to live in towns or villages, to have parents with lower academic attainment, and to have mothers who smoked. During the follow-up period, 17,854 children were lost to follow-up by the fifteenth survey (at 15 years of age). Children lost to follow-up tended to have younger mothers, mothers who smoked, and mothers with lower academic attainment (Table S2).

Table 1 Demographic characteristics of children included in the analysis at age 1.5 years by BMI status at age 15 years (N = 26,778).

Categorical aggregate changes in each BMI status group

Figure 2A shows the results of a retrospective analysis whereby we calculated the percentages of children in the four BMI categories (underweight, normal weight, overweight, or obese) every year during childhood according to their BMI group at age 15 years. Children with normal weights at age 15 years mostly maintained normal weights throughout childhood. Although 83.1% of the children who were obese at age 15 years had normal weights at age 1.5 years, the proportion of overweight or obese children increased annually, with a large percentage of children becoming obese after age 13 years. Figure 2B shows the results of a prospective analysis whereby we calculated the proportion of children in each BMI category at each survey (from ages 1.5–14 years) who subsequently became underweight, normal weight, overweight, or obese at age 15 years. Overall, 31.0% of 7-year-old obese children had normal weights at age 15 years. The proportion of overweight/obese children who returned to normal weights by age 15 years gradually decreased, and markedly decreased after the age of 12 years. Only a small proportion of underweight/normal weight children in earlier surveys became overweight/obese at 15 years of age.

Figure 2
figure 2

Annual categorical body mass index (BMI) changes by BMI category at age 15 years. (Panel A) Retrospective tracking of BMI status during childhood (age 18 months to 14 years) according to BMI status at age 15 years: underweight (A), normal weight (B), overweight (C) and obese (D). (Panel B). Prospective tracking of annual BMI status [underweight (A), normal weight (B), overweight (C), and obese (D)] from childhood (age 18 months to 14 years) to adolescence (age 15 years). BMI status categorization was based on the WHO definitions (under 5 years: overweight ≥ 2 SD, obese ≥ 3 SD; over 5 years: overweight ≥ 1 SD, obese ≥ 2 SD).

BMI status and BMI changes during childhood

The average trajectories of BMI Z-scores for boys and girls are shown in Fig. 3A. These trajectories depict the fixed effects component using mixed effects models with natural cubic splines. The average BMI Z-score trajectories of children with normal weights at age 15 years remained stable around 0 throughout childhood, whereas children who were overweight/obese at age 15 years already had relatively high BMI Z-scores by 1.5 years of age. The average trajectory for BMI Z-scores in participants who were overweight/obese at age 15 years showed a continuous increase in both genders throughout childhood, with a greater slope during puberty. Children who were underweight at age 15 years already had relatively low BMI SD scores at 1.5 years of age and, in contrast to the trajectory for participants who were obese at age 15 years, showed a marked decline in slope after puberty. Comparing the average trajectories of annual change in BMI Z-scores (Fig. 3B), participants of both genders who were obese at age 15 years showed a less pronounced dip around age 5 years than the other groups, a continuous increase in BMI Z-scores across ages, and a greater slope at puberty. By contrast, the average trajectory of annual change in BMI Z-score in participants who were overweight at age 15 years was similar to that of participants who had normal weights at age 15 years, albeit with relatively larger changes in the overweight group compared with the normal weight group.

Figure 3
figure 3

Dynamics of BMI Z-scores (A) and annual BMI Z-score changes (B) through age 15 years by gender.

Impact of adiposity rebound timing

We compared adiposity rebound timing by BMI status at age 15 years overall and by gender (Table 2 and Table S3). Adiposity rebound occurred earlier in participants who were overweight/obese at age 15 years (prior to age 4.5 years) than in those who had normal weights at age 15 years. Moreover, more than 95% of participants who were obese at age 15 years, had experienced adiposity rebound before 2.5 years of age. In contrast, adiposity rebound tended to occur later in participants who were underweight at age 15 years.

Table 2 Timing of adiposity rebound and BMI status at age 15 years.

BMI status at age 15 years and APHV

Of the 26,778 participants included in the analysis, we excluded eight children whose annual height gain was never measured (i.e., no two consecutive responses). We used mixed-effects models with natural cubic regression splines to calculate the fixed-effects portion of the trajectory of annual height gain for participants of each obesity status at age 15 years, by gender (Fig. 4). Among boys, the APHV occurred earliest in participants who were obese at age 15 years, followed by those who were overweight, normal weight, and underweight at age 15 years. A similar trend was observed for girls with no marked differences between those who were obese and overweight at age 15 years.

Figure 4
figure 4

Dynamics of annual height growth (cm) through age 15 years.

Discussion

In the present study, we delineated the BMI trajectories leading to obesity in adolescents and examined associations between BMI status in childhood and obesity at age 15 years using data from a large birth cohort of all Japanese children born during specific weeks of 2001. The role of annual height gain on BMI trajectories in children and adolescents was also evaluated.

Our data regarding changes in BMI during childhood are partially consistent with the findings of a German population-based study examining BMI trends from childhood to adolescence (age 15–18 years). Mandy et al. reported that BMI acceleration (i.e., a rapid increase in BMI) during childhood increased the risk of obesity in adolescence and that almost 90% of children who were obese at 3 years of age remained overweight or obese in adolescence35. In the present study, we found that adiposity rebound timing occurred earlier in participants of both genders who were overweight/obese at age 15 years; this difference was especially marked in those who were obese at age 15 years. The population average trajectories of BMI Z-score change among participants who were obese at age 15 years also showed an increase in BMI Z-scores over time, with no dip observed at preschool age in either gender. In our study, only 17.5% of children who were obese at age 5.5 years remained obese at the age of 15 years, and more than half of participants who were obese at age 15 years were overweight at age 13 years. Analysis of the population average trajectory for participants of each BMI status at age 15 years showed that unlike those who had normal weights at age 15 years, participants of both genders who were obese at age 15 years maintained relatively high BMI Z-scores throughout childhood, with an increasing trend over time and an increasing slope during puberty. This rapid increase in BMI Z-score during adolescence (age 14–15 years) was not observed in a previous German study. Unlike some prior studies, we included participants from a large nationwide population-based study for whom data were collected annually. In contrast with a previous German study, in which data were available for 13 or more time points in only 1% of participants, in our study the majority of children (91.74%) had responses for more than 12 surveys and 11,093 children (41.53%) had responses for all 15 surveys between the ages of 6 months and 15 years. On the basis of these comprehensive data, we were able to model BMI trajectories using multilevel models with natural cubic splines and depict the increases in BMI characteristic of obese adolescents.

Racial differences between study populations may explain some of the discrepancy observed between studies of BMI trajectories and obesity. For example, a follow-up study conducted in northern China identified a subgroup of children with a rapid increase in BMI after the onset of puberty36 These findings suggested that Chinese adolescents with overweight or obesity experienced BMI acceleration at two time points: at preschool ages and during puberty. Although several studies of BMI trajectories have included Asian participants, few studies have examined childhood BMI trajectories by BMI status in adolescence in a large cohort of children followed from birth until puberty. Liang et al. assessed the BMI trajectories of children aged 2–18 years using group-based trajectory modeling with random sampling from five cohorts in China. Their study mainly focused on social factors related to obesity and BMI trajectories could not be compared by BMI status in adolescence because the study included children from various ages and backgrounds. Haga et al. followed 1644 children born during an 8-year period in a district of Japan until age 12 years and identified five latent class patterns in boys and six latent class patterns in girls using latent class growth modeling. However, few large studies have longitudinally tracked BMI from birth to adolescence. The methodology used in our study would be expected to be less susceptible to cohort effects because children were born around the same time9.

In study of BMI, Sheila et al. assessed the influence of height gain on early adiposity rebound. BMI during puberty is expected to be affected by height gain. Several studies have shown that early adiposity rebound indicates faster growth, more advanced development, and earlier puberty37,38,39,40. In fact, puberty has been reported to occur earlier in obese individuals, and differences in the timing of puberty may have accentuated the increase in BMI Z-scores observed after age 13 years in obese children10,11. Adolescence, characterized by changes in body composition, physical fitness, and decreased insulin sensitivity during puberty, is a critical period for preventing the onset and continuation of obesity throughout the lifespan41,42,43. Ohlsson et al. showed that increased BMI through puberty and adolescence, but not in childhood, was associated with risk of adult stroke44. Further long-term studies are needed to assess the impact of BMI acceleration in adolescence on obesity and disease risk; at the time of BMI assessment, differences in acceleration of height growth based on childhood BMI status should be considered. Therefore, we analyzed the trajectory of annual height gain in this study. The population average trajectories for annual height gain by BMI status at age 15 years showed that APHV occurred earlier in participants who were obese/overweight at age 15 years and later in participants who were underweight at age 15 years compared with those who had normal weights at age 15 years. This phenomenon may partially explain why BMI Z-score trajectories in adolescence diverge by BMI status at age 15 years.

To date, few studies have considered the role of height when examining BMI trajectories10,11,12, especially in studies of Asian children16,17. Japan has been noted as a country with a rapidly growing obesity epidemic. We expect that our report will provide valuable insights for the prevention of obesity18.

Our study had several limitations. First, information on maternal history of obesity was unavailable. Since individual genetic predisposition and dietary habits can affect the risk of obesity45,46, future studies that include these data may identify additional group traits contributing to adolescent obesity. Second, we did not consider fat mass index and focused only on BMI, which may have resulted in misclassification of adiposity rebound timing47. However, this misclassification would likely be non-differential and bias effect estimates toward the null. Third, information on height and weight was obtained on the basis of parental reports rather than clinical measurement, which may have introduced measurement errors. Self-reported BMI may overestimate BMI in underweight individuals and underestimate BMI in overweight/obese individuals48. Fourth, some participants were lost to follow-up, which may have introduced selection bias. Children lost to follow-up (who tended to have younger mothers, mothers who smoked, and mothers with lower academic attainment) may have been at higher risk for overweight/obesity, and thus selection bias might have reduced the number of overweight/obese children in our study. Finally, we targeted Japanese children, which might limit generalizability to other populations.

In conclusion, our study using data from a Japanese national birth cohort showed that obese adolescents experienced early adiposity rebound timing and maintained serial BMI Z-score increases throughout childhood, with a greater slope during puberty. An earlier peak in height gain during puberty may have contributed to the observed patterns of BMI change.