Association Between Mode of Delivery and Body Mass Index at 4-5 Years in White British and Pakistani Children: the Born in Bradford Birth Cohort

DOI: https://doi.org/10.21203/rs.3.rs-117961/v1

Abstract

Background Globally, it is becoming more common for pregnant women to deliver by caesarean section (CS). In 2012, 21% of births in the UK were CS, surpassing the recommended prevalence of 15%. Concerns have been raised regarding potential unknown consequences of this mode of delivery.

Childhood adiposity is also an increasing concern. Previous research provides inconsistent conclusions on the association between CS and childhood adiposity. More studies are needed to investigate the consequences of CS in different populations and ethnicities. Therefore, this study investigates the association between mode of delivery and BMI, in children of 4-5 years and if this differs between White British (WB) and Pakistani ethnicities, in Bradford UK.

Methods Data were obtained from the Born in Bradford (BiB) cohort, which recruited pregnant women at the Bradford Royal Infirmary, between 2007-2010. For these analyses, a sub-sample (n=6410) of the BiB cohort (n=13858) was used.

Linear regression models determined the association between mode of delivery (vaginal or CS) and BMI z-scores at 4-5 years. Children were categorised as underweight/healthy weight, overweight and obese, and logistic regression models determined the odds of adiposity. Effect modification by ethnicity was also explored.

Results Multivariable analysis found no evidence for a difference in BMI z-score between children of CS and vaginal delivery (0.005 kg/m2, 95% CI= -0.062–0.072, p=0.88). Neither was there evidence of CS affecting the odds of being overweight (OR=1.05, 95% CI=0.86–1.28, p=0.65), or obese (OR=0.98, 95% CI=0.74–1.29, p=0.87). There was no evidence that ethnicity was an effect modifier of these associations (p=0.97).

Conclusion Having CS, compared to a vaginal delivery, was not associated with greater adiposity in children of 4-5 years in this population. Concerns over CS increasing adiposity in children are not supported by the findings reported here using the BiB study population, of both WB and Pakistani families. 

Background

Delivery by caesarean section (CS) is increasing globally. Using data from 150 countries, longitudinal analysis suggests that CS represent 18.6% of all births (1). In the UK, CS rates rose from 11–21%, between 1988 and 2012 (2).

CS is increasing due to research suggesting protective effects against fetal death (3) and to avoid adverse impacts of macrosomia in obese and diabetic mothers (4). Also CS is sometimes perceived as more convenient, less painful and more profitable for private hospitals (4).

However, the rise in CS has aroused alarm due to the lack of knowledge on the short- and long-term risks. With this uncertainty in mind, the World Health Organisation recommended that the CS rate should not surpass 15% (5). Some evidence suggests those who have undergone CS have just over twice the odds of severe maternal morbidity, compared to those experiencing vaginal deliveries (3). CS has also been associated with other complications, such as a higher risk of immune and metabolic disorders in children (6), and offspring overweight and obesity (710). The latter complication will be investigated in this report.

Overweight and obesity in England was prevalent in 22.6% of children aged 4–5 years old, in 2018/2019; more specifically 23.1% in WB children and 19.9% in Pakistani children (11).

In the first six months of life, the colonisation and diversity of gut microbiota is associated with the mode of delivery (12). Those born vaginally have a higher abundance of Bifidobacteria and Bacteroides then those born by CS (12). These bacteria genera have a protective effect against being overweight as they are well equipped to obtain nutrients from breast milk oligosaccharides (13). Additionally, there is evidence to suggest the gut microbiota of a child born by CS is more abundant in Staphylococcus aureus, which has been associated with the development of obesity (13, 14). It is important to note that guidelines endorse the use of prophylactic antibiotics for women undergoing CS, to prevent wound infection (15).

However, in children over the age of six months, there is very weak evidence of an association between mode of delivery and gut microbiota (12, 16), suggesting that the protective effect of vaginal deliveries against adiposity attenuates through early childhood. This conflicts with some evidence of an association between mode of delivery and BMI found in adult life (17).

Previous research presents mixed results. One systematic review concluded children delivered by CS had higher odds of being overweight or obese at 0–8 years (pooled odds ratio from 10 studies = 1.32, 95% CI = 1.15–1.51) (7). Another review also determined CS children to be at higher risk of being obese at 2–18 years (pooled risk ratio from 19 studies = 1.34, 95% CI = 1.18–1.51) (8).

Nine other studies provide evidence to suggest delivering by CS increases the risk of child adiposity (1826). However, three studies have found no evidence of differences in child BMI between CS and vaginal deliveries (2729).

A further search found one study to address the effects of ethnicity on the primary association of mode of delivery and child adiposity. The study found differing race-specific effects of CS with body size at 2 years between African American and non-African American mothers (25). In children of African American mothers, CS was associated with a significantly higher odds of obesity, whereas no association was found in children of non-African American mothers (25). It was suggested that ethnic differences in the developing gut microbiome or epigenetic structure, could be the cause of the effect modification (25).

Rationale For This Study

There is limited published research on the direct association of mode of delivery and child BMI, at the age of children starting school (4–5 years old). This association has not been investigated in UK South Asian mothers compared to WB mothers.

The aim of this study was to determine if there is any association between mode of delivery (CS and vaginal delivery) and BMI at 4–5 years of age, in the BiB cohort, and if this differs between ethnicities (WB and Pakistani).

Methods

Study design

Born in Bradford (BiB) is a longitudinal multi-ethnic birth cohort study. BiB aims to investigate parent and child wellbeing by examining physiological, environmental and genetic factors in the City of Bradford (30). Bradford is situated in the north of England; it is ethnically diverse and has high levels of socio-economic deprivation. BiB recruitment occurred from September 2007 to December 2010. Women who attended the Bradford Royal Infirmary at 26-28 weeks gestation for a routine glucose-tolerance test, which is offered to all women booked to give birth in Bradford, were invited to join the study and 87% of those approached agreed to participate. The BiB population is broadly representative of the maternal population in Bradford (30). 12453 pregnant women gave consent to be involved with the study, a total of 13776 pregnancies were reported, which gave rise to 13858 children. A baseline questionnaire was conducted at recruitment. The mother self-reported most variables, including ethnicity, socioeconomic indicators, alcohol-related and smoking habits. School nurse teams took anthropometric measurements of children in reception class (aged 4-5) as part of the National Child Measurement Programme (NCMP). Anthropometric measurements were also obtained from Primary Care and Child Health Records at this age when NCMP data were missing (82.2% from NCMP, 10.2% from Primary Care, 7.6% from Child Health Records). A subgroup of mothers was followed up for data on breastfeeding, at 6, 12, 18, 24 and 36 months. This subgroup was part of BiB1000, a nested cohort of the BiB prospective birth cohort (31).

Study population

From the total 13858 children enrolled in BiB, 7448 were excluded due to not meeting the inclusion criteria of having mode of delivery, BMI at 4-5 years, singleton birth and being WB or Pakistani. Mothers self-reporting an ethnicity other than WB or Pakistani were excluded as they were a very heterogeneous ethnic group (429 = Indian, 288 = White Other, 253 = Bangladeshi, 226 = Black, 108 = Mixed-White and Black, 61 = Mixed-White and South Asian, 309 = other). This left a sample size of 6410 children (Figure 1). The original BiB cohort had a similar distribution of child BMI and mode of delivery to the study population.

Exposure (mode of delivery)

Exposed participants were children who were delivered by elective or emergency CS. Unexposed participants were children who were delivered vaginally; including normal, forceps and ventouse extraction deliveries.

Outcome (child BMI)

BMI values, recorded by school nurse teams, were transformed to a standardised measure (z-scores). The z-scores were calculated using the LMS method. This is prepared via an Excel spreadsheet, which can be obtained online for free (32). The LMS growth application includes access to a 1990 UK (UK90) reference population. Using this reference, each individual is assigned a z-score which adjusts for age, sex and the BMI distribution for skewness (33). The UK90 reference group is recommended for population monitoring and clinical assessment in children aged four years and over (34). It serves as an anchor for comparison; it is used by the National Child Measurement Programme and has been used for other BiB studies (35). Children with BMI z-scores above the 85th percentile were classified as overweight, and those above the 95th percentile as obese (36).

Sample size calculations

Sample size calculations were conducted in OpenEpi. This study had a power of 99%, determined from a post-hoc power calculation using the parameters from this study (vaginal to CS ratio of 3.71 and prevalence of childhood overweight or obesity in vaginal births at 14.8%), and the odds ratio from a previous study (odds ratio of 2.10 (95% CI 1.36-3.23) of obesity in children aged 7 years, by CS compared to vaginal birth) (9).

Statistical analysis

All analyses were conducted in Stata/IC 15.1. Figures were produced using RStudio version 1.3.1056. Variables that had good evidence (chi-squared tests, ANOVA, and judging correlation to have approximately p<0.05) to suggest they had an association with both mode of delivery and z-scores, as well as not being on the causal pathway, met the criteria to be potential confounders.

Multivariable analysis

The forward selection approach was used to create regression models. Potential confounders were added individually according to their effect size. The covariate was retained in the model if there was an appreciable (10%) difference in effect size of mode of delivery on z-score.

The final multivariable regression model assessed the association between mode of delivery and BMI z-score. Preliminary analysis confirmed the assumptions of the regression were met; z-scores were normally distributed and lacked collinearity. BMI z-scores were also categorised and logistic regression models performed to obtain odd ratios for being overweight and obese in children delivered by CS.

Effect modification

Potential effect modification was judged by stratifying the final model by ethnicity to observe the separate association of mode of delivery on z-score in WB and Pakistani children. A formal test for effect modification was also conducted; a likelihood ratio test compared the final model with a model which also included an interaction term between mode of delivery and ethnicity.

Results

Descriptive results

Table 1 and 2 summarise the baseline characteristics of the study population stratified by mode of delivery and ethnicity. In this study, 21.3% (n=1361) of babies were delivered by CS.  There were more Pakistani mothers (54.6%, n=3502) than WB mothers (45.4%, n=2908). Amongst Pakistani mothers, 19.8% had CS, whereas the CS prevalence among WB mothers was higher at 23.0%. Most children were underweight/ healthy (84.7%), and fewer were overweight (10.0%) or obese (5.2%) (Figure 2). The mean BMI z-score was slightly higher in CS deliveries (0.32) than vaginal deliveries (0.22) at age 4-5 years. Furthermore, the mean BMI z-score was higher in WB children (0.43) than Pakistani children (0.08).

Table 1. Baseline characteristics, of the study population, stratified by mode of delivery.

   

All

Vaginal

(n= 5049)

Caesarean

(n= 1361)

P value

   

n

%

n

%

n

%

Child BMI categorised

(n= 6410)

 

Underweight/

Healthy

5432

84.74

4303

85.22

1129

82.95

0.039

 

 

 

Overweight

643

10.03

490

9.70

153

11.24

Obese

335

5.23

256

5.07

79

5.80

Child BMI z-score

(n= 6410)

Mean

0.24

0.22

0.32

0.003

SD

1.12

1.11

1.15

Ethnicity

(n= 6410)

White British

2908

45.37

2239

44.35

669

49.16

0.002

 

Pakistani

3502

54.63

2810

55.65

692

50.84

Maternal age (years)

(n= 6410)

Mean

27.47

27.12

28.77

<0.0001

 

SD

5.61

5.53

5.73

Maternal BMI (kg/m2)

(n= 6135)

Mean

26.09

25.67

27.67

<0.0001

 

SD

5.69

5.46

6.26

Maternal education

(n= 6395)

 

 

 

 

<5 GCSE equivalent

1502

23.49

1230

24.41

272

20.04

<0.001

 

 

 

 

 5 GCSE equivalent

2155

33.70

1718

34.10

437

32.20

A-level equivalent

876

13.70

695

13.80

181

13.34

Higher than A-level

1420

22.20

1058

21.00

362

26.68

Foreign unknown/other

442

6.91

337

6.69

105

7.74

Maternal job status

(n= 6400) 

Currently employed

2707

42.30

2029

40.24

678

49.93

<0.001

 

 

Previously employed

1885

29.45

1517

30.09

368

27.10

Never employed

1808

28.25

1496

29.67

312

22.97

Maternal house tenure

(n= 6396)

 

 

 

 

Owns outright

1014

15.85

807

16.02

207

15.25

0.055

 

 

 

 

Mortgage

3236

50.59

2518

49.97

718

52.91

Private landlord

964

15.07

792

15.72

172

12.68

Social housing

667

10.43

522

10.36

145

10.69

Rent free/other

515

8.05

400

7.94

115

8.47

Maternal benefits received

(n= 6385)

Yes

2733

42.80

2237

44.51

496

36.50

<0.001

 

No

3652

57.20

2789

55.49

863

63.50

Maternal drinking of alcohol during pregnancy or 3 months before

(n= 6396)

Yes

2065

32.29

1570

31.17

495

36.42

0.001

 

 

No

4328

67.67

3465

68.79

863

63.50

Don't remember

3

0.05

2

0.04

1

0.07

Maternal smoking during pregnancy

(n= 6396)

Yes

1066

16.67

856

16.99

210

15.45

0.176

 

No

5330

83.33

4181

83.01

1149

84.55

Breastfeeding

(n= 971)

Ever

701

72.19

557

72.53

144

70.94

0.653

 

Never

270

27.81

211

27.47

59

29.06

Parity

(n= 6137) 

Primiparous

2388

38.91

1802

37.22

586

45.22

<0.001

 

Multiparous

3749

61.09

3039

62.78

710

54.78

Maternal gestational diabetes

(n= 6404)

Yes

489

7.64

342

6.78

147

10.83

<0.001

 

No

5915

92.36

4705

93.22

1210

89.17

Child gender

(n= 6410)

Male

3266

50.95

2548

50.47

718

52.76

0.134

 

Female

3144

49.50

2501

49.53

643

47.24

Child birthweight (g)

(n= 6410)

Mean

3238.10

3244.10

3215.84

0.089

 

SD

544.78

503.61

675.73

Gestational period (days)

(n= 6410)

Mean

276.86

277.67

273.90

<0.0001

 

SD

12.22

11.09

15.34

P values to provide the level of statistical evidence on the difference between mode of delivery; obtained from chi-squared tests or ANOVA, where appropriate. Abbreviations: n, sample size; BMI, body mass index; SD, standard deviation.

 

Table 2. Baseline characteristics, of the study population, stratified by ethnicity.

   

All

White British

(n= 2908)

Pakistani

(n= 3502)

P value

   

n

%

n

%

n

%

Child BMI categorised

(n= 6410)

 

Underweight/

Normal

5432

84.74

2435

83.73

2997

85.58

0.041

 

Overweight

643

10.03

335

11.52

308

8.79

Obese

335

5.23

138

4.75

197

5.63

Child BMI z-score

(n= 6410)

Mean

0.24

0.43

0.08

<0.0001

 

SD

1.12

0.97

1.21

Maternal age (years)

(n= 6410)

Mean

27.47

27.03

27.84

<0.0001

 

SD

5.61

6.09

5.15

Maternal BMI (kg/m2)

(n= 6135)

Mean

26.09

26.80

25.50

<0.0001

 

SD

5.69

5.97

5.39

Maternal education

(n= 6395)

 

 

 

 

<5 GCSE equivalent

1502

23.49

577

19.87

925

26.50

<0.001

 

 

 

 

 5 GCSE equivalent

2155

33.70

1028

35.40

1127

32.28

A-level equivalent

876

13.70

458

15.77

418

11.97

Higher than A-level

1420

22.20

560

19.28

860

24.63

Foreign unknown/other

442

6.91

281

9.68

161

4.61

Maternal job status

(n= 6400) 

Currently employed

2707

42.30

1924

66.21

783

22.41

<0.001

 

 

Previously employed

1885

29.45

775

26.67

1110

31.77

Never employed

1808

28.25

207

7.12

1601

45.82

Maternal house tenure

(n= 6396)

 

 

 

 

Owns outright

1014

15.85

119

4.10

895

25.63

<0.001

 

 

 

 

Mortgage

3236

50.59

1519

52.31

1717

49.17

Private landlord

964

15.07

654

22.52

310

8.88

Social housing

667

10.43

451

15.53

216

6.19

Rent free/other

515

8.05

161

5.54

354

10.14

Maternal benefits received

(n= 6385)

Yes

2733

42.80

1063

36.73

1670

47.84

<0.001

 

No

3652

57.20

1831

63.27

1821

52.16

Maternal drinking of alcohol during pregnancy or 3 months before

(n= 6396)

Yes

2065

32.29

2055

70.72

10

0.29

<0.001

 

 

No

4328

67.67

848

29.18

3480

99.71

Don't remember

3

0.05

3

0.10

0

0.00

Maternal smoking during pregnancy (n= 6396)

Yes

1066

16.67

961

33.08

105

3.01

<0.001

 

No

5330

83.33

1944

66.92

3386

96.99

Breastfeeding

(n= 971)

Yes

701

72.19

270

63.08

431

79.37

<0.001

 

No

270

27.81

158

36.92

112

20.63

Parity

(n= 6137) 

Primiparous

2388

38.91

1351

48.11

1037

31.15

<0.001

 

Multiparous

3749

61.09

1457

51.89

2292

68.85

Maternal gestational diabetes

(n= 6404)

Yes

489

7.64

141

4.85

348

9.95

<0.001

 

No

5915

92.36

2764

95.15

3151

90.05

Child gender

(n= 6410)

Male

3266

50.95

1487

51.13

1779

50.80

0.789

 

Female

3144

49.05

1421

48.87

1723

49.20

Child birthweight (g)

(n= 6410)

Mean

3238.10

3357.60

3138.86

<0.0001

 

SD

544.78

550.53

519.56

Gestational period (days)

(n= 6410)

Mean

276.86

277.63

276.22

<0.0001

 

SD

12.22

12.51

11.94

P values to provide the level of statistical evidence on the difference between White British and Pakistani ethnic groups; obtained from chi-squared tests or ANOVA, where appropriate. Abbreviations: n, sample size; BMI, body mass index; SD, standard deviation.


Mother’s booking BMI

Mothers who gave birth by CS had a higher mean BMI at booking of 2 kg/m2 compared to mothers giving birth vaginally. The mean BMI of all mothers in this study population was 26.09 kg/m2, which is categorised as overweight (37). WB mothers had a slightly higher BMI (+1.30 kg/m2) than Pakistani women.

Family sociodemographic factors

Mothers who had CS often had achieved a higher level of education than mothers with vaginal deliveries (26.7% of CS mothers, 21.0% of vaginal mothers). More Pakistani women achieved higher than A-level qualifications compared to WB women (24.6 vs 19.3%).

Most mothers were not currently employed (67.7%). Current employment was more common in women who had CS births, compared to those with vaginal births (49.9% and 40.2%, respectively) and also more common in WB women (66.2%), compared to Pakistani women (22.4%).  Additionally, a higher proportion of mothers having vaginal deliveries (44.5%) received benefits than those having CS deliveries (36.5%).  Receiving benefits was more common in Pakistani mothers (47.8%) than WB mothers (36.7%).

Gestational factors

The mean age of mothers who had a CS was 28.8 years old, which was 1.7 years older than those who gave birth vaginally. Also, mothers having vaginal deliveries were more likely to be multiparous (62.8%) compared to mothers having CS deliveries (54.8%). The difference in gestation period between CS and vaginal deliveries was minimal.

Only 10 out of 3490 Pakistani mothers (0.3%) drank alcohol during pregnancy or 3 months before, whereas 70.7% of WB mothers reported alcohol consumption. Further to this, mothers having CS were marginally more likely to have drunk alcohol during pregnancy or 3 months before (36.4% of CS mothers and 31.2% of vaginal mothers). More WB women reportedly smoked during pregnancy than Pakistani women (33.1% of WB women, 3.0% of Pakistani women).

Gestational diabetes was slightly more prevalent in mothers having CS (10.8%) compared to those having vaginal deliveries (6.8%). Additionally, prevalence was higher in Pakistani mothers than WB mothers (10.0% and 4.9% respectively).

Child factors

Children born by CS had a mean birthweight 28.3 g lower than vaginal births. Irrespective of mode of delivery, children with WB mothers had a higher mean birthweight than those with Pakistani mothers (3357.6 g and 3138.9 g, respectively).

Missing data

Approximately 85% of the study population have missing data on breastfeeding (n= 5439 missing). Maternal parity (n= 273 missing) and maternal BMI (n= 275 missing) are also missing data, but on about 4% of the study population. There was no evidence of selection bias based on the distribution of missing data.      

Multivariable analysis

The unadjusted linear regression calculated the predicted difference in z-score between mode of delivery, the z-score being higher with CS (n=6410, difference= 0.103; 95% CI= 0.035–0.170) (Figure 3). The adjusted model calculated the predicted mean difference in z-score between mode of delivery, controlling for all factors which met the confounding criteria (ethnicity, maternal BMI, maternal job status and maternal drinking of alcohol during pregnancy or 3 months before), the z-score being higher with CS (n= 6115, difference= 0.005; 95% CI= -0.062–0.072) (Figure 3).

The proportion of variance in z-scores explained by the mode of delivery was 7.49% for the adjusted model (adjusted R2 value). The F ratio (71.76) shows how much variability the model can explain relative to how much it cannot explain. The standard error (0.034) did not differ between unadjusted and adjusted models, suggesting an absence of collinearity.

The adjusted logistic regression models with the outcome of overweight and obesity obtained odds ratios and confidence intervals of no strong statistical support for a difference in odds between mode of delivery (Figure 4). 

Effect modification

Stratified analysis suggests the adjusted association of mode of delivery on BMI z-score is similar, irrespective of ethnicity (Figure 3). Additionally, there was no evidence of effect modification from the likelihood ratio test (p= 0.97). When a test for effect modification was performed on the categorised z-scores, similar results were obtained. There was weak evidence of effect modification by ethnicity on the association between mode of delivery and overweight (p=0.14), and no evidence for effect modification by ethnicity on the association between mode of delivery and obesity (p=0.79).

Discussion

In this cohort study, it was found that undergoing a CS was not associated with an increased risk of overweight and obesity in children, and there was no difference between ethnic groups.

Mothers who had undergone a CS were generally of higher socio-economic status than those who had vaginal deliveries; CS women were more educated, more likely to be currently employed, more likely to have a stable housing situation and less likely to be receiving benefits. Also, mothers who experienced CS had baseline characteristics to suggest they had poorer health than mothers giving birth vaginally; CS mothers had a higher mean BMI, were more likely to drink alcohol during pregnancy or 3 months before and more prevalent gestational diabetes. The distribution of alcohol drinking varies vastly between ethnicities, this is most likely due to religious beliefs (30). This explains the very low prevalence of alcohol drinking and avoidance of smoking amongst Pakistani women.

The linear regression for the adjusted model offers no evidence for a difference in BMI z-score between children born via CS and vaginal deliveries. The low adjusted R2 value suggests there are other variables which have an influence on the primary association. The adjusted logistic regression models also suggest no evidence for children delivered by CS having different odds of being overweight or obese, compared to children of vaginal deliveries.

The stratified analysis and formal test for effect modification both indicate there is no evidence that the association between mode of delivery and children’s BMI z-scores varied by ethnicity.

As discussed in the introduction, previous studies have varied interpretations. Two leading systematic reviews suggest there is evidence that CS increases child BMI (7, 8). However, there were several studies which found no statistical association between mode of delivery and child BMI. The findings from this project are compatible with the later studies mentioned.

Two out of three studies conducted in the UK concluded there was no ‘statistical significant’ difference in risk of childhood overweight or obesity between modes of delivery, at 3 years old (38) and 5 years old (27).

Furthermore, maternal BMI explained most of the observed association in this study and was hence the main confounding factor. All previous studies looking at the association between mode of delivery and BMI, also adjusted for this factor. Additionally, none of the previous studies identified any confounders that had a biologically plausible reason to include them in the unadjusted model.

However, this study also differs with previous research. The other UK study found that CS increased the odds of being overweight or obese, at 7 years old (18). This was a study which used data from the Avon Longitudinal Study of Parents and Children (ALSPAC); participants were recruited from the Avon area if they were born in 1991–1992.

Several confounders (child gender, gestational factors and child feeding patterns) were adjusted for in the ALSPAC study but did not meet the confounding criteria to be adjusted for in this study. There were also inconsistencies with other factors adjusted for in this study compared to previous studies, such as not adjusting for antibiotics during pregnancy (19). Different factors could have met the confounding criteria in previous studies due to their population type, for example, by having a different BMI distribution as the children were leaner.

Adjusting for ethnicity was not seen in previous research in the UK. Due to the large proportion of Pakistani women in this BiB study, there was sufficient power to investigate differences between WB and Pakistani ethnic groups, whereas this would not be possible in studies like ALSPAC. As previous studies did not adjust for ethnicity, other variables could have acted as confounders. Overall, the differences in study design, study population and confounding adjustments could explain the inconsistent conclusions reached.

The large sample size used in this study allowed sufficient power to identify any meaningful differences in association between BMI z-scores of two different modes of delivery. Additionally, consistent statistical methodology with previous studies was used and there was minimal recruitment bias due to the BiB study having a high recruitment rate of 87%.

As CS and vaginal deliveries are very different procedures, in theory, there was no opportunity for this to be incorrectly recorded. Hence no information bias, in the form of differential misclassification, should have occurred. Furthermore, observer bias would not arise when recording the child’s BMI, as nurses taking anthropometric measurements at ages 4–5 were blinded to information regarding the child’s mode of delivery.

There is evidence to suggest BMI measurements systematically underestimate childhood adiposity (39). This has also been specifically investigated in South Asians with evidence to suggest that BMI additionally appears to underestimate adiposity in this ethnic group. Despite South Asians being generally smaller and lighter, they seem to have greater relative fatness compared to white European populations (40).

Most of the data on covariates were collected in the baseline questionnaire, completed by the mother. As data were self-reported, information bias in the form of differential misclassification could have occurred which would tend the results to overestimate or underestimate the true association. An example would be smoking as this is a likely factor to be underreported. Underreporting could underestimate the association between mode of delivery and smoking, which would have led to it not being adjusted for in the final analysis.

There may be residual confounding which is obscuring the true effect of mode of delivery on child BMI. The low adjusted R2 value implies other factors could have an influence on the association, therefore suggesting factors which were not included in the analysis explained some of the association. These could be factors such as amount of exercise, feeding pattern of the child or breastfeeding. Adding the effect of breastfeeding on childhood adiposity as a covariable in the multivariable analysis would have been desirable, but this was not possible due to the missing data on breastfeeding. Approximately 85% of the study population had missing data on breastfeeding due to the data being collected in a subgroup of women who participated in the BiB1000 as described in the study design section.

Additionally, the generalisability of the results is limited to the BiB population, however, similarities with other multi-ethnic UK cities are likely to be seen.

Conclusion

Overall CS was not associated with an increased risk of overweight or obesity in children aged 4–5 years in Bradford. Neither was there a difference in association seen between White British or Pakistani children. To our knowledge, this is the first study to assess this association, between these ethnicities, at this age.

As CS deliveries are becoming more common globally and health concerns have been raised, the results from this study, combined with similar studies, should be informative to prospective parents and healthcare advisors.

Data collection within the BiB cohort should be continued to provide more reliable estimates of adiposity and to allow investigation at older ages. This will enable examination of whether any association exists at subsequent ages between mode of delivery and later life adiposity.

As there is some uncertainty around how well BMI represents child adiposity, the use of body fat centile curves should be explored instead. To do this, data on fat mass analysed using a DXA scanner would be needed. A DXA scanner is an extremely accurate method for analysing body composition, could be used as a gold standard in the population. The sensitivity and specificity of BMI against the DXA scanner in the population can be obtained and incorporated in the interpretation.

Further research could also investigate if the method of fetus extraction acts as an effect modifier on the association between mode of delivery and BMI. Vaginal deliveries can occur by: natural vaginal birth, forceps assistance or by ventouse techniques. CS can occur by: emergency, elective or semi-elective delivery. These studies would need to have a large sample size to power the subgroup analyses.

Additionally, if the mechanism for any potential association is related to the developing gut microbiota, then more studies could focus on the differences between gut microbiota stratified by mode of delivery and by ethnicity. This is likely to involve a genetic approach when looking at the differences between the WB and Pakistani population (16).

Abbreviations

ALSPAC, Avon Longitudinal Study of Parents and Children; BiB, Born in Bradford; BMI, body mass index; CI, confidence intervals; CS, caesarean section; DXA, Dual-energy X-ray absorptiometry; MoD, mode of delivery; NCMP, National Child Measurement Programme; OR, odds ratio; WB, White British.

Declarations

Ethics approval and consent to participate

All methods were carried out in accordance with relevant guidelines and regulations. Informed consent was obtained from all participants or, if subjects are under 18, from a parent or legal guardian. BiB ethical approval was granted by Bradford Research Ethics Committee (Ref 07/H1302/112). Participating women consented at recruitment to access to their routine primary and secondary health care records. Ethical approval for the current study was received from the Ethics Committee at the London School of Hygiene and Tropical Medicine.

Consent for publication

Not applicable.

Availability of data and materials

Data from the BiB study is available to researchers following approval from the Executive Committee (https://borninbradford.nhs.uk/research/how-to-access-data/)

Competing interests

The authors declare that they have no competing interests.

Funding

Born in Bradford is supported by the Wellcome Trust [101597]; the UK Medical Research Council (MRC) and UK Economic and Social Science Research Council (ESRC) [MR/N024391/1]; British Heart Foundation [CS/16/4/32482]; the National Institute for Health Research under its Applied Research Collaboration Yorkshire and Humber [NIHR200166] and the NIHR Clinical Research Network which provided research delivery support for this study. Lucy Pembrey was supported by the Wellcome Trust (grant numbers: 083521/Z/07/Z, 083521/Z/07/A), the AsthmaPhenotypes Study Grant; European Research Council under the European Union’s Seventh Framework Programme (FP7/2007-2013)/ ERC grant agreement no. 668954 and National Institute for Health Research, Health Technology Assessment grant (no. 16/150/06).

Authors’ contributions

ER was responsible for analysing the data, interpreting the results and writing the first draft. LP had input into design of the study, analysis plan, and interpreting the results. JW and GS conceived the original research idea, were responsible for facilitating data extraction from the Born in Bradford cohort. All authors contributed to further iterations of the draft. All authors approved the final version.

Acknowledgments

Born in Bradford is only possible because of the enthusiasm and commitment of the Children and Parents in BiB. We are grateful to all the participants, health professionals and researchers who have made Born in Bradford happen.

References

  1. Betrán AP, Ye J, Moller A-B, Zhang J, Gülmezoglu AM, Torloni MR. The increasing trend in caesarean section rates: global, regional and national estimates: 1990–2014. PloS One. 2016;11(2):e0148343.
  2. Rose A, Raja EA, Bhattacharya S, Black M. Intervention thresholds and cesarean section rates: A time-trends analysis. Acta Obstet Gynecol Scand. 2018;97(10):1257–1266.
  3. Villar J, Carroli G, Zavaleta N, Donner A, Wojdyla D, Faundes A, et al. Maternal and neonatal individual risks and benefits associated with caesarean delivery: multicentre prospective study. Bmj. 2007;335(7628):1025.
  4. Mi J, Liu F. Rate of caesarean section is alarming in China. The Lancet. 2014;383(9927):1463–1464.
  5. Appropriate technology for birth. Lancet Lond Engl. 1985 Aug 24;2(8452):436–7.
  6. Dominguez-Bello MG, De Jesus-Laboy KM, Shen N, Cox LM, Amir A, Gonzalez A, et al. Partial restoration of the microbiota of cesarean-born infants via vaginal microbial transfer. Nat Med. 2016;22(3):250.
  7. Li HT, Zhou YB, Liu JM. The impact of cesarean section on offspring overweight and obesity: a systematic review and meta-analysis. Int J Obes. 2013;37(7):893–899.
  8. Kuhle S, Tong OS, Woolcott CG. Association between caesarean section and childhood obesity: a systematic review and meta-analysis. Obes Rev. 2015;16(4):295–303.
  9. Huh SY, Rifas-Shiman SL, Zera CA, Edwards JWR, Oken E, Weiss ST, et al. Delivery by caesarean section and risk of obesity in preschool age children: a prospective cohort study. Arch Dis Child. 2012;97(7):610–616.
  10. Mueller NT, Zhang M, Hoyo C, Østbye T, Benjamin-Neelon SE. Does cesarean delivery impact infant weight gain and adiposity over the first year of life? Int J Obes. 2019;43(8):1549–1555.
  11. NCMP and Child Obesity Profile - PHE [Internet]. [cited 2020 Nov 24]. Available from: https://fingertips.phe.org.uk/profile/national-child-measurement-programme/data#page/7/gid/8000011/pat/126/par/E47000002/ati/101/are/E08000016/iid/20601/age/200/sex/4/cid/4/page-options/ine-ao-1_ine-vo-0_ine-yo-1:2014:-1:-1_ine-ct-17_ine-pt-0
  12. Rutayisire E, Huang K, Liu Y, Tao F. The mode of delivery affects the diversity and colonization pattern of the gut microbiota during the first year of infants’ life: a systematic review. BMC Gastroenterol. 2016;16(1):86.
  13. Kalliomäki M, Carmen Collado M, Salminen S, Isolauri E. Early differences in fecal microbiota composition in children may predict overweight. Am J Clin Nutr. 2008;87(3):534–538.
  14. Dominguez-Bello MG, Costello EK, Contreras M, Magris M, Hidalgo G, Fierer N, et al. Delivery mode shapes the acquisition and structure of the initial microbiota across multiple body habitats in newborns. Proc Natl Acad Sci. 2010;107(26):11971–11975.
  15. Gholitabar M, Ullman R, James D, Griffiths M. Caesarean section: summary of updated NICE guidance. Bmj. 2011;343:d7108.
  16. Vinding RK, Sejersen TS, Chawes BL, Bønnelykke K, Buhl T, Bisgaard H, et al. Cesarean delivery and body mass index at 6 months and into childhood. Pediatrics. 2017;139(6).
  17. Darmasseelane K, Hyde MJ, Santhakumaran S, Gale C, Modi N. Mode of delivery and offspring body mass index, overweight and obesity in adult life: a systematic review and meta-analysis. PloS One. 2014;9(2):e87896.
  18. Blustein J, Attina T, Liu M, Ryan AM, Cox LM, Blaser MJ, et al. Association of caesarean delivery with child adiposity from age 6 weeks to 15 years. Int J Obes. 2013;37(7):900–906.
  19. Mueller NT, Whyatt R, Hoepner L, Oberfield S, Dominguez-Bello MG, Widen EM, et al. Prenatal exposure to antibiotics, cesarean section and risk of childhood obesity. Int J Obes. 2015;39(4):665–670.
  20. Portela DS, Vieira TO, Matos SM, de Oliveira NF, Vieira GO. Maternal obesity, environmental factors, cesarean delivery and breastfeeding as determinants of overweight and obesity in children: results from a cohort. BMC Pregnancy Childbirth. 2015;15(1):94.
  21. Veile A, Kramer KL. Childhood body mass is positively associated with cesarean birth in Y ucatec M aya subsistence farmers. Am J Hum Biol. 2017;29(2):e22920.
  22. Yuan C, Gaskins AJ, Blaine AI, Zhang C, Gillman MW, Missmer SA, et al. Association between cesarean birth and risk of obesity in offspring in childhood, adolescence, and early adulthood. JAMA Pediatr. 2016;170(11):e162385–e162385.
  23. Rutayisire E, Wu X, Huang K, Tao S, Chen Y, Tao F. Cesarean section may increase the risk of both overweight and obesity in preschool children. BMC Pregnancy Childbirth. 2016;16(1):338.
  24. Pluymen LP, Smit HA, Wijga AH, Gehring U, De Jongste JC, Van Rossem L. Cesarean delivery, overweight throughout childhood, and blood pressure in adolescence. J Pediatr. 2016;179:111–117.
  25. Cassidy-Bushrow AE, Wegienka G, Havstad S, Levin AM, Lynch SV, Ownby DR, et al. Race-specific association of caesarean-section delivery with body size at age 2 years. Ethn Dis. 2016;26(1):61.
  26. Mueller NT, Mao G, Bennet WL, Hourigan SK, Dominguez-Bello MG, Appel LJ, et al. Does vaginal delivery mitigate or strengthen the intergenerational association of overweight and obesity? Findings from the Boston Birth Cohort. Int J Obes. 2017;41(4):497–501.
  27. Black M, Bhattacharya S, Philip S, Norman JE, McLernon DJ. Planned cesarean delivery at term and adverse outcomes in childhood health. Jama. 2015;314(21):2271–2279.
  28. Carrillo-Larco RM, Miranda JJ, Bernabé-Ortiz A. Delivery by caesarean section and risk of childhood obesity: analysis of a Peruvian prospective cohort. PeerJ. 2015;3:e1046.
  29. Morgen CS, Ängquist L, Baker JL, Andersen A-MN, MIchaelsen KF, Sørensen TI. Prenatal risk factors influencing childhood BMI and overweight independent of birth weight and infancy BMI: a path analysis within the Danish National Birth Cohort. Int J Obes. 2018;42(4):594–602.
  30. Wright J, Small N, Raynor P, Tuffnell D, Bhopal R, Cameron N, et al. Cohort profile: the Born in Bradford multi-ethnic family cohort study. Int J Epidemiol. 2013;42(4):978–991.
  31. Bryant M, Santorelli G, Fairley L, West J, Lawlor DA, Bhopal R, et al. Design and characteristics of a new birth cohort, to study the early origins and ethnic variation of childhood obesity: the BiB1000 study. Longitud Life Course Stud. 2013;4(2):119–135.
  32. LMSgrowth | Health for all Children [Internet]. [cited 2020 Nov 24]. Available from: https://www.healthforallchildren.com/shop-base/shop/software/lmsgrowth/
  33. Cole TJ, Freeman JV, Preece MA. Body mass index reference curves for the UK, 1990. Arch Dis Child. 1995;73(1):25–29.
  34. NHS N. National Obesity Observatory on behalf of the Public Health Observatories in England. A simple guide to classifying body mass index in children, 2011. 2011.
  35. Santorelli G, Petherick ES, Wright J, Wilson B, Samiei H, Cameron N, et al. Developing prediction equations and a mobile phone application to identify infants at risk of obesity. PLoS One. 2013;8(8):e71183.
  36. Reilly JJ, Dorosty AR. Epidemic of obesity in UK children. The Lancet. 1999;354(9193):1874–1875.
  37. NICE (UK C for PHE at, Care (UK NCC for P. Obesity: the prevention, identification, assessment and management of overweight and obesity in adults and children. 2006;
  38. Weng SF, Redsell SA, Nathan D, Swift JA, Yang M, Glazebrook C. Estimating Overweight Risk in Childhood From Predictors During Infancy. Pediatrics. 2013 Aug 1;132(2):e414–21.
  39. McCarthy HD, Ellis SM, Cole TJ. Central overweight and obesity in British youth aged 11–16 years: cross sectional surveys of waist circumference. Bmj. 2003;326(7390):624.
  40. Nightingale CM, Rudnicka AR, Owen CG, Cook DG, Whincup PH. Patterns of body size and adiposity among UK children of South Asian, black African–Caribbean and white European origin: Child Heart And health Study in England (CHASE Study). Int J Epidemiol. 2011;40(1):33–44.