DNA methylation study of fetus genome through a genome-wide analysis

DNA methylation is a crucial epigenetic modification of the genome which is involved in embryonic development, transcription, chromatin structure, X chromosome inactivation, genomic imprinting and chromosome stability. Consistent with these important roles, DNA methylation has been demonstrated to be required for vertebrate early embryogenesis and essential for regulating temporal and spatial expression of genes controlling cell fate and differentiation. Further studies have shown that abnormal DNA methylation is associated with human diseases including the embryonic development diseases. We attempt to study the DNA methylation status of CpG islands in fetus related to fetus growth and development. GeneChip® Human Tiling 2.0R Array set is used for analysis of methylated DNA in a whole-genome wide in 8 pairs amniotic fluid and maternal blood DNA samples. We found 1 fetus hypermethylation DNA markers and 4 fetus hypomethylation DNA markers though a Genome-wide analysis. These DNA markers all found to be associated with the critical genes for fetus growth and development (SH2D3C gene, EML3 gene, TRIM71 gene, HOXA3 gene and HOXA5 gene). These genes can be used as a biomarker for association studying of embryonic development, pathological pregnancy and so on. The present study has provided new and fundamental insights into the roles that DNA methylation has in embryonic development and in the pathological pregnancy.

Background DNA methylation, the most important and common part of epigenetic, regulates gene expression by modified DNA without the changing of the genomic sequence. DNA methylation is a process in which a methyl group is covalently bound to a 5-position of cytosine in the context of a cytosine-guanine dinucleotide (CpG) in the mammalian genome [1]. Methylation of CpG-rich promoters is used by mammals to prevent transcriptional initiation, to ensure the inactivation of X chromosome and the silence of imprinted genes and parasitic DNAs.
Recent studies show that methylation also affects the stability and structure of chromosomes. And the abnormal hypomethylation can cause many human diseases including many kinds of cancer [2,3]. However, the potential roles of methylation in tissue-specific gene expression and in the regulation of CpG-poor promoters are less well established [4]. On the other hand, Xiang H et al. estimate that more than 0.11% of genomic cytosines are methylcytosines, all of which probably occur in CG dinucleotides. CG methylation is substantially enriched in gene bodies and positively correlated with gene expression levels, suggesting it has a positive role in gene transcription [5].
Recently, researchers have conducted extensive research on the mechanism of DNA methylation and its role in embryonic development. DNA methylation has been demonstrated to be required for vertebrate early embryogenesis and can regulate temporal and spatial expression of genes controlling cell fate and differentiation which is a complicated and predetermined developmental program. Moreover, DNA methylation provides necessary direction for the multitude of changes that are required to proceed from a fertilized oocyte to a fully developed adult animal [6,7]. Sperm and eggs carry distinctive epigenetic modifications that are adjusted by reprogramming after fertilization. The paternal genome in a zygote undergoes active DNA demethylation before the first mitosis [8]. During the early phases of mammalian development, DNA methylation is extensively reprogrammed [9]. The maintenance of DNA methylation and normal methylation level is necessary for normal embryonic development and tissue-specific differentiation. The changes of DNA methylation will result in abnormal expression of gene and produce various phenotypes such as pregnancy failure, congenital defects and acquired diseases. For the design of targeted therapies to abnormal DNA methylation, there is important theoretical and clinical significance to studying and understanding the status of DNA methylation in the development of embryo.
Methylated Microarray analysis is efficient technology platform for looking for epigenetic methylation markers. With the completion of many whole-genome sequences, new types of genome-wide experiments are possible. Tiling Arrays offer a physical readout of a genome and can be used as a discovery tool for mapping sites of DNA interaction in chromatin immunoprecipitation (ChIP) experiments, understanding global epigenomic changes like methylation. Currently, most studies use various animal models to study the methylation status during embryonic development. In this study, GeneChip® Human Tiling 2.0R Array set is used to look for methylation markers in amniotic fluid samples in an early development of human embryo. Moreover, clone sequencing is used to test and verify the results of Methylated Microarray analysis. And Real-time quantitative PCR (RTQ-PCR) is used to verify the function of hypermethylation and hypomethylation in the identified genomic regions.

Study samples
Eight women with 18-24 weeks of euploid pregnancies admitted to Prenatal Diagnosis Center of Henan Provincial People' hospital between March 2013 and May 2013 were recruited. 10 ml maternal peripheral blood and 20 ml amniotic fluid were collected from the 8 participants. All samples were frozen at −80°C before analysis. The pregnant information about the participants was shown in Table 1. All the participants were not complicated with any diseases and provided their written informed consent for the collection of the samples and subsequent analysis. The investigation was conducted in accordance with humane and ethical research principles of Henan Provincial People' hospital, China. This study was approved by the Ethics Committee of Henan Provincial People's Hospital, China.

Exclusion of maternal DNA
The fifteen autosomal STR loci and a Amelogenin locus are amplified to determine whether there is maternal DNA in fetus amniotic fluid by using the PowerPlex® 16 system kit in a multiplex amplification reaction system following manufacturer's instructions using 25 μl reactions containing 1.0 μl (0.5-2 ng) genomic DNA, 5.0 μl Power-Plex® 16

Clone sequencing
Clone sequencing strategy has been used to verify the selected methylation markers by microarray analysis method in a whole-genome wide (Sangon Biotech Incorporation, Shanghai, China). Before the clone sequencing was conducted, bisulfite conversion of the sample genomic DNA was done by using EpiTect Bisulfite kit (Qiagen Incorporation, Valencia, CA, Spain) according to the manufacturer's procedure. Nine fetus hypermethylation markers and twelve fetus hypomethylation markers were tested. The primer sequences were shown in Table 2.

Real-time quantitative PCR
Real-time quantitative PCR was used to verify the function of hypermethylation and hypomethylation in the identified genomic regions. Total RNA was prepared from maternal peripheral blood and amniotic fluid samples (400 mg) using TRIzol® Reagent (Life Technologies  Incorporation

Results and discussion
DNA Methylation is a crucial epigenetic modification involved in embryonic development, transcription, chromation structure, X chromosome inactivation, gene imprinting, chromosome stability and occurrence of different kinds of tumors [10][11][12][13][14][15]. In recent years, researchers have conducted extensive studies on the mechanism of DNA methylation and its role in embryonic development using different biological models, such as human and mouse cell lines, gene knockout mice, zebra fish and so on [6,7,[16][17][18]. These studies all show that DNA methylation is important for regulating temporal and spatial expression of genes for controlling cell fate and differentiation. In some other studies, DNA methylation has been demonstrated to be required for vertebrate early embryogenesis [19][20][21]. Here we collect eight pairs of amniotic fluid and maternal blood samples for selecting the hypermethylation and hypomethylation DNA markers. The fetus and mother DNA samples are all amplified by using the PowerPlex® 16 system kit in a multiplex amplification reaction system and there is no maternal DNA in fetus amniotic fluid. DNA methylation analysis of the 8 pairs amniotic fluid and maternal blood DNA samples are conducted by using GeneChip® Human Tiling 2.0R Array set. In this study, from the fetus and mother groups we have got  the values of the DNA methylation levels for thousands of methylation markers respectively. In order to get fetus methylation markers which are significant different from that of mother in DNA methylation level, we have compared the values of the DNA methylation levels of the methylation markers of fetus with that of mother by using MAT calculation (p < 0.0001). Scatter plots are constructed to show the results of the comparison using R Software Packages 3.0.1. Two of the Scatter plots were showed in Figure 1. As shown in Figure 1, red solid triangles which represent the DNA methylation level of fetus in this chromosomal location gathered together, blue open circles which represent the DNA methylation level of mother in this chromosomal location gathered together, and red solid triangle groups are separated from blue open circle groups. If the Scatter plots results are the same as Figure 1, we can determine that there must be significant different in DNA methylation level from that of fetus and mother. According to the above principle we select 59 fetus hypermethylation DNA markers and 56 fetus hypomethylation DNA markers which all show significant differences between fetus and mothers. To further understand the sequence structure and function of the 59 fetus hypermethylation DNA markers and 56 fetus hypomethylation DNA markers, we check their associated genes where these markers locate in (http://www.ncbi.nlm.nih.gov/). We have finally identified 9 fetus hypermethylation DNA markers and 12 fetus hypomethylation DNA markers and their associated genes according to their function (showed in Table 4).
To test and verify the 21 selected methylation DNA markers by microarray analysis method in a wholegenome wide, we used clone sequencing strategy. One of the clone sequencing results was shown in Figure 2 If these 7 genes are really regulated by specific methylation status in the fetus, the mRNA expression of ERG and SH2D3C should be lower in the amniotic fluid, whereas those of EML3, TRIM71, HOXA3, HOXA5, and HOXA11 are higher in the amniotic fluid than in the maternal peripheral blood cells. To verify the function of the 7 genes, we did the real-time quantitative PCR (RTQ-PCR) comparison between amniotic fluid samples and maternal peripheral blood cells on these 7 genes. The results of the RTQ-PCR were consistent with the clone sequencing results. However, Analysis of Variance (ANOVA) showed that only SH2D3C was reduced by 0.47 ± 0.33 folds (p < 0.01) in the amniotic fluid, and EML3, TRIM71, HOXA3, and HOXA5 were enhanced by 1.31 ± 0.72, 2.57 ± 1.37, 13.07 ± 9.45 and 2.38 ± 2.69 folds (p < 0.01) respectively in the amniotic fluid than in the maternal peripheral blood cells. In order to clarify whether the differences between mother and fetus were due to fetus gender differences, the experimental data were grouped according to fetus gender. The Analysis of Variance (ANOVA) results showed that there were no significant differences between male and female fetuses in methylation status. So we conclude that the differences in methylation status between mother and fetus have nothing to do with gender differences in our study. Marker FH7 was associated with SH2D3C gene which encoded an adaptor protein and was a member of cytoplasmic protein family involved in cell migration. Until now, Research on the function of SH2D3C gene was still very rare. Marker FL2 was associated with EML3 gene which was a nuclear microtubule-binding protein required for the correct alignment of chromosomes in metaphase [22]. Marker FL3 was associated with TRIM71 gene which cooperated with microRNAs to repress Cdkn1a expression and promote embryonic stem cell proliferation [23,24]. HOX genes called homeobox genes encoded the class of transcription factors were found in clusters A, B, C, and D on four separate chromosomes. The expression of these genes was spatially and temporally regulated during embryonic development. Marker FL1 was associated with HOXA3 gene which was a member of HOX gene family and was part of the A cluster on chromosome 7. It encoded a DNAbinding transcription factor which may regulate gene expression, morphogenesis and differentiation. Marker FL7 was associated with HOXA5 gene which also encoded the class of transcription factors. As HOXA3, the expression of this gene was spatially and temporally regulated during embryonic development. Methylation of this gene may result in the loss of its expression. Since the encoded protein up-regulates the tumor suppressor p53, this protein may play an important role in tumorigenesis. The previous studies showed HOX genes were involved in stem cell differentiation [25]. High-density association study of the HOX gene for volumetric BMD at the femoral neck and lumbar spine among older men was conducted and proved that HOXA gene regions were associated with both femoral neck and lumbar spine BMD [26]. However, HOX genes has not been given enough attention that they can regulate gene expression, morphogenesis, and differentiation and were involved in the regulation of uterine development and required for female fertility.

Conclusions
In this study, we found 1 fetus hypermethylation DNA markers and 4 fetus hypomethylation DNA markers which all associated with the critical genes for fetus growth and development. Our results may provide valuable information for future research of fetus growth and development, and even can give some inspiration of the studies of the pathological pregnancy, such as Unexplained Recurrent Spontaneous Abortion (URSA) and early pregnant failure which generally occur at about 20 weeks in the embryonic development [27].