Search

Patterns of transcriptional activity are encoded in our genome through regulatory elements such as promoters or enhancers that, paradoxically, contain similar assortments of sequence-specific transcription factor (TF) binding sites1-3. Knowledge of how these sequence motifs encode multiple, often overlapping, gene expression programs is central to understanding gene regulation and how mutations in non-coding DNA manifest in disease4,5. Here, by studying gene regulation from the perspective of individual transcription start sites (TSSs), using natural genetic variation, perturbation of endogenous TF protein levels and massively parallel analysis of natural and synthetic regulatory elements, we show that the effect of TF binding on transcription initiation is position dependent. Analysing TF-binding-site occurrences relative to the TSS, we identified several motifs with highly preferential positioning. We show that these patterns are a combination of a TF's distinct functional profiles-many TFs, including canonical activators such as NRF1, NFY and Sp1, activate or repress transcription initiation depending on their precise position relative to the TSS. As such, TFs and their spacing collectively guide the site and frequency of transcription initiation. More broadly, these findings reveal how similar assortments of TF binding sites can generate distinct gene regulatory outcomes depending on their spatial configuration and how DNA sequence polymorphisms may contribute to transcription variation and disease and underscore a critical role for TSS data in decoding the regulatory information of our genome.

The prevalent one-dimensional alignment of genomic signals to a reference landmark is a cornerstone of current methods to study transcription and its DNA-dependent processes but it is prone to mask potential relations among multiple DNA elements. We developed a systematic approach to align genomic signals to multiple locations simultaneously by expanding the dimensionality of the genomic-coordinate space. We analyzed transcription in human and uncovered a complex dependence on the relative position of neighboring transcription start sites (TSSs) that is consistently conserved among cell types. The dependence ranges from enhancement to suppression of transcription depending on the relative distances to the TSSs, their intragenic position, and the transcriptional activity of the gene. Our results reveal a conserved hierarchy of alternative TSS usage within a previously unrecognized level of genomic organization and provide a general methodology to analyze complex functional relationships among multiple types of DNA elements.

Accurate transcription start site (TSS) annotations are essential for understanding transcriptional regulation and its role in human disease. Gene collections such as GENCODE contain annotations for tens of thousands of TSSs, but not all of these annotations are experimentally validated nor do they contain information on cell type-specific usage. Therefore, we sought to generate a collection of experimentally validated TSSs by integrating RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression (RAMPAGE) data from 115 cell and tissue types, which resulted in a collection of approximately 50 thousand representative RAMPAGE peaks. These peaks are primarily proximal to GENCODE-annotated TSSs and are concordant with other transcription assays. Because RAMPAGE uses paired-end reads, we were then able to connect peaks to transcripts by analyzing the genomic positions of the 3' ends of read mates. Using this paired-end information, we classified the vast majority (37 thousand) of our RAMPAGE peaks as verified TSSs, updating TSS annotations for 20% of GENCODE genes. We also found that these updated TSS annotations are supported by epigenomic and other transcriptomic data sets. To show the utility of this RAMPAGE rPeak collection, we intersected it with the NHGRI/EBI genome-wide association study (GWAS) catalog and identified new candidate GWAS genes. Overall, our work shows the importance of integrating experimental data to further refine TSS annotations and provides a valuable resource for the biological community.

The general transcription factor TFIIH contains three ATP-dependent catalytic activities. TFIIH functions in nucleotide excision repair primarily as a DNA helicase and in Pol II transcription initiation as a dsDNA translocase and protein kinase. During initiation, the XPB/Ssl2 subunit of TFIIH couples ATP hydrolysis to dsDNA translocation facilitating promoter opening and the kinase module phosphorylates Pol II to facilitate the transition to elongation. These functions are conserved between metazoans and yeast; however, yeast TFIIH also drives transcription start-site scanning in which Pol II scans downstream DNA to locate productive start-sites. The ten-subunit holo-TFIIH from S. cerevisiae has a processive dsDNA translocase activity required for scanning and a structural role in scanning has been ascribed to the three-subunit TFIIH kinase module. Here, we assess the dsDNA translocase activity of ten-subunit holo- and core-TFIIH complexes (i.e. seven subunits, lacking the kinase module) from both S. cerevisiae and H. sapiens. We find that neither holo nor core human TFIIH exhibit processive translocation, consistent with the lack of start-site scanning in humans. Furthermore, in contrast to holo-TFIIH, the S. cerevisiae core-TFIIH also lacks processive translocation and its dsDNA-stimulated ATPase activity was reduced ~5-fold to a level comparable to the human complexes, potentially explaining the reported upstream shift in start-site observed in vitro in the absence of the S. cerevisiae kinase module. These results suggest that neither human nor S. cerevisiae core-TFIIH can translocate efficiently, and that the S. cerevisiae kinase module functions as a processivity factor to allow for robust transcription start-site scanning.

Using the Cap Analysis of Gene Expression (CAGE) technology, the FANTOM5 consortium provided one of the most comprehensive maps of transcription start sites (TSSs) in several species. Strikingly, ~72% of them could not be assigned to a specific gene and initiate at unconventional regions, outside promoters or enhancers. Here, we probe these unassigned TSSs and show that, in all species studied, a significant fraction of CAGE peaks initiate at microsatellites, also called short tandem repeats (STRs). To confirm this transcription, we develop Cap Trap RNA-seq, a technology which combines cap trapping and long read MinION sequencing. We train sequence-based deep learning models able to predict CAGE signal at STRs with high accuracy. These models unveil the importance of STR surrounding sequences not only to distinguish STR classes, but also to predict the level of transcription initiation. Importantly, genetic variants linked to human diseases are preferentially found at STRs with high transcription initiation level, supporting the biological and clinical relevance of transcription initiation at STRs. Together, our results extend the repertoire of non-coding transcription associated with DNA tandem repeats and complexify STR polymorphism.

Characterizing transcription start sites is essential for understanding the regulatory mechanisms that control gene expression. Recently, a new bovine genome assembly (ARS-UCD1.2) with high continuity, accuracy, and completeness was released; however, the functional annotation of the bovine genome lacks precise transcription start sites and contains a low number of transcripts in comparison to human and mouse. By using the RAMPAGE approach, this study identified transcription start sites at high resolution in a large collection of bovine tissues. We found several known and novel transcription start sites attributed to promoters of protein-coding and lncRNA genes that were validated through experimental and in silico evidence. With these findings, the annotation of transcription start sites in cattle reached a level comparable to the mouse and human genome annotations. In addition, we identified and characterized transcription start sites for antisense transcripts derived from bidirectional promoters, potential lncRNAs, mRNAs, and pre-miRNAs. We also analyzed the quantitative aspects of RAMPAGE to produce a promoter activity atlas, reaching highly reproducible results comparable to traditional RNA-seq. Coexpression networks revealed considerable use of tissue-specific promoters, especially between brain and testicle, which expressed several genes in common from alternate loci. Furthermore, regions surrounding coexpressed modules were enriched in binding factor motifs representative of each tissue. The comprehensive annotation of promoters in such a large collection of tissues will substantially contribute to our understanding of gene expression in cattle and other mammalian species, shortening the gap between genotypes and phenotypes.

Human genome stability requires efficient repair of oxidized bases, which is initiated via damage recognition and excision by NEIL1 and other base excision repair (BER) pathway DNA glycosylases (DGs). However, the biological mechanisms underlying detection of damaged bases among the million-fold excess of undamaged bases remain enigmatic. Indeed, mutation rates vary greatly within individual genomes, and lesion recognition by purified DGs in the chromatin context is inefficient. Employing super-resolution microscopy and co-immunoprecipitation assays, we find that acetylated NEIL1 (AcNEIL1), but not its non-acetylated form, is predominantly localized in the nucleus in association with epigenetic marks of uncondensed chromatin. Furthermore, chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) revealed non-random AcNEIL1 binding near transcription start sites of weakly transcribed genes and along highly transcribed chromatin domains. Bioinformatic analyses revealed a striking correspondence between AcNEIL1 occupancy along the genome and mutation rates, with AcNEIL1-occupied sites exhibiting fewer mutations compared to AcNEIL1-free domains, both in cancer genomes and in population variation. Intriguingly, from the evolutionarily conserved unstructured domain that targets NEIL1 to open chromatin, its damage surveillance of highly oxidation-susceptible sites to preserve essential gene function and to limit instability and cancer likely originated ∼500 million years ago during the buildup of free atmospheric oxygen.

Sustained, drug-free control of HIV-1 replication is naturally achieved in less than 0.5% of infected individuals (here termed 'elite controllers'), despite the presence of a replication-competent viral reservoir1. Inducing such an ability to spontaneously maintain undetectable plasma viraemia is a major objective of HIV-1 cure research, but the characteristics of proviral reservoirs in elite controllers remain to be determined. Here, using next-generation sequencing of near-full-length single HIV-1 genomes and corresponding chromosomal integration sites, we show that the proviral reservoirs of elite controllers frequently consist of oligoclonal to near-monoclonal clusters of intact proviral sequences. In contrast to individuals treated with long-term antiretroviral therapy, intact proviral sequences from elite controllers were integrated at highly distinct sites in the human genome and were preferentially located in centromeric satellite DNA or in Krüppel-associated box domain-containing zinc finger genes on chromosome 19, both of which are associated with heterochromatin features. Moreover, the integration sites of intact proviral sequences from elite controllers showed an increased distance to transcriptional start sites and accessible chromatin of the host genome and were enriched in repressive chromatin marks. These data suggest that a distinct configuration of the proviral reservoir represents a structural correlate of natural viral control, and that the quality, rather than the quantity, of viral reservoirs can be an important distinguishing feature for a functional cure of HIV-1 infection. Moreover, in one elite controller, we were unable to detect intact proviral sequences despite analysing more than 1.5 billion peripheral blood mononuclear cells, which raises the possibility that a sterilizing cure of HIV-1 infection, which has previously been observed only following allogeneic haematopoietic stem cell transplantation2,3, may be feasible in rare instances.

G-quadruplexes are nucleic acid secondary structures formed by a stack of square planar G-quartets. G-quadruplexes were implicated in many biological functions including telomere maintenance, replication, transcription, and translation, in many species including humans and plants. For wheat, however, though it is one of the world's most important staple food, no G-quadruplex studies have been reported to date. Here, we computationally identify putative G4 structures (G4s) in wheat genome for the first time and compare its distribution across the genome against five other genomes (human, maize, Arabidopsis, rice, and sorghum). We identified close to 1 million G4 motifs with a density of 76 G4s/Mb across the whole genome and 93 G4s/Mb over genic regions. Remarkably, G4s were enriched around three regions, two located on the antisense and one on the sense strand at the following positions: 1) the transcription start site (TSS) (antisense), 2) the first coding domain sequence (CDS) (antisense), and 3) the start codon (sense). Functional enrichment analysis revealed that the gene models containing G4 motifs within these peaks were associated with specific gene ontology (GO) terms, such as developmental process, localization, and cellular component organization or biogenesis. We investigated genes encoding MADS-box transcription factors and showed examples of G4 motifs within critical regulatory regions in the VRN-1 genes in wheat. Furthermore, comparison with other plants showed that monocots share a similar distribution of G4s, but Arabidopsis shows a unique G4 distribution. Our study shows for the first time the prevalence and possible functional roles of G4s in wheat.

Metazoan development requires the robust proliferation of progenitor cells, the identities of which are established by tightly controlled transcriptional networks1. As gene expression is globally inhibited during mitosis, the transcriptional programs that define cell identity must be restarted in each cell cycle2-5 but how this is accomplished is poorly understood. Here we identify a ubiquitin-dependent mechanism that integrates gene expression with cell division to preserve cell identity. We found that WDR5 and TBP, which bind active interphase promoters6,7, recruit the anaphase-promoting complex (APC/C) to specific transcription start sites during mitosis. This allows APC/C to decorate histones with ubiquitin chains branched at Lys11 and Lys48 (K11/K48-branched ubiquitin chains) that recruit p97 (also known as VCP) and the proteasome, which ensures the rapid expression of pluripotency genes in the next cell cycle. Mitotic exit and the re-initiation of transcription are thus controlled by a single regulator (APC/C), which provides a robust mechanism for maintaining cell identity throughout cell division.

MeCP2 plays a multifaceted role in gene expression regulation and chromatin organization. Interaction between MeCP2 and methylated DNA in the regulation of gene expression is well established. However, the widespread distribution of MeCP2 suggests it has additional interactions with chromatin. Here we demonstrate, by both biochemical and genomic analyses, that MeCP2 directly interacts with nucleosomes and its genomic distribution correlates with that of H3K27me3. In particular, the methyl-CpG-binding domain of MeCP2 shows preferential interactions with H3K27me3. We further observe that the impact of MeCP2 on transcriptional changes correlates with histone post-translational modification patterns. Our findings indicate that MeCP2 interacts with genomic loci via binding to DNA as well as histones, and that interaction between MeCP2 and histone proteins plays a key role in gene expression regulation.

Although promoter-associated CpG islands have been established as targets of DNA methylation changes in cancer, previous studies suggest that epigenetic dysregulation outside the promoter region may be more closely associated with transcriptional changes. Here we examine DNA methylation, chromatin marks, and transcriptional alterations to define the relationship between transcriptional modulation and spatial changes in chromatin structure. Using human papillomavirus-related oropharyngeal carcinoma as a model, we show aberrant enrichment of repressive H3K9me3 at the transcriptional start site (TSS) with methylation-associated, tumor-specific gene silencing. Further analysis identifies a hypermethylated subtype which shows a functional convergence on MYC targets and association with CREBBP/EP300 mutation. The tumor-specific shift to transcriptional repression associated with DNA methylation at TSSs was confirmed in multiple tumor types. Our data may show a common underlying epigenetic dysregulation in cancer associated with broad enrichment of repressive chromatin marks and aberrant DNA hypermethylation at TSSs in combination with MYC network activation.

The potential impact of structural variants includes not only the duplication or deletion of coding sequences, but also the perturbation of noncoding DNA regulatory elements and structural chromatin features, including topological domains (TADs). Structural variants disrupting TAD boundaries have been implicated both in cancer and developmental disease; this likely occurs via "enhancer hijacking," whereby removal of the TAD boundary exposes enhancers to new target transcription start sites (TSSs). With this functional role, we hypothesized that boundaries would display evidence for negative selection. Here we demonstrate that the chromatin landscape constrains structural variation both within healthy humans and across primate evolution. In contrast, in patients with developmental delay, variants occur remarkably uniformly across genomic features, suggesting a potentially broad role for enhancer hijacking in human disease.

CCCTC-binding factor (CTCF) plays a key role in the formation of topologically associating domains (TADs) and loops in interphase. During mitosis TADs are absent, but how TAD formation is dynamically controlled during the cell cycle is not known. Several contradicting observations have been made regarding CTCF binding to mitotic chromatin using both genomics- and microscopy-based techniques. Here, we have used four different assays to address this debate. First, using 5C, we confirmed that TADs and CTCF loops are readily detected in interphase, but absent during prometaphase. Second, ATAC-seq analysis showed that CTCF sites display greatly reduced accessibility and lose the CTCF footprint in prometaphase, suggesting loss of CTCF binding and rearrangement of the nucleosomal array around the binding motif. In contrast, transcription start sites remain accessible in prometaphase, although adjacent nucleosomes can also become repositioned and occupy at least a subset of start sites during mitosis. Third, loss of site-specific CTCF binding was directly demonstrated using CUT&RUN. Histone modifications and histone variants are maintained in mitosis, suggesting a role in bookmarking of active CTCF sites. Finally, live-cell imaging, fluorescence recovery after photobleaching, and single molecule tracking showed that almost all CTCF chromatin binding is lost in prometaphase. Combined, our results demonstrate loss of CTCF binding to CTCF sites during prometaphase and rearrangement of the chromatin landscape around CTCF motifs. This, combined with loss of cohesin, would contribute to the observed loss of TADs and CTCF loops during mitosis and reveals that CTCF sites, key architectural cis-elements, display cell cycle stage-dependent dynamics in factor binding and nucleosome positioning.

Although recent studies provide evidence for a common genetic basis between complex traits and Mendelian disorders, a thorough quantification of their overlap in a phenotype-specific manner remains elusive. Here, we have quantified the overlap of genes identified through large-scale genome-wide association studies (GWASs) for 62 complex traits and diseases with genes containing mutations known to cause 20 broad categories of Mendelian disorders. We identified a significant enrichment of genes linked to phenotypically matched Mendelian disorders in GWAS gene sets; of the total 1,240 comparisons, a higher proportion of phenotypically matched or related pairs (n = 50 of 92 [54%]) than phenotypically unmatched pairs (n = 27 of 1,148 [2%]) demonstrated significant overlap, confirming a phenotype-specific enrichment pattern. Further, we observed elevated GWAS effect sizes near genes linked to phenotypically matched Mendelian disorders. Finally, we report examples of GWAS variants localized at the transcription start site or physically interacting with the promoters of genes linked to phenotypically matched Mendelian disorders. Our results are consistent with the hypothesis that genes that are disrupted in Mendelian disorders are dysregulated by non-coding variants in complex traits and demonstrate how leveraging findings from related Mendelian disorders and functional genomic datasets can prioritize genes that are putatively dysregulated by local and distal non-coding GWAS variants.

RNA-sequencing (RNA-seq) is a powerful technique to investigate and quantify entire transcriptomes. Recent advances in the field have made it possible to explore the transcriptomes of single cells. However, most widely used RNA-seq protocols fail to provide crucial information regarding transcription start sites. Here we present a protocol, Tn5Prime, that takes advantage of the Tn5 transposase-based Smart-seq2 protocol to create RNA-seq libraries that capture the 5' end of transcripts. The Tn5Prime method dramatically streamlines the 5' capture process and is both cost effective and reliable. By applying Tn5Prime to bulk RNA and single cell samples, we were able to define transcription start sites as well as quantify transcriptomes at high accuracy and reproducibility. Additionally, similar to 3' end-based high-throughput methods like Drop-seq and 10× Genomics Chromium, the 5' capture Tn5Prime method allows the introduction of cellular identifiers during reverse transcription, simplifying the analysis of large numbers of single cells. In contrast to 3' end-based methods, Tn5Prime also enables the assembly of the variable 5' ends of the antibody sequences present in single B-cell data. Therefore, Tn5Prime presents a robust tool for both basic and applied research into the adaptive immune system and beyond.

DNA sequence signals in the core promoter, such as the initiator (Inr), direct transcription initiation by RNA polymerase II. Here we show that the human Inr has the consensus of BBCA+1BW at focused promoters in which transcription initiates at a single site or a narrow cluster of sites. The analysis of 7678 focused transcription start sites revealed 40% with a perfect match to the Inr and 16% with a single mismatch outside of the CA+1 core. TATA-like sequences are underrepresented in Inr promoters. This consensus is a key component of the DNA sequence rules that specify transcription initiation in humans.

DNA sequence signals in the core promoter, such as the initiator (Inr), direct transcription initiation by RNA polymerase II. Here we show that the human Inr has the consensus of BBCA+1BW at focused promoters in which transcription initiates at a single site or a narrow cluster of sites. The analysis of 7678 focused transcription start sites revealed 40% with a perfect match to the Inr and 16% with a single mismatch outside of the CA+1 core. TATA-like sequences are underrepresented in Inr promoters. This consensus is a key component of the DNA sequence rules that specify transcription initiation in humans.

The Cancer Genome Atlas (TCGA) network study of 12 cancer types (PanCancer 12) revealed frequent mutation of TP53, and amplification and expression of related TP63 isoform ΔNp63 in squamous cancers. Further, aberrant expression of inflammatory genes and TP53/p63/p73 targets were detected in the PanCancer 12 project, reminiscent of gene programs comodulated by cREL/ΔNp63/TAp73 transcription factors we uncovered in head and neck squamous cell carcinomas (HNSCCs). However, how inflammatory gene signatures and cREL/p63/p73 targets are comodulated genome wide is unclear. Here, we examined how the inflammatory factor tumor necrosis factor-α (TNF-α) broadly modulates redistribution of cREL with ΔNp63α/TAp73 complexes and signatures genome wide in the HNSCC model UM-SCC46 using chromatin immunoprecipitation sequencing (ChIP-seq). TNF-α enhanced genome-wide co-occupancy of cREL with ΔNp63α on TP53/p63 sites, while unexpectedly promoting redistribution of TAp73 from TP53 to activator protein-1 (AP-1) sites. cREL, ΔNp63α and TAp73 binding and oligomerization on NF-κB-, TP53- or AP-1-specific sequences were independently validated by ChIP-qPCR (quantitative PCR), oligonucleotide-binding assays and analytical ultracentrifugation. Function of the binding activity was confirmed using TP53-, AP-1- and NF-κB-specific REs or p21, SERPINE1 and IL-6 promoter luciferase reporter activities. Concurrently, TNF-α regulated a broad gene network with cobinding activities for cREL, ΔNp63α and TAp73 observed upon array profiling and reverse transcription-PCR. Overlapping target gene signatures were observed in squamous cancer subsets and in inflamed skin of transgenic mice overexpressing ΔNp63α. Furthermore, multiple target genes identified in this study were linked to TP63 and TP73 activity and increased gene expression in large squamous cancer samples from PanCancer 12 TCGA by CircleMap. PARADIGM inferred pathway analysis revealed the network connection of TP63 and NF-κB complexes through an AP-1 hub, further supporting our findings. Thus, inflammatory cytokine TNF-α mediates genome-wide redistribution of the cREL/p63/p73, and AP-1 interactome, to diminish TAp73 tumor suppressor function and reciprocally activate NF-κB and AP-1 gene programs implicated in malignancy.

Nuclear pore complexes (NPCs) emerged as nuclear transport channels in eukaryotic cells ∼1.5 billion years ago. While the primary role of NPCs is to regulate nucleo-cytoplasmic transport, recent research suggests that certain NPC proteins have additionally acquired the role of affecting gene expression at the nuclear periphery and in the nucleoplasm in metazoans. Here we identify a widely expressed variant of the transmembrane nucleoporin (Nup) Pom121 (named sPom121, for "soluble Pom121") that arose by genomic rearrangement before the divergence of hominoids. sPom121 lacks the nuclear membrane-anchoring domain and thus does not localize to the NPC. Instead, sPom121 colocalizes and interacts with nucleoplasmic Nup98, a previously identified transcriptional regulator, at gene promoters to control transcription of its target genes in human cells. Interestingly, sPom121 transcripts appear independently in several mammalian species, suggesting convergent innovation of Nup-mediated transcription regulation during mammalian evolution. Our findings implicate alternate transcription initiation as a mechanism to increase the functional diversity of NPC components.

We report that the mammalian 5-methylcytosine (5mC) oxidase Tet3 exists as three major isoforms and characterized the full-length isoform containing an N-terminal CXXC domain (Tet3FL). This CXXC domain binds to unmethylated CpGs, but, unexpectedly, its highest affinity is toward 5-carboxylcytosine (5caC). We determined the crystal structure of the CXXC domain-5caC-DNA complex, revealing the structural basis of the binding specificity of this domain as a reader of CcaCG sequences. Mapping of Tet3FL in neuronal cells shows that Tet3FL is localized precisely at the transcription start sites (TSSs) of genes involved in lysosome function, mRNA processing, and key genes of the base excision repair pathway. Therefore, Tet3FL may function as a regulator of 5caC removal by base excision repair. Active removal of accumulating 5mC from the TSSs of genes coding for lysosomal proteins by Tet3FL in postmitotic neurons of the brain may be important for preventing neurodegenerative diseases.

Bromodomain and extra terminal domain (BET) proteins are important epigenetic regulators facilitating the transcription of genes in chromatin areas linked to acetylated histones. JQ1, a BET protein inhibitor, has antiproliferative activity against many cancers, mainly through inhibition of c-MYC and upregulation of p21. In this research, we investigated the use of JQ1 for human osteosarcoma (OS) treatment. JQ1 significantly inhibited the proliferation and survival of OS cells inducing G1 cell cycle arrest, premature senescence, but little effect on apoptosis. Interestingly, c-MYC protein levels in JQ1-treated cells remained unchanged, whereas the upregulation of p21 protein was still observable. Although effective in vitro, JQ1 alone failed to reduce the size of the MNNG/HOS xenografts in immunocompromised mice. To overcome the resistance of OS cells to JQ1 treatment, we combined JQ1 with rapamycin, an mammalian target of rapamycin (mTOR) inhibitor. JQ1 and rapamycin synergistically inhibited the growth and survival of OS cells in vitro and in vivo. We also identified that RUNX2 is a direct target of bromodomain-containing protein 4 (BRD4) inhibition by JQ1 in OS cells. Chromatin immunoprecipitation (ChIP) showed that enrichment of BRD4 protein around RUNX2 transcription start sites diminished with JQ1 treatment in MNNG/HOS cells. Overexpression of RUNX2 protected JQ1-sensitive OS cells from the effect of JQ1, and siRNA-mediated inhibition of RUNX2 sensitized the same cells to JQ1. In conclusion, our findings suggest that JQ1, in combination with rapamycin, is an effective chemotherapeutic option for OS treatment. We also show that inhibition of RUNX2 expression by JQ1 partly explains the antiproliferative activity of JQ1 in OS cells.

Divergent transcription, in which reverse-oriented transcripts occur upstream of eukaryotic promoters in regions devoid of annotated genes, has been suggested to be a general property of active promoters. Here we show that the human basal RNA polymerase II transcriptional machinery and core promoter are inherently unidirectional and that reverse-oriented transcripts originate from their own cognate reverse-directed core promoters. In vitro transcription analysis and mapping of nascent transcripts in HeLa cells revealed that sequences at reverse start sites are similar to those of their forward counterparts. The use of DNase I accessibility to define proximal promoter borders revealed that about half of promoters are unidirectional and that unidirectional promoters are depleted at their upstream edges of reverse core promoter sequences and their associated chromatin features. Divergent transcription is thus not an inherent property of the transcription process but rather the consequence of the presence of both forward- and reverse-directed core promoters.

Accurate gene model annotation of reference genomes is critical for making them useful. The modENCODE project has improved the D. melanogaster genome annotation by using deep and diverse high-throughput data. Since transcriptional activity that has been evolutionarily conserved is likely to have an advantageous function, we have performed large-scale interspecific comparisons to increase confidence in predicted annotations. To support comparative genomics, we filled in divergence gaps in the Drosophila phylogeny by generating draft genomes for eight new species. For comparative transcriptome analysis, we generated mRNA expression profiles on 81 samples from multiple tissues and developmental stages of 15 Drosophila species, and we performed cap analysis of gene expression in D. melanogaster and D. pseudoobscura. We also describe conservation of four distinct core promoter structures composed of combinations of elements at three positions. Overall, each type of genomic feature shows a characteristic divergence rate relative to neutral models, highlighting the value of multispecies alignment in annotating a target genome that should prove useful in the annotation of other high priority genomes, especially human and other mammalian genomes that are rich in noncoding sequences. We report that the vast majority of elements in the annotation are evolutionarily conserved, indicating that the annotation will be an important springboard for functional genetic testing by the Drosophila community.

Rare copy number variants (CNVs) disrupting ASTN2 or both ASTN2 and TRIM32 have been reported at 9q33.1 by genome-wide studies in a few individuals with neurodevelopmental disorders (NDDs). The vertebrate-specific astrotactins, ASTN2 and its paralog ASTN1, have key roles in glial-guided neuronal migration during brain development. To determine the prevalence of astrotactin mutations and delineate their associated phenotypic spectrum, we screened ASTN2/TRIM32 and ASTN1 (1q25.2) for exonic CNVs in clinical microarray data from 89 985 individuals across 10 sites, including 64 114 NDD subjects. In this clinical dataset, we identified 46 deletions and 12 duplications affecting ASTN2. Deletions of ASTN1 were much rarer. Deletions near the 3' terminus of ASTN2, which would disrupt all transcript isoforms (a subset of these deletions also included TRIM32), were significantly enriched in the NDD subjects (P = 0.002) compared with 44 085 population-based controls. Frequent phenotypes observed in individuals with such deletions include autism spectrum disorder (ASD), attention deficit hyperactivity disorder (ADHD), speech delay, anxiety and obsessive compulsive disorder (OCD). The 3'-terminal ASTN2 deletions were significantly enriched compared with controls in males with NDDs, but not in females. Upon quantifying ASTN2 human brain RNA, we observed shorter isoforms expressed from an alternative transcription start site of recent evolutionary origin near the 3' end. Spatiotemporal expression profiling in the human brain revealed consistently high ASTN1 expression while ASTN2 expression peaked in the early embryonic neocortex and postnatal cerebellar cortex. Our findings shed new light on the role of the astrotactins in psychopathology and their interplay in human neurodevelopment.

Exome sequence analysis of affected individuals from two families with autosomal-dominant inheritance of coloboma identified two different cosegregating heterozygous nonsense mutations (c.370C>T [p.Arg124*] and c. 1066G>T [p.Glu356*]) in YAP1. The phenotypes of the affected families differed in that one included no extraocular features and the other manifested with highly variable multisystem involvement, including hearing loss, intellectual disability, hematuria, and orofacial clefting. A combined LOD score of 4.2 was obtained for the association between YAP1 loss-of-function mutations and the phenotype in these families. YAP1 encodes an effector of the HIPPO-pathway-induced growth response, and whole-mount in situ hybridization in mouse embryos has shown that Yap1 is strongly expressed in the eye, brain, and fusing facial processes. RT-PCR showed that an alternative transcription start site (TSS) in intron 1 of YAP1 and Yap1 is widely used in human and mouse development, respectively. Transcripts from the alternative TSS are predicted to initiate at codon Met179 relative to the canonical transcript (RefSeq NM_001130145). In these alternative transcripts, the c.370C>T mutation in family 1305 is within the 5' UTR and cannot result in nonsense-mediated decay (NMD). The c. 1066G>T mutation in family 132 should result in NMD in transcripts from either TSS. Amelioration of the phenotype by the alternative transcripts provides a plausible explanation for the phenotypic differences between the families.

Genome-wide maps of DNase I hypersensitive sites (DHSs) reveal that most human promoters contain perpetually active cis-regulatory elements between -150 bp and +50 bp (-150/+50 bp) relative to the transcription start site (TSS). Transcription factors (TFs) recruit cofactors (chromatin remodelers, histone/protein-modifying enzymes, and scaffold proteins) to these elements in order to organize the local chromatin structure and coordinate the balance of post-translational modifications nearby, contributing to the overall regulation of transcription. However, the rules of TF-mediated cofactor recruitment to the -150/+50 bp promoter regions remain poorly understood. Here, we provide evidence for a general model in which a series of cis-regulatory elements (here termed 'cardinal' motifs) prefer acting individually, rather than in fixed combinations, within the -150/+50 bp regions to recruit TFs that dictate cofactor signatures distinctive of specific promoter subsets. Subsequently, human promoters can be subclassified based on the presence of cardinal elements and their associated cofactor signatures. In this study, furthermore, we have focused on promoters containing the nuclear respiratory factor 1 (NRF1) motif as the cardinal cis-regulatory element and have identified the pervasive association of NRF1 with the cofactor lysine-specific demethylase 1 (LSD1/KDM1A). This signature might be distinctive of promoters regulating nuclear-encoded mitochondrial and other particular genes in at least some cells. Together, we propose that decoding a signature-based, expanded model of control at proximal promoter regions should lead to a better understanding of coordinated regulation of gene transcription.

Class switch DNA recombination (CSR) crucially diversifies Ab biologic effector functions. 14-3-3γ specifically binds to the 5'-AGCT-3' repeats in the IgH locus switch (S) regions. By interacting directly with the C-terminal region of activation-induced cytidine deaminase (AID), 14-3-3γ targets this enzyme to S regions to mediate CSR. In this study, we showed that 14-3-3γ was expressed in germinal center B cells in vivo and induced in B cells by T-dependent and T-independent primary CSR-inducing stimuli in vitro in humans and mice. Induction of 14-3-3γ was rapid, peaking within 3 h of stimulation by LPSs, and sustained over the course of AID and CSR induction. It was dependent on recruitment of NF-κB to the 14-3-3γ gene promoter. The NF-κB recruitment enhanced the occupancy of the CpG island within the 14-3-3γ promoter by CFP1, a component of the COMPASS histone methyltransferase complex, and promoter-specific enrichment of histone 3 lysine 4 trimethylation (H3K4me3), which is indicative of open chromatin state and marks transcription-competent promoters. NF-κB also potentiated the binding of B cell lineage-specific factor E2A to an E-box motif located immediately downstream of the two closely-spaced transcription start sites for sustained 14-3-3γ expression and CSR induction. Thus, 14-3-3γ induction in CSR is enabled by the CFP1-mediated H3K4me3 enrichment in the promoter, dependent on NF-κB and sustained by E2A.

Bartonella quintana is a vector-borne bacterial pathogen that causes fatal disease in humans. During the infectious cycle, B. quintana transitions from the hemin-restricted human bloodstream to the hemin-rich body louse vector. Because extracytoplasmic function (ECF) sigma factors often regulate adaptation to environmental changes, we hypothesized that a previously unstudied B. quintana ECF sigma factor, RpoE, is involved in the transition from the human host to the body louse vector. The genomic context of B. quintana rpoE identified it as a member of the ECF15 family of sigma factors found only in alphaproteobacteria. ECF15 sigma factors are believed to be the master regulators of the general stress response in alphaproteobacteria. In this study, we examined the B. quintana RpoE response to two stressors that are encountered in the body louse vector environment, a decreased temperature and an increased hemin concentration. We determined that the expression of rpoE is significantly upregulated at the body louse (28°C) versus the human host (37°C) temperature. rpoE expression also was upregulated when B. quintana was exposed to high hemin concentrations. In vitro and in vivo analyses demonstrated that RpoE function is regulated by a mechanism involving the anti-sigma factor NepR and the response regulator PhyR. The ΔrpoE ΔnepR mutant strain of B. quintana established that RpoE-mediated transcription is important in mediating the tolerance of B. quintana to high hemin concentrations. We present the first analysis of an ECF15 sigma factor in a vector-borne human pathogen and conclude that RpoE has a role in the adaptation of B. quintana to the hemin-rich arthropod vector environment.

Eukaryotic transcription initiation requires the assembly of general transcription factors into a pre-initiation complex that ensures the accurate loading of RNA polymerase II (Pol II) at the transcription start site. The molecular mechanism and function of this assembly have remained elusive due to lack of structural information. Here we have used an in vitro reconstituted system to study the stepwise assembly of human TBP, TFIIA, TFIIB, Pol II, TFIIF, TFIIE and TFIIH onto promoter DNA using cryo-electron microscopy. Our structural analyses provide pseudo-atomic models at various stages of transcription initiation that illuminate critical molecular interactions, including how TFIIF engages Pol II and promoter DNA to stabilize both the closed pre-initiation complex and the open-promoter complex, and to regulate start--initiation complexes, combined with the localization of the TFIIH helicases XPD and XPB, support a DNA translocation model of XPB and explain its essential role in promoter opening.

5'- and 3'-untranslated regions (UTRs) are important regulators of gene expression and play key roles in disease progression and susceptibility. The 5'-UTR of the fragile X mental retardation 1 (FMR1) gene contains a CGG repeat element that is expanded (>200 CGG repeats; full mutation) and methylated in fragile X syndrome (FXS), the most common form of inherited intellectual disability (ID) and known cause of autism. Significant phenotypic involvement has also emerged in some individuals with the premutation (55-200 CGG repeats), including fragile X-associated premature ovarian insufficiency (FXPOI) in females, and the neurodegenerative disorder, fragile X-associated tremor/ataxia syndrome (FXTAS), in older adult carriers. Here, we show that FMR1 mRNA in human and mouse brain is expressed as a combination of multiple isoforms that use alternative transcriptional start sites and different polyadenylation sites. Furthermore, we have identified a novel human transcription start site used in brain but not in lymphoblastoid cells, and have detected FMR1 isoforms generated through the use of both canonical and non-canonical polyadenylation signals. Importantly, in both human and mouse, a specific regulation of the UTRs is observed in brain of FMR1 premutation alleles, suggesting that the transcript variants may play a role in premutation-related pathologies.

BackgroundThe Salmonella PreA/PreB two-component system (TCS) is an ortholog of the QseBC TCS of Escherichia coli. In both Salmonella and E. coli, this system has been shown to affect motility and virulence in response to quorum-sensing and hormonal signals, and to affect the transcription of the Salmonella enterica serovar Typhimurium (S. Typhimurium) pmrAB operon, which encodes an important virulence-associated TCS.ResultsTo determine the PreA/PreB regulon in S. Typhimurium, we performed DNA microarrays comparing the wild type strain and various preA and/or preB mutants in the presence of ectopically expressed preA (qseB). These data confirmed our previous findings of the negative effect of PreB on PreA gene regulation and identified candidate PreA-regulated genes. A proportion of the activated loci were previously identified as PmrA-activated genes (yibD, pmrAB, cptA, etc.) or were genes located in the local region around preA, including the preAB operon. The transcriptional units were defined in this local region by RT-PCR, suggesting three PreA activated operons composed of preA-preB, mdaB-ygiN, and ygiW-STM3175. Several putative virulence-related phenotypes were examined for preAB mutants, resulting in the observation of a host cell invasion and slight virulence defect of a preAB mutant. Contrary to previous reports on this TCS, we were unable to show a PreA/PreB-dependent effect of the quorum-sensing signal AI-2 or of epinephrine on S. Typhimurium with regard to bacterial motility.ConclusionThis work further characterizes this unorthadox OmpR/EnvZ class TCS and provides novel candidate regulated genes for further study. This first in-depth study of the PreA/PreB regulatory system phenotypes and regulation suggests significant comparative differences to the reported function of the orthologous QseB/QseC in E. coli.

BackgroundDespite an intense interest in the biological functions of the phosphoinositide 3-kinase (PI3K) signalling enzymes, little is known about the regulation of PI3K gene expression. This also applies to the leukocyte-enriched p110delta catalytic subunit of PI3K, an enzyme that has attracted widespread interest because of its role in immunity and allergy.Principal findingsWe show that p110delta expression is mainly regulated at the transcriptional level. In fibroblasts, lymphocytes and myeloid cells, p110delta gene transcription appears to be constitutive and not subject to acute stimulation. 5'RACE experiments revealed that p110delta mRNA transcripts contain distinct upstream untranslated exons (named exon -1, -2a, -2b, -2c and -2d), which are located up to 81 kb upstream of the translational start codon in exon 1. The levels of all the different p110delta transcripts are higher in leukocytes compared to non-leukocytes, with the p110delta transcript containing exon -2a most abundantly expressed. We have identified a highly conserved transcription factor (TF) binding cluster in the p110delta gene which has enhanced promoter activity in leukocytes compared to non-leukocytes. In human, this TF cluster is located immediately upstream of exon -2a whilst in mouse, it is located within exon -2a.ConclusionThis study identifies a conserved PIK3CD promoter region that may account for the predominant leukocyte expression of p110delta.