|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 283, Issue 27, 18612-18620, July 4, 2008
Genomic Matrix Attachment Region and Chromosome Conformation Capture Quantitative Real Time PCR Assays Identify Novel Putative Regulatory Elements at the Imprinted Dlk1/Gtl2 Locus*
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
13 kb upstream of the Gtl2 gene. This region is highly methylated on the paternal chromosome and unmethylated on the maternal chromosome (15). This methylation of the IG-DMR is therefore associated with activity from Dlk1 and repression of Gtl2. The pattern of differential methylation originates from the gametes and is maintained in all cell types during the entire life of the individual (primary imprinting mark). The significance of the IG-DMR was demonstrated by its deletion in the mouse (16). In the embryo, the deletion of the IG-DMR (
IG-DMR) on the maternally inherited chromosome leads to loss of imprinting (activation) of the normally repressed paternally expressed genes (Dlk1, Rtl1, and Dio3) and to the repression of the maternally expressed noncoding RNAs. However, the latter effect is not observed in the placenta to the same extent (14, 16). Furthermore, upon paternal transmission of the deletion, imprinting is unaltered. These findings demonstrate that the unmethylated IG-DMR acts as a repressor for imprinted protein coding genes on the maternal chromosome over a region of
1 Mb. To better understand how such long range regulation may occur at this locus, we investigated MAR activities and performed chromosome conformation capture (3C-qPCR) assays (17) in both the embryonic and the postnatal mouse liver.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Quantitative Analysis of MAR/HRS Assays
In each HRS assay, target DNA sequences were quantified by real time PCR using a SYBR Green mix (18) and a LightCycler apparatus (Roche Applied Science; software version 3.5). Amplifications were performed as previously described (3), and the enrichments were calculated as the ratio of the amount of target DNA in the MAR fraction versus the loop fraction. To standardize each assay, the ratio values were normalized against the ratio obtained for a negative control (NC) located 41 kb upstream of the Gtl2 gene and lacking predicted MAR features (see Fig. 2B). For primer sequences used for quantifications, please see supplemental data.
Quantification of the Allele Proportions
To assess the allelic specificity of the IG-DMR attachment, we took advantage of its parental origin-specific methylation pattern and the presence of restriction sites for methylation-sensitive/dependent enzymes. The following real time PCR primers were designed on either side of a restriction site for the methylation-dependent McrBC enzyme, which corresponds to a HhaI restriction site that was previously reported to be fully methylated on the paternal allele and unmethylated on the maternal (15): 5'-ACTCCTGGAGTGAGGGAAGG-3' and 5'-CAGCTAACCTGAGCTCCATGC-3'. The relative allele proportions in the DNA samples were determined by real time PCR quantifications performed before and after McrBC digestion. The reliability of the method was first demonstrated by using genomic DNA samples of wild-type mice (expected IG-DMR methylation level: 50%) and of
IG-DMR mice (16) upon both maternal (expected IG-DMR methylation, 100%) and paternal transmission (expected IG-DMR methylation, 0%) (data not shown).
To assess the allelic specificity of the Gtl2-MAR attachment, we performed MAR assays on nuclei from reciprocal hybrid mice obtained by mating (C57BL/6JxCBA) F1 mice (Mus musculus domesticus) with the JF1 strain (Mus musculus molossinus). The following primers for real time PCR were designed on either side of a BglII polymorphic restriction site absent in M. musculus domesticus: 5'-GCAGTTCTTTGCTCCACTGAAG-3' and 5'-ATCTTGACACGGGGAAGAG-3'. Amplifications were equally efficient on M. musculus molossinus and M. musculus domesticus alleles (i.e. no PCR biases were detected) (data not shown).
The relative proportions of both alleles were deduced from quantifications of DNA digested or undigested with BglII (see Refs. 2 and 3).
To assess MAR4 allelic attachment, a primer pair, specific to M. musculus molossinus, was designed: 5'-GTGCCTCCTCATCTTCTCCC-3' and 5'-GTGCGCTCACTGGGTGACAA-3'. This primer pair does not significantly amplify MAR4 from the M. musculus domesticus alleles. The allelic attachment of MAR4 was deduced from the M. musculus molossinus MAR4 attachment in MAR assays performed on reciprocal hybrid mice.
Luciferase Reporter Assay
DNA fragments containing the MAR4 and F9 sequences were PCR amplified on genomic DNA with the following primer pairs: MAR4 5'-AAGAGGGTGCCTCAGTCTGG-3' and 5'-TTGGGTTCCACCCTCAGTGG-3'; F95'-GGAATGCCACCCACAAGGAG-3' and 5'-TCAGCCTCCGGAACAGCTAG-3'. The resulting 2248- and 2720-bp fragments were cloned in both directions in the pDrive vector (Qiagen) and subsequently cloned into the "pGL3promotor" vectors (firefly luciferase under the SV40 promoter) using XhoI and KpnI restriction sites. MAR4 constructs were then subjected to mutagenesis by using the QuikChange site-directed mutagenesis kit (Stratagene) with the following primers: 5'-GCGAAGATCTCCCAGACGTTGTGCTAAGACGAGG-3' and 5'-CCTCGTCTTAGCACAACGTCTGGGAGATCTTCGC-3'. The resulting
MAR4 constructs contain a 2169-bp insert in the sense or antisense directions with a 79-bp deletion that removes the A/T-rich MAR4 sequence from the constructs. Primary bipotential mouse embryonic liver cells, derived from E14.5 embryos, were transfected with these reporter plasmids, together with the Renilla luciferase (Rluc) control reporter vector pRL-TK vector (Promega) using the jetPEI transfection reagent (PolyPlus Transfection) according to the supplier's protocol. Luciferase activity was determined with a dual luciferase reporter assay system (Promega).
3C-qPCR Assays
3C-qPCR assays were performed as previously described (17). Briefly, the assays were performed on nucleus preparations extracted from E15.5/16.5 or D30 mouse liver samples. 1 x 107 nuclei were formaldehyde cross-linked, digested using EcoRI, and ligated as previously described (17). DNA was then decross-linked (overnight at 65 °C) and purified by classical phenol extraction procedures. An additional digestion using BamHI was then performed on each sample to allow more accurate quantifications (see Ref. 17). Real time PCR was performed on a LighCycler apparatus (Roche Applied Science) using the TaqMan technology (QuantiTect Probe PCR Master Mix; Qiagen). We used a 5'FAM-3'BHQ1 oligonucleotidic probe (MWG-Bioctech). For primer sequences used for quantifications, please see the supplemental data. 3C-qPCR data processing was as follows.
Normalization 1: GAPDH Normalization ("Loading Control")
As previously described (17), Ct obtained for each chimerical ligation fragments were processed using parameters of a standard curve (slope and intercept) to obtain quantification values that were normalized to a "GAPDH loading control."
Normalization 2: Noise Band Determination and Normalization of Data to Basal Interaction Level
This normalization compensates for experimental variations and allows comparison between different 3C assays; it replaces the Pdhb or Ercc3 normalizations usually used, (see Ref. 17). The procedures below (Parts A, B, and C) should be followed independently for each sample.
Part A: Removal of "Deviant" Experimental Points—Important: removal can only be done when at least three independent 3C assays have been performed for a given sample. In the present study, three independent assays were performed for both the E16.5 liver and the D30 liver samples. When only one or two independent assays have been performed, select all points and go to Part B of the procedure.
For Step A.1, for each experimental point, calculate the Log of the values (v) normalized to GAPDH (Normalization 1), [Log v].
For Step A.2, for each fragment (fx), calculate the mean [m(fx)] of [Log v].
For Step A.3, for each fragment (fx), calculate the standard deviation [sd(fx)] of [Log v].
For Step A.4, for each experimental points corresponding to a given fragment (fx), calculate x(fx) = m(fx) – (sd(fx)*k), where k is the tolerance factor that we systematically fix at 1.05.
For Step A.5, for each experimental point corresponding to a given fragment (fx), calculate y(fx) = m(fx) + (sd(fx)*k), where k is the tolerance factor that we systematically fix at 1.05.
For Step A.6, for each fragment (fx), select all experimental points for which x(fx) < [Log v] < y(fx). Remove all other points.
Part B: Determination of the Basal Interaction Level—For Step B.1, for each experimental point selected above (Part A) corresponding to a given fragment (fx), calculate the mean of values (v) normalized to GAPDH (normalization 1) [m'(fx)].
For Step B.2, for each experimental points selected above (Part A) corresponding to a given fragment (fx), calculate the standard deviation of the values normalized to GAPDH v (normalization 1) [sd'(fx)].
For Step B.3, calculate the mean [M] of all the [m'(fx)] values.
For Step B.4, calculate the mean [SD] of all [sd'(fx)] values.
For Step B.5, for each fragment (fx), select all experimental points for which m'(fx) < M – (SD).
For Step B.6, calculate the mean [M1] of all the [m'(fx)] values of the points selected in Step B.5.
For Step B.7, calculate the standard deviation [SD1] of all the [m'(fx)] values of the points selected in Step B.5.
For Step B.8, for each fragment (fx), select all experimental points for which m'(fx) > M1 + (SD1).
For Step B.9, calculate the mean [M2] of all the [m'(fx)] values of the points selected in Step B.8.
For Step B.10, calculate the standard deviation [SD2] of all the [m'(fx)] values of the points selected in Step B.8.
The value of M2 is the raw basal interaction level that we then use to normalize our data (see Part C, below). M2 ± SD2 is the raw noise band.
Part C: Normalization to the Basal Interaction Level and Determination of the Noise Band—For Step C.1, for each fragment (fx), calculate the normalized mean M(fx) = m'(fx)/M2.
For Step C.2, for each fragment (fx), calculate the normalized standard deviation SD(fx) = sd'(fx)/M2.
For Step C.3, calculate the normalized basal interaction level = M2/M2 = 1.
For Step C.4, calculate the normalized noise band (NB) = SD2/M2. The basal interaction level ± NB is the normalized noise band.
For Step C.5, make a graph showing the distribution of the M(fx) values ("relative cross-linking frequencies") as a function of the distance (in kb) between the fx fragments and the "constant fragment."
| RESULTS |
|---|
|
|
|---|
160 kb) than the equivalent mouse locus (
110 kb), both the number and relative positions of the predicted MARs relative to the IG-DMR and the Dlk1 gene are conserved between both species (Fig. 1). Altogether, the situation is reminiscent of MARs previously identified at the Igf2/H19 locus (3).
|
DNA fragments containing the predicted MARs were then investigated. The fragments containing the MAR1, MAR2, MAR3, and the Gtl2-MAR were found to be enriched in the MAR fraction in both the E15.5 and 30-day-old mouse liver (Fig. 2B). We can thus conclude that their attachment is not linked to gene expression. In contrast, MAR4 attachment appears strongly correlated with gene expression (p = 0.0384). Indeed, its attachment level, which is very high in the E15.5 liver when the genes are strongly expressed, falls to background level in the 30-day-old mouse liver, when the Dlk1 and Gtl2 genes are repressed.
We also assessed the enrichment of DNA fragments containing the previously reported conserved sequences (cs9, 10/11, 12, 13, and 14) as well as the differentially methylated regions (Dlk1-DMR, IG-DMR, and Gtl2-DMR) (15, 19). The conserved sequences are not or are only weakly enriched into the MAR fraction, with exception of the cs10/11 fragment, which displays an attachment pattern similar to the neighboring MAR4 (high attachment in the E15.5 liver; p = 0.0384). The attachment pattern of the Dlk1-DMR fragment, which essentially corresponds to exonic sequences, seems to follow the Dlk1 gene transcription pattern. Therefore, as previously shown for other actively transcribed regions (20), attachment of the Dlk1-DMR fragment may thus be driven by the transcription machinery itself. The fragment containing the Gtl2-DMR (F12) was also found enriched in the nuclear matrix fraction. However, the Gtl2-DMR corresponds with the Gtl2 promoter, and it contains several conserved sequences (cs16 to 19), and therefore it is difficult to hypothesize which feature is responsible for the nuclear matrix attachment.
The IG-DMR Is Retained within the Nuclear Matrix Fraction—Surprisingly, the fragment containing the IG-DMR was also found to be strongly and persistently attached in the mouse liver (Fig. 2B). After performing experiments at a higher resolution, we were able to isolate the attachments of the IG-DMR and of the neighboring Gtl2-MAR into two distinct restriction fragments separated by several other weakly attached fragments (data not shown).
|
MAR4 construct) (Fig. 3). We conclude that the 79-bp A/T-rich MAR4 fragment contains intrinsic directional activity that favors transcription. MAR4 and Gtl2-MAR Attachments Occur on Both Parental Chromosomes—Because MAR4 attachment correlates with expression and because both Dlk1 and Gtl2 genes are monoallelically expressed, it was of interest to determine whether one or both parental copies is recovered from the nuclear matrix. To address this, we used M. musculus domesticus x M. musculus molossinus hybrid mice and designed two real time quantitative PCR primer pairs. The first one is able to detect MAR4 independently of the genotype, i.e. from both domesticus and molossinus alleles (dom+mol), whereas the second primer pair is specific for the molossinus alleles (mol) (Fig. 4A). Only the (dom+mol) primer pair showed enrichment into the MAR DNA fraction. We conclude that genomic MAR assays retain the MAR4 sequence exclusively from the M. musculus domesticus chromosome. This domesticus-specific finding is most likely a genetic background-specific phenomenon, because this attachment is similar in a (dom x mol) cross compared with the reverse cross (mol x dom) (Fig. 4A). We therefore conclude that attachment of the domesticus MAR4 allele can occur on both chromosomes and does not depend on its parental origin.
We next examined the allelic attachment of the Gtl2-MAR. Here, no allele-specific primer pair could be successfully designed; however, a BglII restriction site specific for the M. musculus molossinus allele was identified. We thus used a previously described technique that allows determination of allele ratios using such polymorphic restriction sites and the real time quantitative technique (2). After assessing for the absence of any PCR bias (data not shown), we showed that the M. musculus domesticus allele was preferentially retained in our assays (Fig. 4B). However, again, this preferential attachment was found to be independent of its parental origin because both the (dom x mol) and (mol x dom) crosses display similar allele ratios (Fig. 4B). We conclude that the Gtl2-MAR is equally attached on both parental chromosomes.
Parental Imprinting Influences IG-DMR Attachment—Finally, it was particularly important to determine the relative allelic MAR attachment level of the IG-DMR sequence because this sequence displays allele-specific features characteristic of imprinting control regions. This sequence is hypermethylated on the paternal chromosome and completely unmethylated on the maternal copy (15). We used this imprinted allele-specific methylation pattern as a means to determine the parental origin of the alleles. We applied the DNA methylation-dependent McrBC restriction enzyme and determined the ratio between the methylated versus unmethylated alleles that correspond to the allele ratio (2). Therefore, for this experiment, classical inbred mouse strains (M. musculus domesticus) could be used instead of hybrid mice, thus precluding any potential PCR or genetic background related bias. We first checked that, in the loop fraction of genomic MAR assays, this region possesses a mean methylation level of 50%, as a result of its monoallelic methylation pattern (Fig. 4C, left panel). We next examined the MAR fraction and found that the methylated paternal allele is strongly enriched in this fraction (76.2% ±3.1) (Fig. 4C, right panel). Because the signal observed in Fig. 2B for the IG-DMR sequence is four times higher than the control sequence (NC), the remaining signal observed on the maternal allele (23.8% ±3.1) can be entirely attributed to nonspecific background. Therefore, the IG-DMR sequence is attached exclusively on the methylated paternal allele.
|
|
15–20 kb on each side of the MAR4, prevented detection of any specific interactions that may occur close to this sequence. However, from 20 to 80 kb away from MAR4, long range interaction frequencies appear to stabilize to a basal level that corresponds to random interactions between the fragment containing the MAR4 and other fragments throughout the locus. Therefore, for each sample (E16.5 or D30 liver), a mean basal interaction level covering the whole locus, was calculated and fixed to a value of 1. The mean of the experimental variations of the points used to calculate this value was then determined and used to define a "noise band," within which cross-linking frequencies can be considered as reflecting random interactions (see "Experimental Procedures"). Local peaks observed outside the side effect zone and above the noise band represent favored (nonrandom) interactions. Weak interactions peaks were found between the restriction fragment containing the MAR4 and fragments –6, –3, and +7 in the E16.5 liver and fragments –10, –6, –3, and +9 in the D30 liver. Interestingly, MAR4 interacts with the 5' part of the IG-DMR (fragment +7) in E16.5 liver, but not in D30 liver. Conversely, interactions with fragments –10 and +9 (containing the Gtl2-MAR) are only observed in D30 liver. However, in the absence of any functional data about these unexpected interactions, their potential roles remain unclear. Finally, no interaction was found with the fragments containing the Dlk1/Gtl2 promoters (fragments –8 and +10, respectively). Furthermore, in the E16.5 liver, interaction frequencies with these fragments (as well as with fragment –10) even fall below the noise band (see inset in Fig. 5). This means that the frequencies of random interactions between the MAR4 sequence and the Dlk1/Gtl2 promoters are lower than those observed for the other fragments of the locus. Noticeably, this effect is not observed in D30 liver when the genes are repressed, and it may thus simply reflect local steric impairment around the promoters because of the high transcriptional activity in E16.5 liver.
|
IG-DMR deletion (16). When the deletion was paternally inherited (Fig. 6A), no interaction peak was detected between the fragments containing the MAR4 and the remaining maternally inherited IG-DMR (fragment 7); the interaction level falls to the basal level of random interaction. In contrast, when the deletion was maternally inherited, the interaction is kept to the same level as observed for wild-type mice. The measured value of the relative cross-linking frequency for fragment 7 is 4 for both wild-type and maternal heterozygous mice (Fig. 6B). We conclude that the MAR4/IG-DMR interaction occurs exclusively on the methylated paternal allele. Upon maternal inheritance of the deletion, novel interaction peaks are also observed between the MAR4 and fragments 6 and 9 (containing the Gtl2-MAR). Because such interactions are not detected upon paternal inheritance or even in the wild-type mouse, we conclude that, when maternally inherited, the deletion chromosome adopts a peculiar genomic organization where the MAR4 is not only interacting with the IG-DMR but also with other regions including the Gtl2-MAR.
| DISCUSSION |
|---|
|
|
|---|
MAR sequences have been operationally defined in in vitro assays as being able to attach to a purified nuclear matrix or scaffold. Therefore, most analysis of genomic MAR assays have focused on classical A/T-rich MAR sequences. However, after having adapted the sensitive real time PCR technology to the analysis of genomic MAR assays, it now clearly appears that not only such classical A/T-rich MAR sequences can be retained but also other very important regulatory sequences like enhancers (3) or G/C-rich differentially methylated regions/imprinting control regions (this work). In genomic MAR assays, such genomic sequences have the potential to be trapped (in a tissue- or developmental stage-specific manner) into the so-called nuclear matrix upon high salt treatment of nuclei and can thus be recovered. Therefore, we propose a strict and wider operational definition and to rename such sequences the HRSs.
|
Altogether, these findings lead us to conclude that, just like the DNase sensitivity assays, the HRS assays reveal epigenetic features of the chromatin. Because the only feature in common with all the HRS identified so far is that they appear to act at a long distance (several tens to a few hundred kb), we propose that the HRS assays most likely reflect epigenetic features involved in higher order chromatin organization of the mammalian genome. Because of technical limitations, the organization of the mammalian genome at that scale remains largely unknown. However, a recent technological breakthrough allowed the identification of chromatin loops of several hundred kb in the mammalian genome (21). Such chromatin loops appear to be linked to gene expression, and therefore it was suggested that specific higher order chromatin architectures should be associated with gene activity or repression. Furthermore, it becomes increasingly clear that A/T-rich sequences and factors binding such regions play important roles in tissue-specific higher order organization (22). Therefore, MARs, as other HRS, appear as versatile regulatory elements that could combine with differentially methylated regions and other imprinting control regions to confer cell-specific higher order chromatin organization and to control monoallelic expression (3, 23). As proposed previously, it may well be that such higher order chromatin architecture corresponds to a real "genome format" that specifies the transcriptional status of the genes in a tissue- and developmental stage-specific manner (24).
Interestingly, we also identified a 79-bp A/T-rich sequence that we called MAR4 as a novel putative regulatory element that possesses an intrinsic activity that favors transcription. This element is located 9136 bp upstream from a conserved dodecamer motif that contains a single nucleotide polymorphism causing the Callipyge phenotype in sheep (29, 33). Interestingly, the Callipyge mutation was recently shown to enhance bidirectional long range Dlk1-Gtl2 intergenic transcription in cis (25). Therefore, because the novel MAR4 activity appears to be stronger when the MAR4 is inserted in the antisense direction (with respect to the endogenous locus) in a reporter construct, its action may potentially favor transcription of antisense transcripts. Antisense transcripts have been found at numerous loci, and they play important roles in several imprinting mechanisms (26–28); however, such antisense RNAs are very weakly expressed and, although they have been described in the sheep (CLPG1 transcript) (25, 29), they remain to be identified at the Dlk1/Gtl2 locus in mouse. We also show that the MAR4 sequence does not physically interact with the Dlk1/Gtl2 promoters. This element could therefore be similar to the DMR2/MAR2 sequence that maintains high Igf2 transcription levels on mouse chromosome 7 (30). The DMR2/MAR2 was not found to interact with the Igf2 promoters but was rather proposed to act through a paternal-specific interaction with the imprinting control region (31) and endodermic enhancers (3). Similarly, we could propose that the MAR4 element, which interacts with the IG-DMR on the paternal allele (Figs. 5 and 6), may act at the endogenous Dlk1/Gtl2 locus by recruiting distant enhancers. However, none of the 3C assays performed at the MAR4 sequence revealed interaction peaks that may correspond to such interactions. Therefore, we speculate that sequestering of the IG-DMR by MAR4 may help to ensure the activity of the Dlk1 gene on the paternal chromosome. Alternatively, paternally specific MAR4/IG-DMR interaction might contribute to the repression of MAR4-mediated transcriptional activation of the Gtl2 gene on that chromosome, thus contributing to the process of genomic imprinting at this locus. Inactivation of MAR4 by homologous recombination in the mouse is now required to prove the functional role of this element in vivo.
| FOOTNOTES |
|---|
The on-line version of this article (available at http://www.jbc.org) contains supplemental Tables 1 and 2. ![]()
1 Supported by Association pour la Recherche contre le Cancer Fellowship JR/MLD/MDV-P05/2 and P06/2. ![]()
2 To whom correspondence should be addressed: IGMM, UMR5535 CNRS-UMII, IFR122, 1919, Route de Mende, 34293 Montpellier Cedex 5, France. Tel.: 33-467-61-36-84; Fax: 33-467-04-02-31; E-mail: forne{at}igmm.cnrs.fr.
3 The abbreviations used are: MAR, matrix attachment region; 3C, chromosome conformation capture; qPCR, quantitative real time PCR; DMR, differentially methylated region; IG, intergenic germ line-derived; HRS, high salt recovered sequence; NC, negative control; En, embryonic day n; Dn, postnatal day n; GAPDH, glyceraldehyde-3-phosphate dehydrogenase. ![]()
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| All ASBMB Journals | Molecular and Cellular Proteomics |
| Journal of Lipid Research | ASBMB Today |