|
|
|
|
RESEARCH PAPER
1 Department of Medical Oncology, Dana Farber Cancer Institute; Departments of Medicine and Genetics, Harvard Medical School, Boston, Massachusetts 02115, USA; 2 Department of Statistics, Stanford University, Stanford, California 94305, USA
| Abstract |
|---|
|
|
|---|
[Keywords: Histone modifications; knock-out; mSin3 complex; mSin3A; transcriptional regulation; tumorigenesis]
Received December 6, 2004; revised version accepted May 20, 2005.
These biological processes are affected through histone-modifying enzymes whose activities are targeted to specific sites on the DNA via their association with corepression scaffolding molecules that can be recruited by sequence-specific binding proteins such as transcription factors. The mammalian Sin3 proteins, mSin3A and mSin3B, represent prototypical corepressor molecules that were discovered initially through their direct interaction with the sequence-specific transcriptional repressors of the Myc family network, Mad1 and Mxi1 (Ayer et al. 1995
; Schreiber-Agus et al. 1995
; Laherty et al. 1997
) and later shown to associate with HDAC1 and HDAC2 within a large multiprotein complex, the Sin3/HDAC corepressor complex (for review, see Schreiber-Agus and DePinho 1998
; Ayer 1999
; Knoepfer and Eisenman 1999; Silverstein and Ekwall 2004
). Subsequent efforts have also revealed that the mSin3/HDAC proteins are recruited by an expanding array of factors involved primarily in transcriptional regulation including p53 (Murphy et al. 1999
), Ikaros (Koipally et al. 1999
), REST (Huang et al. 1999
), and E2F4 (Rayman et al. 2002
), contributing to repression of their target genes.
mSin3A and mSin3B proteins share extensive similarity within their four paired amphipathic helix (PAH) protein-interaction domains. The presence of two highly related mSin3 homologs in mammals and other species (Dang et al. 1999
; Silverstein and Ekwall 2004
), suggests that each mSin3 homolog serves separable functions. The mSin3 corepressor complex's central enzymatic activity is imparted through its Class I HDAC constituents, HDAC1 and HDAC2. The maintenance of mSin3-associated HDAC activity is dependent on another integral component, mSds3, whose ortholog in Saccharomycese cerevisiae is epistatic to the yeast sin3p and HDAC orthologs (Vannier et al. 1996
; Alland et al. 2002
). Additional core mSin3/HDAC proteins include the retinoblastoma-associated proteins, RbAp46 and RbAp48, which are thought to stabilize interactions with the mSin3 complex and the nucleosome, mSin3-associated proteins (SAPs: SAP18, SAP30, SAP130, and SAP180), RBP1, p33ING1b, and BRMS1, which have been isolated through biochemical means and remain the focus of ongoing functional investigations (for review, see Ayer 1999
; Kuzmichev et al. 2002
; Lai et al. 2001
; Skowyra et al. 2001
; Fleischer et al. 2003
; Meehan et al. 2004
).
As initially reported for the Mad family, the mSin3/HDAC complex also drives the repression activity of numerous other sequence-specific transcription factors whose list has expanded considerably in the past few years. Among those factors that recruit the mSin3/HDAC complex, several are involved in cell survival and apoptosis, thereby linking this complex activity to cancer-relevant processes. For instance, tumor-suppressor proteins, including p53 (Murphy et al. 1999
), pRb (Brehm et al. 1998
; Luo et al. 1998
; Magnaghi-Jaulin et al. 1998
; Lai et al. 2001
), and Menin (Kim et al. 2003
) have been shown to repress their target genes through interaction with the mSin3/HDAC complex. Moreover, aberrant recruitment of this complex by altered transcription factors has been shown to be pathogenetic in several human cancers, such as Acute Promyelocytic Leukemia and Acute Myeloid Leukemia (for review, see Melnick and Licht 2002
).
In addition to its role in transcriptional repression, new activities for the mSin3/HDAC complex are continually uncovered in diverse model organisms. For example, S. cerevisiae sin3p has been shown to participate in the transcriptional activation (as opposed to repression) of MAPK Hog1 target genes upon osmotic stress (De Nadal et al. 2004
), and recent links have been forged between S. cerevisiae sin3 and DNA replication (Aparicio et al. 2004
). S. cerevisiae and Caenorhabditis elegans Sin3 molecules have also been implicated in the DNA damage-repair process in a transcription-independent manner (Pothof et al. 2003
; Jazayeri et al. 2004
). Finally, several Sin3-associated molecules have been shown to participate in activities other than transcriptional regulation including a role for the Sin3 complex in centromere function and genomic stability in both yeast and mammals. Along these lines, genetic inactivation of Pst1 (the homolog of mSin3A in fission yeast) or mSds3 (an essential component of the mSin3/HDAC complex in mouse cells) results in loss of specific chromatin modifications at centromeric or pericentromeric loci, and ultimately, impaired chromosomal segregation (David et al. 2003
; Nakayama et al. 2003
; Silverstein et al. 2003
).
Here, against this backdrop of significant knowledge of the mSin3 corepressor complex, we have sought to better define the biological roles of mSin3A in normal mammalian development, the growth and survival of primary and cancer cells, and the regulation of transcriptional networks. To this end, a conditional mSin3A knock-out model system, biochemical analyses, and computational methods were used in an integrated manner to understand the biological impact of zygotic and somatic deletion of mSin3A in vivo and in cultured primary and transformed cells.
| Results |
|---|
|
|
|---|
Gene targeting and chimera formation were used to generate a conditional mSin3A allele (mSin3AL) in the mouse germline (Fig. 1A,B) that retained wild-type function, as evidenced by normal expression of mSin3A and the lack of a phenotype in mice and derivative cells that were heterozygous or homozygous for mSin3AL (Fig. 1C; data not shown). Moreover, mSin3A+/ mice, generated from E2A-Cre and mSin3AL/+ matings, displayed normal genotype distribution patterns and were outwardly healthy and grossly indistinguishable from mSin3A+/+ littermates up to 1.5 yr of age. While mSin3A+/ intercrosses readily generated mSin3A/ blastocysts, nullizygotes were not present on embryonic day 6.5 (E6.5) and thereafter, indicating that mSin3A serves an essential role in early development circa at peri-implantation (Fig. 1D).
The early lethal phenotype associated with zygotic homozygous deletion of mSin3A prompted somatic conditional mSin3A deletion studies in primary cell cultures. To that end, Cre-encoding retroviruses were used to infect mSin3AL/L mouse embryonic fibroblasts (MEFs), resulting in steady and complete depletion of mSin3A protein only in the mSin3AL/L MEFs over a 5-d period (Fig. 1C). While the well-recognized modest Cre-induced stasis was evident in mSin3A+/+ MEFs, Cre expression in mSin3AL/L MEFs was associated with profound growth arrest and increased apoptosis, as revealed by microscopic inspection, TUNEL assay, and by FACS profiles with a sub-G1 peak (Fig. 2; Supplementary Fig. 1). In addition to a G2/M arrest, the BrdU-PI FACS profiles of mSin3A/ cultures revealed a subpopulation with S-phase DNA content, yet no BrdU incorporation (Fig. 2C, arrow)a profile consistent with unscheduled DNA replication with subsequent triggering of an S-phase checkpoint.
|
Distinctive roles for mSin3A and mSds3 in chromosome biology
Previous yeast studies have suggested a strong epistatic relationship between Sin3 and its associated HDAC and Sds3 components. In mammals, the mSin3A/ and mSds3/ phenotypes appeared to be somewhat overlapping, as inferred by early embryonic lethality and cellular growth arrest and apoptosis. In particular, both mSin3A/ and mSds3/ cells showed G2/M arrest profiles, which raised the possibility that mSin3A may also be involved in chromosome dynamics as reported recently for mSds3 (David et al. 2003
). As reported previously (David et al. 2003
), mSds3/ cells exhibited an increase in multinucleated cells and rampant aneuploidy due to impaired chromosome segregation. In contrast, mSin3A/ MEFs showed no increase in multinucleated cells and possessed relatively normal chromosome counts (Fig. 3A; data not shown). Reminiscent of mSds3/ cells, >75% of the mSin3A/ cells exhibited diffuse HP1
signal throughout the nucleus, rather than the normal concentrated signal within DAPI-stained heterochromatin structures (Fig. 3C).
Despite these abnormal HP1
patterns (and consistent with normal chromosome counts), immunofluorescence studies revealed proper heterochromatin localization of methylated Lys 9 of histone H3 (meH3K9) and Lys 20 at histone H4 (meH4K20), as well as hypoacetylated Lys 5 and Lys 12 of histone H4 (AcH4K5 and AcH4K12) (Fig. 3B,C; data not shown). Additionally, no differences were observed in overall acetylation status of Lys 5, Lys 8, Lys 12, and Lys 16 of histone H4, indicating that mSin3A is not involved in global deacetylation of histone H4 (data not shown). Notably, this mSin3A/ heterochromatin profile is highly reminiscent of Orc2 (origin of replication 2)-depleted cells that similarly show a DNA replication block, G2/M arrest, and loss of HP1
pericentric heterochromatin localization, yet lack of chromosomal instability (Prasanth et al. 2004
). Accordingly, Orc2 failed to localize to pericentric heterochromatin in mSin3A/ cells (Fig. 3D), pointing to a link between Orc2 and mSin3A in the proper targeting of HP1
to pericentric heterochromatin. However, as Orc2 disappears from pericentric heterochromatic during late G2/M, one cannot exclude the possibility that Orc2 is not present at pericentric heterochromatin as mSin3A/ cells are constrained at the late G2/M phase. In summary, mSin3A and mSds3 deficiencies differentially impact on cell cycle kinetics, as evidenced by the specific link of mSin3A to DNA replication processes and the seminal role of Sds3 in chromosome segregation. These differences in chromosome dynamics suggest that either mSds3 functions independently of mSin3A, or that mSin3A and B homologs serve overlapping functions in these specific cellular processes.
|
The cell cycle arrest and apoptosis phenotypes of mSin3A deficiency prompted an in-depth genetic analysis of the p53p21Cip1 pathway, particularly in the light of previous reports revealing (1) mSin3Ap53 physical association (Murphy et al. 1999
), (2) modulation of p53 acetylation and transcriptional activity upon genotoxic stress (Sakaguchi et al. 1998
), and (3) capacity of HDAC inhibitors or genetic deletion of Hdac1 to induce the p53 transcriptional target p21Cip1 (Lagger et al. 2002
, 2003
; for review, see Johnstone and Licht 2003
). In line with these observations, we observed consistent induction of p21Cip1 expression upon Cre-mediated deletion of mSin3A (Fig. 4A,B). Indeed, chromatin immunoprecipitation (ChIP) experiments revealed the presence of mSin3A on the promoter of p21Cip1 (data not shown). While mSin3AL/L MEFs infected with Cre retroviruses showed no significant additional increase in p21Cip1 protein levels, independent experiments revealed hyperacetylation of p53 on Lys 317 and Lys 379, upon treatment with doxorubicin (Fig. 4B). These data strongly suggest that mSin3A regulates the p53p21Cip1 axis on two levels; it represses p21Cip1 transcription by regulating acetylation of histones in the promoter of p21Cip1 and regulates p53 activity via deacetylation upon genotoxic stress.
To substantiate the mSin3Ap53p21Cip1 interrelationship on the genetic level, mSin3A+/+p53/, mSin3AL/+p53/ and mSin3AL/p53/ MEFs were infected with either Cre or control retrovirus, and subsequently analyzed for DNA damage-induced p21Cip1 induction and associated effects on proliferation and cell cycle. p21Cip1 induction was found to be robust in mSin3A/p53+/+ cells and totally abolished in mSin3A/p53/ cells (Fig. 4A), indicating that p53 is responsible for genotoxin-induced p21Cip1 up-regulation in the setting of mSin3A deficiency. Despite this, p53 deficiency did not alleviate the mSin3A/ phenotypes of poor cell proliferation, reduced S-phase and increased G2/M fractions, block in DNA replication, and marked apoptosis (Fig. 4CE; Supplementary Fig. 3A). Additionally, functional inactivation of both the p53 and pRb pathways by means of expression of SV40-Large T Antigen (T-Ag) did not overcome the cell lethal impact of mSin3A deficiency (Supplementary Fig. 3B).
|
The mSin3A transcriptome
Since mSin3 proteins serve integral roles in the regulation of gene expression (Knoepfer and Eisenman 1999) and mSin3A loss exerts broad biological effects (this study), transcriptional profiling and in silico pathway analysis were utilized to understand more precisely these diverse actions of mSin3A. In time course studies, the 4-hydroxytamoxifen (4-OHT)-inducible Cre/Estrogen Receptor (CreERT2) system was used, as it affords tight temporal control of Cre activity in mSin3AL/L MEFs. As shown in Figure 6A, 4-OHT induction of the CreERT2-transduced mSin3AL/L resulted in a cellular phenotype that was similar to that brought about by constitutive Cre retrovirus in mSin3AL/L MEFs (cf. Fig. 2A), although cell cycle arrest and apoptosis were less pronounced, as observed with cells expressing constitutive Cre-recombinase (Supplementary Fig. 4). On the basis of the kinetics showing steady depletion of mSin3A mRNA and protein levels (Fig. 6B), a comparative time-course transcriptional profile of CreERT2-transduced mSin3A+/+ and mSin3AL/L MEFs was performed at times 0, 24, 48, 72, and 96 h following 4-OHT treatment. This approach allows the monitoring of differential gene expression upon mSin3A deletion before the onset of the mSin3A phenotype, and more readily enables discrimination between direct transcriptional effects due to mSin3A ablation versus nonspecific secondary transcriptional changes accompanied by the aggravated cell cycle arrest and apoptosis. dChip software (Li and Wong 2001
) was used to analyze the microarray data generated with Affymetrix M430A oligo DNA microarrays.
An Analysis of Variance (ANOVA) procedure was performed to identify genes with nonconstant expression values between the two genetic backgrounds (mSin3A+/+ vs. mSin3A/). Imposing a p-value threshold of <0.0005, 1308 genes were determined as having significant nonconstant expression. In concordance with the depletion kinetics of the mSin3A protein, most of these genes showed differential expression starting at 48 h after addition of 4-OHT. It is worth noting that there were no obvious cellular phenotypes at this 48-h time point, and thus it is more likely that the observed transcriptional changes are more proximal to mSin3A function as opposed to secondary events and terminal cellular collapse (Supplementary Fig. 4). The majority (977) of these differentially expressed genes were up-regulated (Fig. 6C; Supplementary Table 1; Supplementary Information 1, 2). A subset of the genes showing altered expression on the Affymetrix expression profiles were verified by real-time RTPCR analysis on independent generated RNA samples from mSin3AL/L MEFs expressing Cre-ERT2 (Fig. 7A). As shown in Figure 7A, whereas mSin3A RNA levels declined progressively at 48 and 72 h of 4-OHT treatment, increased mRNA levels for mSin3B, Mxi1, Sap30, Ing1, Arid4b, Peg3, Set, p130, cyclin E, Odc1, Msh6, Brca1, 53Bp1, Stmn1, and Aldh6a1 were already markedly increased at 48 h, findings consistent with a direct link to mSin3A regulation. Correspondingly, ChIP analysis of mSin3B, Mxi1, Ing1, Arid4b, Peg3, Set, p130, Caspase 3, 53Bp, Pmm1, and Aldh6a1 showed that mSin3A physically localizes to the proximal promoters of these genes (Fig. 7B).
|
Several up-regulated genes functionally cluster as genes involved in G2/M transition and cytokinesis including 53Bp1, p130, Gadd45-
, Mcm4, Mcm5, Mcm7, Securin, Set, Tacc3, and Cdc2A. None of these genes have been reported to be specifically up-regulated in G2 or mitosis, suggesting that their up-regulation is not merely a consequence of the G2/M cell cycle arrest imposed by mSin3A deficiency (Whitfield et al. 2002
). Therefore, the fact that some of these genes, when overexpressed, give rise to a G2/M arrest (e.g., Gadd45
and 53Bp1) provides a possible explanation for the observed G2/M cell cycle arrest in mSin3A/ cells. Similarly, the apoptosis associated with mSin3A depletion could relate in part to transcriptional up-regulation of Caspase-3, Ing1, Peg3, Pdcd5, Trim27, Mcl1, Pin1, and sphingosine phosphate lyase 1, all encoding pro-apoptotic activities (Nicholson et al. 1995
; Helbing et al. 1997
; Reiss et al. 2004
).
|
A prominent group of genes present in the mSin3A transcriptome participate in DNA-repair processes such as Rpa1, Msh6, Rad17, Lig1, Lig3, Apex1, and Nsep1. These observations are in line with work in yeast supporting a role for Sin3 in the nonhomologous end-joining pathway (Jazayeri et al. 2004
). In this regard, there is an up-regulation of several DNA repair genes involved in the homologous recombination pathway including Brca1, 53Bp1, and Rad51. The up-regulation of homologous recombination genes in a context of impaired NHEJ activity raises the intriguing possibility that mSin3A may be involved in coordinating the relative levels of NHEJ versus homologous recombination repair activities in the cell.
Consistent with findings in Drosophila, mSin3A deletion resulted in the up-regulation of genes involved in cytosolic and mitochondrial energy-generating pathways (Pdk1, Mrsp23, Oghd, Adac9, Bdh, Got2, Hadh2, Cs, Fxc1, Prdx3), strengthening the provocative role for mSin3A in the regulation of mitochondrial respiratory activity (Pile et al. 2003
; see below). Finally, mSin3A regulates the expression of several genes involved in DNA replication (Rpa1, Lig1, Mcm4, Mcm5, Mcm7, Nap1l1) providing a basis for the observed defects in DNA replication upon deletion of mSin3A, and strengthening the case for a direct role for mSin3A in the DNA replication process.
Identification of known and novel mSin3A transcriptional nodes
Enrichment of known Myc, E2F, and p53 targets in the mSin3A transcriptome prompted us to search the regulatory regions of all genes regulated by mSin3A for additional Myc, E2F, and p53 target genes, as well as for unanticipated links other transcription factors and their binding sites. In the case of E2F, 66 mSin3A/ up-regulated genes contained E2F-binding sites within 500 bp upstream or downstream of the transcriptional start site. Compared with a random collection of 1161 genes, there was a 1.95-fold enrichment (4.40% vs. 8.61%; p = 0.044) in E2F-binding sites in the collection of induced genes upon deletion of mSin3A. Scanning for p53-binding elements (which are typically positioned at greater distances from the transcriptional start site), 64 mSin3A/ up-regulated genes were found to contain p53-binding sites within 5000 bp of the transcriptional start site in either upstream or downstream direction. A 1.85-fold enrichment was observed for p53-binding sites in the mSin3A/ up-regulated genes compared with the random genes (8.47% vs. 4.57%, p = 0.046). Finally, a search for Myc target genes by scanning the promoters of the up-regulated genes for E-box sequences identified 189 genes and 1.34-fold enrichment (24.6% vs. 18.33%, p = 0.183) for E-box sequences over a set of random genes. The relatively low enrichment for E-box sequences is likely due to the reduced information present in E-box sequences compared with the longer E2F- and p53-binding sequences, as well as the known ability of the Myc oncoprotein to regulate a significant fraction of the genome, as revealed by studies in human and Drosophila (Fernandez et al. 2003
; Orian et al. 2003
). Together, these findings show that gene expression profiling and in silico promoter analysis can lead to the identification of E2F-, Myc-, and p53-responsive targets.
|
|
Next, we conducted an unbiased search for cis-elements that are overrepresented in the promoters of validated mSin3A target genes: mSin3B, Mxi1, Ing1, Arid4b, Sap30, 53Bp1, Peg3, Set, Caspase 3, p130, and Aldh6a1 (i.e., genes confirmed by ChIP analysis and real-time RTPCR) (Fig. 7A,B). DNA sequences spanning 500 bp on each side of the mSin3A ChIP-positive genome fragment were overlaid across genomes of multiple species to identify conserved regions (Siepel and Hauser 2005). These conserved regions were subsequently analyzed for DNA-binding motifs using the MEME motif-finding program (Bailey and Elkan 1994
) and compared with known transcription factor-binding motifs using the Match Program in Transfac (Matys et al. 2003
). As shown in Figure 7C, two cis-regulatory elements were identified that are similar to DNA-binding motifs for the transcriptional activator Stat, whose DNA-binding activity has been shown to be acetylation dependent (Wang et al. 2005
; Yuan et al. 2005
) and Falz (also know as FAC1/BPTF/Nurf301), a nucleosome remodeling factor that is used in the ISWI ATP-dependent chromatin remodeling complex NURF (Bowser et al. 1995
; Jordan-Sciutto et al. 1999
; Jones et al. 2000
; Xiao et al. 2001
; Badenhorst et al. 2002
). Together, these data suggest that the combined computational and ChIP analyses provides a means for the validation of known mSin3A cofactors and an entry point for the identification of novel mSin3A cofactors.
| Discussion |
|---|
|
|
|---|
Despite up-regulation of mSin3A partners, as well as ubiquitous coexpression of the highly related mSin3B homolog, the loss of mSin3A function exacted profound physiological consequences, underscoring the essential and separable roles of the mSin3 homologs. That mSin3A and mSin3B possess nonoverlapping functions is in line with the capacity of mSin3A and mSin3B to occupy different promoters in mammalian cells (Rayman et al. 2002
), as well as the assignment of distinct functions for the two yeast Sin3 orthologs (Nakayama et al. 2003
; Silverstein et al. 2003
). The nature and extent of the common and distinct functions of mSin3A and mSin3B will require additional direct genetic study. Along these lines, it will be paramount to determine whether there is either functional redundancy among mSin3A and mSin3B in the establishment of functional centromeric heterochromatin, or whether there is indeed a unique mSin3A/B-independent function for mSds3 in this aspect of chromosome biology.
At the same time, although mSds3 appears to be critical for maintenance of HDAC-associated enzymatic activity in the mSin3A and mSin3B corepressor complex (Alland et al. 2002
), it is important to note that mSin3A/ cells exhibit phenotypes that are not observed in the mSds3/ cells. Specifically, mSin3A/ MEF cultures harbor a subpopulation containing S-phase DNA content with no ongoing DNA replication (no BrdU incorporation). Given that this DNA content anomaly cannot be attributed to chromosomal instability and that there is prominent representation of DNA replication genes in the mSin3A transcriptome, we propose that mSin3A plays a direct role in regulation of DNA replication in mammalian cells. This hypothesis is in line with previous studies in S. cerevisiae, showing a role for Sin3 and Rpd3 (ortholog of HDAC1) in the activation of late origins of replication independent relating to their role in transcriptional repression (Aparicio et al. 2004
). Additionally, in Drosophila, activation of origin of replication is dependent on Rpd3 as evidenced by the loss of Rpd3 function, resulting in genome-wide hyperacetylation, genomic replication, and redistribution of the origin binding protein ORC2 (Aggarwal and Calvi 2004
). Together, these collective insights have shaped our view of the mSin3 corepressor complex. Rather than the existence of a monolithic collection of components, the data are consistent with the ability of the mSin3 complex to operate in a more modular fashion, such that various transcription factors and other DNA-binding elements may recruit only a subset of the mSin3/HDAC components. This speculation is consistent with emerging observation of the existence of several biochemically distinct forms of the mSin3 multiprotein complex (Fleischer et al. 2003
).
In the context of cancer, a strong consensus has emerged that modulators of chromatin architecture are intimately linked to key cancer-relevant pathways. Particularly relevant to this issue is the previously documented physical and biochemical interactions between the mSin3A corepressor and the p53 tumor suppressor. Here, our genetic analysis has substantiated the ability of mSin3A to modulate p53 activation in the setting of DNA damage, as evidenced by Lys 317 and Lys 379 hyperacetylation and increased p53 protein levels upon loss of mSin3A function. Thus, we propose that an important function of the mSin3A/HDAC complex is the regulation of p53 acetylation status and cell cycle control following genotoxin-induced injury. In addition, it is worth noting that several deacetylases and HDAC-containing complexes have shown the capacity to deacetylate p53, including the NuRD complex (Luo et al. 2000
), all Class I HDACs (Juan et al. 2000
), and the class III deacetylase SIRT1 (Luo et al. 2001
; Vaziri et al. 2001
; Cheng et al. 2003
). These observations, together with the findings of this report, raise the interesting possibility that multiple protein complexes, and therefore, multiple biochemical signaling pathways, converge on the p53 protein to regulate its activity via acetylation in various physiological settings.
In silico promoter scanning revealed significant enrichment of Myc, E2F, and p53 cis-elements in the mSin3A transcriptome. It is interesting to note, however, that not all Myc, E2F, and p53 target genes are up-regulated, which implies that mSin3A might regulate only a subset of these targets, that mSin3A may only engage these transcription factors under certain physiological conditions or cell types, and/or that mSin3B may serve overlapping functions on a subset of these gene targets. One important caveat to these suppositions is the fact that only a limited region of the promoter sequences of the mSin3A transcriptome genes was scanned, and it is possible that Myc, E2F, and p53 cis-regulatory elements were not captured in our analyses, thereby resulting in an underestimation of the mSin3A target genes coregulated by Myc, p53, and E2F. The identification of the Stat and Falz cis-regulatory elements in a subset of the mSin3A target genes shows that mSin3A uses various DNA-binding transcription factors to confer selectivity and specificity in transcriptional regulation. Recently, it was shown that acetylation of Stat3 is critical for the protein to form stable dimers required for cytokine-stimulated DNA-binding and transcriptional regulation and to enhance transcription of cell growth-related genes (Wang et al. 2005
; Yuan et al. 2005
). These studies, together with our finding that mSin3A can repress Stat target gene transcription, suggests a model in which mSin3A-associated HDAC activity inhibits Stat3 dimerization and neutralizes its transcriptional potential. Various studies have implied collaborative interactions between different chromatin remodeling complexes, including mSin3A/HDAC and SWI/SNF (Sternberg et al. 1987
; Sif et al. 2001
; Kuzmichev et al. 2002
). In this context, our finding that mSin3A regulates gene transcription through Falz appears relevant, as Falz is a component of the ISWI ATP-dependent nucleosome remodeling complex NURF. Falz contains a bromodomain, which are known to bind acetylated lysines in histones and nonhistone proteins (for review, see Yang 2004
), suggesting that Falz might recruit mSin3A/HDAC complexes to such targets to deacetylate specific promoter regions or modify the transcriptional activity of nonhistone proteins. A human NURF complex was shown to contain, beside FALZ, also RbAp46 and RbAp48, two proteins that are often found in complex with mSin3A (Barak et al. 2003
). Notably, Nurf301 can act as a negative regulator of the JAK/STAT signaling pathway, downstream of the STAT proteins, suggesting that mSin3A may regulate gene transcription through both DNA-binding transcription factors (Badenhorst et al. 2002
). Along these lines, our computational and ChIP analyses reveal that the mSin3A target gene Sap30 contains Statas well as Falz-binding motifs in its promoter.
In summary, the integrated genetic, biochemical, and computational approaches of this study has enabled us to gain insight into the diverse functions of mSin3A. The application of this integrated approach and conditional somatic inactivation may find utility in the analysis of multifunctional complexes serving cell essential functions. In the case of mSin3A, this study confirmed the essential role of mSin3A, and by extension, the mSin3A corepressor complex in numerous biological processes relevant to normal physiology and to cancer, and served to uncover known and novel cofactors regulating various biochemical processes (e.g., DNA replication, nuclesomal remodeling, etc.).
The observations of this report suggest a general strategy for assessment of the other components of the mSin3 corepressor complex in normal cells and in cancer.
| Materials and methods |
|---|
|
|
|---|
A conditional mSin3A-targeting construct was generated by flanking mSin3A exon 4 with a 5' loxP-Frt-PGK-neomcyin expression cassette (Aguirre et al. 2003
) and a 3' loxP site. Upon electroporation of 129/SvOla embryonic stem (ES) cells and selection at 200 µg/mL G418, homologous recombinants were identified by Southern blot analysis. Correctly targeted ES clones were injected into C57Bl/6 blastocysts and germline transmitting chimeric mice were bred to CAGG-Flpe (Rodriguez et al. 2000
) and EIIa-Cre transgenic mice (Lakso et al. 1996
) to generate the mSin3AL and mSin3A alleles, respectively. Mice were genotyped by Southern analysis and multiplex PCR (primers and conditions are available upon request).
Preparation of MEFs, cell culture, and retroviral infection
MEFs were isolated from 13.5-d-post-coitum (dpc) embryos using standard procedures (Dannenberg et al. 2000
) and cultured in DMEM supplemented with 10% fetal bovine serum (FBS), 0.1 mM
-mercapto-ethanol, and antibiotics. Ecotropic retroviral supernatants were produced by cotransfection of 293T packaging cells with 6 µg viral DNA and 6 µg pCL vector (Naviaux et al. 1996
) using lipofectamin 2000 (Invitrogen) according to the manufacturer's instructions. Viral supernatant was collected and filtered through a 0.45-µm filter 48, 56, and 72 h post-transfection. Low-density cultures of early passage MEFs were infected two times for at least 6 h with viral supernatants containing 4 µg/mL polybrene (Sigma) with an 8-h interval. After a recovery period of 36 h, the MEFs were selected for at least 48 h using medium containing 2.5 µg/mL puromycin (Sigma) before the cells were subjected to any of the analyses reported in this study.
To study p53 acetylation, MEFs were infected with either control or Cre-expressing viruses and treated with doxorubicin 72 h post-infection and processed as described in Cheng et al. (2003
).
Proliferation curves, cell cycle, and apoptosis analyses
Proliferation curves were carried out as described in Dannenberg et al. (2000
). Cell cycle analysis was performed by incubating MEFs on 10-cm tissue culture plates with 10 µM 5-Bromodeoxy uridine (Sigma) for 1 h and processed as described in Dannenberg et al. (2000
). Assessment of apoptosis was performed by growing MEFs on eight-well chamber slides for 2448 h, and subsequently analyzed using an In Situ Cell Detection Kit (Roche) according to the manufacturer's instructions. Quantification of apoptotic cells was performed by counting the number of TUNEL-positive nuclei in 100 DAPI-stained nuclei in three different fields at 100x magnification. All experiments were performed at 57 d after retroviral infection with pBabepuro or pBabe-Cre virus.
Metaphase spreads and immunofluoresence
Metaphase spreads from MEFs were generated as described in David et al. (2003
). Immunofluoresence analysis of MEFs was performed 57 d after retroviral infection with pBabe-puro or pBabe-Cre virus. MEFs were cultured for 2448 h on eight-well chamber slides, fixed with 2% fresh paraformaldehyde (Pierce) and permeabilized with 0.1% Triton X-100. Antibodies against the following proteins were used: AcK12H4, AcK5H4, meK9H3, meK20H4 (Upstate), HP1
(Chemicon). For the Orc2 (BD Biosciences)-staining, Triton extraction was performed prior to fixation as described in Prasanth et al. (2004
). Immunofluorescent images were acquired using either a Zeiss Axiovert 200-M inverted microscope or a CoolSnap HQ CCD camera (Photometrics) driven by the Slidebook 4.0 b 2.5 software (intelligent Imaging Innovations). Pictures were acquired with 1 x 1 binning and a 0.1-mm Z-step. Deconvolution was performed by the nearest neighbor method with 20 iterations.
Protein analyses
Protein levels were determined by Western blot analyses using routine protocols. Antibodies against mSin3A (K-20), mSin3B (AK-12), p130 (C-20), Cdk4 (C-22), and p21Cip1 (C-19) were obtained from Santa Cruz, HDAC1 from Upstate, p53 antibody (Ab-7) from Oncogene Science, and p53AcK317 and p53AcK379 from Trevigen. Peroxidase-conjugated goat anti-rabbit and goat anti-mouse secondary antibodies were from Pierce. To detect anti-p53 antibodies, we used peroxidase-conjugated bovine anti-sheep antibodies from Santa Cruz.
Generation and culturing of mSin3AL/Lp53/ sarcomas and lymphomas
mSin3A+/+p53/ and mSin3AL/Lp53/ mice were generated by mSin3AL/+p53/ intercrosses, and subsequently allowed to age. On average, at 1015 wk of age, these mice developed (angio)sarcoma and T-cell lymphomas. Single-cell suspensions of sarcomas and lymphomas were prepared and taken into culture. Sarcomas were maintained in culture using exactly the same conditions as for MEFs. Lymphomas were cultured in culture flasks in RPMI medium (Gibco), supplemented with 10% FBS, 0.1 mM
-mercapto ethanol, and antibiotics.
Gene expression profiling
Subconfluent mSin3A+/+ and mSin3AL/L MEFs cultures on 10-cm culture dishes were twice infected with retroviruses expressing the CreERT2 recombinase and selected for puromycin (2.5 µg/mL) for the rest of the experiment. After reaching confluency, the cells were passaged onto 15-cm culture dishes and selected for an additional 2 d. Subsequently, each plate was split into 5 x 3 10-cm culture dishes. For five consecutive days, 4-hydroxy tamoxifen was added to a new set of three 10-cm dishes at a final concentration of 100 mM. RNA and protein were isolated at the end of the experiment to document mSin3A status. Two 10-cm plates for each time point were used to isolate RNA using Trizol (Invitrogen) and an RNeasy mini kit (Qiagen). Gene expression profiling was performed on the Affymetrix M430A platform.
Real time RTPCR
Real-time PCR was performed using the Quantitec SYBR kit from Qiagen according to the manufacturer's prescription. Primer sequences are available on request. RNA was isolated from mSin3AL/L MEFs expressing CreERT2 and treated with Tamoxifen for either 0, 48, or 72 h using Rneasy kits (Qiagen). To calculate RNA quantities, we established standard curves with serial dilutions of known RNA concentrations. One-hundred nanograms of total RNA was used for each RTPCR (Qiagen). All RTPCRs were performed in triplicate.
ChIP and identification of cis-regulatory elements
ChIP experiments were performed using confluent early passage (passage 4) primary MEFs (pMEFs). For each ChIP assay, we used 3 x 107 cells and followed a published protocol (Ren et al. 2002
). Antibodies used for ChIP are anti-mSin3A (K-20; Santa Cruz), anti-SV40 LargeT (gift of J. DeCaprio, Dana-Farber Cancer Institute, Boston, MA), anti-acetylated histone H4 (Upstate). PCR primers were designed for each kilobase of a 5-kb promoter fragment, and are available on request. Information on the methodology used in the identification of novel cis-regulatory elements in promoters of mSin3A target genes can be found in the Supplemental Material.
Bioinformatic analyses of the gene expression profiling data
Analysis of the raw gene expression profiling was performed using dChip software (Li and Wong 2001
). In order to identify potential E2f-binding, sites in the promoters of 1182 probe sets that showed significant (p = 0.0005) up-regulation upon mSin3A deletion, we mapped the genes onto the genome and determined the transcriptional start site (TSS). Promoters as large as 500 bp upstream and 500 bp down stream of the TSS were scanned with a known E2F matrix for potential E2f-binding sites. Only the two best hits per promoter were recorded for location and score. p53-binding sites in the same promoters were found by following the same procedure, except that the scan region for the p53 matrix was expanded to 5000 bp upstream and 5000 bp down stream of the TSS. E-box sequences were identified in a region spanning 1000 bp upstream and 200 bp downstream of the TSS. Sequences for the E2F, p53, Myc, and AP-1-binding sites used to scan the promoter regions of the genes in the mSin3A transcriptome are documented in Transfac (Matys et al. 2003
). To identify gene networks within the mSin3A transcriptome, we used the Ingenuity Pathways Knowledge Base (http://www.ingenuity.com).
| Acknowledgments |
|---|
|
|
|---|
| Footnotes |
|---|
Article and publication are at http://www.genesdev.org/cgi/doi/10.1101/gad.1286905.
A study by Cowley et al. (2005
) using an independently generated mSin3A conditional mouse shows similar effects of mSin3A deficiency on embryonic development, MEF proliferation, and chromosomal dynamics.
4 These authors equally contributed to this work. ![]()
3 Present address: Department of Pharmacology, New York University, School of Medicine, New York, NY, 10016 USA. ![]()
E-MAIL ron_depinho{at}dfci.harvard.edu; FAX (617) 632-6069. ![]()
| References |
|---|
|
|
|---|
Aguirre, A.J., Bardeesy, N., Sinha, M., Lopez, L., Tuveson, D.A., Horner, J., Redston, M.S., and DePinho, R.A. 2003. Activated Kras and Ink4a/Arf deficiency cooperate to produce metastatic pancreatic ductal adenocarcinoma. Genes & Dev. 17: 31123126.
Alland, L., David, G., Shen-Li, H., Potes, J., Muhle, R., Lee, H.C., Hou, H.J., Chen, K., and DePinho, R.A. 2002. Identification of mammalian Sds3 as an integral component of the Sin3/histone deacetylase corepressor complex. Mol. Cell. Biol. 22: 27432750.
Aparicio, J.G., Viggiani, C.J., Gibson, D.G., and Aparicio, O.M. 2004. The Rpd3Sin3 histone deacetylase regulates replication timing and enables intra-S origin control in Saccharomyces cerevisiae. Mol. Cell. Biol. 24: 47694780.
Ayer, D.E. 1999. Histone deacetylases: Transcriptional repression with siners and nurds. Trends Cell. Biol. 9: 193198.[CrossRef][Medline]
Ayer, D.E., Lawrence, Q.A., and Eisenman, R.N. 1995. Mad-Max transcriptional repression is mediated by ternary complex formation with mammalian homologs of yeast repressor Sin3 Cell 80: 767776.[CrossRef][Medline]
Badenhorst, P., Voas, M., Rebay, I., and Wu, C. 2002. Biological functions of the ISWI chromatin remodeling complex NURF. Genes & Dev. 16: 31863198.
Bailey, T.L. and Elkan, C. 1994. Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc. Int. Conf. Intell. Syst. Mol. Biol. 2: 2836.[Medline]
Barak, O., Lazzaro, M.A., Lane, W.S., Speicher, D.W., Picketts, D.J., and Shiekhattar, R. 2003. Isolation of human NURF: A regulator of Engrailed gene expression. EMBO J. 22: 60896100.[CrossRef][Medline]
Bowser, R., Giambrone, A., and Davies, P. 1995. FAC1, a novel gene identified with the monoclonal antibody Alz50, is developmentally regulated in human brain. Dev. Neurosci. 17: 2037.[Medline]
Brehm, A., Miska, E.A., McCance, D.J., Reid, J.L., Bannister, A.J., and Kouzarides, T. 1998. Retinoblastoma protein recruits histone deacetylase to repress transcription. Nature 391: 597601.[CrossRef][Medline]
Cheng, H.L., Mostoslavsky, R., Saito, S., Manis, J.P., Gu, Y., Patel, P., Bronson, R., Appella, E., and Alt, F.W. 2003. Developmental defects and p53 hyperacetylation in Sir2 homolog (SIRT1)-deficient mice. Proc. Natl. Acad. Sci. 100: 1079410799.
Cowley, S.M., Iritani, B.M., Mendrysa, S.M., Xu, T., Cheng, P.F., Yada, J., Liggitt, H.D., and Eisenman, R.N. 2005. The mSin3A chromatin modifying complex is essential for embryogenesis and T cell development. Mol. Cell. Biol. (in press).
Dang, V.D., Benedik, M.J., Ekwall, K., Choi, J., Allshire, R.C, and Levin, H.L. 1999. A new member of the Sin3 family of corepressors is essential for cell viability and required for retroelement propagation in fission yeast. Mol. Cell. Biol. 19: 23512365.
Dannenberg, J.-H., van Rossum, A., Schuijff, L., and te Riele, H. 2000. Ablation of the retinoblastoma gene family deregulates G(1) control causing immortalization and increased cell turnover under growth-restricting conditions. Genes & Dev. 14: 30513064.
David, G., Turner, G.M., Yao, Y., Protopopov, A., and DePinho, R.A. 2003. mSin3-associated protein, mSds3, is essential for pericentric heterochromatin formation and chromosome segregation in mammalian cells. Genes & Dev. 17: 23962405.
De Nadal, E., Zapater, M., Alepuz, P.M., Sumoy, L., Mas, G., and Posas, F. 2004. The MAPK Hog1 recruits Rpd3 histone deacetylase to activate osmoresponsive genes. Nature 427: 370374.[CrossRef][Medline]
Fan, M., Rhee, J., St-Pierre, J., Handschin, C., Puigserver, P., Lin, J., Jaeger, S., Erdjument-Bromage, H., Tempst, P., and Spiegelman, B.M. 2004. Suppression of mitochondrial respiration through recruitment of p160 myb binding protein to PGC-1
: Modulation by p38 MAPK. Genes & Dev. 18: 278289.
Fernandez, P.C., Frank, S.R., Wang, L., Schroeder, M., Liu, S., Greene, J., Cocito, A., and Amati, B. 2003. Genomic targets of the human c-Myc protein. Genes & Dev. 17: 11151129.
Fischle, W., Wang, Y., and Allis, C.D. 2003. Histone and chromatin cross-talk. Curr. Opin. Cell. Biol. 15: 172183.