Skip to main content

The impact of cesarean delivery on infant DNA methylation

Abstract

Background

Mounting evidence suggests that cesarean delivery may have a long-lasting effect on infant health. But the underlying mechanisms remain unclear. This study aims to examine whether cesarean delivery on maternal request without any medical indications (CDMR) impacts DNA methylation status in the umbilical cord blood of the infant.

Methods

A cross-sectional study was conducted in Shanghai, China. A total of 70 CDMR and 70 vaginal deliveries (VD) were recruited in 2012. The cord blood DNA methylation status was measured in 30 CDMR and 30 VD newborns using Illumina Infinium Human Methylation 450 K BeadChip. To validate the results, the cord blood DNA methylation status was measured in another 40 CDMR and 40 VD newborns using targeted bisulfite sequencing assay. A total of 497 CpG sites from 40 genes were included in the analysis.

Results

A total of 165 differentially methylated positions (DMPs) exhibited differences in DNA methylation by 10% or more between the CDMR and VD groups, many of which were related to the development of the immune system. Based on the targeted bisulfite sequencing assay, 16 genes (16/22, 72.7%) had higher methylation level in the CDMR group than the VD group. Among them, 5 genes were related to the immune system. After considering the estimation of cell type proportions, there was few significant differences in DNA methylation between CDMR and VD groups.

Conclusions

The DMPs identified between CDMR and VD groups might be largely explained by the cell type proportions. Further studies are needed to examine DNA methylation in each cell type separately.

Peer Review reports

Background

Birth is probably the most dramatic experience in the whole life. After labor for several hours, a totally dependent fetus must rely on his body systems to survive in the outside world. Labor, therefore, may not only be a physical process of birth but also a physiological process to prepare the fetus to survive independently [1]. Pre-labor cesarean delivery (CD) interrupts this natural process and, therefore, might have important health impacts. However, our understanding in this aspect is very limited.

China is one of the countries with the highest CD rates in the world [2]. Moreover, almost half of CDs were due to maternal preference or other social factors instead of clinical reasons [3]. The overuse of CD has raised public health concerns [4]. More and more evidence indicate that CD may have long-term impacts on not only maternal but also child health. For example, epidemiologic studies showed that CD may increase risks of childhood asthma and other allergic disorders [5,6,7], and obesity [8]. But the biological mechanisms remain unclear.

Among possible biological links between CD and childhood disorders, physiological effects of stress caused by strong uterine contraction during labor may play an important role. Normal epigenetic status may also be altered in CDs compared with vaginal deliveries (VD) [9,10,11]. The epigenetic alterations not only change the expression of the functional genome and mediate adaptations to a dynamic environment but also retain its stability for the cells lifetime, even after divisions [12, 13]. Therefore, DNA methylation may be one of the possible mechanisms that explain why the mode of delivery may have long-term effects on the offspring’s health.

Our study aimed to investigate the association between CD and DNA methylation status in the cord blood. To avoid any potential confounding effects of maternal and/or fetal complications on the association, we selected women who had CD due to maternal request without any medical indication (CDMR) as the exposed group and normal VD as the non-exposed group.

Methods

Study population

In 2012, we recruited 70 women who had CDMR and 70 women who had VD from Xinhua Hospital affiliated to the Shanghai Jiao Tong University School of Medicine in Shanghai, China. Women who had CDMR and VD were matched on their age at delivery (within 5 years), pre-pregnancy body mass index (BMI) (within 1 kg/m2), infant birthweight (within 300 g), and gestational weeks at delivery (within 1 weeks). A research nurse introduced the project to the women and a written informed consent was obtained prior to delivery. The ethics approval was obtained from the Ethics Committee of Xinhua Hospital affiliated to the Shanghai Jiao Tong University School of Medicine. The women were asked of a wide array of questions such as age, education, occupation, family income, alcohol use, smoking and exposure to second hand smoke during pregnancy. At birth, umbilical cord blood samples were taken and newborn anthropometric measures were recorded. Detailed obstetric information was extracted from medical records after delivery.

Sample collection and preparation

Blood samples were collected at delivery by double clamping the umbilical cord at the placental end and the infant end. The whole blood was centrifuged at 4 °C (10 min, 1500 g), then buffy coat was collected in BD EDTA tubes (QIAamp DNA Mini Kit; Qiagen, Hilden, Germany) and stored at − 80 °C until DNA extraction.

Genome-wide, locus-specific DNA methylation analysis in cord blood

The 30 CDMR and individually matched 30 VD samples were randomly selected for methylation analysis. Genomic DNA was isolated from 200 μL whole blood for each sample. Then the bisulfite treatment of genomic DNA was performed with an EZ DNA Methylation kit (Zymo, Irvine, CA) in accordance with the protocol. Quantitative DNA methylation were measured using Illumina Infinium Human Methylation 450 K BeadChip (Illumina, San Diego, CA). More than 485,000 cytosine phosphate-guanine (CpG) sites were included in each chip. The arrays were scanned using an Illumina iScan with initial quality review using GenomeStudio software. The level of methylation was expressed as a “beta” value (β-value), ranging from 0 (no cytosine methylation) to 1 (complete cytosine methylation). β-values are reported as percentages. Quality control measures were routinely taken and the details were described elsewhere [14]. CpG sites located at known single-nucleotide polymorphisms and on the X and Y chromosomes (list provided by the manufacturer) were discarded to avoid potential confounding by single-nucleotide polymorphisms and gender.

Differential DNA methylation positions (DMPs) were identified using the GenomeStudio methylation module v1.8 [15]. DMPs were defined as those that exhibited a 10% or greater difference in the DNA methylation between the CDMR and VD groups at p < 0.05 [11]. By using the annotation provided by Illumina, we summarized the DMP information over genes by computing the number of probes for each gene, the number of DMPs, and how many DMPs were hypermethylated or hypomethylated. In addition, gene ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were used to represent methylated genes. To investigate the functional relevance of the selected DMPs, we used the online DAVID bioinformatics resource to test for enrichment in gene ontology biological processes, and we used the Benjamini-Hochberg procedure to control for false discoveries [16].

Moreover, we estimated the proportions of different cell types using FlowSorted.CordBloodCombined.450 k package in R. This package was derived from 263 umbilical blood cell references. Details of this package and the accompanied data are shown elsewhere [17].

Targeted bisulfite sequencing assay

To validate the results of the methylation chip, we used the rest 40 CDMR and individually matched 40 VD samples. Targeted bisulfite sequencing assay DNA extraction and bisulfite conversion were performed as previously described. Considering the genes function and the design of upstream and downstream primers, we selected 40 genes to validate the results of 450 K BeadChips. Then the primers were carefully designed in order to detect them in a panel (Additional file Table S1). The net-PCR was performed firstly to amplify the targeted DNA sequence. Then, the target genes were sequenced by Illumina Hiseq 2000. All samples achieved a mean coverage of 600 base-pairs. More CpG sites around the target CpG sites were found in the validation experiments. Finally, 497 CpG sites were included in the validation assay. The methylation level of each CpG site was calculated as the percentage of the methylated cytosines over total tested cytosines. The average methylation level was calculated using methylation levels of all measured CpG sites within the gene.

Identification of CpG sites associated with diseases

The comparative toxicogenomics database (CTD, http://ctdbase.org/) was used to find integrated chemical-gene, chemical-disease, and gene-disease interactions to generate expanded networks and predict novel associations [18]. We used the database to analyze relationships between CpG sites and diseases. Then, the relationships between CpG sites and diseases were identified.

Statistical analysis

Socio-demographic characteristics were compared between CDMR and VD using median values with ranges. When analyzing the effects of covariates on DNA methylation, beta-values were converted to M-values. The methylation status was compared between CDMR and VD using t-test and multiple linear regression model. In addition, we applied the empirical Bayes smoothing to the standard errors using the R Bioconductor limma pipeline to detect DMPs on M-value. For multiple testing, FDR correction was performed. A significance level of 0.05 is used in the analyses. SAS version 9.4 for Windows (SAS Institute, Inc., Cary, NC) was used to compare subjects’ characteristics while R was used to analyze the locus specific methylation.

Results

Table 1 shows that there were no significant differences in maternal characteristics, gender distribution or birthweight between the VD and CDMR groups except that the VD group had slightly more multipara than the CDMR group.

Table 1 The characteristics of the study group in the DNA methylation analysis by delivery mode

A total of 180 DMPs exhibited differences in DNA methylation by 10% or more between the two groups. After 15 DMPs located on the sex chromosomes were excluded, 165 DMPs still remained (P < 0.05) (Additional file Table S2). Among them, 121 (73.3%) were associated with known genes, and the majority of the DMPs were in gene bodies. Among them, 61.2% (101 of 165) were found to be hypermethylated in the CDMR as compared to VD (Fig. S1 in Supplemental Material). After the conversion of Beta-values to M-values, we used multiple linear regression to validate the present results. The p values and FDR values were shown in Table S4 in the Supplemental Material.

The estimate proportions of cell type in our samples were shown in Table S5 in the Supplemental Material. These cell types contained B cells, CD4 T cells, CD8 T cells, granulocytes, monocytes, natural killer cells (NK), and nucleated red blood cells (nRBC). According to the estimate, the majority of cord blood was granulocyte, followed by CD4 T cells. Multiple linear regression model adjusting for the estimate proportions of cell type was used to avoid potential confounding effect of the cell type (Additional file Table S4).

A total of 20 GO terms were enriched based on the 450 K analysis (Table 2). Many of them have been reported to be associated with the immune system. Furthermore, the DMPs were overlapped with 10 KEGG pathways that are potentially associated with the immune system (Table 3).

Table 2 Significantly enriched gene ontology (GO) terms of differentially methylated genes in VD vs CDMR group based on the 450 K analysis
Table 3 Significantly enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) terms of differentially methylated genes in VD vs CDMR group based on the 450 K analysis

Additionally, we applied the empirical Bayes smoothing t using the R Bioconductor limma pipeline to detect DMPs on M-value. With the threshold of FDR (< 0.01) and a difference in the mean β value |Δβ| ≥ 5% [19], 1208 DMPs were identified after adjusting the covariates. Under the threshold difference in the mean β value |Δβ| > 10%, 4 DMPs were identified and shown in Table S6.

In the targeted bisulfite sequencing assay, the mean coverage of sequencing was 600 base-pair, 497 CpGs around the selected 40 DMPs were identified finally. There were 212 CpG sites around 25 DMPs were different between the CDMR and VD groups. The methylation level of 148 CpG sites (148/212, 69.8%) was higher in the CDMR group than in the VD group. The methylation level of 64 CpG sites (64/212, 30.2%) was lower in the CDMR group than in the VD group. The average methylation level of all CpG sites was calculated as the methylation level of the genes (Additional file Table S3). A total of 22 (55%) genes showed differences between CDMR and VD. Six genes (6/22, 27.3%) had lower methylation status in CDMR group than in VD group, and 16 genes (16/22, 72.7%) had higher methylation level in CDMR group than VD group (Table 4). Among them, five genes were related to the immune system, including phosphatidylinositol- 4,5- bisphosphate 3- kinase catalytic subunit delta (PIK3CD), CXXC finger protein 5 (CXXC5), lymphotoxin β receptor (LTBR), Spi-1 proto-oncogene (SPI1) and serine protease inhibitor B9 (SERPINB9).

Table 4 Differentially methylated genes in VD vs CDMR group

We also used the CTD database to analyze relationships between the above genes and diseases. Table 5 showed that the genes were related to immune system diseases. All the 5 genes were related to inflammation. In addition, PIK3CD and SPI1 were also related to different types of lymphoma.

Table 5 Relationship to immune system diseases related to CpG sites in the validation assay based on the CTD database

Discussion

Our study showed that infants born by CDMR had more locus-specific DNA hypermethylation in genes that were related to the immune system. However, the differences might be largely explained by the cell type proportions.

Evidence has indicated that CD may be associated with the development of various disorders in childhood. Short-term risks of CD include altered immune development, an increased likelihood of allergy, atopy, and asthma, and reduced intestinal gut microbiome diversity [20]. At 6 months of age, children born by CD had a significantly higher BMI [21]. Continuously, CD was associated with offspring obesity until age 20 to 28 years [8]. Besides, CD was weakly associated with otitis media in childhood, especially in the first 7 years of life [22]. In China, children born by CDMR prior to 39 gestational weeks were identified with an increased risk of emotional and behavioral problems at preschool age [23]. Moreover, emergency or planned CD was consistently associated with a modestly increased risk of autism spectrum disorder (ASD) from gestational weeks of 36 to 42 when compared with vaginal delivery in four Nordic countries and Western Australia [24].

Epigenetic changes have been shown to be essential during hematopoiesis by contributing to the formation of stable, heritable gene expression patterns. Aberrant DNA methylation has been associated with several immune deficiencies and autoimmune disorders [25]. A number of studies have added to the growing body of evidence supporting significant associations of epigenetic regulation of gene expression and asthma and allergy phenotypes [26]. DNA hypermethylation has also been identified to be related to more than six types of autoimmune diseases, including rheumatoid arthritis, Sjogren’s syndrome, Graves’ disease, Type 1 diabetes, and Crohn’s disease [27].

A previous study found that DNA-methylation in fetal leucocytes was higher in CDMR births than that in VD [9]. It was also reported that methylation of single genes (e.g., ELA2 and IRF1) was significantly elevated in elective CD than in VD [10]. These results were consistent with our findings. Our study extended and refined evidence in this area. After the GO enrichment, we found 20 GO terms, most of which were related to the immune system. In the present pathway analysis, the DMPs were overlapped with 10 KEGG pathways, including asthma and several autoimmunity diseases.

Moreover, several genes involved in the normal function of the immune system were differentially methylated in VD comparing to CDMR. Based on the validation analysis, the different methylation status was also found in PIK3CD, CXXC5, LTBR, SPI1 and SERPINB9 genes. Of them, three genes were involved in the maturation and normal function of B cells, including PIK3CD, CXXC5 and SPI1. The PIK3CD, one of the phosphoinositide 3-kinase (PI3K) catalytic subunits, has emerged as a critical mediator of multiple B cell functions [28]. PIK3CD was confirmed to participate in the proliferation of B cells and activation in humoral immunity, mediating the survival pathways of T cells and B cells [29]. Based on the previous study, CXXC5 was also involved in the regulation of Toll-like receptors (TLRs) response and antibody production in B cells [30]. Additionally, SPI1 encoded an ETS-domain transcription factor that activated gene expression during myeloid and B-lymphoid cell development [31].

The other 2 genes participated in the process of immune response. LTBR was one of the subunits of lymphotoxin, which was a member of the tumor necrosis factor (TNF) superfamily [32]. The LTBR activation induced the NFκB pathways, which leads to inflammatory gene transcription and nuclear translocation, had a major role in lymphoid organogenesis, immune response and initiation of inflammation [33, 34]. Moreover, SERPINB9 was considered as the main inhibitor for the cytotoxic effector protein granzyme B and caspase-1, which was involved in the immune response [35, 36].

These changes in the above gene expression levels had profound implications for immune function development in the neonates. The immunity abnormity could lead to a series of chronic inflammation diseases and autoimmune diseases, including viral encephalitis, hepatitis B, IgA nephropathy, type 1 diabetes, lymphoblastic leukemia, and cancers [30, 37,38,39,40].

However, the mechanisms on how different delivery mode affects different DNA methylation status remained unclear. A plausible explanation was that the long and extreme fetal stress during labor and his physiological responses might program DNA methylation and gene expression [10]. It was also reported that fetal exposure to glucocorticoids may change the epigenetic state, affecting both DNA methylation and H3K9 acetylation, resulting in significant variations in genome-wide promoter methylation [41]. Consequently, the mode of delivery may affect the banking quality in terms of the number of fetal cells [42].

We also observed that cell type proportions in the umbilical cord blood of CDMR and VD infants explained a large portion of variation in analyses of DNA methylation patterns. When we selected DMPs on M-value, the number of DMPs showing significant DNA methylation associated with delivery mode was reduced from 1208 to 2 single CpG when controlling for blood cell type proportions. After the estimation of cell types in cord blood using FlowSorted.CordBloodCombined.450 k, significant differences were found in CD8T cells, CD4T cells, B cells and Granulocytes between the two groups. More CD8T cells, CD4T cells and B cells were found in CDMR group, while more Granulocytes were found in VD group. It was reported that neonates born by caesarean delivery had lower T and NK lymphocyte counts than neonates born by vaginal delivery. Within the T-lymphocyte population, helper and naive T lymphocyte counts decreased in neonates born by caesarean delivery [43]. Significant associations were also found between basophilic granulocyte counts, plasmacytoid dendritic cells and delivery mode [44].

The present study used cord blood, which were mixed-cell samples. The previous studies using blood samples from both adults and newborns illustrated that a large number of DMPs could be influenced by the cellular composition of the blood sample [45, 46]. Similar results were also found in gene expression study that used adult blood samples to define the cell type specific gene expression patterns [47]. Purification of the specific cell type would be the most reliable approach for eliminating the interindividual variability. However, a large volume of blood and specialized technology would be required to obtain sufficient number of various cell types.

Our study has several strengths. First, only women who had CDMR were selected as the exposed group in our study, to avoid any potential effects of maternal and/or fetal complications. This study is the first to examine the effect of CDMR on DNA methylation status in the umbilical cord blood of the infants. Second, the cord blood DNA methylation status was measured using commercial DNA methylation chip and targeted bisulfite sequencing assay, which provided reliable assessments. Third, the women who had CDMR and VD were matched on the basic characteristics, to avoid the potential confounding effects.

Nevertheless, our study had limitations. The sample size was still small, which may have limited our ability to identify more differentially methylated positions. In addition, there may be many crucial DMPs not included in the commercial chip. Finally, the cell proportions in cord blood were not well considered in the present analysis as our stored samples were no longer suitable for cell typing. Although we have used FlowSorted.CordBloodCombined.450 k to estimate the proportion of cells, further methylation analysis in specific cell types would be more reliable. Further studies focusing on the specific pathway with larger sample size should be performed.

Conclusions

Our study suggested that DNA methylation patterns were altered in different delivery modes, while the differences might be largely explained by the cell type proportions. Further studies investigating the DNA methylation status of the specific genes may be essential to understand the long-lasting effect of the delivery mode on the offspring’s health.

Availability of data and materials

The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

CDMR:

Cesarean deliveries on maternal request

VD:

Vaginal deliveries

DMPs:

Differentially methylated positions

CD:

Cesarean delivery

BMI:

Body mass index

CpG:

Cytosine phosphate-guanine

GO:

Gene ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

HLA:

Human leukocyte antigen

PI3K:

Phosphoinositide 3-kinase

TCRs:

T cell antigen receptors

BCRs:

B cell antigen receptors

TLRs:

Toll-like receptors

LTBR:

Lymphotoxin β receptor

TNF:

Tumor necrosis factor

SERPINB9:

Serine protease inhibitor B9

References

  1. 1.

    Benfield RD, Newton ER, Tanner CJ, Heitkemper MM. Cortisol as a biomarker of stress in term human labor: physiological and methodological issues. Biol Res Nurs. 2014;16(1):64–71. https://doi.org/10.1177/1099800412471580.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

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

    Article  PubMed  Google Scholar 

  3. 3.

    Jun Zhang LY, Meikle S, Zheng J, Sun W, Li Z. Cesarean delivery on maternal request in Southeast China. Obstet Gynecol. 2008;111(5):6–1082. https://doi.org/10.1097/AOG.0b013e31816e349e.

    Article  Google Scholar 

  4. 4.

    Betran AP, Torloni MR, Zhang JJ, Gülmezoglu AM, the WHO Working Group on Caesarean Section, Aleem HA, et al. WHO statement on caesarean section rates. BJOG. 2016;123(5):667–70. https://doi.org/10.1111/1471-0528.13526.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Cuppari C, Manti S, Salpietro A, Alterio T, Arrigo T, Leonardi S, et al. Mode of delivery and risk for development of atopic diseases in children. Allergy Asthma Proc. 2015;36(5):344–51. https://doi.org/10.2500/aap.2015.36.3870.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Wu P, Feldman AS, Rosas-Salazar C, James K, Escobar G, Gebretsadik T, et al. Relative importance and additive effects of maternal and infant risk factors on childhood asthma. PLoS One. 2016;11(3):e0151705. https://doi.org/10.1371/journal.pone.0151705.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Brandão HV, Vieira GO, de Oliveira VT, et al. Increased risk of allergic rhinitis among children delivered by cesarean section: a cross-sectional study nested in a birth cohort. BMC Pediatr. 2016;16(1):1–8. https://doi.org/10.1186/s12887-016-0594-x.

    Article  Google Scholar 

  8. 8.

    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:e162385. https://doi.org/10.1001/jamapediatrics.2016.2385 Association Between Cesarean Birth and Risk of Obesity in Offspring in Childhood, Adolescence, and Early Adulthood.

  9. 9.

    Schlinzig T, Johansson S, Gunnar A, Ekström TJ, Norman M. Epigenetic modulation at birth - altered DNA-methylation in white blood cells after caesarean section. Acta Paediatr. 2009;98(7):1096–9. https://doi.org/10.1111/j.1651-2227.2009.01371.x.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Franz MB, Poterauer M, Elhenicky M, Stary S, Birner P, Vinatzer U, et al. Global and single gene DNA methylation in umbilical cord blood cells after elective caesarean: a pilot study. Eur J Obstet Gynecol Reprod Biol. 2014;179:121–4. https://doi.org/10.1016/j.ejogrb.2014.05.038.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Almgren M, Schlinzig T, Gomez-Cabrero D, et al. Cesarean delivery and hematopoietic stem cell epigenetics in the newborn infant: implications for future health? Am J Obstet Gynecol. 2014;211(5):502 e1–8. https://doi.org/10.1016/j.ajog.2014.05.014.

    Article  PubMed  Google Scholar 

  12. 12.

    Reik W, Dean W. Epigenetic reprogramming: Back to the beginning. Nature. 2002;420(6912):127. https://doi.org/10.1038/420127a.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Dennis C. Epigenetics and disease: altered states. Nature. 2003;421(6924):686–8. https://doi.org/10.1038/421686a.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Wilhelm-Benartzi CS, Koestler DC, Karagas MR, Flanagan JM, Christensen BC, Kelsey KT, et al. Review of processing and analysis methods for DNA methylation array data. Br J Cancer. 2013;109(6):1394–402. https://doi.org/10.1038/bjc.2013.496.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Cheishvili D, Parashar S, Mahmood N, et al. Identification of an Epigenetic Signature of Osteoporosis in Blood DNA of Postmenopausal Women. J Bone Miner Res. 2018;33(11):1980–9. https://doi.org/10.1002/jbmr.3527 published Online First: 2018/06/21.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57. https://doi.org/10.1038/nprot.2008.211 published Online First: 2009/01/10.

    CAS  Article  Google Scholar 

  17. 17.

    Gervin K, Salas LA, Bakulski KM, et al. Systematic evaluation and validation of reference and library selection methods for deconvolution of cord blood DNA methylation data. Clin Epigenetics. 2019;11(1):125. https://doi.org/10.1186/s13148-019-0717-y published Online First: 2019/08/29.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Davis AP, Grondin CJ, Johnson RJ, et al. The Comparative Toxicogenomics Database: update 2017. Nucleic Acids Res. 2017;45(D1):D972–D78. https://doi.org/10.1093/nar/gkw838 published Online First: 2016/09/22.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Husquin LT, Rotival M, Fagny M, et al. Exploring the genetic basis of human population differences in DNA methylation and their causal impact on immune gene regulation. Genome Biol. 2018;19(1):222. https://doi.org/10.1186/s13059-018-1601-3 published Online First: 2018/12/20.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Sandall J, Tribe RM, Avery L, Mola G, Visser GHA, Homer CSE, et al. Short-term and long-term effects of caesarean section on the health of women and children. Lancet. 2018;392(10155):1349–57. https://doi.org/10.1016/s0140-6736(18)31930-5.

    Article  PubMed  Google Scholar 

  21. 21.

    Masukume G, McCarthy FP, Baker PN, et al. Association between caesarean section delivery and obesity in childhood: a longitudinal cohort study in Ireland. BMJ Open. 2019;9(3):e025051. https://doi.org/10.1136/bmjopen-2018-025051 published Online First: 2019/03/18.

    Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Hartley M, Woolcott CG, Langley JM, et al. Birth by Caesarean section and otitis media in childhood: a retrospective cohort study. Sci Rep. 2020;10(1):5219. https://doi.org/10.1038/s41598-020-62229-y published Online First: 2020/04/07.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Huang K, Yan S, Wu X, et al. Elective caesarean section on maternal request prior to 39 gestational weeks and childhood psychopathology: a birth cohort study in China. BMC Psychiatry. 2019;19(1):22. https://doi.org/10.1186/s12888-019-2012-z published Online First: 2019/01/16.

    Article  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Yip BHK, Leonard H, Stock S, et al. Caesarean section and risk of autism across gestational age: a multi-national cohort study of 5 million births. Int J Epidemiol. 2017;46(2):429–39. https://doi.org/10.1093/ije/dyw336 published Online First: 2016/12/27.

    Article  PubMed  Google Scholar 

  25. 25.

    Suarez-Alvarez B, Rodriguez RM, Fraga MF, et al. DNA methylation: a promising landscape for immune system-related diseases. Trends Genet. 2012;28(10):506–14. https://doi.org/10.1016/j.tig.2012.06.005 published Online First: 2012/07/25.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Davidson EJ, Yang IV. Role of epigenetics in the development of childhood asthma. Curr Opin Allergy Clin Immunol. 2018;18(2):132–8. https://doi.org/10.1097/ACI.0000000000000429 published Online First: 2018/02/02.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Ehrlich M. DNA hypermethylation in disease: mechanisms and clinical relevance. Epigenetics. 2019;14(12):1141–63. https://doi.org/10.1080/15592294.2019.1638701 published Online First: 2019/07/10.

    Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Pauls SD, Lafarge ST, Landego I, Zhang T, Marshall AJ. The phosphoinositide 3-kinase signaling pathway in normal and malignant B cells: activation mechanisms, regulation and impact on cellular functions. Front Immunol. 2012;3:224. https://doi.org/10.3389/fimmu.2012.00224.

    Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Delgado P, Cubelos B, Calleja E, Martínez-Martín N., Ciprés A., Mérida I., Bellas C., Bustelo X.R., Alarcón B. Essential function for the GTPase TC21 in homeostatic antigen receptor signaling. Nat Immunol 2009;10(8):880–8. doi: https://doi.org/10.1038/ni.1749, 888.

  30. 30.

    Ma S, Wan X, Deng Z, Shi L, Hao C, Zhou Z, et al. Epigenetic regulator CXXC5 recruits DNA demethylase Tet2 to regulate TLR7/9-elicited IFN response in pDCs. J Exp Med. 2017;214(5):1471–91. https://doi.org/10.1084/jem.20161149.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Zhang JA, Mortazavi A, Williams BA, Wold BJ, Rothenberg EV. Dynamic transformations of genome-wide epigenetic marking and transcriptional control establish T cell identity. Cell. 2012;149(2):467–82. https://doi.org/10.1016/j.cell.2012.01.056.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Rennert PD, Browning JL, Mebius R, Mackay F, Hochman PS. Surface lymphotoxin alpha/beta complex is required for the development of peripheral lymphoid organs. J Exp Med. 1996;184(5):1999–2006. https://doi.org/10.1084/jem.184.5.1999.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    S B, SC L. Functions of NF-kappaB1 and NF-kappaB2 in immune cell biology. Biochem J. 2004;382(2):393–409. https://doi.org/10.1042/BJ20040544.

    Article  Google Scholar 

  34. 34.

    Piao W, Xiong Y, Famulski K, Brinkman CC, Li L, Toney N, et al. Regulation of T cell afferent lymphatic migration by targeting LTbetaR-mediated non-classical NFkappaB signaling. Nat Commun. 2018;9(1):3020. https://doi.org/10.1038/s41467-018-05412-0.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Rizzitelli A, Meuter S, Vega Ramos J, Bird CH, Mintern JD, Mangan MSJ, et al. Serpinb9 (Spi6)-deficient mice are impaired in dendritic cell-mediated antigen cross-presentation. Immunol Cell Biol. 2012;90(9):841–51. https://doi.org/10.1038/icb.2012.29.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Mangan MS, Vega-Ramos J, Joeckel LT, Mitchell AJ, Rizzitelli A, Roediger B, et al. Serpinb9 is a marker of antigen cross-presenting dendritic cells. Mol Immunol. 2017;82:50–6. https://doi.org/10.1016/j.molimm.2016.12.011.

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    Wang X, Lei D, Ding J, Liu S, Tao L, Zhang F, et al. A DNA-methylated sight on autoimmune inflammation network across RA, pSS, and SLE. J Immunol Res. 2018;2018:4390789–13. https://doi.org/10.1155/2018/4390789.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  38. 38.

    Lucas CL, Chandra A, Nejentsev S, Condliffe AM, Okkenhaug K. PI3Kdelta and primary immunodeficiencies. Nat Rev Immunol. 2016;16(11):702–14. https://doi.org/10.1038/nri.2016.93.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Liu Y. A global immune gene expression signature for human cancers. Oncotarget. 2019;10(20):1993-2005. https://doi.org/10.18632/oncotarget.26773 published online first 2019/03/08.

  40. 40.

    Rvd B, Meeldijk J, Jongeneel L, et al. Reduced serpinB9-mediated caspase-1 inhibition can contribute to autoinflammatory disease. Oncotarget. 2016;7(15):19265-71. https://doi.org/10.18632/oncotarget.8086 published online first 2016/03/15.

  41. 41.

    Crudo A, Suderman M, Moisiadis VG, Petropoulos S, Kostaki A, Hallett M, et al. Glucocorticoid programming of the fetal male hippocampal epigenome. Endocrinology. 2013;154(3):1168–80. https://doi.org/10.1210/en.2012-1980.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Yang H, Loutfy MR, Mayerhofer S, Shuen P. Factors affecting banking quality of umbilical cord blood for transplantation. Transfusion. 2011;51(2):284–92. https://doi.org/10.1111/j.1537-2995.2010.02826.x.

    Article  PubMed  Google Scholar 

  43. 43.

    Duijts L, Bakker-Jonges L, Labout J, et al. Perinatal stress influences lymphocyte subset counts in neonates. The generation R study. Pediatr Res. 2008;63(3):292–8. https://doi.org/10.1203/PDR.0b013e318163a29f.

    Article  PubMed  Google Scholar 

  44. 44.

    Lura MP, Gorlanova O, Muller L, et al. Response of cord blood cells to environmental, hereditary and perinatal factors: A prospective birth cohort study. PLoS One. 2018;13(7):e0200236. https://doi.org/10.1371/journal.pone.0200236 published Online First: 2018/07/07.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Lam LL, Emberly E, Fraser HB, et al. Factors underlying variable DNA methylation in a human community cohort. Proc Natl Acad Sci U S A. 2012;109(Suppl 2):17253–60. https://doi.org/10.1073/pnas.1121249109 published Online First: 2012/10/10.

    Article  PubMed  PubMed Central  Google Scholar 

  46. 46.

    Braid SM, Okrah K, Shetty A, et al. DNA Methylation Patterns in Cord Blood of Neonates Across Gestational Age: Association With Cell-Type Proportions. Nurs Res. 2017;66(2):115–22. https://doi.org/10.1097/NNR.0000000000000210 published Online First: 2017/01/27.

    Article  PubMed  PubMed Central  Google Scholar 

  47. 47.

    Palmer C, Diehn M, Alizadeh AA, et al. Cell-type specific gene expression profiles of leukocytes in human peripheral blood. BMC Genomics. 2006;7:115. https://doi.org/10.1186/1471-2164-7-115 published Online First: 2006/05/18.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgments

We thank all the participants of the present study, and the workers in the Ministry of Education-Shanghai Key Laboratory of Children’s Environmental Health.

Funding

The design of this study was supported by the Shanghai Municipal Health and Family Planning Commission (GWIII-26; GWIV-26; 15GWZK0401) and the Xinhua Hospital Early Life Plan. The measurement and collection of data in the present study was supported by grants from National Natural Science Foundation of China (81530086; 81273091; 81803244; 81703235).

Author information

Affiliations

Authors

Contributions

LH, QC and JZ designed the study. QC analyzed the data, and drafted the manuscript. QC, LH, YZ and XW collected the data. QC, YM, YG and YL performed the experiments. All authors provided critical comments and input in the paper. JZ extensively edited the manuscript and takes the overall responsibility for the study. All authors have reviewed the final version of the manuscript and agreed to the submission.

Corresponding author

Correspondence to Jun Zhang.

Ethics declarations

Ethics approval and consent to participate

The ethics approval was obtained from the Ethics Committee of Xinhua Hospital affiliated to the Shanghai Jiao Tong University School of Medicine. Informed consent was obtained from the all participants. This study was conducted in accordance with the principles of the Helsinki Declaration.

Consent for publication

Not applicable.

Competing interests

The authors report no conflict of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1 Table S1

. Information of primers used in targeted bisulfite sequencing assay. Table S2. Differentially methylated loci between vaginal delivery and CDMR group from the Methylation Chip. Table S3. Differentially methylated loci between vaginal delivery and CDMR group from targeted bisulfite sequencing assay. Table S4. The comparison of M-values between VD and CDMR group using multiple linear regression. Table S5. The estimate proportions of cell type in cord blood. Table S6. Differentially methylated CpG sites between VD and CDMR group from the Methylation Chip analyzed by M values, under the threshold difference in the mean β value |Δβ| ≥ 10%. Fig. S1. A volcano plot of lg-transformed P values vs differences in DNA methylation (b-value) between VD and CDMR group.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Chen, Q., Ming, Y., Gan, Y. et al. The impact of cesarean delivery on infant DNA methylation. BMC Pregnancy Childbirth 21, 265 (2021). https://doi.org/10.1186/s12884-021-03748-y

Download citation

Keywords

  • Cesarean deliveries on maternal request
  • Cord blood
  • DNA methylation
  • Perinatal epidemiology
\