GC-rich Motifs Research Articles

IGS sequences in Cestrum present AT- and GC-rich conserved domains, with strong regulatory potential for 5S rDNA.

The 35S and 5S ribosomal DNA (rDNA) organized in thousands of copies in genomes, have been widely used in numerous comparative cytogenetic studies. Nevertheless, several questions related to the diversity and organization of regulatory motifs in 5S rDNA remain to be addressed. The 5S rDNA unit is composed of a conserved 120bp length coding region and an intergenic spacer (IGS) containing potential regulatory motifs (Poly-T, AT-rich and GC-rich) differing in number, redundancy and position along the IGS. The Cestrum species (Solanaceae) have large genomes (about 10pg/1C) and conserved 2n = 16 karyotypes. Strikingly, these genomes show high diversity of heterochromatin distribution, variability in 35S rDNA loci and the occurrence of B chromosomes. However, the 5S rDNA loci are highly conserved in the proximal region of chromosome 8. Comparison of seventy-one IGS sequences in plants revealed several conserved motifs with potential regulatory function. The AT- and GC-rich domains appeared highly conserved in Cestrum chromosomes. The 5S genic and the GC-rich IGS probe produced FISH signals in both A (pair 8) and B chromosomes. The GC-rich domain presented a strong potential for regulation because it may be associated with CpG islands organization, as well as to hairpin and loop organization. Another interesting aspect was the ability of AT- and GC-rich motifs to produce non-heterochromatic CMA/DAPI signals. While the length of the 5S rDNA IGS region varied in size between the Cestrum species, the individual sequence motifs seem to be conserved suggesting their regulatory function. The most striking feature was the conserved GC-rich domain in Cestrum, which is recognized as a signature trait of the proximal region of chromosome pair 8.

Inhibition of cancer cell proliferation by adenosine triphosphate-triggered codelivery system of p53 gene and doxorubicin

The codelivery of drugs and genes by stimuli-responsive nanocarriers is a promising strategy for achieving an effective cancer treatment. In this study, an adenosine triphosphate (ATP)-responsive nanosystem was constructed to codeliver doxorubicin (DOX) and p53 gene based on an ATP-triggered aptamer and polyethyleneimine (PEI). In this system, DOX interacts with the GC-rich motif of duplex formed by the aptamer and its cDNA sequence. Then, a ternary nanocomplex DOX-Duplex/PEI/p53 was constructed using cationic carrier PEI25K for facilitating the intracellular delivery and release of p53 gene and DOX. The DOX-Duplex/PEI/p53 nanocomplex was found to possess an efficient anticancer effect, which was attributed to the ability of the system to trigger cell apoptosis and meanwhile block the cell cycle at G2 phase. The favorable antiproliferative effect was found to be associated with the rapid DOX and p53 gene release in response to the intracellular ATP concentration and the synergistic effect of therapeutic drug and gene.

Intragenomic heterogeneity of intergenic ribosomal DNA spacers in Cucurbita moschata is determined by DNA minisatellites with variable potential to form non-canonical DNA conformations.

The intergenic spacer (IGS) of rDNA is frequently built of long blocks of tandem repeats. To estimate the intragenomic variability of such knotty regions, we employed PacBio sequencing of the Cucurbita moschata genome, in which thousands of rDNA copies are distributed across a number of loci. The rRNA coding regions are highly conserved, indicating intensive interlocus homogenization and/or high selection pressure. However, the IGS exhibits high intragenomic structural diversity. Two repeated blocks, R1 (300–1250 bp) and R2 (290–643 bp), account for most of the IGS variation. They exhibit minisatellite-like features built of multiple periodically spaced short GC-rich sequence motifs with the potential to adopt non-canonical DNA conformations, G-quadruplex-folded and left-handed Z-DNA. The mutual arrangement of these motifs can be used to classify IGS variants into five structural families. Subtle polymorphisms exist within each family due to a variable number of repeats, suggesting the coexistence of an enormous number of IGS variants. The substantial length and structural heterogeneity of IGS minisatellites suggests that the tempo of their divergence exceeds the tempo of the homogenization of rDNA arrays. As frequently occurring among plants, we hypothesize that their instability may influence transcription regulation and/or destabilize rDNA units, possibly spreading them across the genome.

GM604 regulates developmental neurogenesis pathways and the expression of genes associated with amyotrophic lateral sclerosis

BackgroundAmyotrophic lateral sclerosis (ALS) is currently an incurable disease without highly effective pharmacological treatments. The peptide drug GM604 (GM6 or Alirinetide) was developed as a candidate ALS therapy, which has demonstrated safety and good drug-like properties with a favorable pharmacokinetic profile. GM6 is hypothesized to bolster neuron survival through the multi-target regulation of developmental pathways, but mechanisms of action are not fully understood.MethodsThis study used RNA-seq to evaluate transcriptome responses in SH-SY5Y neuroblastoma cells following GM6 treatment (6, 24 and 48 h).ResultsWe identified 2867 protein-coding genes with expression significantly altered by GM6 (FDR < 0.10). Early (6 h) responses included up-regulation of Notch and hedgehog signaling components, with increased expression of developmental genes mediating neurogenesis and axon growth. Prolonged GM6 treatment (24 and 48 h) altered the expression of genes contributing to cell adhesion and the extracellular matrix. GM6 further down-regulated the expression of genes associated with mitochondria, inflammatory responses, mRNA processing and chromatin organization. GM6-increased genes were located near GC-rich motifs interacting with C2H2 zinc finger transcription factors, whereas GM6-decreased genes were located near AT-rich motifs associated with helix-turn-helix homeodomain factors. Such motifs interacted with a diverse network of transcription factors encoded by GM6-regulated genes (STAT3, HOXD11, HES7, GLI1). We identified 77 ALS-associated genes with expression significantly altered by GM6 treatment (FDR < 0.10), which were known to function in neurogenesis, axon guidance and the intrinsic apoptosis pathway.ConclusionsOur findings support the hypothesis that GM6 acts through developmental-stage pathways to influence neuron survival. Gene expression responses were consistent with neurotrophic effects, ECM modulation, and activation of the Notch and hedgehog neurodevelopmental pathways. This multifaceted mechanism of action is unique among existing ALS drug candidates and may be applicable to multiple neurodegenerative diseases.

Curcumin ameliorates PRMT5-MEP50 arginine methyltransferase expression by decreasing the Sp1 and NF-YA transcription factors in the A549 and MCF-7 cells.

The protein arginine methyltransferase 5 (PRMT5) and its catalytic partner methylosome protein MEP50 (WDR77) catalyse the mono- and symmetric di-methylation of selective arginines in various histones and non-histone target proteins. It has emerged as a crucial epigenetic regulator in cell proliferation and differentiation; which also reported to be overexpressed in many forms of cancers in humans. In this study, we aimed to assess the modulations in the expression of this enzyme upon exposure to the well-studied natural compound from the spice turmeric, curcumin. We exposed the lung and breast cancer cell lines(A549 and MCF-7) to curcumin (2 and 20μM) and observed a highly significant inhibitory effect on the expression of both PRMT5 and MEP50. The level of symmetrical dimethylarginine (SDMA) in multiple proteins, and more specifically, the H4R3me2s mark (which predominates in GC-rich motifs in nucleosomal DNA) was also diminished significantly. We also found that curcumin significantly reduced the level and enrichment of the transcription factors Sp1 and NF-YA which shares their binding sites within the GC-rich region of the PRMT5 proximal promoter. Furthermore, the involvement of both PKC-p38-ERK-cFos and AKT-mTOR signalling was observed in reducing the Sp1 and NF-YA expression by curcumin. Therefore, we propose curcumin decreased the expression of PRMT5 in these cells by affecting at least these two transcription factors. Altogether, we report a new molecular target of curcumin and further elucidation of this proposed mechanism through which curcumin affects the PRMT5-MEP50 methyltransferase expression might be explored for its therapeutic application.

Supramolecular aptamer nano-constructs for receptor-mediated targeting and light-triggered release of chemotherapeutics into cancer cells

Platforms for targeted drug-delivery must simultaneously exhibit serum stability, efficient directed cell internalization, and triggered drug release. Here, using lipid-mediated self-assembly of aptamers, we combine multiple structural motifs into a single nanoconstruct that targets hepatocyte growth factor receptor (cMet). The nanocarrier consists of lipidated versions of a cMet-binding aptamer and a separate lipidated GC-rich DNA hairpin motif loaded with intercalated doxorubicin. Multiple 2′,6′-dimethylazobenzene moieties are incorporated into the doxorubicin-binding motif to trigger the release of the chemotherapeutics by photoisomerization. The lipidated DNA scaffolds self-assemble into spherical hybrid-nanoconstructs that specifically bind cMet. The combined features of the nanocarriers increase serum nuclease resistance, favor their import into cells presumably mediated by endocytosis, and allow selective photo-release of the chemotherapeutic into the targeted cells. cMet-expressing H1838 tumor cells specifically internalize drug-loaded nanoconstructs, and subsequent UV exposure enhances cell mortality. This modular approach thus paves the way for novel classes of powerful aptamer-based therapeutics.

Structural basis for genome wide recognition of 5-bp GC motifs by SMAD transcription factors

Smad transcription factors activated by TGF-β or by BMP receptors form trimeric complexes with Smad4 to target specific genes for cell fate regulation. The CAGAC motif has been considered as the main binding element for Smad2/3/4, whereas Smad1/5/8 have been thought to preferentially bind GC-rich elements. However, chromatin immunoprecipitation analysis in embryonic stem cells showed extensive binding of Smad2/3/4 to GC-rich cis-regulatory elements. Here, we present the structural basis for specific binding of Smad3 and Smad4 to GC-rich motifs in the goosecoid promoter, a nodal-regulated differentiation gene. The structures revealed a 5-bp consensus sequence GGC(GC)|(CG) as the binding site for both TGF-β and BMP-activated Smads and for Smad4. These 5GC motifs are highly represented as clusters in Smad-bound regions genome-wide. Our results provide a basis for understanding the functional adaptability of Smads in different cellular contexts, and their dependence on lineage-determining transcription factors to target specific genes in TGF-β and BMP pathways.

In-silico analysis of cis-acting regulatory elements of pathogenesis-related proteins of Arabidopsis thaliana and Oryza sativa.

Pathogenesis related (PR) proteins are low molecular weight family of proteins induced in plants under various biotic and abiotic stresses. They play an important role in plant-defense mechanism. PRs have wide range of functions, acting as hydrolases, peroxidases, chitinases, anti-fungal, protease inhibitors etc. In the present study, an attempt has been made to analyze promoter regions of PR1, PR2, PR5, PR9, PR10 and PR12 of Arabidopsis thaliana and Oryza sativa. Analysis of cis-element distribution revealed the functional multiplicity of PRs and provides insight into the gene regulation. CpG islands are observed only in rice PRs, which indicates that monocot genome contains more GC rich motifs than dicots. Tandem repeats were also observed in 5’ UTR of PR genes. Thus, the present study provides an understanding of regulation of PR genes and their versatile roles in plants.

Identification of Plagl1/Zac1 binding sites and target genes establishes its role in the regulation of extracellular matrix genes and the imprinted gene network.

PLAGL1/ZAC1 undergoes parental genomic imprinting, is paternally expressed, and is a member of the imprinted gene network (IGN). It encodes a zinc finger transcription factor with anti-proliferative activity and is a candidate tumor suppressor gene on 6q24 whose expression is frequently lost in various neoplasms. Conversely, gain of PLAGL1 function is responsible for transient neonatal diabetes mellitus, a rare genetic disease that results from defective pancreas development. In the present work, we showed that Plagl1 up-regulation was not associated with DNA damage-induced cell cycle arrest. It was rather associated with physiological cell cycle exit that occurred with contact inhibition, growth factor withdrawal, or cell differentiation. To gain insights into Plagl1 mechanism of action, we identified Plagl1 target genes by combining chromatin immunoprecipitation and genome-wide transcriptomics in transfected cell lines. Plagl1-elicited gene regulation correlated with multiple binding to the proximal promoter region through a GC-rich motif. Plagl1 target genes included numerous genes involved in signaling, cell adhesion, and extracellular matrix composition, including collagens. Plagl1 targets also included 22% of the 409 genes that make up the IGN. Altogether, this work identified Plagl1 as a transcription factor that coordinated the regulation of a subset of IGN genes and controlled extracellular matrix composition.

Short DNA sequence patterns accurately identify broadly active human enhancers

BackgroundEnhancers are DNA regulatory elements that influence gene expression. There is substantial diversity in enhancers’ activity patterns: some enhancers drive expression in a single cellular context, while others are active across many. Sequence characteristics, such as transcription factor (TF) binding motifs, influence the activity patterns of regulatory sequences; however, the regulatory logic through which specific sequences drive enhancer activity patterns is poorly understood. Recent analysis of Drosophila enhancers suggested that short dinucleotide repeat motifs (DRMs) are general enhancer sequence features that drive broad regulatory activity. However, it is not known whether the regulatory role of DRMs is conserved across species.ResultsWe performed a comprehensive analysis of the relationship between short DNA sequence patterns, including DRMs, and human enhancer activity in 38,538 enhancers across 411 different contexts. In a machine-learning framework, the occurrence patterns of short sequence motifs accurately predicted broadly active human enhancers. However, DRMs alone were weakly predictive of broad enhancer activity in humans and showed different enrichment patterns than in Drosophila. In general, GC-rich sequence motifs were significantly associated with broad enhancer activity, and consistent with this enrichment, broadly active human TFs recognize GC-rich motifs.ConclusionsOur results reveal the importance of specific sequence motifs in broadly active human enhancers, demonstrate the lack of evolutionary conservation of the role of DRMs, and provide a computational framework for investigating the logic of enhancer sequences.

Inhibition of cell proliferation through an ATP-responsive co-delivery system of doxorubicin and Bcl-2 siRNA.

Herein, DNA duplex was constructed through the hybridization of adenosine triphosphate (ATP)-responsive aptamer and its cDNA in which GC-rich motif could be used to load doxorubicin (DOX), and then, cationic polymer PEI25K was used as a carrier to simultaneously condense DOX-Duplex and Bcl-2 siRNA to prepare the ternary nanocomplex polyethylenimine (PEI)/DOX-Duplex/siRNA. The ATP concentration gradient between the cytosol and extracellular environment could achieve the stable loading of DOX in duplex and the rapid drug release in an ATP-responsive manner. Using human prostate tumor cell line PC-3 as a model, an obvious induction of cell proliferation could be detected with a cell viability of 53.3%, which was stronger than single cargo delivery, indicating the synergistic effect between these two components. The enhanced anti-proliferative effect of ternary nanocomplex could be attributed to the improved induction of cell apoptosis in a mitochondria-mediated pathway and cell-cycle arrest at the G2 phase. Overall, the ATP-responsive nanocarrier for co-delivering DOX and Bcl-2 siRNA has been demonstrated to be a smart delivery system with favorable anti-proliferative effect, especially for solving the multidrug resistance of tumors.

An ATP-Responsive Codelivery System of Doxorubicin and MiR-34a To Synergistically Inhibit Cell Proliferation and Migration.

Establishing stimulus-responsive nanosystems for the codelivery of anticancer drug and oligonucleotide is a promising strategy in cancer treatment owing to the combination of chemotherapy and gene therapy in a synergistic manner. Herein, an ATP aptamer and its cDNA sequence were first hybridized to produce the duplex, into which chemotherapeutic agent doxorubicin (DOX) interacted through the GC-rich motif of duplex, and PEI25K was then employed as a carrier to condense the DOX-loading duplex and miR-34a to construct the ternary nanocomplex PEI/DOX-Duplex/miR-34a. The nanocomplex exhibited a favorable drug release profile through the response to high concentration of ATP in the cytosol. The ATP-responsive delivery system was demonstrated to possess higher antiproliferative effect (cell viability of <40%) than the single cargo delivery, which could be attributed to the synergistic induction of cell apoptosis and cell cycle arrest from DOX and miR-34a. Furthermore, wound healing and Transwell assay elucidated the higher antimigration effect of ternary nanocomplex than DOX-Duplex or miR-34a delivery. Overall, the combinatorial delivery of DOX and miR-34a through an ATP-responsive manner could trigger the rapid release of cargoes in the cytosol and enhance the inhibition of cell proliferation and migration through the synergistic manner of these two components.

A GC-rich sequence feature in the 3' UTR directs UPF1-dependent mRNA decay in mammalian cells.

Up-frameshift protein 1 (UPF1) is an ATP-dependent RNA helicase that has essential roles in RNA surveillance and in post-transcriptional gene regulation by promoting the degradation of mRNAs. Previous studies revealed that UPF1 is associated with the 3′ untranslated region (UTR) of target mRNAs via as-yet-unknown sequence features. Herein, we aimed to identify characteristic sequence features of UPF1 targets. We identified 246 UPF1 targets by measuring RNA stabilization upon UPF1 depletion and by identifying mRNAs that associate with UPF1. By analyzing RNA footprint data of phosphorylated UPF1 and two CLIP-seq data of UPF1, we found that 3′ UTR but not 5′ UTRs or open reading frames of UPF1 targets have GC-rich motifs embedded in high GC-content regions. Reporter gene experiments revealed that GC-rich motifs in UPF1 targets were indispensable for UPF1-mediated mRNA decay. These findings highlight the important features of UPF1 target 3′ UTRs.

Gene Evolutionary Trajectories and GC Patterns Driven by Recombination in Zea mays.

Recombination occurring during meiosis is critical for creating genetic variation and plays an essential role in plant evolution. In addition to creating novel gene combinations, recombination can affect genome structure through altering GC patterns. In maize (Zea mays) and other grasses, another intriguing GC pattern exists. Maize genes show a bimodal GC content distribution that has been attributed to nucleotide bias in the third, or wobble, position of the codon. Recombination may be an underlying driving force given that recombination sites are often associated with high GC content. Here we explore the relationship between recombination and genomic GC patterns by comparing GC gene content at each of the three codon positions (GC1, GC2, and GC3, collectively termed GCx) to instances of a variable GC-rich motif that underlies double strand break (DSB) hotspots and to meiocyte-specific gene expression. Surprisingly, GCx bimodality in maize cannot be fully explained by the codon wobble hypothesis. High GCx genes show a strong overlap with the DSB hotspot motif, possibly providing a mechanism for the high evolutionary rates seen in these genes. On the other hand, genes that are turned on in meiosis (early prophase I) are biased against both high GCx genes and genes with the DSB hotspot motif, possibly allowing important meiotic genes to avoid DSBs. Our data suggests a strong link between the GC-rich motif underlying DSB hotspots and high GCx genes.

E2F coregulates an essential HSF developmental program that is distinct from the heat-shock response.

Heat-shock factor (HSF) is the master transcriptional regulator of the heat-shock response (HSR) and is essential for stress resilience. HSF is also required for metazoan development; however, its function and regulation in this process are poorly understood. Here, we characterize the genomic distribution and transcriptional activity of Caenorhabditis elegans HSF-1 during larval development and show that the developmental HSF-1 transcriptional program is distinct from the HSR. HSF-1 developmental activation requires binding of E2F/DP to a GC-rich motif that facilitates HSF-1 binding to a heat-shock element (HSE) that is degenerate from the consensus HSE sequence and adjacent to the E2F-binding site at promoters. In contrast, induction of the HSR is independent of these promoter elements or E2F/DP and instead requires a distinct set of tandem canonical HSEs. Together, E2F and HSF-1 directly regulate a gene network, including a specific subset of chaperones, to promote protein biogenesis and anabolic metabolism, which are essential in development.

Characterization and utilization of microsatellites in the Coffea canephora genome to assess genetic association between wild species in Kenya and cultivated coffee

Coffee is an important beverage crop in the world and has a significant contribution to Kenya’s economy. Here, we analyzed the genome-wide distribution of microsatellites in the Coffea canephora genome. A total of 159,041 SSRs were identified, with an overall density of 308 SSRs per Mb. Tetra-nucleotide repeats are the most abundant, accounting for 32 % of the total SSRs. AT-rich motifs are dominant across all SSR repeat units, while GC-rich motifs were generally rare. A set of 100 SSRs was selected to amplify 96 coffee accessions, including 10 wild accessions collected from Mt. Marsabit (Kenya). Of these SSRs, 33 % generated clear polymorphic bands among all tested accessions, with an average of 3.9 alleles per SSR locus. Wild coffee species from Mt. Marsabit showed a close genetic similarity with cultivated accessions in Kenya, suggesting that the wild species in Mt. Marsabit played an important role in the domestication of cultivated coffee in Kenya. Significantly low pairwise genetic divergence was observed between cultivated and wild accessions in Kenya, suggesting a relatively narrow level of genetic basis among coffee germplasm in Kenya. In addition, cultivated and wild coffee accessions in Kenya show a great divergence from those in other countries. Our results not only provide molecular tools for genetic studies in coffee but are also helpful for conservation and coffee breeding programs in Kenya.

Hepatic DGAT2 gene expression is regulated by the synergistic action of ChREBP and SP1 in HepG2 cells

ABSTRACTDiacylglycerol acyltransferase 2 (DGAT2) catalyzes the final step of triglyceride synthesis and plays a critical role in the development of fatty liver disease and insulin resistance. The expression of DGAT2 is mostly induced by dietary carbohydrate in the mammalian liver; however, the transcription factors that regulate DGAT2 expression have yet to be identified. In this study, we investigated the molecular mechanism underlying the glucose-induced transcriptional regulation of human DGAT2 in HepG2 cells. Human DGAT2 expression was increased by glucose in HepG2 cells. Transfection studies of the DGAT2 promoter-luciferase reporter construct in vitro and in silico analysis identified two glucose-responsive regions in the DGAT2 promoter. Each region contains one ChREBP/MLX (carbohydrate response element binding protein/Max-like protein) binding site (ChoRE, −539 to −551 and −6067 to −6083) and one specificity protein 1 (SP1) binding site (GC-rich motif, −556 to −572 and −6016 to −6032). Mutational analysis showed that both ChREBP/MLX and SP1 sites are required for the glucose-induced transcription of DGAT2. Gel shift assays and chromatin immunoprecipitation assays demonstrated that ChREBP and SP1 bind directly to ChoRE and the GC-rich motif, respectively. High glucose promoted the recruitment of ChREBP to ChoRE, whereas SP1 was recruited to the GC-rich motif even under low-glucose conditions. These data demonstrate that both ChREBP and SP1 are key factors to regulate the expression of human DGAT2 by glucose.

Development of genic SSR marker resource from RNA-Seq data in Dendrocalamus latiflorus

Bamboo is one of the most essential forest resources. Microsatellites have important role in genome, and are preferred over other molecular markers. For the first time, we exploited public transcriptome data of Dendrocalamus latiflorus for identification and characterization of microsatellite markers in the species. Frequency of SSRs was found to be 1 per 4 Kb. Tri-repeats were found to be the most abundant (47.7 %). Among them, GC rich repeat motifs namely, CCG/CGG (24.18 %) and AGG/CCT (22.29 %) were predominant. 12,028 potential genic SSR markers were identified, of which, 92 % were class II type. A set of 21 functionally relevant markers were successfully validated in 5 accessions of D. latiflorus. Twenty markers were found polymorphic, with number of alleles and polymorphism information content averaging 4.8 and 0.452, respectively. Interestingly, 100 % markers showed cross transferability across 14 related bamboo species, and were utilized for establishing phylogenetic relationships. Genic SSR markers developed herein would assist in genetic diversity, evolutionary and association mapping studies in bamboo.

Ribosome A and P sites revealed by length analysis of ribosome profiling data.

The high-throughput sequencing of nuclease-protected mRNA fragments bound to ribosomes, a technique known as ribosome profiling, quantifies the relative frequencies with which different regions of transcripts are translated. This technique has revealed novel translation initiation sites with unprecedented scope and has furthered investigations into the connections between codon biases and translation rates. Yet the location of the codon being decoded in ribosome footprints is still unknown, and has been complicated by the recent observation of footprints with non-canonical lengths. Here we show how taking into account the variations in ribosome footprint lengths can reveal the ribosome aminoacyl (A) and peptidyl (P) site locations. These location assignments are in agreement with the proposed mechanisms for various ribosome pauses and further enhance the resolution of the profiling data. We also show that GC-rich motifs at the 5′ ends of footprints are found in yeast, calling into question the anti-Shine-Dalgarno effect's role in ribosome pausing.

Zinc finger independent genome-wide binding of Sp2 potentiates recruitment of histone-fold protein Nf-y distinguishing it from Sp1 and Sp3.

Transcription factors are grouped into families based on sequence similarity within functional domains, particularly DNA-binding domains. The Specificity proteins Sp1, Sp2 and Sp3 are paradigmatic of closely related transcription factors. They share amino-terminal glutamine-rich regions and a conserved carboxy-terminal zinc finger domain that can bind to GC rich motifs in vitro. All three Sp proteins are ubiquitously expressed; yet they carry out unique functions in vivo raising the question of how specificity is achieved. Crucially, it is unknown whether they bind to distinct genomic sites and, if so, how binding site selection is accomplished. In this study, we have examined the genomic binding patterns of Sp1, Sp2 and Sp3 in mouse embryonic fibroblasts by ChIP-seq. Sp1 and Sp3 essentially occupy the same promoters and localize to GC boxes. The genomic binding pattern of Sp2 is different; Sp2 primarily localizes at CCAAT motifs. Consistently, re-expression of Sp2 and Sp3 mutants in corresponding knockout MEFs revealed strikingly different modes of genomic binding site selection. Most significantly, while the zinc fingers dictate genomic binding of Sp3, they are completely dispensable for binding of Sp2. Instead, the glutamine-rich amino-terminal region is sufficient for recruitment of Sp2 to its target promoters in vivo. We have identified the trimeric histone-fold CCAAT box binding transcription factor Nf-y as the major partner for Sp2-chromatin interaction. Nf-y is critical for recruitment of Sp2 to co-occupied regulatory elements. Equally, Sp2 potentiates binding of Nf-y to shared sites indicating the existence of an extensive Sp2-Nf-y interaction network. Our results unveil strikingly different recruitment mechanisms of Sp1/Sp2/Sp3 transcription factor members uncovering an unexpected layer of complexity in their binding to chromatin in vivo.