Like most places in the world, the total fertility rate in Brazil is low, and continuing to decline. The average number of births per woman was 2·3 in the year 2000 and 1·9 in 2010( 1 , 2 ). Recent studies suggest that the consequent changes in family structure may influence obesity risk in adolescence( Reference Haugaard, Ajslev and Zimmermann 3 , Reference Wang, Sekine and Chen 4 ). Some have observed that only children have higher BMI, greater fat mass, and are more likely to be overweight or obese( Reference Chen and Escarce 5 – Reference Guedes, Rocha and Silva 7 ). The influence of birth order on the occurrence of obesity has also been investigated, but the findings are inconsistent in previous studies. While some authors showed higher BMI and greater fat mass between firstborn adolescents( Reference Celi, Bini and De Giorgi 8 – Reference Stettler, Tershakovec and Zemel 11 ), others have not observed such associations( Reference Gopinath, Baur and Burlutsky 6 , Reference Li, Goran and Kaur 12 , Reference Wells, Hallal and Reichert 13 ). Studies have also shown that a larger number of siblings is associated with lower prevalence of obesity in adolescence( Reference Haugaard, Ajslev and Zimmermann 3 , Reference Wang, Sekine and Chen 4 ), regardless of whether siblings are younger or older( Reference Wang, Sekine and Chen 4 ).
There are both social and biological explanations of how birth order and/or number of siblings might influence adolescent body composition. Lower-birth-order infants tend to be smaller at birth than later-born infants( Reference Karn and Penrose 14 – Reference Wilcox, Chang and Johnson 16 ) and more likely to experience catch-up growth, a pattern of growth associated with obesity risk( Reference Wells, Hallal and Reichert 13 ). However, lower-birth-order children also tend to have fewer siblings, which might reduce their opportunities for playing games and other physical activities( Reference Hallal, Wells and Reichert 17 ). It is thus critically important to better understand the relative importance of both birth order and number of siblings, a challenge few previous studies have addressed.
To help fill this gap, we used data from the 1993 Birth Cohort of Pelotas, Brazil, which included detailed information on family structure, to evaluate the association between both birth order and number of siblings with body composition in 18-year-old adolescents.
Methods
Study design and sample
Data are from a population-based birth cohort study located in Pelotas, Brazil. Between 1 January and 30 December, 1993, all maternity hospitals in the city were visited daily and 5265 births from woman living in the city were recorded. Of this, 5249 agreed to take part in the longitudinal study. Mothers and their infants have since been followed up on numerous occasions. Visits to the full cohort took place at 11, 15 and 18 years of age. Topic-specific sub-studies were conducted at the ages of 4, 6, 9, 11, 13 and 18 years. More detailed information about the study can be found in specific methodological publications( Reference Goncalves, Assuncao and Wehrmeister 18 , Reference Victora, Araujo and Menezes 19 ).
To be included in the sample used in this analysis, individuals must have participated in the follow-ups at ages 11, 15 and 18 years, and have information on the perinatal exposures (birth order), at age 15 years (number of younger siblings and total siblings) and on the outcomes at age 18 years (BMI, fat mass and fat-free mass). Individuals who at 18 years of age were pregnant or suspected of being pregnant, in a cast, or using a wheelchair were excluded from the sample.
Measurements
Body composition at 18 years of age was assessed by fat mass index (fat mass in kg divided by height in m2) and fat-free mass index (fat-free mass in kg divided by height in m2), both collected by air displacement plethysmography (BOD POD®), using the Siri equation. These indices include the height of the individual in the calculations, and thus improving interpretation in adolescents with different heights( Reference VanItallie, Yang and Heymsfield 20 ). Moreover, the BMI (weight in kg divided by height in m2) was also assessed at the age of 18 years.
For the examination, the young individuals were attired in a top and shorts made of average compression spandex, and a silicone cap with good grip on the head.
Height was measured using an aluminium stadiometer with size of 2 m and precision of 1 mm; and weight was measured with an electronic scale connected to the BOD POD®. All measurements were taken by trained interviewers.
Data collection on major exposures was performed by asking the mother two questions. For data on birth order, we asked ‘How many times have you ever been pregnant including this pregnancy?’, and to collect information on number of younger siblings, we asked ‘Have you had any children after (name)? How many?’ Then, the two responses were combined to build the total siblings variable.
The first question was made during the perinatal follow-up, and the second at the 15-year-old follow-up.
Analytical methods
The main exposures used in the analyses were birth order (dummy-variable coded as first born, 0; second born, 1; third born, 2; and fourth born or more, 3), number of younger siblings and total siblings (each dummy-variable coded for 0, 1, 2, and more than 3 siblings). The outcomes BMI, fat mass index and fat-free mass index were assessed as continuous variables in kg/m2 and standardised as z-scores.
For controlling possible confounding factors, the following variables were included in the analysis: family income (in minimum wages), maternal education (in completed years), presence of the father (yes/no), maternal skin colour (white/black/other), maternal age at birth (in years), maternal height (in cm), gestational weight gain (in g), pregnancy smoking consumption (yes/no), pregnancy alcohol consumption (yes/no), gestational systemic hypertension (yes/no) and gestational diabetes (yes/no).
At the age of 18 years, 4563 adolescents were located, of whom 4106 were interviewed and 3974 adolescents were included in our analysis with complete data on body composition. Those who completed the interviews, added to those known to have died, represented 81·3 % of the original cohort. Of those located, 127 (2·3 %) refused to participate in the study and 330 (7·2 %) were considered losses, 196 were found living in other cities and were not interviewed. The maternal height variable had the maximum percentage of unknown observations (7·9 %).
Descriptive analysis of the exposures and outcomes were performed, presenting the absolute and relative frequencies of categorical variables, and measures of central tendency and dispersion for continuous variables.
Later, linear regression models were used to estimate associations between exposures (birth order, number of younger siblings and total siblings) and outcomes (BMI, fat mass index and fat-free mass index). We report crude estimates, as well as those adjusted for the potential confounders described earlier. All analyses were sex stratified.
To evaluate the relative importance of birth order and number of siblings, we compared the fit of three different models based on their respective R 2, adjusted R 2, Akaike's information criterion and mean squared error. In model 1, only birth order and number of younger siblings were included as independent variables. In model 2, we included these variables plus potential confounders (Fig. 1). Model 3 included the total siblings (which is equal to the birth order plus the number of younger siblings, less one), and the confounding variables of model 2. The purpose of this third model was to evaluate if the position of the adolescent in the family influences the body composition, or whether this association is rather due to the number of siblings. The model that presented the best fit was the one with higher values of R 2 and adjusted R 2, and lower values of Akaike's information criterion and mean squared error.
Ethical considerations
The present study is part of the 18-year-old follow-up in the 1993 birth cohort, titled ‘Early and contemporary influences on body composition, human capital, mental health and complex chronic diseases precursors in the 1993 birth cohort of Pelotas, Brazil’, which was approved by the Ethics Committee of the Medicine School of the Federal University of Pelotas in the official letter numbered 05/11.
All participants signed an informed consent form before the procedures, interviews and examinations in both follow-ups. At the age of 18 years, body composition examinations did not present health risks to participants; nevertheless, girls were always asked about the possibility of pregnancy. Pregnant women and those suspected of being pregnant did not perform the tests.
Results
Characteristics of sample studied are shown in Table 1. Most of the adolescents were the firstborns (39·8 %), and at the age of 15, 30·8 % had one sibling. Regarding the presence of the father, 11·9 % of adolescents had no father living in the house when they were born. In relation to the mothers, one-third reported smoking during pregnancy (32·9 %) and a small proportion reported drinking alcohol during pregnancy (5·2 %). The average maternal education at birth was about 7 (sd 3·5) years and maternal age at birth was 26 (sd 6·4) years on average. The median family income was 2·6 times (interquartile range 1·5–4·6) the minimum wage.
MMW, monthly minimum wages.
* Maximum percentage of unknown observations: (n 326; 7·9 %) for the maternal height variable.
Table 2 shows the mean and standard deviations of BMI, fat mass index and fat-free mass index among male and female adolescents, and in the total population. These averages were 23·4 (sd 4·5) kg/m2, 6·1 (sd 3·9) kg/m2 and 17·3 (sd 2·5) kg/m2, respectively. Although BMI was the same for both sexes, girls had an average fat mass index almost two times higher than boys, while the mean fat-free mass index was higher in males. The prevalence of overweight and obesity in the adolescents, using the BMI-for-age WHO criterion, was 17·2 (95 % CI 16·0, 18·3) % and 10·1 (95 % CI 9·2, 11·0) %, respectively.
Table 3 presents the fit of the three models analysed, for each of the three outcomes investigated. In general, model 3 showed the best fit for each outcome, and was more parsimonious than model 2. While models 1 and 2 indicated that both birth order and number of younger siblings were associated with the outcome (data not shown), the better fit of model 3 led us to conclude that the total siblings was the key influence on fat mass index, fat-free mass index and BMI, rather than the position that the adolescents occupy among siblings.
AIC, Akaike's information criterion; MSE, mean squared error.
* Birth order and total siblings adjusted one for each other.
† Birth order and younger siblings adjusted for confounding variables: family income, maternal education, presence of the father, maternal skin colour, maternal age at birth, maternal height, gestational weight gain, pregnancy smoking consumption, pregnancy alcohol consumption, gestational systemic hypertension and gestational diabetes.
‡ Total siblings adjusted for confounding variables: family income, maternal education, presence of the father, maternal skin colour, maternal age at birth, maternal height, gestational weight gain, pregnancy smoking consumption, pregnancy alcohol consumption, gestational systemic hypertension and gestational diabetes.
The distributions of BMI, fat mass index and fat-free mass index, across total siblings, by sex are presented in Fig. 2. Boys who had three or more siblings had lower BMI mean (22·6, 95 % CI 22·3, 22·9) compared to those of only children (24·4, 95 % CI 23·8, 25·1). Fat mass index was associated to total siblings in both sexes. Adolescents with three or more siblings had lower fat mass index mean than those who had no siblings. Girls who had three or more siblings had higher fat-free mass index mean (15·7, 95 % CI 15·5, 15·8) than those who had none (15·3, 95 % CI 15·1, 15·5).
Table 4 shows the crude and adjusted analyses of the association between total siblings and the outcomes studied according to sex. After adjustment for possible confounders (family income, maternal education, presence of the father, maternal skin colour, maternal age at birth, maternal height, gestational weight gain, pregnancy smoking consumption, pregnancy alcohol consumption, gestational systemic hypertension and gestational diabetes), the total siblings remained associated with fat mass index in both sexes with an inverse linear trend in males (P< 0·001). Boys who had three or more siblings showed a decrease of 0·37 z-scores (95 % CI − 0·52, − 0·23) in fat mass index compared to those who had no siblings. Fat-free mass index was associated with the total siblings only in girls (P= 0·035). Those with three or more siblings showed a 0·06 increase in z-score (95 % CI − 0·04, 0·17) compared to girls who had none. Moreover, the total siblings remained inversely related to BMI in males (P< 0·001). Those boys with three or more siblings had a reduction of 0·39 z-scores in BMI when compared to an only child (95 % CI − 0·55, − 0·22).
* For family income, maternal education, presence of the father, maternal skin colour, maternal age at birth, maternal height, gestational weight gain, pregnancy smoking consumption, pregnancy alcohol consumption, gestational systemic hypertension and gestational diabetes.
† Linear trend test.
Discussion
An important result of the present study was the association between number of total siblings and fat mass index in adolescents of both sexes. The smaller the number of total siblings, the greater the fat mass index, even after adjustment for possible confounding factors. Furthermore, it was shown that the number of total siblings remained directly associated with the fat-free mass index in girls, even after adjustment for potential confounders. However, it is important to consider both associations in the girls were weak, although significant.
Similar results were found in a study in India with women only, which showed a negative correlation between number of siblings and both fat mass and fat-free mass. However, the Indian study was not adjusted for potential confounding factors and body composition was assessed by skinfold thickness( Reference Ghosh and Bandyopadhyay 21 ). Studies that assess the body composition of adolescents using reference measurement methods are still rare in the literature, which hinders the comparison of results.
Concerning BMI, it was shown that this variable remained inversely associated with the number of total siblings in males in the adjusted analysis. One possible explanation for this association is that the greater number of siblings leads to an increase in family size, which has an inverse effect on family income and, consequently, on obesity. According to a study conducted in Brazil, there is a direct association between income and obesity in adolescents of both sexes( 22 ).
Corroborating the results of the present study, a longitudinal study conducted in the USA showed that only children had higher mean BMI compared to those who had at least two brothers, even after adjustment for possible confounding factors( Reference Chen and Escarce 5 ). Other studies analysing overweight and/or obesity in adolescents, using different methods( Reference Cole, Bellizzi and Flegal 23 – 25 ), observed inverse associations between number of siblings and overweight and/or obesity, with only children having a higher prevalence of overweight and obesity when compared with individuals who had siblings( Reference Haugaard, Ajslev and Zimmermann 3 , Reference Wang, Sekine and Chen 4 , Reference Gopinath, Baur and Burlutsky 6 , Reference Guedes, Rocha and Silva 7 , Reference Mushtaq, Gull and Shahid 26 ). In contrast, Hesketh and colleagues( Reference Hesketh, Qu and Tomkins 27 ) showed no association between number of siblings and overweight in Chinese adolescents. Results of stratified analysis by sex are presented for one study only( Reference Wang, Sekine and Chen 4 ), which hinders the comparison of findings.
Some limitations in the present study should be highlighted. First, as with any longitudinal study, potential selection biases due to loss to follow-up are a substantial limitation. To help address this concern, we compared the sample participants with the original participants examined in 1993. Adolescents with worse socioeconomic and nutritional profiles were slightly less likely to be followed up. Socioeconomically intermediate participants were more likely to located as compared with very poor or very rich adolescents (81·8 v. 75·6 and 76·1 %, respectively). Furthermore, participants whose mothers had no schooling were less likely to be followed up compared to those whose mothers had 9 or more years of schooling (69·4 v. 77·5 %). However, the magnitude of such differences is modest, therefore minimising the likelihood of bias( Reference Goncalves, Assuncao and Wehrmeister 18 ).
Another limitation is the lack of information on the number of siblings of adolescents at 18 years of age, which limits us to the information collected at the 15-year-old follow-up. Finally, the measurements of body composition are only available at age 18 years; therefore, we were not able to investigate the development of body composition across the entire life course.
Strengths of the present study include its prospective design, the use of air displacement plethysmography (BOD POD®) as a method for assessing body composition, and the high follow-up rate at age 18 years (81·3 %) ensuring the representativeness of the sample despite some small differences between study participants and those lost to follow-up.
In conclusion, only children of both sexes are more likely to have higher fat mass index compared to those who have siblings. Girls who have no siblings present the lowest levels of fat-free mass, while only boy children have greater BMI.
Given the reduction in the fertility rate in Brazil and in the world( 1 , 2 , 28 , 29 ) and the larger number of couples choosing to have an only child( 29 ), our research further highlights the need for prevention of obesity or fat mass accumulation in adolescents who are only children.
Acknowledgements
All colleagues and workers who contributed to the study included since its 1st year.
The 1993 birth cohort study is currently supported by the Wellcome Trust through the programme titled ‘Major Awards for Latin America on Health Consequences of Population Change’.
Contributions of authors are as follows: F. O. M. was responsible for the literature review, results analysis, data interpretation and writing of the paper; M. C. F. A., F. C. B. and D. L. D. helped with the overall concept of the paper, the design of the analysis, and writing of the paper; A. A. S., C. L. M. and A. J. D. B. contributed with the methodology issues and the data analysis. All authors discussed the results and implications and commented on the manuscript at all stages.
The authors declare that they have no conflicts of interest.