- Open Access
Placental surface area mediates the association between FGFR2 methylation in placenta and full-term low birth weight in girls
Clinical Epigeneticsvolume 10, Article number: 39 (2018)
Fibroblast growth factor receptor 2 (FGFR2) gene encodes a protein of the fibroblast growth factor receptor family. FGFR2 gene expression is associated with the regulation of implantation process of placenta which plays a vital role in fetal growth. DNA methylation is widely known as a mechanism of fetal growth. However, it is unclear whether and how DNA methylation of FGFR2 gene in the placenta is associated with full-term low birth weight. This case-control study aims to explore the links between FGFR2 methylation in placenta and full-term low birth weight and to further examine the mediation effect of placental surface area on this association.
We conducted analyses for each of the five valid CpG sites at FGFR2 in 165 mother-baby pairs (86 FT-LBW vs. 79 FT-NBW) and found that per one standard deviation increase in the DNA methylation of CpG 11 at FGFR2 was associated with 1.64-fold higher risk of full-term low birth weight (OR = 1.64, 95% CI = [1.07, 2.52]) and 0.18 standard deviation decrease in placental surface area (β = − 0.18; standard error = 0.08, p = 0.02). The mediation effect of placental surface area on the association between DNA methylation and full-term low birth weight was significant in girls (OR = 1.38, 95% CI = [1.05, 1.80]) but not in boys. The estimated mediation proportion was 48.38%.
Our findings suggested that placental surface area mediated the association between DNA methylation of FGFR2 in placenta and full-term low birth weight in a sex-specific manner. Our study supported the importance of placental epigenetic changes in placental development and fetal growth.
Low birth weight (LBW), defined as the newborn birth weight less than 2500 gram (g) , is one of the most common adverse birth outcomes that affects about 7.2% of pregnancies in China  and as much as 15–20% worldwide . LBW not only increases morbidity and mortality in infants but is also associated with increased susceptibility to chronic non-communicable diseases in adulthood, such as type 2 diabetes, obesity, hypertension, stroke, and cardiovascular diseases [4,5,6,7,8,9,10]. LBW is the result of either preterm birth or intrauterine growth retardation (IUGR), or both . To distinguish the IUGR-driven LBW from prematurity LBW, in this study, we only focused on full-term low birth weight (FT-LBW), the LBW with gestational age from 37 to 42 weeks .
The placenta is a temporary organ that supports fetal growth and development . LBW can be caused by placental development insufficiency  which manifests in different ways including lower placental weight [15,16,17], shorter breadth of the placental surface , and less placental surface area. Placental surface area reflects the invasion of trophoblasts into the maternal decidua and the interface development , which affects mother-fetal substance transmission and fetal growth . However, the molecular mechanism of placental development is still poorly understood.
Fibroblast growth factors (FGFs) are critical in the regulation of implantation process of the placenta, including trophoblast differentiation, hormone production , and migration . Additionally, upregulated FGF2 signaling pathways participate in the regulation of human placental artery endothelial cell proliferation and angiogenesis [23,24,25,26,27]. FGFs carry out these regulation functions by binding to their receptors (FGFRs), and one of these receptors is the fibroblast growth factor receptor 2 (FGFR2), a high-affinity receptor of most FGFs . Animals and human studies have demonstrated that FGFR2 expresses on the membrane of trophoblast [28,29,30]. Furthermore, inhibited FGFR2 expression causes decrease in trophoblast formation and delayed the timing of trophoblast outgrowth [31, 32]. Therefore, it is worth looking more closely into the associations of FGFR2 with placental development and fetal growth.
DNA methylation is a well-documented mechanism of fetal growth . The importance of DNA methylation lies in the fact that DNA methylation participates in regulating gene expression ; also, DNA methylation was proved to be a critical component of fetal programming and long-term onset diseases . LBW has been found to be associated with placental DNA methylation on a general level or of target genes (hypomethylated IGF2 , HUS1B , and hypermethylated HSD11B2, WNT2, and AHRR [38, 39]). Studies have also suggested the vital role of DNA methylation in placental development [40,41,42]; however, the associations between DNA methylation of specific genes and placental development remain largely unknown. Hence, it is valuable to determine the association between DNA methylation of placenta itself and development of the placenta, and fetal growth, also.
We hypothesized that FGFR2 DNA methylation changes might contribute to maldevelopment and poor implantation of the placenta, followed by fetal growth disruption. Here, this study aims to investigate the association of placental FGFR2 methylation with placental development and FT-LBW and to investigate whether placenta mediates the association between placental FGFR2 methylation and FT-LBW.
Study design and subjects
A case-control study design was employed. FT-LBW (< 2500 g at birth) mother-baby pairs were defined as cases, while those with birth weight in the range of 2500 to 4000 g and with a gestational age of 37 to 42 weeks (full-term normal birth weight, FT-NBW) were regarded as controls. During September 2009 and March 2011, 86 FT-LBW and 79 FT-NBW mother-baby pairs were enrolled at the maternity and child health care hospitals of two cities (Foshan and Shenzhen) in China. The cases and controls were individually matched for gestational age (± 7 days), pre-pregnant BMI (± 1.5), parity (the same), and newborn gender (the same). We excluded the subjects if (a) the mothers had heart disease, hepatitis, kidney disease, hypertension (gestational hypertension), diabetes mellitus (gestational diabetes mellitus), hyperthyreosis, anemia, or tuberculosis; (b) the mothers used hazardous drugs to the fetus; (c) the mothers had alcohol abuse; (d) the mothers had multiple gestation; (e) the newborns had hereditary disease or congenital malformation.
We surveyed the pregnant women face-to-face using questionnaires to inquire about their demographic characteristics, alcohol and tobacco use, environmental tobacco smoke (ETS) exposure during pregnancy, medical history, and reproductive history. We collected information of last menstrual, pre-pregnancy body mass index (BMI), parity, placental weight and size, birth weight, birth length, head circumference (HC), and newborn gender from the medical records of mothers and newborns. We estimated gestational age at delivery by subtracting the last menstrual period (LMP) date from the delivery date or by the ultrasound if the ultrasound estimation differed from the LMP estimation by more than 10 days.
Newborn birth weight was measured using SECA baby weight balance in a scale interval of 5 g, and birth length and HC were measured by tape in a scale interval of 1 mm. The placental surface of the mother side was recognized as an ellipse. The maximal diameter of the surface and its perpendicular bisecting diameter were considered as major axis and minor axis. The placental surface area was calculated using the formula for the area of the ellipse: major axis × minor axis × π/4.
Tissue sample collection
Placental tissues were collected immediately after delivery. Approximately 1 cm3 in size samples were cut from the placental maternal side without the surficial membrane, in the middle of the line between the umbilical cord insertion point and the edge of the placenta, avoiding calcification points. The tissue samples were washed in cold physiological saline until it became pale, then cut into small pieces and stored in a − 80 °C freezer.
DNAs from 165 maternal side placental tissue samples were isolated using TIANamp Genomic DNA Kit (TIANGEN Biotech, Beijing, China). DNA quality was assessed by Thermo Scientific™ NanoDrop 2000 (Thermo Scientific™, San Jose) and DNA gel electrophoresis (Major science). The DNA molecular size of all samples were more than 10 Kb, OD260/280 ≥ 1.8, OD260/230 ≥ 1.5, and the nucleic acid concentration of all samples was in 0.2–2.0 μg/μL.
DNA methylation measurement using Sequenom MassARRAY EpiTYPER
We measured the DNA methylation of CpG sites in the 5′UTR region of FGFR2 gene using Sequenom MassARRAY EpiTYPER approach (Sequenom, San Diego, CA, USA). This region located in chr10:123355182-123355644 (University of California, Santa Cruz, GRCh37/hg19 assembly). The amplicons of this region were generated by polymerase chain reaction using the sodium bisulphite-converted DNAs (Zymo Research, Orange, CA) as the templates. The primers for this region were designed using Epidesigner. The forward and backward primers were aggaagagagTATTGGGGTTTAGATTTTAGGAAGG and cagtaatacgactcactatagggagaaggctCAAAATACAAAAAAATTTTCTACCTCT. DNA methylation level was quantified by Honor Tech company (Beijing, China) using Sequenom MassARRAY platform  (San Diego, CA, USA). The methylation calls were performed by the EpiTYPER software v1.0 (Sequenom, San Diego, CA, USA), which generates quantitative results for each CpG site or an aggregate of multiple CpG sites. In total, our target region contains 16 CpG sites, and 15 out of the 16 sites could be measured independently except CpG 13, because its fragment molecular weight is too large and exceeds the detectable range. The molecular weight of CpG 3- and CpG 4-containing fragments overlap, which reduces the ability to resolve CpG 3 and CpG 4 methylation levels independently; thus, the DNA methylation of CpG 3.4 was estimated based on the average of the two CG-containing fragments. For details about this region, see Additional file 1: Tables S1 and S2. DNA methylation level (β) of each CpG site was calculated as the ratio of the intensity of methylated cytosines over the sum of the intensity of methylated and unmethylated cytosines. DNA methylation level ranges from 0 (unmethylated) to 1 (methylated).
To reduce batch effects, we tested half case and half control samples within each batch. To verify the internal replication, one random sample per batch was tested as a repeat. The difference of the methylation level between repeated samples was less than 5% across all CpG sites. To evaluate the efficiency of bisulphite conversion, one external unmethylated region of INS gene was tested in every batch as a conversion control. The bisulphite conversion efficiency of all samples was greater than 98%. In the samples’ quality control, we eliminated samples having a poor readout for each CpG. Then, we also removed the samples having DNA methylation level lower than 5% or higher than 95% for each CpG, due to the poor accuracy resulted from the detection limitation of Sequenom MassARRAY. In the following CpG quality control, we excluded the CpGs with less than 75% valid samples. Hence, our final analysis included CpG 2 of 130 samples (65 cases; 65 controls), CpG 3.4 of 160 samples (82 cases; 78 controls), CpG 5 of 138 samples (73 cases; 65 controls), CpG 7 of 133 samples (69 cases; 64 controls), and CpG 11 of 124 samples (69 cases; 55 controls). The flowchart of quality control is shown in Additional file 1: Figure S1.
Continuous variables corresponding to normal distribution were described using mean and standard deviation. Categorical variables were described with proportions. Student’s t test and chi-squared test were used to compare the distributions of characteristics of mothers, placentas, and newborns.
We fitted a series of logistic regressions to explore the association of FT-LBW with DNA methylation and placental surface area, adjusting for maternal age, education, family monthly income, ETS exposure during pregnancy, gestational age, and newborn’s sex. To test the effect of prenatal ETS exposure on the findings, we conducted a sensitivity analysis in which the ETS did not include in the model as a covariate. Moreover, we performed another sensitivity analysis in which we excluded the mother-child pairs exposed to prenatal alcohol to test if the alcohol use influences the associations between DNA methylation levels of FGFR2 and FT-LBW. Due to lack of the placental cell components’ measurements, we failed to adjust for cell type composition in our primary analysis. However, we generated an epigenome-wide dataset for a subset (n = 26) of the samples in our present study using Infinium Human Methylation 450K BeadChip (Illumina, San Diego, CA), following the standard manufacturer’s protocols. We displayed the process of quality control of the high-throughput data in Additional file 2, supplementary method. Surrogate variable analysis (SVA) approach  was used to identify the unmeasured sources of variation in the high-throughput DNA methylation data, which can confound the association between DNA methylation and FT-LBW. To test if our findings were driven by the unmeasured variations, in the subset samples, we compared the effect sizes of the individual CpG sites involved in our current study in the model adjusted for age, gender, and surrogate variables to the raw model adjusted for age and gender. The surrogated variables might capture the variations contributed by the cell type composition. Hence, we took the variations in the cell type composition into consideration. Additionally, we replicated the analysis and the comparison using the epigenome-wide 450k data. We performed the replication for only CpG 7 (cg25052156) that is the only CpG site measured independently by both EpiTYPER and Methylation 450K approaches.
We explored the association between DNA methylation and placental surface area using several multiple linear regression models, adjusting for all the confounders above. We also performed these models stratified by newborn’s sex, adjusting for the same confounders except newborn’s sex. To make the coefficients comparable across equations in mediation analysis, we used the standard score, which was calculated using z = (x − μ)/σ, of the placental surface area and DNA methylation in all the models. According to the method of Kenny and Sobel test, we considered the mediation being present when (a) the DNA methylation of any CpG sites at FGFR2 gene was correlated with the placental surface area; (b) the DNA methylation of the CpG site(s) in (a) was/were correlated with FT-LBW without adjusting for the placental surface area; (c) the placental surface area was correlated with FT-LBW; and (d) the association between the DNA methylation of CpG site(s) in (a) and FT-LBW became weaker when the model was additionally adjusted for placental surface area. Using the method provided by VanderWeele , we calculated the direct and indirect effect of DNA methylation on FT-LBW; also, we estimated the proportion of the mediation effect of the placental surface area . Mediation analysis was performed using “medflex” package in R 3.4.0, and the other analyses were conducted using SPSS 21.0 (SPSS Inc., Chicago, IL, USA). We applied Bonferroni correction and considered a P value less than 0.025 to be statistically significant.
Characteristics of all qualified subjects are presented in Table 1. On average, the mothers of FT-LBW newborns have lower education levels, lower family monthly income, and higher exposure to environmental tobacco smoke during pregnancy. None of the mothers have smoked ever or during pregnancy. As expected, the FT-LBW newborns had shorter body length and head circumference; also, their placentas were smaller and lighter.
Associations of DNA methylation of FGFR2 in placenta with placental surface area and the risk of FT-LBW
Comparing to FT-NBW, DNA methylation levels of CpG 3.4, CpG 5, CpG 7, and CpG 11 were higher in FT-LBW, but CpG 2 methylation level was lower in FT-LBW (Fig. 1). The correlation coefficients (r) among CpG 2, CpG 3.4, CpG 5, CpG 7, and CpG 11 are shown in Additional file 1: Figure S2. We defined two blocks of dependent CpGs according to r. One block included CpG 2 and CpG 3.4 (r = 0.66), while the other block included CpG 5, CpG 7, and CpG 11 (r = 0.36 − 0.52). The correlation between the two blocks is negligible (r < 0.30) . Hence, we performed Bonferroni correction for two independent tests. Comparing FT-LBW to FT-NBW, only the DNA methylation of CpG 11 at FGFR2 was associated with an increased risk of FT-LBW (OR = 1.48; 95% CI = [1.01, 2.18]). When we further adjusted for the confounders, this association remained statistically significant with an effect size slightly greater (OR = 1.64; 95% CI = [1.07, 2.52]) (Table 2). Namely, per standard deviation increase in DNA methylation of CpG 11 in placenta was associated with 1.64-fold risk of FT-LBW. On average, FT-LBW babies had 4% higher methylation level of CpG 11 in placenta than FT-NBW babies (Table 2). In the stratification by newborn’s sex, this association was present and significant in girls (OR = 1.92; 95% CI = [1.07, 3.43]) but weaker and non-significant in boys (OR = 1.14; 95% CI = [0.52, 2.48]) (Table 2). In other words, per standard deviation increase in DNA methylation of CpG 11 in placenta was associated with 1.92-fold risk of FT-LBW in girls, but this association was not observed in boys. The mean methylation of CpG 11 was 6% higher in FT-LBW baby girls than in FT-NBW. Two sensitivity analyses demonstrated that the associations between CpG methylation of FGFR2 and FT-LBW were robust in the model unadjusted for ETS (Additional file 1: Table S3) and in the subjects who had no prenatal alcohol exposure (Additional file 1: Table S4). We found an association between DNA methylation of CpG sites and placental surface area only at CpG 11 (β = − 0.18; 95% CI = [− 0.34, − 0.03], p = 0.02) (see the correlation scatter plots in Additional file 1: Figure S3). In the further sex-specific analysis, this association was found only in girls (β = − 0.21; 95% CI = [− 0.41, − 0.01]; p = 0.045) but not in boys (β = − 0.04; 95% CI = [− 0.30, 0.22]; p = 0.77). Placental surface area was significantly associated with FT-LBW (OR = 0.43; 95% CI, 0.28–0.68). This association was found significant in both girls (OR = 0.46; 95% CI = [0.26, 0.79]) and boys (OR = 0.38; 95% CI = [0.17, 0.87].
The association between the DNA methylation of CpG 11 and FT-LBW became weaker and non-significant when the model was additionally adjusted for the placental surface area in girls, compared to the model without adjustment of the placental surface area (Fig. 2a). However, in boys, the CpG 11 methylation was not linked to FT-LBW (Fig. 2b). Therefore, the mediation effect was present in girls and was not detected in boys, and the mediation analysis was only conducted in girls. The total effect was decomposed into a direct effect (OR = 1.41, 95% CI = [0.90, 2.22]) and an indirect effect (mediation effect) (OR = 1.38, 95% CI = [1.05, 1.80]) which attributed the effect of DNA methylation on FT-LBW through placental surface area (Fig. 3). The estimated mediation proportion was 48.38%.
In this study, higher methylation levels at of one novel CpG site (CpG 11, Chr10: 123355343) of FGFR2 gene in the mother-side placental tissue were found in full-term low-birth compared to full-term normal-weight newborns. We also found that this association was more apparent in girls rather than in boys. Additionally, higher methylation level of CpG 11 was associated with lower placental surface area which was also confirmed to be associated with higher risk of FT-LBW. Moreover, the placental surface area played a role in mediating the association between CpG 11 DNA methylation and FT-LBW in girls.
Our findings demonstrated the association between increased placental FGFR2 DNA methylation and FT-LBW. The link between DNA methylation of FGFR family and fetal growth has also been shown in previous studies. Haworth’s epigenome-wide study identified that DNA methylation of cg18566515 (not covered in our study) at FGFR2 significantly increased in cord blood of newborns with birth weight centile greater than 85th compared to the medium birth weight centile (40th–60th) group. However, their study did not compare low birth weight centile (< 15th) group to the 40th–60th birth weight group . In another epigenome-wide study, the DNA methylation level of cg15791248 at FGFR1, an important paralog of FGFR2 , in cord blood was found negatively correlated with birth weight centile . An animal study also demonstrated the role of FGFR2 DNA methylation in growth rate . Our study provides additional evidence in placental tissue of the link between FGFR family DNA methylation and fetal growth. However, the previous studies of the association between genome-wide DNA methylation in placenta and low birth weight did not found DNA methylation alternations of FGFR2 [50, 52,53,54,55]. One possible reason for this inconsistency could lie in the coverage of the epigenome-wide experimental approaches. CpG 11 (Chr10: 123355343, UCSC GRCh37/hg19) found in our study was covered by neither Illumina HumanMethylation 27K nor 450K BeadChip; thus, it was impossible to detect CpG 11 in the previous epigenome-wide studies. Another possible reason might be the different definitions of the outcome. The previous studies were performed on birth weight (consisting of the premature and full term) or small for gestational age (SGA). In contrast, this study only focused on the full-term birth weight.
To explore the genetic variation’s influence on the association between DNA methylation of the target region of FGFR2 and FT-LBW, we identified 156 methylation quantitative trait loci (mQTLs) in relation to CpG 3 (cg22633036) and CpG 7 (cg25052156) through performing a look-up in an online catalog of mQTLs (http://www.mqtldb.org) . None of the 156 identified mQTLs was associated with birth weight according to an online genome-wide association study (GWAS) catalog (https://www.ebi.ac.uk/gwas) [57,58,59,60], which suggested that the associations of CpG 3 (cg22633036) and CpG 7 (cg25052156) DNA methylation with low birth weight might not be confounded by shared genetic variations. However, it is worth to mention that the mQTLs in the online catalog were identified in cord blood which may not reflect the specific mQTL status in placental tissue; moreover, we cannot identify the mQTLs related to the other CpG sites uncovered by the Illumine 450k methylation array. We further explored the DNA methylation level distribution of all valid CpGs and did not find a very frequent pattern that is an important characteristic of DNA methylation distribution influenced by genetic variations . Therefore, the associations between DNA methylation modifications and FT-LBW were less likely driven by the genetic variations. More precise evaluation method and the process of the genetic variation’s effect should be developed and performed in the future study.
To explore the unmeasured variation’s influence on our findings, we estimated five significant surrogate variables in the subset samples of our present study using epigenome-wide methylation data and SVA approach and compared the effect size of the individual CpG sites in the model adjusted for age, gender, and surrogate variables to the raw model only adjusted for age and gender (Additional file 1: Table S5). We found that the effect sizes of CpG 2 and CpG 5 on FT-LBW in the adjusted models (adjusted for age, gender, and surrogate variables) changed as large as 151 and 29%, respectively, compared to those in model 0 (unadjusted for surrogate variables); however, the effect sizes of CpG 3.4, CpG 7, and CpG 11 on FT-LBW in the adjusted models changed less than 10% (2%, 10%, 6%), compared to those in model 0. It suggested that the associations of CpG 2 and CpG 5 with FT-LBW might be strongly driven by placental unmeasured variations while indicated a much weaker effect of the unmeasured variations on the associations of CpG 3.4, CpG 7, and CpG 11 with FT-LBW. This supplementary analysis made us more confident of our major positive finding (CpG 11). In the replication, we found that the effect size of cg25052156 (450k data) was much similar to the effect size of CpG 7 (EpiTYPER data) in both raw model 0 and adjusted model 3 and the directions of the effects were the same (Additional file 1: Table S5), which indicated consistency between two measurements. The effect size of cg25052156 on FT-LBW in model 3 increased 13% compared to that in model 0, which was also similar to the corresponding changes of CpG 7. One thing to note is that the evaluation of the unmeasured variation effect is based on very small samples that have different characteristic distributions from the entire population (Additional file 1: Table S6), which might not adequately reflect the real situation of our original study.
This study, for the first time, elucidated the association between the DNA methylation of FGFR2 in placenta and placental characteristics. Existing biological evidence can help explain this association. The FGFR2 gene expresses in placenta and upregulates trophoblast formation [31, 32]. If gene expression of FGFR2 is inhibited, the embryo implantation and maternal-fetal interface formation will be interrupted , which can be reflected by reduced placental surface area . Furthermore, FGFR2 affects placental angiogenesis and artery endothelial cell proliferation by inducing the FGF2 signal pathway [23,24,25,26,27]. Reduced trophoblast formation and peripheral villous vascularization lead to further intrauterine growth restriction . It is consistent with these biological mechanisms that we observed the mediation effect of placental surface area on the association between FGFR2 methylation and FT-LBW in girls. Although we did not measure the association between FGFR2 methylation and gene expression due to lack of RNA samples, another study of thyroid cancer documented an inverse correlation between the DNA methylation of 5′UTR region at FGFR2 and gene expression . The region we measured in our study is also in the 5′UTR region, and this region locates in a predicted gene expression regulatory element (enhancer) region with a strong signal in Layered H3K27Ac . This evidence highlights the probability of FGFR2 DNA methylation as a regulation mechanism involved in FGFR2 expression and underscores the importance of higher methylation of CpG 11 at FGFR2 observed in our study relative to placental development and fetal growth.
It is well known that the placenta plays a key role in fetal growth. As a proxy for placental development, the placental surface area was found associated with fetal growth . In our study, per one standard deviation increase in the placental surface area was associated with 0.45-fold lower risk of FT-LBW (OR = 0.45 for girls, OR = 0.37 for boys). In other studies, a 1-cm increase in placental axis was associated with 20.9–28.1 g higher birth weight . Larger uteroplacental interface surface area in the first trimester associated with higher birth weight centile . The surface area of the placenta at term was also highly and positively correlated with birth weight . Our findings were in line with these studies.
The associations above between FGFR2 DNA methylation in placenta, placental surface area, and FT-LBW suggested a mediation role of the placental surface area. In our study, we found that the association between FGFR2 DNA methylation and FT-LBW became non-significant after being additionally adjusted for placental surface area, which proved the mediation effect. Placental surface area mediated 48.38% association between FGFR2 DNA methylation and FT-LBW in girls in our study.
We found the associations of FGFR2 methylation with the placental surface area and FT-LBW only in girls but not in boys. After stratifying for sex, we observed greater effect size of FGFR2 methylation on FT-LBW in girls (OR = 1.92) compared to all (OR = 1.64). The sex-specific DNA methylation difference was also observed in other studies of birth weight. For example, in Kippler’s epigenome-wide study , TSH7DA (cg07846874) was found inversely associated with birth weight in girls but not in boys. On the contrary, cg19119945 was only observed in boys. The sex-specific DNA methylation changes with low birth weight were also detected in animal models [70, 71]. However, to our knowledge, none of the studies has characterized the sex-specific DNA methylation of FGFR2 to compare with our findings. Our findings suggested that DNA methylation of CpG 11 at FGFR2 has a sex-specific effect on placenta and fetal growth, but this might have been a question of low statistical power due to the limited sample size in boys.
Our study has several strengths. The outcome of this study focused specifically on full-term low birth weight which eliminated the effect of prematurity on birth weight and directly reflects the effect of intrauterine growth retardation. Secondly, we matched cases and controls for several characteristics which may influence both DNA methylation and birth weight. Thirdly, the mediation analysis method used in this study is specific for the binary outcome and continuous mediator, which allows analyzing the effect estimated with linear regression and logistic regression on the same scale.
Several limitations of the study should be addressed. Firstly, FGFR2 gene expression regulates placental development mostly in the first trimester ; however, we collected placenta samples at delivery. Thus, we were unable to observe the pattern of placental FGFR2 DNA methylation in the early and sensitive stage during pregnancy. Secondly, due to lack of the measures of placental cell type compositions and appropriate computational estimation methods for individual target DNA methylation, we failed to adjust the effects of cell type compositions on the results in this study. Thirdly, although the placental tissue samples were washed in physiological saline water to remove the blood tissue, the maternal DNA contamination effect may not be eliminated completely. Additionally, the placental surface area is not the exclusive proxy of uteroplacental interface development, but the density of vessel and villi are also important metrics. More placental morphology measurements should be performed in future studies. Furthermore, due to lack of RNA data, we were unable to estimate the association between DNA methylation of CpG 11 at FGFR2 and FGFR2 gene expression. The function examination of FGFR2 methylation modification should be conducted in future studies.
Our findings suggested that placental surface area mediates the association between DNA methylation of FGFR2 in placenta and FT-LBW in a sex-specific manner. Our study supported the importance of placental epigenetic changes in placental development and fetal growth. The findings extend our understanding of the epigenetic mechanism of intrauterine growth.
WHO. P07 disorders related to short gestation and low birth weight in icd-10. 2016. http://www.icd10data.com/ICD10CM/Codes/P00-P96/P05-P08/P07-/P07.10. Accessed 08/31/2017.
Chen Y, Wu L, Zhang W, Zou L, Li G, Fan L. Delivery modes and pregnancy outcomes of low birth weight infants in China. J Perinatol. 2016;36(1):41–6. https://doi.org/10.1038/jp.2015.137.
WHO. Global nutrition targets 2025: low birth weight policy brief. 2014. http://www.who.int/nutrition/publications/globaltargets2025_policybrief_lbw/en/. Accessed 08/31/2017.
McNamara BJ, Gubhaju L, Chamberlain C, Stanley F, Eades SJ. Early life influences on cardio-metabolic disease risk in aboriginal populations—what is the evidence? A systematic review of longitudinal and case-control studies. Int J Epidemiol. 2012;41(6):1661–82. https://doi.org/10.1093/ije/dys190.
Pfab T, Slowinski T, Godes M, Halle H, Priem F, Hocher B. Low birth weight, a risk factor for cardiovascular diseases in later life, is already associated with elevated fetal glycosylated hemoglobin at birth. Circulation. 2006;114(16):1687–92. https://doi.org/10.1161/CIRCULATIONAHA.106.625848.
Zarrati M, Shidfar F, Razmpoosh E, Nezhad FN, Keivani H, Hemami MR, et al. Does low birth weight predict hypertension and obesity in schoolchildren? Ann Nutr Metab. 2013;63(1–2):69–76. https://doi.org/10.1159/000351869.
Silverwood RJ, Pierce M, Hardy R, Sattar N, Whincup P, Ferro C, et al. Low birth weight, later renal function, and the roles of adulthood blood pressure, diabetes, and obesity in a British birth cohort. Kidney Int. 2013;84(6):1262–70. https://doi.org/10.1038/ki.2013.223.
Lund LK, Vik T, Lydersen S, Lohaugen GC, Skranes J, Brubakk AM, et al. Mental health, quality of life and social relations in young adults born with low birth weight. Health Qual Life Outcomes. 2012;10:146. https://doi.org/10.1186/1477-7525-10-146.
Walter EC, Koepsell TD, Chien JW. Low birth weight and respiratory hospitalizations in adolescence. Pediatr Pulmonol. 2010; https://doi.org/10.1002/ppul.21398.
Main KM, Jensen RB, Asklund C, Hoi-Hansen CE, Skakkebaek NE. Low birth weight and male reproductive function. Horm Res. 2006;65(Suppl 3):116–22. https://doi.org/10.1159/000091516.
Wardlaw TM, World Health Organization, UNICEF, editors. Low birthweight: country, regional and global estimates. Geneva : New York: WHO ; UNICEF; 2004.
Fedrick J, Adelstein P. Factors associated with low birth weight of infants delivered at term. Br J Obstet Gynaecol. 1978;85(1):1–7.
Novakovic B, Gordon L, Robinson WP, Desoye G, Saffery R. Glucose as a fetal nutrient: dynamic regulation of several glucose transporter genes by DNA methylation in the human placenta across gestation. J Nutr Biochem. 2013;24(1):282–8. https://doi.org/10.1016/j.jnutbio.2012.06.006.
Valero De Bernab J, Soriano T, Albaladejo R, Juarranz M, Calle MAE, Mart Nez D, et al. Risk factors for low birth weight: a review. Eur J Obstet Gyn R B. 2004;116(1):3–15. https://doi.org/10.1016/j.ejogrb.2004.03.007.
Niu Z, Xie C, Wen X, Tian F, Ding P, He Y, et al. Placenta mediates the association between maternal second-hand smoke exposure during pregnancy and small for gestational age. Placenta. 2015;36(8):876–80. https://doi.org/10.1016/j.placenta.2015.05.005.
Roland MC, Friis CM, Voldner N, Godang K, Bollerslev J, Haugen G, et al. Fetal growth versus birthweight: the role of placenta versus other determinants. PLoS One. 2012;7(6):e39324. https://doi.org/10.1371/journal.pone.0039324.
Brett KE, Ferraro ZM, Yockell-Lelievre J, Gruslin A, Adamo KB. Maternal-fetal nutrient transport in pregnancy pathologies: the role of the placenta. Int J Mol Sci. 2014;15(9):16153–85. https://doi.org/10.3390/ijms150916153.
Alwasel SH, Abotalib Z, Aljarallah JS, Osmond C, Al OS, Harrath A, et al. The breadth of the placental surface but not the length is associated with body size at birth. Placenta. 2012;33(8):619–22. https://doi.org/10.1016/j.placenta.2012.04.015.
Roberts DJ, Post MD. The placenta in pre-eclampsia and intrauterine growth restriction. J Clin Pathol. 2008;61(12):1254–60. https://doi.org/10.1136/jcp.2008.055236.
Capellini I. The evolutionary significance of placental interdigitation in mammalian reproduction: contributions from comparative studies. Placenta. 2012;33(10):763–8. https://doi.org/10.1016/j.placenta.2012.07.004.
Powers CJ, McLeskey SW, Wellstein A. Fibroblast growth factors, their receptors and signaling. Endocr Relat Cancer. 2000;7(3):165–97.
Yang QE, Giassetti MI, Ealy AD. Fibroblast growth factors activate mitogen-activated protein kinase pathways to promote migration in ovine trophoblast cells. Reproduction. 2011;141(5):707–14. https://doi.org/10.1530/REP-10-0541.
Feng L, Liao WX, Luo Q, Zhang HH, Wang W, Zheng J, et al. Caveolin-1 orchestrates fibroblast growth factor 2 signaling control of angiogenesis in placental artery endothelial cell caveolae. J Cell Physiol. 2012;227(6):2480–91. https://doi.org/10.1002/jcp.22984.
Mata-Greenwood E, Liao WX, Zheng J, Chen DB. Differential activation of multiple signalling pathways dictates eNOS upregulation by FGF2 but not VEGF in placental artery endothelial cells. Placenta. 2008;29(8):708–17. https://doi.org/10.1016/j.placenta.2008.05.005.
Wang K, Song Y, Chen DB, Zheng J. Protein phosphatase 3 differentially modulates vascular endothelial growth factor- and fibroblast growth factor 2-stimulated cell proliferation and signaling in ovine fetoplacental artery endothelial cells. Biol Reprod. 2008;79(4):704–10. https://doi.org/10.1095/biolreprod.108.068957.
Zheng J, Wen Y, Song Y, Wang K, Chen DB, Magness RR. Activation of multiple signaling pathways is critical for fibroblast growth factor 2- and vascular endothelial growth factor-stimulated ovine fetoplacental endothelial cell proliferation. Biol Reprod. 2008;78(1):143–50. https://doi.org/10.1095/biolreprod.107.064477.
Javerzat S, Franco M, Herbert J, Platonova N, Peille AL, Pantesco V, et al. Correlating global gene regulation to angiogenesis in the developing chick extra-embryonic vascular system. PLoS One. 2009;4(11):e7856. https://doi.org/10.1371/journal.pone.0007856.
Kunath T, Yamanaka Y, Detmar J, MacPhee D, Caniggia I, Rossant J, et al. Developmental differences in the expression of FGF receptors between human and mouse embryos. Placenta. 2014;35(12):1079–88. https://doi.org/10.1016/j.placenta.2014.09.008.
Pfarrer C, Weise S, Berisha B, Schams D, Leiser R, Hoffmann B, et al. Fibroblast growth factor (fgf)-1, fgf2, fgf7 and fgf receptors are uniformly expressed in trophoblast giant cells during restricted trophoblast invasion in cows. Placenta. 2006;27(6–7):758–70. https://doi.org/10.1016/j.placenta.2005.06.007.
Baczyk D, Dunk C, Huppertz B, Maxwell C, Reister F, Giannoulias D, et al. Bi-potential behaviour of cytotrophoblasts in first trimester chorionic villi. Placenta. 2006;27(4–5):367–74. https://doi.org/10.1016/j.placenta.2005.03.006.
Ozawa M, Yang QE, Ealy AD. The expression of fibroblast growth factor receptors during early bovine conceptus development and pharmacological analysis of their actions on trophoblast growth in vitro. Reproduction. 2013;145(2):191–201. https://doi.org/10.1530/REP-12-0220.
Arman E, Haffner-Krausz R, Chen Y, Heath JK, Lonai P. Targeted disruption of fibroblast growth factor (FGF) receptor 2 suggests a role for FGF signaling in pregastrulation mammalian development. Proc Natl Acad Sci U S A. 1998;95(9):5082–7.
Januar V, Desoye G, Novakovic B, Cvitic S, Saffery R. Epigenetic regulation of human placental function and pregnancy outcome: considerations for causal inference. Am J Obstet Gynecol. 2015;213(4 Suppl):S182–96. https://doi.org/10.1016/j.ajog.2015.07.011.
Moore LD, Le T, Fan G. DNA methylation and its basic function. Neuropsychopharmacology. 2013;38(1):23–38. https://doi.org/10.1038/npp.2012.112.
Waterland RA, Michels KB. Epigenetic epidemiology of the developmental origins hypothesis. Annu Rev Nutr. 2007;27:363–88. https://doi.org/10.1146/annurev.nutr.27.061406.093705.
Bouwland-Both MI, van Mil NH, Stolk L, Eilers PH, Verbiest MM, Heijmans BT, et al. DNA methylation of IGF2DMR and H19 is associated with fetal and infant growth: the generation R study. PLoS One. 2013;8(12):e81731. https://doi.org/10.1371/journal.pone.0081731.
Rumbajan JM, Yamaguchi Y, Nakabayashi K, Higashimoto K, Yatsuki H, Nishioka K, et al. The HUS1B promoter is hypomethylated in the placentas of low-birth-weight infants. Gene. 2016;583(2):141–6. https://doi.org/10.1016/j.gene.2016.02.025.
Xiao X, Zhao Y, Jin R, Chen J, Wang X, Baccarelli A, et al. Fetal growth restriction and methylation of growth-related genes in the placenta. Epigenomics-Uk. 2016;8(1):33–42. https://doi.org/10.2217/epi.15.101.
Tian F, Hivert M, Wen X, Xie C, Niu Z, Fan L, et al. Tissue differences in DNA methylation changes at AHRR in full term low birth weight in maternal blood, placenta and cord blood in Chinese. Placenta. 2017;52:49–57. https://doi.org/10.1016/j.placenta.2017.02.017.
Bianco-Miotto T, Mayne BT, Buckberry S, Breen J, Rodriguez LC, Roberts CT. Recent progress towards understanding the role of DNA methylation in human placental development. Reproduction. 2016;152(1):R23–30. https://doi.org/10.1530/REP-16-0014.
Kohan-Ghadr HR, Kadam L, Jain C, Armant DR, Drewlo S. Potential role of epigenetic mechanisms in regulation of trophoblast differentiation, migration, and invasion in the human placenta. Cell Adhes Migr. 2016;10(1–2):126–35. https://doi.org/10.1080/19336918.2015.1098800.
Bashir ST, Gastal MO, Tazawa SP, Tarso SG, Hales DB, Cuervo-Arango J, et al. The mare as a model for luteinized unruptured follicle syndrome: intrafollicular endocrine milieu. Reproduction. 2016;151(3):271–83. https://doi.org/10.1530/REP-15-0457.
Ehrich M, Nelson MR, Stanssens P, Zabeau M, Liloglou T, Xinarianos G, et al. Quantitative high-throughput analysis of DNA methylation patterns by base-specific cleavage and mass spectrometry. Proc Natl Acad Sci U S A. 2005;102(44):15785–90. https://doi.org/10.1073/pnas.0507816102.
Leek JT, Johnson WE, Parker HS, Jaffe AE, Storey JD. The SVA package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 2012;28(6):882–3. https://doi.org/10.1093/bioinformatics/bts034.
Vanderweele TJ, Vansteelandt S. Odds ratios for mediation analysis for a dichotomous outcome. Am J Epidemiol. 2010;172(12):1339–48. https://doi.org/10.1093/aje/kwq332.
Wang W, Zhang B. Assessing natural direct and indirect effects for a continuous exposure and a dichotomous outcome. J Stat Theory Pract. 2016;10(3):574–87. https://doi.org/10.1080/15598608.2016.1203843.
Mukaka MM. Statistics corner: a guide to appropriate use of correlation coefficient in medical research. Malawi Med J. 2012;24(3):69–71.
Haworth KE, Farrell WE, Emes RD, Ismail KM, Carroll WD, Hubball E, et al. Methylation of the FGFR2 gene is associated with high birth weight centile in humans. Epigenomics-Uk. 2014;6(5):477–91. https://doi.org/10.2217/epi.14.40.
Katoh M. FGFR2 abnormalities underlie a spectrum of bone, skin, and cancer pathologies. J Invest Dermatol. 2009;129(8):1861–7. https://doi.org/10.1038/jid.2009.97.
Fryer AA, Emes RD, Ismail KM, Haworth KE, Mein C, Carroll WD, et al. Quantitative, high-resolution epigenetic profiling of CpG loci identifies associations with cord blood plasma homocysteine and birth weight in humans. Epigenetics-Us. 2011;6(1):86–94. https://doi.org/10.4161/epi.6.1.13392.
Hu Y, Xu H, Li Z, Zheng X, Jia X, Nie Q, et al. Comparison of the genome-wide DNA methylation profiles between fast-growing and slow-growing broilers. PLoS One. 2013;8(2):e56411. https://doi.org/10.1371/journal.pone.0056411.
Hillman SL, Finer S, Smart MC, Mathews C, Lowe R, Rakyan VK, et al. Novel DNA methylation profiles associated with key gene regulation and transcription pathways in blood and placenta of growth-restricted neonates. Epigenetics-Us. 2015;10(1):50–61. https://doi.org/10.4161/15592294.2014.989741.
Gordon L, Joo JE, Powell JE, Ollikainen M, Novakovic B, Li X, et al. Neonatal DNA methylation profile in human twins is specified by a complex interplay between intrauterine environmental and genetic factors, subject to tissue-specific influence. Genome Res. 2012;22(8):1395–406. https://doi.org/10.1101/gr.136598.111.
Turan N, Ghalwash MF, Katari S, Coutifaris C, Obradovic Z, Sapienza C. DNA methylation differences at growth related genes correlate with birth weight: a molecular signature linked to developmental origins of adult disease? BMC Med Genet. 2012;5:10. https://doi.org/10.1186/1755-8794-5-10.
Wilhelm-Benartzi CS, Houseman EA, Maccani MA, Poage GM, Koestler DC, Langevin SM, et al. In utero exposures, infant growth, and DNA methylation of repetitive elements and developmentally related genes in human placenta. Environ Health Persp. 2012;120(2):296.
Gaunt TR, Shihab HA, Hemani G, Min JL, Woodward G, Lyttleton O, et al. Systematic identification of genetic influences on methylation across the human life course. Genome Biol. 2016;17(1) https://doi.org/10.1186/s13059-016-0926-z.
Horikoshi M, Beaumont RN, Day FR, Warrington NM, Kooijman MN, Fernandez-Tajes J, et al. Genome-wide associations for birth weight and correlations with adult disease. Nature. 2016;538(7624):248–52. https://doi.org/10.1038/nature19806.
Horikoshi M, Yaghootkar H, Mook-Kanamori DO, Sovio U, Taal HR, Hennig BJ, et al. New loci associated with birth weight identify genetic links between intrauterine growth and adult height and metabolism. Nat Genet. 2013;45(1):76–82. https://doi.org/10.1038/ng.2477.
Urbanek M, Hayes MG, Armstrong LL, Morrison J, Lowe LP, Badon SE, et al. The chromosome 3q25 genomic region is associated with measures of adiposity in newborns in a multi-ethnic genome-wide association study. Hum Mol Genet. 2013;22(17):3583–96. https://doi.org/10.1093/hmg/ddt168.
Metrustry SJ, Edwards MH, Medland SE, Holloway JW, Montgomery GW, Martin NG, et al. Variants close to ntrk2 gene are associated with birth weight in female twins. Twin Res Hum Genet. 2014;17(4):254–61. https://doi.org/10.1017/thg.2014.34.
Daca-Roszak P, Pfeifer A, Żebracka-Gala J, Rusinek D, Szybińska A, Jarząb B, et al. Impact of SNPs on methylation readouts by Illumina Infinium HumanMethylation450 BeadChip Array: implications for comparative population studies. BMC Genomics. 2015;16(1) https://doi.org/10.1186/s12864-015-2202-0.
Red-Horse K, Zhou Y, Genbacev O, Prakobphol A, Foulk R, McMaster M, et al. Trophoblast differentiation during embryo implantation and formation of the maternal-fetal interface. J Clin Invest. 2004;114(6):744–54. https://doi.org/10.1172/JCI22991.
Chen C, Bajoria R, Aplin JD. Decreased vascularization and cell proliferation in placentas of intrauterine growth—restricted fetuses with abnormal umbilical artery flow velocity waveforms. Am J Obstet Gynecol. 2002;187(3):764–9. https://doi.org/10.1067/mob.2002.125243.
Kondo T, Zheng L, Liu W, Kurebayashi J, Asa SL, Ezzat S. Epigenetically controlled fibroblast growth factor receptor 2 signaling imposes on the RAS/BRAF/mitogen-activated protein kinase pathway to modulate thyroid cancer progression. Cancer Res. 2007;67(11):5461–70. https://doi.org/10.1158/0008-5472.CAN-06-4477.
UCSC: H3k27ac mark (often found near active regulatory elements) on 7 cell lines from encode. 2017. https://genome.ucsc.edu/cgi-bin/hgTrackUi?hgsid=509764167_hBFLTAE3Vuhnfv0O6RMuAiUE2hXp&c=chr5&g=wgEncodeRegMarkH3k27ac. Accessed 5 Sept 2017.
Salafia CM, Zhang J, Charles AK, Bresnahan M, Shrout P, Sun W, et al. Placental characteristics and birthweight. Paediatr Perinat Epidemiol. 2008;22(3):229–39. https://doi.org/10.1111/j.1365-3016.2008.00935.x.
Suri S, Muttukrishna S, Jauniaux E. 2D-ultrasound and endocrinologic evaluation of placentation in early pregnancy and its relationship to fetal birthweight in normal pregnancies and pre-eclampsia. Placenta. 2013;34(9):745–50. https://doi.org/10.1016/j.placenta.2013.05.003.
Sivarao S, Vidyadaran MK, Jammal ABE, Zainab S, Goh YM, Ramesh KN. Weight, volume and surface area of placenta of normal pregnant women and their relation to maternal and neonatal parameters in Malay, Chinese and Indian ethnic groups. Placenta. 2002;23(8–9):691–6. https://doi.org/10.1053/plac.2002.0817.
Kippler M, Engstrom K, Mlakar SJ, Bottai M, Ahmed S, Hossain MB, et al. Sex-specific effects of early life cadmium exposure on DNA methylation and implications for birth weight. Epigenetics-Us. 2013;8(5):494–503. https://doi.org/10.4161/epi.24401.
Ke X, McKnight RA, Caprau D, O'Grady S, Fu Q, Yu X, et al. Intrauterine growth restriction affects hippocampal dual specificity phosphatase 5 gene expression and epigenetic characteristics. Physiol Genomics. 2011;43(20):1160–9. https://doi.org/10.1152/physiolgenomics.00242.2010.
Baserga M, Kaur R, Hale MA, Bares A, Yu X, Callaway CW, et al. Fetal growth restriction alters transcription factor binding and epigenetic mechanisms of renal 11beta-hydroxysteroid dehydrogenase type 2 in a sex-specific manner. Am J Physiol Regul Integr Comp Physiol. 2010;299(1):R334–42. https://doi.org/10.1152/ajpregu.00122.2010.
We acknowledge the support of Foshan and Shenzhen Maternity and Child Health Care Hospitals. This work was also supported by the scholarship from Chinese Scholarship Council (CSC).
The study was funded by the National Natural Science Foundation of China (Grant Numbers 30872164 and 81172758) and Guangdong Science and Technology Plan (2016A020218014).
Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
The Institutional Review Board of the public health school of the Sun Yat-Sen University approved the study, and every single pregnant woman was recruited with informed consent.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Quality control flowchart. Figure S2. The correlation of DNA methylation across all the valid CpG sites in valid samples. The correlation coefficients with the gray background are not statistically significant (P > 0.05). Figure S3. The scatter plots of placental surface area versus DNA methylation of CpG sites at FGFR2. The coefficient and P value were given by regression of placental surface area on DNA methylation using multiple linear regression model adjusted for maternal age, education, family monthly income, ETS exposure during pregnancy, gestational age and sex. (DOCX 294 kb)
Includes the supplementary method. (DOCX 25 kb)