GTC
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE ADVANCED SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Genes to Cells (2007) 12, 1123-1132. doi:10.1111/j.1365-2443.2007.01120.x
© 2007 Blackwell Publishing or its licensors

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Sakamoto, H.
Right arrow Articles by Shiota, K.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Sakamoto, H.
Right arrow Articles by Shiota, K.

Cell type-specific methylation profiles occurring disproportionately in CpG-less regions that delineate developmental similarity

Hideki Sakamoto1,2,a,b, Masako Suzuki1,a,c, Tetsuya Abe1, Tohru Hosoyama1, Emi Himeno1, Satoshi Tanaka1, John M Greally2, Naka Hattori1, Shintaro Yagi1 and Kunio Shiota1,*

1 Cellular Biochemistry, Animal Resource Sciences/Veterinary Medical Sciences, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo 113-8657, Japan
2 Departments of Medicine (Hematology) and Molecular Genetics, Albert Einstein College of Medicine, Bronx, NY, USA


    Abstract
 Top
 Abstract
 Introduction
 Results
 Discussion
 Experimental procedures
 References
 
Our previous studies using restriction landmark genomic scanning (RLGS) defined tissue- or cell-specific DNA methylation profiles. It remains to be determined whether the DNA sequence compositions in the genomic contexts of the NotI loci tested by RLGS influence their tendency to change with differentiation. We carried out 3834 methylation measurements consisting of 213 NotI loci in the mouse genome in 18 different tissues and cell types, using quantitative real-time PCR based on a VIRTUAL IMAGE RLGS database. Loci were categorized as CpG islands or other, and as unique or repetitive sequences, each category being associated with a variety of methylation categories. Strikingly, the tissue-dependently and differentially methylated regions (T-DMRs) were disproportionately distributed in the non-CpG island loci. These loci were located not only in 5'-upstream regions of genes but also in intronic and non-genic regions. Hierarchical clustering of the methylation profiles could be used to define developmental similarity and cellular phenotypes. The results show that distinctive tissue- and cell type-specific methylation profiles by RLGS occur mostly at NotI sites located at non-CpG island sequences, which delineate developmental similarity of different cell types. The finding indicates the power of NotI methylation profiles in evaluating the relatedness of different cell types.


    Introduction
 Top
 Abstract
 Introduction
 Results
 Discussion
 Experimental procedures
 References
 
DNA methylation is associated with transcriptional silencing and chromatin condensation, involved in various phenomena such as tissue-specific gene expression and tumorigenesis in mammals (Bird & Wolffe 1999; Jones & Takai 2001). In mammalian genomes, CpGs are found in dense clusters called CpG islands, which are frequently located at the 5'-ends of housekeeping and tissue-specific genes (Bird 1987; Gardiner-Garden & Frommer 1987; Larsen et al. 1992). Restriction landmark genomic scanning (RLGS) using NotI, a methylation-sensitive restriction enzyme recognizing the 8-bp GCGGCCGC motif, has been employed to define tissue-specific DNA methylation profiles (Ohgane et al. 1998). Interestingly, several instances of tissue-specific methylation of CpG islands have been found to be inversely correlated with gene expression at those loci (Imamura et al. 2001; Newell-Price et al. 2001; Pao et al. 2001; Futscher et al. 2002). These data suggest that CpG islands are more complex than merely representing constitutively unmethylated loci. Methylation of the deoxycytidine of CpG dinucleotides is heritable from parent to daughter cells, but alters with mammalian development and differentiation, with DNA methylation profiles revealing numerous tissue- or cell type-dependently and differentially methylated regions (T-DMRs) (Shiota 2004). These profiles can discriminate cell types but also measure the similarity between them by defining as an "epigenetic distance," the number of differential RLGS spots (Shiota et al. 2002).

Recent progress in understanding genome sequence information has raised questions about these gene- and CpG island-centric views. Genes without CpG islands have been also reported recently to have tissue-specific methylation in 5'-flanking regions (Cho et al. 2001; Hattori et al. 2004b; Nishino et al. 2004; Ko et al. 2005). While a previous study has estimated that more than 90% of NotI sites represent CpG islands in mammalian genome (Lindsay & Bird 1987), the recent completion of mouse genome sequencing reveals that 33% of NotI sites distributes in non-CpG islands (Fazzari & Greally 2004). Meanwhile, previous RLGS-based studies have assumed that repetitive elements are rarely detected in RLGS screens (Shiota et al. 2002; Matsuyama et al. 2003), as repetitive elements are generally methylated to silence their transposable activity (Yoder et al. 1997). However, even some repetitive elements can be differentially methylated between tissues (Chalitchagorn et al. 2004; Burden et al. 2005), and analysis of genome sequence database also showed that 36% of NotI loci represent repetitive elements in the mouse genome (Fazzari & Greally 2004).

We have previously defined epigenetic distance as number of differentially methylated RLGS loci, and found that this distance reflects similarities in cell differentiation lineages (Shiota et al. 2002). We have previously analyzed the methylation profiles of randomly chosen NotI sites from RLGS studies, defining the genomic source of the differentially methylated loci using Virtual Image RLGS (Vi-RLGS), which is an in silico RLGS simulation based on genome sequence predictions of restriction enzyme digestions (Matsuyama et al. 2003; Hattori et al. 2004a). In the present study, to examine the distribution of loci in the DNA methylation status, we analyzed the methylation status of 213 potentially informative NotI loci identified using Vi-RLGS of 18 types of mouse tissues and cells by methylation-sensitive quantitative real-time PCR and characterized the T-DMRs. In addition, we analyzed the methylation profiles of T-DMRs by agglomerative hierarchical clustering to evaluate the correlation between cellular phenotype and epigenetic status.


    Results
 Top
 Abstract
 Introduction
 Results
 Discussion
 Experimental procedures
 References
 
Sampling of the NotI sites represented in the RLGS profiles

For sampling non-redundant NotI loci identified by RLGS, the Vi-RLGS software was used (Matsuyama et al. 2003; Hattori et al. 2004a). The sequences of 213 non-redundant Vi-RLGS spots, each of which represents one NotI-PstI fragment, were sampled from 2170 spots of the Vi-RLGS image. It has been reported that the mouse genome has 6057 NotI sites that are distributed among every chromosome in the genome (Fazzari & Greally 2004). The 213 NotI loci we used represent every chromosome except for the Y chromosome, where there are very few NotI sites (Supplementary Table S1). Of the 213 sampled NotI sites, some were located in CpG islands or non-CpG island sequences, and could be categorized further as unique or repetitive sequences (Fig. 1A). Repetitive elements were underrepresented compared to the total distribution of NotI loci because the redundant repetitive NotI-PstI fragments, which appear as aggregated large spots in the Vi-RLGS profile, were omitted from the sampling. We also analyzed the number of CpGs and (G + C) content in 500-bp sequences surrounding each of the 213 NotI loci (Fig. 1B). These 213 NotI loci are obviously distributed in a wide range of CpG densities (0.8–16.6 CpGs/100 bp) and (G + C) content (35.4%–79.4%), and do not merely represent a biased sampling from the genome with uniform base compositional contexts.


Figure 1
View larger version (17K):
[in this window]
[in a new window]

 
Figure 1  Characteristics of the total 6057 and the selected 213 NotI loci in the mouse genome. (A) Distribution of NotI loci by CpG island and by unique or repetitive sequence. The information regarding the total population of NotI sites in the mouse genome is summarized from our previous study (Fazzari & Greally 2004). (B) The CpG density of flanking sequences of the sampled NotI loci. CpG density and (G + C)%, calculated based on 500-bp windows centered at NotI sites, are plotted for 213 NotI sites.

 
The methylation status of the selected NotI sites define cell- and tissue-specific DNA methylation profiles

Genomic DNA was extracted from 18 types of tissues and cells of male C57BL/6 mice [brain, placenta, kidney, testis, sperm, liver, spleen, pituitary, skeletal muscle (soleus and gastrocnemius), brown adipocytes, brown preadipocytes, white adipocytes, white preadipocytes, ES cells (stem and differentiated conditions) and EG cells (stem and differentiated conditions)], and subsequently digested by NotI. Analysis of the DNA methylation status on the unbiased sample of 213 NotI loci was performed by quantitative real-time PCR. The methylation level at each locus was expressed as percentage of the proportion of the amount of undigested DNA in the NotI-treated genomic DNA sample to that in the NotI-untreated sample.

To examine if the methylation profiles of the Vi-RLGS-based NotI loci actually represent similarities in cell types as well as RLGS profile, agglomerative hierarchical clustering was applied for analysis (Fig. 2). The data are shown in a matrix format, representing methylation ratio by blue–yellow color scale. Euclidean distance (the square root of the sum of the squared differences in each dimension) was applied to reflect the range of methylation ratios. The distance was measured using the group average method (average link). The euclidean distance matrix was determined from the 213 dimensional vectors of DNA methylation ratios based on the total sample of NotI loci. The methylation profiles of tissues and cells (horizontal colored bars), consisting of various levels of methylation at the NotI sites, showed differences reflecting the phenotypic diversity among the tissues and cells. The clustering of the profiles revealed a remarkably ordered pattern, and shows that many but not all of NotI loci are tissue-dependently and differentially methylated, whereas many loci showed hypomethylated status universally in all of the sampled tissues and cells. The data also show a greater relatedness between tissues and cell types having common developmental pathways or similar phenotypes. Germ-line cells (sperm and testis) showed strong similarity in methylation profiles, creating one small cluster separate from other somatic tissues and cells. Stem ES and EG cells, both of which have comparable pluripotent phenotypes, formed a discrete cluster as well. Two skeletal muscle tissues sampled from different parts of the body (soleus and gastrocnemius) also showed a similarity compared to other tissues. The remarkable cluster within somatic tissues represents the adipose tissue group. Each of combinations of adipocytes and preadipocytes (white adipocytes and white preadipocytes/brown adipocytes and brown preadipocytes) formed discrete pairs of clusters. These clusters fell into a group based solely on similarities in patterns of methylation. The adipocyte cluster, skeletal muscle cluster and kidney, derived commonly from mesodermal lineage, were grouped into one large cluster. Pituitary, the majority of the sample represents anterior pituitary, is grouped not with brain, an ectodermal tissue, but with liver, a tissue with a similar endodermal origin. Differentiated ES and EG cells, consisting of multiple cell types, also formed a cluster. Spleen and placenta, inherently heterogeneous groups of cells, form a discrete cluster together with differentiated ES and EG cells. These results demonstrated that the clustering of methylation profiles of unbiased sampled Vi-RLGS NotI loci can reveal similarities in developmental pathways and functional phenotypes within and between stem, germ and somatic cells.


Figure 2
View larger version (24K):
[in this window]
[in a new window]

 
Figure 2  Delineation of developmental similarity by hierarchical clustering of methylation profiles of 213 NotI loci. Depicted is the methylation matrix consisting of 213 NotI loci (row) and 18 tissue or cell types (column). The top blue–yellow color scale bar refers to the degree of methylation (percent). Gray indicates missing or excluded data. A dendrogram obtained from agglomerative hierarchical clustering analysis is shown at the top of the panel. Relatedness of tissues or cells is represented by the tree whose branch lengths represent the degree of similarity.

 
The genomic context around NotI sites influences the variability of their methylation levels between cells and tissues

A total of 3834 methylation analyses demonstrate a non-uniform distribution of methylation status (Fig. 3). Whereas two-thirds (2647) of the measurements showed ratios indicating less than 10% methylation; hyper-methylation and intermediate methylation values were also observed (Fig. 3A). When analyzed in terms of CpG island status, the majority of the measurements at CpG islands (93%) show less than 10% methylation, but there also existed hypermethylated and intermediately methylated NotI loci (Fig. 3B). Non-CpG island loci show a bimodal distribution of methylation values, although intermediate methylation was also detected. Both unique and repetitive sequences also showed a range of methylation values (Fig. 3C). The results demonstrate that the diverse methylation status at NotI sites reflect information about the genomic contexts of these loci.


Figure 3
View larger version (19K):
[in this window]
[in a new window]

 
Figure 3  Distribution of DNA methylation status determined by total 3834 real-time PCR measurements. (A) The distribution of total methylation measurements shows a bimodal distribution of hypomethylation and hypermethylation, with variable intermediate methylation status. (B) Distribution of DNA methylation status by CpG island categorization. The majority of CpG island loci showed hypomethylation, whereas non-CpG island loci showed a variable distribution. (C) Distribution of DNA methylation status in unique sequences and repetitive elements. Both unique sequence loci and repetitive element loci showed diverse distribution.

 
In order to observe variability of methylation status between tissues or cells, the standard deviations of methylation values at individual loci are depicted (Fig. 4). The loci were aligned in the rank order of standard deviation of methylation status among tissues to illustrate the variability of methylation status. Whereas unique and repetitive sequences demonstrated similar distribution of standard deviation in methylation, non-CpG island loci showed higher variability in methylation status than CpG island loci. This raises the possibility that variability of methylation is affected by CpG density. We measured the number of CpGs in the 500-bp sequences surrounding each of the 213 sampled NotI loci and correlated these values with the variability of methylation status at each NotI locus (Fig. 5). We chose a 500-bp window because of a previous report stating that methylated and unmethylated CpGs of 14 rodent and human genes can be distinguished by the CpG density based on this window size (Matsuo et al. 1993). The variability of methylation status is indicated as standard deviation of methylation ratio. Figure 5 depicts the correlation between variability of DNA methylation and CpG density of the surrounding 500 bp. The result clearly shows that, in spite of their distribution in the diverse genomic contexts, the loci showing the larger differences in methylation status between tissues or cells are disproportionately distributed in less-CpG-dense regions.


Figure 4
View larger version (18K):
[in this window]
[in a new window]

 
Figure 4  Variability of methylation status of 213 NotI loci in 18 types of tissues or cells. NotI loci were arranged in the order of decreasing standard deviation in methylation status. (A) The standard deviation of methylation status of the total sample of 213 NotI loci. (B) The standard deviation of methylation status of CpG island/non-CpG island loci. Non-CpG island loci have generally a larger standard deviation of methylation status than CpG island loci. (C) The standard deviation of methylation status of unique sequences and repetitive elements.

 

Figure 5
View larger version (14K):
[in this window]
[in a new window]

 
Figure 5  Average methylation ratio, standard deviation of methylation status and CpG density of the 213 NotI loci. Legends show average methylation ratio of one NotI locus among different tissue or cell types, and the length of vertical line of each legend shows the standard deviation of methylation ratio of the locus.

 
The DNA methylation status matrix suggested that clustering of T-DMRs represents the developmental relatedness. In order to identify such T-DMRs, we extracted the minimum subset of loci required to produce the same topology of the dendrogram. We chose the loci in the rank order of standard deviation of methylation status among tissues. A subset of 92 NotI loci with standard deviation values exceeding 4.0% constituted the minimum required to recreate the tree topology. In order to understand the link between genomic context and the variability of methylation status of these loci, we extracted information about the relative positions of these T-DMRs to genes using the UCSC genome browser <http://genome.ucsc.edu>. While more than half of the T-DMRs were associated with genic regions, the majority were located within gene bodies rather than upstream or overlapping the first exon (Table 1).


View this table:
[in this window]
[in a new window]

 
Table 1  The relative position of T-DMRs responsible for the clustering to genes/transcripts
 

    Discussion
 Top
 Abstract
 Introduction
 Results
 Discussion
 Experimental procedures
 References
 
The DNA methylation profile is unique to each cell type (Shiota & Yanagimachi 2002; Shiota 2004). Tissue-specific differential RLGS spots have been identified (Shiota et al. 2002), with subsequent identification and confirmation of 150 corresponding genomic loci using Vi-RLGS (Song et al. 2005). In this study, we show that the application of Vi-RLGS to identify the genomic loci corresponding to differently methylated RLGS spots followed by hierarchical clustering analysis reveals that profiles consisting of limited number of loci alone could be used to discriminate tissues and cells on the basis of similarities in developmental pathways or functional phenotypes. We found diversity in the genomic contexts of NotI loci in these profiles, but the tissue-specific methylated loci occur disproportionately at sites flanked by low CpG density.

The disproportionate distribution of tissue-specific methylation in CpG-poor regions is consistent with previous reports about T-DMRs. CpG islands were formerly believed to maintain an unmethylated status so as to enable them to express associated genes ubiquitously (Yen et al. 1984; Antequera & Bird 1993; Cross & Bird 1995). T-DMRs have been identified at the edges of the CpG islands where the density of CpGs is lower than the centers of CpG islands (Imamura et al. 2001). This is consistent with the finding of T-DMRs in non-CpG islands, where the CpG density is lower; therefore, CpG density may be a key contributor to variability of methylation status. It has been reported that CpG-poor promoters are more prone to demethylation occurring in cancerous cells (Hatada et al. 2006). The greater flexibility in methylation states of regions less dense in CpGs may contribute to the diversity of methylation states that underlies heterogeneity of cell types.

One of the unsolved questions in typing of cells is the issue of the differentiation of white versus brown adipocytes from their precursors. It is still unknown whether brown and white adipocytes derive from the same preadipocytic precursors or whether they arise independently from undetermined mesenchymal stem cells (Rosen et al. 2000). Our finding that the close but distinguishable clusters of methylation patterns of these two types of adipocytes gives novel insights into understanding similarity in their origin and characteristics. No rigorous morphologic distinction between white and brown adipocytes in developing tissues is possible; even ultrastructural studies using electron microscopy cannot detect marked differences (Poissonnet et al. 1988). Molecular-based studies have reported several signaling factors, transcription factors and cell cycle regulators in the distinctive control pathways of white versus brown adipocyte differentiation (Puigserver et al. 1998; Cederberg et al. 2001; Picard et al. 2002; Hansen et al. 2004). Microarray analysis of gene expression reveals that white preadipocytes and brown preadipocytes have very similar gene expression profiles; even uncoupling protein 1 (UCP-1), a gene responsible for mitochondrial thermogenesis considered to be a marker for brown adipocytes (Klaus et al. 1995), can be seen in white preadipocytes and adipocytes (Boeuf et al. 2001). The present study demonstrated that each of the preadipocytes show greater similarity in DNA methylation profiles to their derived mature adipocytes rather than to the other types of preadipocyte. This is consistent with the belief that preadipocytes isolated separately from the white and brown adipose tissues are already committed to their distinct adipocyte lineages. On the other hand, both of the clusters formed a single large cluster, possibly reflecting the commonality in origins of these two adipocyte lineages relative to the other cell types tested.

The loci showing cell- or tissue-specific methylation were located not only 5'-upstream of the first exon but also in intragenic, intergenic or non-genic regions, raising the issue of the functional significance of these sites. It has also been reported in an RLGS-based study that almost half of the NotI loci that were methylated in a tissue-dependent manner were not located in 5'-promoter CpG islands (Song et al. 2005). These non-promoter T-DMRs may have more complex cis-regulatory properties rather than simple methyl-CpG-mediated gene silencing mechanisms. For example, Lorincz et al. have reported that intragenic CpG methylation forms a chromatin structure that causes reduced efficiency of RNA polymerase II elongation (Lorincz et al. 2004). It has also been demonstrated that the differentially methylated regions of the imprinted Igf2/H19 domain contribute to various levels of higher-order chromatin structure such as chromatin loops and insulator formation (Bell & Felsenfeld 2000; Hark et al. 2000; Murrell et al. 2004). Notably, there are several reports that cells having common developmental pathways have similarities in higher-order chromatin structure and spatial organization of chromosomes (Parada et al. 2004; Crawford et al. 2006). The link between DNA methylation profiles reflecting similarity in development and differentiation and these kinds of higher-order regulatory mechanisms of chromatin and chromosomes should be a logical next step in elucidating higher-order gene regulatory processes.

Numerous repetitive elements also formed tissue- or cell-type-specific DNA methylation profiles which distinguish tissue and cell types, suggesting a possible biological significance in development and differentiation. Although it is well recognized that repetitive elements are usually methylated to silence transposable activity (Yoder et al. 1997), their methylation levels can differ between some tissues (Chalitchagorn et al. 2004; Burden et al. 2005). It has been previously reported that hypomethylation and activation of intracisternal A-particle (IAP)-LTR occurs in spermatogenesis (Dupressoir & Heidmann 1996). The differential methylation of the individual transposable elements may modify their transposable activity and contribute to the generation of diversity in cellular differentiation.

The present study confirmed and extended the concept that epigenetic distance based on differential methylation of RLGS spots reflects cellular similarity and developmental origins (Shiota et al. 2002). The correlation between similarity in DNA methylation profiles and relatedness in cell types and origins might be applied for diagnosis of aberrant cells in diseases. It is well known that aberrant DNA methylation is a potential cause of tumorigenesis (Jones 2003; Ushijima 2005). As for cancer cells, several studies have reported that cells with similar origins or phenotypes show relatedness in RLGS profiles (Costello et al. 2000; Feltus et al. 2003). Epigenetic distance in DNA methylation profiles can be an indicator not only of commonality in developmental pathways but also of abnormality of cells, even in non-cancerous cellular alteration in diseases.

In summary, this study describes the DNA methylation profiles for a number of loci using Vi-RLGS, defining T-DMRs responsible for patterns of cytosine methylation that reflect tissue and cell type similarity, thereby allowing their developmental origins to be defined epigenetically. A common characteristic to the T-DMRs in diverse genomic contexts was their distinctive flanking CpG density. These data show how T-DMRs can be used to generate a novel means of specifying tissue and cell phenotypes, allowing cell lineages to be traced and functional characters (including aberrancy such as neoplasia) to be defined.


    Experimental procedures
 Top
 Abstract
 Introduction
 Results
 Discussion
 Experimental procedures
 References
 
Preparation of tissues and cells

ES and EG cells were cultured as previously described (Matsui et al. 1992; Matise et al. 2000). Spontaneous differentiation of ES and EG cells was induced by culturing in DMEM (Invitrogen, Carlsbad, CA) + 15% FCS (Biowest, Miami, FL) for 4–6 days in a monolayer culture after removal of LIF and feeder cells. These stem cells were derived from the same background strain (C57BL/6) as the male mice from which tissues were dissected. All tissues except for pituitary and sperm [brain, placenta, kidney, testis, liver, spleen, skeletal muscle (soleus and gastrocnemius), interscapular brown fat and epididymal white fat] were sampled from three to five mice. The pituitary gland was sampled from 27 mice, and sperm were obtained from cauda epididymides of 66 mice. Isolation of adipocytes and its precursors (stromal–vascular fraction) from each type of adipose tissue was performed essentially as previously described (Nechad et al. 1983). Briefly, each tissue was incubated in the presence of 0.2% (w/v) collagenase in 4 ml DMEM. After incubation, the tissue debris was removed by filtration through a 250-µm nylon screen. Cells were fractionated by the centrifugation at 700 g for 10 min. The mature adipocytes were separated by their floating to the surface of the cell suspension whereas the stromal–vascular fraction was isolated in the pellet. This separation process was repeated three times, and all the three collected samples were pooled for DNA extraction.

Preparation of genomic DNA

Genomic DNA was extracted as described previously (Ohgane et al. 1998). Briefly, cells or tissues were lysed in lysis buffer [150 mM EDTA, 10 mM Tris–HCl pH 8.0, 1% SDS, 100 µg/mL proteinase K (Merck, Darmstadt, Germany)]. The lysate was incubated at 55 °C for 20 min. Following two extractions in phenol/chloroform/isoamyl alcohol (50:49:1), the genomic DNA was precipitated in ethanol, and dissolved in 50 µL of TE buffer (10 mM Tris–HCl, 1 mM EDTA, pH 7.6).

Identification of 213 candidate informative NotI sites

Virtual image restriction landmark genomic scanning (Vi-RLGS) software, which simulates 2-D electrophoresis of DNA digested in silico by the combination of NotI-PvuII-PstI, using genome sequence information obtained from the UCSC genome browser [March 2005 (mm5) assembly, build 33 assembly by NCBI], was used to identify non-redundant NotI sites (Matsuyama et al. 2003). The redundant repetitive NotI-PstI fragments, which are plotted as aggregated large spots in the Vi-RLGS image, were removed from consideration. The Vi-RLGS simulation identified 2170 spots in the field corresponding to the observed experimental RLGS profile. The field was divided into 7 x 6 matrix areas, and the sequence information of 213 non-redundant spots was sampled from the matrix to sample 13–45 spots from each vertically divided area and 26–41 spots from each horizontally divided area (Supplementary Fig. S1 and Supplementary Table S2).

Methylation-sensitive quantitative real-time PCR

To analyze the methylation status of the 213 candidate informative NotI sites sampled from Vi-RLGS, methylation-sensitive quantitative real-time PCR was performed as previously described (Hattori et al. 2004a). Genomic DNA was digested by PstI, and the aliquot was treated subsequently with NotI. Primer sets were designed to amplify the genomic fragments including the NotI site. Ten nanograms of the genomic DNA treated with or without NotI was analyzed by real-time PCR with the primers. The amount of undigested DNA both in NotI-treated and -untreated genomic DNA was calculated by real-time PCR with SYBR Green I Universal PCR Master Mix by using ABI Prism 7500 Sequence Detection System (Applied Biosystems, Foster City, CA) according to the manufacturer's protocol. The methylation level at each locus was expressed as percentage of the proportion of the amount of undigested DNA in NotI-treated DNA compared to that in the untreated sample. The initial DNA amount in the reaction mix was normalized using primers which amplify genomic regions without NotI sites at the c-met and 5-lipo2 loci. The average of two independent assays was used for data analysis. The primers used in the analysis, methylation data and real-time PCR data are listed in Supplementary Tables S3–S5, respectively.

Data analysis

Two-dimensional agglomerative hierarchical clustering using the average linkage method was applied (Eisen et al. 1998). Euclidean distance was applied for the distance matrices. For the analysis, the publicly available program MeV was utilized <http://www.tigr.org/software/tm4/> (Saeed et al. 2003).

The genomic contexts of each of the NotI sites were retrieved from UCSC genome browser <http://genome.ucsc.edu> using the tracks within the browser. Overlap with CpG island or non-CpG island sequences was categorized using the CpG island track (Gardiner-Garden & Frommer 1987). The information of repetitive and low complexity sequence was retrieved from the RepeatMasker track (Smit, AFA, Hubley, R & Green, P. RepeatMasker Open-3.0. 1996–2004 <http://www.repeatmasker.org>) using RepBase (Jurka 2000). The relative position to genes (upstream, exon and intron) were defined as positive when the Refseq genes or mRNAs appeared within 5 kbp of a NotI site.


    Acknowledgements
 
This work was supported by the Program for Promotion of Basic Research Activities for Innovative Biosciences (PROBRAIN), Grant-in-aid for Scientific Research (15080202) and Special Coordination Funds for Promoting Science and Technology (040600000729), Ministry of Education, Culture, Sports, Science and Technology (K.S.), and Research Fellowships of the Japan Society for the Promotion of Science for Young Scientists (H.S.).


    Footnotes
 
Communicated by: Yo-ichi Nabeshima

aThese authors contributed equally to the work. Back

bPresent address: Pfizer Japan Inc., Shinjuku Bunka Quint, 3-22-7 Yoyogi, Tokyo 151-8589, Japan. Back

cPresent address: Department of Medicine (Hematology), Albert Einstein College of Medicine, Bronx, NY, 1046 USA. Back

* Correspondence: E-mail: ashiota{at}mail.ecc.u-tokyo.ac.jp


    References
 Top
 Abstract
 Introduction
 Results
 Discussion
 Experimental procedures
 References
 
Antequera, F. & Bird, A. (1993) Number of CpG islands and genes in human and mouse. Proc. Natl. Acad. Sci. USA 90, 11995–11999.[Abstract/Free Full Text]

Bell, A.C. & Felsenfeld, G. (2000) Methylation of a CTCF-dependent boundary controls imprinted expression of the Igf2 gene. Nature 405, 482–485.[CrossRef][Medline]

Bird, A.P. (1987) CpG islands as gene markers in the vertebrate nucleus. Trends Genet. 3, 342–347.[CrossRef]

Bird, A.P. & Wolffe, A.P. (1999) Methylation-induced repression—belts, braces, and chromatin. Cell 99, 451–454.[CrossRef][Medline]

Boeuf, S., Klingenspor, M., Van Hal, N.L., Schneider, T., Keijer, J. & Klaus, S. (2001) Differential gene expression in white and brown preadipocytes. Physiol. Genomics 7, 15–25.[Abstract/Free Full Text]

Burden, A.F., Manley, N.C., Clark, A.D., Gartler, S.M., Laird, C.D. & Hansen, R.S. (2005) Hemimethylation and non-CpG methylation levels in a promoter region of human LINE-1 (L1) repeated elements. J. Biol. Chem. 280, 14413–14419.[Abstract/Free Full Text]

Cederberg, A., Grønning, L.M., Ahrén, B., Taskén, K., Carlsson, P. & Enerbäck, S. (2001) FOXC2 is a winged helix gene that counteracts obesity, hypertriglyceridemia, and diet-induced insulin resistance. Cell 106, 563–573.[CrossRef][Medline]

Chalitchagorn, K., Shuangshoti, S., Hourpai, N., Kongruttanachok, N., Tangkijvanich, P., Thong-ngam, D., Voravud, N., Sriuranpong, V. & Mutirangura, A. (2004) Distinctive pattern of LINE-1 methylation level in normal tissues and the association with carcinogenesis. Oncogene 23, 8841–8846.[CrossRef][Medline]

Cho, J.H., Kimura, H., Minami, T., Ohgane, J., Hattori, N., Tanaka, S. & Shiota, K. (2001) DNA methylation regulates placental lactogen I gene expression. Endocrinology 142, 3389–3396.[Abstract/Free Full Text]

Costello, J.F., Fruhwald, M.C., Smiraglia, D.J. et al. (2000) Aberrant CpG-island methylation has non-random and tumour-type-specific patterns. Nat. Genet. 24, 132–138.[CrossRef][Medline]

Crawford, G.E., Holt, I.E., Whittle, J., et al. (2006) Genome-wide mapping of DNase hypersensitive sites using massively parallel signature sequencing (MPSS). Genome Res. 16, 123–131.[Abstract/Free Full Text]

Cross, S.H. & Bird, A.P. (1995) CpG islands and genes. Curr. Opin. Genet. Dev. 5, 309–314.[CrossRef][Medline]

Dupressoir, A. & Heidmann, T. (1996) Germ line-specific expression of intracisternal A-particle retrotransposons in transgenic mice. Mol. Cell. Biol. 16, 4495–4503.[Abstract/Free Full Text]

Eisen, M.B., Spellman, P.T., Brown, P.O. & Botstein, D. (1998) Cluster analysis and display of genome-wide expression patterns. Proc. Natl. Acad. Sci. USA 95, 14863–14868.[Abstract/Free Full Text]

Fazzari, M.J. & Greally, J.M. (2004) Epigenomics: beyond CpG islands. Nat. Rev. Genet. 5, 446–455.[CrossRef][Medline]

Feltus, F.A., Lee, E.K., Costello, J.F., Plass, C. & Vertino, P.M. (2003) Predicting aberrant CpG island methylation. Proc. Natl. Acad. Sci. USA 100, 12253–12258.[Abstract/Free Full Text]

Futscher, B.W., Oshiro, M.M., Wozniak, R.J., Holtan, N., Hanigan, C.L., Duan, H. & Domann, F.E. (2002) Role for DNA methylation in the control of cell type specific maspin expression. Nat. Genet. 31, 175–179.[CrossRef][Medline]

Gardiner-Garden, M. & Frommer, M. (1987) CpG islands in vertebrate genomes. J. Mol. Biol. 196, 261–282.[CrossRef][Medline]

Hansen, J.B., Jørgensen, C., Petersen, R.K., Hallenborg, P., De Matteis, R., Bøye, H.A., Petrovic, N., Enerbäck, S., Nedergaard, J., Cinti, S., te Riele, H. & Kristiansen, K. (2004) Retinoblastoma protein functions as a molecular switch determining white versus brown adipocyte differentiation. Proc. Natl. Acad. Sci. USA 101, 4112–4117.[Abstract/Free Full Text]

Hark, A.T., Schoenherr, C.J., Katz, D.J., Ingram, R.S., Levorse, J.M. & Tilghman, S.M. (2000) CTCF mediates methylation-sensitive enhancer-blocking activity at the H19/Igf2 locus. Nature 405, 486–489.[CrossRef][Medline]

Hatada, I., Fukasawa, M., Kimura, M., Morita, S., Yamada, K., Yoshikawa, T., Yamanaka, S., Endo, C., Sakurada, A., Sato, M., Kondo, T., Horii, A., Ushijima, T. & Sasaki, H. (2006) Genome-wide profiling of promoter methylation in human. Oncogene 25, 3059–3064.[CrossRef][Medline]

Hattori, N., Abe, T., Hattori, N., Suzuki, M., Matsuyama, T., Yoshida, S., Li, E. & Shiota, K. (2004a) Preference of DNA methyltransferases for CpG islands in mouse embryonic stem cells. Genome Res. 14, 1733–1740.[Abstract/Free Full Text]

Hattori, N., Nishino, K., Ko, Y.G., Hattori, N., Ohgane, J., Tanaka, S. & Shiota, K. (2004b) Epigenetic control of mouse Oct-4 gene expression in embryonic stem cells and trophoblast stem cells. J. Biol. Chem. 279, 17063–17069.[Abstract/Free Full Text]

Imamura, T., Ohgane, J., Ito, S., Ogawa, T., Hattori, N., Tanaka, S. & Shiota, K. (2001) CpG island of rat sphingosine kinase-1 gene: tissue-dependent DNA methylation status and multiple alternative first exons. Genomics 76, 117–125.[CrossRef][Medline]

Jones, P.A. (2003) Epigenetics in carcinogenesis and cancer prevention. Ann. NY Acad. Sci. 983, 213–219.[CrossRef][Medline]

Jones, P.A. & Takai, D. (2001) The role of DNA methylation in mammalian epigenetics. Science 293, 1068–1070.[Abstract/Free Full Text]

Jurka, J. (2000) Repbase update: a database and an electronic journal of repetitive elements. Trends Genet. 16, 418–420.[CrossRef][Medline]

Klaus, S., Ely, M., Encke, D. & Heldmaier, G. (1995) Functional assessment of white and brown adipocyte development and energy metabolism in cell culture. Dissociation of terminal differentiation and thermogenesis in brown adipocytes. J. Cell Sci. 108, 3171–3180.[Abstract]

Ko, Y.G., Nishino, K., Hattori, N., Arai, Y., Tanaka, S. & Shiota, K. (2005) Stage-by-stage change in DNA methylation status of Dnmt1 locus during mouse early development. J. Biol. Chem. 280, 9627–9634.[Abstract/Free Full Text]

Larsen, F., Gundersen, G., Lopez, R. & Prydz, H. (1992) CpG islands as gene markers in the human genome. Genomics 13, 1095–1107.[CrossRef][Medline]

Lindsay, S. & Bird, A.P. (1987) Use of restriction enzymes to detect potential gene sequences in mammalian DNA. Nature 327, 336–338.[CrossRef][Medline]

Lorincz, M.C., Dickerson, D.R., Schmitt, M. & Groudine, M. (2004) Intragenic DNA methylation alters chromatin structure and elongation efficiency in mammalian cells. Nat. Struct. Mol. Biol. 11, 1068–1075.[CrossRef][Medline]

Matise, M.P., Auerbach, W. & Joyner, A.L. (2000) Production of targeted embryonic stem cell clones. In: Gene Targeting: A Practical Approach (ed. A.L. Joyner), pp. 102–132, New York: Oxford University Press.

Matsui, Y., Zsebo, K. & Hogan, B.L. (1992) Derivation of pluripotential embryonic stem cells from murine primordial germ cells in culture. Cell 70, 841–847.[CrossRef][Medline]

Matsuo, K., Clay, O., Takahashi, T., Silke, J. & Schaffner, W. (1993) Evidence for erosion of mouse CpG islands during mammalian evolution. Somat. Cell Mol. Genet. 19, 543–555.[CrossRef][Medline]

Matsuyama, T., Kimura, M.T., Koike, K., Abe, T., Nakano, T., Asami, T., Ebisuzaki, T., Held, W.A., Yoshida, S. & Nagase, H. (2003) Global methylation screening in the Arabidopsis thaliana and Mus musculus genome: applications of virtual image restriction landmark genomic scanning (Vi-RLGS). Nucleic Acids Res. 31, 4490–4496.[Abstract/Free Full Text]

Murrell, A., Heeson, S. & Reik, W. (2004) Interaction between differentially methylated regions partitions the imprinted genes Igf2 and H19 into parent-specific chromatin loops. Nat. Genet. 36, 889–893.[CrossRef][Medline]

Nechad, M., Kuusela, P., Carneheim, C., Bjorntorp, P., Nedergaard, J. & Cannon, B. (1983) Development of brown fat cells in monolayer culture. I. Morphological and biochemical distinction from white fat cells in culture. Exp. Cell Res. 149, 105–118.[CrossRef][Medline]

Newell-Price, J., King, P. & Clark, A.J. (2001) The CpG island promoter of the human proopiomelanocortin gene is methylated in nonexpressing normal tissue and tumors and represses expression. Mol. Endocrinol. 15, 338–348.[Abstract/Free Full Text]

Nishino, K., Hattori, N., Tanaka, S. & Shiota, K. (2004) DNA methylation-mediated control of Sry gene expression in mouse gonadal development. J. Biol. Chem. 279, 22306–22313.[Abstract/Free Full Text]

Ohgane, J., Aikawa, J., Ogura, A., Hattori, N., Ogawa, T. & Shiota, K. (1998) Analysis of CpG islands of trophoblast giant cells by restriction landmark genomic scanning. Dev. Genet. 22, 132–140.[CrossRef][Medline]

Pao, M.M., Tsutsumi, M., Liang, G., Uzvolgyi, E., Gonzales, F.A. & Jones, P.A. (2001) The endothelin receptor B (EDNRB) promoter displays heterogeneous, site specific methylation patterns in normal and tumor cells. Hum. Mol. Genet. 10, 903–910.[Abstract/Free Full Text]

Parada, L.A., McQueen, P.G. & Misteli, T. (2004) Tissue-specific spatial organization of genomes. Genome Biol. 5, R44.[CrossRef][Medline]

Picard, F., Géhin, M., Annicotte, J., Rocchi, S., Champy, M.F., O'Malley, B.W., Chambon, P. & Auwerx, J. (2002) SRC-1 and TIF2 control energy balance between white and brown adipose tissues. Cell 111, 931–941.[CrossRef][Medline]

Poissonnet, C.M., LaVelle, M. & Burdi, A.R. (1988) Growth and development of adipose tissue. J. Pediatr. 113, 1–9.[CrossRef][Medline]

Puigserver, P., Wu, Z., Park, C.W., Graves, R., Wright, M. & Spiegelman, B.M. (1998) A cold-inducible coactivator of nuclear receptors linked to adaptive thermogenesis. Cell 92, 829–839.[CrossRef][Medline]

Rosen, E.D., Walkey, C.J., Puigserver, P. & Spiegelman, B.M. (2000) Transcriptional regulation of adipogenesis. Genes Dev. 14, 1293–1307.[Free Full Text]

Saeed, A.I., Sharov, V., White, J. et al. (2003) TM4: a free, open-source system for microarray data management and analysis. Biotechniques 34, 374–378.[Medline]

Shiota, K. (2004) DNA methylation profiles of CpG islands for cellular differentiation and development in mammals. Cytogenet. Genome Res. 105, 325–334.[CrossRef][Medline]

Shiota, K., Kogo, Y., Ohgane, J., Imamura, T., Urano, A., Nishino, K., Tanaka, S. & Hattori, N. (2002) Epigenetic marks by DNA methylation specific to stem, germ and somatic cells in mice. Genes Cells 7, 961–969.[Abstract]

Shiota, K. & Yanagimachi, R. (2002) Epigenetics by DNA methylation for development of normal and cloned animals. Differentiation 69, 162–166.[CrossRef][Medline]

Song, F., Smith, J.F., Kimura, M.T., Morrow, A.D., Matsuyama, T., Nagase, H. & Held, W.A. (2005) Association of tissue-specific differentially methylated regions (TDMs) with differential gene expression. Proc. Natl. Acad. Sci. USA 102, 3336–3341.[Abstract/Free Full Text]

Ushijima, T. (2005) Detection and interpretation of altered methylation patterns in cancer cells. Nat. Rev. Cancer 5, 223–231.[CrossRef][Medline]

Yen, P.H., Patel, P., Chinault, A.C., Mohandas, T. & Shapiro, L.J. (1984) Differential methylation of hypoxanthine phosphoribosyltransferase genes on active and inactive human X chromosomes. Proc. Natl. Acad. Sci. USA 81, 1759–1763.[Abstract/Free Full Text]

Yoder, J.A., Walsh, C.P. & Bestor, T.H. (1997) Cytosine methylation and the ecology of intragenomic parasites. Trends Genet. 13, 335–340.[CrossRef][Medline]

Received: 20 March 2007
Accepted: 28 June 2007




This article has been cited by other articles:


Home page
Proc. Natl. Acad. Sci. USAHome page
M. Oda, S. Tanaka, Y. Yamazaki, H. Ohta, M. Iwatani, M. Suzuki, J. Ohgane, N. Hattori, R. Yanagimachi, T. Wakayama, et al.
From the Cover: Establishment of trophoblast stem cell lines from somatic cell nuclear-transferred embryos
PNAS, September 22, 2009; 106(38): 16293 - 16297.
[Abstract] [Full Text] [PDF]


Home page
Hum Mol GenetHome page
R. A. Waterland, R. Kellermayer, M.-T. Rached, N. Tatevian, M. V. Gomes, J. Zhang, L. Zhang, A. Chakravarty, W. Zhu, E. Laritsky, et al.
Epigenomic profiling indicates a role for DNA methylation in early postnatal liver development
Hum. Mol. Genet., August 15, 2009; 18(16): 3026 - 3038.
[Abstract] [Full Text] [PDF]


Home page
Genome ResHome page
S. Yagi, K. Hirabayashi, S. Sato, W. Li, Y. Takahashi, T. Hirakawa, G. Wu, N. Hattori, N. Hattori, J. Ohgane, et al.
DNA methylation profile of tissue-dependent and differentially methylated regions (T-DMRs) in mouse promoter regions demonstrating tissue-specific gene expression
Genome Res., December 1, 2008; 18(12): 1969 - 1978.
[Abstract] [Full Text] [PDF]


Home page
J BiochemHome page
K. Yamagata
Capturing Epigenetic Dynamics During Pre-implantation Development Using Live Cell Imaging
J. Biochem., March 1, 2008; 143(3): 279 - 286.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Sakamoto, H.
Right arrow Articles by Shiota, K.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Sakamoto, H.
Right arrow Articles by Shiota, K.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE ADVANCED SEARCH TABLE OF CONTENTS