rs10732516 polymorphism at the IGF2/H19 locus associates with genotype-specific effects on placental DNA methylation and birth weight of newborns conceived by assisted reproductive technology
Clinical Epigenetics volume 10, Article number: 80 (2018)
Assisted reproductive technology (ART) has been associated with low birth weight of fresh embryo transfer (FRESH) derived and increased birth weight of frozen embryo transfer (FET)-derived newborns. Owing to that, we focused on imprinted insulin-like growth factor 2 (IGF2)/H19 locus known to be important for normal growth. This locus is regulated by H19 imprinting control region (ICR) with seven binding sites for the methylation-sensitive zinc finger regulatory protein (CTCF). A polymorphism rs10732516 G/A in the sixth binding site for CTCF, associates with a genotype-specific trend to the DNA methylation. Due to this association, 62 couples with singleton pregnancies derived from FRESH (44 IVF/18 ICSI), 24 couples from FET (15 IVF/9 ICSI), and 157 couples with spontaneously conceived pregnancies as controls were recruited in Finland and Estonia for genotype-specific examination. DNA methylation levels at the H19 ICR, H19 DMR, and long interspersed nuclear elements in placental tissue were explored by MassARRAY EpiTYPER (n = 122). Allele-specific changes in the methylation level of H19 ICR in placental tissue (n = 26) and white blood cells (WBC, n = 8) were examined by bisulfite sequencing. Newborns’ (n = 243) anthropometrics was analyzed by using international growth standards.
A consistent trend of genotype-specific decreased methylation level was observed in paternal allele of rs10732516 paternal A/maternal G genotype, but not in paternal G/maternal A genotype, at H19 ICR in ART placentas. This hypomethylation was not detected in WBCs. Also genotype-specific differences in FRESH-derived newborns’ birth weight and head circumference were observed (P = 0.04, P = 0.004, respectively): FRESH-derived newborns with G/G genotype were heavier (P = 0.04) and had larger head circumference (P = 0.002) compared to newborns with A/A genotype. Also, the placental weight and birth weight of controls, FRESH- and FET-derived newborns differed significantly in rs10732516 A/A genotype (P = 0.024, P = 0.006, respectively): the placentas and newborns of FET-derived pregnancies were heavier compared to FRESH-derived pregnancies (P = 0.02, P = 0.004, respectively).
The observed DNA methylation changes together with the phenotypic findings suggest that rs10732516 polymorphism associates with the effects of ART in a parent-of-origin manner. Therefore, this polymorphism should be considered when the effects of environmental factors on embryonic development are studied.
Although the results of assisted reproductive technology (ART) in Western countries are impressive and the children born are generally healthy, they have been associated with increased risk of adverse perinatal outcome . Especially, an increased risk of low birth weight and preterm birth have been observed in in vitro fertilized (IVF) singleton pregnancies compared to natural conception [2,3,4]. Furthermore, there has been a suspicion of a higher frequency of imprinting disorders, such as Angelman, Beckwith-Wiedemann, and Silver-Russell syndromes [5, 6]. On the other hand, increased birth weight has been associated with newborns derived from frozen embryo transfer (FET) [7,8,9]. The reason for the differences in perinatal outcome is unclear, but it could be explained by parental characteristics, subfertility, or gonadotrophin stimulation of the ovaries [10,11,12]. Furthermore, some adverse effect of laboratory procedures involving use of culture media, prolonged culturing of the embryos and freezing/thawing methods has not been excluded .
The procedures of IVF and intracytoplasmic sperm injection (ICSI) are performed in the beginning of embryonic development, which is a period of epigenetic reprogramming. During this dynamic period of cell divisions, the epigenetic marks are erased and then established again . The adequate methylation profiles are needed for normal embryonic development, and indeed, it has been recently shown that many developmentally important transcription factors display preference for sequences containing DNA methylation . Altered levels of DNA methylation have been observed not only in different human tissues derived from IVF/ICSI pregnancies [15,16,17] but also in mouse [18, 19], suggesting that IVF protocol, even without infertility can alter the epigenome. Although the results are inconsistent, theoretically IVF could affect the epigenetic reprogramming of early embryo and consequently influence the perinatal outcome.
Owing to the low birth weight associated with IVF and increased birth weight of newborns derived from FET, we focused on the imprinted insulin-like growth factor 2 (IGF2)/H19 locus on chromosome 11p15.5. These two genes are expressed in parent-of-origin manner; IGF2, a major driver of growth, is expressed from paternal allele  and non-coding, negative growth controller H19 from maternal allele . Allele-specific gene expression is needed for normal placental and embryonic growth. The locus is regulated by allele-specific DNA methylation at the H19 imprinting control region (H19 ICR) locating between the genes, as well as H19 promoter region (H19 DMR) and three differentially methylated regions (DMR0, DMR1, and DMR2) at the IGF2 (Fig. 1). H19 ICR contains seven binding sites for a methylation-sensitive, zinc-finger protein CCCTC-binding factor (CTCF). These binding factors organize chromatin contacts and have a critical role in the establishment and maintenance of imprinting . According to mouse studies, unmethylated H19 ICR sequence on the maternal allele enables binding of the CTCF protein, which is required to prevent enhancers from acting on maternal IGF2, thus, repressing its expression .
Imprinting disorders have shown the importance of the adequate H19 ICR methylation: hypomethylation results in downregulation of IGF2 and biallelic expression of H19, leading to a growth restriction disorder, Silver-Russell syndrome. By contrast, hypermethylation of H19 ICR leads to overexpression of IGF2, downregulation of H19, and consequently fetal over-growth known as Beckwith-Wiedemann syndrome [24, 25]. Hypomethylation of the sixth CTCF binding site (CTCF6) at the H19 ICR has been previously associated with placental cells , buccal epithelium cells , and cord blood mononuclear cells  in human pregnancies conceived by IVF or ICSI. Also, increased inter- and intra-individual variation in allele-specific DNA methylation and decreased IGF2 and H19 expression have been observed in placental tissue of in vitro conceived children .
In our recent study, we observed a single nucleotide polymorphism rs10732516 G/A in CTCF6, which associated with the distinct DNA methylation profiles of H19 ICR in human placenta . Moreover, when the samples were divided in four groups according to the genotype (rs10732516 G/G, paternal G/maternal A, paternal A/maternal G, and A/A), we observed decreased methylation level in alcohol-exposed placentas of paternal A/maternal G genotype, but not in paternal G/maternal A genotype. Surprisingly, alcohol exposure associated with decreased head circumference in all genotypes except A/A, in which increased head circumference was observed . Interestingly, previous studies have shown parent-of-origin associations between birth weight and polymorphisms rs4929984 and rs2071094, which both are in linkage to rs10732516 [28, 29].
Although prenatal alcohol exposure and IVF are very different environmental factors, they both have been associated with growth-restricted phenotype of newborns. To examine if there is similar genotype-specific decreased DNA methylation level at the IGF2/H19 locus caused by ART, we collected placental tissue from fresh embryo transfer (FRESH) and frozen embryo transfer (FET)-derived pregnancies of Finnish and Estonian couples (Table 1; Additional file 1: Table S1). We compared them to placentas of spontaneous, naturally conceived pregnancies. We explored the methylation levels of H19 ICR and H19 DMR in placenta. Owing to the previously detected ART-associated changes in global DNA methylation levels , also long interspersed nuclear elements (LINE-1) were examined. In addition to placenta, we collected umbilical cord blood to explore if we could see similar changes in DNA methylation in both extra embryonic placental cells and embryonic blood cells of the newborns. Potential genotype-specific effects of ART on the newborns’ phenotype were studied by using international growth standards .
Significant differences between mothers in studied groups (controls, FRESH- and FET-derived pregnancies) were observed in age (P < 0.0001, two-way ANOVA) and parity (P = 0.003, two-way ANOVA), but not in maternal BMI (Table 1). Finnish and Estonian mothers differed significantly in age (P = 0.005, two-way ANOVA) and parity (P = 0.04, two-way ANOVA); however, the interaction effect was not significant (P = 0.8, P = 1, respectively). There was a significant difference in gestational age between the study groups (P = 0.002, two-way ANOVA), as well as between the Finnish and Estonian study populations (P = 0.01, two-way ANOVA), although the interaction effect was not significant (P = 0.7). Furthermore, ART had no effect on the 5 min Apgar score which is used to evaluate the vitality of the newborn at birth.
Birth weight, birth length, and head circumference were examined for FRESH and FET derived as well as control newborns using international growth standards . The placental weights differed significantly between the controls, FRESH- and FET-derived newborns (P = 0.04, two-way ANOVA) (Table 1, Fig. 2a). Although there was a significant difference in the placental weights between Finnish and Estonian newborns (P < 0.001, two-way ANOVA), the interaction effect was not significant (P = 0.4) and the data was combined. Placentas of FET-derived pregnancies were heavier compared to FRESH-derived placentas and control placentas (P < 0.001, P = 0.01, respectively, Bonferroni post hoc). Also FRESH-derived placentas were lighter compared to controls (P = 0.01, Bonferroni post hoc) (Fig. 2a). According to the international growth standards, the standard deviations (SDs) of birth weight (Fig. 2b) or head circumference did not differ significantly between the Finnish and Estonian newborns or between the studied groups. Birth length differed significantly between Finnish and Estonian newborns (P = 0.001, two-way ANOVA), although the interaction effect was not significant (P = 0.9) and no difference could be observed between the study groups (P = 0.2). We did not observe significant differences between sexes, when all samples were compared.
DNA methylation profiles at H19 ICR, H19 DMR, and LINE-1
To explore the potential association between ART and placental DNA methylation changes, we compared the methylation levels of H19 imprinting control region (ICR) and H19 differentially methylated region (DMR) by EpiTYPER (Sequenom). We did not observe differences between ART and control placentas (Additional file 2: Table S1). We also determined the effects of ART on global methylation level in placenta by examining methylation in LINE-1 by EpiTYPER. However, we did not observe any significant alterations in global methylation level either (Additional file 2: Table S1).
Genotype-specific DNA methylation at H19 ICR and H19 DMR by EpiTYPER
Owing to the genotype-specific DNA methylation profiles of CTCF6 at H19 ICR , we divided our samples into four groups according to the genotype: rs10732516 G/G, paternal G/maternal A (patG/matA), paternal A/maternal G (patA/matG), and A/A. The allele frequencies of this polymorphism are almost equal in Finnish population (G = 0.47, A = 0.53) , and there were no differences in the prevalence of rs10732516 genotypes between controls and ART-derived samples in this study (X2(3) = 5.52, P = 0.138, chi-square test).
We compared first the genotype-specific methylation levels of placental CTCF6 at H19 ICR and H19 DMR by EpiTYPER to explore potential effects of ART. We did not see any genotype-specific differences between control and ART samples at H19 ICR (Additional file 2: Table S1). At the H19 DMR, we observed increased methylation level in CpG_3 and CpG_16 units in A/A genotype of ART samples (nominal P values: P = 0.03 and P = 0.05, respectively, Student’s t test), but changes were not significant after Bonferroni multiple testing correction.
Genotype-specific DNA methylation at H19 ICR by bisulfite sequencing
We also compared genotype-specific methylation levels of CTCF6 at H19 ICR between control and ART placentas by traditional bisulfite sequencing. To discern maternal and paternal alleles, we used only heterozygous samples (patG/matA and patA/matG). We observed a bias in PCR product: hypomethylated maternal allele of patA/matG genotype was amplified more efficiently compared to hypermethylated paternal allele. Owing to that, we counted the average methylation percentages separately for both alleles and then calculated the total methylation level for each CpG sites (CpG_1-27). We observed similar, but much more prominent common genotype-specific methylation profiles in placenta as we detected by EpiTYPER (Additional file 2: Table S1 and S2).
When comparing genotype-specific DNA methylation within heterozygotes (patG/matA and patA/matG) controls to ART samples, we observed decreased methylation level at sites CpG_1-3, CpG_5, CpG_14, and CpG_24 patA/matG genotype in the ART placentas (nominal P values 0.008, 0.02, 0.001, 0.013, 0.013, and 0.029, respectively, Mann-Whitney) (Fig. 3d). Instead of hypomethylation, we observed increased methylation level at site CpG_26 in patG/matA genotype (nominal P value = 0.041, Mann-Whitney). However, changes in methylation level were not significant after multiple testing correction. We did not see similar trend of decreased methylation in the patA/matG genotype by EpiTYPER method (Additional file 2: Table S1 and S2), which could be explained by the amplification bias in PCR.
Allele-specific DNA methylation at H19 ICR
We next assessed the allele-specific methylation levels of CTCF6 at H19 ICR in heterozygous genotypes (patG/matA and patA/matG) in placenta by bisulfite sequencing. When comparing the methylation levels of paternal and maternal alleles separately in control and ART samples, we observed consistently decreased methylation level in paternal allele of patA/matG genotype at sites CpG_1-5, CpG_14, CpG_24, and CpG_25 in ART placentas (nominal P values 0.013, 0.013, 0.001, 0.029, 0.013, 0.003, 0.005, and 0.029, respectively, Mann-Whitney) (Fig. 3e, Additional file 2: Table S2). Conversely, increased methylation level was observed at site CpG_26 in the maternal allele of patA/matG genotype (nominal P value = 0.005, Mann Whitney) (Fig. 3f), and both increased and decreased methylation at sites CpG_24 and CpG_26, respectively, in the paternal allele of patG/matA genotype (nominal P value = 0.026 and 0.015, respectively, Mann-Whitney) (Fig. 3b). However, changes in methylation levels were not significant after multiple testing correction.
To see if similar decreased methylation level in patA/matG genotype can be seen also in the blood, we examined white blood cells (WBCs) of newborns’ umbilical cord blood from the same ART-derived pregnancies. However, we did not observe similar decreased methylation level in the paternal allele of patA/matG genotype in ART-derived WBCs as we saw in placental tissue (Fig. 4, Additional file 2: Table S3). Conversely, a subtle but consistent increased methylation level in ART-derived WBCs was detected. The methylation level of CpG_4 site in the paternal allele of ART samples was clearly increased (P = 0.03, Mann-Whitney), but the difference is not significant after multiple testing correction.
Genotype-specific phenotypes of newborns
Finally, we assessed genotype-specific phenotypes of newborns by using international growth standards. Genotype-specific examination revealed differences in FRESH-derived newborns’ birth weight and head circumference (P = 0.04 and P = 0.004, respectively, one-way ANOVA). FRESH-derived newborns with G/G genotype were heavier (P = 0.04, Bonferroni post hoc) and had larger head circumference (P = 0.002, Bonferroni post hoc) compared to newborns with A/A genotype (Fig. 5a, b). We did not observe significant differences in birth length. Genotype-specific differences were not compared in the FET-derived newborns since the sample size was too low in the heterozygous genotypes. We also saw that the placental weight and birth weight differed significantly between controls, FRESH-derived and FET-derived newborns in the A/A genotype (P = 0.024 and P = 0.006, one-way ANOVA) (Fig. 5a, c). Both the placentas and newborns of FET-derived pregnancies were heavier than FRESH-derived pregnancies (P = 0.02 and P = 0.004, respectively, Bonferroni post hoc). We did not see similar differences between the groups in the G/G genotype.
Owing to the decreased birth weight and decreased methylation level at the H19 ICR associated previously with IVF treatments, we focused on the imprinted IGF2/H19 locus, which is crucial for normal placental and embryonic growth. In our previous study, a polymorphism rs10732516 at this locus associated with a genotype-specific trend in placental DNA methylation and head circumference of prenatally alcohol-exposed newborns . Due to the growth-restricted phenotype in both prenatal alcohol exposure and IVF, we explored if these rather different environmental factors could associate with similar changes. We observed consistently decreased DNA methylation at the sixth binding site (CTCF6) of H19 ICR in ART-derived placentas, which is consistent with previous ART study . Interestingly, the decreased methylation level was detected only in the paternal allele in rs10732516 patA/matG genotype of two studied heterozygous genotypes, and alteration was even more profound than in the alcohol-exposed placentas in our previous study. This suggests that the effect of ART on DNA methylation in placenta is genotype-specific.
We did not observe similar changes in the methylation level of paternal allele in WBCs from cord blood of ART-derived newborns. This could be explained by the more advanced differentiation stage of the extraembryonic trophoblast cells compared to inner cell mass, the location of the cells in the blastocyst, or the better DNA methylation repairing and maintaining mechanisms of the embryonic cells. However, our result is consistent with the earlier study, where decreased methylation in this same region was detected only in mononuclear cells, not in WBCs . This suggests that these genotype-specific changes in methylation have not occurred in blood cells in the early embryonic development or they are not fixed in all cell types.
We are aware of some limitations in this study. Traditional bisulfite sequencing is useful only for allele-specific examination of heterozygous genotypes, and thus, the methylation information about homozygous genotypes is lacking. Furthermore, we did not observe similar trend of decreased methylation by EpiTYPER, indicating that as a PCR-based method, it is not a convenient method to detect relatively small but consistent allele specific methylation changes in this specific imprinted region. Also, based on these results, it is not possible to see if decreased methylation at the H19 ICR is caused by ART or infertility.
The number of FET samples was too low to explore the phenotypic effects on each genotype, but the genotype-specific variation in the birth weight and head circumference of FRESH-derived newborns suggest that the polymorphism could associate with the growth. Genotype-specific examination also revealed that in the rs10732516 A/A genotype, the placental weight and birth weight (SD) of controls, FRESH- and FET-derived newborns, differed from each other. Interestingly, we did not see the same in the G/G genotype. Both the placentas and newborns of FET-derived pregnancies were the heaviest in the A/A genotype. This is consistent with the increased head circumference of prenatally alcohol-exposed children  as well as the strongest growth phenotype of infantile hemangiomas , which both associate with this specific genotype.
Both genotype-specific methylation profiles and phenotypic findings suggest that rs10732516 polymorphism associates with the effects of ART in a parent-of-origin manner. The polymorphism locates on the binding sequence of CTCF protein, and allele A deletes a CpG binding site for a methyl group. Whether the A allele affects slightly the binding efficiency of CTCF protein and consequently makes A/A genotype particularly sensitive to environmental factors, it needs to be clarified in functional studies. More studies are also needed to find out if changes in this locus have occurred already in the very beginning of the embryonic development, in the period of epigenetic reprogramming, and the causality of the alterations: could the genotype-specific changes in DNA methylation affect the gene expression and thus the phenotype of developing embryo. Owing to the genotype-specific methylation changes at the H19 ICR in ART-derived placentas and previous associations between ART and imprinting disorders, it would be interesting to find out if the prevalence of imprinting disorders associates with the rs10732516 G/A polymorphism.
Study design and sample collection
Couples applied to fertilization treatment in the Reproductive Medicine Unit of Helsinki University Central Hospital, Finland or Fertility Clinic of the Family Federation of Finland or Tartu University Hospital were recruited to this study. IVF or ICSI have been used in the treatments. The conception has been done using fresh embryo transfer (FRESH) or frozen embryo transfer (FET). Placental and cord blood samples from 47 Finnish cases (29 FRESH: 23 IVF/6 ICSI and 18 FET: 12 IVF/6 ICSI) are collected during years 2013–2017 and 39 Estonian cases (33 FRESH: 21 IVF/12 ICSI and 6 FET: 3 IVF/3 ICSI) 2016–2017. Spontaneously conceived 100 Finnish controls have been collected during years 2013–2015 in Helsinki University Central Hospital, Finland , and 57 Estonian controls in Tartu University Hospital. All the samples were Caucasian origin, from Finnish and Estonian newborns. Sample information and variation between Finnish and Estonian samples are shown in Table 1 and Additional file 1: Table S1.
The placental biopsies (1 cm3) and umbilical cord blood samples were collected immediately after delivery. The placental biopsies were collected from the fetal side of placenta within a radius of 2–3 cm from the umbilical cord, rinsed in cold 1× PBS and stored in RNAlater® (Thermo Fisher Scientific, Vilnius, Lithuania) at − 80 °C. White blood cells (WBCs) were extracted as soon as possible, at latest 16 h after birth (Additional file 3: Protocol S1).
Birth weight (g), birth length (cm), and head circumference (cm) were examined for both Finnish and Estonian newborns using international growth standards, the Fentom Preterm Growth Chart by PediTools (http://peditools.org/), in which the gestational age at birth and sex are considered when calculating the SD (z-score) of birth measures . This chart has also been used previously for full-term deliveries [34, 35]. Measures deviating more than ± 2 SDs are commonly considered abnormal.
Placental genomic DNA was extracted by commercial QIAamp Fast DNA Tissue Kit (Finnish samples, Qiagen, Valencia, CA, USA) or PureLink Genomic DNA Kit (Estonian samples, Invitrogen, Life Technologies, USA). The extractions were done from one to four pieces (on average from three pieces) of placental tissue. The total DNA methylation levels of H19 ICR (CTCF6), H19 DMR, and LINE-1 regions in placental samples were measured by MassARRAY EpiTYPER (SEQUENOM Inc.) based on matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry. First, DNA (1000 ng) was bisulfite converted (EZ-96 DNA Methylation™ kit, Zymo Research, Irvine, CA, USA) and PCR was performed in three independent 10 μl reactions using HotStar PCR kit (Qiagen, Valencia, CA, USA) according to manufacturer’s instructions. Primers for the regions of interest were obtained from previous publications [36, 37] (Additional file 3: Table S2). The EpiTYPER measurements were done for pooled PCR reactions. Altogether, 60 controls and 62 ART-derived (48 FRESH: 33 IVF/15 ICSI and 14 FET: 8 IVF/6 ICSI) Finnish and Estonian placental samples were analyzed by EpiTYPER.
To find out the allele-specific methylation profiles and to confirm the EpiTYPER results as well as genotypes, the CTCF6 at H19 ICR of heterozygous ART-derived placental samples with patG/matA genotype (4 FRESH: 4 IVF, and 2 FET: 1 IVF/1 ICSI) and patA/matG genotype (7 FRESH: 4 IVF/3 ICSI, and 1 FET: IVF), and eight WBC samples with patA/matG genotype (4 controls and 4 ART-derived samples: 3 FRESH: 1 IVF/2 ICSI, and 1 FET: IVF) were subjected to bisulfite sequencing. All the samples were from Finnish newborns. The control placental samples had been published previously . Due to heterozygosity and imprinting, the paternal and maternal alleles could be distinguished. Two separate bisulfite conversions were performed for 500 ng of genomic DNA (EZ DNA Methylation™ kit, Zymo Research, Irvine, CA, USA) and pooled afterwards. To avoid possible PCR bias, three independent 20 μl PCR reactions (HotStar PCR kit, Qiagen, Valencia, CA, USA) were performed per sample. Primers were obtained from previous publication and allowed to detect the polymorphism in units CpG_17,18,19,20  (Additional file 2: Table S3). PCR reactions were gel isolated, and the three reactions of each sample were pooled and purified using NucleoSpin Gel and PCR Clean-up Kit (Macherey-Nagel, Düren, Germany). The purified PCR fragments were ligated into pGEM®-T Easy Vector (Promega, Madison, WI, USA) and cloned by standard protocol. The recombinant-DNA clones were purified using NucleoSpin® Plasmid EasyPure kit (Macherey-Nagel, Düren, Germany) according to manufacturer’s instructions. Fifty to eighty clones of each individual were sequenced. The sequences were analyzed by BIQ Analyzer  excluding the clones with lower than 90% conversion rate from the dataset.
Placental samples were genotyped by Sanger sequencing. According to our sequencing analyses, the heterozygous samples were able to distinguish from each other due to uneven amplification and hence different signal levels of the alleles. The peak of rs10732516 A in sequence of patA/matG genotype was lower compared to patG/matA genotype. One 20 μl PCR reaction was performed for each sample using commercial HotStar PCR kit (Qiagen Valencia, CA, USA) with 100–300 ng of template DNA. Primers were designed by using NCBI/Primer Blast (Additional file 3: Table S2). PCR products were purified with SAP treatment (FastAP Thermosensitive Alkaline Phosphatase (1 U/μL), Thermo Scientific, Waltham, MA, USA) according to manufacturer’s instructions. The genotypes of samples that were analyzed by EpiTYPER could be confirmed by detecting genotype-specific fragmentation and distinct methylation levels in unit CpG_10 of H19 ICR. The methylation level in patG/matA was ~ 0.80, in G/G ~ 0.30, in patA/matG ~ 0.02, and in A/A there was no value.
Statistical analyses were conducted using either SPSS software for Windows version 22.0 (NY, USA) or GraphPad Prism 7 software (GraphPad Software, Inc., La Jolla, CA, USA). All data are expressed as the mean with ±SD for a normal distribution of variables. Samples were divided into four groups according to the genotype and the chi-square test was used to compare the prevalence of the rs1072516 in the control and ART samples. The non-parametric Mann–Whitney test was used to compare the methylation level of CpG sites analyzed by bisulfite sequencing. Student’s t test was used to compare CpG units analyzed by EpiTYPER. In the methylation analysis, the nominal P value was considered significant when < 0.05 and Bonferroni correction was used for multiple testing correction. Two-way ANOVA, followed by Bonferroni post hoc test when significant, was used to identify the differences among the study groups as well as to eliminate the interaction effect if significant differences between the Finnish and Estonian newborns were observed.
Assisted reproductive technology
Zinc-finger protein CCCTC-binding factor
The sixth binding site for CTCF
Differentially methylated region
Frozen embryo transfer
Fresh embryo transfer
Imprinting control region
Intracytoplasmic sperm injection
Insulin-like growth factor 2
In vitro fertilization
Long interspersed nuclear elements
White blood cell
Qin JB, Sheng XQ, Wu D, Gao SY, You YP, Yang TB, Wang H. Worldwide prevalence of adverse pregnancy outcomes among singleton pregnancies after in vitro fertilization/intracytoplasmic sperm injection: a systematic review and meta-analysis. Arch Gynecol Obstet. 2017;295:285–301.
Jackson RA, Gibson KA, Wu YW, Croughan MS. Perinatal outcomes in singletons following in vitro fertilization: a meta-analysis. Obstet Gynecol. 2004;103:551–63.
Pandey S, Shetty A, Hamilton M, Bhattacharya S, Maheshwari A. Obstetric and perinatal outcomes in singleton pregnancies resulting from IVF/ICSI: a systematic review and meta-analysis. Hum Reprod Update. 2012;18:485–503.
Malchau SS, Loft A, Henningsen AK, Nyboe Andersen A, Pinborg A. Perinatal outcomes in 6,338 singletons born after intrauterine insemination in Denmark, 2007 to 2012: the influence of ovarian stimulation. Fertil Steril. 2014;102:1110–6.
Maher ER, Brueton LA, Bowdin SC, Luharia A, Cooper W, Cole TR, Macdonald F, Sampson JR, Barratt CL, Reik W, et al. Beckwith-Wiedemann syndrome and assisted reproduction technology (ART). J Med Genet. 2003;40:62–4.
Lazaraviciute G, Kauser M, Bhattacharya S, Haggarty P, Bhattacharya S. A systematic review and meta-analysis of DNA methylation levels and imprinting disorders in children conceived by IVF/ICSI compared with children conceived spontaneously. Hum Reprod Update. 2014;20:840–52.
Pelkonen S, Koivunen R, Gissler M, Nuojua-Huttunen S, Suikkari AM, Hyden-Granskog C, Martikainen H, Tiitinen A, Hartikainen AL. Perinatal outcome of children born after frozen and fresh embryo transfer: the Finnish cohort study 1995–2006. Hum Reprod. 2010;25:914–23.
Wennerholm UB, Henningsen AK, Romundstad LB, Bergh C, Pinborg A, Skjaerven R, Forman J, Gissler M, Nygren KG, Tiitinen A. Perinatal outcomes of children born after frozen-thawed embryo transfer: a Nordic cohort study from the CoNARTaS group. Hum Reprod. 2013;28:2545–53.
Maheshwari A, Pandey S, Amalraj Raja E, Shetty A, Hamilton M, Bhattacharya S. Is frozen embryo transfer better for mothers and babies? Can cumulative meta-analysis provide a definitive answer? Hum Reprod Update. 2018;24:35–58.
Pinborg A, Loft A, Romundstad LB, Wennerholm UB, Söderström-Anttila V, Bergh C, Aittomäki K. Epigenetics and assisted reproductive technologies. Acta Obstet Gynecol Scand. 2016;95:10–5.
Sunkara SK, La Marca A, Seed PT, Khalaf Y. Increased risk of preterm birth and low birthweight with very high number of oocytes following IVF: an analysis of 65 868 singleton live birth outcomes. Hum Reprod. 2015;30:1473–80.
Kamath MS, Kirubakaran R, Mascarenhas M, Sunkara SK. Perinatal outcomes after stimulated versus natural cycle IVF: a systematic review and meta-analysis. Reprod BioMed Online. 2018;36:94–101.
Reik W, Dean W, Walter J. Epigenetic reprogramming in mammalian development. Science. 2010;293:1089–93.
Yin Y, Morgunova E, Jolma A, Kaasinen E, Sahu B, Khund-Sayeed S, Das PK, Kivioja T, Dave K, Zhong F, et al. Impact of cytosine methylation on DNA binding specificities of human transcription factors. Science. 2017;356:eaaj2239.
Nelissen EC, Dumoulin JC, Daunay A, Evers JL, Tost J, van Montfoort AP. Placentas from pregnancies conceived by IVF/ICSI have a reduced DNA methylation level at the H19 and MEST differentially methylated regions. Hum Reprod. 2013;28:1117–26.
Loke YJ, Galati JC, Saffery R, Craig JM. Association of in vitro fertilization with global and IGF2/H19 methylation variation in newborn twins. J Dev Orig Health Dis. 2015;6:115–24.
Castillo-Fernandez JE, Loke YJ, Bass-Stringer S, Gao F, Xia Y, Wu H, Lu H, Liu Y, Wang J, Spector TD, et al. DNA methylation changes at infertility genes in newborn twins conceived by in vitro fertilisation. Genome Med. 2017;9:1–15.
Mann MR, Lee SS, Doherty AS, Verona RI, Nolen LD, Schultz RM, Bartolomei MS. Selective loss of imprinting in the placenta following preimplantation development in culture. Development. 2004;131:3727–35.
Market-Velker BA, Fernandes AD, Mann MR. Side-by-side comparison of five commercial media systems in a mouse model: suboptimal in vitro culture interferes with imprint maintenance. Biol Reprod. 2010;83:938–50.
DeChiara TM, Robertson EJ, Efstratiadis A. Parental imprinting of the mouse insulin-like growth factor II gene. Cell. 1991;64:849–59.
Gabory A, Ripoche MA, Le Digarcher A, Watrin F, Ziyyat A, Forne T, Jammes H, Ainscough JF, Surani MA, Journot L, et al. H19 acts as a trans regulator of the imprinted gene network controlling growth in mice. Development. 2009;136:3413–21.
Phillips JE, Corces VG. CTCF: master weaver of the genome. Cell. 2009;137:1194–211.
Hark AT, Schoenherr CJ, Katz DJ, Ingram RS, Levorse JM, Tilghman SM. CTCF mediates methylation-sensitive enhancer-blocking activity at the H19/Igf2 locus. Nature. 2000;405:486–9.
Soejima H, Higashimoto K. Epigenetic and genetic alterations of the imprinting disorder Beckwith–Wiedemann syndrome and related disorders. J Hum Genet. 2013;58:402–9.
Gicquel C, Rossignol S, Cabrol S, Houang M, Steunou V, Barbu V, et al. Epimutation of the telomeric imprinting center region on chromosome 11p15 in Silver-Russell syndrome. Nat Genet. 2005;37:1003–7.
Turan N, Katari S, Gerson LF, Chalian R, Foster MW, Gaughan JP, Coutifaris C, Sapienza C. Inter- and intra-individual variation in allele-specific DNA methylation and gene expression in children conceived using assisted reproductive technology. PLoS Genet. 2010;6:e1001033.
Marjonen H, Kahila H, Kaminen-Ahola N. rs10732516 polymorphism at the IGF2/H19 locus associates with a genotype-specific trend in placental DNA methylation and head circumference of prenatally alcohol-exposed newborns. Hum Reprod Open. 2017;(3):hox014.
Adkins RM, Somes G, Morrison JC, Hill JB, Watson EM, Magann EF, Krushkal J. Association of birth weight with polymorphisms in the IGF2, H19, and IGF2R genes. Pediatr Res. 2010;68:429–34.
Petry CJ, Seear RV, Wingate DL, Acerini CL, Ong KK, Hughes IA, Dunger DB. Maternally transmitted foetal H19 variants and associations with birth weight. Hum Genet. 2011;130:663–70.
Ghosh J, Coutifaris C, Sapienza C, Mainigi M. Global DNA methylation levels are altered by modifiable clinical manipulations in assisted reproductive technologies. Clin Epigenetics. 2017;6:14.
Fenton TR, Nasser R, Eliasziw M, Kim JH, Bilan D, Sauve R. Validating the weight gain of preterm infants between the reference growth curve of the fetus and the term infant. BMC Pediatr. 2013;13:92.
Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, Korbel JO, Marchini JL, McCarthy S, McVean GA, et al. Genomes Project C: a global reference for human genetic variation. Nature. 2015;526:68–74.
Schultz B, Yao X, Deng Y, Waner M, Spock C, Tom L, Persing J, Narayan D. A common polymorphism within the IGF2 imprinting control region is associated with parent of origin specific effects in infantile hemangiomas. PLoS One. 2015;10:e0113168.
Gilbert-Diamond D, Emond JA, Baker ER, Korrick SA, Karagas MR. Relation between in utero arsenic exposure and birth outcomes in a cohort of mothers and their newborns from New Hampshire. Environ Health Perspect. 2016;124:1299–307.
Vanker A, Barnett W, Workman L, Nduru PM, Sly PD, Gie RP, Zar HJ. Early-life exposure to indoor air pollution or tobacco smoke and lower respiratory tract illness and wheezing in African infants: a longitudinal birth cohort study. Lancet Planet Health. 2017;1:e328–36.
Ollikainen M, Smith KR, Joo EJ, Ng HK, Andronikos R, Novakovic B, Abdul Aziz NK, Carlin JB, Morley R, Saffery R, et al. DNA methylation analysis of multiple tissues from newborn twins reveals both genetic and intrauterine components to variation in the human neonatal epigenome. Hum Mol Genet. 2010;19:4176–88.
Wang L, Wang F, Guan J, Le J, Wu L, Zou J, Zhao H, Pei L, Zheng X, Zhang T. Relation between hypomethylation of long interspersed nucleotide elements and risk of neural tube defects. Am J Clin Nutr. 2010;91:1359–67.
Coolen MW, Statham AL, Gardiner-Garden M, Clark SJ. Genomic profiling of CpG methylation and allelic specificity using quantitative high-throughput mass spectrometry: critical evaluation and improvements. Nucleic Acids Res. 2007;35:e119.
Bock C, Reither S, Mikeska T, Paulsen M, Walter J, Lengauer T. BiQ analyzer: visualization and quality control for DNA methylation data from bisulphite sequencing. Bioinformatics. 2005;21:4067–8.
We thank all participants in Finland and Estonia for their invaluable contributions to the study and research nurses Teija Karkkulainen and Riikka Hiltunen.
The study was funded by Faculty of Medicine, University of Helsinki, and Helsinki University Hospital, Estonian Ministry of Education and Research (IUT34-16), Enterprise Estonia (EU48695), Horizon 2020 innovation program (WIDENLIFE, 692065), European Union’s FP7 Marie Curie Industry-Academia Partnerships and Pathways funding (IAPP, SARM, EU324509), and MSCA-RISE-2015 project MOMENDO (691058).
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional files.
Ethics approval and consent to participate
The study was approved by the Ethics Committee of Helsinki University Central Hospital (386/13/03/03/2012 and 285/13/03/03/2013) and Research Ethics Committee of the University of Tartu (256/M-17). Informed consent was obtained from all participants.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Information of ART samples. (PDF 165 kb)
Table S1. DNA methylation levels of H19 ICR (CTCF6), H19 DMR and LINE-1 in control and ART-derived placentas by EpiTYPER method. Methylation average values with SDs (±) of CpG units are presented. Gray boxes present significant methylation level difference at CpG sites in H19 DMR between ART and control samples within A/A genotype (nominal p-value < 0.05, Student’s t-test). Table S2. Total, rs10732516 genotype-specific (patG/matA and patA/matG) and allele-specific (paternal and maternal) DNA methylation levels at H19 ICR (CTCF6) in control and ART-derived placentas by traditional bisulfite sequencing. Gray boxes present significant methylation level difference at CpG sites between control and ART samples within genotypes and alleles (nominal p-value < 0.05, Mann-Whitney). Total methylation levels: total means of both genotypes and both alleles are included. Genotype-specific methylation levels: total means of both alleles are included. Allele-specific methylation levels: maternal and paternal alleles are presented separately. Table S3. rs10732516 patA/matG genotype- and allele-specific DNA methylation levels at H19 ICR1 (CTCF6) in controls and ART-derived newborns’ white blood cells (WBCs) in cord blood by traditional bisulfite sequencing. Genotype-specific methylation levels: total means of both alleles are included. Allele-specific methylation levels: maternal and paternal alleles are presented separately. (PDF 145 kb)
Protocol S1. Extraction of total white blood cells (WBC) from EDTA-tube. Table S2. Primers (A) and PCR reactions (B). (PDF 306 kb)
About this article
Cite this article
Marjonen, H., Auvinen, P., Kahila, H. et al. rs10732516 polymorphism at the IGF2/H19 locus associates with genotype-specific effects on placental DNA methylation and birth weight of newborns conceived by assisted reproductive technology. Clin Epigenet 10, 80 (2018). https://doi.org/10.1186/s13148-018-0511-2