Previous Article | Next Article ![]()
Molecular and Cellular Biology, November 2007, p. 7871-7885, Vol. 27, No. 22
0270-7306/07/$08.00+0 doi:10.1128/MCB.00990-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Department of Biology, University of Texas at San Antonio, San Antonio, Texas 78249,1 Department of Cellular and Structural Biology, University of Texas Health Science Center at San Antonio, San Antonio, Texas 782292
Received 5 June 2007/ Returned for modification 23 July 2007/ Accepted 8 September 2007
|
|
|---|
|
|
|---|
Among the many germ cell-specific genes expressed during spermatogenesis is a set encoding isozymes that function during glycolysis, primarily in spermatozoa. These include the Ldhc, Gapdhs, Hk1s, Pgm2, Pdha2, and Pgk2 genes, which encode sperm-specific isozymes of lactate dehydrogenase, glyceraldehyde-3-phosphate dehydrogenase, hexokinase, phosphoglycerate mutase, pyruvate dehydrogenase E1 alpha, and phosphoglycerate kinase, respectively (13, 17, 20, 42). In mammals, two genes encode the glycolytic enzyme phosphoglycerate kinase, the X-linked Pgk1 gene and the autosomal Pgk2 gene (58). The Pgk1 gene is widely expressed in all somatic cells, oogenic cells, and premeiotic spermatogenic cells, whereas the Pgk2 gene is expressed exclusively in meiotic spermatocytes and postmeiotic spermatids, where expression of the Pgk1 gene is repressed by meiotic sex chromosome inactivation (MSCI) (57) and by the subsequent repressive effects of postmeiotic sex chromatin (PMSC) (38, 45, 60).
The promoter of the Pgk2 gene is believed to have derived from that of the Pgk1 gene as part of a retroposition event that led to the genesis of the Pgk2 gene early during mammalian evolution (12, 40). Thus, although the Pgk2 promoter was likely initially identical to the Pgk1 promoter, it appears to have diverged, especially in eutherian mammals, to facilitate tissue-specific expression. Evidence favoring a common ancestry of these two promoters includes the absence of a TATA box and the presence of GC and CAAT box elements in both (34, 35). Although this "core promoter" region from the Pgk2 gene is able to direct transcription of a downstream reporter gene in transient transfection assays (12), this portion alone is insufficient to direct expression of a reporter gene in vivo in transgenic mice (48). Indeed, a minimum of an additional 42 bp of enhancer sequence immediately upstream of the Pgk2 core promoter, along with the core promoter itself, is required to direct testis-specific expression in vivo (62). The acquisition of tissue-specific function at this enhancer region, along with the loss in the Pgk2 promoter of the CpG island present in the Pgk1 promoter, represents key differences that we previously suggested may have facilitated the evolution of testis-specific regulation of the Pgk2 gene (34, 35). We previously used a combination of analyses in vitro and in vivo to determine that spermatogenesis-specific expression of the Pgk2 gene is regulated by a 5'-flanking region that includes enhancer and core promoter elements that bind testis-specific and ubiquitous transcription factors, respectively (21, 48, 62). We also showed that demethylation of a domain encompassing this regulatory region precedes decondensation of chromatin, factor binding, and initiation of transcription by several days (3, 4, 22, 28, 29, 63).
The mammalian spermatogenic cell lineage is well characterized with respect to the sequential appearance of specific spermatogenic cell types, permitting a detailed study of developmental processes (37, 39). Here we have taken advantage of this system to examine in vivo the developmental dynamics of protein-DNA interactions that regulate the activation of the Pgk2 gene and to identify the key factors involved. By characterizing expression patterns during spermatogenesis of the genes encoding these factors, we can suggest those most likely to be the key regulators of tissue-specific, cell type-specific, and developmental stage-specific transcription of the Pgk2 gene during spermatogenesis. Our results provide the first demonstration of a spermatogenesis-specific gene that is regulated by the testis-specific homeodomain factor PBX4 and form the basis for a model of gene activation whereby this "pioneer transcription factor" (8, 15, 16, 50) initiates an ordered cascade of events that includes remodeling of chromatin and binding of other transcription factors to initiate transcription of the Pgk2 gene. Finally, our results facilitate a functional comparison of the promoter of the testis-specific Pgk2 gene with that of the ubiquitously expressed Pgk1 gene that allows us to suggest key molecular changes that accompanied the evolution of tissue-specific regulation of the Pgk2 gene.
|
|
|---|
60 dpp) mice. Testicular spermatozoa were isolated from the testes of adult mice by dissociation of tissue and sonication to lyse all other cells (38). The purity of each cell population was determined by morphological examination under phase optics. Purities of spermatogonial cell types and juvenile spermatocytes were
85%, while those of adult spermatocytes, spermatids, and spermatozoa were
95%. Spleen cells from adult CD-1 mice were used as a somatic control. Spleens were mechanically teased open with hypodermic needles in phosphate-buffered saline, and the resulting cell suspensions were filtered through 100-µm nylon mesh (Millipore, Billerica, MA) and then treated with Tris-buffered ammonium chloride (0.017 M Tris-HCl, pH 7.6, 0.84% NH4Cl) to lyse contaminating erythrocytes. Dimethyl sulfate and piperidine treatments of DNA for in vivo footprinting and ligation-mediated PCR (LM-PCR). In vivo dimethyl sulfate (DMS) treatment was performed with tissues or enriched populations of specific cell types, as previously described by Hornstra and Yang (24). The double-stranded linker (made from oligonucleotides LP1 and LP2 [Table 1]) used for annealing to the piperidine-cleaved blunt ends of DMS-treated genomic DNA has been described previously (43, 44). One pmol of gene-specific primer 1 (Table 1) was annealed to 1 µg of treated DNA and extended using Vent (exo–) DNA polymerase (New England Biolabs, Ipswich, MA). The reaction was performed for 1 cycle at 94°C for 5 min followed by 1 min at each of the following temperatures, in succession: 57°C, 55°C, 54°C, 53°C, 52°C, 51°C, and 50°C. Following the annealing steps, extension was performed at 72°C for 20 min. The double-stranded linker (150 pmol) was ligated to the genomic DNA fragments in an overnight ligation at 16°C with T4 DNA ligase (New England Biolabs). Following ligation, genomic DNA fragments were purified by ethanol precipitation and then amplified by using primer LP2 and a gene-specific primer 2 (Table 1). The PCR was performed as follows: 94°C for 12 min, followed by 20 cycles of 94°C for 30 s, 61°C for 2 min, and 76°C for 30 s. Following this amplification, 0.1 pmol of radiolabeled gene-specific primer 3 (Table 1) was added to each tube, and an additional five cycles were run, consisting of 94°C for 1 min, 69°C for 1 min, and 76°C for 2 min. The amplification products were purified by ethanol precipitation, resolubilized in formamide loading buffer, and denatured by boiling for 5 min. They were then loaded, along with an aliquot of treated naked DNA (G ladder), onto a 6% polyacrylamide gel and subjected to electrophoresis at 55 W of constant power for 3 h. Each gel was then dried and exposed overnight to an imaging screen (Bio-Rad, Hercules, CA) and visualized on a Molecular Imager FX (Bio-Rad).
|
View this table: [in a new window] |
TABLE 1. Primers for in vivo footprinting/ligation-mediated PCR
|
Real-time qPCR for ChIP assay. qPCR analysis was carried out with DNA precipitated by ChIP using SYBR Premix Ex Taq (Takara Bio USA, Madison, WI) on a Chromo4 real-time PCR detection system (Bio-Rad). Each PCR assay was performed in triplicate. PCR primers used are listed in Table 2. Relative amounts of DNA in the input and bound fractions were determined from cycle threshold (CT) values, using a standard curve generated from DNA of known concentrations. Values for enrichment were calculated as the average from at least two independent ChIP experiments and multiple independent PCR analyses with each antibody. Results were normalized to control values from rabbit IgG arbitrarily set at 1.0.
|
View this table: [in a new window] |
TABLE 2. Primers for PCRs used with ChIP assays
|
|
View this table: [in a new window] |
TABLE 3. Primers for quantitative real-time RT-PCRs
|
|
|
|---|
![]() View larger version (16K): [in a new window] |
FIG. 1. Quantitative RT-PCR analysis of Pgk gene expression during spermatogenesis. Expression of the Pgk1 (A) and Pgk2 (B) genes during spermatogenesis and in representative somatic tissues was quantified by real-time RT-PCR. The relative expression of target mRNA was calculated from the target CT values and ß-actin mRNA CT values, using the standard curve method. Each error bar indicates standard error of the mean. Populations of the specific spermatogenic cell types examined included the following: Type A, type A spermatogonia; Type B, type B spermatogonia; Prelep, preleptotene spermatocytes; Lep/Zyg, leptotene plus zygotene spermatocytes; Ea Pachs, early (juvenile) pachytene spermatocytes; Pachs, pachytene spermatocytes; Round, round spermatids; Elong/RBs, elongated spermatids plus residual bodies.
|
![]() View larger version (56K): [in a new window] |
FIG. 2. In vivo footprinting of the Pgk2 promoter during spermatogenesis in the mouse. (A) Schematic representation of putative regulatory elements in the Pgk2 gene promoter. The transcription start site is represented as a bent arrow. The core promoter includes a GC box and a CAAT box (open boxes). Immediately upstream from the core promoter is enhancer region E1/E2 and the E3/E4 subregions that have been shown to bind up to four anonymous factors. These sites are represented, respectively, as oval light and dark gray boxes. The location of amplimers analyzed by LM-PCR for in vivo footprinting are depicted as horizontal lines with numbers representing nucleotide positions relative to the transcription start site. Amplimers shown above and below the line represent analysis of upper and lower DNA strands, respectively. (B through E) DMS modification and LM-PCR were used to in vivo footprint the core promoter and enhancer regions of the Pgk2 promoter in purified populations of mouse spermatogenic cells. Open circles adjacent to bands represent protected cleavage sites, whereas filled circles represent enhanced cleavage sites. Numbers represent nucleotide positions relative to the transcriptional initiation site. Vertical lines represent the location of known or putative factor-binding elements. Prim. type A, primitive type A spermatogonia; type A, type A spermatogonia; type B, type B spermatogonia; Lep/Zyg, leptotene plus zygotene spermatocytes; early pachs, early (juvenile) pachytene spermatocytes; pachytene, pachytene spermatocytes; round, round spermatids; testis sperm, testicular sperm.
|
Because transcription of the Pgk2 gene is repressed in somatic cells, we expected to detect few protein-DNA interactions at the Pgk2 promoter in splenocytes other than those representing putative repressors. Indeed, we detected a binding activity at position –115 that extends to –110 and that overlaps the E3/E4 region, consistent with our previous EMSA data indicating binding of a single factor at the E3/E4 region in nonexpressing somatic cells (Fig. 2D) (21). This suggests that a repressor is constitutively bound to the Pgk2 enhancer in nonexpressing somatic cells; however, the identity of this repressor remains unknown. We also detected an enhancement at –12 in splenocytes, but we suspect that this represents nonspecific binding since it was detected in splenocytes, spermatocytes, and spermatids but not in spermatogonia.
Active transcription of the Pgk2 gene was previously demonstrated in primary spermatocytes and round spermatids by nuclear runoff experiments (29). We detected a variety of distinct footprints at the Pgk2 promoter site in vivo in spermatogenic cells. In the core promoter, a faint footprint first became evident at the GC box in preleptotene spermatocytes (extending from position –30 to –24) and became much more distinct in later spermatocytes (Fig. 2B). This footprint consisted of an enhancement at position –24 and several protections at –25, –27, –29, and –30 and developed gradually throughout the progression of first meiotic prophase, coinciding with the increase in Pgk2 transcriptional activity shown in Fig. 1B. In round spermatids, the footprint at the GC box diminished slightly (the protection at –29 disappeared) (Fig. 2B). This may indicate the beginning of a gradual unloading of a GC box-bound factor associated with changes in chromosomal architecture as histones are replaced by transition proteins and subsequently protamines during spermiogenesis (49).
The CAAT box (from position –105 to –100) lacks any guanosines within the putative binding sequence on the upper strand of the Pgk2 core promoter, such that factor binding in vivo cannot be assessed directly on this strand by using DMS as the footprinting agent (Fig. 2C). However, the guanosines (–105, –104, and –100) on the bottom strand also showed no significant footprint in spermatocytes and spermatids (Fig. 2E). Therefore, although a binding activity at the CAAT box was detected in vitro in EMSA experiments (21), we were not able to confirm a protein-DNA interaction at this site in vivo, even in cell types that were actively transcribing the Pgk2 gene.
In spermatocytes and round spermatids, protein-DNA interactions were also detected in the Pgk2 enhancer region, specifically in the E3/E4 region at about the same stage as factor binding developed in the core promoter (Fig. 2D). This was indicated by the appearance, initially in preleptotene spermatocytes, of an enhancement at position –115 and a protection at –110. Subsequently, in a mixture of leptotene and zygotene spermatocytes, an additional protection became visible at –119. These factor-binding activities persisted through the round spermatid stage. We conclude that this footprint is distinct from that seen in this same region in somatic splenocytes for two reasons. First, the footprint detectable in spermatocytes extends from –110 to –119 and, as such, is more expansive than the footprint detected in splenocytes that extends from –110 to only –115. Second, the footprint detected in splenocytes indicative of a putative repressor is not observed in spermatogonia, indicating that the putative repressor becomes dissociated from the E3/E4 region prior to activation and is subsequently replaced by an activator that binds to this same region in spermatocytes and spermatids.
While we were able to detect a protein-DNA interaction at the E3/E4 region of the Pgk2 enhancer in expressing spermatocytes and spermatids in vivo, we detected no binding activity at the E1/E2 region in these same cell types (Fig. 2D). Thus, as with the CAAT box, we found that certain sequence elements that demonstrated the capacity to bind nuclear factors in reconstitution EMSA experiments in vitro are, in fact, not bound by proteins in vivo. This demonstrates the need to confirm results obtained in vitro with analyses performed in vivo.
Interestingly, even though the GC box is bound in spermatocytes and spermatids by a ubiquitously expressed factor, we detected no evidence of protein binding at this site in any of the spermatogonial cell types we examined (Fig. 2B). This suggests that the Pgk2 promoter exists in a condensed (repressed) chromatin state in spermatogonia, such that factors that are already present at this stage are prevented from accessing their cognate binding sites. Similarly, we observed no evidence of factor binding in the Pgk2 enhancer region in spermatogonia (Fig. 2D). However, we did detect several apparent protein-DNA interactions in the form of protections in the region between the CAAT box and the GC box, specifically on the lower strand in primitive type A and type A spermatogonia (at positions –100, –93, –80, –74, –70, –56, and –40) and, to a lesser extent, in type B spermatogonia (only at positions –74, –56, and –40) (Fig. 2E). The pattern of these footprints suggests a periodicity of about 20 bp that could represent a phased nucleosome, but we have no other data to confirm this. If a phased nucleosome is present in this position in spermatogonia, then it apparently becomes displaced coincident with the subsequent initiation of active transcription of the Pgk2 gene in spermatocytes and spermatids.
Postmeiotic spermatogenic cells undergo the unique differentiative process of spermiogenesis that includes, among other things, significant chromatin condensation as histones associated with DNA in earlier spermatogenic cells are ultimately replaced by protamines (49). As a result, global transcription ceases, presumably due to restricted access of the transcriptional machinery to promoter sequences. Therefore, we expected to observe a reduction in protein-DNA interactions associated with the Pgk2 promoter in spermatozoa. Indeed, the protections and enhancements seen at the GC box and enhancer region in expressing spermatocytes and spermatids largely disappeared in testicular spermatozoa, wherein we observed a footprint pattern similar to that seen in nonexpressing somatic cells (Fig. 2B to D).
Identification by ChIP analysis of factors regulating Pgk gene transcription in vivo. The GC box is well characterized as an SP/KLF family binding site and is a common cis regulatory element present in the promoters of several genes expressed in a developmentally specific manner in differentiating mouse germ cells (11, 25, 54, 56, 61). The SP/KLF family of transcription factors is united by a particular combination of three conserved C2H2 zinc fingers that form the sequence-specific DNA-binding domain (31). Within the SP/KLF family of transcription factors, SP1 and SP3 are ubiquitously expressed in mammalian cells, including spermatogenic cells (56) and so were primary candidates for the factors responsible for the footprint at the Pgk2 GC box in spermatocytes and spermatids.
To examine the association of the SP1 and SP3 factors with GC boxes in the Pgk1 and Pgk2 promoters in vivo, we conducted ChIP experiments on cells from adult mouse testis and spleen. In contrast to the Pgk2 promoter, which contains a single GC box, the murine Pgk1 promoter has five GC box elements located within a CpG island (as revealed by TFSEARCH [http://www.cbrc.jp/research/db/TFSEARCH.html]). We observed an association of both SP1 and SP3 with the Pgk1 promoter in somatic splenocytes (Fig. 3A and B). Similarly, in adult mouse testis, which contains both somatic and germ cells, we also detected both SP1 and SP3 in association with the Pgk1 promoter (Fig. 3A and B). In contrast, the Pgk2 promoter showed only background levels of association with either factor in splenocytes (where this gene is not expressed) but showed a significant association with SP3, but not SP1, in the testis (where this gene is expressed) (Fig. 3A and B).
![]() View larger version (26K): [in a new window] |
FIG. 3. In vivo binding of transcription factors to the Pgk1 and Pgk2 promoters. ChIP analyses were performed with testicular germ cells and splenocytes, using antibodies specific for SP1 (A), SP3 (B), PBX1/2/3/4 (C), PBX1/2/3 (D), PREP1 (E), CREB (F), CREM (G), and CBF-B/NF-YA (H). Precipitated DNA was subjected to real-time PCR with the use of specific primers for the Pgk1 and Pgk2 promoters (see below), and relative enrichment of amplified regions in immunoprecipitated DNA (bound fractions; BO) was determined and compared to that in input (IP) DNA. All results were normalized to values for normal rabbit IgG, with that level arbitrarily set to 1.0. Each error bar indicates standard error of the mean. All PCR primers used are listed in Table 2. The PCR primer set named "Pgk1 CAAT box" was used to analyze the binding of all transcription factors examined in the Pgk1 promoter. This primer set amplified the "NF1-like region" (from position –180 to –77) (14, 47), including the consensus binding site for PBX factors and the CAAT box element (the binding site for CBF-B/NF-YA) in the Pgk1 promoter. Because of the size of the sonication fragments, amplification of this region was also used to examine the binding of SP factors to the nearby GC boxes at positions –66 and –254 in the Pgk1 promoter. Similarly, the PCR primer sets named "Pgk2 E1-E4," "Pgk2 CAAT box," and "Pgk2 GC box" (Table 2) were used to selectively amplify specific regions in the Pgk2 promoter, including a region containing the E1/E2-plus-E3/E4 enhancer sites, a region containing the CAAT box, and a region containing the GC box, respectively. The "Pgk2 GC box" primer set was also used to analyze the binding of CREB and/or CREM to the putative CRE site (+333). Deduction of specificity of binding of each factor to each element was based on the fact that, except where noted otherwise, each putative binding element represented the only candidate consensus binding sequence for each factor within at least 1 kb in each direction in each promoter.
|
We used ChIP to investigate the interactions in vivo between PBX4 and PREP1 and the Pgk2 enhancer elements E3/E4. Because an antibody that solely recognizes PBX4 was not commercially available, we performed ChIP assays using an anti-PBX1/2/3/4 antibody, which recognizes all four isoforms of the PBX family, and an anti-PBX1/2/3 antibody that recognizes the first three but not the fourth isoform (Fig. 3C and D). The anti-PBX1/2/3/4 antibody precipitated Pgk1 sequences from both spleen and testis and Pgk2 sequences from testis only. The anti-PBX1/2/3 antibody also precipitated Pgk1 sequences at high levels from spleen but only at very low levels from testis. Pgk2 sequences were not precipitated from either spleen or testis, using the anti-PBX1/2/3 antibody. This suggests that PBX4 is the sole PBX isoform bound to the Pgk2 gene and is the predominant isoform bound to the Pgk1 gene in the testis. We also confirmed that PREP1 localizes to the Pgk2 promoter in the testis, but not in spleen, and to the Pgk1 gene in both spleen and testis (Fig. 3E). These associations support the contention that PBX factors function as activators in conjunction with the coactivator PREP1 to stimulate transcription from Pgk promoters.
Although they were not indicated by our footprinting experiments, two other candidate regulators of Pgk transcription were the cyclic AMP-responsive element modulator (CREM) and its family member the cyclic AMP-responsive element binding protein (CREB). CREM is a transcription factor known to be critically required for spermatogenesis (26). A microarray study of gene expression in testes of CREM-deficient mice revealed a 2.8-fold decrease in the level of Pgk2 transcripts (http://www.dkfz.de/tbi_old/crem/index.html) (6). Indeed, half CRE sites (5'-TGACG-3' and/or 5'-CGTCA-3') are found at –1632 and +333 of the murine Pgk2 gene (http://natural.salk.edu/CREB/) (64). CREB is known to function in glucose homeostasis, for which the glycolytic enzyme phosphoglycerate kinase, encoded by the Pgk1 and Pgk2 genes, is crucial (33). Previous studies showed that the human PGK1 promoter is occupied by CREB in primary hepatocytes and in the human embryonic kidney cell line HEK293T (64).
We examined the association of CREM and CREB with the Pgk1 and Pgk2 promoters in testis and spleen. We confirmed the binding of CREB to the Pgk1 promoter in both spleen and testis, but no significant association of CREB with the Pgk2 promoter was observed in either tissue (Fig. 3F). However, CREM, a master regulator of spermatogenesis (26), was found to be associated with both the Pgk1 and Pgk2 promoters in the testis (Fig. 3G). Thus, both CREB and CREM are bound to the Pgk1 gene in the testis, but only CREM is bound to the Pgk2 gene in this same tissue. There are both activating and repressing isoforms of CREM (26). Because a microarray study of gene expression during spermatogenesis in a CREM-deficient mouse revealed a 2.8-fold decrease in expression of the Pgk2 gene (6), we assume there are normally activating isoforms of CREM bound to this gene in spermatocytes and spermatids in vivo.
Finally, we examined the association of CBF-B/NF-YA in both the Pgk1 and Pgk2 promoters, because it is the most likely factor to bind to a CAAT box (55). Pfeifer et al. (46) and Dai et al. (14) used genomic sequencing to detect a footprint in vivo at the CAAT box in the actively expressed human PGK1 gene, and our ChIP analysis revealed an association of CBF-B/NF-YA with the mouse Pgk1 promoter in somatic (spleen) cells where the Pgk1 gene is actively expressed (Fig. 3H). However, as noted above, our in vivo footprinting analysis of the Pgk2 gene in spermatocytes and spermatids revealed no protein-DNA interaction at the Pgk2 CAAT box (Fig. 2C and E). This is consistent with our ChIP analysis that revealed no association of CBF-B/NF-YA with the Pgk2 promoter in testicular cells (Fig. 3H).
Spermatogenic stage-specific binding of transcription factors to the Pgk2 gene. Having identified the key transcriptional regulators of the Pgk2 gene, we next sought to determine the developmental kinetics by which these factors interact with the Pgk2 promoter during spermatogenesis in vivo. We performed ChIP assays for PBX4, SP3, and CREM on highly enriched populations of specific mouse spermatogenic cell types at stages prior to, during, and following the period of active transcription of the Pgk2 gene. As shown in Fig. 1B, transcription of the Pgk2 gene begins at the onset of meiosis, and levels of Pgk2 transcript then increase steadily throughout spermatogenesis. On the other hand, the X-linked Pgk1 gene is expressed in all somatic cells and in premeiotic spermatogonia, but transcription of this gene ceases in primary spermatocytes due to MSCI (57) and continues to be repressed in spermatids due to the effects of PMSC (Fig. 1A) (45, 60). We found that PBX4, SP3, and CREM are all recruited to the Pgk2 promoter predominantly at the onset of meiosis and then remain present throughout meiosis and the subsequent differentiation of spermatogenic cells to the round spermatid stage (Fig. 4A to C, left side), coinciding exactly with the pattern of active transcription of the Pgk2 gene during spermatogenesis.
![]() View larger version (27K): [in a new window] |
FIG. 4. Developmental stage- and cell type-specific recruitment of transcription factors to the Pgk1 and Pgk2 promoters. ChIP analyses were performed with relatively pure populations of specific spermatogenic cell types prior to, during, and following the period of active transcription of the Pgk2 gene, using antibodies specific for PBX1/2/3/4 (A), SP3 (B), and CREM (C). Precipitated DNA was subjected to real-time PCR with the use of specific primers for the Pgk2 (left side panels) and Pgk1 (right side panels) promoters, and enrichment of immunoprecipitated DNA (bound fractions, BO) relative to input (IP) DNA was calculated in each case. All results were normalized to values for normal rabbit IgG, with that level arbitrarily set to 1.0. Each error bar indicates standard error of the mean. Type A, type A spermatogonia; Type B, type B spermatogonia; Prelep, preleptotene spermatocytes; Lep/Zyg, leptotene plus zygotene spermatocytes; Ea Pachs, early (juvenile) pachytene spermatocytes; Pachs, pachytene spermatocytes; Round, round spermatids.
|
Expression profile of transcriptional activators of the Pgk2 gene during spermatogenesis. Our factor- and stage-specific ChIP assays revealed the identities of key factors regulating transcription of the Pgk2 gene but did not reveal which of these factors might be the primary regulator(s) of the tissue- and stage-specific transcription of this gene during spermatogenesis. To gain further insight into this aspect, we analyzed expression profiles during spermatogenesis of Pbx4 and Prep1, along with other members of the Pbx (Pbx1, Pbx2, Pbx3) and Meis (Meis1, Meis2, Meis3, and Prep2) families of genes encoding transcription factors. Using real-time RT-PCR, we examined total RNA from highly enriched populations of specific spermatogenic cell types, as well as from whole adult testis, liver, spleen, and brain tissues (Fig. 5). The relative expression levels of Pbx1a and Pbx1b, alternative splicing variants of transcripts from the Pbx1 gene, were consistently low at all stages of spermatogenesis and in all tissues examined, except for the expression level of Pbx1b in the brain (Fig. 5A and B), while the Pbx2 and Pbx3 genes were expressed at various levels throughout spermatogenesis, with the Pbx2 expression level showing a decrease during meiosis (Fig. 5C and D). However, the testis-specific Pbx4 gene showed a tightly regulated expression pattern during spermatogenesis, initiating transcription in meiotic spermatocytes and continuing in postmeiotic spermatids (Fig. 5E), mirroring the initial expression profile of the Pgk2 gene (Fig. 1B). Meis1 was also expressed at various levels during spermatogenesis (Fig. 5F), while the Meis2 and Meis3 genes were repressed in male germ cells, although these genes were expressed at high levels in the brain (Fig. 5G and H). Expression of the Prep1 and Prep2 genes was upregulated in the testis. Prep1 showed increasing expression beginning at the onset of meiosis in a pattern similar to that of the Pbx4 gene (Fig. 5I), whereas Prep2 was expressed almost exclusively in postmeiotic spermatids (Fig. 5J). These findings are complemented by the results of previous studies in which PBX4 protein was detected in a subnuclear domain in pachytene spermatocytes (59) and the PREP1 protein was detected in the adult testis (19). Because the expression patterns of the Pbx4 and Prep1 genes mirror that of Pgk2 and because our ChIP analysis demonstrates that these factors associate with the Pgk2 promoter in an appropriate stage-specific pattern, we suggest that a complex of PBX4 and PREP1 is likely to be the primary factor involved in initiating stage- and cell type-specific transcription of the Pgk2 gene during spermatogenesis.
![]() View larger version (26K): [in a new window] |
FIG. 5. qRT-PCR analysis of expression during spermatogenesis of genes encoding PBX and MEIS factors. Expression of different members of the Pbx family of genes, Pbx1a (A), Pbx1b (B), Pbx2 (C), Pbx3 (D), and Pbx4 (E), and the Meis family of genes, Meis1 (F), Meis2 (G), Meis3 (H), Prep1 (I), and Prep2 (J), was quantitated by real-time RT-PCR. The relative expression of each target mRNA was calculated from the target CT values and the 18S rRNA CT values, using the standard curve method. Each error bar indicates standard error of the mean. Type A, type A spermatogonia; Type B, type B spermatogonia; Prelep, preleptotene spermatocytes; Lep/Zyg, leptotene plus zygotene spermatocytes; Ea Pachs, early (juvenile) pachytene spermatocytes; Pachs, pachytene spermatocytes; Round, round spermatids; Elong/RBs, elongated spermatids plus residual bodies.
|
|
|
|---|
Although transcription of the Pgk2 gene does not begin until the meiotic spermatocyte stage (Fig. 1B), we observed unique footprints within the Pgk2 promoter region in premeiotic spermatogonia (Fig. 2E). These footprints, located between the CAAT and GC boxes in the Pgk2 core promoter, were observed for premeiotic spermatogonia but not for subsequent spermatocytes or spermatids and were most prominent in primitive type A spermatogonia isolated from prepuberal males at 6 dpp and in type A spermatogonia isolated at 8 dpp and less prevalent in differentiating type B spermatogonia, also isolated from males at 8 dpp. Because these footprints were observed prior to but not during active transcription of the Pgk2 gene, we envision two potential explanations for what they represent. One possibility, supported by the
20-bp periodicity of these footprints, is that they represent the presence of a single, phased nucleosome positioned in this region just prior to the activation of transcription that then becomes displaced to facilitate factor binding required for the initiation of transcription. A second possibility is that these footprints could represent the presence of chromatin remodeling complexes involved in initiating critical changes in chromatin structure to potentiate the initiation of transcription of this gene in spermatocytes, although detectable decondensation of chromatin and formation of a DNase I hypersensitive site in the Pgk2 promoter do not occur until late in the spermatogonia type B stage (28, 29).
A distinct set of footprints was detected uniquely in those spermatogenic cell types in which the Pgk2 gene is known to be actively transcribed (Fig. 1B) (29, 38). Thus, in meiotic spermatocytes and postmeiotic round spermatids, we observed unique footprints at the GC box and in the E3/E4 enhancer region (Fig. 2B and D). Based on our ChIP assays, we report here, for the first time, the putative identity of the transcription factors responsible for initiation of transcription of the Pgk2 gene during spermatogenesis. We find that the SP3 factor is preferentially bound to the GC box in the Pgk2 core promoter in spermatocytes and spermatids and that the PBX4 factor is bound to the E3/E4 region in the Pgk2 enhancer in these same cells (Fig. 3 and 4). Mammalian PBX factors have been shown to independently dimerize with either HOX proteins or MEIS/PREP proteins (41). Trimeric complexes encompassing all three homeoproteins, PBX-HOX-MEIS/PREP, have also been characterized (41). Our ChIP assays detected PREP1 at the Pgk2 promoter in the testis (Fig. 3). Therefore, initiation of Pgk2 transcription appears to be driven by a complex involving the PBX4 and PREP1 factors bound to the E3/E4 enhancer region, and this complex presumably interacts with other factors bound to the Pgk2 promoter, including SP3 and CREM, to form an enhanceosome-like structure that directs the initiation of transcription of the Pgk2 gene (62). Interestingly, the activity of certain tissue-specific genes activated by PBX-containing complexes was previously shown to be regulated by a balance, or switch, between activator and repressor complexes (5, 51). The alternate protein-DNA interactions we observed at the E3/E4 enhancer in somatic cells, where the Pgk2 gene is repressed, and spermatocytes and spermatids, where the Pgk2 gene is actively transcribed, could represent the function of such a switch mechanism.
Our developmental analysis demonstrated that occupation of the Pgk2 enhancer site by the PBX4-PREP1 complex first occurs coincidentally with the onset of meiosis in primary spermatocytes and that the initial occupation of the GC box by SP3 in the Pgk2 core promoter and the association of CREM with the Pgk2 gene also occur at approximately this same time. This is despite the fact that SP3 and CREM are present in earlier, premeiotic spermatogonia, as are other members of the PBX and MEIS transcription factor families. This suggests that prior to the early primary spermatocyte stage, the Pgk2 promoter is inaccessible to binding by most transcription factors, consistent with our previous observations that chromatin in this region remains relatively condensed in spermatogonia and then decondenses in spermatocytes (28, 29). The stage-specific expression patterns during spermatogenesis of the genes encoding the PBX4 and PREP1 factors suggest that the initial formation of this complex directly regulates the coincident, stage-specific transcriptional activation of the Pgk2 gene.
In this regard, it is interesting to note that somatic isoforms of the PBX factor family have been implicated in the activation of muscle-specific gene expression (8, 15, 27). In muscle, PBX forms a complex with MEIS and penetrates repressive chromatin in the myogenin promoter to recruit MYOD, which then interacts with histone acetyltransferases to promote the acetylation of histones and the association of SWI/SNF complexes that then initiate chromatin remodeling, and this, in turn, facilitates access for stable binding of additional activator complexes to the myogenin promoter (8, 15). Indeed, it has been suggested that PBX-related complexes act as "pioneer" transcription factors to enter repressed (condensed) chromatin and initiate remodeling to facilitate initiation of transcription of tissue-specific genes (8, 15, 16, 50). Our observations of molecular changes in the Pgk2 promoter at the beginning of the spermatocyte stage, including the disappearance of a potential phased nucleosome (Fig. 2E), decondensation of chromatin (28, 29), and changes in histone modifications (H. Yoshioka, J. L. Hornecker, and J. R. McCarrey, unpublished data), all support the contention that chromatin remodeling is indeed induced in the Pgk2 promoter at this stage. The contention that binding of the PBX4-PREP1 complex to the Pgk2 enhancer is critical for the activation of Pgk2 transcription is further supported by our earlier observations that the E3/E4 region is required to direct the testis-specific expression in vivo of Pgk2 promoter/CAT reporter transgenes (48, 62) and that this region is also required to direct changes in histone modifications in the Pgk2 promoter (H. Yoshioka et al., unpublished data). Thus, these observations are consistent with a model in which the PBX4-PREP1 complex binds to the Pgk2 promoter and recruits histone-modifying enzymes and chromatin remodeling factors to facilitate decondensation of chromatin, reconfiguration of nucleosomes, and binding of transcription factors including SP3 and CREM to initiate spermatogenesis-specific transcription of the Pgk2 gene.
Our in vivo analysis of protein-DNA interactions in specific spermatogenic cell types allowed us to examine the developmental sequence of events leading to the initiation of transcription of the Pgk2 gene. Previous developmental studies of molecular events at the promoters of the HNF-4
(23), the collagenase (32), the IFN-ß (1), and the chicken lysozyme (30) genes suggest that the order in which histone modifications, chromatin remodeling, and factor binding normally occur varies among different tissue-specific genes undergoing transcriptional activation. Our results indicate that recruitment of the PBX4-PREP1 complex, SP3, and CREM to the Pgk2 promoter occurs at approximately the same time as chromatin decondensation (28) and the formation of a DNase I hypersensitive site (29), as well as changes in histone modifications (H. Yoshioka et al., unpublished data), all of which occur coincident with the spermatogonium-spermatocyte transition. Importantly, however, the resolution of our developmental analyses in vivo was limited to days rather than minutes or hours. Thus, we cannot distinguish between truly simultaneous events and sequential events that occur in rapid succession. Nevertheless, we have shown that these parameters all undergo significant changes at this stage and none (other than demethylation of DNA [22]) is preset significantly prior to the stage at which transcription of the Pgk2 gene is first initiated.
As noted in the introduction, the intronless Pgk2 gene is believed to have arisen in the mammalian genome as a retroposed copy of an intron-containing Pgk1-like progenitor gene early during mammalian evolution (12, 40). Transcription of the newly formed Pgk2 gene was most likely facilitated by the inclusion of a portion of the promoter from the Pgk1-like progenitor gene as part of the retroposed sequence (34-36). Thus, we previously speculated that the Pgk2 gene initially carried a housekeeping promoter sequence that directed ubiquitous low-level transcription and that this regulatory sequence has subsequently undergone a transition to direct, high-level, testis-specific transcription, at least in eutherian mammals (34-36). A comparison of mammalian Pgk1 and Pgk2 promoter sequences is shown in Fig. 6. One potentially important distinction not specifically shown in this figure is that the Pgk1 promoter is embedded in a CpG island, whereas the Pgk2 promoter has lost the density of CpG dinucleotides required to be considered an island and has retained only dispersed CpG dinucleotides (35). This difference is likely responsible in part for the fact that the Pgk1 promoter region typically remains constitutively hypomethylated (except on the inactive X chromosome) (47), whereas the Pgk2 promoter undergoes highly regulated demethylation in a tissue-, stage-, and cell type-specific manner uniquely in spermatogenic cells (22).
![]() View larger version (27K): [in a new window] |
FIG. 6. Pgk promoter sequences in mammals. Comparisons are shown among promoter sequences for the Pgk1 (A) and Pgk2 (B) genes from various mammalian species. Alignment of sequences upstream from the translational start codon (ATG, dashed rectangle) of each Pgk gene was done using CLUSTALW (http://align.genome.jp/). Identical nucleotides are represented by asterisks, and dashes indicate a gap. The numbers to the right of each sequence are relative to the transcription start site, which is designated +1. Highly conserved sequences for which potential binding functions are known are surrounded by rectangles. (C) An alignment is shown of the NF1-like region from the human and mouse Pgk1 genes and the E3/E4 region from the human and mouse Pgk2 genes and the potential relationship among these regions. Asterisks indicate sequence identity, dashes indicate a lack of identity, and vertical lines indicate partial identity.
|
A key distinction between the Pgk1 and Pgk2 promoters is found in the region just upstream from the CAAT box in each. In the Pgk1 promoter, there is a region of 41 bp that has been termed an "NF1-like" binding sequence and has been shown to bind factors found in somatic cells (14, 46). This sequence is relatively well conserved among the four mammalian Pgk1 promoters, as shown in Fig. 6A, but is not found intact in the Pgk2 promoter. Rather, in the equivalent position, the Pgk2 gene shows a different conserved element of 18 bp that we previously termed the E3/E4 region (Fig. 6B) and showed by EMSA to have the capacity to bind one or more testis-specific factors (21). The results presented in this study indicate that this is the site of binding of the testis-specific PBX4-PREP1 complex. Interestingly, the E3/E4 enhancer element in the Pgk2 promoter shows some similarity to the 5' half of the NF1-like element in the Pgk1 promoter (Fig. 6C). We assume that this portion of the Pgk1 promoter is responsible for binding somatic members of the PBX family, as indicated by precipitation of the Pgk1 promoter by the PBX1/2/3 antibody in splenocytes (Fig. 3C and D). However, the extensive conservation of sequence in the E3/E4 region among mammalian Pgk2 genes suggests that it has undergone selection for optimal binding of the PBX4-PREP1 complex.
This information allows us to speculate that the Pgk1-like promoter included in the original Pgk2 retroposon may have diverged to facilitate its transition from a housekeeping promoter to its present status as a tissue-specific promoter which functions uniquely in spermatogenic cells. It is noteworthy that the differences in protein-DNA interactions associated with each of the Pgk promoters are based largely on the preferential binding of different members of similar families of transcription factors. In particular, it would appear that the evolution of testis-specific members of the PBX activator and MEIS/PREP coactivator families, along with the optimization of a binding sequence to preferentially bind the testis-specific isoform of the PBX activator family (PBX4), have combined to impose testis-specific regulation onto the Pgk2 retroposon in eutherian mammals. As genomic sequences for additional mammalian species become available, it will be of interest to recapitulate this evolutionary process to determine if, as we suspect, the appearance of testis-specific members of the PBX and MEIS factor families coevolved in mammals with the appearance of an optimal binding site for this complex in the Pgk2 promoter.
In conclusion, we have identified key transcription factors that are associated with the actively transcribed Pgk2 gene in spermatocytes and spermatids, including CREM, SP3, and the testis-specific PBX4 activator and PREP1 coactivator. We have exploited the highly accessible spermatogenic cell lineage to characterize the developmental kinetics of associations between these factors and the Pgk2 promoter. Finally, we have suggested how the Pgk2 promoter may have evolved from a Pgk1-like housekeeping promoter to the tightly regulated, testis-specific promoter now associated with this gene in eutherian mammals.
This work was supported by NIH grant HD46637 to J.R.M. K.T.P. was supported by NIH MBRS-RISE grant GM60655.
Published ahead of print on 17 September 2007. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»