Molecular and Cellular Biology, August 2003, p. 5475-5488, Vol. 23, No. 16
0270-7306/03/$08.00+0 DOI: 10.1128/MCB.23.16.5475-5488.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Philippe Arnaud,1 Emma Gordon,1 Wendy Dean,1 Elizabeth A. Coar,1 Christine M. Williamson,2 Robert Feil,3 Jo Peters,2 and Gavin Kelsey1*
Developmental Genetics Programme, The Babraham Institute, Cambridge CB2 4AT,1 MRC Mammalian Genetics Unit, Harwell, Didcot, Oxfordshire OX11 0RD, United Kingdom,2 Institute of Molecular Genetics, CNRS, UMR-5535, 34293 Montpellier, France3
Received 20 December 2002/ Returned for modification 7 April 2003/ Accepted 20 May 2003
|
|
|---|
|
|
|---|
One of the first imprinted effects described was identified from uniparental inheritance of the distal region of chromosome (Chr) 2 in the mouse (7). Maternal and paternal duplications of this region were found to cause striking and superficially opposite neonatal phenotypes, with behavioral and morphological effects. Through the use of a number of reciprocal translocations, the region responsible for the imprinted phenotypes was narrowed down to an
7-Mb interval (36, 53), and by methylation-sensitive representational difference analysis, we subsequently identified a complex imprinted cluster at the Gnas locus (24, 37). Gnas encodes the stimulatory G-protein subunit Gs
. In addition to the coding transcript for Gs
, the locus was found to comprise two imprinted transcripts: Nesp expressed from the maternal allele (which codes for the chromogranin-like neuroendocrine secretory protein NESP55 [21]) and Gnasxl expressed from the paternal allele (which codes for XL
S, a variant Gs
that has a large noncanonical amino-terminal domain [23]). These two transcripts arise from alternative upstream promoters, and both transcripts are spliced to exon 2 of Gnas and contain downstream exons in common with Gnas. Additional complexity of the locus has emerged from identification of a noncoding transcript, Nespas, which runs antisense to Nesp (27, 54, 55), and an alternative noncoding first exon for Gnas with paternal-specific expression (29). The human GNAS locus has a very similar organization (17-19, 28).
Imprinting of GNAS had been implicated from the different clinical manifestations of inactivating mutations of Gs
, which cause the autosomal dominant disorder Albright's hereditary osteodystrophy (8, 52). Maternally inherited mutations in GNAS are associated with multihormone resistance, a condition referred to as pseudohypoparathyroidism type 1a (PHP1a), because patients present with renal resistance to parathyroid hormone. Tissue-specific imprinting in both humans and mice has subsequently been described, with exclusive or prominent expression of the maternal allele in sites such as proximal renal tubules, brown and white adipose tissue, and the pituitary and thyroid glands (13, 16, 30, 61). Imprinting in these target tissues accounts for some of the endocrine anomalies (52). The Nesp and Gnasxl promoters, in contrast, display monoallelic expression at all sites in which they are expressed (18, 19, 27, 37).
The cis-acting elements that control imprinting at the Gnas and GNAS clusters and the mechanisms by which monoallelic expression of the various promoters is executed are not known. Imprinted control regions (ICRs) defined by deletion analysis at other loci coincide with differentially methylated regions (DMRs), where the methylation state of the two parental alleles differs markedly (11, 48, 56, 60) and where methylation of one allele is laid down in the respective germ line (41, 45, 50, 59). Female germ line methylation at imprinted loci depends upon the DNA methyltransferases Dnmt3a and Dnmt3b and the related protein Dnmt3l (6, 15), but the sequence features that specify DMRs for methylation in either germ line are not known. Direct repeats often found within or adjacent to DMRs have been implicated on the basis that such repeats are not found in CpG islands of nonimprinted genes (32) or nonimprinted homologues in other species (33, 34). Such direct repeats have been noted at the human GNAS locus in the NESP55 and XL
S exons (18, 19). Ultimately, differential methylation, in concert with specific chromatin organization at ICRs, is translated into monoallelic expression of linked promoters in somatic tissues by a variety of mechanisms (39, 43).
Three DMRs have been identified at Gnas and GNAS. The Nesp DMR has paternal methylation, while the Gnasxl DMR and a DMR covering Gnas exon 1A (also referred to as exon A/B) have maternal methylation (18, 19, 24, 28, 29, 37). The exon 1A DMR has been shown to be a gametic methylation mark in the mouse (29); the equivalent human region may also be a primary DMR (22). The control of tissue-specific imprinting of the Gs
-coding transcript may reside in the exon 1A region, as patients with hormone resistance in the absence of the other features of Albright's hereditary osteodystrophy (a condition known as PHP1b) almost invariably display altered methylation at exon 1A (3, 28). Whether this DMR controls imprinting of the entire complex locus is not clear: PHP1b patients may or may not also show altered methylation at the NESP55 and XL
S DMRs. Here, we present a characterization of the epigenetic properties of the mouse Gnas locus as a means of pinpointing potential ICRs and predicting their possible modes of actions. We have mapped the extent of differential methylation at the Nesp and Nespas/Gnasxl DMRs, examined gross chromatin organization as revealed by DNase I sensitivity, and identified a second germ line DMR.
|
|
|---|
Collection of gametes and early embryos. Oocytes were obtained from juvenile F1 (C57BL/6J x CBA/Ca) mice; morulae and blastocysts were obtained from an F1 x F1 cross, except where indicated. Oocytes were collected from superovulated immature females, as described by Hogan et al. (20). Mature spermatozoa were isolated from epididymis of adult CBA/Ca mice.
ES cells.
Embryonic stem (ES) cells used in this study have been described previously (9). Hybrid ES cell line SF1-1 was obtained from F1 x Mus spretus hybrid blastocysts created by in vitro fertilization. Monoparental ES cell lines used were AG-A (androgenetic) and PR-8 (parthenogenetic). Cells were cultured on gelatin-coated flasks (0.1% gelatin) with feeder cells (
-irradiated primary embryonic fibroblasts) at 37°C under 5% CO2, in Iscove's modified Dulbecco's modified Eagle medium supplemented with 15% fetal calf serum, recombinant mouse leukemia inhibitory factor (20 ng/ml), penicillin (50 U/ml), streptomycin (50 µg/ml), 0.1 mM ß-mercaptoethanol, 1x modified Eagle medium, nonessential amino acids, and 2 mM glutamine. Prior to harvesting, ES cells were passaged onto gelatin-coated flasks in the absence of feeders to reduce their contribution to the final cell pellet.
Southern analysis of methylation.
Embryos (12.5 days postcoitum [dpc]) with uniparental partial disomy for distal Chr 2 were generated by standard methods of intercrossing reciprocal translocation heterozygotes and have been described before (24, 37). DNAs (10 µg per reaction) were digested with the enzymes indicated, together with Bsh1236I, Hin6I, HpaII, or MspI, resolved by electrophoresis on 1% agarose-TAE gels (TAE is 40 mM Tris-acetate, 10 mM EDTA), and transferred by capillary blotting onto charged nylon membranes. Probes were restriction fragments subcloned into pBluescript II KS(+) (Stratagene) from genomic phage and cosmids for Nesp and Gnasxl (24). Hybridizations were performed with gel-purified probes labeled with [
-32P]dCTP (ICN) by random priming.
Bisulfite sequence analysis.
Oocytes (200 to 600), morulae (5 to 20), or blastocysts (5 to 8) were resuspended in 32.5 µl of a solution containing 10 µg of glycogen, 1 mM sodium dodecyl sulfate, and proteinase K (280 µg/ml) and were incubated for 90 min at 37°C and then for 15 min at 95°C in a thermocycler. The resulting DNA lysate was denatured by addition of 1.1 µl of 10 N NaOH and incubation at 50°C for 15 min. For bisulfite treatment, 200 µl of
4 M sodium bisulfite, pH 5.0 (final concentration,
3.5 M; Sigma); 1.5 µl of 75 mM hydroquinone (final concentration, 0.5 mM; Sigma); and 5 µg of glycogen were added, and DNA incubated at 55°C for 4 h. Desalting was carried out using the QIAquick PCR purification kit (Qiagen), and eluted DNA (in 50 µl Tris-HCl, pH 7.5) was desulfonated by treatment with 1.6 µl of 10 N NaOH. DNA was ethanol precipitated and resuspended in H2O (5 µl per 100 cell equivalents). A nested primer strategy was used to amplify bisulfite-treated oocyte and early embryo DNA. PCR, cloning, and sequencing were performed as previously described (44). Primer sequences are available on request. Prior to cloning, PCR products were tested for full conversion and methylation status by pilot digestion with appropriate restriction enzymes.
Isolation of nuclei and DNase I sensitivity analysis.
Tissues for isolation of nuclei were obtained from (C57BL/6J x M. spretus) mice and the backcross offspring from (C57BL/6J x M. spretus) females to C57BL/6J males (the latter were genotyped by PCR for the presence of M. spretus alleles at D2Mit22 and D2Mit74). Nuclei were isolated from frozen tissues (brain, liver or kidney) after disruption under liquid nitrogen and homogenization, as described elsewhere (25). For ES cells, 5 x 107 to 5 x 108 cells were harvested for preparation of nuclei and treated as previously described (25). DNase I digestion of nuclei was performed immediately after isolation. Nuclei (aliquots of
107 suspended in 200 µl of 60 mM KCl, 15 mM NaCl, 5 mM MgCl2, 0.1 mM EGTA, 15 mM Tris-HCl [pH 7.5], 0.5 mM dithiothreitol, 0.3 M sucrose, 5% [vol/vol] glycerol) were treated with DNase I (Roche grade I) at 0 to 750 U/ml at 25°C for 10 min. Digestion was stopped by addition of a solution containing 200 µl of 20 mM EDTA, 1% (wt/vol) sodium dodecyl sulfate, and proteinase K (200 µg/ml) and treatment at 50°C for 16 h. DNA was purified by extraction with phenol-chloroform-isoamylalcohol (25:24:1) and chloroform-isoamylalcohol (24:1), precipitated with ethanol, and resuspended in 10 mM Tris-HCl-1 mM EDTA (pH 8.0). For detection of DNase I-hypersensitive sites (HSSs), 20 µg of each treated DNA was digested with appropriate restriction enzymes, electrophoresed, blotted and hybridized as above. Probes were PCR fragments labeled directly with [
-32P]dCTP (ICN) by random priming (primer details available on request).
|
|
|---|
![]() View larger version (50K): [in a new window] |
FIG. 1. The Nesp-Gnasxl domain of the Gnas imprinted cluster. (A) Schematic overview of the mouse Gnas locus. Exons of the Nesp, Gnasxl and Gnas transcripts are shown above the line (for simplicity, not all Gnas exons are shown): coding regions are filled, noncoding regions open. Exons for the Nespas antisense transcript are shown below the line. Promoters are indicated by horizontal arrows, with those maternally expressed (mat) shown above the line and those paternally expressed (pat) shown below the line. The Nesp and Gnasxl exons are 48.8 and 34.0 kb upstream of Gnas exon 2, respectively, onto which they are both spliced. (B) Sequence properties of the Nesp-Gnasxl domain. A graphical output of sequence AJ251761 analyzed by CPGPLOT is shown. CpG islands (CGIs) are identified as boxes.
|
S domain. The human XL
S exon contains similar CG-rich repeats (Table 1) (18). The Nesp or Gnasxl repeats do not comprise reiterations of good matches to transcription factor binding site motifs, in particular, multiple sites for the insulator and boundary factors CTCF and YY1 are absent. Two less reiterated tandem repeats (a 16-mer at 19363 to 19399 and a 24-mer at 21631 to 21687), which do not contribute CpG dinucleotides, are located downstream of the Gnasxl exon and CpG islands.
![]() View larger version (29K): [in a new window] |
FIG. 2. Extent of differential methylation across the Nesp-Gnasxl domain. (A) The restriction fragments analyzed on Southern blots are indicated by the horizontal lines, and the respective restriction enzyme sites are abbreviated as follows: B, BamHI; E, Eco32I; H, HindIII; and R, EcoRI. The probes used are shown as striped boxes. Below, the methylation status of HpaII sites in 12.5-dpc embryo DNAs is summarized. Each vertical line represents a single HpaII site (in a few cases, two inseparable sites). Methylation is given separately on the maternal (mat) and paternal (pat) alleles, where filled circles represent fully methylated, open circles represent unmethylated, and paired circles represent partially methylated. Those sites grouped in square brackets are all methylated or unmethylated, as indicated. For sites without symbols, methylation state could not be determined. Methylation status was also assayed for many Hin6I and Bsh1236I sites, which did not differ appreciably from that for HpaII, but is not shown for the sake of clarity. The methylation summary includes results from difference product clones (24). The AJ251761 sequence analyzed in Fig. 2B extends between the HindIII and EcoRI sites marked (*); the outer BamHI (B) sites are situated at 82392 and 36822 in the sequence of mouse BAC RP23-439H2 (AL593857). (B) Representative Southern blots showing digests of 12.5-dpc MatDp(dist2) (Mat) and PatDp(dist2) (Pat) embryo DNAs hybridized with the probes indicated. DNAs are digested with the restriction enzymes indicated above each blot in combination with no other enzyme (lane 1), HpaII (lane 2), MspI (lane 3), Hin6I (lane 4), or Bsh1236I (lane 5).
|
|
View this table: [in a new window] |
TABLE 1. Tandem repeats in Nesp and Gnasxl DMRs
|
4.4 kb, covering the two Nesp exons and the CpG islands. Downstream of the Nesp DMR, CpG density declines (Fig. 1B) and an
4.0-kb region contains sites with full or partial methylation on both maternal and paternal alleles. Further downstream, there is a transition region in which the paternal allele shows less complete methylation before the very low level of methylation that characterizes the DMR (probe M1 on a BamHI digest; Fig. 2B). Methylation on the maternal allele is extensive (probe X6 on a EcoRI digest [Fig. 2B]), spanning
6.6 kb, including both Nespas and Gnasxl promoters and the CpG island region. Downstream of the Gnasxl exon, CpG density falls away and the methylation pattern is complex, with partially methylated as well as unmethylated sites (Fig. 2B; 12.3-kb BamHI fragment analyzed with probe X2); however, there is no marked difference between MatDp(dist2) and PatDP(dist2) DNAs. Moreover, we found no differential methylation further downstream until the DMR at Gnas exon 1A (29) which, as expected, showed maternal methylation in this material (data not shown). In conclusion, we have mapped the extents of the DMRs, as present in midgestation embryos, find them each to extend several kilobases and to coincide with the regions of greatest CpG density. Extensive germ line methylation mark at the Nespas-Gnasxl DMR. A hallmark of an ICR is that distinct methylation patterns are established in male and female gametes, and differential methylation is maintained in the zygote and during embryonic development. At Gnas, exon 1A has been shown by bisulfite genomic sequencing to be contained in a region methylated in oocyte DNA and unmethylated in sperm DNA (29). Whether this gametic DMR controls imprinting of the entire locus, including the Nesp-Gnasxl domain, is unclear. The report from Liu et al. (29) found no evidence of gametic methylation marks in the Nesp and Gnasxl DMRs, but only very limited regions were examined. We have undertaken a more extensive characterization of these regions.
Methylation in gametes and preimplantation embryos was determined by sequencing PCR products obtained from bisulfite-treated DNAs (Fig. 3). At the Nespas-Gnasxl DMR, analysis focused on two regions represented in four PCR products. PCR products b through d examined regions corresponding to prominent HSSs present specifically on the paternal allele in ES cells (see below), product d additionally examined a highly conserved region at the putative Nespas promoter (17, 27, 54), PCR product e covers the putative promoter for Gnasxl (Williamson et al., unpublished data). Bisulfite-modified sequences revealed a high level of methylation in oocyte DNA and very low levels in sperm DNA for each of these regions. In addition, sequences from preimplantation embryo DNAs (morulae or blastocysts) showed equal proportions of methylated and unmethylated sequences, suggestive of the maintenance of the oocyte and sperm derived methylation patterns after fertilization. To assign parental allele origin of the methylation present in preimplantation embryo DNAs, a C57BL/6J (B6) versus CBA/Ca (CBA) single-nucleotide polymorphism was identified in the Nespas promoter/conserved region (PCR product d). Bisulfite analysis made on DNA from morulae resulting from a B6 x CBA cross confirmed that methylated molecules were derived exclusively from the maternal allele (Fig. 3; PCR product d).
![]() View larger version (48K): [in a new window] |
FIG. 3. Germ line and early embryo methylation of the Nespas-Gnasxl DMR. The features of the regions analyzed are indicated at the top. The Gnasxl exon is depicted as a bar with coding portion in black, untranslated in grey and direct repeats indicated by arrowheads; the Nespas exon 1 is an open box. The location of a MER DNA transposon is indicated by the stippled box. The approximate positions of HSSs present specifically on the paternal allele in ES cells are shown by vertical arrows (Fig. 6). Nucleotide positions are given according to sequence AJ251761. The extents of the PCR products sequenced after bisulfite modification of DNAs are represented by open bars labeled a to f. Methylation status is given below. Each line of circles represents an individual sequence molecule, with each circle corresponding to a separate CpG. Methylated CpGs are indicated by filled circles, nonmethylated CpGs by open circles. For PCR product d, the dot represents the position of a single-nucleotide polymorphism (5'-GGTCGG-3' to 5'-GGTCTG-3') found between C57BL/6J and CBA/Ca, which results in loss of the indicated CpG in the CBA/Ca sequence. Sperm DNA is from CBA/Ca and oocyte DNA from (C57BL/6J x CBA/Ca)F1s. Morula (M) and blastocyst (B) DNAs were obtained from an F1 x F1 cross, except for PCR product d, for which morulae were from a C57BL/6J x CBA/Ca cross, with maternal B6 allele sequences and paternal CBA allele sequences identified as mat and pat, respectively
|
For Nesp, because the promoter region had been analyzed previously (29), we obtained bisulfite sequences for two elements potentially able to attract de novo methylation: the CpG-containing direct repeats in the Nesp exon; and a B1 element (58) 2 kb downstream of Nesp. These regions were found to be unmethylated or partially methylated in sperm DNA and, where examined, were unmethylated in oocytes and morulae (data not shown). Therefore, we confirm that the Nesp DMR does not have the properties of a methylation imprint mark.
In conclusion, bisulfite sequence analysis revealed an extensive gametic methylation mark at the Nespas-Gnasxl DMR, covering >3.2 kb. The upstream extent of the methylation mark coincides with the boundary of the CpG-rich region, with the extent of the somatic DMR, and maps close to prominent ES cell-specific DNase I HSSs described below.
Investigating the chromatin organization of the Nesp-Gnasxl domain. As a further indication of elements likely to regulate imprinting of the locus, we examined chromatin organization, as revealed by hypersensitivity to DNase I in isolated nuclei. Nuclei were prepared from tissues from adult mice and from ES cell lines, as a representation of the inner cell mass of preimplantation embryos. Mice were (B6 x M. spretus)F1 hybrids or backcross offspring from F1 hybrid females to B6 males (which hereafter we refer to as M. spretus x B6 for simplicity), which provided restriction fragment length polymorphisms (RFLPs) to distinguish maternally derived and paternally derived alleles. The choice of tissues included those in which Nesp and Gnasxl are expressed in a significant proportion of cells, i.e., brain, and essentially nonexpressing tissues liver and kidney. Analysis of DNA methylation in these adult tissues indicated that the Nesp and Gnasxl DMRs exist essentially as in midgestation embryos (data not shown). The ES cells used were of three types and have been used previously for investigating chromatin at imprinted loci (9, 25). SF1-1 is a (B6 x CBA) x M. spretus hybrid cell line with paternal M. spretus alleles; PR-8 is derived from diploid parthenogenetic embryos which only contain oocyte-derived chromosomes; and AG-A is from androgenetic embryos which contain only sperm-derived chromosomes. The Nesp DMR was hypomethylated on both alleles in SF1-1, unmethylated in PR-8, and largely unmethylated in AG-A (Fig. 4B); the Nespas-Gnasxl DMR was unmethylated in AG-A, methylated in PR-8, and showed the appropriate differential methylation in SF1-1 (see Fig. 6). This pattern, in keeping with the finding of the bisulfite sequence analysis, is consistent with the Nespas-Gnasxl DMR being gametic in origin, while the Nesp DMR becomes established after implantation. Despite these methylation patterns, reverse transcription-PCR assays showed that Nesp was expressed specifically from the maternal allele in the hybrid ES cells SF1-1 and not detected in AG-A ES cells, and Gnasxl was expressed in AG-A ES cells and from the paternal allele in the hybrid cells (data not shown).
![]() View larger version (47K): [in a new window] |
FIG. 4. DNase I hypersensitivity analysis of the Nesp region. (A) Analysis in nuclei from adult mouse kidney and brain, from (B6 x M. spretus) hybrids, in which the paternal allele is of M. spretus origin, and backcross (labeled spretus x B6) in which the maternal allele is M. spretus. Nuclei in the following lanes were digested with DNase I at the indicated concentrations: lane 1, 0 U/ml; lane 2, 50 U/ml; lane 3, 200 U/ml; lane 4, 400 U/ml; and lane 5, 750 U/ml. Purified DNA was digested with DraI and electrophoresed, and Southern blots were hybridized with probes D1 or D2. Maternal (M) and paternal (P) alleles are distinguished as a DraI RFLP. Points on the left of each blot represent DNA markers formed by DraI (D)-cut DNA digested with EcoRI (R), Eco32I (E), or XbaI (X). The M. spretus-specific DraI fragment is marked (*D). Location of thesereference sites is given in C. DNase I cleavages are indicated by arrows to the right of the blots, weaker cleavages in grey, more prominent sites by the labeled black arrows. (B) Analysis in nuclei from ES cell lines (see text for description of cell lines). The key is the same as that for panel A, except that lanes marked H and M are untreated DNAs digested with DraI plus HpaII or MspI. (C) Interpretation of DNase I HSSs. The DraI fragment analyzed is shown, with nucleotide positions according to sequence AJ251761. Nesp exons are given as filled boxes, with coding portions in black and untranslated regions in grey; Nespas exons are given as open boxes. The locations of probes D1 and D2 are given as striped boxes. The approximate positions of HSSs are indicated by the vertical arrows and represent a summation of mapping with probes D1 and D2 in A and B and from additional blots (not shown).
|
![]() View larger version (40K): [in a new window] |
FIG. 6. DNase I analysis of the Nespas promoter region in ES cells. (A) Analysis of nuclei from SF1-1 (F1 x M. spretus) ES cells, using an EcoRI RFLP. Nuclei were digested with DNase I (as in Fig. 4), purified DNA was digested with EcoRI and electrophoresed, and Southern blots were hybridized with probes R1 and R2. Maternal (M) and paternal (P) alleles are indicated to the right of the blot. Points on the left of the blot represent DNA markers formed by EcoRI (R) cut DNA separately digested with HindIII (H), SacI (S) or TaqI (T). The M. spretus-specific EcoRI fragment is given as *R. Location of these reference sites is shown in D. DNase I cleavages are indicated by arrows to the right of each blot, weaker cleavages in grey, more prominent sites in black. Lanes marked H and M are untreated DNAs digested with EcoRI plus HpaII or MspI. (B) Similar analysis in nuclei from adult mouse brain (B6 x M. spretus hybrid). (C) Analysis using probe M1 in BamHI digested DNAs, for fine mapping of the ES cell-specific HSSs. Restriction site markers are BamHI (B), HindIII (H), NcoI (N), Eco32I (E) and TaqI (T). (D) Interpretation of DNase I HSSs. The start of the Gnasxl exon is represented as a filled box with 5' untranslated region in grey, Nespas exon 1 as the open box. Locations of probes M1, R1 and R2 are given as striped boxes. The DNase I sites mapped from these and additional blots (not shown) are indicated using the convention above. HSSs on the maternal allele are labeled M, on the paternal allele P.
|
Imprinted chromatin features at the Nespas-Gnasxl DMR: prominent chromatin features specific to ES cells. The Nespas-Gnasxl region was analyzed in adult tissues using a ScaI RFLP in ScaI-XbaI digests (Fig. 5A). By using probe SX2 upstream of the polymorphic ScaI site, it was apparent that HSSs are present near the Nespas and Gnasxl promoter regions (sites IX and X) specifically on the paternal allele, and that this hypersensitivity is detected in three tissues examined (brain and liver are shown in Fig. 5A). In addition, one HSS upstream of the Nespas promoter (site VI) is present on both alleles. Hybridization with downstream probe SX1 illustrates the greater DNase I sensitivity of the paternal versus the maternal allele, in addition, multiple DNase I cleavages are detected across the Gnasxl exon region with this probe, which we assume to be a property of the unmethylated paternal allele.
![]() View larger version (55K): [in a new window] |
FIG. 5. DNase I hypersensitivity analysis of the Nespas-Gnasxl region in mouse tissues. (A) Analysis in nuclei from adult mouse brain and liver, from (B6 x M. spretus) hybrids, in which the paternal allele is M. spretus, and (M. spretus x B6) backcross, with maternal M. spretus. Nuclei were digested with DNase I (as in Fig. 4), purified DNA was digested with ScaI-XbaI and electrophoresed, and Southern blots were hybridized with probe SX2 or SX1. Maternal (M) and paternal (P) alleles are distinguished as a ScaI RFLP. Points on the left of the blots represent DNA markers formed by ScaI (Sc)-cut DNA digested with XbaI (X), EcoRI (R), HindIII (H), or Eco32I (E). The M. spretus-specific ScaI fragment is given (*Sc). Location of these reference sites is shown in panel B. DNase I cleavages are indicated by arrows to the right of each blot, weaker cleavages in grey, more prominent sites in black. (B) Interpretation of DNase I HSSs. The ScaI-XbaI fragment analyzed is shown. The Gnasxl exon is given as the black box, with 5' untranslated region in grey; the Nespas exon 1 is shown as the open box. The locations of probes SX1 and SX2 are given as striped boxes. The DNase I HSSs mapped from these and additional blots (not shown) are indicated using the convention above. HSSs on the maternal or paternal allele are labeled M or P, respectively.
|
1.5 kb upstream of the exon in all ES cells, as well as in adult brain, and was interpreted as a constitutive, biparental site. |
|
|---|
The Gnas cluster contains two gametic methylation imprints.
As a first expectation, it might be assumed that a single ICR could suffice for the imprinting of a compact cluster such as Gnas. At more extended imprinted gene clusters, there are examples of more than a single germ line DMR/ICR. The distal Chr 7 imprinting cluster in mouse, and the homologous 11p15.5 region in humans associated with the Beckwith-Wiedemann syndrome, each divide into two domains with separate ICRs: at H19, to regulate imprinting of Igf2 (48); and at KvDMR1, to regulate at least six linked imprinted transcripts (11). The former is a paternally methylated and the latter a maternally methylated gametic mark (50, 59). These elements are separated by
750 kb, and human patient data as well as knockout work in the mouse indicate that the two elements operate as distinct ICRs to regulate nonoverlapping sets of imprinted transcripts (39). Equally, however, there are other imprinted regions that appear to be regulated by a single ICR associated with a single gametic methylation mark. At Igf2r an intronic DMR is required for monoallelic expression not only of the Igf2r promoter but also two genes 110 to 155 kb downstream (62).
Our finding that the Nespas-Gnasxl DMR is a gametic imprint, together with the earlier identification of a germ line DMR at Gnas exon 1A (29), therefore predicts that the locus could contain two ICRs and is divided between separate domains regulated by independent imprinting mechanisms. This possibility will need to be tested by targeting experiments in the mouse. In the meantime, support for this contention comes from imprinting anomalies of human GNAS encountered in the disorder PHP1b. PHP1b results from loss of GS
expression in those tissues in which expression is strictly from the maternal allele. A consistent finding in PHP1b is that the GNAS exon 1A DMR is unmethylated on both alleles (3, 28). Loss of imprinted methylation at exon 1A is accompanied in most patients (sporadic and familial) by normal monoallelic methylation at the NESP55 and XL
S/AS DMRs (AS is the human equivalent to Nespas), indicating that methylation at the NESP55-XL
S/AS domain can be set independently of events at exon 1A. In a few patients, however, biallelic methylation at NESP55 is seen, with or without loss of maternal methylation of XL
S/AS (3, 28), implying that in some circumstances epigenotype can be regulated in concert across the whole locus.
A second informative human condition is hydatidiform mole. Complete hydatidiform moles are normally sporadic androgenetic conceptuses, which develop without an oocyte-derived chromosome complement. In biparental complete hydatidiform mole, molar pregnancies are recurrent and the disorder is thought to arise from a defect in setting up the methylation of imprinted genes in the mother's germ line. In biparental complete hydatidiform moles analyzed to date, both the GNAS exon 1A DMR and AS DMR were fully unmethylated (therefore biallelically), a finding compatible with a germ line origin of either or both their normal methylation patterns, while the NESP55 DMR was fully methylated (22). It was interpreted that NESP55 is a secondary DMR, whose methylation is dependent upon absence of methylation at AS. This is consistent with our findings in mouse gametes, as well as in ES cells in which methylation was present on maternal alleles at Nespas-Gnasxl but Nesp was unmethylated.
Properties of ICRs and methylation signals. Germ line DMRs are CpG rich elements which fulfill the criteria for CpG islands, except for the unique property of acquiring methylation in the male or female germ line. What features are responsible for this methylation and for germ line selectivity? The Nespas-Gnasxl DMR seems to be typical of many gametic imprints in its association with direct repeats, although the repeat region itself was not methylated in all oocytes. Such repeats have been proposed to attract methylation de novo through possible formation of unusual DNA structures (32, 40). Members of the Dnmt3 family are required for methylation at DMRs in oocytes (6, 15), but whether they respond to such repeats has not been shown. The most persuasive evidence for the involvement of direct repeats comes from Rasgrf1, whose paternal-specific expression is controlled by a remote DMR whose sperm-derived methylation depends upon a direct repeat block (60). In contrast, a GC-rich repeat in H19 appears to be dispensable for imprinting (38, 49). These pertain to paternal methylation, and for maternally methylated DMRs decisive tests through targeting have not been reported, although the direct repeat region of the U2af1-rs1 DMR appears to be superfluous (46). In the context of a synthetic transgene construct, the Igf2r/Air DMR can act as an imprinted methylation signal, this activity residing in the direct repeat region (40). However, the Igf2r/Air DMR imprints very inefficiently as a transgene in its own right (42), suggesting that elements in addition to the direct repeats are necessary for attracting or maintaining methylation allele-specifically. Direct repeats may facilitate the spread of methylation from surrounding regions undergoing methylation de novo in the germ lines. At the Nespas-Gnasxl DMR, as well as those at U2af1-rs1, Rasgrf1 and Grb10 (1, 41), methylation was present immediately upstream of the DMR in both germ lines. It will be important to track the time of appearance of methylation during germ cell development to see whether methylation of the DMR occurs in concert with that of the surrounding sequences, which would favor a spreading model. Alternatively, the spread of methylation may be limited in one germ line by factors bound to the DNA. The 5' border of the Nesp-Gnasxl DMR is associated with prominent HSSs on the sperm-derived allele in ES cells, but whether such HSSs are also present specifically during male gametogenesis and could provide an impediment to methylation remains to be shown.
At most ICRs studied to date, parental-allele-specific methylation is accompanied by parental-allele-specific chromatin features. These include DNase I HSSs likely to reflect the binding of methylation-sensitive nonhistone proteins on the unmethylated allele, methyl-cytosine binding proteins on the unmethylated allele, as well as differences in histone modifications (10, 12). Only in the case of the Angelman syndrome region is there evidence for an ICR demarcated by differential chromatin organization in the absence of differential methylation (35). Both DMRs in the Nesp-Gnasxl domain were associated with HSSs in adult tissues, most of which were present constitutively on the unmethylated allele. We looked specifically for HSSs in ES cells, because these cells represent an early embryonic state, a stage at which the maintenance of differential methylation against the genome-wide changes in methylation is critical, and because the availability of androgenetic and parthenogenetic ES cells allowed us to examine maternal and paternal alleles separately. Furthermore, studies at H19 revealed that pivotal elements in the imprinting of the locus are recognized as prominent HSSs in these cells (25, 47). Strong HSS were indeed identified at the Nespas-Gnasxl DMR, which flanked the putative start site for Nespas, and were present specifically on paternally derived, unmethylated alleles. It will be an important future aim to determine the identity of the protein-DNA interactions at these sites.
Models for imprinted expression of Nesp and Gnasxl. Collecting our observations together, we consider two models for the reciprocal imprinting of the Nesp and Gnasxl promoters (Fig. 7). In the first model, the Nespas-Gnasxl DMR is a gametic mark and represents a boundary, in the manner of the H19 insulator (4, 14). The essential elements of the model are: the Nesp and Gnasxl promoters have similar expression profiles (predominantly neuroendocrine tissues) (21, 23); the germ line DMR covers the Gnasxl promoter and the boundary element; the boundary is predicted to operate in a methylation-sensitive fashion. On the paternal allele, binding of insulator factors occurs which limit the action of downstream enhancers (yet to be characterized) to the unmethylated Gnasxl promoter. The Nesp promoter is thus quiescent and methylation ensues secondary to promoter inactivity. On the maternal allele, the Gnasxl promoter is silenced by methylation, the boundary fails to establish and downstream enhancers are free to interact with the unmethylated Nesp promoter. Relevant to this model, particularly prominent HSSs at the Nespas/Gnasxl DMR were found in ES cells but not adult tissues. This might suggest that the predicted boundary is functional in early embryonic cells, but becomes redundant once methylation at Nesp and Gnasxl is firmly set up so that differential enhancer access is controlled directly by the robust methylation and accompanying chromatin changes at the promoters. The methylation-sensitive boundary elements of the H19 and Peg3 imprinted genes comprise reiterated binding sites for the multifunctional DNA binding factors CTCF and YY1, respectively (4, 14, 26). In the Nesp-Gnasxl domain, we did not find similar clustered binding sites mapping within repeat arrays or within the HSSs. Isolated imperfect matches to CTCF and YY1 binding motifs were found throughout the mouse Nesp-Gnasxl domain, but these were not conserved in the human sequence, except for a pair of putative CTCF binding sites at the Nesp and NESP55 promoter regions (data not shown). Therefore, the nature of the factors at a hypothetical boundary would need to be investigated further. A second model draws analogies from the Igf2r locus and regulation by the antisense transcript Air (56). In this scenario, the Nespas-Gnasxl DMR is a unidirectional, cis-acting silencer on the unmethylated paternal allele, being the start site for the paternally expressed Nespas transcript antisense to Nesp. Expression of Nespas prevents expression of the Nesp promoter in cis, possibly by organization of a repressive chromatin structure or by inducing methylation. One difficulty with the antisense model is the imperfect concordance between sites of Nespas and Nesp expression (2, 27). These studies examined midgestation embryos or adult tissues, however, and do not exclude a model in which Nespas expression is required specifically at early stages to establish monoallelic expression and initiate permanent silencing of the Nesp promoter (via methylation and chromatin changes), after which Nespas may be redundant. Functional tests using gene targeting and other assays will need to be done to differentiate among these and alternative models.
![]() View larger version (23K): [in a new window] |
FIG. 7. Epigenetic features and models of reciprocal imprinted expression of Nesp and Gnasxl. (A) Summary of epigenetic features. The Nesp, Nespas and Gnasxl exons are depicted as before. Above the line, the approximate locations of DNase I HSSs are marked by the vertical arrows, with thickness of arrow reflecting prominence of HSS, and are shown for maternal (mat) and paternal (pat) alleles in adult tissues and ES cells. Below the line, methylation status is summarized, in both gametes and embryonic tissues, with methylated regions (+) and unmethylated regions (-) marked. (B) Enhancer-boundary model for monoallelic expression of Nesp and Gnasxl. Methylated promoters are indicated (mmm). The open circle represents the position of hypothetical enhancers able to control the Nesp and Gnasxl promoters, and the diamond a methylation-sensitive boundary. Active promoters are indicated by the black arrows. (C) Antisense model. Expression of the Nespas antisense is indicated by the wavy line. See text for description of the models.
|
We are sincerely grateful to J. Richard Chaillet (Rangos Research Center, Pittsburgh) for advice in setting up bisulfite sequencing analysis from mouse germ cells. We thank Simon Ball (Harwell) for M. spretus stocks and Nick Allen (Babraham) for uniparental ES cell lines. We also thank Rachel Smith for helpful comments on the manuscript.
P.A. is supported by a Marie Curie Individual Fellowship from the European Community Programme in Human Potential (under contract HPMF-CT-2001-01122); C.C. and E.A.C. were supported by studentships from the MRC and BBSRC. G.K. is a senior fellow of the MRC.
Present address: Johns Hopkins University, Baltimore, MD 21231. ![]()
|
|
|---|
gene GNAS1 in the pathogenesis of acromegaly. J. Clin. Investig. 107:R31-R36.
s is a new type of G protein. Nature 372:804-809.[Medline]
gene: predominant maternal origin of transcription in human thyroid gland and gonads. J. Clin. Endocrin. Metab. 87:4736-4740.
-subunit (Gs
) knockout mice is due to tissue-specific imprinting of the Gs
gene. Proc. Natl. Acad. Sci. USA 95:8715-8720.
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»