Associations Between Delivery Mode and Early Childhood Body Mass Index Z-Score Trajectories: A Retrospective Analysis of 2,685 Children From Mothers Aged 18 to 35 Years at Delivery

Objective: To investigate the association between cesarean delivery (CD) and trajectory patterns of age- and sex-specific body mass index (BMI) z-score in early childhood.

Methods: A retrospective cohort study was conducted among 2,685 children whose maternal age at the time of birth was between 18 and 35 years, and birth data and anthropometric measurement data during their ages 3–60 months were collected. A group-based trajectory modeling approach was used to identify distinct BMI z-score trajectories, and multinomial logistic regressions were applied to estimate the associations among CD (both elective and non-elective combined), elective and non-selective CD, and BMI z-score trajectory classes.

Results: Of the 2,685 participants, 46.5% (N = 1,248) were born by vaginal delivery (VD), 20.7% (N = 556) by elective CD, and 32.8% (N = 881) by non-elective CD. Five BMI z-score trajectory patterns were identified, and they were “increasing from moderate to high” (10.1%, n = 270), “increasing from mild to moderate” (34.2%, n = 919), “increasing from low to high” (10.5%, n = 283), “stable mild” (30.1%, n = 808), and “stable low” (15.1%, n = 405) groups. Compared with children delivered by VD, those who delivered by CD (both elective and non-elective combined), elective CD, and non-elective CD were associated with the “increasing from moderate to high” trajectory [odds ratio (OR) = 1.61, 95% confidence interval (CI): 1.13–2.29; OR = 1.64, 95%CI: 1.06–2.54; and OR = 1.59, 95%CI: 1.05–2.39, respectively] and were also associated with the “increasing from low to high” trajectory (OR = 1.60, 95%CI: 1.17–2.19, OR = 1.75, 95%CI: 1.16–2.63; and OR = 1.53, 95%CI: 1.00–2.34, respectively).

Conclusion: Both elective and non-elective CD were associated with the risk of accelerated weight gain in early childhood.

Introduction

The cesarean delivery (CD) rate has increased significantly in the past decades worldwide (1). Ranking as one of the countries with the highest CD rate, China has witnessed a more rapid increase in CD rate, which increased from 28.8% in 2008 to 36.7% in 2018 with an annual percentage change of 1.8% (2). CD could be life-saving for both the mothers and their offspring if appropriately performed and medically indicated, which should only be performed when a vaginal delivery (VD) places the mothers or newborns at risk because of medical reasons (3). However, many Chinese mothers chose a CD without any medical necessities for the reasons of fearing pain by VD, perceiving CD safer than VD, choosing a lucky day as babies' birthdate, and obstetricians' beliefs or attitudes toward CD (4, 5). Given the high rate of CD, especially elective CD, the potential adverse impact of CD on mother and offspring has been concerning, and in comparison with VD, higher risk of maternal death or infection, neonatal respiratory morbidity, and adverse effects to subsequent pregnancies have been identified (6). In addition, CD may result in altered neonatal gut microbiota and thus pre-dispose offspring to an increased risk of childhood overweight and obesity (7).

Childhood overweight and obesity have become a global public health concern, and numerous studies have shown that there is a distinct rise in overweight and obesity prevalence in both urban and rural regions of China (8). Childhood overweight and obesity are associated with short-term risk of cardiovascular risk factors [such as elevated blood pressure (9), adverse lipid profiles (10)]. Noteworthily, there is accumulating evidence suggesting that childhood obesity tends to track into adulthood and subsequently increase the long-term risk of adult cardiovascular and metabolic diseases (11, 12).

Some studies have investigated the association between CD and childhood obesity. Although some studies have reported that CD is significantly associated with risk of childhood overweight/obesity or higher body mass index (BMI) (13–15), heterogeneous results presenting null associations have also been observed (16). In addition, studies exploring the association between CD subtypes (elective and non-elective CD) and indicators of childhood obesity also reported inconsistent results (17, 18). Furthermore, most existing studies focused on BMI in a single or limited number of measurements, ignoring its dynamic changes during a specific period.

In response to these literature gaps, we designed a retrospective cohort study to investigate the association between CD as well as its subtypes (elective and non-elective CD) and trajectory patterns of age- and sex-specific BMI z-score in children aged 3–60 months, which could overcome the limitation of static analysis and increase comparability across ages.

Methods

Study Population

We designed this retrospective cohort study in Huangdao Maternal and Child Health Care Hospital, which provides prenatal care, delivery, postnatal, and child health care service for all residents in Huangdao district, Qingdao, China. This study is based on the National Basic Public Health Service Project (http://www.nbphsp.org.cn/jbgw/06rt/), which was launched by the Chinese government in 2009. In this project, one of the nine service categories was providing health services for Chinese children aged 0–6 years for free, including the establishment of maternal and children health records, post-partum visits, physical examinations, and health consultation. Health providers involved in this project were all trained according to national standards and required to provide elaborate health services for children who were at ages of 3, 6, 8, 12, 18, 24, 30, 36, 48, 60, and 72 months.

Participants were singleton offspring born at term in this hospital between January 1, 2012, and December 31, 2013, and had height/length and weight repeatedly measured 4 to 10 times during the following 5 years. Inclusion criteria were (1) singleton offspring born with 37–42 completed weeks of gestation and (2) mother's age at birth between 18 and 35 years. Exclusion criteria were (1) mothers who were already diagnosed with types 1 or 2 diabetes before the index pregnancy or (2) offspring conceived by in vitro fertilization, born with congenital disabilities or postnatal diseases that could interfere with body composition development, or (3) offspring with missing covariate data. The study was conducted in accordance with the Declaration of Helsinki, and the protocol and the use of the database were approved by the Ethics Committee of the Huangdao Maternal and Child Health Care Hospital.

Exposure and Covariates

Delivery mode (VD, elective CD, or non-elective CD) data and the following individual-level covariates were collected from the birth records: maternal age at birth, maternal education level, parity, hypertensive disorders of pregnancy, gestational diabetes mellitus, gestational age at birth, and birth weight. Maternal BMI was measured at the first prenatal visit, and BMI status was defined by using Asian BMI cut points (underweight: BMI <18.5, normal: 18.5 ≤ BMI <23.0, and overweight/obesity: BMI ≥ 23.0) (19).

Outcomes

Age- and sex-specific BMI z-score trajectory patterns were the outcomes of interest, and anthropometric measurements were performed by trained health providers for participants dressed in light indoor clothing without shoes. Recumbent length for children aged <24 months and standing height for children aged ≥24 months were measured to the nearest 0.1 cm with an infantometer and a stadiometer, respectively. Body weight was measured to the nearest 0.1 kg using a digital scale. All of the instruments were calibrated before measuring, and both length/height and weight were measured twice; then, the mean values were recorded. BMI was calculated as weight (kilogram) divided by height/length (square meter), and age- and sex-specific BMI z-score was calculated using the “igrowup”package for SAS based on World Health Organization Child Growth Standards (20).

Statistical Analysis

A group-based trajectory modeling approach was applied to identify trajectories of BMI z-score from aged 3 months through 60 months. This approach has been widely used in other childhood populations across longitudinal data to identify several discrete classes and individuals in each class following similar trends of change (21). To fully reflect the potential BMI z-score trajectories, only an individual having BMI measured at least four time points out of ages of 3, 6, 8, 12, 18, 24, 30, 36, 48, and 60 months could be included in this analysis. Analyses were conducted via PROC TRAJ macro in SAS (22), and the time scale was age at anthropometric measurement. Repeated trajectory analyses were conducted to identify the latent classes by setting the number of groups from 2 to 7, and model fit was assessed based on the following criteria: Bayesian Information Criterion, the mean posterior probability for group membership >0.7, and clinical plausibility and the interpretability of the classes (23, 24). Once the optimal number was determined, linear, quadratic, and cubic terms were tested to determine the appropriate shape of each trajectory (22). Then, participants were assigned to the trajectory group for which they had the greatest posterior predictive probability.

Differences in maternal and offspring characteristics among the identified BMI z-score trajectory groups were examined using one-way analysis of variance or chi-square-test. To examine the association between delivery mode and BMI z-score trajectory patterns, crude and adjusted multinomial logistic regression models were used to report OR with 95% CI. We constructed three models: model 1 was the crude model, model 2 adjusted for gestational age at birth and birth weight, and model 3 adjusted for gestational age at birth, birth weight, maternal age at birth, maternal education level, maternal BMI at the first prenatal visit, parity, hypertensive disorders of pregnancy, and gestational diabetes mellitus. We further conducted a stratification analysis by child sex to test whether there is a sex-specific effect on BMI z-score trajectory patterns. All statistical analyses were conducted using SAS version 9.4 (SAS Institute Inc., Cary, NC, USA).

Results

Study Population Characteristics

A total number of 3,405 birth records in the years 2012 and 2013 were retrieved from the birth records, and 2,685 offspring were included in this study according to the inclusion and exclusion criteria. The characteristics of the studied mother–child pairs are shown in Table 1 and the flow chart in Figure 1. In this analysis population, 46.5% (N = 1,248) of the children were born by VD, 20.7% (N = 556) by elective CD, and 32.8% (N = 881) by non-elective CD. The mean maternal age at birth was 28.2 years, and the prevalence of maternal obesity/overweight was 29.8%. In addition, 52.3% (N = 1,404) of the children were boys, and the mean gestational age and birth weight were 39.1 weeks and 3.06 kg, respectively. We compared the characteristics of the excluded children and the children having <4 anthropometric measures with the characteristics of the included children, respectively (Supplementary Table 1). Excluded children tended to be from older mothers and having lower birth weight, and children having <4 anthropometric measures tended to be from mothers with lower education level.

TABLE 1

Table 1. Participant characteristics for each of the five different BMI z-score trajectory groups.

FIGURE 1

Figure 1. Flow chart describing the selection process of subjects for the present study.

Group-Based Trajectory Modeling

Children had a median of six repeated anthropometric measurements (interquartile range: 5–8 times) from 3 to 60 months of age. Numbers of children who had anthropometric measurements at the age of 3, 6, 8, 12, 18, 24, 30, 36, 48, and 60 months were 2,136, 2,029, 1,767, 1,938, 1,713, 1,725, 1,608, 1,627, 1,534, and 1,489, respectively. The five-group model proved optimal for describing BMI z-score trajectory patterns from age 3 months through 60 months (Figure 2). We labeled the five trajectories as “increasing from moderate to high” (10.1%, n = 270), “increasing from mild to moderate” (34.2%, n = 919), “increasing from low to high” (10.5%, n = 283), “stable mild” (30.1%, n = 808), and “stable low” (15.1%, n = 405) groups. Among the five groups, significant differences in proportions of delivery mode, maternal BMI status, maternal education level, gestational diabetes mellitus, parity, offspring sex, and means of gestational age and birth weight were observed (Table 1).

FIGURE 2

Figure 2. Trajectories of age- and sex-specific body mass index z-score for children aged 3 to 60 months.

Association of Delivery Mode With Body Mass Index z-Score Trajectory Groups

In the unadjusted model, compared with children delivered by VD, those who delivered by CD were significantly associated with the two trajectories characterized by accelerated BMI z-score gain (for “increasing from moderate to high” group: OR = 1.73, 95%CI: 1.29–2.32, for “increasing from low to high” group: OR = 1.89, 95%CI: 1.41–2.53), and results were similar for both elective and non-elective CD (Table 2). In addition, the non-elective CD was significantly associated with the “stable low” group (OR = 1.34, 95%CI: 1.03–1.79). After adjustment for potential confounding factors representing maternal and offspring characteristics, compared with children delivered by VD, those who delivered by CD (both elective and non-elective combined), elective CD, and non-elective CD were still significantly associated with the risk of having the trajectory of “increasing from moderate to high” group (OR = 1.61, 95%CI: 1.13–2.29; OR = 1.64, 95%CI: 1.06–2.54; and OR = 1.59, 95%CI: 1.05–2.39, respectively) and the trajectory of “increasing from low to high” group (OR = 1.60, 95%CI: 1.17–2.19; OR = 1.75, 95%CI: 1.16–2.63, respectively; OR = 1.53, 95%CI: 1.01–2.34, respectively). In the stratification analysis by child sex (Supplementary Table 2), the effect estimates were similar to the main results, revealing no sex-specific effect on BMI z-score trajectory patterns.

TABLE 2

Table 2. Crude and multivariable-adjusted odds ratios (OR) and 95% confidence intervals (CI) for offspring BMI z-score trajectories associated with cesarean delivery as well as its subtypes vs. vaginal delivery.

Discussion

In this retrospective cohort study, we identified five distinct BMI z-score trajectories during early childhood and found significant associations of CD and its subtypes with two trajectories characterized by accelerated BMI z-score gain. Our findings suggest that newborns delivered by CD, irrespective of subtypes, would have an accelerated weight growth pattern during early childhood. To the best of our knowledge, this is the first study to explore the association between delivery mode and BMI z-score trajectory patterns during early childhood. This study provides new insights into the origins of overweight and obesity in early childhood life and emphasizes the importance of reducing the overused cesarean section.

In this population, the prevalence of maternal obesity/overweight is 29.8%, which is similar to the prevalence of obesity/overweight in women aged 20–39 years from the 2011 China Health and Nutrition Survey (26.3%) (25). A total of 53.5% of children were delivered by CD during 2012–2013 in this hospital from an eastern coastal district, which was higher than the mean CD rate in China in 2013 (34.6%) (26). However, there was distinct geographic variation in CD rates across China, and the CD rates of Shanghai and Hangzhou (two eastern coastal cities) were also larger than 50% during 2012–2013 (26). All of these regions ranked in the highest tertile economic level in China, and the expanded access to hospital care and financial incentives for physicians to perform cesareans may lead to a higher CD rate (26).

Our association analyses showed that children delivered by CD were associated with an increased risk of having trajectories characterized by accelerated weight gain, although the effect estimates attenuated distinctly after adjustment for maternal and offspring variables. This result is similar to previous studies reporting positive associations between CD and larger childhood BMI or risk of overweight/obesity (13–15). However, two cohort studies revealed inconsistent results. Pei et al. reported a significant association between CD and childhood obesity risk only at 2 years (OR = 1.68, 95% CI: 1.10–2.58) but not at age 6 or 10 years in a German cohort (27). Vinding et al. (28) reported a higher mean BMI at 6 months (β = 0.37, 95% CI: 0.14–0.60) in children delivered by CD but not at 5 or 13 years in two Danish cohorts. In contrast, one meta-analysis found that CD was associated with overweight and obesity not only among children (3–8 years) but also in adolescence (9–18 years) and adult (≥19 years) populations (29).

With regard to the effects of different subtypes of the CD on offspring overweight and obesity, existing studies also showed discordant findings. Both elective and non-elective CDs were found to be associated with obesity in children aged 3–6 years from China (18) and in children aged 7 years from the United States (17). However, Cai et al. (30) reported that only elective CD, but not non-elective CD, was associated with risk of overweight in Singaporean children aged 12 months, and Smithers et al. (16) found elective CD was not associated with BMI z-score among Australian children aged 3–6 years. Differences in age groups of subjects and definition of outcomes of interest in analytic models among the studies mentioned earlier may explain the discrepant results. Moreover, smaller sample size and CD and obesity rates in some studies may also reduce the power to detect a significant association between CD and outcomes.

The group-based trajectory modeling approach in our study makes it possible to explore the potential heterogeneity of BMI z-score changes during early childhood, which is an essential difference against the studies mentioned earlier. Several previous studies have identified the role of stable high BMI or rapid BMI gain during childhood contributing to increased risk for obesity in adulthood, suggesting the potential ability of childhood BMI trajectories in predicting health outcomes in later life (23, 31). These findings underline the necessity to explore influencing factors leading to those seemingly unhealthy BMI trajectories.

The exact mechanisms of CD leading to obesity or accelerated BMI gain during childhood is still unclear, but alterations in the gut microbiota via different delivery modes are a potential biological mechanism speculated by many researchers. The fetus and the intrauterine environments are believed to be sterile, and thus, the intestinal microbiome developed beginning at birth may be dependent on whether the newborn's first contact was with the vaginal canal or skin microbiota (17). Several studies have reported differences in gut microbiota between overweight or obese children with children in a normal weight status (32, 33). Bacteroides genera, which could protect against overweight by influencing gut energy harvest from the diet (12), were found to be lower in offspring delivered by CD compared with their VD counterparts (34). In addition, mothers who experienced a cesarean operation generally have higher odds of breastfeeding initiation failure or shorter breastfeeding duration, which has been proven to be a risk factor for childhood obesity (35). Therefore, lower Bacteroides genera and poorer maternal breastfeeding practices attributed to CD may be associated with a higher risk of abnormally accelerated growth trajectory. However, the lack of these data impeded us from exploring the exact mechanisms. In summary, the fat and glucose metabolism of offspring may be modified due to the earlier mentioned influencing factors either independently or in conjunction with each other, followed by disturbing their appetite regulation or metabolic control, and thus, the BMI growth trajectory pattern will be altered (27).

This study has several strengths. The retrospective cohort design with sufficient repeated anthropometric measurements and the group-based trajectory modeling approach enabled us to classify children into diverse growth trajectories. Besides, BMI z-score other than BMI was used to model the growth trajectory, which allowed for estimating how individuals tracked according to their rank order within BMI distributions and exempting the need for age and sex adjustment and also providing the most meaningful measure to investigate the relationship between delivery mode and change in weight status (36).

Limitations

First, the development of obesity in childhood is a multifactorial phenomenon, but data on some confounding factors, such as maternal passive smoking, maternal breastfeeding pattern, child daily diet structure, physical and sedentary activity, and family income level, are absent in this study (35, 37, 38). However, several previous studies presenting both unadjusted and adjusted results revealed that adjustment for some environmental (antenatal active or passive smoking), breastfeeding (duration of breastfeeding), and economic (family income) factors had only a slight influence on the estimated effects (18, 30). Therefore, the impact of these unmeasured confounding factors on our results may be minimal. Second, we did not model trajectories for boys and girls separately, given that the growth curves were already adjusted by sex. Third, we included mothers only aged 18–35 years at the time of birth, and older mothers were not included in this study. However, advanced maternal age is defined as childbearing in a woman more than 35 years of age, and this is a well-established risk factor for neonatal morbidity and future childhood development (39, 40). Therefore, we excluded women >35 years of age at delivery aiming to control the complicated effects of advanced maternal age. Lastly, this study collected samples only from one hospital in a Chinese coastal district, and ~90% of the mothers were primiparous; thus, we cannot rule out the possibility that our results might be biased and not be generalizable to a population in other regions. Further prospective multiple-center studies collecting more information on maternal and offspring characteristics, such as diet, feeding practices, and environmental and economic factors that may impact the infant/child's growth trajectory, are warranted.

Conclusion

The control of CD rate, especially the elective CD rate, should be considered as a priority goal for reducing childhood obesity and subsequently decreasing the risk of cardiovascular and metabolic diseases in later life. Communities, antenatal care sections, and hospitals should promote effective health education for pregnant women.

Data Availability Statement

The data analyzed in this study is subject to the following licenses/restrictions: The data underlying the study are not available to the public as the project does not meet the budget requiring data sharing. Requests to access these datasets should be directed to Xiao Jia, amlheGlhbzIwMTkmI3gwMDA0MDsxNjMuY29t.

Ethics Statement

The studies involving human participants were reviewed and approved by The Ethics Committee of the Huangdao Maternal and Child Health Care Hospital.

Author Contributions

LZ, WQ, and XJ conceptualized the study and drafted the article. LZ, LH, and RD collected all data with initial pre-processing. HL, ZZ, and LZ analyzed data. ZZ and LZ also made contributions to revising the manuscript. All authors approved the final submitted and published versions.

Funding

This study was sponsored by the Foundation of Shandong Provincial Science and Technology Development (Nos. 2018GSF118229 and 2019GSF108204).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We acknowledge the dedication of the nursing and medical staff of the participating centers and the study participants, without whom this study could not have been completed.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fped.2020.598016/full#supplementary-material

References

1. Betrán AP, Ye J, Moller AB, Zhang J, Gülmezoglu AM, Torloni MR. The increasing trend in caesarean section rates: global, regional and national estimates: 1990–2014. PLoS ONE. (2016) 11:e0148343. doi: 10.1371/journal.pone.0148343

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Li HT, Hellerstein S, Zhou YB, Liu JM, Blustein J. Trends in cesarean delivery rates in China, 2008–2018. JAMA. (2020) 323:89–91. doi: 10.1001/jama.2019.17595

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Zheng YY, Pubu ZM, Ge Y, Bianba ZM, Xu XH. Prevalence and causes of cesarean births among women residing at high altitude in Lhasa, Tibet: a retrospective observational study of 7365 women. Taiwan J Obstet Gynecol. (2019) 58:380–84. doi: 10.1016/j.tjog.2019.01.030

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Sun N, Yin X, Qiu L, Yang Q, Gong Y. Chinese obstetricians' attitudes, beliefs, and clinical practices related to cesarean delivery on maternal request. Women Birth. (2020) 33:e67–71. doi: 10.1016/j.wombi.2019.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Mi J, Liu F. Rate of caesarean section is alarming in China. Lancet. (2014) 383:1463–64. doi: 10.1016/s0140-6736(14)60716-9

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Yu Y, Zhang X, Sun C, Zhou H, Zhang Q, Chen C. Reducing the rate of cesarean delivery on maternal request through institutional and policy interventions in Wenzhou, China. PLoS ONE. (2017) 12:e0186304. doi: 10.1371/journal.pone.0186304

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Tun HM, Bridgman SL, Chari R, Field CJ, Guttman DS, Becker AB, et al. Roles of birth mode and infant gut microbiota in intergenerational transmission of overweight and obesity from mother to offspring. JAMA Pediatr. (2018) 172:368–77. doi: 10.1001/jamapediatrics.2017.5535

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Song Y, Agardh A, Ma J, Li L, Lei Y, Stafford RS, et al. National trends in stunting, thinness and overweight among Chinese school-aged children, 1985–2014. Int J Obes (Lond). (2019) 43:402–11. doi: 10.1038/s41366-018-0129-7

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Dong Y, Jan C, Zou Z, Dong B, Wang Z, Yang Z, et al. Effect of overweight and obesity on high blood pressure in Chinese children and adolescents. Obesity (Silver Spring). (2019) 27:1503–12. doi: 10.1002/oby.22562

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Brzezinski M, Metelska P, Mysliwiec M, Szlagatys-Sidorkiewicz A. Lipid disorders in children living with overweight and obesity- large cohort study from Poland. Lipids Health Dis. (2020) 19:47. doi: 10.1186/s12944-020-01218-6

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Weihrauch-Bluher S, Schwarz P, Klusmann JH. Childhood obesity: increased risk for cardiometabolic disease and cancer in adulthood. Metabolism. (2019) 92:147–52. doi: 10.1016/j.metabol.2018.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Geserick M, Vogel M, Gausche R, Lipek T, Spielau U, Keller E, et al. Acceleration of BMI in early childhood and risk of sustained obesity. N Engl J Med. (2018) 379:1303–12. doi: 10.1056/NEJMoa1803527

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Chu S, Zhang Y, Jiang Y, Sun W, Zhu Q, Liu S, et al. Cesarean section and risks of overweight and obesity in school-aged children: a population-based study. QJM. (2018) 111:859–65. doi: 10.1093/qjmed/hcy195

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Goldani MZ, Barbieri MA, da Silva AA, Gutierrez MR, Bettiol H, Goldani HA. Cesarean section and increased body mass index in school children: two cohort studies from distinct socioeconomic background areas in Brazil. Nutr J. (2013) 12:104. doi: 10.1186/1475-2891-12-104

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Pluymen LP, Smit HA, Wijga AH, Gehring U, De Jongste JC, Van Rossem L. Cesarean delivery, overweight throughout childhood, and blood pressure in adolescence. J Pediatr. (2016) 179:111–17.e3. doi: 10.1016/j.jpeds.2016.08.059

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Smithers LG, Mol BW, Jamieson L, Lynch JW. Cesarean birth is not associated with early childhood body mass index. Pediatr Obes. (2017) 12(Suppl. 1):120–24. doi: 10.1111/ijpo.12180

CrossRef Full Text | Google Scholar

17. Mueller NT, Whyatt R, Hoepner L, Oberfield S, Dominguez-Bello MG, Widen EM, et al. Prenatal exposure to antibiotics, cesarean section and risk of childhood obesity. Int J Obes (Lond). (2015) 39:665–70. doi: 10.1038/ijo.2014.180

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Rutayisire E, Wu X, Huang K, Tao S, Chen Y, Tao F. Cesarean section may increase the risk of both overweight and obesity in preschool children. BMC Pregnancy Childbirth. (2016) 16:338. doi: 10.1186/s12884-016-1131-5

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Choo V. WHO reassesses appropriate body-mass index for Asian populations. Lancet. (2002) 360:235. doi: 10.1016/s0140-6736(02)09512-0

PubMed Abstract | CrossRef Full Text | Google Scholar

20. World Health Organization. WHO Child Growth Standards SAS igrowup package. Geneva, Switzerland. (2011). Available online at: https://www.who.int/childgrowth/software/en/ (accessed January 16, 2020).

Google Scholar

21. Mattsson M, Maher GM, Boland F, Fitzgerald AP, Murray DM, Biesma R. Group-based trajectory modelling for BMI trajectories in childhood: a systematic review. Obes Rev. (2019) 20:998–1015. doi: 10.1111/obr.12842

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Jones BL, Nagin DS, Roeder K. A SAS procedure based on mixture models for estimating developmental trajectories. Sociol Methods Res. (2001) 29:374–93. doi: 10.1177/0049124101029003005

CrossRef Full Text | Google Scholar

23. Hu Z, Tylavsky FA, Han JC, Kocak M, Fowke JH, Davis RL, et al. Maternal metabolic factors during pregnancy predict early childhood growth trajectories and obesity risk: the CANDLE study. Int J Obes (Lond). (2019) 43:1914–22. doi: 10.1038/s41366-019-0326-z

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Haga C, Kondo N, Suzuki K, Sato M, Ando D, Yokomichi H, et al. Developmental trajectories of body mass index among Japanese children and impact of maternal factors during pregnancy. PLoS ONE. (2012) 7:e51896. doi: 10.1371/journal.pone.0051896

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Mi YJ, Zhang B, Wang HJ, Yan J, Han W, Zhao J, et al. Prevalence and secular trends in obesity among Chinese adults, 1991–2011. Am J Prev Med. (2015) 49:661–69. doi: 10.1016/j.amepre.2015.05.005

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Li HT, Luo S, Trasande L, Hellerstein S, Kang C, Li JX, et al. Geographic variations and temporal trends in cesarean delivery rates in China, 2008–2014. JAMA. (2017) 317:69–76. doi: 10.1001/jama.2016.18663

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Pei Z, Heinrich J, Fuertes E, Flexeder C, Hoffmann B, Lehmann I, et al. Cesarean delivery and risk of childhood obesity. J Pediatr. (2014) 164:1068–73.e2. doi: 10.1016/j.jpeds.2013.12.044

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Vinding RK, Sejersen TS, Chawes BL, Bonnelykke K, Buhl T, Bisgaard H, et al. Cesarean delivery and body mass index at 6 months and into childhood. Pediatrics. (2017) 139:e20164066. doi: 10.1542/peds.2016-4066

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Li HT, Zhou YB, Liu JM. The impact of cesarean section on offspring overweight and obesity: a systematic review and meta-analysis. Int J Obes (Lond). (2013) 37:893–9. doi: 10.1038/ijo.2012.195

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Cai M, Loy SL, Tan KH, Godfrey KM, Gluckman PD, Chong YS, et al. Association of elective and emergency cesarean delivery with early childhood overweight at 12 months of age. JAMA Netw Open. (2018) 1:e185025. doi: 10.1001/jamanetworkopen.2018.5025

CrossRef Full Text | Google Scholar

31. Peneau S, Giudici KV, Gusto G, Goxe D, Lantieri O, Hercberg S, et al. Growth Trajectories of body mass index during childhood: associated factors and health outcome at adulthood. J Pediatr. (2017) 186:64–71.e1. doi: 10.1016/j.jpeds.2017.02.010

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Indiani C, Rizzardi KF, Castelo PM, Ferraz LFC, Darrieux M, Parisotto TM. Childhood obesity and firmicutes/bacteroidetes ratio in the gut microbiota: a systematic review. Child Obes. (2018) 14:501–09. doi: 10.1089/chi.2018.0040

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Karvonen AM, Sordillo JE, Gold DR, Bacharier LB, O'Connor GT, Zeiger RS, et al. Gut microbiota and overweight in 3-year old children. Int J Obes (Lond). (2019) 43:713–23. doi: 10.1038/s41366-018-0290-z

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Rutayisire E, Huang K, Liu Y, Tao F. The mode of delivery affects the diversity and colonization pattern of the gut microbiota during the first year of infants' life: a systematic review. BMC Gastroenterol. (2016) 16:86. doi: 10.1186/s12876-016-0498-0

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Eny KM, Chen S, Anderson LN, Chen Y, Lebovic G, Pullenayegum E, et al. Breastfeeding duration, maternal body mass index, and birth weight are associated with differences in body mass index growth trajectories in early childhood. Am J Clin Nutr. (2018) 107:584–92. doi: 10.1093/ajcn/nqx081

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Smith AJ, O'Sullivan PB, Beales DJ, de Klerk N, Straker LM. Trajectories of childhood body mass index are associated with adolescent sagittal standing posture. Int J Pediatr Obes. (2011) 6:e97–106. doi: 10.3109/17477166.2010.530664

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Koning M, Hoekstra T, de Jong E, Visscher TL, Seidell JC, Renders CM. Identifying developmental trajectories of body mass index in childhood using latent class growth (mixture) modelling: associations with dietary, sedentary and physical activity behaviors: a longitudinal study. BMC Public Health. (2016) 16:1128. doi: 10.1186/s12889-016-3757-7

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Patel R, Tilling K, Lawlor DA, Howe LD, Hughes RA, Bogdanovich N, et al. Socioeconomic differences in childhood BMI trajectories in Belarus. Int J Obes (Lond). (2018) 42:1651–60. doi: 10.1038/s41366-018-0042-0

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Thompson JA. The risks of advancing parental age on neonatal morbidity and mortality are U- or J-shaped for both maternal and paternal ages. BMC Pediatr. (2020) 20:453. doi: 10.1186/s12887-020-02341-0

CrossRef Full Text | Google Scholar

40. Casteleiro A, Paz-Zulueta M, Parás-Bravo P, Ruiz-Azcona L, Santibañez M. Association between advanced maternal age and maternal and neonatal morbidity: a cross-sectional study on a Spanish population. PLoS ONE. (2019) 14:e0225074. doi: 10.1371/journal.pone.0225074

PubMed Abstract | CrossRef Full Text | Google Scholar