- Open Access
Association of four imprinting disorders and ART
Clinical Epigenetics volume 11, Article number: 21 (2019)
Human-assisted reproductive technologies (ART) are a widely accepted treatment for infertile couples. At the same time, many studies have suggested the correlation between ART and increased incidences of normally rare imprinting disorders such as Beckwith-Wiedemann syndrome (BWS), Angelman syndrome (AS), Prader-Willi syndrome (PWS), and Silver-Russell syndrome (SRS). Major methylation dynamics take place during cell development and the preimplantation stages of embryonic development. ART may prevent the proper erasure, establishment, and maintenance of DNA methylation. However, the causes and ART risk factors for these disorders are not well understood.
A nationwide epidemiological study in Japan in 2015 in which 2777 pediatrics departments were contacted and a total of 931 patients with imprinting disorders including 117 BWS, 227 AS, 520 PWS, and 67 SRS patients, were recruited. We found 4.46- and 8.91-fold increased frequencies of BWS and SRS associated with ART, respectively. Most of these patients were conceived via in vitro fertilization (IVF) and intracytoplasmic sperm injection (ICSI), and showed aberrant imprinted DNA methylation. We also found that ART-conceived SRS (ART-SRS) patients had incomplete and more widespread DNA methylation variations than spontaneously conceived SRS patients, especially in sperm-specific methylated regions using reduced representation bisulfite sequencing to compare DNA methylomes. In addition, we found that the ART patients with one of three imprinting disorders, PWS, AS, and SRS, displayed additional minor phenotypes and lack of the phenotypes. The frequency of ART-conceived Prader-Willi syndrome (ART-PWS) was 3.44-fold higher than anticipated. When maternal age was 37 years or less, the rate of DNA methylation errors in ART-PWS patients was significantly increased compared with spontaneously conceived PWS patients.
We reconfirmed the association between ART and imprinting disorders. In addition, we found unique methylation patterns in ART-SRS patients, therefore, concluded that the imprinting disorders related to ART might tend to take place just after fertilization at a time when the epigenome is most vulnerable and might be affected by the techniques of manipulation used for IVF or ICSI and the culture medium of the fertilized egg.
Human-assisted reproductive technologies (ART) are becoming increasingly common due to late marriage and improvements in medical technology in developed countries, including Japan . The majority of studies published previously have suggested that babies conceived after ART have increased incidences of normally rare imprinting disorders such as Beckwith-Wiedemann syndrome (BWS), Angelman syndrome (AS), Prader-Willi syndrome (PWS), and Silver-Russell syndrome (SRS) [2,3,4,5,6,7]. On the other hand, some studies found no relation between ART and imprinting disorders [8, 9]. Imprinting disorders are caused by genetic defects or epigenetic mutations (DNA methylation); i.e., aberrant DNA methylation of differentially methylated regions (DMRs) that regulate allele-specific expression of imprinted genes . The relationship between ART and aberrant genomic imprinting is still unclear. However, many studies have suggested that ovarian stimulation, culture media used for gametes and early embryos, in vitro fertilization (IVF) and intracytoplasmic sperm injection (ICSI) manipulations, and the freezing and thawing of embryos may prevent the proper establishment and maintenance of genomic imprinting and cause imprinting disorders [11, 12].
In 2009, we performed a nationwide epidemiological study in Japan to determine the frequencies of four imprinting disorders (BWS, AS, PWS, and SRS) and found that the frequencies of ART-conceived BWS (ART-BWS) and ART-SRS patients were nearly tenfold higher than anticipated . In addition, the majority of ART-SRS and ART-BWS patients showed aberrant methylation not only in domains responsible for imprinting disorders but also in additional multiple imprinted loci. Other reports have also described multi-locus imprinting disturbances after ART, supporting our hypothesis [14,15,16]. Furthermore, in addition to aberrant DNA methylation patterns at multiple imprinted loci, there are also mosaic methylation errors. We assumed that the imprinting errors occurred after fertilization rather than in the gametes. However, it is unknown who is predisposed to the imprinting disorders and which factor(s) of ART cause them.
In this study, we conducted another Japanese nationwide epidemiological study of the imprinting disorders in 2015 and reconfirmed the tight association between ART and imprinting disorders. Next, we studied new patients with SRS and found significant differences in the genome-wide DNA methylation between the ART patients and spontaneously conceived (Sp) patients.
Relationship between the imprinting disorders and ART
We conducted a Japanese nationwide epidemiological study of four imprinting disorders to determine associations with ART. We contacted a total of 2777 pediatric departments of all hospitals that were identified based on a hospital lists. In all, 1957 departments (70.5%) responded to the first-stage survey questionnaire. We inquired about detailed clinical information and examined the relationship with ART using a second-stage survey ensuring the exclusion of duplicates and a total of 931 patients with imprinting disorders were recruited. The number, age, and gender distributions of patients are given in Additional file 1. The numbers of patients with the four imprinting disorders slightly increased year by year, especially for BWS and SRS. However, there were no differences in the gender ratios of the imprinting disorders. We ascertained the frequencies of ART-BWS, ART-AS, ART-PWS, and ART-SRS via a questionnaire sent to doctors (Table 1). The observed/expected ratios for four ART-imprinting disorders were calculated by considering the live birth rate after ART (1.34%) from 1985 through 2015 in Japan (http://plaza.umin.ac.jp/~jsog-art/). For BWS and SRS, the number of ART patients were 4.46- and 8.91-fold higher than anticipated, accounting for 6.0% (7/117) of BWS and 11.9% (8/67) of SRS patients, respectively. These proportions were similar to those we found in 2009 . The number of ART-PWS patients was 3.44-fold higher than anticipated, but no increase in ART-AS patients was observed compared with that anticipated (1.8%, 4/227).
In the second questionnaire, we examined the results of genetic and DNA methylation tests (Table 2 and Additional file 2). In both ART-BWS and ART-SRS patients, DNA methylation error rates were 100% (4/4 and 5/5, respectively). Additionally, that in ART-BWS was significantly higher than in spontaneously conceived BWS (Sp-BWS) patients (p = 0.027). Nearly one-half of ART-BWS (42.9%, 3/7) and the majority of ART-SRS (75%, 6/8) patients received IVF or ICSI treatment (Additional file 3). The rate of maternal uniparental disomy of chromosome 15 (UPD(15)mat) in ART-PWS patients was significantly increased compared with that of Sp-PWS patients (p = 0.001). Next, we classified PWS patients into two categories based on maternal ages ≤ 37 years and ≥ 38 years and compared ART patients with Sp patients since the aneuploidy rate in the human embryo dramatically increases when the maternal age ≥ 38 years . The rate of DNA methylation errors in ART-PWS patients was significantly increased compared with that in Sp-PWS patients for maternal age ≤ 37 years, whereas the rate of UPD(15)mat in ART-PWS patients was significantly increased compared with that in Sp-PWS patients with maternal age ≥ 38 years (p = 0.021, Table 3 and Additional file 4). There was no correlation between maternal age and UPD in BWS, AS, and SRS patients (Additional file 5).
Phenotypic differences in patients conceived with ART
To investigate whether disease phenotypes in patients could be altered by ART, we compared the clinical features in detail. There were no major differences overall between ART patients and Sp patients for the four diseases, but there were some significant phenotypic differences (Table 4). For ART-SRS patients, the frequency of heart malformation was significantly higher (25.0%, 2/8; 3.4%, 2/59, p = 0.036) than in Sp-SRS patients. In addition, the frequency of mental retardation in ART-SRS was higher (62.5%, 5/8; 32.2%, 19/59, p = 0.124) than in Sp-SRS patients. Tumorigenesis and diabetes occurred in only one ART-SRS patient. For ART-AS patients, the frequencies of hyposomnia and ataxic gait were significantly higher (100%, 4/4; 46.2%, 103/223, p = 0.048) and lower (0%, 0/4; 60.1%, 134/223, p = 0.027), respectively, than in Sp-AS patients. For ART-PWS patients, the frequencies of dyschromatosis (91.7%, 22/24; 69.0%, 342/496, p = 0.020) and acromicria (87.5%, 21/24; 58.3%, 289/496, p = 0.005) were significantly higher and the frequency of obesity was lower (14.3%, 2/14; 49.4%, 205/415, p = 0.012) than in Sp-PWS patients. For BWS, there were some phenotypical differences between ART-BWS and Sp-BWS patients, but these were not significant.
Genome-wide methylation in ART-SRS patients
To investigate the genome-wide DNA methylation changes in ART-SRS patients, we performed reduced representation bisulfite sequencing (RRBS) to produce DNA methylomes derived from the peripheral blood of five ART-SRS patients with H19/IGF2 intergenic (IG)-DMR hypomethylation, five Sp-SRS patients with H19/IGF2 IG-DMR hypomethylation, and ten Sp-normal children (control) who were randomly selected. We obtained an average of 3,242,341 CpG cytosines per sample and analyzed 2,069,685 CpG cytosines covered in the autosomes of all samples (Additional file 6). Global evaluations of the DNA methylation levels of all CpG cytosines in the genome among individuals of the three groups were very similar (R > 0.975), and the correlation coefficients between individuals within the ART-SRS and Sp-SRS groups were greater than between individuals within the control group (Additional file 7). Among the mean methylation levels of the various genomic regions, the most methylated was in the ART-SRS group, followed in order by the Sp-SRS and control groups. The mean methylation levels of the gene bodies, exons, introns, intergenic regions, CpG island (CGI) shores, short interspersed nuclear elements (SINEs), long interspersed nuclear elements (LINEs), long terminal repeats (LTRs), and DNA repeats in the ART-SRS and Sp-SRS groups were significantly increased compared to the control group. In addition, the mean methylation levels of the CGI shelves and simple repeats in the ART-SRS group were significantly increased compared to the control group (Additional file 8).
Characterization of DMVs in ART-SRS patients
We annotated genomic regions and then identified the DNA methylation variation (DMVs) showing absolute methylation changes ≥ 7.5% and statistically significant differences (false discovery rate (FDR) < 0.05) in the ART-SRS and Sp-SRS groups compared with the control group (Fig. 1). The numbers of DMVs in the ART-SRS group were significantly increased compared with those in the Sp-SRS group in the promoter, gene body, CGI, CGI shore, CGI shelf, SINE, LTR, and simple repeat regions. In addition, the number of methylated DMVs was larger than that of demethylated DMVs in all regions. On the other hand, we found that there were no more DMVs in the ART-SRS and Sp-SRS groups than in the control group for SINE-VNTR-Alu (SVA).
To investigate genes related to the characteristics of SRS, we next focused on promoter regions. A total of 83 promoters (corresponding to 79 genes) showed DMVs; 37 promoters (36 genes) with DMVs only in the ART-SRS group, 14 promoters (13 genes) with DMVs only in the Sp-SRS group, and 32 promoters (30 genes) with DMVs in both the ART-SRS and Sp-SRS groups (Fig. 2a and Additional file 9). These promoters with DMVs were distributed in almost all autosomes (Fig. 2b). We used the Database for Annotation, Visualization and Integrated Discovery (DAVID) to interpret the biological meanings of lists of ART-SRS-specific, Sp-SRS-specific, and all genes with DMVs but there were no significant gene ontology (GO) terms. However, some genes with DMVs might be associated with clinical phenotypes of SRS; therefore, we describe these genes below (see “Discussion” section).
Germline-specific methylation patterns of identified DMVs in ART-SRS patients
To investigate the characteristics of DMVs identified in ART-SRS patients, we classified the DMVs into four categories: sperm-specific methylated regions ≥ 80% methylated in sperm and ≤ 20% methylated in oocytes, oocyte-specific methylated regions ≤ 20% methylated in sperm and ≥ 80% methylated in oocytes, both hypermethylated regions ≥ 80% methylated in both sperm and oocytes, and both hypomethylated regions ≤ 20% methylated in both sperm and oocytes, according to our previously reported data . The proportion of DMVs occupying sperm-specific methylated regions was the highest in the promoter, gene body, CGI, CGI shelf, SINE, LTR, DNA repeat, and simple repeat regions (Fig. 2c and Additional file 10). The proportion of DMVs occupying both hypermethylated regions often followed the proportion of DMVs occupying sperm-specific methylated regions. In addition, these DMVs were incompletely (mosaic) methylated. These data suggested that the sperm-specific methylated regions were more likely to show DMVs just after fertilization in ART. Therefore, we concluded that DMVs might be affected by the techniques of the manipulation used for IVF or ICSI and the culture medium of the fertilized egg.
Using a genome-wide approach, we demonstrated the characteristics of the DNA methylation in ART-SRS patients. They were as follows: (1) the numbers and locations of DMVs were larger and more widespread, respectively. However, some domains such as imprinted regions, SVAs, and retrotransposons had few DMVs. (2) The DMVs tended to more frequently be methylated than demethylated events. (3) The DMVs occurred more frequently in sperm-specific methylated regions. These results may suggest that ART caused the DNA methylation alterations. Major methylation dynamics take place during cell development and the preimplantation stages of embryonic development . ART may prevent the proper erasure, establishment, and maintenance of DNA methylation. Recently, we and others reported that the human paternal genome is globally demethylated after fertilization, which is due to ten–eleven translocation (TET)-mediated active demethylation [18, 20,21,22]. On the other hand, the maternal genome is demethylated to a much lesser extent in human blastocysts, which was different from the mouse. The pattern of incomplete DMVs we observed in ART-SRS patients might indicate that the changes occurred after fertilization rather than in the gamete. Furthermore, the increased frequency of DMVs in the paternal-specific methylation domain suggested that they occurred just after fertilization.
Previously, we examined the DNA methylation status of 22 germline DMRs (gDMRs) located within the imprinted loci in an ART-BWS and five ART-SRS patients compared to those of Sp-SRS patients . We observed aberrant DNA methylation patterns in the majority of ART patients (1) at multiple imprinted loci, (2) with both maternal and paternal gDMRs, (3) with both hypermethylation and hypomethylation events, and (4) mosaic methylation errors. These aberrant methylation patterns suggested that imprinting defects caused by the global demethylation after fertilization rather than in the gamete might be involved. Others also reported similar results [23, 24]. In this study, we hypothesized that the demethylation and maintenance mechanisms in the preimplantation stage might fail in ART and lead to DMVs and imprinting disorders such as BWS and SRS. Many reports have demonstrated the possibility that the process of ART might affect epigenetic alterations. For example, the phenome of large offspring syndrome has been described in bovine and sheep IVF . Use of sheep embryo cultures in vitro has led to the birth of overweight male offspring with phenotypic similarity to BWS and associated with loss of imprinted methylation at IGF2R . Some studies have shown that exposure of mouse embryos to different culture conditions can alter the gene expression and DNA methylation, which could result in abnormal development [27, 28]. Most recent studies have shown an association of human preimplantation embryos after ART with incorrect DNA methylation at gDMRs [29,30,31]. Surprisingly, day 3 embryos (76%) and blastocysts (50%) after ART had a high frequency of imprinted methylation errors at three gDMRs .
In this study, the frequencies of ART-BWS and ART-SRS were still high and the frequency of ART-PWS patients increased more than threefold compared to our 2009 survey . Mussa et al. reported that ART entails a tenfold increased risk of BWS . We reconfirmed that ART increases risk of not only BWS but also SRS and PWS. The frequency of paternal uniparental disomy of chromosome 11 (UPD(11)pat) in Sp-BWS increased approximately two times compared with previous data . Sasaki et al. reported that the frequency of hypomethylation in KCNQ1OT1:transcriptional start site (TTS)-DMR in Japanese BWS patients was lower than that in North American/European BWS patients . In addition, we compared the frequency of UPD(11)pat in Sp-BWS patients stratified according to maternal age but there was no correlation between UPD(11)pat and maternal age. Based on the above, the increased frequency of UPD(11)pat in Sp-BWS might be due to race rather than maternal age. Matsubara et al. reported that advanced maternal age affected the onset of PWS through UPD(15)mat, resulting from a disjunction error during meiosis I in oocytes . In stratified analyses of the pathogenesis of PWS according to maternal age, the rate of DNA methylation errors in ART-PWS patients was significantly increased compared with Sp-PWS patients when maternal age was ≤ 37 years, while the rate of UPD(15)mat in ART-PWS patients was significantly increased compared with that in Sp-PWS patients with maternal age ≥ 38 years. These results suggested that the increasing incidence of ART-PWS might be affected not only by maternal age but also DNA methylation altered by ART. The major pathogenesis of AS is deletion or mutation of the UBA3A gene on chromosome 15q11.2-q13 and in total the paternal UPD and imprinting defects, which might be associated with maternal age and ART, are less than 5% of the pathogenesis of AS . Therefore, it may not be possible to observe the increased frequency of ART-AS in this study.
We showed that phenotypic differences between ART patients and Sp patients were largely unreported, while changes to phenotypes might be altered by the frequency and the degree of epimutations caused by ART. Lim et al. reported that ART-BWS patients had a significantly lower frequency of exomphalos and higher risk of non-Wilms’ tumor neoplasia . However, we found no major differences in clinical features between the ART-BWS and Sp-BWS patients. The ART-SRS patients having defects at additional loci and DMVs other than at the domain responsible for that disorder displayed additional minor phenotypes, heart malformation, tumorigenesis, diabetes, and mental retardation. For mental retardation, Zhu et al. reported that children born after ART had a slight delay in motor and cognitive development compared with children born to infertile couples who conceived naturally, and it suggested that ART may increase the risk of mental retardation . The higher frequency of mental retardation in ART-SRS patients compared with Sp-SRS patients may be affected by the risk of ART being added to the risk of SRS itself. Since it is unlikely that methylomes in peripheral blood precisely reflect phenotypes of imprinting disorders, we might not be able to clarify the relationship between phenotypes and epimutations. In addition, it may be that DMVs alone have not reached the threshold for the appearance of clinical phenotypes.
Of the 79 genes with DMVs, 19 candidate genes might be involved in clinical symptoms of SRS. TSPAN1, LIMS2, SCHIP1, HOXA5, ABALON, and AATBC genes might be involved in the growth restriction phenotype of SRS. The SCHIP1 gene is associated with facial size and shape in African children  and Schip1-knockout mice have skeletal, craniofacial, and digestive problems . HOXA5 regulates gene expression, morphogenesis and differentiation, and upregulates the tumor suppressor p53 . ABALON and AATBC are involved in apoptosis, which plays an important role in embryonic development . ELEN-AS1, NPM2, GSTO2, H19, KCNQ1OT1, MEG3, and BANP might be involved in tumorigenesis. H19, KCNQ1OT1, and MEG3 are also known imprinted genes that function as tumor suppressors. The main clinical phenotype of SRS is growth retardation. A tumor is a type of abnormal and excessive growth of cells. These are opposite growth disturbances. Nevertheless, we observed tumorigenesis in ART-SRS patients in our epidemiological study. In addition, at least three papers have reported tumorigenesis in SRS patients [41,42,43]. Therefore, we focused on genes related to tumorigenesis. SLC25A25 and TCF7L2 might be involved in hypoglycemia. SLC25A25 encodes an ATP-Mg/Pi inner mitochondrial membrane solute transporter and Slc25a25-knockout mice are metabolically inefficient . The Wnt signaling pathway regulator TCF7L2 is associated with an increased risk of type 2 diabetes and negative regulators of hepatic gluconeogenesis . RIPOR2 might be involved in defective hearing since it encodes a plasma membrane-associated protein of hair cell stereocilia that is essential for hearing . CMYA5, PRKAG2, and CTF1 might be involved in heart malformation. CMYA5 is associated with left ventricular hypertrophy . PRKAG2 is associated with glycogen-storage cardiomyopathy  and CTF1 induces cardiac myocyte hypertrophy in vitro .
There are some limitations in this study. First, we could not recruit a sufficient number of ART patients with four imprinting disorders to compare with Sp patients for clinical phenotypes and pathogeneses and find any novel appearance for the imprinting disorders in the ART patients. However, there might be differences between the two groups that were not dealt with in our questionnaires. Second, there may be additional strong associations with imprinting disorders and ART because the details about whether or not ART treatments were employed are unknown for half of the patients with imprinting disorders. Third, we cannot provide the reasons for the use of ART for ART-associated patients because we did not inquire about the reason for using ART in the questionnaire. Fourth, we analyzed genome-wide methylation levels in peripheral blood using RRBS technology, which can only cover approximately 10% of genomic CpG cytosines. However, a strength of our study is that the sample group consisted of only Japanese subjects. The Japanese are a relatively homogeneous people geographically and traditionally, because Japan is a small island country and due to the Japanese temperament. Furthermore, genetic differences might affect clinical phenotypes and pathogeneses of imprinting disorders because they are present in the gene-specific DNA methylation levels at birth .
We demonstrated increased frequencies of ART-BWS, ART-SRS, and ART-PWS compared with the live birth rate after ART from 1985 through 2015 in Japan and the characteristics of incomplete and broader DMVs in ART-SRS patients. We reconfirmed the association between ART and three imprinting disorders (BWS, PWS, and SRS). This is perhaps not surprising given the major epigenetic events that take place during early development at a time when the epigenome is most vulnerable. The numbers of children with imprinting disorders born after ART treatment remain limited and the vast majority of children born after such treatments are healthy. On the other hand, it is also well known that epigenetic alterations can increase the risks for various diseases later in life such as cardiovascular disease, type 2 diabetes, and hypertension . Further large-scale, long-term studies will be needed to determine the associations with epimutations and phenotypes of children conceived after ART.
Materials and methods
Nationwide search of four congenital imprinting disorders
The protocol was almost the same as that we used previously and was established by the Research Committee on the Epidemiology of Intractable Diseases. The nationwide survey used a two-stage postal method . First, a survey was conducted to estimate the number of individuals with any of the four imprinting diseases: BWS, SRS, PWS, and AS. In the following stage, a second survey was used to identify the clinico-epidemiological features of these syndromes. We previously reported the questionnaires in detail .
Reduced representation bisulfite sequencing
Genomic DNA was obtained from whole blood using the standard extraction protocol and treated with PureLink RNase A (Thermo Fisher Scientific, Waltham, MA, USA). The DNA was used for reduced representation bisulfite sequencing (RRBS) library construction as previously described with modifications . Briefly, 20 ng of genomic DNA was digested with 20 units of MspI enzymes (NEB, Beverly, MA, USA) in 15 μl reactions with 1x NEBuffer 2 at 37 °C for 16 h and 80 °C for 20 min. The digested products were filled in, and an adenosine was added with 10 units of Klenow Fragment (3′ → 5′ exo−, NEB), 4 μM dGTP, 4 μM dCTP, 40 μM dATP, and 1× NEBuffer 2 to a final volume of 17 μl. The end-repaired/dA-tailing reaction was incubated at 30 °C for 20 min, 37 °C for 20 min, and 75 °C for 20 min. The end-repaired fragments were immediately ligated to 1 μM pre-annealed methylated adaptors using 2000 cohesive end units of T4 DNA Ligase (NEB) in 20 μl reactions with 1× Ligase Reaction Buffer. The ligation reactions were incubated at 16 °C for 16 h and 65 °C for 20 min. The ligated products were separated on 3% NuSieve 3:1 agarose gel (Lonza Japan, Tokyo, Japan) by electrophoresis to isolate 150–350 bp fragments and size-selected fragments were extracted and purified using a MinElute Gel Extraction kit (Qiagen, Valencia, CA, USA) and 20 μl of warmed EB buffer. Then bisulfite conversion was performed using an EpiTect Bisulfite Kit (Qiagen). To incorporate the sample-specific index and amplify the library, PCR was performed with KAPA HiFi HotStart Uracil+ ReadyMix (Kapa Biosystems, Woburn, MA). The PCR amplification was carried out as follows: an initial denaturation at 95 °C for 2 min, 14 cycles of 98 °C for 20 s, 65 °C for 30 s and 72 °C for 30 s, and a final 1-min extension at 72 °C. The RRBS libraries were purified twice using Agencourt AMPure XP beads (Beckman Coulter, Brea, CA, USA), quantified with KAPA Library Quantification Kit (Kapa Biosystems), and sequenced on an HiSeq 2500 platform (Illumina, CA, USA) with 101-bp single-end reads using a TruSeq SR Cluster Kit v3-cBot-HS and TruSeq SBS Kit v3-HS (Illumina). Sequenced reads were processed using an Illumina standard base-calling pipeline (v1.8.2) and the index and adapter sequences were removed. After trimming the first and last four bases, the reads were aligned to Human Genome assembly (GrCh37/hg19) using Bismark (v.0.9.0) with default parameters . Reads that showed < 90% bisulfite conversion were filtered to remove those that resulted from incomplete bisulfite-converted molecules. The methylation level of each cytosine was calculated using the Bismark methylation extractor. We analyzed only CpG cytosines covered with ≥ 5 reads.
Annotations of genomic regions
Annotations of Refseq genes and repeat sequences were downloaded from the UCSC Genome Browser (https://genome.ucsc.edu/). Refseq genes shorter than 300 bp (encoding microRNAs or small nucleolar RNAs) were excluded from our analyses. Promoters were defined as regions 500 bp upstream and downstream from the transcription start sites of the Refseq transcripts. For calculation of the mean methylation levels, we analyzed the genomic regions (promoters, exons, introns, intergenic regions, imprinted gDMRs, CGIs, CGI shores, and CGI shelves) containing ≥ 10 CpG cytosines and SINEs, LINEs, LTRs, DNA repeats, SVAs, and simple repeats containing ≥ 3 CpG cytosines with sufficient coverage for calculation of the methylation levels. Regions and names of the 21 imprinted gDMRs were defined as previously reported .
Statistical analyses were performed using R or JMP Pro 13.1.0 (SAS Institute, Cary, NC, USA). Statistical significance among the clinical characteristics was determined using Fisher’s exact test or the Mann–Whitney U test at p < 0.05. Correlations of all CpG cytosines among individual ART-SRS patients, Sp-SRS patients, and ten spontaneously conceived normal children (control) were calculated using the Pearson correlation coefficient. Statistical significance in the mean methylation levels of various genomic regions among the three groups was determined using Steel’s test at p < 0.05. Using the Mann–Whitney U test with the Benjamini-Hochberg correction, we selected the DNA methylation variations (DMVs) showing absolute methylation changes ≥ 7.5% and a FDR < 0.05. Fisher’s exact test was used to calculate the p values for the number of DMVs. Fisher’s exact test with the Benjamini-Hochberg procedure was used to calculate the p values for enriched GO terms in DAVID analysis .
Assisted reproductive technologies
Database for Annotation, Visualization and Integrated Discovery
DNA Data Bank of Japan
Differentially methylated region
DNA methylation variation
False discovery rate
Germline differentially methylated region
Intracytoplasmic sperm injection
In vitro fertilization
Long interspersed nuclear element
Long terminal repeat
Reduced representation bisulfite sequencing
Short interspersed nuclear element
Transcriptional start site
Paternal uniparental disomy of chromosome 11
Maternal uniparental disomy of chromosome 15
Variable number of tandem repeat
Dyer S, Chambers GM, de Mouzon J, Nygren KG, Zegers-Hochschild F, Mansour R, Ishihara O, Banker M, Adamson GD. International committee for monitoring assisted Reproductive technologies world report: assisted reproductive technology 2008, 2009 and 2010. Hum Reprod. 2016;31(7):1588–609.
Lucifero D, Chaillet JR, Trasler JM. Potential significance of genomic imprinting defects for reproduction and assisted reproductive technology. Hum Reprod Update. 2004;10(1):3–18.
Chiba H, Hiura H, Okae H, Miyauchi N, Sato F, Sato A, Arima T. DNA methylation errors in imprinting disorders and assisted reproductive technology. Pediatr Int. 2013;55(5):542–9.
Matsubara K, Murakami N, Fukami M, Kagami M, Nagai T, Ogata T. Risk assessment of medically assisted reproduction and advanced maternal ages in the development of Prader-Willi syndrome due to UPD(15)mat. Clin Genet. 2016;89(5):614–9.
Cortessis VK, Azadian M, Buxbaum J, Sanogo F, Song AY, Sriprasert I, Wei PC, Yu J, Chung K, Siegmund KD. Comprehensive meta-analysis reveals association between multiple imprinting disorders and conception by assisted reproductive technology. J Assist Reprod Genet. 2018;35(6):943–52.
Mussa A, Molinatto C, Cerrato F, Palumbo O, Carella M, Baldassarre G, Carli D, Peris C, Riccio A, Ferrero GB. Assisted Reproductive Techniques and Risk of Beckwith-Wiedemann Syndrome. Pediatrics. 2017;140(1):e20164311.
Uk A, Collardeau-Frachon S, Scanvion Q, Michon L, Amar E. Assisted reproductive technologies and imprinting disorders: results of a study from a French congenital malformations registry. Eur J Med Genet. 2018;61(9):518–23.
Lidegaard O, Pinborg A, Andersen AN. Imprinting diseases and IVF: Danish national IVF cohort study. Hum Reprod. 2005;20(4):950–4.
Doornbos ME, Maas SM, McDonnell J, Vermeiden JP, Hennekam RC. Infertility, assisted reproduction technologies and imprinting disturbances: a Dutch study. Hum Reprod. 2007;22(9):2476–80.
Butler MG. Genomic imprinting disorders in humans: a mini-review. J Assist Reprod Genet. 2009;26(9–10):477–86.
Hiura H, Okae H, Chiba H, Miyauchi N, Sato F, Sato A, Arima T. Imprinting methylation errors in ART. Reproductive Med Biol. 2014;13(4):193–202.
Uyar A, Seli E. The impact of assisted reproductive technologies on genomic imprinting and imprinting disorders. Curr Opin Obstet Gynecol. 2014;26(3):210–21.
Hiura H, Okae H, Miyauchi N, Sato F, Sato A, Van De Pette M, John RM, Kagami M, Nakai K, Soejima H, et al. Characterization of DNA methylation errors in patients with imprinting disorders conceived by assisted reproduction technologies. Hum Reprod. 2012;27(8):2541–8.
Lim D, Bowdin SC, Tee L, Kirby GA, Blair E, Fryer A, Lam W, Oley C, Cole T, Brueton LA, et al. Clinical and molecular genetic features of Beckwith-Wiedemann syndrome associated with assisted reproductive technologies. Hum Reprod. 2009;24(3):741–7.
Tee L, Lim DH, Dias RP, Baudement MO, Slater AA, Kirby G, Hancocks T, Stewart H, Hardy C, Macdonald F, et al. Epimutation profiling in Beckwith-Wiedemann syndrome: relationship with assisted reproductive technology. Clin Epigenetics. 2013;5(1):23.
Eggermann T, Heilsberg AK, Bens S, Siebert R, Beygo J, Buiting K, Begemann M, Soellner L. Additional molecular findings in 11p15-associated imprinting disorders: an urgent need for multi-locus testing. J Mol Med (Berl). 2014;92(7):769–77.
Fragouli E, Alfarawati S, Spath K, Jaroudi S, Sarasa J, Enciso M, Wells D. The origin and impact of embryonic aneuploidy. Hum Genet. 2013;132(9):1001–13.
Okae H, Chiba H, Hiura H, Hamada H, Sato A, Utsunomiya T, Kikuchi H, Yoshida H, Tanaka A, Suyama M, et al. Genome-wide analysis of DNA methylation dynamics during early human development. PLoS Genet. 2014;10(12):e1004868.
Messerschmidt DM, Knowles BB, Solter D. DNA methylation dynamics during epigenetic reprogramming in the germline and preimplantation embryos. Genes Dev. 2014;28(8):812–28.
Guo H, Zhu P, Yan L, Li R, Hu B, Lian Y, Yan J, Ren X, Lin S, Li J, et al. The DNA methylation landscape of human early embryos. Nature. 2014;511(7511):606–10.
Smith ZD, Chan MM, Humm KC, Karnik R, Mekhoubad S, Regev A, Eggan K, Meissner A. DNA methylation dynamics of the human preimplantation embryo. Nature. 2014;511(7511):611–5.
Zhu P, Guo H, Ren Y, Hou Y, Dong J, Li R, Lian Y, Fan X, Hu B, Gao Y, et al. Single-cell DNA methylome sequencing of human preimplantation embryos. Nat Genet. 2018;50(1):12–9.
Rossignol S, Steunou V, Chalas C, Kerjean A, Rigolet M, Viegas-Pequignot E, Jouannet P, Le Bouc Y, Gicquel C. The epigenetic imprinting defect of patients with Beckwith-Wiedemann syndrome born after assisted reproductive technology is not restricted to the 11p15 region. J Med Genet. 2006;43(12):902–7.
Bliek J, Verde G, Callaway J, Maas SM, De Crescenzo A, Sparago A, Cerrato F, Russo S, Ferraiuolo S, Rinaldi MM, et al. Hypomethylation at multiple maternally methylated imprinted regions including PLAGL1 and GNAS loci in Beckwith-Wiedemann syndrome. Eur J Human Genet. 2009;17(5):611–9.
Young LE, Sinclair KD, Wilmut I. Large offspring syndrome in cattle and sheep. Rev Reprod. 1998;3(3):155–63.
Young LE, Fernandes K, McEvoy TG, Butterwith SC, Gutierrez CG, Carolan C, Broadbent PJ, Robinson JJ, Wilmut I, Sinclair KD. Epigenetic change in IGF2R is associated with fetal overgrowth after sheep embryo culture. Nat Genet. 2001;27(2):153–4.
Doherty AS, Mann MR, Tremblay KD, Bartolomei MS, Schultz RM. Differential effects of culture on imprinted H19 expression in the preimplantation mouse embryo. Biol Reprod. 2000;62(6):1526–35.
Khosla S, Dean W, Brown D, Reik W, Feil R. Culture of preimplantation mouse embryos affects fetal development and the expression of imprinted genes. Biol Reprod. 2001;64(3):918–26.
Chen SL, Shi XY, Zheng HY, Wu FR, Luo C. Aberrant DNA methylation of imprinted H19 gene in human preimplantation embryos. Fertil Steril. 2010;94(6):2356–8. 8 e1
Shi X, Chen S, Zheng H, Wang L, Wu Y. Abnormal DNA methylation of imprinted loci in human preimplantation embryos. Reprod Sci. 2014;21(8):978–83.
White CR, Denomme MM, Tekpetey FR, Feyles V, Power SG, Mann MR. High frequency of imprinted methylation errors in human preimplantation embryos. Sci Rep. 2015;5:17311.
Brioude F, Kalish JM, Mussa A, Foster AC, Bliek J, Ferrero GB, Boonen SE, Cole T, Baker R, Bertoletti M, et al. Expert consensus document: clinical and molecular diagnosis, screening and management of Beckwith-Wiedemann syndrome: an international consensus statement. Nat Rev Endocrinol. 2018;14(4):229–49.
Sasaki K, Soejima H, Higashimoto K, Yatsuki H, Ohashi H, Yakabe S, Joh K, Niikawa N, Mukai T. Japanese and North American/European patients with Beckwith-Wiedemann syndrome have different frequencies of some epigenetic and genetic alterations. Eur J Hum Genet. 2007;15(12):1205–10.
Matsubara K, Murakami N, Nagai T, Ogata T. Maternal age effect on the development of Prader-Willi syndrome resulting from upd(15)mat through meiosis 1 errors. J Hum Genet. 2011;56(8):566–71.
Buiting K, Williams C, Horsthemke B. Angelman syndrome—insights into a rare neurogenetic disorder. Nat Rev Neurol. 2016;12(10):584–93.
Zhu JL, Basso O, Obel C, Hvidtjorn D, Olsen J. Infertility, infertility treatment and psychomotor development: the Danish National Birth Cohort. Paediatr Perinat Epidemiol. 2009;23(2):98–106.
Cole JB, Manyama M, Kimwaga E, Mathayo J, Larson JR, Liberton DK, Lukowiak K, Ferrara TM, Riccardi SL, Li M, et al. Genomewide association study of African children identifies association of SCHIP1 and PDE8A with facial size and shape. PLoS Genet. 2016;12(8):e1006174.
Schmahl J, Raymond CS, Soriano P. PDGF signaling specificity is mediated through multiple immediate early genes. Nat Genet. 2007;39(1):52–60.
Jeannotte L, Gotti F, Landry-Truchon K. Hoxa5: a key player in development and disease. J Dev Biol. 2016;4(2):1–18.
Brill A, Torchinsky A, Carp H, Toder V. The role of apoptosis in normal and abnormal embryonic development. J Assist Reprod Genet. 1999;16(10):512–9.
Draznin MB, Stelling MW, Johanson AJ. Silver-Russell syndrome and craniopharyngioma. J Pediatr. 1980;96(5):887–9.
Weiss GR, Garnick MB. Testicular cancer in a Russell-Silver dwarf. J Urol. 1981;126(6):836–7.
Chitayat D, Friedman JM, Anderson L, Dimmick JE. Hepatocellular carcinoma in a child with familial Russell-Silver syndrome. Am J Med Genet. 1988;31(4):909–14.
Anunciado-Koza RP, Zhang J, Ukropec J, Bajpeyi S, Koza RA, Rogers RC, Cefalu WT, Mynatt RL, Kozak LP. Inactivation of the mitochondrial carrier SLC25A25 (ATP-Mg2+/Pi transporter) reduces physical endurance and metabolic efficiency in mice. J Biol Chem. 2011;286(13):11659–71.
Ip W, Shao W, Chiang YT, Jin T. The Wnt signaling pathway effector TCF7L2 is upregulated by insulin and represses hepatic gluconeogenesis. Am J Phys Endocrinol Metab. 2012;303(9):E1166–76.
Diaz-Horta O, Subasioglu-Uzak A, Grati M, DeSmidt A, Foster J 2nd, Cao L, Bademci G, Tokgoz-Yilmaz S, Duman D, Cengiz FB, et al. FAM65B is a membrane-associated protein of hair cell stereocilia required for hearing. Proc Natl Acad Sci U S A. 2014;111(27):9864–8.
Nakagami H, Kikuchi Y, Katsuya T, Morishita R, Akasaka H, Saitoh S, Rakugi H, Kaneda Y, Shimamoto K, Ogihara T. Gene polymorphism of myospryn (cardiomyopathy-associated 5) is associated with left ventricular wall thickness in patients with hypertension. Hypertens Res. 2007;30(12):1239–46.
Arad M, Maron BJ, Gorham JM, Johnson WH Jr, Saul JP, Perez-Atayde AR, Spirito P, Wright GB, Kanter RJ, Seidman CE, et al. Glycogen storage diseases presenting as hypertrophic cardiomyopathy. N Engl J Med. 2005;352(4):362–72.
Pennica D, King KL, Shaw KJ, Luis E, Rullamas J, Luoh SM, Darbonne WC, Knutzon DS, Yen R, Chien KR, et al. Expression cloning of cardiotrophin 1, a cytokine that induces cardiac myocyte hypertrophy. Proc Natl Acad Sci U S A. 1995;92(4):1142–6.
Adkins RM, Krushkal J, Tylavsky FA, Thomas F. Racial differences in gene-specific DNA methylation levels are present at birth. Birth Defects Res Part A Clin Mol Teratol. 2011;91(8):728–36.
Bateson P, Barker D, Clutton-Brock T, Deb D, D'Udine B, Foley RA, Gluckman P, Godfrey K, Kirkwood T, Lahr MM, et al. Developmental plasticity and human health. Nature. 2004;430(6998):419–21.
Kobayashi N, Okae H, Hiura H, Chiba H, Shirakata Y, Hara K, Tanemura K, Arima T. Genome-scale assessment of age-related DNA methylation changes in mouse spermatozoa. PLoS One. 2016;11(11):e0167127.
Krueger F, Andrews SR. Bismark: a flexible aligner and methylation caller for bisulfite-Seq applications. Bioinformatics. 2011;27(11):1571–2.
Court F, Tayama C, Romanelli V, Martin-Trujillo A, Iglesias-Platas I, Okamura K, Sugahara N, Simon C, Moore H, Harness JV, et al. Genome-wide parent-of-origin DNA methylation analysis reveals the intricacies of human imprinting and suggests a germline methylation-independent mechanism of establishment. Genome Res. 2014;24(4):554–69.
Huang d W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
We would like to thank all the individuals, their families, and medical doctors who participated in this study. We would also like to thank Dr. K. Nakayama, Dr. R. Funayama, Ms. M. Tsuda, Ms. M. Kikuchi, Ms. M. Nakagawa, and Mr. K. Kuroda for technical assistance.
This work was supported by the Japan Agency for Medical Research and Development (AMED: JP17ek0109132h0003 to TA) and Japan Society for the Promotion of Science (KAKENHI: 17H04335 to TA).
Availability of data and materials
All sequence data have been deposited in the DDBJ Japanese Genotype-phenotype Archive under the accession number DRA007187.
Ethics approval and consent to participate
The study was approved by the Ethics Committee of the Tohoku University School of Medicine (2015-1-130, 2015-1-735, 2017-1-880, 2018-1-81). All samples were obtained after receiving written informed consent. All experiments handling human blood cells were conducted in accordance with the principles of the Declaration of Helsinki.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The numbers and age distributions of patients with the four imprinting diseases. (a) BWS. (b) AS. (c) PWS. (d) SRS. The vertical axis shows the number of patients and the horizontal axis shows age. (JPG 1026 kb)
Flowchart showing the molecular testing in four imprinted disorders. (a) BWS. (b) AS. (c) PWS. (d) SRS. The new BWS consensus score  and the Netchine-Harbison clinical scoring system (NH-CSS) were used for the diagnoses of BWS and SRS, respectively. NH-CSS, Netchine-Harbison clinical scoring system; UPD, uniparental disomy; GOM, gain of methylation; LOM, loss of methylation. (JPG 1285 kb)
Clinical characterization of ART-patients. OI/TI, ovulation induction/timed intercourse; IUI, intrauterine insemination; IVF, in vitro fertilization; ICSI, intracytoplasmic sperm injection. N.A. indicates that data was not available. Clinical features were as follows, 1. Macroglossia, 2. Earlobe creases, 3. Umbilical hernia, 4. Hemihypertrophy, 5. Exomphalos, 6. Exophthalmos, 7. Hepatomegaly, 8. Nephromegaly, 9. Ocular hypertelorism, 10. Cryptorchism, 11. Occlusal interference, 12. Mental retardation, 13. Epilepsy, 14. Dysphasia, 15. Dyschromatosis, 16. Ictal laughter, 17. Prognathism, 18. Hyposomnia, 19. Convulsions, 20. Microcephaly, 21. Hypotonia, 22. Feeding difficulties, 23. Almond-shaped eyes, 24. Short stature, 25. Triangular mouse, 26. Acromicria, 27. Bulimia, 28. Obesity, 29. Hypogonadism, 30. Diabetes, 31. Gastrointestinal injury, 32. Cardiac failure, 33. Failure to thrive, 34. Triangular-shaped face, 35. Body asymmetry, 36. Clinodactyly of the fifth fingers, 37. Sweating, 38. Heart malformation. 39. Tumorigenesis. (XLSX 12 kb)
Frequency of different prevalence of UPD(15)mat and DNA methylation error in ART-patients and Sp-patients with PWS according to maternal age. All data were obtained from the questionnaire. DNA methylation error indicates gain of methylation at SNRPN-DMR. UPD(15)mat, maternal uniparental disomy of chromosome 15. (XLSX 11 kb)
Frequency of different pathogeneses in ART-patients and Sp-patients with BWS, AS and SRS stratified according to maternal age. (a) BWS. (b) AS. (c) SRS. The numbers and percentages of patients with chromosomal abnormalities, gene mutations and methylation abnormalities were obtained from a questionnaire. For BWS, UPD and gene indicate paternally uniparental disomy of chromosome 11 and CDKN1C, and methylation errors include both gain of methylation at H19/IGF2 IG-DMR and loss of methylation (LOM) at KCNQ1OT1:TSS-DMR, respectively. For AS, UPD and gene indicate paternally uniparental disomy of chromosome 15 and UBE3A, respectively. For SRS, UPD and methylation error indicate maternally uniparental disomy of chromosome 7 and LOM at H19/IGF2 IG-DMR, respectively. UPD, uniparental disomy; LOM, loss of methylation. (XLSX 11 kb)
Sequencing information. (XLSX 10 kb)
Correlations of the methylation levels of all CpG cytosines covered in all samples. (JPG 832 kb)
Methylation levels of various genomic features in autosomes. Mean methylation levels (%) of CpG cytosines in the promoter, gene body, exon, intron, intergenic region, CGI, CGI shore, CGI shelf, SINE, LINE, LTR, DNA repeat element, SVA and simple repeat. Data are shown as mean ± standard deviation. Since the numbers of ART-SRS and Sp-SRS groups were small, we did not calculate p-values between two groups. SINE, short interspersed nuclear element; LINE, long interspersed nuclear element; LTR, long terminal repeat element; SVA, SINE-VNTR-Alu. (XLSX 11 kb)
Lists of DMVs in various genomic regions. N.A. indicates that data was not available. (XLSX 1005 kb)
DMVs in ART-SRS and Sp-SRS patients classified based on methylation of gametes. Sperm-specific methylated regions were ≥ 80% methylated in sperm and ≤ 20% methylated in oocytes, oocyte-specific methylated regions were ≤ 20% methylated in sperm and ≥ 80% methylated in oocytes, both hypermethylated regions were ≥ 80% methylated in both sperm and oocytes, both hypomethylated regions were ≤ 20% methylated in both sperm and oocytes according to our previously reported data . N.A. indicates that data was not available. (JPG 1269 kb)
About this article
Cite this article
Hattori, H., Hiura, H., Kitamura, A. et al. Association of four imprinting disorders and ART. Clin Epigenet 11, 21 (2019) doi:10.1186/s13148-019-0623-3
- Imprinting disorders
- Assisted reproductive technologies (ART)
- Silver-Russell syndrome (SRS)
- Nationwide epidemiological study
- DNA methylome
- DNA methylation variations (DMVs)