The ets-Related Transcription Factor GABP Directs Bidirectional Transcription

Approximately 10% of genes in the human genome are distributed such that their transcription start sites are located less than 1 kb apart on opposite strands. These divergent gene pairs have a single intergenic segment of DNA, which in some cases appears to share regulatory elements, but it is unclear whether these regions represent functional bidirectional promoters or two overlapping promoters. A recent study showed that divergent promoters are enriched for consensus binding sequences of a small group of transcription factors, including the ubiquitous ets-family transcription factor GA-binding protein (GABP). Here we show that GABP binds to more than 80% of divergent promoters in at least one cell type. Furthermore, we demonstrate that GABP binding is correlated and associated with bidirectional transcriptional activity in a luciferase transfection assay. In addition, we find that the addition of a strict consensus GABP site into a set of promoters that normally function in only one direction significantly increases activity in the opposite direction in 67% of cases. Our findings demonstrate that GABP regulates the majority of divergent promoters and suggest that bidirectional transcriptional activity is mediated through GABP binding and transactivation at both divergent and nondivergent promoters

A functional analysis of the CREB signaling pathway using HaloCHIP-chip and high throughput reporter assays

<p>Abstract</p> <p>Background</p> <p>Regulation of gene expression is essential for normal development and cellular growth. Transcriptional events are tightly controlled both spatially and temporally by specific DNA-protein interactions. In this study we finely map the genome-wide targets of the CREB protein across all known and predicted human promoters, and characterize the functional consequences of a subset of these binding events using high-throughput reporter assays. To measure CREB binding, we used HaloCHIP, an antibody-free alternative to the ChIP method that utilizes the HaloTag fusion protein, and also high-throughput promoter-luciferase reporter assays, which provide rapid and quantitative screening of promoters for transcriptional activation or repression in living cells.</p> <p>Results</p> <p>In analysis of CREB genome-wide binding events using a comprehensive DNA microarray of human promoters, we observe for the first time that CREB has a strong preference for binding at bidirectional promoters and unlike unidirectional promoters, these binding events often occur downstream of transcription start sites. Comparison between HaloCHIP-chip and ChIP-chip data reveal this to be true for both methodologies, indicating it is not a bias of the technology chosen. Transcriptional data obtained from promoter-luciferase reporter arrays also show an unprecedented, high level of activation of CREB-bound promoters in the presence of the co-activator protein TORC1.</p> <p>Conclusion</p> <p>These data suggest for the first time that TORC1 provides directional information when CREB is bound at bidirectional promoters and possible pausing of the CREB protein after initial transcriptional activation. Also, this combined approach demonstrates the ability to more broadly characterize CREB protein-DNA interactions wherein not only DNA binding sites are discovered, but also the potential of the promoter sequence to respond to CREB is evaluated.</p

Functional analysis of transcription factor binding sites in human promoters

BACKGROUND: The binding of transcription factors to specific locations in the genome is integral to the orchestration of transcriptional regulation in cells. To characterize transcription factor binding site function on a large scale, we predicted and mutagenized 455 binding sites in human promoters. We carried out functional tests on these sites in four different immortalized human cell lines using transient transfections with a luciferase reporter assay, primarily for the transcription factors CTCF, GABP, GATA2, E2F, STAT, and YY1. RESULTS: In each cell line, between 36% and 49% of binding sites made a functional contribution to the promoter activity; the overall rate for observing function in any of the cell lines was 70%. Transcription factor binding resulted in transcriptional repression in more than a third of functional sites. When compared with predicted binding sites whose function was not experimentally verified, the functional binding sites had higher conservation and were located closer to transcriptional start sites (TSSs). Among functional sites, repressive sites tended to be located further from TSSs than were activating sites. Our data provide significant insight into the functional characteristics of YY1 binding sites, most notably the detection of distinct activating and repressing classes of YY1 binding sites. Repressing sites were located closer to, and often overlapped with, translational start sites and presented a distinctive variation on the canonical YY1 binding motif. CONCLUSIONS: The genomic properties that we found to associate with functional TF binding sites on promoters -- conservation, TSS proximity, motifs and their variations -- point the way to improved accuracy in future TFBS predictions

Genome-Wide Discovery of Drug-Dependent Human Liver Regulatory Elements

Inter-individual variation in gene regulatory elements is hypothesized to play a causative role in adverse drug reactions and reduced drug activity. However, relatively little is known about the location and function of drug-dependent elements. To uncover drug-associated elements in a genome-wide manner, we performed RNA-seq and ChIP-seq using antibodies against the pregnane X receptor (PXR) and three active regulatory marks (p300, H3K4me1, H3K27ac) on primary human hepatocytes treated with rifampin or vehicle control. Rifampin and PXR were chosen since they are part of the CYP3A4 pathway, which is known to account for the metabolism of more than 50% of all prescribed drugs. We selected 227 proximal promoters for genes with rifampin-dependent expression or nearby PXR/p300 occupancy sites and assayed their ability to induce luciferase in rifampin-treated HepG2 cells, finding only 10 (4.4%) that exhibited drug-dependent activity. As this result suggested a role for distal enhancer modules, we searched more broadly to identify 1,297 genomic regions bearing a conditional PXR occupancy as well as all three active regulatory marks. These regions are enriched near genes that function in the metabolism of xenobiotics, specifically members of the cytochrome P450 family. We performed enhancer assays in rifampin-treated HepG2 cells for 42 of these sequences as well as 7 sequences that overlap linkage-disequilibrium blocks defined by lead SNPs from pharmacogenomic GWAS studies, revealing 15/42 and 4/7 to be functional enhancers, respectively. A common African haplotype in one of these enhancers in the GSTA locus was found to exhibit potential rifampin hypersensitivity. Combined, our results further suggest that enhancers are the predominant targets of rifampin-induced PXR activation, provide a genome-wide catalog of PXR targets and serve as a model for the identification of drug-responsive regulatory elements

Identification, Replication, and Functional Fine-Mapping of Expression Quantitative Trait Loci in Primary Human Liver Tissue

The discovery of expression quantitative trait loci (“eQTLs”) can help to unravel genetic contributions to complex traits. We identified genetic determinants of human liver gene expression variation using two independent collections of primary tissue profiled with Agilent (n = 206) and Illumina (n = 60) expression arrays and Illumina SNP genotyping (550K), and we also incorporated data from a published study (n = 266). We found that ∼30% of SNP-expression correlations in one study failed to replicate in either of the others, even at thresholds yielding high reproducibility in simulations, and we quantified numerous factors affecting reproducibility. Our data suggest that drug exposure, clinical descriptors, and unknown factors associated with tissue ascertainment and analysis have substantial effects on gene expression and that controlling for hidden confounding variables significantly increases replication rate. Furthermore, we found that reproducible eQTL SNPs were heavily enriched near gene starts and ends, and subsequently resequenced the promoters and 3′UTRs for 14 genes and tested the identified haplotypes using luciferase assays. For three genes, significant haplotype-specific in vitro functional differences correlated directly with expression levels, suggesting that many bona fide eQTLs result from functional variants that can be mechanistically isolated in a high-throughput fashion. Finally, given our study design, we were able to discover and validate hundreds of liver eQTLs. Many of these relate directly to complex traits for which liver-specific analyses are likely to be relevant, and we identified dozens of potential connections with disease-associated loci. These included previously characterized eQTL contributors to diabetes, drug response, and lipid levels, and they suggest novel candidates such as a role for NOD2 expression in leprosy risk and C2orf43 in prostate cancer. In general, the work presented here will be valuable for future efforts to precisely identify and functionally characterize genetic contributions to a variety of complex traits

Multispecific Antibody Development Platform Based on Human Heavy Chain Antibodies

Heavy chain-only antibodies (HCAbs) do not associate with light chains and their VH regions are functional as single domains, forming the smallest active antibody fragment. These VH regions are ideal building blocks for a variety of antibody-based biologics because they tolerate fusion to other molecules and may also be attached in series to construct multispecific antibodies without the need for protein engineering to ensure proper heavy and light chain pairing. Production of human HCAbs has been impeded by the fact that natural human VH regions require light chain association and display poor biophysical characteristics when expressed in the absence of light chains. Here, we present an innovative platform for the rapid development of diverse sets of human HCAbs that have been selected in vivo. Our unique approach combines antibody repertoire analysis with immunization of transgenic rats, called UniRats, that produce chimeric HCAbs with fully human VH domains in response to an antigen challenge. UniRats express HCAbs from large transgenic loci representing the entire productive human heavy chain V(D)J repertoire, mount robust immune responses to a wide array of antigens, exhibit diverse V gene usage and generate large panels of stable, high affinity, antigen-specific molecules

Serum response factor binding sites differ in three human cell types

The serum response factor (SRF) is essential for embryonic development and maintenance of muscle cells and neurons. The mechanism by which this factor controls these divergent pathways is unclear. Here we present a genome-wide view of occupancy of SRF at its binding sites with a focus on those that vary with cell type. We used chromatin immunoprecipitation (ChIP) in combination with human promoter microarrays to identify 216 putative SRF binding sites in the human genome. We performed independent quantitative PCR validation at over half of these sites that resulted in 146 sites we assert to be true binding sites at over 90% confidence. Nearly half of the sites are bound by SRF in only one of the three cell types we tested, providing strong evidence for the diverse roles for SRF in different cell types. We also explore possible mechanisms controlling differential binding of SRF in these cell types by assaying cofactor binding, DNA methylation, histone methylation, and histone acetylation at a subset of sites bound preferentially in smooth muscle cells. Although we did not see a strong correlation between SRF binding and epigenetics modifications, at these sites, we propose that SRF cofactors may play an important role in determining cell-dependent SRF binding sites. ELK4 (previously known as SAP-1 [SRF-associated protein-1]) is ubiquitously expressed. Therefore, we expected it to occupy sites where SRF binding is common in all cell types. Indeed, 90% of SRF sites also bound by ELK4 were common to all three cell types. Together, our data provide a more complete understanding of the regulatory network controlled by SRF

Transcriptional regulation and binding of heat shock factor 1 and heat shock factor 2 to 32 human heat shock genes during thermal stress and differentiation

Transcription of mammalian heat shock genes can be regulated by heat shock factors (HSF) 1 and 2. Although it has been shown previously that these factors respond to distinct stimuli, a broad analysis of the induction and function of these factors in living cells has not been performed. In our study, we assayed binding of human HSF1 and HSF2 at the promoters of 32 genes identified through LocusLink as heat shock genes in response to elevated temperature and hemin-induced differentiation in human K562 erythroleukemic cells using the chromatin immunoprecipitation technique. We also measured the induced expression of these genes under these 2 conditions. We found that 17 of the 32 genes were transcriptionally induced during heat shock, and HSF1 binding was detected at 15 of the 17 promoters. Nearly all the genes induced by heat shock were also induced to a lesser degree during hemin treatment. However, some genes were induced significantly more during hemin treatment than during heat shock. A new finding is that HSF1 and HSF2 bind to the same targets, but HSF1 binding is activated more by heat than by hemin treatment, and HSF2 binding is only activated by hemin treatment and not by heat. This technology also identified previously unknown HSF1 binding sites near genes that were previously shown to be heat inducible that may contribute to gene-specific regulation

Comprehensive analysis of transcriptional promoter structure and function in 1% of the human genome

Transcriptional promoters comprise one of many classes of eukaryotic transcriptional regulatory elements. Identification and characterization of these elements are vital to understanding the complex network of human gene regulation. Using full-length cDNA sequences to identify transcription start sites (TSS), we predicted more than 900 putative human transcriptional promoters in the ENCODE regions, representing a comprehensive sampling of promoters in 1% of the genome. We identified 387 fragments that function as promoters in at least one of 16 cell lines by measuring promoter activity in high-throughput transient transfection reporter assays. These positive functional results demonstrate widespread use of alternative promoters. We show a strong correlation between promoter activity and the corresponding endogenous RNA transcript levels, providing the first experimental quantitative estimate of promoter contribution to gene regulation. Finally, we identified functional regions within a randomly selected subset of 45 promoters using deletion analyses. These experiments showed that, on average, the sequence -300 to -50 bp of the TSS positively contributes to core promoter activity. Interestingly, putative negative elements were identified -1000 to -500 bp upstream of the TSS for 55% of genes tested. These data provide the largest and most comprehensive view of promoter function in the human genome

Identification and Functional Analysis of Human Transcriptional Promoters

Genomic and full-length cDNA sequences provide opportunities for understanding human gene structure and transcriptional regulatory elements. The simplest regulatory elements to identify are promoters, as their positions are dictated by the location of transcription start sites. We aligned full-length cDNA clones from the Mammalian Gene Collection to the human genome rough draft sequence to estimate the start sites of more than 10,000 human transcripts. We selected genomic sequence just upstream from the 5′ end of these cDNA sequences and designated these as putative promoters. We assayed the functions of 152 of these DNA fragments, chosen at random from the entire set, in a luciferase-based transfection assay in four human cultured cell types. Ninety-one percent of these DNA fragments showed significant transcriptional activity in at least one of the cell lines, whereas 89% showed activity in at least two of the lines. We analyzed the distributions of strengths of these promoter fragments in the different cell types and identified likely alternative promoters in a large fraction of the genes. These data indicate that this approach is an effective method for predicting human promoters and provide the first set of functional data collected in parallel for a large set of human promoters. [Supplemental material is available online at www.genome.org and http://www-shgc.stanford.edu/myerslab/.