Human Molecular Genetics, 2003, Vol. 12, No. 13 1569-1578
DOI: 10.1093/hmg/ddg167
© 2003 Oxford University Press
Familial aggregation of abnormal methylation of parental alleles at the IGF2/H19 and IGF2R differentially methylated regions
1Fels Institute for Cancer Research and Molecular Biology, Temple University School of Medicine, 3307 North Broad Street, Philadelphia, PA 19140, USA, 2Eccles Institute of Human Genetics, University of Utah, 15 N 2030 E, Salt Lake City, UT 84112, USA, 3Temple University School of Medicine, 3400 North Broad Street, Philadelphia, PA 19140, USA, 4Bucknell University, 701 Moore Avenue, Lewisburg, PA 17837, USA, 5Phillips Academy, 180 Main Street, Andover, MA 01810, USA and 6Department of Pathology and Laboratory Medicine, Temple University School of Medicine, 3307 North Broad Street, Philadelphia, PA 19140, USA
Received March 5, 2003; Accepted April 28, 2003
| ABSTRACT |
|---|
|
|
|---|
Loss of imprinting (LOI) has been observed in many types of human tumors and may be a predisposing event in some colon cancers. LOI is strongly associated with alteration of normal DNA methylation patterns in differentially methylated regions (DMRs) of affected loci but it is not known whether LOI is caused by stochastic, environmental or genetic factors. We have developed a simple, quantitative assay for measurement of allelic methylation ratios based on methylation-sensitive restriction endonuclease digestion of genomic DNA and hot-stop PCR. We examined allelic methylation ratios at DMRs within the IGF2/H19- and IGF2R-loci in a panel of 48 three-generation families. We observed familial clustering of individuals with abnormal methylation ratios at the IGF2/H19 DMR, as well as stability of this trait over a period of nearly two decades, consistent with the possibility that constitutional LOI at this locus is due largely to genetic factors. At the IGF2R DMR, we observed more variability in the allelic methylation ratios over time but also observed familial clustering of abnormal methylation ratios. Overall, our observations at IGF2R suggest that shared genetic factors are responsible for a major fraction of inter-individual variability in parental origin-dependent epigenetic modifications. However, temporal changes also occur in isolated cases, as well as within multiple individuals in the same family, indicating that environmental factors may also play a role.
| INTRODUCTION |
|---|
|
|
|---|
Genomic imprinting is a differential epigenetic modification of chromosomes that may result in the parent-of-origin specific transcription of some genes. DNA methylation (reviewed in 1), as well as various histone modifications (methylation, acetylation, phosphorylation and ubiquitination) (reviewed in 2,3), are likely to be important components of the somatically heritable imprint.
The great majority of transcriptionally imprinted genes show allelic differences in DNA methylation at so-called differentially methylated regions (DMRs). Loss of imprinting (LOI), or relaxation of imprinting (4,5), in which a gene transcribed normally from only one parent's allele is transcribed from both alleles, has been observed in both pediatric and adult cancers (reviewed in 6), as well as normal tissues (7). There is a strong association between LOI and altered DNA methylation patterns at DMRs.
A number of studies suggest that there is inter-individual variability in allele-specific expression of some imprinted genes. For example, the IGF2R gene, although differentially methylated on maternal and paternal alleles, is transcribed from both alleles in the vast majority of the population (8). However, a subset of individuals express predominantly or only the maternal allele (9). The WT1 gene has been reported to be transcribed from both alleles in the majority of placentae but placentae from some individuals express only the maternal allele (10). Further evidence for inter-individual variability in imprinting comes from the observation of significant biallelic expression of the insulin-like growth factor II gene (IGF2) in about 10% of normal individuals (11). LOI at IGF2 also occurs in tumor tissue of 31% (12) to 44% (7) of sporadic colon cancer cases. Surprisingly, LOI is not thought to be a tumor-specific event in these cases, but is likely to be constitutional because both alleles of IGF2 are expressed in the normal colonic mucosa of 90% of these patients and also in the lymphocytes of one-third of these same individuals (7).
In the present study, we have attempted to determine whether stochastic, environmental or genetic factors are responsible for the potentially large population of normal individuals (11) who have altered genome imprints. We have developed an assay based on methylation-sensitive restriction endonuclease digestion of genomic DNA and hot-stop PCR. We used this assay to examine inter-individual variability in methylation at CpG sites within the IGF2/H19 DMR and the IGF2R DMR among a population of 680 individuals from 48 three-generation families. We also examined intra-individual variability in this phenotype by comparing peripheral blood DNA samples taken nearly two decades apart from a subpopulation of the same individuals as a way of detecting the influence of environmental or stochastic factors on parental origin-dependent epigenetic modifications over this time period.
| RESULTS |
|---|
|
|
|---|
Familial aggregation of individuals with abnormal methylation ratio at the IGF2/H19-DMR
Because we wished to screen a population of moderate size (680 individuals) for LOI using archival material, we developed an assay based on the correlation between parent of origin-specific transcription and methylation of CpG sites within DMRs (Fig. 1). Briefly, CpG sites within DMRs on the paternal allele at the IGF2/H19 locus are normally methylated, while those on the maternal allele are normally unmethylated (13). There are several DMRs at this locus but the most consistent observations indicating a role in the control of transcription of the IGF2 and H19 genes involve a CpG island located in a 5 kb region centromeric to the H19 gene (14). This is also the region that contains seven different binding sites for the CTCF protein (15) and the methylation status of the sixth binding site was found to be most consistently associated with the transcriptional status of both IGF2 and H19 (16). We used a single nucleotide polymorphism (a C/T polymorphism within a CfoI site, Fig. 1A) and pedigree analysis to identify maternal and paternal alleles of heterozygous individuals and a methylation-sensitive restriction endonuclease (MluI and, in selected cases, MaeII) to determine the methylation status of specific CpG sites within the DMR. If all paternal alleles are methylated and all maternal alleles are unmethylated at these sites in a sample of genomic DNA, then all maternal alleles should be cleaved by MluI while all paternal alleles will remain uncleaved. Amplification of the region by PCR using primers that flank the MluI site should amplify only paternal alleles (identified by post-PCR cleavage with CfoI). Amplification of maternal alleles indicates resistance to cleavage by MluI. This may occur as a result of methylation of the CpG site within the MluI recognition sequence (the principle upon which the assay is based), mutation of the MluI site or technical artifact. The latter two possibilities may be distinguished from the first by DNA sequencing, assay reproducibility and use of additional methylation-sensitive restriction endonucleases.
|
We identified 318 individuals (46.8%) who were heterozygous at the CfoI polymorphism (Fig. 1A) within the IGF2/H19 DMR and were able to determine unequivocally the parental origins of the alleles in 163 cases. Because previous studies have indicated that LOI is a quantitative trait (7), i.e. expression of the maternal IGF2 allele is variable and may not occur to the same extent in all individuals identified as having LOI, we calculated the ratio between the DNA methylation levels on maternal and paternal (M/P) alleles as an indicator of imprinting status (see Materials and Methods). A ratio of zero corresponds to exclusive methylation of the MluI site on the paternal allele, while a ratio of one signifies methylation of this site on an equal number of maternal and paternal alleles.
The distribution of individual M/P methylation ratios (Fig. 2A) shows that the great majority of individuals have a methylated paternal allele, with no detectable methylation of the maternal allele. (Because none of the 163 individuals for whom we could determine parental origin of alleles by pedigree analysis had a maternal allele that was more methylated than the paternal allele, we have assumed that the less methylated allele is maternal in the remainder of cases compiled in Fig. 2A.) However, 21 individuals (6.6%) have M/P ratios greater than 0.1. We obtained similar ratios after digestion of genomic DNA with either MluI or MaeII, indicating that methylation of the maternal allele in these individuals is not restricted to the MluI site. We excluded the possibility that a subpopulation of lymphocytes from these individuals carried a mutation at the MluI or MaeII sites by DNA sequencing of the amplified maternal allele (Fig. 2C). An additional 28 individuals exhibit lower, but reproducible, levels of methylation on the maternal allele (M/P ratios between 0.05 and 0.1). While some of these individuals occur as isolated cases in some families, a significant proportion are found within families in which other members also have methylation of maternal alleles. Nine out of the 21 individuals (42.8%) with M/P ratios over 0.1 are found within two families (Fig. 3). In one of these families, six of 10 informative individuals have M/P ratios over 0.1 (for H0: six abnormal individuals occur among 10 informative individuals by chance, given 6.6% population frequency,
2 with Yates correction=20.2, P<0.001). When a lower level of maternal allele methylation is considered (an M/P ratio over 0.05), the familial aggregation of individuals is even more pronounced, with 34 out of 49 (69.4%) individuals having at least one relative with this phenotype. Moreover, not every individual in these families is informative, and the true number of cases in which only a single family member is affected may be overestimated. We also note that 10 of the 15 individuals who have no relatives with methylation of the maternal allele are in the first generation of their pedigrees. In these cases, no sibling, and only a single offspring (either the mother or the father of the third generation), is available for analysis.
|
|
Familial aggregation of individuals with altered differential methylation at the IGF2R DMR
Although the human IGF2R gene is normally expressed from both alleles, differential methylation of maternal and paternal alleles is maintained in the human, as it is in the mouse (17), and some humans do appear to have transcriptional imprinting of IGF2R (9). We designed an assay (Fig. 1B) to examine relative levels of methylation at CpG sites on maternal and paternal alleles within the DMR of IGF2R (located in the second intron) using the same principles upon which the IGF2/H19 assay was based. We used a C/T polymorphism within an MspI site to distinguish heterozygous individuals, pedigree analysis to determine maternal and paternal alleles and cleavage of genomic DNA by a methylation-sensitive restriction endonuclease followed by PCR using primers that flank the cleavage site to distinguish methylation of CpG sites on maternal and paternal alleles.
We observed 165 individuals (24.3%) who were informative for the MspI polymorphism found within the IGF2R DMR and we identified maternal and paternal alleles for 112 individuals by pedigree analysis. Although maternal and paternal alleles were distinguishable in almost all individuals on the basis of expected preferential methylation of the maternal allele, the distribution of paternal/maternal (P/M) methylation ratios at IGF2R DMR shows that, unlike the IGF2/H19 DMR, most individuals have a measurable level of methylation at the CpG within the XhoI cleavage site on the paternal allele (P/M>0.1, Fig. 4). P/M ratios were similar for both XhoI and NotI restriction endonuclease digestion of genomic DNA (correlation coefficient for paired measurements r=0.715) indicating, again, that the methylation of CpG sites within the paternal DMR is not restricted to a single site.
|
Because the distribution in Figure 4A is two-tailed (unlike the distribution of P/M methylation ratios at IGF2/H19 DMR), we examined the distribution of individuals from each tail among the 48 families. Because the vast majority of individuals have some methylation of the paternal allele, this modification seems unlikely to inhibit transcription of the normally transcribed paternal allele. However, it is possible that individuals who have little or no methylation of the paternal allele are those individuals in which the maternal allele is predominately or exclusively transcribed, as reported by Xu et al. (9). We observed 19 individuals (11.5%) with P/M ratios under 0.1. Fourteen of these 19 individuals are distributed among only five families, with five individuals being found in a single family (Fig. 5) (for H0: five abnormal individuals occur among seven informative by chance, given 11.5% population frequency,
2 with Yates correction=10.5, P<0.01).
|
We also identified five individuals in which the CpG site being examined was methylated on approximately equal numbers of maternal and paternal alleles (P/M>0.8). Two of these individuals are siblings within the same family (pedigree not shown).
Familial aggregation of individuals with abnormal methylation of IGF2/H19 and IGF2R DMRs; genetics or environment?
The familial clustering of individuals with abnormal allelic methylation ratios may be explained by the action of shared environmental factors or shared genetic factors. We have attempted to distinguish between these possibilities by genetic linkage analysis in those families in which multiple members are affected and by examining the temporal stability of the allelic methylation ratio over a period of nearly two decades.
Genetic linkage analysis of the abnormal methylation ratio phenotype
The methylation of CpG sites within DMRs is one component of the overall chromatin structure of imprinted regions of chromosomes. Abnormal methylation of normally unmethylated CpG sites could result from some aspect of the structure of the particular chromosome in which it is found (an effect in cis) or result from the action of modifying genetic factors that do not reside in the vicinity of the DMR (an effect in trans). We may gain supporting evidence for the involvement of cis-acting factors by determining whether family members who share abnormal allelic methylation ratios also share the same chromosome haplotype in the vicinity of the DMR.
The assays used in this study identify abnormal methylation of the maternal IGF2/H19 DMR and abnormal methylation of the paternal IGF2R DMR. Therefore, we determined which maternal chromosome 11p15.5 was inherited by individuals in families in which more than one individual had an abnormal methylation ratio at the IGF2/H19 DMR and which paternal chromosome 6q26 was inherited by individuals in families in which more than one individual had an abnormal methylation ratio at the IGF2R DMR (see Materials and Methods).
We identified only one family (data not shown) in which all individuals with abnormal IGF2/H19 DMR methylation ratios shared the same maternal chromosome 11p15.5 haplotype. Cis-acting genetic factors may be involved in the abnormal methylation of maternal IGF2/H19 alleles in this family. However, some normal individuals in this family inherited the same maternal chromosome 11p15.5 as did those individuals with abnormal methylation ratio so that there were no families in which there was perfect concordance between methylation ratio phenotype and chromosome 11p15.5 haplotype. There were no families in which all individuals with abnormal methylation ratios at the IGF2R DMR shared the same paternal chromosome 6q26 haplotype.
We have not attempted to perform linkage analysis of candidate trans-acting genetic factors because the number of potential candidate genes is large and the likelihood of detecting spurious associations within any one of the families is high. However, the chromosome haplotype analyses we have performed indicate that if familial clustering of the abnormal methylation ratio trait is due to genetic factors, these factors must often act in trans.
Temporal stability of allelic methylation ratio
In the absence of a hypothesis on the specific identity of an environmental agent that might influence abnormal methylation of DMRs (and associated hypotheses on how and when such an agent might act), we sought to obtain evidence on whether any environmental factor might be involved in the generation of abnormal methylation ratios by comparing the allelic methylation ratios of the same individual at two different times. This approach is based on the premise that, if an environmental factor causes an alteration of the allelic methylation ratio, then one predicts a change in this phenotype between the time of establishment of the original, normal ratio and the time at which the environmental factor exerts its effect.
To test the stability of the allelic methylation ratio over time, two sets of peripheral blood DNA samples taken from the same individuals (the first taken between 1982 and 1985 and the second in 2001) were analyzed. Both old and new samples were available for 133 individuals informative for the CfoI polymorphism at IGF2/H19 DMR and for 75 individuals informative for the MspI polymorphism at IGF2R DMR.
At the IGF2/H19 DMR, only three individuals sampled at both times (individual nos 12, 100 and 125 in Fig. 6A) had M/P allelic methylation ratios of greater than 0.1 and these individuals had a ratio of greater than 0.1 at both sampling times. In fact, there was no case in which the methylation ratio at the IGF2/H19 DMR changed substantially between the two sampling times (Fig. 6A). Because the level of methylation on the maternal allele at the IGF2/H19 DMR locus was very close to zero for the majority of individuals in both old and new DNA samples, we did not attempt to assign statistical significance to small differences between individual paired results. For the data set, as a whole, a paired t-test indicates that there is no significant difference between the two time points (P=0.1051).
|
At the IGF2R DMR, we observed eight cases in which a change in methylation ratio [only changes greater than three times the level of assay reproducibility (±0.067, see Materials and Methods) were considered as significant for this purpose] moved an individual from a normal to an abnormal category or vice versa. One individual (individual no. 43 in Fig. 6B) had a P/M methylation ratio of 0.2 in the first sample and this ratio was nearly one in the second sample. Two unrelated individuals (nos 32 and 38 in Fig. 6B) had P/M methylation ratios of less than 0.1 in the first sample and ratios of 0.26 and 0.37, respectively, in the second sample. Five individuals (individual nos 52, 54, 55, 57, and 61 in Fig. 6B) exhibited potentially significant changes in the opposite direction, i.e. attaining a P/M methylation ratio of less than 0.1 in the second sample. Three of these individuals are from the same family (a mother and two offspring, pedigree not shown), suggesting a common environmental effect on this phenotype or a genetic predisposition to the observed change.
Overall, we observed very similar results for the two sets of measurements (correlation coefficient for paired measurements r=0.7821) and the slope of the linear regression line (P/M methylation ratio at time 1 versus P/M methylation ratio at time 2) is 0.848±0.079, near the expected value of 1.0 if no change in methylation ratio occurred in the interval between collection of the two samples. These data indicate that, although changes in allelic methylation ratio can be observed, this phenotype is temporally stable in the majority of the population.
Bisulfite sequencing analysis of individuals with abnormal pattern of DNA methylation at IGF2/H19-DMR
DNA sequence analysis of bisulfite treated DNA (18,19) has become the standard technique for determination of the methylation status of CpG sites in imprinted regions. The technique has the advantage that all CpG sites are amenable to analysis but has the disadvantage that it is labor-intensive and relatively expensive to use as a screening tool. In addition, the number of individual DNA molecules that must be sequenced to demonstrate a significant difference between individuals with a M/P allelic ratio of zero and individuals with a M/P allelic ratio of 0.1, for example, is much greater than the 1020 molecules normally assayed in this procedure. These considerations led us to attempt this method only in those families in which we observed abnormal methylation ratios using the methylation-sensitive restriction endonuclease-based method. As expected, we observed DNA molecules with an abnormal mixed methylation pattern for many of these individuals (Fig. 7). However, we also analyzed a series of three single nucleotide polymorphisms present in the sequenced region of the individual shown in Figure 7, as well as additional abnormal and informative individuals, in order to confirm the parental origin of the abnormal alleles. We found that the majority of these mixed methylation pattern molecules result from maternal/paternal heteroduplex products that are created during the PCR step of the assay and then repaired, randomly, by the bacterial mismatch repair system during the cloning procedure. We confirmed the potential for artifactual origin of molecules with a mixed methylation pattern by mixing the DNAs of individuals who are homozygous for different SNP haplotypes and detecting individual molecules with both mixed SNP haplotypes and mixed maternal/paternal methylation patterns. Because we were unable to establish conditions in which such heteroduplex molecules could be completely eliminated [including techniques designed to maintain favorable primer : template ratios (20) and techniques for eliminating heteroduplexes using enzymatic digestion with T7 endonuclease I (21); data not shown], we did not pursue this method of analysis further.
|
| DISCUSSION |
|---|
|
|
|---|
Characterization of variability in imprinting phenotype using allele-specific methylation assays
We have examined the relative level of methylation of CpG sites on maternal and paternal alleles within differentially methylated regions of the IGF2/H19 and IGF2R loci as a surrogate measure of imprinting. The distributions of these quantitative measures among informative individuals indicates that approximately 6.5% of the population has a significant fraction (>10%) of peripheral blood lymphocytes that have methylation of CpG sites on the maternal IGF2/H19 allele. We note that this estimate is similar to the fraction of the population designated as having LOI (10%) after measuring transcription of maternal and paternal IGF2 alleles (11). The distribution of P/M allelic methylation at IGF2R is more variable than that at IGF2/H19. Approximately 11% of the population may be described as having abnormally low levels of methylation on the paternal allele (P/M allele ratios of <0.1) while a smaller fraction of individuals (about 3%) have P/M allelic methylation ratios greater than 0.8. Strong familial clustering of individuals from the tails of both distributions is observed (Figs 3 and 5), indicating that shared genetic factors or shared environmental factors are responsible for the majority of abnormal allelic methylation ratios (data summarized in Table 1).
|
We attempted to obtain evidence that allelic methylation ratios may be affected by environmental factors by examining samples taken from the same individuals nearly 20 years apart. We did not detect significant changes at the IGF2/H19 DMR in any of the 133 individuals examined. Although we cannot eliminate the possibility that shared environmental factors may have acted before the first peripheral blood lymphocyte sample was taken, we suggest that familial clustering of abnormal methylation ratios at the IGF2/H19 DMR is most likely due to shared genetic factors. At the IGF2R DMR, significant changes in methylation ratio did occur in eight of 75 informative individuals. Three of these changed individuals are from the same family, indicating that this acquired phenotype can also cluster in families and such familial clustering may be due to shared environmental factors or a genetic predisposition to temporal or environmentally induced changes.
We note that the restriction endonuclease/PCR-based assay used in our experiments provides a ratio of the number of maternal/paternal alleles at which a particular CpG site is methylated rather than the absolute number of DNA molecules on which a CpG site is methylated/unmethylated. For example, a sample in which neither the maternal nor the paternal IGF2/H19 DMR was methylated in 90% of cells but in which the maternal IGF2/H19 DMR was methylated in 10% of the remaining 10% of cells would give the same ratio as a sample in which the paternal DMR was methylated in all cells and the maternal DMR was methylated in 10% of all cells. These two possibilities may be distinguished by bisulfite sequencing because the vast majority of DNA molecules, be they paternal or maternal, will have the unmethylated pattern in the first case. We did not encounter such a situation (or the reverse) in any of the 10 abnormal or four normal samples that we examined by bisulfite sequencing. Nearly equal numbers of molecules with methylated and unmethylated patterns were found (disregarding molecules arising from repair of heteroduplexes). Corroborating evidence in favor of the second possibility was obtained by real-time PCR (see Materials and Methods) on several additional individuals, comparing methylation sensitive restriction endonuclease cleaved DNA with uncleaved DNA (data not shown).
Bisulfite sequencing assay versus methyl-sensitive restriction endonuclease assay
Bisulfite sequencing has become the standard method for single-molecule and single-base resolution analysis of DNA methylation (18,19). The main advantage of this method is the ability to examine the methylation of many CpG sites on a single DNA molecule. However, this method has some limitations (reviewed in 22,23), especially in situations in which inter-cellular or inter-individual variability itself is of interest. The most obvious disadvantage is the practical limitation on the number of chromosomes that can be analyzed. Sequencing of 20 DNA molecules from each individual, for example, yields an estimate of epigenetic variability based on a sample of only 10 cells.
An additional problem is the generation of heteroduplex artifacts during the amplification of mixed templates using universal primers. After bisulfite treatment a high diversity of template molecules will occur as the result of different genomic DNA methylation patterns between the parental chromosomes as well as differences in methylation of specific sites between cells. The presence of native single-nucleotide polymorphisms (SNPs) and any errors induced by the DNA polymerase will increase further the potential for heteroduplexes. When such a heteroduplex molecule is cloned, the host's mismatch repair system can convert it into a single hybrid sequence (reviewed in 24). As the repair enzymes cannot identify a parent strand, either strand is chosen as a template for the synthesis of the repaired strand. The repaired sequence will be a mosaic, i.e. composed of portions of the two parent heterologous sequences. Several protocols to reduce the incidence of mosaic molecules have been described (20,21). Such protocols rely on maintaining a primer : template ratio that favors synthesis of new molecules over hybridization of old strands, or post-synthesis enzymatic cleavage of heteroduplexes. However, it is difficult to eliminate completely the cloning of heteroduplex molecules and the utility of such data depends on the significance of these molecules to the interpretation of results. A mosaic pattern may be misinterpreted as abnormal hypermethylation or hypomethylation of some portion of a parental chromosome (Fig. 7). In our case, we could not tolerate the generation of heteroduplex artifacts because they occurred at levels greater than or equal to the fraction of endogenous molecules targeted for assay (those paternal IGF/H19 alleles and maternal IGF2R alleles that become methyalted at a particular site, see Figs 2 and 4).
The main limitation of methyl-sensitive restriction endonuclease assays is the ability to analyze only one or a few CpG dinucleotides in a region. On the other hand, this method is likely to provide a better quantitative estimate of the ratio of methylation between the two parental alleles because a much larger number of chromosomes are analyzed. In a sample of 100 ng DNA, the methylation status of the selected CpG site is assayed on over 15 000 homologous pairs of a particular chromosome, compared with the larger number of CpG sites (but generally less than 50) that are analyzed on the five to 10 homologous chromosome pairs usually assessed by bisulfite-sequencing. The incorporation of the hot-stop PCR technique circumvents the difficulty of analysis of heteroduplex molecules because only molecules synthesized during the last amplification cycle are radio-labeled (25).
LOI, parent of origin dependent methylation and cancer
The presence of LOI at IGF2 in both tumor cells and normal tissues raises the possibility that constitutional LOI represents a risk factor for cancer (7). Although use of the term LOI implies a qualitative difference between individuals in whom there is parental-origin-dependent transcription of only one allele and individuals in whom there is transcription of both alleles, such differences are generally quantitative (7). Individuals characterized as having LOI exhibit a relaxation of allele-specific methylation with partial methylation on the normally unmethylated allele. This methylation defect at IGF2/H19 is similar to what we found in the majority of individuals at the IGF2R DMR and one might argue that this circumstance may explain the normal biallelic expression of this gene. However, we also identified a number of individuals with nearly exclusive methylation of the maternal IGF2R allele. In these cases, one might hypothesize maternal allele-specific transcription, as described by Xu et al. (9). In this regard, the imprinting status of IGF2R is a potentially interesting variable in cancer susceptibility. A number of reports support its activity as a tumor suppressor gene in some types of tissues, including liver (26), breast (27) and kidney (28) and IGF2R is expressed from only the maternal allele in 50% of Wilms' tumors (28), consistent with the hypothesis that the transcriptional imprinting of this gene is associated with tumorigenesis.
Assays of allele-specific methylation exploit the fact that differential DNA methylation of maternal and paternal alleles is one of the hallmarks of genome imprinting (reviewed in 29). Many imprinted loci show allele-specific DNA methylation, and in the absence of the maintenance methyltransferase, or in the presence of methyltransferase inhibitors, the imprinting of several genes is disrupted (30,31). A major advantage of allele-specific methylation assays is that the methylation phenotype of an individual can be scored in a DNA sample. A similar surrogate measure of the transcriptional status of the androgen receptor gene has been used extensively as an indicator of X-chromosome inactivation in the cells of human females (32,33). A disadvantage of surrogate measures of imprinting is the perception that methylation differences, in and of themselves, do not reflect real imprinting unless they are also accompanied by parental allele-specific transcription (29). However, such differences must be at the root of the imprinting process, regardless of whether or not the process results in transcriptional differences between maternal and paternal alleles (34). The identification of inter-individual variability in the establishment or maintenance of epigenetic marks, and determining the causes of this variability, is an important step in determining the role of imprinting in evolution, development and cancer.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Subjects
DNA samples obtained from unfractionated nucleated peripheral blood cells (not lymphoblast cell lines) from the Salt Lake City collection of CEPH/Utah pedigrees between 1982 and 1984 (680 individuals from 48 three-generation families) were studied. For 25 of these families (288 individuals) a second set of DNA samples was obtained during 2001. All subjects gave informed consent under University of Utah I.R.B. approved protocol number 6090-96.
Methylation analysis using methyl-sensitive restriction endonucleases MluI, MaeII, NotI and XhoI
The upstream H19 sequence used in this study is available from GenBank (accession number AF125183), as is the sequence of IGF2R (GenBank accession number AF069333). In order to analyze allele-specific methylation we screened the DNA samples for two SNPs: a C/T polymorphism (Fig. 1A) recognized by CfoI at the IGF2/H19-DMR (close to the sixth binding site for CTCF) (16) and a C/T polymorphism recognized by MspI at the IGF2R-DMR (Fig. 1B). Parental origin of alleles of heterozygous individuals was determined by pedigree analysis.
A 100 ng sample of genomic DNA from informative individuals was digested overnight at 37°C with an excess of a methyl-sensitive restriction endonuclease: MluI or MaeII for the IGF2/H19-DMR and XhoI or NotI for the IGF2R-DMR respectively. In each experiment we used as control individuals who were homozygous for C alleles, homozygous for T alleles, as well as a mix of 50 ng genomic DNA from an individual homozygous for the C allele and 50 ng genomic DNA from an individual homozygous for the T allele. Digested DNA was then ethanol precipitated and amplified in a hot-stop PCR assay using the following primers: F3, 5'-GAGATGGGAGGAGATACTAGG-3' and R3, 5'-GTCAGTTCAGTAAAAGGCTGG-3' for the IGF2/H19-DMR, and F4, 5'-GACAGACCAGGATTAGCGTG-3' and R4, 5'-GGACATGCACTTCTCCCGTG-3' for IGF2R-DMR. Briefly, hot-stop PCR involves the addition, after 30 cycles, of a primer that has been radiolabeled at the 5' position and allowing synthesis to take place for a final cycle (25). The labeled primer for the last cycle was one of the same primers used in the first 30 cycles. However, competition between labeled and unlabeled primer did not appear to be a significant problem because the concentration of the unlabeled primer was substantially reduced during the first 30 cycles and the labeled primer was added at the original concentration of the unlabelled primer (0.4 µM).
PCR products were digested with the enzyme used for identification of the parental origin of the alleles, the fragments were separated on 5% polyacrylamide gels and the intensity of the bands (alleles) was quantified using a Fuji BAS 2000 phosphorimager (33). In the case of controls, no uncleaved C alleles were detected in any C/C homozygous individuals (indicating that the CfoI and MspI enzymes cleaved the PCR products with >99% efficiency) and the mixture of C homozygote/T homozygote resulted in a 1 : 1 ratio of C/T intensity (data not shown).
All the assays were performed three times, and the average values were used. At the IGF2R-DMR, all samples were assayed using both XhoI and NotI, while at the IGF2/H19-DMR, MaeII was used only for those samples with abnormal patterns of methylation after MluI digestion and for a limited number of samples with a normal pattern of methylation.
Quantification of the methylation level using real-time PCR
Samples of 100 ng of genomic DNA were digested under the same conditions as mentioned above using MluI, XhoI or NotI restriction endonucleases. As controls, samples of 100 ng genomic DNA were digested using DraI which cleaves outside the region of interest. For both the IGF2/H19-DMR and IGF2R DMR we selected representative samples with normal and abnormal methylation pattern in the previous method. After ethanol precipitation, PCR reactions were set up using 10 ng of DNA as template in the LightCycler DNA Master SYBR Green I kit (Roche), according to the manufacturer's instructions. The efficiency of the methyl-sensitive restriction enzymes was quantified as the ratio of fluorescence between the PCR products from DNA samples digested with MluI, XhoI or NotI enzymes and PCR products from the same samples digested with DraI restriction enzyme. The results were similar using the two methods and suggest that abnormal ratios are the result of increased methylation of the normally unmethylated alleles (data not shown).
Bisulfite-sequencing analysis at IGF2/H19 DMR
A 2.5 µg aliquot of genomic DNA was denatured at 37°C for 20 min using freshly prepared 3 M sodium hydroxide (final concentration of 0.3 M). Denatured DNA was subsequently incubated at 55°C with monosulfite solution (2 M sodium metabisulfite and 1 mM hydroquinone) in the dark, overnight (19). After purification with the QIAquick PCR purification kit (Qiagen), the converted DNA was used as a template for a nested PCR. The first set of primers was F1, 5'-GTGTTTGATTTA-TTTTAGGGTGTATTG-3' and R1, 5'-CCAACAACACAAAAATCCTAAACCA-AAA-3' and the second was F2, 5'-GGAATAATGAGGTGTTTTAGTTTTA-3' and R2, 5'-CTAACCACTTAAAACTAAAAAAATC-3'. The PCR products were separated in 1.5% agarose gels, purified using the QIAEX II gel extraction kit (Qiagen) and subsequently subcloned into a TA Cloning vector (Invitrogen) according to the manufacturer's instructions. The DNA from 1520 individual clones was sequenced for each individual.
Haplotype analysis
In order to analyze the concordance between a shared haplotype and familial aggregation of individuals with abnormal methylation ratios, we determined the genotypes of these individuals for several SNPs adjacent to the IGF/H19 and IGF2R loci. For chromosome 11 we analyzed seven polymorphisms, in the following physical order: D11S2071 (0.29 Mb from the telomere), the HRAS1 minisatellite, the A/C polymorphism at exon 5 of the H19 gene (recognized by RsaI), the A/G polymorphism at exon 5 of H19 gene (recognized by AluI), the C/T polymorphism at IGF2/H19-DMR (recognized by CfoI), the A/G polymorphism at exon 9 of IGF2 gene (recognized by ApaI enzyme), and the VNTR minisatellite at the 5' flanking region of insulin gene (2.27 Mb from the telomere). For chromosome 6 we analyzed the C/T polymorphism in intron 2 of IGF2R (recognized by MspI) and the C/T polymorphism at intron 35 of IGF2R gene (recognized by MspI). The genotypes of some individuals at the insulin VNTR and HRAS1 minisatellite loci were obtained from the CEPH database. Haplotypes were established by pedigree analysis.
| ACKNOWLEDGEMENTS |
|---|
We are grateful to Andrew Feinberg and Hengmi Cui for discussion and communication of results prior to publication, to Laurie Mecham, Tina Varvil and Joe Vaughan for administrative and technical assistance. We would like to extend our sincere thanks to all family members who participated in the Utah Genetic Reference Project. Thanks also to Andreas P. Peiffer, MD, PhD, UGRP Medical Director, and Melissa M. Dixon, UGRP Study Coordinator. This research was supported by a grant from the National Institutes of Health (NIH R21ES/CA11607) (to C.S.) and by a Public Health Services research grant to the Huntsman General Clinical Research Center at the University of Utah, grant number M01-RR00064 from the National Center for Research Resources. It was also supported by generous gifts from the W.M. Keck Foundation and from the George S. and Delores Doré Eccles Foundation.
| FOOTNOTES |
|---|
* To whom correspondence should be addressed. Tel: +1 2157077373; Fax: +1 2157071454; Email: sapienza{at}temple.edu
| REFERENCES |
|---|
|
|
|---|
-
Bird, A. (2002) DNA methylation patterns and epigenetic memory. Genes Dev., 16, 621.
[Free Full Text] - Li, E. (2002) Chromatin modification and epigenetic reprogramming in mammalian development. Nat. Rev. Genet., 3, 662673.[CrossRef][Web of Science][Medline]
-
Jenuwein, T. and Allis, C.D. (2001) Translating the histone code. Science, 293, 10741780.
[Abstract/Free Full Text] - Rainier, S., Johnson, L.A., Dobry, C.J., Ping, A.P., Grundy, P.E. and Feinberg, A.P. (1993) Relaxation of imprinted genes in human cancer. Nature, 362, 747749.[CrossRef][Medline]
- Ogawa, O., Eccles, M.R., Szeto, J., McNoe, L.A., Yun, K., Maw, M.A., Smith, P.J. and Reeve, A.E. (1993) Relaxation of insulin-like growth factor II gene imprinting implicated in Wilms' tumor. Nature, 362, 749751.[CrossRef][Medline]
-
Falls, J.G., Pulford, D.J., Wylie, A.A. and Jirtle, R.L. (1999) Genomic imprinting: implications for human disease. Am. J. Pathol., 154, 635647.
[Abstract/Free Full Text] - Cui, H., Horon, I.L., Ohlsson, R., Hamilton, S.R. and Feinberg, A.P. (1998) Loss of imprinting in normal tissue of colorectal cancer patients with microsatellite instability. Nat. Med., 4, 12761280.[CrossRef][Web of Science][Medline]
- Kalscheuer, V.M., Mariman, E.C., Schepens, M.T., Rehder, H. and Ropers, H.H. (1993) The insulin-like growth factor type-2 receptor gene is imprinted in the mouse but not in humans. Nat. Genet., 5, 7478.[CrossRef][Web of Science][Medline]
- Xu, Y., Goodyer, C.G., Deal, C. and Polychronakos, C. (1993) Functional polymorphism in the parental imprinting of the human IGF2R gene. Biochem. Biophys. Res. Commun., 197, 747754.[CrossRef][Web of Science][Medline]
- Nishiwaki, K., Niikawa, N. and Ishikawa, M. (1997) Polymorphic and tissue-specific imprinting of the human Wilms tumor gene, WT1. Jpn. J. Hum. Genet., 42, 205211.[Medline]
- Sakatani, T., Wei, M., Katoh, M., Okita, C., Wada, D., Mitsuya, K., Meguro, M., Ikeguchi, M., Ito, H., Tycko, B. and Oshimura, M. (2001) Epigenetic heterogeneity at imprinted loci in normal populations. Biochem. Biophys. Res. Commun., 283, 11241130.[CrossRef][Web of Science][Medline]
-
Nakagawa, H., Chadwick, R.B., Peltomaki, P., Plass, C., Nakamura, Y. and de La Chapelle, A. (2001) Loss of imprinting of the insulin-like growth factor II gene occurs by biallelic methylation in a core region of H19-associated CTCF-binding sites in colorectal cancer. Proc. Natl Acad. Sci. USA, 98, 591596.
[Abstract/Free Full Text] - Vu, T.H., Li, T., Nguyen, D., Nguyen, B.T., Yao, X.M., Hu, J.F. and Hoffman, A.R. (2000) Symmetric and asymmetric DNA methylation in the human IGF2-H19 imprinted region. Genomics, 64, 132143.[CrossRef][Web of Science][Medline]
-
Srivastava, M., Hsieh, S., Grinberg, A., Williams-Simons, L., Huang, S.P. and Pfeifer, K. (2000) H19 and Igf2 monoallelic expression is regulated in two distinct ways by a shared cis acting regulatory region upstream H19. Genes Dev., 14, 11861195.
[Abstract/Free Full Text] - Hark, A.T., Schroenherr, C.J., Katz, D.J., Ingram, R.S., Levorse, J.M. and Tilghman, S.M. (2000) CTCF mediates methylation-sensitive enhancer-blocking activity at the H19/Igf2 locus. Nature, 405, 486489.[CrossRef][Medline]
-
Takai, D., Gonzales, F.A., Tsai, Y.C., Thayer, M.J. and Jones, P.A. (2001) Large scale mapping of methylcytosines in CTCF-binding sites in the human H19 promoter and aberrant hypomethylation in human bladder cancer. Hum. Mol. Genet., 10, 26192626.
[Abstract/Free Full Text] -
Smrzka, O.W., Fae, I., Stoger, R., Kurzbauer, R., Fisher, G.F., Henn, T., Weith, A. and Barlow, D.P. (1995) Conservation of a maternal-specific methylation signal at the human IGF2R locus. Hum. Mol. Genet., 4, 19451952.
[Abstract/Free Full Text] -
Frommer, M., McDonald, L.E., Millar, D.S., Collis, C.M., Watt, F., Grigg, G.W., Molloy, P.L. and Paul, C.L. (1992) A genomic sequencing protocol that yields a positive display of 5-methylcytosine residues in individual DNA strands. Proc. Natl Acad. Sci. USA, 89, 18271831.
[Abstract/Free Full Text] -
Clark, S.J., Harrison, J., Paul, C.L. and Frommer, M. (1994) High sensitivity mapping of methylated cytosines. Nucl. Acids Res., 22, 29902997.
[Abstract/Free Full Text] -
Thompson, J.R., Marcelino, L.A. and Polz, M.F. (2002) Heteroduplexes in mixed-template amplifications: formation, consequence and elimination by reconditioning PCR. Nucl. Acids Res., 30, 20832088.
[Abstract/Free Full Text] - Lowell, J.L. and Klein, D.A. (2000) Heteroduplex resolution using T7 endonuclease I in microbial community analyses. Biotechniques, 28, 676678, 680, 681.[Web of Science][Medline]
-
Grunau, C., Clark, S.J. and Rosenthal, A. (2001) Bisulfite genomic sequencing: systematic investigation of critical experimental parameters. Nucl. Acids Res., 29, e65.
[Abstract/Free Full Text] - Warnecke, P.M., Stirzaker, C., Song, J., Grunau, C., Melki, J.R. and Clark, S.J. (2002) Identification and resolution of artifacts in bisulfite sequencing. Methods, 27, 101107.[CrossRef][Web of Science][Medline]
- Grilley, M., Holmes, J., Yashar, B. and Modrich, P. (1990) Mechanisms of DNA-mismatch correction. Mutat. Res., 236, 253267.[Web of Science][Medline]
- Uejima, H., Lee, M.P., Cui, H. and Feinberg, A.P. (2000) Hot-stop PCR: a simple and general assay for linear quantification of allele ratios. Nat. Genet., 25, 375376.[CrossRef][Web of Science][Medline]
- De Souza, A.T., Hankins, G.R., Washington, M.K., Orton, T.C. and Jirtle, R.L. (1995) M6P/IGF2R gene is mutated in human hepatocellular carcinomas with loss of heterozigosity. Nat. Genet., 11, 447449.[CrossRef][Web of Science][Medline]
- Hankins, G.R., De Souza, A.T., Bentley, R.C., Patel, M.R., Marks, J.R., Iglehart, J.D. and Jirtle, R.L. (1996) M6P/IGF2R receptor: a candidate breast tumor suppressor gene. Oncogene, 12, 20032009.[Web of Science][Medline]
- Xu, Y.Q., Grundy, P. and Polychronakos, C. (1997) Aberrant imprinting of the insulin-like growth factor II receptor gene in Wilms' tumor. Oncogene, 14, 10411046.[CrossRef][Web of Science][Medline]
- Tilgman, S.M. (1999) The sins of the fathers and mothers: genomic imprinting in mammalian development. Cell, 96, 185193.[CrossRef][Web of Science][Medline]
- Li, E., Beard, C. and Jaenisch, R. (1993) Role for DNA methylation in genomic imprinting. Nature, 366, 362365.[CrossRef][Medline]
-
Dao, D., Walsh, C.P., Yuan, L., Gorelov, D., Feng, L., Hensle, T., Nisen, P., Yamashiro, D.J., Bestor, T.H. and Tycko, B. (1999) Multipoint analysis of human chromosome 11p15/mouse distal chromosome 7: Inclusion of H19/IGF2 in the minimal WT2 region, gene specificity of H19 silencing in Wilms' tumorigenesis and methylation hyper-dependence of H19 imprinting. Hum. Mol. Genet., 8, 13371352.
[Abstract/Free Full Text] - Allen, R.C., Zoghbi, H.Y., Moseley, A.B., Rosenblatt, H.M. and Belmont, J.W. (1992) Methylation of HpaII and HhaI sites near the polymorphic CAG repeat in the human androgen-receptor gene correlates with X chromosome inactivation. Am. J. Hum. Genet., 51, 12291239.[Web of Science][Medline]
- Naumova, A.K., Bird, L.M., Slamka, C., Fonseca, M., Verner, A.E., Wang, M., Leppert, M., Morgan, K. and Sapienza, C. (1995) Sex-specific transmission-ratio distorsion of X-chromosome alleles in the offspring of females with nonrandom X-inactivation, Dev. Genet., 17, 198205.[CrossRef][Web of Science][Medline]
-
Pardo-Manuel de Villena, F., de la Casa Esperon, E. and Sapienza, C. (2000) Natural selection and the function of genome imprinting: beyond the silenced minority. Trends Genet., 16, 573579.[CrossRef][Web of Science][Medline]
This article has been cited by other articles:
![]() |
R. Yang, X. Wang, Z. Wu, D. R. Prows, and M. Lin Bayesian model selection for characterizing genomic imprinting effects and patterns Bioinformatics, January 15, 2010; 26(2): 235 - 241. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Suzuki, Y. Li, X. Dong, M. M. Hassan, J. L. Abbruzzese, and D. Li Effect of Insulin-Like Growth Factor Gene Polymorphisms Alone or In Interaction with Diabetes on the Risk of Pancreatic Cancer Cancer Epidemiol. Biomarkers Prev., December 1, 2008; 17(12): 3467 - 3473. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Ito, T. Koessler, A. E.K. Ibrahim, S. Rai, S. L. Vowler, S. Abu-Amero, A.-L. Silva, A.-T. Maia, J. E. Huddleston, S. Uribe-Lewis, et al. Somatically acquired hypomethylation of IGF2 in breast and colorectal cancer Hum. Mol. Genet., September 1, 2008; 17(17): 2633 - 2643. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. T. Bjornsson, M. I. Sigurdsson, M. D. Fallin, R. A. Irizarry, T. Aspelund, H. Cui, W. Yu, M. A. Rongione, T. J. Ekstrom, T. B. Harris, et al. Intra-individual Change Over Time in DNA Methylation With Familial Clustering JAMA, June 25, 2008; 299(24): 2877 - 2883. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Bruce, K. Hannula-Jouppi, C. M. Lindgren, M. Lipsanen-Nyman, and J. Kere Restriction Site-Specific Methylation Studies of Imprinted Genes with Quantitative Real-Time PCR Clin. Chem., March 1, 2008; 54(3): 491 - 499. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. H. Taylor, R. S. Kramer, J. W. Davis, J. Guo, D. J. Duff, D. Xu, C. W. Caldwell, and H. Shi Ultradeep Bisulfite Sequencing Analysis of DNA Methylation Patterns in Multiple Gene Promoters by 454 Sequencing Cancer Res., September 15, 2007; 67(18): 8511 - 8518. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. T. Heijmans, D. Kremer, E. W. Tobi, D. I. Boomsma, and P. E. Slagboom Heritable rather than age-related environmental and stochastic factors dominate variation in DNA methylation of the human IGF2/H19 locus Hum. Mol. Genet., March 1, 2007; 16(5): 547 - 554. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Sparago, S. Russo, F. Cerrato, S. Ferraiuolo, P. Castorina, A. Selicorni, C. Schwienbacher, M. Negrini, G. B. Ferrero, M. C. Silengo, et al. Mechanisms causing imprinting defects in familial Beckwith-Wiedemann syndrome with Wilms' tumour Hum. Mol. Genet., February 1, 2007; 16(3): 254 - 264. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Tost, H. Jammes, J.-M. Dupont, C. Buffat, B. Robert, T.-M. Mignot, F. Mondon, B. Carbonne, U. Simeoni, G. Grange, et al. Non-random, individual-specific methylation profiles are present at the sixth CTCF binding site in the human H19/IGF2 imprinting control region Nucleic Acids Res., November 14, 2006; 34(19): 5438 - 5448. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. W. Santure and H. G. Spencer Influence of Mom and Dad: Quantitative Genetic Models for Maternal Effects and Genomic Imprinting Genetics, August 1, 2006; 173(4): 2297 - 2316. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. D. Anway and M. K. Skinner Epigenetic Transgenerational Actions of Endocrine Disruptors Endocrinology, June 1, 2006; 147(6): s43 - s49. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. A. Waterland, J.-R. Lin, C. A. Smith, and R. L. Jirtle Post-weaning diet affects genomic imprinting at the insulin-like growth factor 2 (Igf2) locus Hum. Mol. Genet., March 1, 2006; 15(5): 705 - 716. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Sandovici, S. Kassovska-Bratinova, J. C. Loredo-Osti, M. Leppert, A. Suarez, R. Stewart, F. D. Bautista, M. Schiraldi, and C. Sapienza Interindividual variability and parent of origin DNA methylation differences at specific human Alu elements Hum. Mol. Genet., August 1, 2005; 14(15): 2135 - 2143. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. A. McCann, Y. Q. Xu, R. Frechette, L. Guazzarotti, and C. Polychronakos The Insulin-Like Growth Factor-II Receptor Gene Is Associated with Type 1 Diabetes: Evidence of a Maternal Effect J. Clin. Endocrinol. Metab., November 1, 2004; 89(11): 5700 - 5706. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Murrell, S. Heeson, W. N. Cooper, E. Douglas, S. Apostolidou, G. E. Moore, E. R. Maher, and W. Reik An association between variants in the IGF2 gene and Beckwith-Wiedemann syndrome: interaction between genotype and epigenotype Hum. Mol. Genet., January 15, 2004; 13(2): 247 - 255. [Abstract] [Full Text] [PDF] |
||||
![]() |
H.T. BJORNSSON, H. CUI, D. GIUS, M.D. FALLIN, and A.P. FEINBERG The New Field Of Epigenomics: Implications for Cancer and Other Common Disease Research Cold Spring Harb Symp Quant Biol, January 1, 2004; 69(0): 447 - 456. [Abstract] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

















