Skip to content

Advertisement

  • Research article
  • Open Access
  • Open Peer Review

Childhood adiposity trajectories are associated with late adolescent blood pressure: birth to twenty cohort

  • 1, 2, 3, 4Email author,
  • 3,
  • 1, 4 and
  • 3
BMC Public HealthBMC series – open, inclusive and trusted201616:665

https://doi.org/10.1186/s12889-016-3337-x

  • Received: 21 April 2016
  • Accepted: 21 July 2016
  • Published:
Open Peer Review reports

Abstract

Background

Elevated blood pressure in childhood is a risk factor for adult hypertension which is a global health problem. Excess adiposity in childhood creates a predisposition to develop adult hypertension. Our aim was to explore distinct sex-specific adiposity trajectories from childhood to late adolescence and examined their association with blood pressure.

Methods

Latent Class Growth Mixture Modeling (LCGMM) on longitudinal data was used to derive sex-specific and distinct body mass index (BMI: kg/m2) trajectories. We studied 1824 black children (boys = 877, girls = 947) from the Birth to Twenty (Bt20) cohort from Soweto, South Africa, and obtained BMI measures at ages 5 through 18 years. Participants with at least two age-point BMI measures, were included in the analysis. Analysis of variance (ANOVA), chi-square test, multivariate linear and standard logistic regressions were used to test study characteristics and different associations.

Results

We identified three (3) and four (4) distinct BMI trajectories in boys and girls, respectively. The overall prevalence of elevated blood pressure (BP) was 34.9 % (39.4 % in boys and 30.38 % in girls). Boys and girls in the early onset obesity or overweight BMI trajectories were more likely to have higher BP values in late adolescence. Compared to those in the normal weight BMI trajectory, girls in early onset obesity trajectories had an increased risk of elevated BP with odds ratio (OR) of 2.18 (95 % confidence interval 1.31 to 4.20) and 1.95 (1.01 to 3.77). We also observed the weak association for boys in early onset overweight trajectory, (p-value = 0.18 and odds ratio of 2.39 (0.67 to 8.57))

Conclusions

Distinct weight trajectories are observed in black South African children from as early as 5 years. Early onset adiposity trajectories are associated with elevated BP in both boys and girls. It is important to consider individual patterns of early-life BMI development, so that intervention strategies can be targeted to at-risk individuals.

Keywords

  • Blood pressure
  • Latent classes
  • Latent class growth mixture modeling
  • Body mass index trajectories
  • Trajectories
  • Hypertension
  • Childhood adiposity
  • Obesity

Background

Hypertension is a global public health problem affecting more than one billion people globally [1]. It is a risk factor for developing cardiovascular disease (CVD) and also contributes to an increase in mortality worldwide. A recent systematic analysis reported that hypertension accounts for about nine million deaths globally every year [2]. In African countries, 25 % of deaths are caused by hypertension [1]. Hypertension is also one of the leading causes of heart attacks, stroke and kidney failure [2]. Results from the Heart of Soweto Study reported a 55 % hypertension prevalence among adult South Africans aged 52.8 years on average [3]. Some studies have looked at elevated blood pressure (BP) in childhood, mostly in rural South Africa, with reported prevalence rates ranging from 1 to 25.9 % [47].

Previous studies have reported that childhood blood pressure weight gain from childhood to adulthood, adulthood obesity [8, 9], childhood and adolescent physical in-activity and certain lifestyle behaviors such as use of tobacco and alcohol consumption [10], are some of the determinant factors of adulthood hypertension. Understanding the life-course progression of adiposity in children is important since childhood adiposity is associated with adult obesity which has been reported to be linked to increased hypertension risk in adults [11]. Kagura and colleagues reported that elevated BP can be observed from childhood in black urban South African children [12]. Most of the studies that reported on childhood obesity and adult hypertension have used cross sectional data. There are few studies that have used longitudinal data to understand the effect of life-course childhood adiposity on late adolescent blood pressure and lack of comprehensive longitudinal data have made it difficult to do such studies in an African setting.

There has been increased interest in the role of growth patterns on BP, it is also still unclear how early childhood adiposity influences later BP [13]. Using longitudinal data will help elucidate this relationship and this study in particular will add further evidence from an African perspective. For the first time in an exclusively African population, we sought to determine whether adiposity progression (from 5 to 18 years of age), focused on using body mass index (BMI) as an adiposity marker, could be used to predict BP among late adolescents.

We use a relatively new method in both biology and epidemiology called Latent Class Growth Mixture Modeling (LCGMM) to identify distinct sex-specific adiposity trajectories. LCGMM groups individuals with the same developmental trajectory in the same class. It also allows for inclusion of variation in several growth parameters both within and between classes [14]. Developmental trajectories of BMI describe individuals’ BMI change over time. One of the advantages for using LCGMM is that adiposity trajectories are identified independently without pre-assumptions of the relative contributions of various lifestyles, genetic and epigenetic factors. Therefore multiple factors may contribute to variation in the information of these adiposity trajectories.

The aims of this study were to: 1) determine the distinct sex-specific patterns of adiposity trajectories in black South Africa children from 5 to 18 years of age, 2) find the prevalence of elevated blood pressure in late adolescence, and 3) to explore associations between these distinct adiposity trajectories to elevated BP in late adolescence.

Methods

Participants

To identify developmental patterns of BMI from 5 to 18 years old and relate them to blood pressure in late adolescence, data from the Birth to Twenty cohort (Bt20) were used. Bt20 is Africa’s largest and longest running longitudinal birth cohort, with 3273 children at time of enrolment. It is focused on the health and development of children born in a South African urban township, namely Soweto in Johannesburg. The cohort comprised of 1682 girls and 1591 boys of which the majority were black children (78.4 %). The participants have been followed up to date using different means of communication. The detailed cohort profile with recruitment and cohort attrition details has been described elsewhere [15]. Only black participants (n = 1824) with weight and height data available for at least two time points between 5 and 18 years were included in this study.

Measures

Anthropometrics

Trained research assistants collected anthropometric measurements. Weight was measured using a digital scale to the nearest 0.1 kg with participants in light clothes and without shoes. A wall-mounted stadiometer (Holtain, UK) was used to measure standing height to the nearest 0.1 cm. Weight and height at 5 years and 7–18 years old were used to calculate BMI (weight (kg)/height (m2)), used as a marker for adiposity at corresponding years.

Blood pressure assessment

Blood pressure (mm Hg) was measured using Omron 6 automated machine (Kyoto, Japan) from 8 years of age onwards and a Dinamap Vital Signs monitor 1846SX (Critikon, USA) was used at 5 years of age. At each assessment, participants’ seated blood pressure was measured three times with a 2 min interval between each measurement. The blood pressure measurements were taken after 5 min of seated rest. The mean average for the second and third right arm readings was recorded for the current analysis. Blood pressure measurements from one visit have been used in blood pressure studies before [12, 13, 16, 17].

The mean arterial pressure (MAP) was calculated from systolic blood pressure (SBP) and diastolic blood pressure (DBP) using the formula; MAP = (SBP + (2 * DBP))/3. We followed a standard procedure as recommended by the fourth National High Blood Pressure Education Program working group on high blood pressure in children and adolescents (NHBPEP) [18] report to classify blood pressure either normotensive (BP less than 90th percentile), prehypertension or hypertension (BP equal or over 90th percentile for sex, age and height) for those participants who were above 17 years but not yet 18 years old the time BP data was collected. For participants who were already 18 years or older we used the cut-offs as recommended in the seventh report of the Joint National Committee on Prevention, Detection, Evaluation, and Treatment of High Blood Pressure [19]. Prehypertension was defined as SBP readings from 120 to 139 mm Hg or a DBP from 80 to 89 mm Hg while hypertensive was defined as SBP readings equal or greater than 140 mm Hg or DBP readings equal or greater than 90 mm Hg. Due to small sample size the prevalence of hypertension, prehypertension and hypertension were combined and defined as elevated blood pressure.

Statistical analysis

Different group based modeling methods such as Latent Class Growth Analysis (LCGA) and Latent Class Growth Mixture Modeling (LCGMM), have been used to identify and classify individuals into different trajectories. LCGA technique was developed by Nagin and colleagues [2023] and is implemented through SAS Proc Traj in the SAS software. There is also a Traj Stata software plugin developed by Jones and Nagin, 2012 [24]. On the other hand LCGMM was developed by Muthén and colleagues and is implemented in Mplus [2528]. We used Mplus and LCGMM since it allows variation in the intercept and slope in both within the class (inter-individual differences) and across classes while SAS Proc Traj allows only across classes variation. Using LCGMM adds more heterogeneity in the model, thus being the more flexible method as such was preferred in the current analysis [14, 22, 25, 26, 2830].

The exploratory analyses to check the patterns of missing data concluded that data was missing at random. Mplus handles missing data by the expectation-maximization algorithm (EM algorithm) with assumption of missing at random (MAR) [28]. This makes use of full information maximum likelihood estimation with missing values (FIML) by using estimation to integrate all available information based on MAR assumptions. This means that we can make full use of all available data in our analysis. This prevents the inclusion of only those participants who have no missing data at all data points, which would subsequently reduce the sample size. So the use of the maximum likelihood (ML) approach implemented in Mplus overcomes potential biases in participants with substantial missing data [25, 26, 28, 31, 32].

We included only those participants with at least two time points of available data for BMI in the current analysis. Those with one time point of available data were excluded as this may have influenced the identification and pattern of the BMI trajectories. On average, seven available data points per participant were used to perform this longitudinal analysis.

Latent Class Growth mixture modeling (LCGMM in Mplus version 7.3 [28] was performed to identify the distinct sex-specific BMI trajectories between 5 and 18 years.

The BMI values were used in this study because they are more sensitive to changes in body composition and have been recommended over BMI z-scores [3335].

Once BMI trajectory group membership was determined for both girls and boys, each respondent was assigned to the class according to the highest probability of belonging to that class. The remainder of analyses was conducted using Stata version 13 [36] to examine additional descriptive characteristics of the classes. We used both standard logistic and multivariate linear regressions to estimate the relationship between latent class membership and blood pressure in late adolescence. Analysis of variance (ANOVA) and chi-square test results were used to assess the differences in different study characteristics among the BMI trajectory classes at 5 % level of significance.

Results

BMI Trajectories: optimal number of BMI trajectories

Using LCGMM, we explored linear, quadratic and cubic slopes to fit the model. The cubic slope coefficients were either not significant or required prohibitive running times and did not converge in all models. This suggests that the cubic time function model did not fit the data well. After applying the criteria as recommended by other authors [14, 25, 26, 29, 30, 3739], only quadratic models involving quadratic time function were further analyzed.

To determine the optimal number of latent classes, we used different model fit indices in conjunction with other criteria as used in previous studies [14, 30, 37, 40, 41]. Firstly, we looked at the three model statistics, the Bayesian Information Criterion (BIC) [42], the Bootstrap Likelihood Ratio Test (BLRT) [43] and Lo, Mendell and Rubin Likelihood Ratio Test (LMR-LRT) [44]. Some studies have applied BIC in variable modeling analyses [41, 42, 45, 46]. It estimates a model to be true using posterior probabilities. A lower BIC would indicate that a model is more likely to be considered as the true model [42]. BIC reduces the false positive rate, hence most individuals will be assigned in their right BMI latent class [47]. We also looked at the entropy value and number of study participants per class. We employed a less restrictive 1 % group membership in this study, which has also been used before in other studies [41, 48]. Lastly we looked at the shape of the BMI trajectories and its clinical interpretability.

Model fit statistics used to come up with optimal number of classes in both boys and girls are shown in Table 1. Applying the criterion discussed above a four-class model fits our data well in girls and a three-class solution has the best classification of BMI trajectory membership in boys.
Table 1

Model fit statistics for quadratic LCGMM comparing solutions for 1 to 7 latent classes

Sex

Latent classes

BICa

LMR-LRTb (P-value)

BLRTc (P-value)

Convergence

Entropy value

Log Likelihood

Number of parameters

Boys

1

23605.4

N/A

N/A

Yes

N/A

−11738.3

19

2

22846.6

0.03

0.00

Yes

0.98

−11345.4

23

3 d

22644.4

0.13

0.67

Yes

0.99

−11230.7

27

4

22544.7

0.53

1.00

Yes

0.98

−11167.3

31

5

22476.5

0.77

0.67

Yes

0.93

−11119.7

35

6

22438.4

0.61

1.00

Yes

0.94

−11087.1

39

7

   

NO

  

43

Girls

1

27949.8

N/A

N/A

Yes

N/A

−13909.8

19

2

27226.2

0.00

0.00

Yes

0.94

−13534.3

23

3

27155.0

0.10

0.00

Yes

0.76

−13485.0

27

4 d

27083.9

0.23

0.00

Yes

0.80

−13435.7

31

5

27026.0

0.45

0.6

Yes

0.82

−13393.1

35

6

26990.9

0.63

1.00

Yes

0.78

−13361.8

39

7

26961.2

0.27

1.00

Yes

0.80

−13333.3

43

LCGMM latent class growth mixture modeling

aBayesian Information Criterion

bLo, Mendell and Rubin Likelihood Ratio Test

cBootstrapped Likelihood Ratio Test, N/A – not applicable

dThe optimal class solutions according to the model fit criteria are shown in bold

In order to facilitate interpretation of the results, Fig. 1a and b show the BMI trajectory classes for 947 girls and 877 boys in relation to the extended International Obesity Task Force (IOFT) [49]. For girls, four BMI trajectory classes can be described as: Trajectory 1 (normal weight: 75.7 %), Trajectory 2 (late onset overweight: 15 %), Trajectory 3 (early onset obesity to overweight: 4.8 %) and lastly Trajectory 4 (early onset obesity to morbidly obese: 4.5 %). For boys, a three-class model fits the data well. These classes can be classified as: Trajectory 1 (normal weight: 92.7 %), Trajectory 2 (early onset overweight to normal: 6 %) and Trajectory 3 (the early onset overweight to obese trajectory: 1.3 %).
Fig. 1
Fig. 1

a BMI Trajectories in girls. BMI latent classes plotted together with Extended International Obesity Task Force (IOTF) Cut-Offs for BMI in girls. b BMI Trajectories in boys. BMI latent classes plotted together with Extended International Obesity Task Force (IOTF) Cut-Offs for BMI in boys

Study characteristics by BMI trajectory classes are shown in Tables 2 and 3. Most of the study characteristics were significantly different among BMI trajectories in girls. In boys, lower maternal education and SBP at 18 years were significantly different among the BMI trajectories and the rest were not.
Table 2

Study characteristics by BMI trajectory in girls: The Birth to Twenty Cohort

 

Normal weight (n = 717)

Trajectory 2 (n = 142)

Trajectory 3 (n = 45)

Trajectory 4 (n = 43)

P value

Study characteristics

    

 Birth weight (kg)

3.0 (0.5)

3.0 (0.5)

3.0 (0.5)

3.3 (0.5)

0.01*

 Height at 18 years (cm)

159.7 (6.2)

158.4 (5.1)

157.3 (6.6)

159.3 (6.1)

0.03*

 SBP at 18 years (mm Hg)

114.3 (9.6)

116.0 (9.5)

117.9 (11.7)

120.1 (11.3)

0.00***

 DBP at 18 years (mm Hg)

71.5 (8.3)

72.8 (8.1)

75.8 (12.7)

75.8 (9.0)

0.00***

 MAP at 18 years (mm Hg)

85.8 (7.9)

87.2 (7.8)

89.8 (11.8)

90.6 (8.9)

0.00***

 Elevated BP n (%)

187 (28.5)

28 (31.1)

17 (43.6)

18 (46.1)

0.03*

For girls: Trajectory 2 = Late Onset Overweight, Trajectory 3 = Early Onset Obese to Overweight, Trajectory 4 = Early Onset Obese to Morbidity Obese

Summary statistics presented as mean (±standard deviation) otherwise stated

*P < 0.05; **P < 0.01; ***P < 0.001

ANOVA and chi square test were conducted to describe the study characteristics differences in BMI trajectories

Percentages may not add up to 100 % due to rounding

Table 3

Study characteristics by BMI trajectory in boys: the birth to twenty cohort

 

Normal weight (n = 813)

Trajectory 2 (n = 53)

Trajectory 3 (n = 11)

P value

Study characteristics

   

 Birth weight (kg)

3.1 (0.5)

3.2 (0.6)

3.2 (0 .7)

0.46

 Height at 18 years (cm)

170.6 (6.7)

172.4 (5.5)

170.4 (7.6)

0.21

 SBP at 18 years (mm Hg)

121.0 (10.8)

124.2 (8.2)

134.5 (16.9)

0.00***

 DBP at 18 years (mm Hg)

71.0 (8.6)

71.9 (8.3)

76.8 (8.7)

0.09

 MAP at 18 years (mm Hg)

87.7 (8.1)

89.4 (7.3)

96.0 (10.3)

0.00**

 Elevated BP n (%)

276 (38.7)

20 (46.5)

6 (60.0)

0.24

For boys: Trajectory 2 = Early Onset Overweight to Normal Weight, Trajectory 3 = Early Onset Overweight to Obese

Summary statistics presented as mean (±standard deviation) otherwise stated

*P < 0.05; **P < 0.01; ***P < 0.001

ANOVA and chi square test were conducted to describe the study characteristics differences in BMI trajectories

Percentages may not add up to 100 % due to rounding

Blood pressure in late adolescence and BMI trajectory classes

Only participants with blood pressure values at 18 years of age were included in this analysis. It included 766 boys and 823 girls. The overall prevalence of elevated blood pressure in late adolescence was 34.91 % (39.43 % boys and 30.38 % girls). Birth weight and height at 18 years were added in the model to determine whether BMI trajectory classes have an effect on SBP, DBP, MAP and elevated BP status as shown in Tables 4 and 5. In girls, Trajectory 3 was significantly associated with SBP, DBP and MAP with mean BP (95 % confidence interval, 95 % CI) of 4.18 mm Hg (1.03–7.33), 4.44 mm Hg (1.65 to 7.22) and 4.35 mm Hg (1.70 to 7.01) compared to normal trajectory, respectively. Trajectory 4 was significantly associated with SBP, DBP and MAP with mean BP of 6.08 mm Hg (2.93 to 9.23), 4.28 mm Hg (1.50 to 7.07) and 4.88 mm Hg (2.23 to 7.54) compared to normal weight trajectory, accordingly. Girls in Trajectories 3 and 4 had a 1.95 fold (1.01 to 3.77) and a 2.18 fold (1.31 to 4.20) increased risk of elevated BP in late adolescence, separately.
Table 4

The association of BMI trajectories, birth weight and height and blood pressure at 18 years for girls

 

Unadjusted

Adjusted

SBP (mm Hg)

  

 BMI Trajectory

  

 Normal Weight

Reference

Reference

 Late Onset Overweight

1.77 (−0.38 to 3.93)a

2.07 (−0.08 to 4.21)

 Early Onset Obese to Overweight

3.66 (0.50 to 6.82)*

4.18 (1.03 to 7.33)**

 Early Onset Obese to Morbidity Obese

5.87 (2.71 to 9.03)**

6.08 (2.93 to 9.23)**

 Birth Weight (kg)

 

−0.61 (−1.99 to 0.78)

 Height (cm) at 18 years

 

0.22 (0.11 to 0.34)**

 R2

0.02

0.04

DBP (mm Hg)

  

 BMI Trajectory

  

 Normal Weight

Reference

Reference

 Late Onset Overweight

1.33 (−0.57 to 3.22)

1.40 (−0.50 to 3.30)

 Early Onset Obese to Overweight

4.29 (1.52 to 7.07)**

4.44 (1.65 to 7.22)**

 Early Onset Obese to Morbidity Obese

4.25 (1.48 to 7.03)**

4.28 (1.50 to 7.07)**

 Birth Weight (kg)

 

−0.05 (−1.27 to 1.18)

 Height (cm) at 18 years

 

0.06 (−0.04 to 0.16)

 R2

0.02

0.02

MAP (mm Hg)

  

 BMI Trajectory

  

 Normal Weight

Reference

Reference

 Late Onset Overweight

1.47 (−0.33 to 3.28)a

1.62 (−0.19 to 3.43)

 Early Onset Obese to Overweight

4.08 (1.43 to 6.73)***

4.35 (1.70 to 7.01)***

 Early Onset Obese to Morbidity Obese

4.79 (2.14 to 7.44)***

4.88 (2.23 to 7.54)***

 Birth Weight (kg)

 

−0.23 (−1.40 to 0.93)

 Height (cm) at 18 years

 

0.12 (0.02 to 0.21)*

 R2

0.03

0.03

Elevated BP at 18 years (Normal BP-reference)

OR

OR

 BMI Trajectory

  

 Normal Weight

Reference

Reference

 Late Onset Overweight

1.13 (0.70 to 1.82)

1.14 (0.71 to 1.84)

 Early Onset Obese to Overweight

1.93 (1.00 to 3.72)*

1.95 (1.01 to 3.77)*

 Early Onset Obese to Morbidity Obese

2.15 (1.11 to 4.12)*

2.18 (1.31 to 4.20)*

 Birth Weight (kg)

 

0.93 (0.67 to 1.27)

 Height (cm) at 18 years

 

1.01 (0.98 to 1.03)

 Pseudo R2

0.01

0.01

SBP systolic blood pressure, DBP diastolic blood pressure, MAP mean arterial pressure, OR odds ratio

Adjusted for birth weight and height at 18 years

*P < 0.05; **P < 0.01; ***P < 0.001

aIntercept (MAP), OR (elevated BP) and 95 % Confidence Interval are presented

Table 5

The association of BMI trajectories, birth weight and height and blood pressure at 18 years for boys

 

Unadjusted

Adjusted

SBP (mm Hg)

  

 BMI Trajectory

  

 Normal Weight

Reference

Reference

 Early Onset Overweight to Normal Weight

3.23 (−0.08 to 6.54)a

2.85 (−0.43 to 6.13)

 Early Onset Overweight to Obese

13.46 (6.75 to 20.17)**

13.45 (6.82 to 20.07)**

 Birth Weight (kg)

 

−1.35 (−2.85 to 0.14)

 Height (cm) at 18 years

 

0.27 (0.15 to 0.38)**

 R2

0.02

0.05

DBP (mm Hg)

  

 BMI Trajectory

  

 Normal Weight

Reference

Reference

 Early Onset Overweight to Normal Weight

0.95 (−1.69 to 3.59)

0.844 (−1.80 to 3.49)

 Early Onset Overweight to Obese

5.77 (0.41 to 11.13)

5.74 (0.39 to 11.10)*

 Birth Weight (kg)

 

−0.53 (−1.7 to 0.68)

 Height (cm) at 18 years

 

0.07 (−0.2 to 0.17)

 R2

0.01

0.01

MAP (mm Hg)

  

 BMI Trajectory

  

 Normal Weight

Reference

Reference

 Early Onset Overweight to Normal Weight

1.71 (−0.77 to 4.19)a

1.51 (−0.96 to 3.99)

 Early Onset Overweight to Obese

8.33 (3.3 to 13.37)***

8.31 (3.30 to 13.32)***

 Birth Weight (kg)

 

−0.80 (−1.93 to 0.33)

 Height (cm) at 18 years

 

0.14 (0.05 to 0.22)**

 R2

0.02

0.03

Elevated BP at 18 years (Normal BP-reference)

OR

OR

 BMI Trajectory

  

 Normal Weight

Reference

Reference

 Early Onset Overweight to Normal Weight

1.38 (0.74 to 2.55)

1.34(0.72 to 2.50)

 Early Onset Overweight to Obese

2.38 (0.66 to 8.49)

2.39 (0.67 to 8.57)

 Birth Weight (kg)

 

0.89 (0.66 to 1.18)

 Height (cm) at 18 years

 

1.02 (1.00 to 1.05)*

 Pseudo R2

0.003

0.01

SBP systolic blood pressure, DBP diastolic blood pressure, MAP mean arterial pressure; OR odds ratio

Adjusted for birth weight and height at 18 years

*P < 0.05; **P < 0.01; ***P < 0.001

aIntercept (MAP), OR (elevated BP) and 95 % Confidence Interval are presented

For boys, Trajectory 3 was significantly associated with SBP, DBP and MAP with mean elevated BP of 13.45 mm Hg (6.82 to 20.07), 5.74 mm Hg (0.30 to 11.10) and 8.31 mm Hg (3.30–13.32) compared to normal weight trajectory respectively. In boys we observed a weak association for Trajectory 3 with a 2.39 fold (0.67 to 8.57) increased risk of elevated BP in late adolescence. The main results did not vary adjusting by birth weight, height at 18 years or both in girls. In boys, only the association between Trajectory 3 and DBP varied after adjusting by both weight and height at 18 years. The association was not significant in the unadjusted model while significant in the adjusted model.

Discussion

There has been a growing interest in studying the heterogeneity in adiposity (BMI) developmental patterns in different populations over past years. Taking advantage of data-driven methods such as LCGMM we explored the distinct developmental patterns of BMI across childhood to late adolescence and its association with late adolescent blood pressure in black South African children. Our results confirm that there is heterogeneity in BMI trajectories in the study sample and that trajectories vary between boys and girls. A four-class BMI trajectory model may best represent heterogeneity in BMI developmental patterns in girls (normal weight, late onset overweight, early onset obesity to overweight and early onset obesity to morbidly obese) and a three-class model among boys (normal weight, early onset overweight to normal and early onset overweight to obese trajectory). These results are consistent with those reported by Hejazi and colleagues in Canadian children aged 2–8 years but with different age group and ancestry population [50]. Similarly, a study by Haga et al. reported that BMI trajectories varied between boys and girls in Japanese children aged between 0 and 12 years [51]. Ventura et al. [46] reported four distinct BMI trajectories in girls only. Few studies have focused on gender differences in BMI trajectories. American, Canadian, Australian and Dutch studies in children reported either three or four-class distinct BMI trajectories in both girls and boys [29, 41, 45, 5257]. Individuals in trajectories characterized by high BMI trends had an increased risk of developing hypertension in late adolescence compared to those in a normal weight trajectory in the current study.

Prevalence of elevated BP in late adolescence was relatively higher in boys (39.43 %) than in girls (30.38 %). This prevalence is consistent with a recent study in urban South African children that reported the prevalence of prehypertension or hypertension ranging between 17.6 and 40.8 % [12]. Another study in youth aged 15 to 24 years in rural Ghana reported levels of prehypertension or hypertension with a prevalence of 37.4 %. The prevalence of prehypertension or hypertension (48.8 vs. 26.6 %) was higher in males than in females [17]. It should be pointed out that the prevalence in the current study was also higher than that previously reported in some studies within and outside South Africa. A study done by Makgae and colleagues in rural South African children aged 6 to 13 years found that 1.0 to 5.8 % of boys and 3.1 to 11.4 % of girls were hypertensive [58]. Moselakgomo et al. assessed BMI and blood pressure among adolescent school children aged 10 to16 years in Limpopo province, South Africa and the results revealed that prevalence of hypertension was between 2.3 and 5.9 % [59]. The following blood pressure prevalence has been reported in other African countries. In Brazzaville, Congo 24.3 % in girls and 16.6 % prevalence of hypertension in boys school children aged 5 to 18 years [60]. In Egyptian adolescents aged 11–19 years revealed that the prevalence rates of prehypertension and hypertension were 5.7 and 4.0 %, respectively [16]. The mechanisms of the observed sex differences in elevated BP in our study sample are not clear but may be attributed to sexual dimorphism as observed in other BP control programs [61] and also due to cardiovascular fetal programming as reported before [61, 62].

Girls in the early onset obesity trajectories were likely to have elevated SBP, DBP and MAP in late adolescence compared to those in a healthy trajectory. They displayed a 118 % increased odds of developing hypertension in late adolescence compared to those in the normal weight trajectory. Boys in early onset overweight to obese trajectory were likely to have a higher SBP, DBP and MAP in late adolescence. For instance boys in this trajectory had a 13 mm Hg mean SBP higher than those in a normal weight trajectory.

One longitudinal study on life-course adiposity and blood pressure in late adolescents has been done in Australian children [63], three of the seven adiposity trajectories were associated with elevated blood pressure at 17 years old of age. We report that being in an early onset obese or overweight trajectory was associated with increased risk of elevated BP in both girls and boys which is in agreement with other studies, although comparative studies were performed in a different population and different age range [63]. Girls in the early onset obesity to overweight BMI trajectory had reduced risk of elevated BP in late adolescence compared to those who were in early onset obesity to morbidity obese trajectory despite the fact that both groups had early onset obesity. This study was not an intervention study but it is important to study the mechanisms behind the weight loss in this group of girls. This could be important in reducing hypertension risk later in life.

The current study has several strengths; firstly, the longitudinal study in black South African children is a representative sample of South African children hence more relevant in understanding the different BMI trajectories in black South African children. Using LCGMM in identifying trajectories is an ideal analytical method since it involves a number of criteria in selecting the best fitting model and it predicts individual class membership from the available data and its ability to deal with missing at random data. The current results therefore clearly show individual heterogeneity in BMI development in boys and girls from five years to late adolescence. Our analysis was stratified according to sex and performed in a South African black population, which gives us an in-depth understanding of the difference in BMI developmental patterns in black boys and girls in South Africa. These results suggest that targeted intervention might be developed for individuals in high-risk trajectories. Limitations of our current study include the fact that the sample used in this study is from Soweto, South Africa and generalizability to other parts of Africa should be treated with caution. Adiposity trajectories explain only 5 % variation in late adolescent blood pressure and we do not know yet what the other factors are, that would explain the remaining variance. We speculate that the role of genetics and epigenetics could be important, as previously reported [64]. The role of genetics and epigenetics on adiposity and growth patterns in African populations is understudied and it calls for further investigation. A recent study by Kagura and colleagues in this cohort reported no association between blood pressure at 18 and alcohol consumption or smoking or both at 18 years in this population [65] and thus we did not include these variables within the models. Furthermore, we examined walking to school as a proxy of physical activity but this variable was not able to different among participants given that most participants are walking extensively during the day to and from school.

Due to the limited sample size we might not be able to capture all trajectories; this also influenced the observation that the high-risk trajectory in boys comprised of very few individuals that might influence association analysis.

Conclusions

In this study, we identified distinct sex-specific trajectories. The early onset obesity or overweight trajectories are associated with elevated blood pressure in late adolescence. These results signify the importance to consider patterns of BMI development, especially at early stage of development, so that prevention strategies may be implemented to target those individuals in high-risk developmental patterns.

Our study has shown that patterns of adiposity could be a preferred predictor of future SBP, DBP, MAP and elevated BP in late adolescents compared to cross-sectional BMI measures. Furthermore, identification of childhood obesity can help with early identification of those at risk of developing elevated BP and other chronic diseases in adulthood.

Abbreviations

BMI, body mass index; DBP, diastolic blood pressure; LCGA: latent class growth analysis; LCGMM: latent class growth mixture modeling; MAP, mean arterial pressure; SBP, systolic blood pressure

Declarations

Acknowledgments

We wish to thank the Bt20 participants for taking part in the study and the Bt20 team for their relentless support throughout the study.

Funding

The Bt20 is supported by the Wellcome Trust (UK) [grant number 092097/Z/10/Z], the South African Medical Research Council, the University of Witwatersrand and DST-NRF Centre of Excellence in Human Development, University of the Witwatersrand, Johannesburg, South Africa funding for the collection, cleaning and storage and management of the data. This work is based on the research supported in part by the Thuthuka Programme of the National Research Foundation of South Africa [Award No. 94007]. The financial assistance of the National Research Foundation (NRF) towards this research is hereby acknowledged. Opinions expressed and conclusions arrived at, are those of the author and are not necessarily to be attributed to the NRF. Norris is supported by the UK MRC/DFID African Research Leader Scheme.

Availability of data and materials

Data on this study is available on request from the Developmental Pathways for Health Research Unit data management department by contacting Prof. Shane A Norris.

Authors’ contributions

RM, JK, ZL and SAN carried out the initial conceptualization and design of the study. RM led the study, conducted the initial analyses and interpreted the data, drafted the initial manuscript, and approved the final manuscript as submitted in close collaboration with JK, ZL and SAN who gave guidance throughout the study from conception to the approval of the final draft, reviewed and edited the data analysis plan and the manuscript for intellectual content. All authors read and approved the final manuscript as submitted.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Prior to the study, ethical approval was obtained from the Human Research Ethics Committee of the University of the Witwatersrand, Johannesburg (Reference number M0101556), participants or their parents/caregivers when the participants were minor gave informed consent at the beginning of each data collection session throughout the study.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

(1)
School of Molecular and Cell Biology, Faculty of Science, University of the Witwatersrand, The Mount, 9 Jubilee Road, Parktown, Johannesburg, 2193, South Africa
(2)
Sydney Brenner Institute for Molecular Biosciences (SBIMB), University of the Witwatersrand, Johannesburg, South Africa
(3)
MRC/Wits Developmental Pathways for Health Research Unit (DPHRU), University of the Witwatersrand, Johannesburg, South Africa
(4)
Division of Human Genetics, School of Pathology, Faculty of Health Sciences, University of the Witwatersrand and National Health Laboratory Service, Johannesburg, South Africa

References

  1. Alwan A. Global status report oSn noncommunicable diseases 2010. Geneva: World Health Organization; 2011.Google Scholar
  2. Lim SS, Vos T, Flaxman AD, Danaei G, Shibuya K, Adair-Rohani H, AlMazroa MA, Amann M, Anderson HR, Andrews KG. A comparative risk assessment of burden of disease and injury attributable to 67 risk factors and risk factor clusters in 21 regions, 1990–2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2013;380(9859):2224–60.View ArticleGoogle Scholar
  3. Sliwa K, Wilkinson D, Hansen C, Ntyintyane L, Tibazarwa K, Becker A, Stewart S. Spectrum of heart disease and risk factors in a black urban population in South Africa (the Heart of Soweto Study): a cohort study. Lancet. 2008;371(9616):915–22.View ArticlePubMedGoogle Scholar
  4. Daniel G, Amusa L, Mhlongo D, Khoza L, Anyanwu F. Elevated blood pressure among rural South African children in Thohoyandou, South Africa. Iran J Public Health. 2013;42(5):489.Google Scholar
  5. Kemp C, Pienaar AE, Schutte AE. The prevalence of hypertension and the relationship with body composition in Grade 1 learners in the North West Province of South Africa. S Afr J Sports Med. 2011;23(4):122–77.View ArticleGoogle Scholar
  6. Monyeki K, Kemper H, Makgae P. The association of fat patterning with blood pressure in rural South African children: the Ellisras Longitudinal Growth and Health Study. Int J Epidemiol. 2006;35(1):114–20.View ArticlePubMedGoogle Scholar
  7. Schutte A, Van Rooyen J, Huisman H, Kruger H, Malan N, De Ridder J. Dietary risk markers that contribute to the aetiology of hypertension in black South African children: the THUSA BANA study. J Hum Hypertens. 2003;17(1):29–35.View ArticlePubMedGoogle Scholar
  8. Bridger T. Childhood obesity and cardiovascular disease. Paediatr Child Health. 2009;14(3):177.PubMedPubMed CentralGoogle Scholar
  9. Christofaro DG, Ritti-Dias RM, Fernandes RA, Polito MD. Andrade SMd, Cardoso JR, Oliveira ARd: High blood pressure detection in adolescents by clustering overall and abdominal adiposity markers. Arq Bras Cardiol. 2011;96(6):465–70.View ArticlePubMedGoogle Scholar
  10. Petkeviciene J, Klumbiene J, Simonyte S, Ceponiene I, Jureniene K, Kriaucioniene V, Raskiliene A, Smalinskiene A, Lesauskaite V. Physical, Behavioural and Genetic Predictors of Adult Hypertension: The Findings of the Kaunas Cardiovascular Risk Cohort Study. 2014.Google Scholar
  11. Singh AS, Mulder C, Twisk JW, Van Mechelen W, Chinapaw MJ. Tracking of childhood overweight into adulthood: a systematic review of the literature. Obes Rev. 2008;9(5):474–88.View ArticlePubMedGoogle Scholar
  12. Kagura J, Adair LS, Musa MG, Pettifor JM, Norris SA. Blood pressure tracking in urban black South African children: birth to twenty cohort. BMC Pediatr. 2015;15(1):78.View ArticlePubMedPubMed CentralGoogle Scholar
  13. Jones A, Charakida M, Falaschetti E, Hingorani AD, Finer N, Masi S, Donald AE, Lawlor DA, Smith GD, Deanfield JE. Adipose and height growth through childhood and blood pressure status in a large prospective cohort study. Hypertension. 2012;59(5):919–25.View ArticlePubMedPubMed CentralGoogle Scholar
  14. Berlin KS, Parra GR, Williams NA. An introduction to latent variable mixture modeling (part 2): longitudinal latent class growth analysis and growth mixture models. J Pediatr Psychol. 2014;39(2):188–203.View ArticlePubMedGoogle Scholar
  15. Richter L, Norris S, Pettifor J, Yach D, Cameron N. Cohort profile: Mandela's children: the 1990 Birth to Twenty study in South Africa. Int J Epidemiol. 2007;36(3):504–11.View ArticlePubMedPubMed CentralGoogle Scholar
  16. Abolfotouh MA, Sallam SA, Mohammed MS, Loutfy AA, Hasab AA. Prevalence of elevated blood pressure and association with obesity in Egyptian school adolescents. Int J Hypertens. 2011;2011:8.Google Scholar
  17. Afrifa–Anane E, Agyemang C, Codjoe SNA, Ogedegbe G, Aikins A-G. The association of physical activity, body mass index and the blood pressure levels among urban poor youth in Accra, Ghana. BMC Public Health. 2015;15(1):1.View ArticleGoogle Scholar
  18. Health UD, Services H. The fourth report on the diagnosis, evaluation, and treatment of high blood pressure in children and adolescents. 2005. 2012.Google Scholar
  19. Chobanian AV, Joint National Committee on Prevention, Detection, Evaluation, and Treatment of High Blood Pressure. National Heart, Lung, and Blood Institute; National High Blood Pressure Education Program Coordinating Committee: Seventh report of the Joint National Committee on Prevention, Detection, Evaluation, and Treatment of High Blood Pressure. Hypertension. 2003;42:1206–52.View ArticlePubMedGoogle Scholar
  20. Nagin D. Group-based modeling of development. Massachusetts: Harvard University Press; 2005.View ArticleGoogle Scholar
  21. Hennig BJ, Fulford AJ, Sirugo G, Rayco-Solon P, Hattersley AT, Frayling TM, Prentice AM. FTO gene variation and measures of body mass in an African population. BMC Med Genet. 2009;10:21.View ArticlePubMedPubMed CentralGoogle Scholar
  22. Kelly T, Yang W, Chen CS, Reynolds K, He J. Global burden of obesity in 2005 and projections to 2030. Int J Obes (Lond). 2008;32(9):1431–7.View ArticleGoogle Scholar
  23. Puoane T, Steyn K, Bradshaw D, Laubscher R, Fourie J, Lambert V, Mbananga N. Obesity in South Africa: The South African Demographic and Health Survey. Obes Res. 2002;10(10):1038–48.View ArticlePubMedGoogle Scholar
  24. Rossouw HA, Grant CC, Viljoen M. Overweight and obesity in children and adolescents: The South African problem. S Afr J Sci. 2012;108(5/6):7.View ArticleGoogle Scholar
  25. Muthen B. Latent variable mixture modeling. In: New developments and techniques in structural equation modeling. 2001. p. 1–33.Google Scholar
  26. Muthén B. The potential of growth mixture modelling. Infant Child Dev. 2006;15(6):623–5.View ArticleGoogle Scholar
  27. Mei H, Chen W, Jiang F, He J, Srinivasan S, Smith EN, Schork N, Murray S, Berenson GS. Longitudinal Replication Studies of GWAS Risk SNPs Influencing Body Mass Index over the Course of Childhood and Adulthood. PLoS One. 2012;7(2):e31470.View ArticlePubMedPubMed CentralGoogle Scholar
  28. Muthén LKaM, B.O. Mplus User's Guide. 1998–2012, Seventh Edition.Google Scholar
  29. Hoekstra T, Barbosa-Leiker C, Koppes LL, Twisk JW. Developmental trajectories of body mass index throughout the life course: an application of Latent Class Growth (Mixture) Modelling. Longitudinal Life Course Stud. 2011;2(3):319–30.Google Scholar
  30. Jung T, Wickrama KAS. An Introduction to Latent Class Growth Analysis and Growth Mixture Modeling. Soc Personal Psychol Compass. 2008;2(1):302–17.View ArticleGoogle Scholar
  31. Steyn K, Fourie JM. Media Data Document Cape Town: Department of Medicine, University of Cape Town, and Chronic Diseases of Lifestyle Unit. 2007.Google Scholar
  32. Enders CK, Bandalos DL. The relative performance of full information maximum likelihood estimation for missing data in structural equation models. Struct Equ Model. 2001;8(3):430–57.View ArticleGoogle Scholar
  33. Berkey CS, Colditz GA. Adiposity in adolescents: change in actual BMI works better than change in BMI z score for longitudinal studies. Ann Epidemiol. 2007;17(1):44–50.View ArticlePubMedGoogle Scholar
  34. Boyer BP, Nelson JA, Holub SC. Childhood body mass Index trajectories predicting cardiovascular risk in adolescence. J Adolesc Health. 2015;56(6):599–605.View ArticlePubMedPubMed CentralGoogle Scholar
  35. Cole TJ, Faith M, Pietrobelli A, Heo M. What is the best measure of adiposity change in growing children: BMI, BMI%, BMI z-score or BMI centile? Eur J Clin Nutr. 2005;59(3):419–25.View ArticlePubMedGoogle Scholar
  36. StataCorp. Stata Statistical Software. College Station: StataCorp LP; 2013. Release 13.Google Scholar
  37. Dahly DL. Growth Mixture Modelling for Life Course Epidemiology. In: Modern Methods for Epidemiology. Berlin: Springer; 2012. p. 223–41.View ArticleGoogle Scholar
  38. Ram N, Grimm KJ. Growth Mixture Modeling: A Method for Identifying Differences in Longitudinal Change Among Unobserved Groups. Int J Behav Dev. 2009;33(6):565–76.View ArticlePubMedPubMed CentralGoogle Scholar
  39. Wang M, Bodner TE. Growth mixture modeling identifying and predicting unobserved subpopulations with longitudinal data. Organ Res Methods. 2007;10(4):635–56.View ArticleGoogle Scholar
  40. Tofighi D, Enders CK. Identifying the correct number of classes in growth mixture models. In: Advances in Latent Variable Mixture Models. Charlotte: Information Age Publishing, Inc; 2008. p. 317–41.Google Scholar
  41. Magee CA, Caputi P, Iverson DC. Identification of distinct body mass index trajectories in Australian children. Pediatric Obesity. 2013;8(3):189–98.View ArticlePubMedGoogle Scholar
  42. Nylund KL, Asparouhov T, Muthén BO. Deciding on the number of classes in latent class analysis and growth mixture modeling: A Monte Carlo simulation study. Struct Equ Model. 2007;14(4):535–69.View ArticleGoogle Scholar
  43. McLachlan G, Peel D. Finite mixture models. Hoboken: Wiley; 2004.Google Scholar
  44. Lo Y, Mendell NR, Rubin DB. Testing the number of components in a normal mixture. Biometrika. 2001;88(3):767–78.View ArticleGoogle Scholar
  45. Garden FL, Marks GB, Simpson JM, Webb KL. Body mass index (BMI) trajectories from birth to 11.5 years: relation to early life food intake. Nutrients. 2012;4(10):1382–98.View ArticlePubMedPubMed CentralGoogle Scholar
  46. Ventura AK, Loken E, Birch LL. Developmental trajectories of girls' BMI across childhood and adolescence. Obesity. 2009;17(11):2067–74.View ArticlePubMedGoogle Scholar
  47. Dziak JJ, Coffman DL, Lanza ST, Li R. Sensitivity and specificity of information criteria. State College: The Methodology Center and Department of Statistics, Penn State, The Pennsylvania State University; 2012.Google Scholar
  48. Jun H-J, Corliss HL, Nichols LP, Pazaris MJ, Spiegelman D, Austin SB. Adult body mass index trajectories and sexual orientation: The Nurses' Health Study II. Am J Prev Med. 2012;42(4):348–54.View ArticlePubMedPubMed CentralGoogle Scholar
  49. Cole T, Lobstein T. Extended international (IOTF) body mass index cut-offs for thinness, overweight and obesity. Pediatric Obesity. 2012;7(4):284–94.View ArticlePubMedGoogle Scholar
  50. Hejazi S, Dahinten VS, Marshall SK, Ratner PA. Developmental pathways leading to obesity in childhood. Health Rep. 2009;20(3):63–9.PubMedGoogle Scholar
  51. Haga C, Kondo N, Suzuki K, Sato M, Ando D, Yokomichi H, Tanaka T, Yamagata Z. Developmental trajectories of body mass index among Japanese children and impact of maternal factors during pregnancy. PLoS One. 2012;7(12):e51896.View ArticlePubMedPubMed CentralGoogle Scholar
  52. Brault M-C, Aimé A, Bégin C, Valois P, Craig W. Heterogeneity of sex-stratified BMI trajectories in children from 8 to 14years old. Physiol Behav. 2015;142:111–20.View ArticlePubMedGoogle Scholar
  53. Jun H-J, Corliss HL, Boynton-Jarrett R, Spiegelman D, Austin SB, Wright RJ. Growing up in a domestic violence environment: relationship with developmental trajectories of body mass index during adolescence into young adulthood. J Epidemiol Comm Health. 2011:jech. 2010.110932.Google Scholar
  54. Li C, Goran MI, Kaur H, Nollen N, Ahluwalia JS. Developmental trajectories of overweight during childhood: role of early life factors. Obesity. 2007;15(3):760–71.View ArticlePubMedGoogle Scholar
  55. Mustillo S, Worthman C, Erkanli A, Keeler G, Angold A, Costello EJ. Obesity and psychiatric disorder: developmental trajectories. Pediatrics. 2003;111(4):851–9.View ArticlePubMedGoogle Scholar
  56. Nonnemaker JM, Morgan‐Lopez AA, Pais JM, Finkelstein EA. Youth BMI trajectories: evidence from the NLSY97. Obesity. 2009;17(6):1274–80.PubMedGoogle Scholar
  57. Pryor LE, Tremblay RE, Boivin M, Touchette E, Dubois L, Genolini C, Liu X, Falissard B, Côté SM. Developmental trajectories of body mass index in early childhood and their risk factors: an 8-year longitudinal study. Arch Pediatr Adolesc Med. 2011;165(10):906–12.View ArticlePubMedGoogle Scholar
  58. Makgae P, Monyeki K, Brits S, Kemper H, Mashita J. Somatotype and blood pressure of rural South African children aged 6–13 years: Ellisras longitudinal growth and health study. Ann Hum Biol. 2007;34(2):240–51.View ArticlePubMedGoogle Scholar
  59. Moselakgomo VK, Toriola AL, Shaw BS, Goon DT, Akinyemi O. Body mass index, overweight, and blood pressure among adolescent schoolchildren in Limpopo province, South Africa. Revista Paulista de Pediatria. 2012;30(4):562–9.View ArticleGoogle Scholar
  60. Ellenga Mbolla BF, Okoko AR, Mabiala Babela JR, Ekouya Bowassa G, Gombet TR, Kimbally-Kaky S-G, Longo-Mbenza B. Prehypertension and Hypertension among Schoolchildren in Brazzaville, Congo. Int J Hypertens. 2014;2014:6.View ArticleGoogle Scholar
  61. Maranon R, Reckelhoff JF. Sex and gender differences in control of blood pressure. Clin Sci. 2013;125(7):311–8.View ArticlePubMedPubMed CentralGoogle Scholar
  62. Kautzky‐Willer A, Handisurya A. Metabolic diseases and associated complications: sex and gender matter! Eur J Clin Invest. 2009;39(8):631–48.View ArticlePubMedGoogle Scholar
  63. Huang R-C, Burrows S, Mori TA, Oddy WH, Beilin LJ. Lifecourse Adiposity and Blood Pressure Between Birth and 17 Years Old. Am J Hypertens. 2015;28:1056–63.View ArticlePubMedGoogle Scholar
  64. Ehret G, Munroe P, Rice K. CARDIoGRAM consortium; CKDGen Consortium; KidneyGen Consortium; EchoGen consortium; CHARGE-HF consortium. Genetic variants in novel pathways influence blood pressure and cardiovascular disease risk. Nature. 2011;478(7367):103–9.View ArticlePubMedGoogle Scholar
  65. Silva CF, Zandoná MR, Vitolo MR, Campagnolo PDB, Rotta LN, Almeida S, Mattevi VS. Association between a frequent variant of the FTO gene and anthropometric phenotypes in Brazilian children. BMC Med Genet. 2013;14(1):34.View ArticlePubMedPubMed CentralGoogle Scholar

Copyright

Advertisement