MacroH2A1 knockdown effects on the Peg3 imprinted domain
© Choo et al; licensee BioMed Central Ltd. 2007
Received: 02 July 2007
Accepted: 31 December 2007
Published: 31 December 2007
MacroH2A1 is a histone variant that is closely associated with the repressed regions of chromosomes. A recent study revealed that this histone variant is highly enriched in the inactive alleles of Imprinting Control Regions (ICRs).
The current study investigates the potential roles of macroH2A1 in genomic imprinting by lowering the cellular levels of the macroH2A1 protein. RNAi-based macroH2A1 knockdown experiments in Neuro2A cells changed the expression levels of a subset of genes, including Peg3 and Usp29 of the Peg3 domain. The expression of these genes was down-regulated, rather than up-regulated, in response to reduced protein levels of the potential repressor macroH2A1. This down-regulation was not accompanied with changes in the DNA methylation status of the Peg3 domain.
MacroH2A1 may not function as a transcriptional repressor for this domain, but that macroH2A1 may participate in the heterochromatin formation with functions yet to be discovered.
The nucleosome, the basic unit of eukaryotic chromatin, consists of a 146-bp DNA wrapped around a histone octamer composed of two copies of core canonical histones, H2A, H2B, H3, and H4. Eukaryotic chromatin also contains a small amount of atypical histones, which have different amino acid sequences than canonical histones, including CENPA, H3.3, H2A.X, H2A.Z, and macroH2A . Among these histone variants, the macroH2A variants display the most unusual protein structure: an N-terminal H2A domain is fused to a large C-terminal non-histone domain. This results in a histone that is nearly three times the size of conventional histones . The N-terminal third of macroH2A (H2A-like) shares 64% sequence identity with H2A, and the remaining two thirds of the protein shows similarity with the domain called 'macro' . MacroH2A is conserved throughout all vertebrate lineages, and two different members of the macroH2A family have been identified to date, macroH2A1 and macroH2A2 [2, 4]. The macro domain functions as a strong transcriptional repressor by inhibiting the initiation step of Pol II transcription and by interfering with p300-dependent histone acetylation . The macro domain is also capable of blocking the chromatin remodeling process mediated by SWI/SNF and ACF [5, 6].
MacroH2A is mainly associated with the heterochromatic regions of chromosomes. This is consistent with the observed repression activity. Notably, early immuno-fluorescence studies revealed the enrichment of macroH2A1 deposition in the inactive X chromosome of mammals [7, 8], which undergoes chromosome-wide repression to balance different gene dosage between female and male . MacroH2A1 is also associated with other regions of heterochromatin, including peri-centromeric regions and senescence-associated heterochromatic foci [10–12]. Another recent study revealed that high levels of macroH2A1 deposition are detected in methylated CpG islands that are located close to promoter regions . In particular, macroH2A1 is highly enriched in the inactive, methylated alleles of Imprinting Control Regions (ICRs), which are shown to be critical for maintaining the imprinting (allele-specific expression) of surrounding genes [14, 15]. Although the mechanism(s) targeting macroH2A1 to the ICRs are not well understood at moment, the allele-specific deposition of macroH2A1 in the ICRs represents another epigenetic marker besides DNA methylation that differentiates two parental alleles.
To characterize potential roles of macroH2A1 in genomic imprinting, the current study has analyzed several imprinted domains with a major focus on the Peg3 domain, an evolutionarily well-conserved domain located in human chromosome 19q13.4/proximal mouse chromosome 7 . This domain contains 6 imprinted genes within a 500-kb genomic interval, including paternally expressed Peg3, Usp29 and Zfp264, and maternally expressed Zim1, Zim2 and Zim3. We have generated several stable cell lines showing low levels of macroH2A1 using Neuro2A cells and analyzed changes in the expression levels of imprinted genes. The results indicate that a subset of imprinted genes, including Peg3 and Usp29, were affected by lowering the cellular levels of the macroH2A1 protein. Furthermore, the expression of these genes were further repressed, rather than de-repressed, in response to the knockdown of the potential repressor macroH2A1. This suggests that macroH2A1 may not function as the dominant repressor for the transcription of imprinted genes.
MacroH2A1 knockdown stable cell lines
We performed chromatin immunoprecipitation (ChIP) experiments using the two macroH2A1 knockdown cell lines to test whether macroH2A1 is indeed depleted in the Imprinting Control Regions (ICRs) of several domains, which show high levels of macroH2A1 deposition . According to several ChIP experiments using anti-macroH2A1 antibody, the deposition levels of macroH2A1 on ICRs in the two stable cell lines were much lower than the levels detected in the positive control Neuro2A cells (Fig. 1B). This confirms the depletion of macroH2A1 in most ICRs in the two macroH2A1-knockdown cell lines. The growth and morphology of these stable cell lines appeared to be normal without any obvious abnormalities.
MacroH2A1 knockdown effects on the transcription of imprinted genes
In the Peg3 domain, the expression of most resident genes except for Zim3 was detected in the Neuro2A-derived control and knockdown cell lines. The expression levels of Peg3 and Usp29 in the two macroH2A1 knockdown cell lines were lower than those in the two control cell lines (Fig. 2A). The results from qRT-PCR also derived a similar conclusion: the expression levels of Peg3 and Usp29 were decreased by 2.8- to 3.0-fold and 1.4- to 2.5-fold, respectively (Fig. 2B). This is an unexpected outcome given the prediction that macroH2A1 is a transcriptional repressor [5, 6]. However, the expression levels of Zim1 and Zfp264 did not deviate significantly from those of the control cells (Student's t test, 0.05 <p < 0.1), indicating that these two genes are probably not changed by macroH2A1 knockdown.
We also analyzed the expression levels of other imprinted genes, including Gtl2/Dlk1, Igf2/H19, Kcnq1/Kcnq1ot1, Nesp/Nespas, GnasXL, Exon1A, and Xist/Tsix (Fig. 2A and see Additional file 1). These genes were analyzed with RT-PCR with a fixed number of cycles, ranging from 30 to 38, after the normalization of cDNAs with the three internal controls (GAPDH, 28S, and β-actin). The expression levels of Gtl2, H19, Kcnq1, Kcnq1ot1, and Tsix were relatively low in the Neuro2A-derived cell lines, and thus we could not analyze the knockdown effects on these genes. The expression levels of Igf2, Nesp and GnasXL in the two knockdown cells were lower than those in the two control cells, indicating potential knockdown effects on these genes. However, due to inconsistent expression levels of these genes in one of the control cells (pSicoR in Fig. 2 and see Additional file 1), we could not confirm the knockdown effects on these genes with statistical significance. In contrast, Xist was highly expressed in Neuro2A-derived cell lines, but showed a similar level of expression among all the cell lines analyzed (Fig. 2A). This indicates that macroH2A1 knockdown had no obvious effect on this locus. The expression levels of DNA methyltransferases (Dnmts) were also analyzed in the macroH2A1 knockdown cells. Dnmt1, Dnmt3a and Dnmt3b were expressed at a similar level in both controls and macroH2A1 knockdown cells, showing that macroH2A1 knockdown had no effect (see Additional file 1). In conclusion, our analyses indicate that the expression levels of two genes located in the Peg3 domain were down-regulated in the macroH2A1 knockdown cells.
MacroH2A1 knockdown effects on the DNA methylation of imprinted genes
To determine the methylation status of promoter regions, we subsequently digested each PCR product with a second series of enzymes that can differentiate methylated and unmethylated DNAs, including Hpa I, Taq I, Hpy 188I, and Eco RV (the top panel in Fig. 3). Digestion of each amplified product with one of these enzymes indicates DNA methylation, whereas non-digestion indicates that the tested DNA is unmethylated. As shown in Fig. 3, all the tested promoter regions within the Peg3 domain derived an equal ratio of methylated to unmethylated DNAs in all the cell lines tested. This is an expected outcome given that only one of two parental alleles of these regions is usually methylated in vivo (thus named DMRs, Differentially Methylated Regions). However, we did not find any difference between the DNA methylation levels of the control and knockdown cell lines, suggesting that DNA methylation status of the promoters located in the Peg3 domain was not changed in the macroH2A1 knockdown cell lines. This is also true for the Xist promoter (Fig. 3). Overall, our DNA methylation analyses indicate that macroH2A1 depletion did not cause any change in the DNA methylation status of imprinted domains.
The current study demonstrates that lowering the macroH2A1 protein levels affects the expression levels of a subset of imprinted genes, including Peg3 and Usp29. Interestingly, the expression of these genes was repressed, rather than de-repressed, in response to the knockdown of the potential repressor macroH2A1. Subsequent analyses also revealed that macroH2A1 knockdown did not change the DNA methylation status of their promoter regions. This suggests that macroH2A1 may not function as a dominant repressor for the transcription of these imprinted genes.
The down-regulation of a subset of imprinted genes, Peg3 and Usp29, in macroH2A1 knockdown cells is an unexpected outcome since several lines of evidence strongly suggest a repressor role for this histone variant [5, 6]. According to the results from a previous study , macroH2A1 is deposited mainly on the inactive allele of the promoter region of Peg3 and Usp29, suggesting a potential role in the inactive allele of this imprinted domain. RNAi-based knockdown indeed depleted macroH2A1 from the inactive allele of Peg3 and Usp29 as shown in Fig. 1B. Despite this depletion, however, the inactive allele of the two genes was still methylated at similar levels between the control and knockdown cells (Fig 3). This suggests that macroH2A1 knockdown does not have any impact on the methylation levels and subsequently silent state of the inactive allele. Nevertheless, the expression levels of Peg3 and Usp29 were down-regulated in the macroH2A1-knockdown cells (Fig. 2). Since macroH2A1 is mainly associated with the inactive allele, the transcription status of which is unaffected in the knockdown cells, the observed down-regulation is likely caused by some changes that have occurred in the active allele. It is, therefore, most likely that the observed effect on Peg3 and Usp29 in the knockdown cells is an indirect outcome involving other unknown components at their active alleles.
Despite the high levels of macroH2A1 deposition in the ICRs of most known imprinted domains , only two imprinted genes, Peg3 and Usp29, are affected in the macroH2A1-knockdown cell lines. This unexpected outcome is also similar to the results from the macroH2A1-deficient mouse . Depletion of macroH2A1 appears to have affected the transcription of only a small number of genes, about 100 genes, although macroH2A1 has been predicted to be a global regulator. Out of 100 genes, some were also down-regulated, indicating another indirect effect similar to that seen in the macroH2A1-knockdown cells. Since YY1 is known to be involved in the transcription of the two imprinted genes, Peg3 and Usp29 [20, 21], we first analyzed the macroH2A1-knockdown effects on the expression levels of YY1 to obtain potential clues regarding the down-regulation of Peg3 and Usp29. However, we did not observe any change in the expression level of this potential candidate (see Additional file 1). We also examined the list of genes that were affected in the macroH2A1-deficient mouse . Among the 100 genes affected by macroH2A1 knockout, the transcription factor USF2 (Upstream Stimulating Factor 2) is one likely candidate that was shown to be down-regulated in the macroH2A1-knockout mouse. This transcription factor deserves to be studied further in the near future, since it is known to interact or compete with YY1 for binding to DNA elements . It is plausible that macroH2A1 knockdown might have down-regulated the transcription of Peg3 and Usp29 through other unknown transcription factors like USF2. Overall, the results described above suggest that macroH2A1 may be involved in the fine tuning of transcriptional control of a subset of genes, but not as a global regulator.
Mammalian imprinted domains are usually co-regulated through a small number of controlling regions, called ICRs [14, 15]. Similarly, within the Peg3 domain, the surrounding regions of Peg3's 1st exon, termed the Peg3-DMR, are predicted to control the whole domain (Fig. 2B). This potential ICR is characterized by high levels of macroH2A1 deposition as well as an unusual tandem array of YY1 transcription factor binding sites [20, 21]. Interestingly, independent knockdown experiments targeting macroH2A1 and YY1 both appear to have resulted in similar global changes in the transcriptional levels of most imprinted genes in this domain. The two imprinted genes of the Peg3 domain displayed a down-regulation in response to the macroH2A1 knockdown, although this down-regulation appears to be an indirect outcome. Also, lowering the cellular level of the YY1 protein has resulted in up-regulation of most of the imprinted genes within the Peg3 domain . Although opposite, these similar global responses are consistent with the prediction that the imprinting and transcription of this domain is linked together and regulated by a shared mechanism, possibly through the Peg3-DMR.
RNAi-based macroH2A1 knockdown experiments in Neuro2A cells changed the expression levels of a subset of genes, including Peg3 and Usp29 of the Peg3 domain. The expression of these genes was down-regulated, but not up-regulated, in response to reduced protein levels of the potential repressor macroH2A1. This down-regulation was not accompanied with changes in the DNA methylation status of the Peg3 domain. Thus, it is likely that macroH2A1 may not function as a transcriptional repressor for this domain, but that macroH2A1 may participate in the heterochromatin formation with functions yet to be discovered.
MacroH2A1 knockdown using siRNA technique
We designed four siRNA constructs for macroH2A1 knockdown. The sequences of these siRNA constructs are as follows: siRNA1, sense strand, 5'-TGTGCTCGCTTCGGCAGCACATATACTGGCCACCCTAAGTATAGGATTCAAGAGATCCTATACTTAGGGTGGCCTTTTTTC-3'; siRNA1, antisense strand, 5'-TCGAGAAAAAAGGCCACCCTAAGTATAGGATCTCTTGAATCCTATACTTAGGGTGGCCAGTATATGTGCTGCCGAAGCGAGCACA-3'; siRNA2, sense strand, 5'-TGTGCTCGCTTCGGCAGCACATATACTGAGACAACAAGAAGGGACGTTCAAGAGACGTCCCTTCTTGTTGTCTCTTTTTTC-3'; siRNA2, antisense strand, 5'-TCGAGAAAAAAGAGACAACAAGAAGGGACGGTCTCTTGAACGTCCCTTCTTGTTGTCTCAGTATATGTGCTGCCGAAGCGAGCACA-3'; siRNA3, sense strand, 5'-TGTGCTCGCTTCGGCAGCACATATACTACTGCTTGGCTCTAGCTGATTCAAGAGATCAGCTAGAGCCAAGCAGTTTTTTTC-3'; siRNA3, antisense strand, 5'-TCGAGAAAAAAACTGCTTGGCTCTAGCTGATCTCTTGAATCAGCTAGAGCCAAGCAGTAGTATATGTGCTGCCGAAGCGAGCACA-3'; siRNA4, sense strand, 5'-TGTGCTCGCTTCGGCAGCACATATACTTGATGAAGAGCTAAACCAGTTCAAGAGACTGGTTTAGCTCTTCATCATTTTTTC-3'; siRNA4, antisense strand, 5'-TCGAGAAAAAATGATGAAGAGCTAAACCAGTCTCTTGAACTGGTTTAGCTCTTCATCAAGTATATGTGCTGCCGAAGCGAGCACA-3'. Double stranded oligonucleotides were inserted into the Hap I and Xho I sites of the pSicoR vector . We purified these construct DNAs from E. coli using the Hispeed plasmid midi kit (Qiagen), and transfected using the GeneJuice transfection reagent (Novagen) according to the manufacturer's protocol. To derive macroH2A1 knockdown cell line, we co-transfected pcDNA 3.1/His/LacZ vector (Invitrogen) containing the neomycin resistant gene. Transfectants were selected by adding G418 (500 μg/ml; Calbiochem) to the culture medium. After transfection of Neuro2A cells with a pool of the four siRNA constructs, we performed PCR-based cloning and DNA sequencing analyses with DNAs derived from each stable cell line to determine which construct(s) was the most potent in terms of macroH2A1 knockdown. Our analyses indicated that the siRNA #3 and #4 constructs were much more potent than the siRNA #1 and #2 constructs.
Histone extraction and Western blot
Histones were extracted with 0.2N HCl from 2 × 106 cells after cell lysis with the triton exaction buffer (0.5% Triton X-100, 2 mM PMSF, 0.02% NaN3 in PBS). Protein concentrations were determined by Bradford assay kit (Pierce). Each sample of histone extracts (20 μg) was separated on 14% SDS-PAGE gels, transferred to the PVDF membrane (Hybond-P; Amersham) using a Mini Trans-Blot transfer Cell (Bio-Rad). Membranes were blocked for 1 hour in Tris-buffered saline (TBS-T) containing 5% skim milk and 0.05% Tween-20, incubated at 4°C overnight with anti-macroH2A1 antibody (Cat. # 07-219; Upstate Biotech.). These membranes were further incubated for 1 additional hour with the secondary antibody linked to horseradish peroxidase (Sigma). The blots were visualized with a Western blot detection system (Intron Biotech) according to the manufacturer's protocol.
RT-PCR and quantitative real time PCR
Total RNAs were isolated from transfectants using Trizol (Invitrogen) and were reverse transcribed using the SuperScript First-Strand Synthesis System (Invitrogen). PCR amplifications were performed with a series of specific primer sets using the Maxime PCR premix kit (Intron Biotech). Quantitative real time PCR was also performed with the iQ SYBR green supermix (Bio-Rad) using the icycler iQ™ multicolor real-time detection system (Bio-Rad). All qRT-PCR were carried out for 40 cycles under the standard PCR conditions. The results derived from qRT-PCR were analyzed based on threshold cycle (Ct) values. Briefly, a ΔCt was first calculated by subtracting the averaged Ct value of two internal controls (GAPDH and 28S) from the averaged Ct value of each gene, and later a ΔΔCt value was calculated by subtracting the two ΔCt values of the targeted gene derived from the Neuro2A-derived knockdown cell and control cell lines. Fold differences were determined by raising 2 to the ΔΔCt power . The information regarding individual primer sequences and PCR conditions are available upon request.
COBRA (Combined bisulfite restriction analysis)
Genomic DNAs were purified from cell lines using DNAzol (Invitrogen), and 2 μg of each genomic DNA was treated for the bisulfite conversion reaction according to the manufacturer's protocol (EZ DNA methylation kit, Zymo Research). The resultant single-stranded DNAs were used as templates for the PCR reaction using specific primers that were designed for the C-to-T converted DNAs. PCR reaction was performed with the Maxime PCR premix kit (Intron Biotech). Each PCR product was digested with a series of restriction enzymes, first to monitor the efficiency of the bisulfite conversion reaction and later to analyze the methylation status of each region. The restriction enzymes for each region are shown in Fig. 3. The oligonucleotide sequences used for this study are available upon request.
Chromatin immunoprecipitation (ChIP) analysis
Chromatin immunoprecipitations were performed according to the protocol provided by Upstate Biotechnology (Upstate Biotech) with slight modification. Briefly, we used 2 × 107 cells from Neuro2A and two macroH2A1 knockdown stable cell lines [2 and #11]. Formaldehyde was added to the culture medium to a final concentration of 1%, and incubated at 37°C for 10 mins. Treated samples were sheared by sonication and immunoprecipitated with anti-macroH2A1 polyclonal antibody (Cat. # 05-457; Upstate Biotech.). Precipitated chromatins were reverse cross-linked and the DNAs were purified by the phenol/chloroform extraction method. Purified DNAs were used as templates for PCR amplification. PCR reactions were carried out for 35 cycles using standard PCR conditions. The resulting PCR products were analyzed by separating on a 1.5% agarose gel. The oligonucleotide sequences used for this study are available upon request.
We thank Jennifer M. Huang for critically reading the manuscript. This research was supported by National Institutes of Health R01 GM66225 (J.K).
- Wolffe AP: Chromatin: structure and function. 1998, Sandiego, California: Academic Press, 3Google Scholar
- Pehrson JR, Fried VA: MacroH2A, a core histone containing a large nonhistone region. Science. 1992, 257 (5075): 1398-1400. 10.1126/science.1529340.PubMedView ArticleGoogle Scholar
- Pehrson JR, Fuji RN: Evolutionary conservation of histone macroH2A subtypes and domains. Nucleic acids research. 1998, 26 (12): 2837-2842. 10.1093/nar/26.12.2837.PubMed CentralPubMedView ArticleGoogle Scholar
- Chadwick BP, Willard HF: Histone H2A variants and the inactive X chromosome: identification of a second macroH2A variant. Human molecular genetics. 2001, 10 (10): 1101-1113. 10.1093/hmg/10.10.1101.PubMedView ArticleGoogle Scholar
- Doyen CM, An W, Angelov D, Bondarenko V, Mietton F, Studitsky VM, Hamiche A, Roeder RG, Bouvet P, Dimitrov S: Mechanism of polymerase II transcription repression by the histone variant macroH2A. Molecular and cellular biology. 2006, 26 (3): 1156-1164. 10.1128/MCB.26.3.1156-1164.2006.PubMed CentralPubMedView ArticleGoogle Scholar
- Angelov D, Molla A, Perche PY, Hans F, Cote J, Khochbin S, Bouvet P, Dimitrov S: The histone variant macroH2A interferes with transcription factor binding and SWI/SNF nucleosome remodeling. Molecular cell. 2003, 11 (4): 1033-1041. 10.1016/S1097-2765(03)00100-X.PubMedView ArticleGoogle Scholar
- Costanzi C, Pehrson JR: Histone macroH2A1 is concentrated in the inactive X chromosome of female mammals. Nature. 1998, 393 (6685): 599-601. 10.1038/31275.PubMedView ArticleGoogle Scholar
- Chadwick BP, Valley CM, Willard HF: Histone variant macroH2A contains two distinct macrochromatin domains capable of directing macroH2A to the inactive X chromosome. Nucleic acids research. 2001, 29 (13): 2699-2705. 10.1093/nar/29.13.2699.PubMed CentralPubMedView ArticleGoogle Scholar
- Huynh N, Dolan B, Lee S, Whitcher JP, Stanley J: Management of phaeniciatic ophthalmomyiasis externa. The British journal of ophthalmology. 2005, 89 (10): 1377-1378. 10.1136/bjo.2005.071597.PubMed CentralPubMedView ArticleGoogle Scholar
- Costanzi C, Stein P, Worrad DM, Schultz RM, Pehrson JR: Histone macroH2A1 is concentrated in the inactive X chromosome of female preimplantation mouse embryos. Development (Cambridge, England). 2000, 127 (11): 2283-2289.Google Scholar
- Grigoryev SA, Nikitina T, Pehrson JR, Singh PB, Woodcock CL: Dynamic relocation of epigenetic chromatin markers reveals an active role of constitutive heterochromatin in the transition from proliferation to quiescence. Journal of cell science. 2004, 117 (Pt 25): 6153-6162. 10.1242/cs.01537.PubMedView ArticleGoogle Scholar
- Zhang R, Poustovoitov MV, Ye X, Santos HA, Chen W, Daganzo SM, Erzberger JP, Serebriiskii IG, Canutescu AA, Dunbrack RL: Formation of MacroH2A-containing senescence-associated heterochromatin foci and senescence driven by ASF1a and HIRA. Developmental cell. 2005, 8 (1): 19-30. 10.1016/j.devcel.2004.10.019.PubMedView ArticleGoogle Scholar
- Choo JH, Kim JD, Chung JH, Stubbs L, Kim J: Allele-specific deposition of macroH2A1 in imprinting control regions. Human molecular genetics. 2006, 15 (5): 717-724. 10.1093/hmg/ddi485.PubMedView ArticleGoogle Scholar
- Brannan CI, Bartolomei MS: Mechanisms of genomic imprinting. Current opinion in genetics & development. 1999, 9 (2): 164-170. 10.1016/S0959-437X(99)80025-2.View ArticleGoogle Scholar
- Spahn L, Barlow DP: An ICE pattern crystallizes. Nature genetics. 2003, 35 (1): 11-12. 10.1038/ng0903-11.PubMed CentralPubMedView ArticleGoogle Scholar
- Kim J, Stubbs L: Rapidly evolving imprinted loci. Encyclopedia of Genetics, Genomics, Proteomics and Bioinformatics. 2005, West Sussex, UK: John Wiley PublishersGoogle Scholar
- Ventura A, Meissner A, Dillon CP, McManus M, Sharp PA, Van Parijs L, Jaenisch R, Jacks T: Cre-lox-regulated conditional RNA interference from transgenes. Proceedings of the National Academy of Sciences of the United States of America. 2004, 101 (28): 10380-10385. 10.1073/pnas.0403954101.PubMed CentralPubMedView ArticleGoogle Scholar
- Xiong Z, Laird PW: COBRA: a sensitive and quantitative DNA methylation assay. Nucleic acids research. 1997, 25 (12): 2532-2534. 10.1093/nar/25.12.2532.PubMed CentralPubMedView ArticleGoogle Scholar
- Changolkar LN, Costanzi C, Leu NA, Chen D, McLaughlin KJ, Pehrson JR: Developmental changes in histone macroH2A1-mediated gene regulation. Molecular and cellular biology. 2007, 27 (7): 2758-2764. 10.1128/MCB.02334-06.PubMed CentralPubMedView ArticleGoogle Scholar
- Kim J, Kollhoff A, Bergmann A, Stubbs L: Methylation-sensitive binding of transcription factor YY1 to an insulator sequence within the paternally expressed imprinted gene, Peg3. Human molecular genetics. 2003, 12 (3): 233-245. 10.1093/hmg/ddg028.PubMedView ArticleGoogle Scholar
- Kim JD, Hinz AK, Bergmann A, Huang JM, Ovcharenko I, Stubbs L, Kim J: Identification of clustered YY1 binding sites in imprinting control regions. Genome research. 2006, 16 (7): 901-911. 10.1101/gr.5091406.PubMed CentralPubMedView ArticleGoogle Scholar
- Breen GA, Jordan EM: Upstream stimulatory factor 2 stimulates transcription through an initiator element in the mouse cytochrome c oxidase subunit Vb promoter. Biochimica et biophysica acta. 2000, 1517 (1): 119-127.PubMedView ArticleGoogle Scholar
- Kim JD, Hinz AK, Choo JH, Stubbs L, Kim J: YY1 as a controlling factor for the Peg3 and Gnas imprinted domains. Genomics. 2007, 89 (2): 262-269. 10.1016/j.ygeno.2006.09.009.PubMed CentralPubMedView ArticleGoogle Scholar
- Biosystems A: Applied Biosystem Bulletin 2. 1997, Foster City, CAGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.