Number of children and body composition in later life among men and women: Results from a British birth cohort study.

BACKGROUND: Although research has found associations between increasing number of children and higher body mass index (BMI), there has been limited research investigating the association with body composition despite abdominal fat being independently associated with cardiovascular and metabolic risk. Most existing research has focussed on women, but investigating the relationship in men can help distinguish biological effects of pregnancy from social pathways related to parenthood.
METHODS: Using the MRC National Survey of Health and Development (NSHD) multiple regression models were applied to test associations between number of children and body composition at age 60-64 years (N = 2229) and body mass index (BMI) and waist circumference (WC) at ages 60-64 and 69 years (N = 2149).
RESULTS: In adjusted models, associations were observed between increasing numbers of children and increasing fat-adjusted lean mass index in women (p = 0.06). Among men, those with children had 0.59kg (95% CI: 0.15 to 1.02) greater lean mass index than those without and fat:lean mass ratio was greater in those with 4+ children because of their slightly higher mean fat mass. Weak evidence of a higher android:gynoid mass ratio among mothers (0.03, 95% CI: 0.00,0.06, p = 0.1) was observed with no associations with fat mass index,android or gynoid fat mass. Increasing parity in women was associated with increasing BMI at 60-64 years and, more strongly, with increasing BMI at 69 years. Associations among men were also observed more clearly with BMI at 69 years.
CONCLUSION: There was little evidence of a consistent association between number of children and body composition in early old age. The strongest associations are observed for lean, rather than fat mass, and in men rather than women, suggesting little evidence of biological effects of pregnancy. The results indicate social pathways associated with parenthood are the likely underlying mechanisms, with suggestion of selection into parenthood among men.

Introduction

Characteristics of family formation such as number of children, age at first birth, age at last birth and birth intervals have been associated with subsequent overweight and obesity [1-7]. The majority of studies have focussed only on women due, at least in part, to the assumption that there is a biological pathway between pregnancy, child birth and obesity. The majority of studies have found parous women to have greater body mass index (BMI), greater risk of obesity [1, 8–10] or faster weight gain trajectories [2, 3, 5, 11] than nulliparous women, with associations also observed between increasing number of children and higher BMI. However, differences in associations do exist according to the age at BMI measurement [5] with few studies investigating the relationship with BMI measured past mid-life, area of residence (non-metropolitan compared to metropolitan) [12], and according to race and ethnicity [3, 5, 10–12]. Social and behavioural factors related to parenting provide an alternative explanation to the biological effects of pregnancy for any such association, as parenthood has been found to be related to poorer health behaviours. Differences have been observed in the diets of parents and those who were not parents [13, 14], for example, parents have been found to eat more saturated fat [14, 15]. Additionally, previous studies have shown that those with 4 or more children are more likely to be smokers, physically inactive and from a manual social class [1]. It has therefore been suggested that parenthood changes the social contexts in which people live in a way that impacts weight gain [11].

Biological pathways, such as reduced exposure to oestrogen or increased levels of relative insulin resistance [16], have been implicated in the association between pregnancy and cardiovascular disease (CVD), with the effect of pregnancy on abdominal obesity being one potential explanatory mechanism [17]. Hence, biological effects of pregnancy may be particularly strong when considering fat distribution and body composition, rather than general markers of adiposity such as BMI. There is increasing evidence that fat mass and android and visceral fat are related to cardiovascular and metabolic risk factors independently of BMI and even among those who are not overweight [18-20]. Few studies have assessed the association between parity and body composition. Existing studies, which have shown an association between parity and increased visceral adiposity measured using computed tomography (CT), have been based on small sample sizes (n≤170) [21, 22]. One of these studies considered the change in adiposity pre- and post-pregnancy [22]. The other study, based on a highly age heterogeneous sample (age 18 to 76), also considered percent body fat measured by dual-energy x-ray absorptiometry (DXA) scans finding no association with parity after adjustment for age, physical activity and smoking [21]. A larger study using DXA scans to measure body composition at a mean age of 59.7 years carried out in an Asian population, did not consider number of children, but found that earlier age of first birth was associated with increased trunk fat mass [23]. However, these findings may not be generalizable to western populations as visceral and abdominal fat distribution for a given BMI is higher in Asian compared to western populations [24].

There has been growing interest in the impact of becoming a father on health, with a few studies suggesting that number of children is associated with obesity and higher BMI in men as well as women [1, 9]. It remains unclear as to whether associations are stronger in men or women as results have so far been inconsistent [2], and to our knowledge there have been few studies investigating associations with body composition in men, with one study in Brazil testing the association with percentage fat mass in both men and women and finding no association for either sexes [25]. Investigating the relationship between parenthood and obesity in men as well as women can also help distinguish a biological effect of pregnancy from a social pathway related to parenthood [8]; similar associations in both sexes would suggest that the underlying mechanisms are social as opposed to biological.

Using data from the Medical Research Council (MRC) National Survey of Health and Development (NSHD), a British birth cohort born in 1946, we investigated the relationship between number of children and body composition measured using DXA scans at ages 60–64 years in men and women. Previous research in NSHD has shown that greater number of children in women, and being a parent in men, was associated with higher BMI between age 36 and 53 years [1]. Therefore we also assess whether such associations remain into old age, investigating the association between number of children and BMI and waist circumference (WC) measured at the same age as body composition (60–64 years) and also at 69 years (where DXA measures were not collected). We assess the extent to which associations can be explained by lifestyle and behavioural factors related to family life, by controlling for childhood cognition, childhood and adult socioeconomic position, education and smoking status.

Materials and methods

The MRC NSHD is a socially stratified sample of singleton births consisting of 2547 women and 2815 men born in one week in March 1946 in England, Scotland and Wales. There have been 24 follow ups including interviews with the mothers, assessments by school doctors and teachers during childhood, and home visits and postal questionnaires during adult life. When study members were aged 60–64 years, 2856 individuals were invited for assessment at one of six clinical research facilities (CRF) or to have a nurse visit them in their own home. Invitations were not sent to those who had died (n = 778), who were living abroad (n = 570), who had previously withdrawn from the study (n = 594), or who had been lost to follow-up (n = 564). Of those invited, a total of 2229 (78%) responded and were assessed, with 1690 (59%) attending a CRF and the remaining 539 being visited at home. At the clinic visit study members underwent (DXA) scans from which measures on body composition were derived and anthropometric measurements. The most recent data collection took place in 2015 when participants were aged 69 years when 2149 were interviewed and examined in their own homes by a team of research nurses where anthropometric measures were taken.

Ethical approval for the age 60–64 data collection was obtained from the Greater Manchester Local Research Ethics Committee and the Scotland A Research Ethics Committee, and for the age 69 data collection from the Queens Square Research Ethics Committee and the Scotland A Research Ethics Committee. Written informed consent was obtained from the study members at each data collection.

Outcome variables

Measures of body composition were obtained at the CRFs using a QDR 4500 Discovery DXA scanner (Hologic Inc, Bedford, MA, USA) whilst the individuals were in a supine position. The measures of body composition used in these analyses were whole body, android and gynoid fat mass and whole body lean mass defined as total mass minus fat and bone mass. For all measures, mass from the head was excluded and measures were converted to kg. All scans were reviewed and centrally analysed by a single operator using APEX 3.1 software (Hologic Inc., Bedford, MA, USA), with use of a the European Spine Phantom at the start and end of the study to ensure cross-calibration between scanners [26, 27]. Fat mass index (kg/m2) and lean mass index (kg/m2) were calculated by dividing the measures by height2. Fat-to-lean mass ratio and android-to-gynoid fat mass ratio were derived.

The anthropometric measures of height, weight, and waist circumference were measured by the nurses using standard protocols at ages 60–64 years at both CRFs and home visits and at 69 years at the home visits. BMI (defined as weight (kg)/height(m)2) was calculated.

Explanatory variable

Number of children was derived from the number of new live births since age 53 years reported at age 60–64 years added to previous reports of number of live births collected throughout adult life at ages of 36, 43 and 53 years. For the purpose of analysis and due to the small numbers with four or more children, all those who reported having four or more children were combined into one group.

Confounding variables

Potential confounders, identified from previous research as being associated with both number of births and adiposity, were childhood cognitive function, childhood and adult socioeconomic position, education and cigarette smoking. Cognitive function in childhood was calculated from the sum of standardised verbal and non-verbal cognitive tests taken at age 8 years. Socioeconomic position was measured by father’s occupational social class when the study member was aged 4 (or 7 or 11 if missing) and by occupational social class of the head of the household at age 53 years. The Register General classification system was used to categorise into I Professional; II Intermediate; IIINM Skilled (Non-Manual); IIIM Skilled (Manual); IV Partly skilled; V Unskilled. Highest educational qualification by age 26 was classified into (1) No Qualifications; (2) Below O-Level; (3) O-Level or Equivalent; (4) A-Level or Equivalent; (5) Degree Level or higher. Smoking status was self-reported at the age of 60–64 years (and at previous ages) and was classified as current smoker, ex-smoker and never smoker.

Statistical analysis

Multivariable regression models were used to examine the associations between number of children (with no children as the reference group) and each DXA measure and BMI and WC at 60–64 years. A test for deviation from a linear trend across groups of number of children was carried out and where significant deviations from linearity were detected (P≤0.05), tests for heterogeneity across groups were carried out. Where there was no evidence of deviation form linearity, number of children was fitted as a continuous variable to obtain a test for linearity. Analyses were carried out separately for men and women and, in addition, sex by number of children interactions were included to test whether the associations differed between men and women in models including both sexes. Models were then adjusted for all potential confounding variables and for lean mass index as the outcome, models were further adjusted for fat mass index due to the adaptive changes in lean mass with increasing fat mass [27]. Similar models were fitted for BMI and WC at age 69 years to assess whether associations were consistent with those at age 60–64 years. Since it has previously been observed that there are differences in outcomes between those with children and those without [1, 28], models were also fitted to test differences between those with at least one child and those with no children and then, in parents only, to test associations with number of children using one child as the reference group. All analysis was carried out using Stata14 (Stata Corp., College Station, TX, USA).

Multiple imputation

A total of 2,331 participants with complete information on number of children and at least one outcome were included in analyses. Since there was some missing information for every confounder, multiple imputation was used in order to maintain the sample size under the assumption that data were missing at random. There were 125 (5.36%) missing and 2,206 (94.63%) complete for education level, 119 (5.10%) missing and 2,212 (94.89%) complete for childhood social class, 294 (12.61%) missing and 2,037 complete (87.39%) for smoking status, 254 (10.89%) missing and 2,077 (89.10%) complete for cognitive score and 43 (1.84%) missing and 2,888 (98.16%) complete for adult social class. There were no significant differences between participants with and without complete information for important confounders, except for waist circumference at 69 years, which was slightly higher in those with missing information, and number of children where those with missing data were more likely to have 4 or more children. A total of 20 imputed data sets were obtained using chained equations and Rubin’s rule was used to combine the estimates from regression models from each of the 20 datasets. Imputation models included the outcome variables and all variables in the adjusted models, as well as additional variables to aid the imputation process (other measures of socioeconomic position, earlier life BMI and waist circumference). For fat mass index, lean mass index and fat-to-lean ratio, 304 (20%) (Male = 152, Female = 152) of the analytic sample had missing data on at least one covariate. For android fat, gynoid fat and android-gynoid ratio, the number was 318 and for BMI and waist circumference aged 60–64, 456 and 457, respectively. At age 69, the number was 348 for both BMI and waist circumference.

Sensitivity analysis

The characteristics of those who attended a CRF, and who therefore had a DXA scan at age 60–64 years, have previously found to differ from those who had a home visit, with indication that they tended to be healthier [29]. Therefore a sensitivity analysis was carried out to examine whether this difference in sample characteristics might have influenced associations in the subgroup who had DXA measures. We compared the difference in mean BMI and WC at age 60–64 years between those who had DXA measures and those who had a home visit (and therefore did not have DXA measures). We then repeated analysis of BMI and WC at age 60–64 years in the sample who had DXA measures.

Results

Characteristics of the sample (N = 2331, 52.2% female) who had at least one of the outcome variables and a valid record of number of live births are presented in Table 1. Mean fat-to-lean ratio, fat mass index and gynoid fat were higher in women compared to men, whilst all other measures were higher in men. Mean BMI was similar in men and women at both ages 60–64 years and 69 years and, as expected, mean WC was higher in men compared to women. The median number of children was 2 for both men and women, with few individuals with 4 or more children (7.27% and 8.46%, respectively).

Table 1

Characteristics of men and women in the analytic sample.

Variable	Male		Female
	Number (%)	Mean (SD)	Number (%)	Mean (SD)
Fat Mass Index at 60–64 years (kg/m2)	701	7.72 (2.33)	787	11.1 (3.47)
Lean Mass Index at 60–64 years (kg/m2)	701	17.42 (1.97)	787	14.18 (1.84)
Fat- Lean Mass ratio at 60–64 years	701	0.44 (0.11)	787	0.77 (0.19)
Android Fat at 60–64 years (kg)	731	2.49 (0.97)	830	2.33 (1.02)
Gynoid Fat at 60–64 years (kg)	731	3.76 (1.01)	830	5.13 (1.48)
Android-Gynoid Ratio at 60–64 years	731	1.18 (0.18)	830	0.92 (0.15)
Body Mass Index at 60–64 years (kg/m2)**	990	27.85 (4.08)	1,117	28.00 (5.51)
Waist Circumference at 60–64 years (cm) **	990	100.68 (11.02)	1,115	92.56 (13.09)
Number of Children*
0	163 (14.63)		133 (10.93)
1	117 (10.50)		142 (11.67)
2	551 (49.46)		563 (46.26)
3	202 (18.13)		276 (22.68)
4+	81 (7.27)		103 (8.46)
Total	1,114		1,217
Confounding Variables
Childhood Cognitive Function (z-score)	985	0.082 (0.83)	1,092	0.12 (0.80)
Childhood Socioeconomic Position
I Professional	80 (7.53)		80 (6.96)
II Intermediate	228 (21.45)		250 (21.76)
IIINM Skilled (Non-Manual)	174 (16.37)		190 (16.54)
IIIM Skilled (Manual);	328 (30.86)		353 (30.72)
IV Partly skilled;	189 (17.78)		216 (18.80)
V Unskilled	64 (6.02)		60 (5.22)
Adult Socioeconomic Position
I Professional	156 (14.30)		108 (9.02)
II Intermediate	453 (41.52)		501 (41.85)
IIINM Skilled (Non-Manual)	114 (10.45)		179 (14.95)
IIIM Skilled (Manual);	273 (25.02)		250 (20.89)
IV Partly skilled;	73 (6.69)		109 (9.11)
V Unskilled	22 (2.02)		50 (4.18)
Education Level
No Qualifications	348 (33.05)		369 (32.00)
Below O-Level	58 (5.51)		110 (9.54)
O-Level or Equivalent	155 (14.72)		305 (26.45)
A-Level or Equivalent	318 (30.20)		303 (26.28)
Degree Level or higher	174 (16.52)		66 (5.27)
Smoking Status
Current Smoker	109 (11.32)		122 (11.36)
Ex-Smoker	580 (60.23)		561 (52.23)
Never Smoked	274 (28.45)		391 (36.41)

*Taken from the number of children reported at age 60–64 years.

** Values for age 69 are as follows: Male BMI 28.15 (s.d 4.55, n = 961), female BMI 28.21 (s.d 5.96, n = 1,047), male WC 100.93 (s.d 12.21, n = 961), female WC 92.39 (s.d 13.91, n = 1,053).

Body composition

Fat mass index and lean mass index were available for 701 men and 787 women and android and gynoid fat for 731 men and 830 women. There was little evidence of a trend across number of children categories for the majority of body composition measures in either men or women (Table 2). Associations between number of children and lean mass index were observed for both sexes. Higher mean lean mass index was observed in men with children compared to those without and there was weak evidence that increasing parity was associated with increasing lean mass index in women (although there was no evidence of a sex interaction). These findings were confirmed in the analysis comparing men with one or more children to those with no children where the mean difference in lean mass index was 0.58kg (95% CI: 0.15 to 1.02, P = 0.008) and among those with at least one child for women (Table 3). Among men, there was also a suggestion of a non-linear relationship between number of children and fat:lean mass ratio where those with 1, 2 or 3 children had lower ratios than men with either 0 or 4+ children (Table 2). Men with 4+ children, despite having greatest lean mass index, also had greatest mean fat mass index, resulting in a ratio which was greater than men with fewer children. The higher fat mass in men with 4+ children was confirmed in analyses including only parents (p = 0.05, Table 3). Associations were hardly affected following adjustment for confounders (Table 2) with slight attenuation of effects for lean mass index and very slight strengthening for the ratio in men. Among women, there was weak evidence that those with children had higher android:gynoid fat mass ratio than those without children (0.03 (95% CI: 0.00–0.06) higher, p = 0.1) (Table 3).

Table 2

Adjusted and unadjusted regression coefficients (mean differences) and 95% confidence intervals (95% CI) for each measure of body composition at 60–64 years and body mass index and waist circumference at 60–64 and 69 years in men and women.

Outcome Variable	N		UnadjustedRegression coefficients (95% CI) comparedwith baseline group of 0 children				Adjusted**Regression coefficients (95% CI) comparedwith baseline group of 0 children
						P value,trend					P value,trend
		1	2	3	4+		1	2	3	4+
Fat Mass Index (kg/m2)
Male	701	0.18(-0.55,0.91)	-0.02(-0.55,0.50)	-0.33(-0.95,0.29)	0.69(-0.11,1.48)	0.8	-0.07(-0.81,0.67)	-0.03(-0.57,0.50)	-0.44(-1.06,0.19)	0.63(-0.18,1.44)	0.9
Female	787	0.66(-0.33, 1.64)	0.13(-0.65,0.91)	0.42(-0.45,1.29)	0.75(-0.33,1.83)	0.3	0.38(-0.62,1.38)	-0.03(-0.82,0.76)	0.23(-0.64,1.10)	0.42(-0.68,1.52)	0.6
Lean Mass Index (kg/m2)
Male	701	0.98 (0.37,1.60)	0.56 (0.11,1.00)	0.56 (0.03,1.08)	1.02 (0.35,1.69)	0.009*	0.84 (0.21,1.46)	0.53(0.07,0.98)	0.46(-0.07,0.99)	1.06 (0.37,1.74)	0.02*
Female	787	0.14(-0.39,0.66)	0.03(-0.39,0.44)	0.34(-0.12,0.80)	0.42(-0.15,0.99)	0.07	0.10(-0.43,0.64)	0.03(-0.39,0.46)	0.33(-0.14,0.80)	0.39(-0.20,0.98)	0.09
Fat: Lean Mass ratio
Male	701	-0.01(-0.05 0.02)	-0.02(-0.04,0.01)	-0.03(-0.06,0.00)	0.01(-0.03,0.05)	0.07*	-0.02(-0.06,0.01)	-0.02(-0.04,0.01)	-0.04(-0.07,-0.01)	0.01(-0.03,0.05)	0.05*
Female	787	0.05 (0.00,0.10)	0.02(-0.03,0.06)	0.02(-0.03,0.06)	0.04(-0.02,0.09)	0.7	0.03(-0.02,0.08)	0.00(-0.04,0.05)	0.00(-0.04,0.05)	0.01(-0.05,0.7)	0.9
Android Fat (kg)
Male	731	0.09(-0.21,0.39)	0.02(-0.19,0.24)	-0.08(-0.33,0.17)	0.24(-0.08,0.56)	0.7	-0.01(-0.31,0.29)	0.00(-0.21,0.22)	-0.13(-0.39,0.12)	0.22(-0.11,0.55)	0.8
Female	830	0.16(-0.12,0.45)	0.07(-0.15,0.30)	0.15(-0.10,0.40)	0.25(-0.06,0.57)	0.2	0.13(-0.16,0.42)	0.06(-0.17,0.29)	0.13(-0.12,0.38)	0.23(-0.09,0.55)	0.2
Gynoid Fat (kg)
Male	731	0.14(-0.17,0.45)	0.02(-0.21,0.24)	-0.01(-0.27,0.26)	0.31(-0.02,0.65)	0.4	0.08(-0.24,0.40)	0.01(-0.22,0.24)	-0.03(-0.30,0.24)	0.31(-0.04,0.66)	0.4
Female	830	0.10(-0.31,0.52)	-0.09(-0.42,0.24)	0.00(-0.36,0.36)	0.20(-0.25,0.65)	0.7	0.04(-0.38,0.47)	-0.11(-0.44,0.23)	0.01(-0.38,0.36)	0.16(-0.30,0.62)	0.7
Android/Gynoid Ratio
Male	731	-0.00(-0.06,0.05)	0.01(-0.03,0.05)	0.00(-0.05,0.05)	0.01(-0.05,0.7)	0.9	-0.01(-0.07,0.05)	0.01(-0.03,0.05)	-0.01(-0.06,0.04)	0.01(-0.05,0.08)	0.9
Female	830	0.03(-0.01,0.07)	0.03 (0.00,0.07)	0.02(-0.01,0.06)	0.03(-0.01,0.08)	0.2	0.02(-0.02,0.07)	0.03 (0.00,0.06)	0.02(-0.01,0.06)	0.03(-0.01,0.08)	0.2
Body Mass Index (kg/m2) 60–64 years
Male	990	0.77(-0.26,1.80)	0.48(-0.27,1.22)	0.36(-0.52,1.25)	1.34 (0.17,2.51)	0.1	0.44(-0.60,1.47)	0.40(-0.34,1.15)	0.18(-0.70,1.06)	1.29 (0.13,2.45)	0.2
Female	1,117	0.44(-0.91,1.80)	0.35(-0.74,1.43)	1.19 (0.00,2.38)	1.43(-0.04,2.91)	0.02	0.15(-1.23,1.51)	0.10(-0.99,1.20)	0.86 (-0.33,2.06)	1.08(-0.41,2.57)	0.05
Waist Circumference (cm) 60-64 years
Male	990	0.99(-1.80,3.78)	0.39(-1.63,2.42)	0.49(-1.91,2.88)	2.87(-0.28,6.02)	0.2	-0.07(-2.84,2.70)	0.16(-1.86,2.18)	0.02(-2.35,2.39)	2.49(-0.65,5.63)	0.3
Female	1,115	1.56(-1.66,4.79)	0.61(-1.97,3.19)	1.30(-1.53,4.13)	4.34 (0.82,7.86)	0.06	0.81(-2.43,4.05)	-0.07(-2.68,2.53)	0.44(-2.41,3.30)	3.57(-0.01,7.12)	0.2
Body Mass Index (kg/m2) 69 years
Male	961	1.46 (0.26,2.65)	0.37(-0.53,1.27)	0.75(-0.29,1.80)	1.92 (0.58,3.26)	0.01*	1.06(-0.13,2.24)	0.30(-0.60,1.19)	0.57(-0.47,1.61)	1.66 (0.32,3.00)	0.07*
Female	1,047	0.23(-1.22,1.68)	0.24(-0.92,1.40)	1.06(-0.20,2.33)	2.19 (0.61,3.78)	0.003	-0.07(-1.52,1.38)	-0.09(-1.26,1.08)	0.74(0.-0.53,2.02)	1.75 (0.15,3.35)	0.01
Waist Circumference (cm) 69 years
Male	961	2.55(-0.66,5.77)	0.57(-1.85,2.99)	1.34(-1.46,4.14)	5.97 (2.37,9.57)	0.006*	1.36 (-1.84,4.57)	0.28(-2.14,2.69)	0.89(-1.90,3.69)	5.12 (1.51,8.73)	0.04*
Female	1,053	-0.18(-3.72,3.36)	0.13(-2.70,2.97)	2.00(-1.10,5.10)	3.70(-0.17,7.58)	0.02	-1.25(-4.80,2.29)	-0.96(-3.82,1.90)	0.91(-2.20,4.03)	2.34(-1.57,6.25)	0.09

*Test for heterogeneity across groups when test for deviation from linearity gives P≤0.05.

** Adjusted model includes childhood cognitive function, childhood and adult socioeconomic position, education and cigarette smoking.

Table 3

Adjusted (adjusted model includes childhood cognitive function, childhood and adult socioeconomic position, education and cigarette smoking) regression coefficients (mean differences) and 95% confidence intervals (95% CI) for each measure of body composition and anthropometric measures aged 60–64 years and 69 years for (a) parents versus non-parents, and (b) number of children for parents only–reference group 1 child.

Outcome Variable	N	Model a–Regression Coefficient (95% CI) comparing those with children to those with no children	P value	N	Model b—Regression coefficients (95% CI) comparedwith baseline group of 1 child, those with no children excluded			P value,trend
					2	3	4+
Fat Mass Index (kg/m2)
Male	701	-0.08 (-0.59,0.44)	0.8	607	0.05 (-0.57,0.66)	-0.38 (-1.06,0.32)	0.68 (-0.17,1.54)	0.05*
Female	787	0.13 (-0.62,0.88)	0.7	691	-0.38 (-1.16,0.40)	-0.12 (-0.97,0.74)	0.03 (-1.04,1.10)	0.6
Lean Mass Index (kg/m2)
Male	701	0.58 (0.15,1.02)	0.008	607	-0.31 (-0.82,0.21)	-0.38 (-0.96,0.20)	0.21 (-0.51,0.92)	0.7
Female	787	0.15 (-0.26,0.55)	0.5	691	-0.06 (-0.48,0.36)	0.24 (-0.22,0.72)	0.29 (-0.29,0.87)	0.09
Fat: Lean Mass ratio
Male	701	-0.02 (-0.04,0.00)	0.1	607	0.01 (-0.02,0.04)	-0.01 (-0.05,0.02)	0.03 (-0.01,0.07)	0.07*
Female	787	0.01 (-0.03,0.05)	0.7	691	-0.03 (-0.07,0.02)	-0.03 (-0.07,0.02)	-0.02 (-0.08,0.04)	0.5
Android Fat (kg)
Male	731	-0.01 (-0.22,0.20)	0.9	635	0.03 (-0.23,0.28)	-0.12 (-0.41,0.16)	0.22 (-0.12,0.57)	0.8
Female	830	0.10 (-0.12,0.32)	0.4	731	-0.07 (-0.29,0.16)	0.01 (-0.24,0.26)	0.09 (-0.22,0.41)	0.4
Gynoid Fat (kg)
Male	731	0.03 (-0.19,0.25)	0.8	635	-0.06 (-0.32,0.21)	-0.11 (-0.41,0.18)	0.23 (-0.14,0.59)	0.4
Female	830	-0.04 (-0.36,0.28)	0.8	731	-0.14 (-0.47,0.19)	-0.04 (-0.40,0.32)	0.11 (-0.34,0.56)	0.4
Android/Gynoid Ratio
Male	731	0.00 (-0.03,0.04)	0.8	635	0.02 (-0.03,0.07)	0.00 (-0.05,0.05)	0.02 (-0.04,0.09)	0.9
Female	830	0.03 (0.00,0.06)	0.1	731	0.00 (-0.03,0.04)	0.00 (-0.04,0.03)	0.01 (-0.04,0.05)	0.9
Body Mass Index (kg/m2) 60–64 years
Male	990	0.43 (-0.28,1.14)	0.2	841	-0.03 (-0.91,0.85)	-0.27 (-1.26,0.71)	0.79 (-0.45,2.04)	0.5
Female	1,117	0.38 (-0.67,1.42)	0.5	995	-0.03 (-1.07,1.02)	0.73 (-0.43,1.90)	0.88 (-0.56,2.33)	0.06
Waist Circumference (cm) 60–64 years
Male	990	0.28 (-1.64,2.20)	0.8	841	0.25 (-2.10,2.60)	-0.03 (-2.61,2.68)	2.36 (-0.96,5.68)	0.3
Female	1,115	0.44 (-2.05,2.93)	0.7	993	-0.85 (-3.31,1.61)	-0.33 (-3.05,2.39)	2.64 (-0.77,6.05)	0.1
Body Mass Index (kg/m2) 69 years
Male	961	0.55 (-0.31,1.41)	0.2	840	-0.74 (-1.71,0.23)	-0.49 (-1.58,0.61)	0.55 (-0.82,1.92)	0.10*
Female	1,047	0.28 (-0.85,1.40)	0.7	934	-0.00 (-1.12,1.12)	0.84 (-0.39,2.07)	1.80 (0.23,3.35)	0.004
Waist Circumference (cm) 69 years
Male	961	0.89 (-1.44,3.22)	0.5	840	-1.06 (-3.65,1.53)	-0.45 (-3.38,2.47)	3.63 (-0.04,7.31)	0.03*
Female	1,053	-0.25 (-2.98,2.47)	0.9	939	0.34 (-2.38,3.06)	2.25 (-0.75,5.24)	3.42 (-0.38,7.22)	0.02

*Test for heterogeneity across groups when test for deviation from linearity gives P≤0.05

The associations with lean mass index for men and women remained after adjustment for fat mass index. (Table 4).

Table 4

Regression models for lean mass index (kg/m2) adjusted for fat mass index (kg/m2).

Lean Mass Index (kg/m2) Adjusted for Fat Mass Index (kg/m2)		Number of Children					P Value
	N	0	1	2	3	4+
Unadjusted Regression Coefficients
Male	701	Reference Group	0.90 (0.39,1.40)	0.57 (0.21,0.93)	0.72 (0.29,1.15)	0.69 (0.14,1.24)	0.003*
Female	787	Reference Group	-0.10 (-0.48,0.29)	-0.02 (-0.33,0.29)	0.19 (-0.15,0.53)	0.15 (-0.27,0.57)	0.1
Adjusted Regression Coefficients
Male	701	Reference Group	0.87 (0.36,1.38)	0.54 (0.17,0.91)	0.67 (0.24,1.11)	0.75 (0.19,1.31)	0.006*
Female	787	Reference Group	-0.04 (-0.43,0.36)	0.04 (-0.27,0.35)	0.25 (-0.10,0.60)	0.24 (-0.20,0.67)	0.07
With Children Versus No Children
Male	701	Reference Group	0.62 (0.27,0.98)**				0.001
Female	787	Reference Group	0.10 (-0.20,0.40)**				0.5
With At Least 1 Child
Male	607	-	Reference Group	-0.33 (-0.75,0.10)	-0.20 (-0.68,0.28)	-0.12 (-0.71,0.47)	0.9
Female	691	-	Reference Group	0.07 (-0.24,0.40)	0.29 (-0.07,0.64)	0.28 (-0.16,0.72)	0.07

* Test for heterogeneity across groups when test for deviation from linearity gives P≤0.05.

** Combines categories 1 to 4+

BMI and WC

BMI at age 60–64 years increased with increasing parity, with the highest regression coefficients among women with 3 or 4+ children. For WC the trend was weaker, although those with 4+ children had the highest mean (Table 2). There is no evidence of a difference between those with children and those without for either BMI and WC with weak increasing trends among parous women. Among men, no evidence of trends across groups were observed, although those with 4+ children had higher mean BMI and WC than other groups. No significant sex interactions were observed for either BMI or WC (Table 2). Linear associations for both outcomes are attenuated in women following adjustment, but remain significant for BMI (p = 0.05, Table 2).

The trends of increasing BMI and WC with increasing number of children in women were stronger at age 69 years than at age 60–64 years. In men, associations were also observed at 69 years with mean levels highest in those with 4+ children (Table 2). Associations were slightly attenuated in adjusted models.

The mean BMI and WC at 60–64 years were lower in men and women who had a DXA scan compared with those who did not (two tailed t-test, P<0.001). Those who had a DXA scan were also of higher socioeconomic position and education level, had a higher childhood cognitive score and were less likely to be a current smoker than those who did not attend. Regression analysis showed weaker associations between number of children and both BMI and WC in the unadjusted and adjusted model for those that attended the DXA compared to the whole sample (S1 Table).

Discussion

There was no clear evidence of a consistent association between number of children and body composition in early old age in either men or women and no evidence of clear sex differences in associations. Lean mass index was higher in men with children than those without children, while men with 4+ children had higher fat-lean ratio than others. In women, there was weak evidence of increasing lean mass index with increasing number of children and higher android:gynoid fat mass ratio in those with compared to without children. Both BMI and WC at ages 60–64 years and, more clearly, at 69 years increased with increasing numbers of children in women, whilst highest mean BMI and WC was observed among men with 4+ children, particularly at 69 years.

Limited previous research has investigated the association between number of children and body composition, particularly in men or at older ages. A study carried out in Brazil investigating the effect of parity on fat mass measure by DXA at mean age of 30.2, similarly found little evidence of an association with number of children in both men and women [25]. The weak association between having children and higher android:gynoid fat mass ratio could be consistent with existing research showing increases in visceral fat measured using computed tomography (CT) scans, although not fat mass measured by DXA, among 14 women who had a child within a 5 year interval compared with 108 who did not [22]. Another small study (n = 170) in an age-heterogeneous group of women from 18–70 years found that after adjustment, increasing number of children was associated with increasing intra-abdominal adipose tissue but not with total fat mass [21]. Other research in a larger sample of women from Korea with body composition measured at older ages (45 years and older) has shown younger age at first birth (those with younger age at first birth tend to have a higher number of children) to be associated with higher total and truncal fat mass measured using DXA [23].

Our finding of associations, particularly among women, between number of children and BMI are more consistent with previous findings. In previous research carried out using the NSHD cohort followed up to 53 years, increasing number of children was found to be associated with higher BMI from age 36 to 53 years in women. In men, there was no difference in BMI according to number of children at age 36 years, but those who had at least one child exhibited faster rates of BMI increase and higher BMI at the age of 53 years than those without children [1]. Together with our current findings, this suggests that the association is established during child rearing and then tracks across adult life and into early old age. A study by Lawlor et al (2002) also found linear associations in men and women for both BMI and WHR, with stronger associations observed in women and attenuation to the null for WHR in men once adjusted for lifestyle risk factors and SEP, but no attenuation in women. Two US studies found no evidence that number of children increased weight gain more in women than in men [9, 10] with one even finding greater weight gain among men when stratified according to SEP and race [10].

To our knowledge no other studies have considered associations with lean mass, the outcome where we observed the strongest associations. That we found higher mean lean mass index in men who had children compared to those who did not and, more weakly, increasing lean mass index with increasing parity in women could indicate a degree of selection into parenthood, for example, according to health or social factors. This effect could be greater for men, and is consistent with previous research in NSHD finding that physical performance, including grip strength which is an indicator of muscle strength, was lower among unmarried and childless men compared to others [28]. An alternative explanation for the greater muscle mass among those with children is that with increasing fat mass, muscle mass also increases to accommodate carrying the extra weight [27], although the associations remained once models were adjusted for fat mass index. It is possible that this selection into parenthood could simultaneously occur with a second process of accumulation of fat associated with high numbers of children given the higher fat: lean mass ratio among men with 4+ children.

Our findings suggest that for women the effects of pregnancy on body composition in older age are likely to be small. While the similar associations in men and women for BMI and WC, which are somewhat attenuated after adjustment, instead suggest that the changing social contexts associated with parenthood for both sexes are likely to be key in explaining those associations. Our research extends previous findings by showing that such associations extend into old age, implying the social contexts associated with parenthood may continue to have an effect into early old age. Health behaviours which might be related to weight and body composition which are established as a result of family life, and may indeed change at different stages of family life, may continue to track over the rest of the life course. It has been found that the diets of parents in western cultures contain more saturated fat than non-parents [14, 15, 30] and an increased intake of sugar sweetened beverages and total energy [30]. Differences between mothers and non-mothers in terms of diet quality have also been demonstrated in Japan, with changes to macronutrients and mineral intake, in particular protein, carbohydrates and Sodium [13]. Physical activity has also been identified as being lower in parents compared to non-parents [30-34]. Some research has found that inactivity is associated with increased parity in women [33] whilst others have found that motherhood itself, as opposed to number of children, is key in explaining reduced levels of physical activity [32].

Barriers to maintaining physical activity such as guilt, family responsibilities, lack of support, scheduling constraints, and work have been identified as reasons for the difference in activity patterns between parents and non-parents [35]. Some studies have even found that identities formed around being a parent contradict positive health behaviours [36] as there is a shift from prioritising individual needs to that of their children and family [35, 36]. Where parents have been able to maintain physical activity after childbirth, exercise levels before and during pregnancy, exercise self-efficacy [37, 38] and the ability to develop strategies to incorporate exercise into family life [35] have been highlighted. Additional possible behaviour changes in parents relate to smoking and drinking, with research focusing solely on mothers finding both increased levels of smoking [31, 33] and drinking [33] compared to non-mothers.

There are some limitations with the current study. Firstly, body composition was measured only at age 60–64 years meaning that associations with the number of children could not be assessed closer to age of childbearing in order to see whether the associations change with increasing age. Moreover, the use of DXA scans do not provide detailed insight to levels of visceral compared to subcutaneous fat which may be more relevant for pregnancy. However, research has shown that DXA is a suitable method for assessing body composition, with high correlations between android fat measured by DXA scans and visceral fat measured by CT scans [39, 40]. Additionally if the effects of number of children on outcomes are small, and greatest at high parity, our study may lack power to investigate the effect of high numbers of children due to few individuals with 4 or more children. However, the study had the power to observe associations with BMI.

There is possibility of selection bias through loss to follow up if those individuals who have dropped out of the study or died differ in their relationship between number of children and the outcome variables. Further, the sample who attended clinics for DXA assessment have previously been shown to be healthier than those who participated in home visits (and who had basic anthropometric measures but not DXA) [29]. The sensitivity analysis provided some evidence to suggest that associations may be weaker in the sample who attended clinics as the association with BMI was weaker in this group.

Our findings are from a sample born in Britain 1946 [41] and although the cohort is nationally representative of the population at the time of recruitment, it does not include immigrants into Britain of the same age and may not be representative of later generations. Research suggests that more recent cohorts are experiencing the obesogenic environment earlier in the life course [42] and childbearing patterns have changed such that on the whole couples have later and fewer births, therefore associations may vary between cohorts. Nevertheless, adjustment of social factors relevant to the cohort under study would be expected to result in a consistent biological effect of pregnancy, if one exists.

The current research is, to the best of our knowledge, unique in that it uses DXA imaging to measure body composition, rather than relying on more general anthropometric measures which cannot distinguish fat from lean mass, in a large western sample and which includes men as well as women. Using DXA imaging alongside anthropometric measures such as BMI and WC is a strength when looking at an older population as changes associated with ageing, such as kyphosis and muscle loss, may result in BMI being an inaccurate measure of fatness in older individuals [43].

Our findings suggest little evidence of a biological effect of pregnancy in women on body composition in older age. Social pathways, contexts and behaviours associated with parenthood are more likely to explain the persisting associations observed with more general measures of adiposity, BMI and WC, and selection into parenthood may explain associations with lean mass, particularly in men. As parents have been recognised as important in the development of health behaviours, such as physical activity engagement [44], in their children, these pathways could also have knock on effect on children’s propensity to overweight. Therefore, public health initiatives should focus on adapting the social contexts of parenthood and the health behaviours in families, as opposed to focusing on the timing and processes of pregnancy.

Sensitivity Analysis—Unadjusted and adjusted regressions for anthropometric measures aged 60–64 comparing full sample (those who had a home visit as well as those who attended the CRF) with only those who attended the CRF and thus also had a DXA scan.

(DOCX)

Click here for additional data file.