#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

DNA Methylation Landscapes of Human Fetal Development


Methylation of DNA is a key epigenetic mark. Adult tissues have highly distinct genome-wide DNA methylation signatures. How these signatures arise during human fetal development is largely unknown. Here, we studied DNA methylation profiles of four tissues (amnion, muscle, adrenal, pancreas) during first and second trimester of human fetal development. Already in the first trimester, a tissue-specific signature was found in each of the tissues. However, during the first and second trimester, a substantial number of genomic regions were found to gain and lose DNA methylation. Genomic regions that gained methylation were associated with the shut-down of developmental processes, while genomic regions that lose methylation were associated with the activation of tissue-specific functions. These findings on the DNA methylation landscape of human fetal development are important as they provide insight into regulatory elements that guide tissue specification and lead to organ functionality.


Published in the journal: . PLoS Genet 11(10): e32767. doi:10.1371/journal.pgen.1005583
Category: Research Article
doi: https://doi.org/10.1371/journal.pgen.1005583

Summary

Methylation of DNA is a key epigenetic mark. Adult tissues have highly distinct genome-wide DNA methylation signatures. How these signatures arise during human fetal development is largely unknown. Here, we studied DNA methylation profiles of four tissues (amnion, muscle, adrenal, pancreas) during first and second trimester of human fetal development. Already in the first trimester, a tissue-specific signature was found in each of the tissues. However, during the first and second trimester, a substantial number of genomic regions were found to gain and lose DNA methylation. Genomic regions that gained methylation were associated with the shut-down of developmental processes, while genomic regions that lose methylation were associated with the activation of tissue-specific functions. These findings on the DNA methylation landscape of human fetal development are important as they provide insight into regulatory elements that guide tissue specification and lead to organ functionality.

Introduction

Methylation of CpG dinucleotides in the mammalian genome is a key epigenetic mark. Adult tissues have highly distinct genome-wide DNA methylation signatures consistent with the regulation of cell differentiation by epigenetic mechanisms [13]. Differences in DNA methylation between tissues have been shown to mark differences between germ layers [4], preferentially at regions with low CpG content [2,5,6], at enhancers [4] and alternative promoters [7,8].

Multiple studies have reported on the reprogramming of the human methylome during preimplantation embryo development [911]. In line with previous data on mice [12], in humans DNA methylation is largely erased after conception, the paternal genome being actively and the maternal genome passively demethylated, to become remethylated with the implantation of the embryo [9,10,13,14]. However, systematic and detailed reports on DNA methylation dynamics during human fetal development remain scarce [15], while such data is key to understand how epigenetic mechanisms drive tissue specification and organ functionality. Current views of fetal DNA methylation dynamics are largely extrapolated from studies on the differentiation of human and mouse cells in vitro [7,1521], and the comparison of differentiated tissues to human induced pluripotent stem cells and human embryonic stem cell lines [15]. An exception is fetal brain development in humans, for which recently reported in vivo data showed significant DNA methylation remodelling [15,22].

Recent developments in technology for interrogating genome-wide DNA methylation at single-nucleotide resolution [23] and detailed functional annotation of the human genome [24,25] provide an opportunity to chart DNA methylation during development and assign biological roles to the regions involved. Taking advantage of these developments, we report on DNA methylation dynamics during human fetal development of one extraembryonic tissue and three organs relevant for complex human diseases. This organ-specific catalogue of DNA methylation during development provides fundamental insights into processes guiding human development, but also into the biological function of non-coding regions, which are emerging as important from genome-wide association studies (GWASs) of complex diseases [26]. In addition, this catalogue may serve as a reference for studies on the role of epigenetic mechanisms in the association between an adverse prenatal environment and adulthood disease [27] since DNA methylation marks may have an heightened sensitivity for environmental perturbations during remodelling [28].

Results

Fetal DNA methylation reflects tissue origin and developmental age

To study DNA methylation dynamics in human fetal development, amnion, muscle, adrenal and pancreas samples of 11 fetuses were obtained at 9, 18 and 22 weeks of gestation (W9, W18 and W22; S1A Fig). Genome-wide DNA methylation was investigated with the Illumina 450k array resulting in data on 452,490 CpG sites after quality control [29] (S1B–S1E Fig). The study included three biological replicates per tissue and time point, except for W22 amnion (n = 2) and W22 pancreas (n = 2) (S2A Fig).

We first assessed differences in overall DNA methylation patterns between time points and tissues using hierarchical clustering based on Euclidean distance (Fig 1A) and multidimensional scaling (MDS) (Figs 1B and S2B). DNA methylation patterns clearly differentiated the four tissue types studied (Figs 1A, 1B and S2B). The amnion, representing an extraembryonic tissue, clustered separately from the three embryonic tissues (Fig 1A). Within the embryonic cluster, all W9 tissues (representative of the first trimester) clustered together, whereas W18-W22 tissues (representative of the second trimester) were present towards the edges of the MDS plot (Figs 1B and S2B). Despite the distinct differences in DNA methylation patterns from the first to second trimester, the total number of hypo-, intermediately and hypermethylated CpGs remained constant across time points and tissues (including the extraembryonic amnion) on both autosomes and, in females, the X chromosome (S2C and S2D Fig). This suggests that the observed differences in the MDS plot were not driven by changes in average levels of DNA methylation, but rather due to tissue- and time-specific changes in DNA methylation.

Fig. 1. Tissue- and time-specific DNA methylation patterns during human fetal development.
Tissue- and time-specific DNA methylation patterns during human fetal development.
(A) Clustering based on Euclidean distance. (B) Multidimensional scaling based on Euclidean distance of the four tissues. (C) Median DNA methylation for each of the four tissues over time with a combined genic and CGI-centric annotation. CGI, CpG island.

To validate our findings, we integrated our data with three previously published Illumina 450k datasets on 10 human fetal tissues [15,22,30]. Hierarchical clustering of all data together (n = 117) confirmed the presence of distinct tissue- and time-specific DNA methylation patterns in fetal tissues (S2E Fig).

The characteristics and biological function of DNA methylation depend on the local CpG content and position relative to genes [31]. We mapped CpG sites (CpGs) to CpG islands (CGIs as defined in the UCSC genome browser; 138,919 CpGs), their shores (±2 kb of CGIs; 103,453 CpGs) and shelves (±2 kb of shores; 42,227 CpGs), and remaining CpG-poor non-CGI regions (157,560 CpGs), and to regions relative to gene locations including distal promoters (-10 kb–-1.5 kb; 21,101 CpGs), proximal promoters (-1.5 kb–+0.5 kb; 171,077 CpGs), gene bodies (+0.5 kb–3’ untranslated region (UTR); 175,062 CpGs), downstream regions (3’ UTR–+5 kb; 8,563 CpGs) and remaining intergenic regions (66,356 CpGs; Fig 1C). CpGs were commonly hypomethylated in CGIs, intermediately methylated in shores and hypermethylated in both shelves and non-CGI regions (Fig 1C). These patterns differed by genic position, e.g. CGI methylation was lowest in proximal promoters and highest in gene bodies. Annotation-specific methylation differences were found between W9 and W22, as CpGs in CGIs and shores tended to increase (e.g. gene body CGIs and distal promoter shores), whereas CpGs in non-CGI regions decreased in methylation (e.g. non-CGI promoters). For a subset of annotations, the amnion showed a slightly different DNA methylation patterns than for embryonic tissues, e.g. for gene body CGIs (Fig 1C). Taken together, our data imply that DNA methylation is highly dynamic during fetal development without affecting the average level of DNA methylation.

Hypomethylation discriminates tissues independent of developmental age

It has been shown that each adult tissue is defined by tissue-specific DNA hypomethylation [15,32,33]. Since the four fetal tissues analysed showed a clear DNA methylation signature that corresponded to separated clusters (Fig 1A and 1B), we investigated whether combinations of tissue-specific DNA hypomethylated CpGs were present irrespective of its developmental stage. To do this, we identified CpGs that were relatively hypomethylated (defined as a DNA methylation difference of > 20%) in each tissue compared to all others throughout the three time points of fetal development investigated. The analysis showed indeed that, independently of the developmental age, each tissue showed a cluster of tissue-specific hypomethylated CpGs (Fig 2A). The early lineage segregation of the amnion was further confirmed by the comparatively large number of CpGs (3,536 CpGs) that were exclusively hypomethylated across amniotic samples. In contrast, the embryonic tissues contained much fewer tissue-specific hypomethylated CpGs (muscle 756 CpGs; adrenal 140 CpGs; pancreas 220 CpGs) reflecting their common origin of the epiblast, that gives rise to all embryonic tissues. Genes mapping (i.e. the nearest gene locus) to the specific hypomethylated CpGs per tissue regardless of the time point (amnion 2372, muscle 548, adrenal 120, pancreas 175) were enriched for biological processes that included GO terms characteristic of amnion, muscle and pancreas development and function (S1 Table).

Fig. 2. Sets of hypomethylated CpG sites are tissue-specific.
Sets of hypomethylated CpG sites are tissue-specific.
(A) Heatmap representing hypomethylated CpGs per tissue, defined by a beta difference of ≥ 0.2 of the studied tissue compared to the other tissues. (B) Combined genic and CGI-centric annotation for the hypomethylated CpGs per tissue represented as the odds ratio (see S3A Fig for odds ratios). CGI, CpG island; DP, distal promoter; DS, downstream; GB, gene body; IG, intergenic; NC, non-CGI; PP, proximal promoter; SHE, shelves; SHO, shores. (C) Enrichment of hypomethylated CpGs in the chromatin state segmentation states for the matching tissues (amnion, fetal muscle, fetal adrenal and adult pancreatic islets; see S3B Fig for odds ratios).

When annotated to genic and CGI-related location, it became evident that tissue-specific hypomethylation was enriched for non-CGI regions (P < 0.0001) and highly depleted at CGIs, in particular when mapping to proximal promoters (P < 0.0001; Figs 2B and S3A). To gain further insight in the biological role of genomic regions displaying tissue-specific hypomethylation, we used chromatin state segmentations for fetal muscle, fetal adrenal, amnion and adult pancreatic islets generated by the Epigenomics Roadmap [25]. Tissue-specific hypomethylation was strongly enriched at enhancers (P < 0.001; Figs 2C and S3B). The functional relevance of those tissue-specific hypomethylated CpGs was further validated by comparison to additional fetal samples [15] and adult somatic tissues [6] from available external datasets (S3C Fig). Intriguingly, a high degree of similarity between tissues sharing the same origin, even into adulthood, was observed (S3C Fig).

We next investigated whether tissue-specific hypomethylated CpGs clustered into Hypomethylated Regions (tHRs, defined as 3 consecutive hypomethylated CpGs within 1kb of each other) [6]. This was the case for amnion, muscle and pancreas (Tables 1 and S2). tHRs comprise robust development-independent epigenetic markers as exemplified by the mapping of pancreatic tHRs to proximal promoters of the nearest genes ACY3, HNF1A, and HNF4A (Table 1 and S3D Fig), genes with a key role in pancreas development [34,35]. Muscle tHRs (S2 Table) mapped to distal elements of transcription factors involved in muscle development (NFATC1) and somitogenesis (UNCX) [36,37] (S3D Fig and Tables 1 and S2). Importantly, we could confirm the tHRs identified with the relatively sparse Illumina 450k array with fetal and adult muscle whole-genome bisulfite sequencing (WGBS) data (S3E Fig) [25]. These data indicate that it is feasible to use combinations of tHRs as tissue-specific and development-independent barcodes.

Tab. 1. Table with the numbers of tHRs and three representative genes per tissue associated with tHRs.
Table with the numbers of tHRs and three representative genes per tissue associated with tHRs.
SM, skeletal muscle; TF, transcription factor.

Distinct roles for gain and loss of DNA methylation in fetal development

We provide evidence for large-scale DNA methylation dynamics between W9 and W22 (DNA methylation difference > 20%) that affected 11.5% of evaluated CpGs (52,134/452,490). Approximately equal numbers of CpGs showed a gain of methylation (GOM) (26,555 CpGs; amnion 5,988; muscle 7,631; adrenal 13,997; pancreas 8,620) and a loss of methylation (LOM) (25,579 CpGs; amnion 10,811; muscle 11,925; adrenal 4,476; pancreas 3,286) (Fig 3A). DNA methylation remodelling occurred predominantly between W9 and W18 and not between W18 and adulthood (Figs 3B and S4A). Intriguingly, the integration and re-analysis of external DNA methylation data of fetal adrenal, brain and liver [6,15,22,30,38] revealed a striking confirmation of DNA methylation dynamics during fetal development. Furthermore, for all embryonic tissues, the DNA methylation levels at W22 were similar to those found in the adult counterpart (Figs 3B and S4A), suggesting that the extent of changes after W22 are limited for these CpGs.

Fig. 3. Gain and loss of DNA methylation during human fetal development.
Gain and loss of DNA methylation during human fetal development.
(A) Heatmap of CpGs with a gain and a loss, respectively, of methylation over time. Gain and loss of methylation was defined as a difference of beta ≥ 0.2 between W9 and W22, and W18 in between. (B) Mean methylation of CpGs with a gain or loss of DNA methylation for fetal tissues and their adult counterpart. (C) Combined genic and CGI-centric annotation for CpGs with a gain or a loss of methylation represented as the odds ratio (see S4B Fig for odds ratios). CGI, CpG island; DP, distal promoter; DS, downstream; GB, gene body; IG, intergenic; NC, non-CGI; PP, proximal promoter; SHE, shelves; SHO, shores. (D) Enrichment of dynamically methylated CpGs in the chromatin segmentation states for the matching tissues (fetal muscle, fetal adrenal, amnion and adult pancreatic islets; see S4C Fig for odds ratios).

GOM CpGs did not show tissue-specific patterns (Fig 3A) and, in line with this observation, often corresponded to genes involved in generic developmental and cellular processes, including embryonic morphogenesis and regulation of transcription (S3 Table). In contrast, the LOM CpGs were highly tissue-specific (Fig 3A) and mapped to genes involved in tissue-specific processes that matched the organ in which the LOM CpGs were identified (S3 Table). CpGs that lost methylation in the amnion mapped, amongst others, to genes that were associated with the regulation of apoptosis and cytoskeleton organization; in the muscle to genes associated with cytoskeleton organization and muscle system processes; in the adrenal to genes associated with regulation of macromolecule metabolism (S3 Table). In the pancreas no significant enrichments were found.

GOM and LOM CpGs differed in their genomic annotation. While GOM CpGs were generally enriched in CGIs and CGI-shores, LOM CpGs were enriched for CGI-shelves and non-CGI regions (Figs 3C and S4B). LOM- and GOM-specific enrichments were also observed for Epigenomics Roadmap chromatin state segmentations. LOM CpGs were strongly enriched for (genic) enhancers and transcribed regions, whereas GOM CpGs were enriched for bivalent and repressed regions and only modestly at enhancers (Figs 3D and S4C). The results underscore the relevance of DNA methylation in enhancer activity, in addition to the well-studied relationship between DNA methylation and promoter activity [39].

We previously reported on transcriptional data of amnion (n = 7), muscle (n = 6), adrenal (n = 3) and pancreas (n = 7) at W9, W18 and W22 [40] and used this data to test the hypothesis that GOM is associated with the epigenetic downregulation of developmental programs and LOM with upregulation of tissue-specific processes. Genes associated with GOM and involved in embryonic morphogenesis (a process enriched for GOM in all tissues) showed a decrease in transcriptional activity from W9 to W22 in all tissues (amnion, muscle, pancreas: P < 0.05; adrenal P = 0.87; Fig 4A). In contrast, genes involved in tissue-specific processes found to be enriched for LOM (S3 Table) increased in transcription from W9 to W22 (P < 0.05; Figs 4B and S4D).

Fig. 4. DNA methylation dynamics are accompanied by changes in gene expression.
DNA methylation dynamics are accompanied by changes in gene expression.
(A) Expression profiles of genes in embryonic morphogenesis near dynamic regions with gain of methylation represented as median with the interquartile range (IQR) [40]. (B) Expression profiles of genes near dynamic regions with loss of methylation grouped by significant, tissue-specific Gene Ontology terms for each of the four tissues from S3 Table represented as median with IQR [40]. Reg., regulation.

Altogether, these findings emphasize that DNA methylation dynamics during human fetal development is associated with the availability to transcription of both general embryonic programs (shutting those down for transcription) as well as tissue-specific developmental programs (making those available for transcription).

Dynamically methylated regions correlate with developmental and tissue-specific genes

From the dynamically methylated CpGs, we identified 2,229 development-related differentially methylated regions (dDMRs, defined as 3 consecutive differentially methylated CpGs within 1kb of each other) undergoing GOM (amnion 185; muscle 530; adrenal 1,065; pancreas 449) and 1,017 undergoing LOM (amnion 388; muscle 482; adrenal 136; pancreas 61; S4 Table). After mapping the dDMRs to the nearest gene locus we observed that the percentage of common genes in the embryonic tissues showing LOM dDMRs was 1.3%, whereas those showing GOM was 10.2% (Fig 5A).

Fig. 5. Association of gain and loss of DNA methylation, DNAse I hypersensitive sites and histone modifications.
Association of gain and loss of DNA methylation, DNAse I hypersensitive sites and histone modifications.
(A) Venn diagram visualizing the overlaps between genes with a gain and a loss of methylation of the three embryonic tissues. (B) Methylation difference between W9 and W22 of MYLK2 in the four tissues. (C) Methylation difference between W9 and W22 of the HOXB cluster. (D) Mean DNAse I hypersensitive (DHS) and histone modifications signals in a 5 kb flanking region of the muscle dDMRs in HSMMs and HSMMtubes. HSMM, human skeletal muscle myoblasts; HSMMtube, human skeletal muscle myotubes.

LOM dDMRs were associated with genes involved in tissue-specific functions such as MYH3 in muscle (muscle contractile protein), MC2R (adrenocorticotropic hormone receptor) in adrenal and PFKFB3 (involved in insulin secretion) in the pancreas (Table 2 and S5A Fig). As an example of tissue specificity encountered in the dDMRs showing a loss of DNA methylation, we zoomed in on the MYLK2 locus, a muscle-specific gene [41]. The methylation of the MYLK2 promoter and first exon in the muscle decreased during development, but increased (or remained constant) in the other organs studied (Fig 5B). Interestingly, GOM dDMRs were associated with (tissue-specific) developmental genes, such as PAX3 in muscle, a key gene in myogenesis [42], and NKX6.1 in pancreas, an important gene in beta-cell development [43] (Table 2 and S5A Fig), but also near well-known developmental genes including the HOXB (Fig 5C) and other HOX clusters (S6A Fig) that play a key role in embryonic patterning and morphogenesis [44].

Tab. 2. Table with six genes per tissue highlighting the tissue specificity of the genes found near dDMRs with a loss of methylation as well as the association of dDMRs with a gain of methylation with tissue-specific developmental genes.
Table with six genes per tissue highlighting the tissue specificity of the genes found near dDMRs with a loss of methylation as well as the association of dDMRs with a gain of methylation with tissue-specific developmental genes.

However, when comparing all identified dDMRs to previously identified adult (tissue-specific) tDMRs using the 450k array [6], about 50% of the GOM dDMRs were not identified as tDMRs in adult tissues (S5B Fig), while 32% and 38%, of the LOM dDMRs in muscle and pancreas, respectively, were unique for those fetal tissues (adult data on adrenal was absent). The persistence into adulthood of subsets of GOM and LOM dDMRs was confirmed using WGBS data for adult muscle [25] (S5C Fig). These results suggest that the study of DNA methylation dynamics in fetal development will identify regions that are remodelled during development and are missed when studying adult tissues only.

We further explored the potential biological validity of the 1,012 muscle dDMRs using ENCODE data [24] on human skeletal muscle myoblasts (HSMMs) and their differentiated derivatives, human skeletal muscle myotubes (HSMMtubes). In HSMMs and HSMMtubes, DNAse I hypersensitive sites (DHSs), which mark genomic regions of open chromatin associated with transcriptional activity, were abundant at LOM dDMRs, particularly in CpG-poor regions (CGI-shelves and non-CGI regions, Fig 5D left). DHSs were depleted at GOM dDMRs in CpG-rich regions (CGIs and CGI-shores, Fig 5D left). Consistent with an increased transcriptional activity, LOM dDMRs were also enriched in myotubes and myoblasts (ENCODE [24]) for histone H3 lysine 4 methylation ((H3K4me1, -me2, -me3) [25,45], and acetylation of histone H3 at lysine 9 and 27 (H3K9ac, H3K27ac), all marks associated with active regulatory regions [25,45] (Fig 5D right). In contrast, these active histone modifications were depleted for GOM dDMRs in CpG-rich regions (Fig 5D right). LOM dDMRs were depleted of H3K9me3 (marking inactive DNA), H3K27me3 (marking Polycomb-repressed regions) and H3K36me3 in HSMMtubes but not in their precursor cells HSMMs (Fig 5D right). A final indication for the functional relevance of the muscle dDMRs was that 124 out of 482 LOM dDMRs significantly overlapped (P < 0.0001) with binding sites of the muscle-specific transcription factor MYOD in HSMMs (188/482 in HSMMtubes), whereas only 7 out of the 530 GOM dDMRs mapped to MYOD binding sites[46] (8/530 in HSMMtubes; S5D Fig).

Discussion

Here, we show that human tissues already exhibit a specific DNA methylation signature as early as W9 of fetal development. In addition, the DNA methylation landscape is subjected to considerable changes from the first to second trimester of gestation as the developing organs gain complexity and functionality. Our study highlights that dynamic DNA methylation is not only an integral part of early preimplantation embryo development and implantation [911], but continues to be a key feature of epigenetic remodelling during human fetal development. While global changes in levels of DNA methylation characterize development until implantation (Fig 6), these are not observed during fetal development. Instead, distinct LOM occurs near tissue-specific genes and GOM occurs near developmental genes in a largely tissue-independent fashion (Fig 6). Our direct assessment of DNA methylation dynamics suggests that a larger proportion of the methylome is remodelled during development than previously thought [3,6,47].

Fig. 6. DNA methylation dynamics during human development.
DNA methylation dynamics during human development.
This illustration depicts the current comprehensive knowledge of DNA methylation during human pre- and postimplantation development. The knowledge about DNA methylation during human preimplantation (left panel) is derived from [912] whereas our study sheds light on postimplantation development (right panel).

Interestingly, the functional relevance of identified dynamic regions was further exemplified by the changes in expression of their nearest genes. While the nearest genes of regions gaining DNA methylation associated with embryonic morphogenesis showed loss of expression, the nearest genes of regions losing DNA methylation showed increased expression over time. In agreement with our observations, LOM of hematopoietic-specific genes has been observed during human hematopoietic differentiation in vitro [48] and have been linked to transcriptional changes in human T-cell development [49]. Moreover, several mouse and human in vitro studies demonstrated that the methylation of developmental genes increases [19,50] and tissue-specific functional genes lose methylation [5,51] during stem and progenitor cell differentiation. Lastly, DNA demethylating agents, such as 5-azacytidine, have been shown to promote stem cell differentiation and maturation of skeletal myotubes in mice [52,53]. Further experimental studies are required to evaluate the mechanistic role of DNA methylation in development.

Although our study reveals general principles of DNA methylation dynamics during human fetal development, it should be noted that a limited number of tissues and individuals was investigated; and that we used a genome-wide method interrogating a relatively small proportion of all CpGs in the human genome. Expansion to more tissues and the application of whole-methylome technologies will lead to a more comprehensive catalogue of regulatory regions. However, by extensive inclusion of external fetal and adult 450k array datasets, we have consolidated our findings. Moreover, the use of external available WGBS data confirmed the results obtained by the 450k array data.

Since we studied organ biopsies similar to previous studies investigating biopsies of human adult tissues [3,15], the methylation profiles we report reflect the average of multiple cell types. The cellular complexity of the organs investigated led to an underestimation of the actual DNA methylation dynamics in individual cell types. This is exemplified by the detection of a considerably larger number of CpGs displaying dynamic methylation in muscle which has an exclusive mesodermal origin in comparison with adrenal and pancreas which are composed by cells originating from two different germ layers (adrenal: mesoderm and ectoderm (neural crest); pancreas: endoderm and mesoderm). However, it is unlikely that the methylation dynamics observed is an epiphenomenon of this cellular complexity instead of being driven by cell differentiation and maturation. This is obvious for genes associated with GOM that appears to be shared across organs to repress general developmental programs during development. In contrast, genes associated with LOM displayed tissue-specific patterns. Their intricate involvement in organ-specific functions was emphasized by tight linkage to biological processes and chromatin states relevant to the organs investigated. Moreover, between W9 and W22, the organs analysed are mainly composed of progenitor cells; perfusion by blood and lymphatic vasculature, and innervation by neural crest cell derivatives still plays a minor role as compared with adult organs (S1A Fig). In the future, single-cell methodology [54,55] will enable comparing single-cell DNA methylomes of the various adult cell types to their fetal progenitor counterparts.

Studies of DNA methylation landscapes of human fetal development may serve as reference in the development of (organoid) differentiation models [56] and, moreover, shed light on potential mechanisms underlying genetic associations and studies in the field of epigenetic epidemiology [57] focussing on the prenatal environment.

Materials and Methods

Ethical statement

The Medical Ethical Committee of the Leiden University Medical Center approved this study (P08.087). Informed consent was obtained on the basis of the Declaration of Helsinki (World Medical Association).

Fetal tissue

Human fetal tissues (amnion, skeletal muscle, adrenal glands, pancreas) at gestational age W9, W18, W22 (S2A Fig) were collected from elective abortion material (vacuum aspiration) without medical indication. In this study, “weeks of gestation” was used as determined by the last menstrual period (LMP). After collection, the material was washed with 0.9% NaCl (Fresenius Kabi, France) and the identified organs were immediately snap-frozen using dry ice and stored at -80°C until further processing. Histology was performed as previously described [58]. The images were taken with an Olympus AX70 microscope (Olympus, Japan) provided with a XC50 digital colour camera (Olympus, Japan).

DNA extraction

Tissues were homogenized with a pestle and lysed overnight at 56°C with proteinase K (600 mAU/ml, Qiagen, Germany) in ATL buffer (Qiagen, Germany). After lysis, residual RNA in the samples was degraded using RNase A (10 mg/μl, Invitrogen, USA). Subsequently, genomic DNA (gDNA) was extracted on the basis of phenol/chloroform. Briefly, lysates were transferred to Phase Lock Heavy Gel 2ml Eppendorf tubes (5PRIME, Germany) and 700 μl of 25:24:1 Phenol/Chloroform/Isoamyl alcohol was added and spun down for 5 minutes. The aqueous phase was transferred to a Phase Lock tube and the latter step was repeated. The aqueous phase was transferred to a new Phase Lock tube and 700 μl 24:1 Chloroform/Isoamyl alcohol was added and spun down for 5 minutes. The aqueous phase was transferred to a Phase Lock tube and the latter step was repeated. The aqueous phase was transferred to a new 2 μl Eppendorf tube (Eppendorf AG, Germany), 70 μl 3M sodium acetate (Ambion, USA) and 1400 μl ice cold 100% ethanol were added. gDNA was precipitated over night at -20°C. Eppendorf tubes were spun down at 4°C for 15 minutes and washed twice with 70% ethanol. After the pellet was dry, gDNA was solubilized in AE buffer (Qiagen, Germany) and stored at 4°C. DNA concentration was determined using the Qubit dsDNA BR Assay Kit on a Qubit 2.0 Fluorometer (Invitrogen, USA). gDNA was bisulfite converted using the EZ-96 DNA methylation kit (Zymo Research, Orange County, USA) with an average input of 600 ng gDNA. Following bisulfite conversion, DNA methylation data was generated using Illumina HumanMethylation450 BeadChip according to the manufacturer’s protocol.

(Pre-) processing of the Illumina 450k BeadChip data

All analyses were performed using R statistics, version 3.0.1. The 65 polymorphic SNP probes featured on the 450k array were used to exclude potential sample mix ups. Data was imported in R using minfi [59] and processed and normalized using a custom pipeline: Arrays were removed if they had a low median intensity, a high background signal or with incomplete bisulfite conversion, but none were excluded (S1B, S1C and S1D Fig). The CpGs on chromosome X and Y were used to confirm sex (S1C Fig). Next, probes with a low bead count (< 3), high detection P-value (> 0.01), and with a low success rate (< 95%), and ambiguously mapped probes [60] were removed. After probe filtering, all arrays contained > 95% of the original number of probes. Background correction and colour correction were applied and the data was quantile normalized (lumi [61]). To adjust for the type I/II bias BMIQ was applied [62]. In our analyses, CpG sites in the sex chromosomes (Y and X) were excluded. High correlation was found between samples from the same time point and tissue but also between time points of the same tissue (S1D Fig). To exclude chromosomal abnormalities, we calculated the copy number aberration based on the signal intensities using the method published by Feber et al. [63] as implemented in the R package ChAMP [64]. From these results, no abnormalities were found in the samples used (S1E Fig).

Bioinformatics analyses

Multidimensional scaling and clustering was performed based on Euclidean distance. For the DNA methylation over time within features, a genic annotation was combined with a CGI-centric annotation as presented before [6] to determine the median methylation per combined feature. P were calculated using quantile regression based on the median (R package quantreg [65]). Figures were made using the R-packages ggplot2 [66] and GenomeGraphs [67].

Tissue-specific hypomethylation: CpGs with a standard deviation ≥ 0.1 within the tissue of interest were discarded. Relative tissue-specific hypomethylation was defined as hypomethylation of the tissue of interest compared to the other tissues, with a difference of ≥ 0.20 in beta value. CpGs sites were selected if a difference was consistent in each of the time points.

Dynamic methylation: CpGs with a high standard deviation ≥ 0.10 within time point/tissues indicative of an instable estimate of DNA methylation were discarded from this analysis. Gain and loss of methylation was defined a gain/loss of ≥ 0.20 between W9 and W22 and W18 in between the two time points (W18 was allowed to be 0.05 lower/higher in beta value than W9/W22 respectively). The CpGs with a gain or a loss of methylation were used for the combined genic/CGI-centric annotation and expressed as an odds ratio. Chromatin state segmentation data [25] were downloaded from the Epigenomics Roadmap Project for fetal muscle (W15 female), fetal adrenal (W13 male), amnion (W16 male) and pancreatic islets (adult) and enrichment of dynamically methylated CpGs was calculated.

DMRs: In both, the relative hypomethylated CpGs (tHRs) and the CpGs with a gain or a loss of methylation (dDMRs), DMRs were called using an algorithm described before [6]. Briefly, DMRs (tHRs and dDMRs) were defined by three consecutive CpGs that matched a criterion (that is, relative hypomethylation or gain/loss of methylation) with a maximum of 1 kb between CpGs and with at highest three CpGs that did not match the criterion.

Gene ontology: Tissue-specific hypomethylated and dynamically methylated CpGs were mapped to their nearest gene (that is to the nearest TSS or TES of a gene) and tested for enrichment of gene ontology terms using DAVID [68]. For the tissue-specific hypomethylation we used a P cut-off of 0.05 on the raw P as the number of CpGs was relatively low. For the dynamically methylated CpGs a FDR cut-off of 0.05 was set as cut-off for enrichment in GO terms. A background set was used containing nearest genes of all CpGs covered on the array.

Gene expression data: Transcriptional data of the four tissues at W9, W18, W22 (amnion: n = 2, 3, 2, muscle: n = 2, 2, 2; adrenal: n = 1, 1, 1; pancreas: n = 3, 2, 2) were used. The counts per million (CPM) expression levels were calculated using the R package edgeR 3.2.4 [69,70]. For the plots, the arithmetic mean of the biological replicates was used and the median of all genes plotted. To access enrichment of up- and downregulation, a probability test was used.

MYOD, DNAse I, histone marks and WGBS: Overlaps between MYOD binding peaks and muscle dDMRs were calculated. To test for significance, we calculated an empirical distribution by performing 20,000 permutations with 482 (gain of methylation) and 530 (loss of methylation) DMR-like regions each and determined the overlap with the MYOD binding sites. DMR-like regions were defined as regions with equal characteristics as dDMRs identified: an inter-CpG distance smaller than 1 kb and an average length of five CpGs per DMR-like region (n~8 x 104 regions). The two-sided P was determined using the empirical distribution.

DNAse I and histone mark data of human skeletal muscle myoblasts (HSMMs) and human skeletal muscle myotubes (HSMMtubes) were downloaded from the ENCODE website [24]. DNAse I hypersensitivity was expressed as the count of DNAse-seq tags. The enrichment of histone marks was expressed as the log2 of the ChIP/input. The total number of reads within the myotubes was different from the total number of reads in the myoblast data and, therefore, the data was normalized. dDMRs were classified as island (CGIs and their shores) or non-island dDMRs, and DNAse-seq tags and histone marks around dDMRs were mapped up to 5 kb up- and downstream.

CpG sites of WGBS data were mapped to hypomethylated and dynamic regions and their 5kb flanking regions. Using a smooth spline, the methylation around the regions was smoothed for the adult and fetal data.

Accession numbers

Methylation data has been deposited in the NCBI’s Gene Expression Omnibus [71] under accession number GSE56515. External datasets that have been used in this manuscript include: fetal and adult DNA methylation data of various tissues from Nazor et al. (Gene Expression Omnibus (GEO) accession number: GSE31848) [15], fetal brain DNA methylation data from Spiers et al. (GEO accession number: GSE58885) [22], fetal liver DNA methylation data from Bonder et al. (GEO accession number: GSE61279) [30], adult DNA methylation data of various tissues from Slieker et al. (GEO accession number: GSE48472) [6], fetal Deep SAGE expression data of the four tissues studied here from Roost et al. (GEO accession number: GSE66302) [40], adult DNA methylation brain data from Pidsley et al.(GEO accession number: GSE61431) [38], WGBS data of fetal and adult muscle generated by the Epigenomics Roadmap consortium (GEO accession numbers: GSM1172596 and GSM1010986), MYOD binding peaks from MacQuarrie et al. (GEO accession numbers: GSM1218849 and GSM1218850) [46].

Supporting Information

Attachment 1

Attachment 2

Attachment 3

Attachment 4

Attachment 5

Attachment 6

Attachment 7

Attachment 8

Attachment 9

Attachment 10


Zdroje

1. Byun HM, Siegmund KD, Pan F, Weisenberger DJ, Kanel G, Laird PW, et al. Epigenetic profiling of somatic tissues from human autopsy specimens identifies tissue- and individual-specific DNA methylation patterns. Hum Mol Genet. 2009;18: 4808–4817. doi: 10.1093/hmg/ddp445 19776032

2. Rakyan VK, Down TA, Thorne NP, Flicek P, Kulesha E, Gräf S, et al. An integrated resource for genome-wide identification and analysis of human tissue-specific differentially methylated regions (tDMRs). Genome Res. 2008;18: 1518–1529. doi: 10.1101/gr.077479.108 18577705

3. Ziller MJ, Gu H, Müller F, Donaghey J, Tsai LT- Y, Kohlbacher O, et al. Charting a dynamic DNA methylation landscape of the human genome. Nature. Nature Publishing Group, a division of Macmillan Publishers Limited. All Rights Reserved.; 2013;500: 477–81. doi: 10.1038/nature12433 23925113

4. Hon GC, Rajagopal N, Shen Y, McCleary DF, Yue F, Dang MD, et al. Epigenetic memory at embryonic enhancers identified in DNA methylation maps from adult mouse tissues. Nat Genet. Nature Publishing Group, a division of Macmillan Publishers Limited. All Rights Reserved.; 2013;45: 1198–206. doi: 10.1038/ng.2746 23995138

5. Nagae G, Isagawa T, Shiraki N, Fujita T, Yamamoto S, Tsutsumi S, et al. Tissue-specific demethylation in CpG-poor promoters during cellular differentiation. Hum Mol Genet. 2011;20: 2710–2721. doi: 10.1093/hmg/ddr170 21505077

6. Slieker RC, Bos SD, Goeman JJ, Bovée JV, Talens RP, van der Breggen R, et al. Identification and systematic annotation of tissue-specific differentially methylated regions using the Illumina 450k array. Epigenetics Chromatin. 2013;6: 26. doi: 10.1186/1756-8935-6-26 23919675

7. Gifford C a, Ziller MJ, Gu H, Trapnell C, Donaghey J, Tsankov A, et al. Transcriptional and epigenetic dynamics during specification of human embryonic stem cells. Cell. Elsevier Inc.; 2013;153: 1149–1163. doi: 10.1016/j.cell.2013.04.037 23664763

8. Maunakea AK, Nagarajan RP, Bilenky M, Ballinger TJ, D’Souza C, Fouse SD, et al. Conserved role of intragenic DNA methylation in regulating alternative promoters. Nature. Nature Publishing Group; 2010;466: 253–257. doi: 10.1038/nature09165 20613842

9. Smith ZD, Chan MM, Humm KC, Karnik R, Mekhoubad S, Regev A, et al. DNA methylation dynamics of the human preimplantation embryo. Nature. Nature Publishing Group, a division of Macmillan Publishers Limited. All Rights Reserved.; 2014;511: 611–615. doi: 10.1038/nature13581 25079558

10. Guo H, Zhu P, Yan L, Li R, Hu B, Lian Y, et al. The DNA methylation landscape of human early embryos. Nature. Nature Publishing Group, a division of Macmillan Publishers Limited. All Rights Reserved.; 2014;511: 606–10. doi: 10.1038/nature13544 25079557

11. Okae H, Chiba H, Hiura H, Hamada H, Sato A, Utsunomiya T, et al. Genome-wide analysis of DNA methylation dynamics during early human development. Oakey RJ, editor. PLoS Genet. 2014;10: e1004868. doi: 10.1371/journal.pgen.1004868 25501653

12. Seisenberger S, Peat JR, Hore TA, Santos F, Dean W, Reik W. Reprogramming DNA methylation in the mammalian life cycle: building and breaking epigenetic barriers. Philos Trans R Soc Lond B Biol Sci. 2013;368: 20110330. doi: 10.1098/rstb.2011.0330 23166394

13. Smith ZD, Chan MM, Mikkelsen TS, Gu H, Gnirke A, Regev A, et al. A unique regulatory phase of DNA methylation in the early mammalian embryo. Nature. Nature Publishing Group; 2012. pp. 339–344. doi: 10.1038/nature10960 22456710

14. Seisenberger S, Andrews S, Krueger F, Arand J, Walter J, Santos F, et al. The Dynamics of Genome-wide DNA Methylation Reprogramming in Mouse Primordial Germ Cells. Mol Cell. 2012;48: 849–862. doi: 10.1016/j.molcel.2012.11.001 23219530

15. Nazor KL, Altun G, Lynch C, Tran H, Harness J V., Slavin I, et al. Recurrent variations in DNA methylation in human pluripotent stem cells and their differentiated derivatives. Cell Stem Cell. 2012;10: 620–634. doi: 10.1016/j.stem.2012.02.013 22560082

16. Bocker MT, Hellwig I, Breiling A, Eckstein V, Ho AD, Lyko F. Genome-wide promoter DNA methylation dynamics of human hematopoietic progenitor cells during differentiation and aging. Blood. 2011;117: e182–9. doi: 10.1182/blood-2011-01-331926 21427290

17. Xie W, Schultz MD, Lister R, Hou Z, Rajagopal N, Ray P, et al. Epigenomic Analysis of Multilineage Differentiation of Human Embryonic Stem Cells. Cell. Elsevier Inc.; 2013;153: 1134–1148. doi: 10.1016/j.cell.2013.04.022 23664764

18. Brunner AL, Johnson DS, Kim SW, Valouev A, Reddy TE, Neff NF, et al. Distinct DNA methylation patterns characterize differentiated human embryonic stem cells and developing human fetal liver. Genome Res. 2009;19: 1044–1056. doi: 10.1101/gr.088773.108 19273619

19. Laurent L, Wong E, Li G, Huynh T, Tsirigos A, Ong CT, et al. Dynamic changes in the human methylome during differentiation. Genome Res. 2010;20: 320–331. doi: 10.1101/gr.101907.109 20133333

20. Bock C, Beerman I, Lien WH, Smith ZD, Gu H, Boyle P, et al. DNA Methylation Dynamics during In Vivo Differentiation of Blood and Skin Stem Cells. Mol Cell. 2012;47: 633–647. doi: 10.1016/j.molcel.2012.06.019 22841485

21. Gilsbach R, Preissl S, Grüning BA, Schnick T, Burger L, Benes V, et al. Dynamic DNA methylation orchestrates cardiomyocyte development, maturation and disease. Nat Commun. Nature Publishing Group; 2014;5: 5288. doi: 10.1038/ncomms6288 25335909

22. Spiers H, Hannon E, Schalkwyk LC, Smith R, Wong CCY, O’Donovan MC, et al. Methylomic trajectories across human fetal brain development. Genome Res. 2015;25: 338–352. doi: 10.1101/gr.180273.114 25650246

23. Bibikova M, Barnes B, Tsan C, Ho V, Klotzle B, Le JM, et al. High density DNA methylation array with single CpG site resolution. Genomics. Elsevier Inc.; 2011;98: 288–295. doi: 10.1016/j.ygeno.2011.07.007 21839163

24. Dunham I, Kundaje A, Aldred SF, Collins PJ, Davis CA, Doyle F, et al. An integrated encyclopedia of DNA elements in the human genome. Nature. Nature Publishing Group, a division of Macmillan Publishers Limited. All Rights Reserved.; 2012. pp. 57–74. doi: 10.1038/nature11247 22955616

25. Ernst J, Kheradpour P, Mikkelsen TS, Shoresh N, Ward LD, Epstein CB, et al. Mapping and analysis of chromatin state dynamics in nine human cell types. Nature. 2011;473: 43–49. doi: 10.1038/nature09906 21441907

26. Hou L, Zhao H. A review of post-GWAS prioritization approaches [Internet]. Frontiers in Genetics. 2013. p. 280. doi: 10.3389/fgene.2013.00280 24367376

27. Waterland RA, Michels KB. Epigenetic epidemiology of the developmental origins hypothesis. Annu Rev Nutr. Annual Reviews; 2007;27: 363–88. Available: http://www.annualreviews.org/doi/abs/10.1146/annurev.nutr.27.061406.093705 17465856

28. Heijmans BT, Tobi EW, Stein AD, Putter H, Blauw GJ, Susser ES, et al. Persistent epigenetic differences associated with prenatal exposure to famine in humans. Proc Natl Acad Sci U S A. 2008;105: 17046–17049. doi: 10.1073/pnas.0806560105 18955703

29. Sandoval J, Heyn HA, Moran S, Serra-Musach J, Pujana MA, Bibikova M, et al. Validation of a DNA methylation microarray for 450,000 CpG sites in the human genome. Epigenetics. 2011;6: 692–702. 21593595

30. Bonder MJ, Kasela S, Kals M, Tamm R, Lokk K, Barragan I, et al. Genetic and epigenetic regulation of gene expression in fetal and adult human livers. BMC Genomics. 2014;15: 860. doi: 10.1186/1471-2164-15-860 25282492

31. Jones PA. Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. Nature Publishing Group, a division of Macmillan Publishers Limited. All Rights Reserved.; 2012;13: 484–92. doi: 10.1038/nrg3230 22641018

32. Colaneri A, Wang T, Pagadala V, Kittur J, Staffa NG, Peddada SD, et al. A Minimal Set of Tissue-Specific Hypomethylated CpGs Constitute Epigenetic Signatures of Developmental Programming. Defossez P-A, editor. PLoS One. 2013;8: e72670. doi: 10.1371/journal.pone.0072670 24069155

33. Song F, Mahmood S, Ghosh S, Liang P, Smiraglia DJ, Nagase H, et al. Tissue specific differentially methylated regions (TDMR): Changes in DNA methylation during development. Genomics. 2009;93: 130–139. doi: 10.1016/j.ygeno.2008.09.003 18952162

34. Oliver-Krasinski JM, Stoffers D a. On the origin of the β cell. Genes and Development. 2008. pp. 1998–2021. doi: 10.1101/gad.1670808 18676806

35. Brink C. Promoter elements in endocrine pancreas development and hormone regulation. Cell Mol Life Sci. 2003;60: 1033–1048. 12861373

36. Leitges M, Neidhardt L, Haenig B, Herrmann BG, Kispert a. The paired homeobox gene Uncx4.1 specifies pedicles, transverse processes and proximal ribs of the vertebral column. Development. 2000;127: 2259–67. 10804169

37. Schulz R a, Yutzey KE. Calcineurin signaling and NFAT activation in cardiovascular and skeletal muscle development. Developmental Biology. 2004. pp. 1–16. 14729474

38. Pidsley R, Viana J, Hannon E, Spiers H, Troakes C, Al-Saraj S, et al. Methylomic profiling of human brain tissue supports a neurodevelopmental origin for schizophrenia. Genome Biol. 2014;15: 483. 25347937

39. Wiench M, John S, Baek S, Johnson TA, Sung M-H, Escobar T, et al. DNA methylation status predicts cell type-specific enhancer activity. EMBO J. 2011;30: 3028–3039. doi: 10.1038/emboj.2011.210 21701563

40. Roost MS, van Iperen L, Ariyurek Y, Buermans HP, Arindrarto W, Devalla HD, et al. KeyGenes, a Tool to Probe Tissue Differentiation Using a Human Fetal Transcriptional Atlas. Stem cell reports. 2015;4: 1112–24. doi: 10.1016/j.stemcr.2015.05.002 26028532

41. Stull JT, Kamm KE, Vandenboom R. Myosin light chain kinase and the role of myosin light chain phosphorylation in skeletal muscle. Archives of Biochemistry and Biophysics. 2011. pp. 120–128. doi: 10.1016/j.abb.2011.01.017 21284933

42. Bismuth K, Relaix F. Genetic regulation of skeletal muscle development. Experimental Cell Research. 2010. pp. 3081–3086. doi: 10.1016/j.yexcr.2010.08.018 20828559

43. Nelson SB, Schaffer AE, Sander M. The transcription factors Nkx6.1 and Nkx6.2 possess equivalent activities in promoting beta-cell fate specification in Pdx1+ pancreatic progenitor cells. Development. 2007;134: 2491–2500. 17537793

44. Heffer A, Pick L. Conservation and variation in Hox genes: how insect models pioneered the evo-devo field. Annu Rev Entomol. Annual Reviews; 2013;58: 161–79. doi: 10.1146/annurev-ento-120811-153601 23317041

45. Zhou VW, Goren A, Bernstein BE. Charting histone modifications and the functional organization of mammalian genomes. Nat Rev Genet. 2011;12: 7–18. doi: 10.1038/nrg2905 21116306

46. MacQuarrie KL, Yao Z, Fong AP, Diede SJ, Rudzinski ER, Hawkins DS, et al. Comparison of genome-wide binding of MyoD in normal human myogenic cells and rhabdomyosarcomas identifies regional and local suppression of promyogenic transcription factors. Mol Cell Biol. 2013;33: 773–84. doi: 10.1128/MCB.00916-12 23230269

47. Lokk K, Modhukur V, Rajashekar B. DNA methylome profiling of human tissues identifies global and tissue-specific methylation patterns. Genome …. 2014; Available: http://www.biomedcentral.com/content/pdf/gb-2014-15-4-r54.pdf

48. Calvanese V, Fernández AF, Urdinguio RG, Suárez-Alvarez B, Mangas C, Pérez-García V, et al. A promoter DNA demethylation landscape of human hematopoietic differentiation. Nucleic Acids Res. 2012;40: 116–131. doi: 10.1093/nar/gkr685 21911366

49. Rodriguez R. Regulation of the transcriptional program by DNA methylation during human αβ T-cell development. Nucleic acids …. 2014; Available: http://nar.oxfordjournals.org/content/early/2014/12/24/nar.gku1340.short

50. Kim M, Park YK, Kang TW, Lee SH, Rhee YH, Park JL, et al. Dynamic changes in DNA methylation and hydroxymethylation when hES cells undergo differentiation toward a neuronal lineage. Hum Mol Genet. 2014;23: 657–667. doi: 10.1093/hmg/ddt453 24087792

51. Miyata K, Miyata T, Nakabayashi K, Okamura K, Naito M, Kawai T, et al. DNA methylation analysis of human myoblasts during in vitro myogenic differentiation: de novo methylation of promoters of muscle-related genes and its involvement in transcriptional down-regulation. Hum Mol Genet. 2015;24: 410–23. doi: 10.1093/hmg/ddu457 25190712

52. Zhou Y, Hu Z. Genome-Wide Demethylation by 5-aza-2’-Deoxycytidine Alters the Cell Fate of Stem/Progenitor Cells. Stem Cell Reviews and Reports. 2014.

53. Hupkes M, Jonsson MKB, Scheenen WJ, van Rotterdam W, Sotoca AM, van Someren EP, et al. Epigenetics: DNA demethylation promotes skeletal myotube maturation [Internet]. The FASEB Journal. 2011. pp. 3861–3872. doi: 10.1096/fj.11-186122 21795504

54. Smallwood SA, Lee HJ, Angermueller C, Krueger F, Saadeh H, Peat J, et al. Single-cell genome-wide bisulfite sequencing for assessing epigenetic heterogeneity. Nat Methods. 2014;11: 817–20. doi: 10.1038/nmeth.3035 25042786

55. Shapiro E, Biezuner T, Linnarsson S. Single-cell sequencing-based technologies will revolutionize whole-organism science. Nat Rev Genet. 2013;14: 618–30. doi: 10.1038/nrg3542 23897237

56. Lancaster MA, Knoblich JA. Organogenesis in a dish: modeling development and disease using organoid technologies. Science. 2014;345: 1247125. doi: 10.1126/science.1247125 25035496

57. Mill J, Heijmans BT. From promises to practical strategies in epigenetic epidemiology. Nat Rev Genet. 2013;14: 585–94. doi: 10.1038/nrg3405 23817309

58. Roost MS, van Iperen L, de Melo Bernardo A, Mummery CL, Carlotti F, de Koning EJ, et al. Lymphangiogenesis and angiogenesis during human fetal pancreas development. Vasc Cell. BioMed Central Ltd; 2014;6: 22. doi: 10.1186/2045-824X-6-22 25785186

59. Hansen KD, Aryee M. minfi: Analyze Illumina’s 450k methylation arrays. R package version 1.6.0. Bioconductor; 2013.

60. Chen Y, Lemire M, Choufani S, Butcher DT, Grafodatskaya D, Zanke BW, et al. Discovery of cross-reactive probes and polymorphic CpGs in the Illumina Infinium HumanMethylation450 microarray. Epigenetics. 2013;8: 203–9. doi: 10.4161/epi.23470 23314698

61. Du P, Kibbe WA, Lin SM. lumi: a pipeline for processing Illumina microarray. Bioinformatics. 2008;24: 1547–8. doi: 10.1093/bioinformatics/btn224 18467348

62. Teschendorff AE, Marabita F, Lechner M, Bartlett T, Tegner J, Gomez-Cabrero D, et al. A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics. 2013;29: 189–96. doi: 10.1093/bioinformatics/bts680 23175756

63. Feber A, Guilhamon P, Lechner M, Fenton T, Wilson GA, Thirlwell C, et al. Using high-density DNA methylation arrays to profile copy number alterations. Genome Biol. 2014;15: R30. doi: 10.1186/gb-2014-15-2-r30 24490765

64. Morris TJ, Butcher LM, Feber A, Teschendorff AE, Chakravarthy AR, Wojdacz TK, et al. ChAMP: 450k Chip Analysis Methylation Pipeline. Bioinformatics. 2014;30: 428–430. doi: 10.1093/bioinformatics/btt684 24336642

65. Koenker R. quantreg: Quantile Regression [Internet]. 2013. Available: http://cran.r-project.org/package=quantreg

66. Wickham H. ggplot2: elegant graphics for data analysis. Springer, New York. Springer, New York; 2009.

67. Durinck S, Bullard J. GenomeGraphs: Plotting genomic information from Ensembl. R package version 1.22.0.

68. Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. Nature Publishing Group; 2009;4: 44–57. doi: 10.1038/nprot.2008.211 19131956

69. Robinson MD, McCarthy DJ, Smyth GK. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2009;26: 139–140. doi: 10.1093/bioinformatics/btp616 19910308

70. Robinson MD, Smyth GK. Small-sample estimation of negative binomial dispersion, with applications to SAGE data. Biostatistics. 2008;9: 321–32. Available: http://biostatistics.oxfordjournals.org/content/9/2/321.short 17728317

71. Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, et al. NCBI GEO: Archive for functional genomics data sets—Update. Nucleic Acids Res. 2013;41: D991–5. doi: 10.1093/nar/gks1193 23193258

Štítky
Genetika Reprodukční medicína

Článek vyšel v časopise

PLOS Genetics


2015 Číslo 10
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

plice
INSIGHTS from European Respiratory Congress
nový kurz

Současné pohledy na riziko v parodontologii
Autoři: MUDr. Ladislav Korábek, CSc., MBA

Svět praktické medicíny 3/2024 (znalostní test z časopisu)

Kardiologické projevy hypereozinofilií
Autoři: prof. MUDr. Petr Němec, Ph.D.

Střevní příprava před kolonoskopií
Autoři: MUDr. Klára Kmochová, Ph.D.

Všechny kurzy
Kurzy Podcasty Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#