This article has Open Peer Review reports available.
DNA methylation study of fetus genome through a genome-wide analysis
- Hong-Dan Wang†1, 2,
- Qiao-Fang Hou†1,
- Qian-Nan Guo1,
- Tao Li1,
- Dong Wu1,
- Xian-Ping Zhang3,
- Yan Chu1,
- Miao He1,
- Hai Xiao1,
- Liang-Jie Guo1,
- Ke Yang1,
- Shi-Xiu Liao1Email author and
- Bo-Feng Zhu2Email author
© Wang et al.; licensee BioMed Central Ltd. 2014
Received: 25 December 2013
Accepted: 10 April 2014
Published: 15 April 2014
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.
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 . 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 . 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 .
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 . During the early phases of mammalian development, DNA methylation is extensively reprogrammed . 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.
General information about the participants
Maternal age (years)
19 week + 3 day
19 week + 6 day
21 week + 5 day
20 week + 5 day
24 week + 2 day
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 PowerPlex® 16 5 × Master Mix, 5.0 μl PowerPlex® 16 5 × Primer Pair Mix and 14.0 μl ddH2O. Thermal cycling was performed using the GeneAmp® PCR System 9700 (Applied Biosystems, Foster City, CA, USA). Temperature cycling conditions for PCR reactions are as follows: denaturation for 11 min at 95°C, and then 1 min at 96°C, followed by 10 cycles for 30 s at 94°C, 30 s at 60°C, 45 s at 70°C, then 20 cycles for 30 s at 90°C, 30 s at 60°C, 45 s at 70°C and a final elongation step at 60°C for 30 min.
The PCR products have been separated and detected by capillary electrophoresis on an ABI 3130 Genetic Analyzer (Applied Biosystems, Foster City, CA, USA). 2 μl PCR product or Allelic Ladder has been mixed with 15 μl Hi-Di formamide and 0.5 μL internal lane standard ILS 600 (Promega, WI, USA). The loading mixture has been denatured at 95°C for 3 min, followed by chilling on ice for 3 min immediately. STR alleles have been analyzed by comparison with kit allelic ladders using the GeneMapper ID 3.2 software (Applied Biosystems, Foster City, CA, USA). Control DNA from cell line 9947A (Promega, Madison, WI, USA) has been genotyped as standard reference in all experiments.
Methylated DNA capture
Methylamp™ Methylated DNA Capture kit (Epigentek Group Incorporation, Farmingdale, New York, USA) has been used to enrich and capture methylated DNA fragments from blood and amniotic fluid samples for use in gene-specific DNA methylation studies on a genome wide scale. According to the manufacturer’s protocol of the kit, Methylated DNA Capture include the following steps: i) DNA isolation, ii) breaking DNA as fragments by sonicator: sonication for 10 s, and pause 30 s, repeat it 15 times at 100 W, check the DNA size by running 1.2% agarose gel, iii) incubating sonicated/digested DNA at 98°C for 3 min and immediately put on ice, iv) bisulfite conversion and methylated DNA immunoprecipitation, v) release of mDNA from antibody complex, vi) capturing of mDNA, vii) elution of mDNA. Following methylated DNA immunoprecipitation (meDIP) processed, the methlation of the DNA can be analyzed using meDIP-ChIP.
DNA methylation analysis
GeneChip® Human Tiling 2.0R Array set (Affymetrix Incorporation, Santa Clara, CA, USA) has been used for analysis of methylated DNA in a whole-genome wide. GeneChip® Fluidics Station 400 and GeneChip® Scanner 3000 7G (Affymetrix Incorporation, Santa Clara, CA, USA) are needed for analysis. According to the Affymetrix chip-on-chip procedures we first amplified immunoprecipitated DNA targets using the following cycle conditions: 14 cycles for 30 s at 95°C, 30 s at 45°C, 30 s at 55°C, 1 min at 72°C, then 14 cycles for 30 s at 95°C, 30 s at 45°C, 30 s at 55°C (for every subsequent cycle add 5 second), 1 min at 72°C, and 4°C for ever. Secondly, we purified the PCR products using MiniElute PCR purification kit (Qiagen Incorporation, Valencia, CA, Spain) and got fragmentation mix, then put up fragmentation mix for Tiling array, following this step we used GeneChip® WT Double-Stranded DNA Terminal Labeling Kit (Affymetrix Incorporation, Santa Clara, CA, USA) to label fragmented Double-Stranded DNA, and used GeneChip® Hybridization, Wash, and Stain Kit (Affymetrix Incorporation, Santa Clara, CA, USA) to hybridize, wash and stain the labeled target DNA according to the procedure of the kit. Finally, we scan the barcode on the chip using appropriate protocol.
The primers for clone sequencing
Amplicon size (bp)
Real-time quantitative PCR
The primers for real-time quantitative PCR
Amplicon size (bp)
Affymetrix® Tiling Analysis Software (TAS) was used for analysis and quality control of the GeneChip Tiling Arrays. TAS analyzed feature intensity data stored in GeneChip® Operating Software (GCOS), output cel files and produced signals and p-values for each genomic position interrogated. Affymetrix Power Tools Software 1.15.1 was used to convert cel files to text files. The methylation regions of fetus group and mother control group were compared by using Modle-based Analysis of Tiling-array (MAT) calculation (p < 0.0001). CisGenome ver 2.0 Software (http://www.biostat.jhsph.edu/~hji/cisgenome/) was used for annotating the gene names which were close to methylation regions. Scatter plots were constructed utilizing R Software Packages 3.0.1. Two sets of results (the fetus and the mohter) of clone sequencing were analyzed by Quantification tool for Methylation Analysis (QUMA) and compared by statistical software SPSS Version 13.0 using Analysis of Variance (ANOVA) method (p < 0.01). The statistical analysis of the RTQ-PCR results was done using ΔΔCt methods. Analysis of Variance (ANOVA) was also used for comparison the two sets of results of RTQ-PCR (p < 0.01).
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–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–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–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.
The information of the 9 fetus hypermethylation DNA markers (FH) and 12 fetus hypomethylation DNA markers (FL)
Mitochondrial DNA D-loop
Mitochondrial DNA D-loop
Free fatty acid receptor 2 (FFAR2)
Erythroblastosis virus E26 oncogene onmcogene homolog (ERG)
Carbohydrate (chondroitin 4) Sulfotransferase 12 (CHST12)
Intermediate filament family orphan 1 (IFFO1)
SH2 domain containing 3C (SH2D3C)
Keratin associated protein 16-1 (KRTAP16-1)
Histocompatibility (minor) HA-1 (HMHA1)
Homeobox A3 (HOXA3)
Echinoderm microtubule associated protein like 3 (EML3)
Tripartite motif containing 71, E3 ubiquitin protein ligase (TRIM71)
Homo sapiens homeobox A11 (HOXA11)
Homo sapiens glucosaminyl (N-acetyl) transferase 2, I-branching enzyme (GCNT2)
Homo sapiens homeobox A2 (HOXA2)
Homeobox A5 (HOXA5)
Nuclear factor I/X (NFIX) (CCAAT-binding transcription factor)
Homo sapiens trichorhinophalangeal syndrome I (TRPS1)
Elongation factor, RNA polymerase II, 2 (ELL2)
Phosphodiesterase 1C, calmodulin-dependent (PDE1C)
Phosphoprotein associated with glycosphingolipid microdomains 1 (PAG1)
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 . 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 DNA-binding 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 . 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 . 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.
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 .
This project was supported by the National Natural Science Foundation of China (NSFC, No. 81170581) and the National Natural Science Foundation of China (NSFC, No. 81373248).
- Geiman TM, Muegge K: DNA methylation in early development. Mol Reprod Dev. 2010, 77 (2): 105-113.PubMedGoogle Scholar
- Takahashi T, Matsuda Y, Yamashita S, Hattori N, Kushima R, Lee YC, Igaki H, Tachimori Y, Nagino M, Ushijima T: Estimation of the fraction of cancer cells in a tumor DNA sample using DNA methylation. PLoS One. 2013, 8 (12): e82302-10.1371/journal.pone.0082302.View ArticlePubMedPubMed CentralGoogle Scholar
- Pinto R, De Summa S, Pilato B, Tommasi S: DNA Methylation and miRNAs regulation in hereditary breast cancer. Curr Mol Med. 2013, In pressGoogle Scholar
- Jones PA, Takai D: The role of DNA methylation in mammalian epigenetics. Science. 2001, 293 (5532): 1068-1070. 10.1126/science.1063852.View ArticlePubMedGoogle Scholar
- Xiang H, Zhu J, Chen Q, Dai F, Li X, Li M, Zhang H, Zhang G, Li D, Dong Y, Zhao L, Lin Y, Cheng D, Yu J, Sun J, Zhou X, Ma K, He Y, Zhao Y, Guo S, Ye M, Guo G, Li Y, Li R, Zhang X, Ma L, Kristiansen K, Guo Q, Jiang J, Beck S, et al: Single base-resolution methylome of the silkworm reveals a sparse epigenomic map. Nat Biotechnol. 2010, 28 (5): 516-520. 10.1038/nbt.1626.View ArticlePubMedGoogle Scholar
- Sen GL, Reuter JA, Webster DE, Zhu L, Khavari PA: DNMT1 maintains progenitor function in self-renewing somatic tissue. Nature. 2010, 463 (7280): 563-567. 10.1038/nature08683.View ArticlePubMedPubMed CentralGoogle Scholar
- Trowbridge JJ, Snow JW, Kim J, Orkin SH: DNA methyltransferase 1 is essential for and uniquely regulates hematopoietic stem and progenitor cells. Cell Stem Cell. 2009, 5 (4): 442-449. 10.1016/j.stem.2009.08.016.View ArticlePubMedPubMed CentralGoogle Scholar
- Gu TP, Guo F, Yang H, Wu HP, Xu GF, Liu W, Xie ZG, Shi L, He X, Jin SG, Iqbal K, Shi YG, Deng Z, Szabó PE, Pfeifer GP, Li J, Xu GL: The role of Tet3 DNA dioxygenase in epigenetic reprogramming by oocytes. Nature. 2011, 477 (7366): 606-610. 10.1038/nature10443.View ArticlePubMedGoogle Scholar
- Borgel J, Guibert S, Li Y, Chiba H, Schübeler D, Sasaki H, Forné T, Weber M: Targets and dynamics of promoter DNA methylation during early mouse development. Nat Genet. 2010, 42 (12): 1093-1100. 10.1038/ng.708.View ArticlePubMedGoogle Scholar
- Zhou F, Tao G, Chen X, Xie W, Liu M, Cao X: Methylation of OPCML promoter in ovarian cancer tissues predicts poor patient survival. Clin Chem Lab Med. 2013, 11: 1-8.Google Scholar
- Suzuki M, Shigematsu H, Shames DS, Sunaga N, Takahashi T, Shivapurkar N, Iizasa T, Frenkel EP, Minna JD, Fujisawa T, Gazdar AF: DNA methylation-associated inactivation of TGFβ-related genes, DRM/Gremlin, RUNX3, and HPP1 in human cancers. Br J Cancer. 2013, 109 (12): 3132-10.1038/bjc.2013.776.View ArticlePubMed CentralGoogle Scholar
- Verma M, Rogers S, Divi RL, Scully SD, Nelson S, Su LJ, Ross S, Pilch S, Winn DM, Khoury MJ: Epigenetic research in cancer epidemiology: trends, opportunities, and challenges. Cancer Epidemiol Biomarkers Prev. 2013, In pressGoogle Scholar
- Wang X, Douglas KC, Vandeberg JL, Clark AG, Samollow PB: Chromosome-wide profiling of X-chromosome inactivation and epigenetic states in fetal brain and placenta of the opossum, Monodelphis domestica. Genome Res. 2013, In pressGoogle Scholar
- Bala Tannan N, Brahmachary M, Garg P, Borel C, Alnefaie R, Watson CT, Thomas NS, Sharp AJ: DNA methylation profiling in X;autosome translocations supports a role for L1 repeats in the spread of X chromosome inactivation. Hum Mol Genet. 2013, In pressGoogle Scholar
- Girardot M, Feil R, Llères D: Epigenetic deregulation of genomic imprinting in humans: causal mechanisms and clinical implications. Epigenomics. 2013, 5 (6): 715-728. 10.2217/epi.13.66.View ArticlePubMedGoogle Scholar
- Hodges E, Molaro A, Dos Santos CO, Thekkat P, Song Q, Uren PJ, Park J, Butler J, Rafii S, McCombie WR, Smith AD, Hannon GJ: Directional DNA methylation changes and complex intermediate states accompany lineage specificity in the adult hematopoietic compartment. Mol Cell. 2011, 44 (1): 17-28. 10.1016/j.molcel.2011.08.026.View ArticlePubMedPubMed CentralGoogle Scholar
- Fang X, Corrales J, Thornton C, Scheffler BE, Willett KL: Global and gene specific DNA methylation changes during zebrafish development. Comp Biochem Physiol B Biochem Mol Biol. 2013, 166 (1): 99-108. 10.1016/j.cbpb.2013.07.007.View ArticlePubMedGoogle Scholar
- Rai K, Nadauld LD, Chidester S, Manos EJ, James SR, Karpf AR, Cairns BR, Jones DA: Zebra fish Dnmt1 and Suv39h1 regulate organ-specific terminal differentiation during development. Mol Cell Biol. 2006, 26 (19): 7077-7085. 10.1128/MCB.00312-06.View ArticlePubMedPubMed CentralGoogle Scholar
- Bulut HE, Ozdemir O, Başimoglu-Koca Y, Korkmaz M, Atalay A: Effects of a DNA demethylating agent–5-azacytidine–on testicular morphology during mouse embryo development. Okajimas Folia Anat Jpn. 1999, 76 (1): 47-53. 10.2535/ofaj1936.76.1_47.View ArticlePubMedGoogle Scholar
- Snyder MD, Pryor JH, Peoples MD, Williamson GL, Golding MC, Westhusin ME, Long CR: 122 suppression of epigenetic modifiers alters the bovine embryonic developmental program during in vitro culture. Reprod Fertil Dev. 2013, 26 (1): 175.View ArticleGoogle Scholar
- Kraushaar DC, Zhao K: The epigenomics of embryonic stem cell differentiation. Int J Biol Sci. 2013, 9 (10): 1134-1144. 10.7150/ijbs.7998.View ArticlePubMedPubMed CentralGoogle Scholar
- Tegha-Dunghu J, Neumann B, Reber S, Krause R, Erfle H, Walter T, Held M, Rogers P, Hupfeld K, Ruppert T, Ellenberg J, Gruss OJ: EML3 is a nuclear microtubule-binding protein required for the correct alignment of chromosomes in metaphase. J Cell Sci. 2008, 121 (Pt 10): 1718-1726.View ArticlePubMedGoogle Scholar
- Lin YC, Hsieh LC, Kuo MW, Yu J, Kuo HH, Lo WL, Lin RJ, Yu AL, Li WH: Human TRIM71 and its nematode homologue are targets of let-7 microRNA and its zebrafish orthologue is essential for development. Mol Biol Evol. 2007, 24 (11): 2525-2534. 10.1093/molbev/msm195.View ArticlePubMedGoogle Scholar
- Chang HM, Martinez NJ, Thornton JE, Hagan JP, Nguyen KD, Gregory RI: Trim71 cooperates with microRNAs to repress Cdkn1a expression and promote embryonic stem cell proliferation. Nat Commun. 2012, 3: 923.View ArticlePubMedPubMed CentralGoogle Scholar
- Chung N, Jee BK, Chae SW, Jeon YW, Lee KH, Rha HK: HOX gene analysis of endothelial cell differentiation in human bone marrow-derived mesenchymal stem cells. Mol Biol Rep. 2009, 36 (2): 227-235. 10.1007/s11033-007-9171-6.View ArticlePubMedGoogle Scholar
- Yerges LM, Klei L, Cauley JA, Roeder K, Kammerer CM, Moffett SP, Ensrud KE, Nestlerode CS, Marshall LM, Hoffman AR, Lewis C, Lang TF, Barrett-Connor E, Ferrell RE, Orwoll ES, Zmuda JM, MrOS Research Group: High-density association study of 383 candidate genes for volumetric BMD at the femoral neck and lumbar spine among older men. J Bone Miner Res. 2009, 24 (12): 2039-2049. 10.1359/jbmr.090524.View ArticlePubMedPubMed CentralGoogle Scholar
- Lambrot R, Xu C, Saint-Phar S, Chountalos G, Cohen T, Paquet M, Suderman M, Hallett M, Kimmins S: Low paternal dietary folate alters the mouse sperm epigenome and is associated with negative pregnancy outcomes. Nat Commun. 2013, 4: 2889.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/7/18/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.