Sites of transcription initiation drive mRNA isoform selection.

The generation of distinct messenger RNA isoforms through alternative RNA processing modulates the expression and function of genes, often in a cell-type-specific manner. Here, we assess the regulatory relationships between transcription initiation, alternative splicing, and 3' end site selection. Applying long-read sequencing to accurately represent even the longest transcripts from end to end, we quantify mRNA isoforms in Drosophila tissues, including the transcriptionally complex nervous system. We find that in Drosophila heads, as well as in human cerebral organoids, 3' end site choice is globally influenced by the site of transcription initiation (TSS). "Dominant promoters," characterized by specific epigenetic signatures including p300/CBP binding, impose a transcriptional constraint to define splice and polyadenylation variants. In vivo deletion or overexpression of dominant promoters as well as p300/CBP loss disrupted the 3' end expression landscape. Our study demonstrates the crucial impact of TSS choice on the regulation of transcript diversity and tissue identity.

Exon-Mediated Activation of Transcription Starts.

The processing of RNA transcripts from mammalian genes occurs in proximity to their transcription. Here, we describe a phenomenon affecting thousands of genes that we call exon-mediated activation of transcription starts (EMATS), in which the splicing of internal exons impacts promoter choice and the expression level of the gene. We observed that evolutionary gain of internal exons is associated with gain of new transcription start sites (TSSs) nearby and increased gene expression. Inhibiting exon splicing reduced transcription from nearby promoters, and creation of new spliced exons activated transcription from cryptic promoters. The strongest effects occurred for weak promoters located proximal and upstream of efficiently spliced exons. Together, our findings support a model in which splicing recruits transcription machinery locally to influence TSS choice and identify exon gain, loss, and regulatory change as major contributors to the evolution of alternative promoters and gene expression in mammals.

Sub-nucleosomal Genome Structure Reveals Distinct Nucleosome Folding Motifs.

Elucidating the global and local rules that govern genome-wide, hierarchical chromatin architecture remains a critical challenge. Current high-throughput chromosome conformation capture (Hi-C) technologies have identified large-scale chromatin structural motifs, such as topologically associating domains and looping. However, structural rules at the smallest or nucleosome scale remain poorly understood. Here, we coupled nucleosome-resolved Hi-C technology with simulated annealing-molecular dynamics (SA-MD) simulation to reveal 3D spatial distributions of nucleosomes and their genome-wide orientation in chromatin. Our method, called Hi-CO, revealed distinct nucleosome folding motifs across the yeast genome. Our results uncovered two types of basic secondary structural motifs in nucleosome folding: α-tetrahedron and ß-rhombus analogous to α helix and ß sheet motifs in protein folding. Using mutants and cell-cycle-synchronized cells, we further uncovered motifs with specific nucleosome positioning and orientation coupled to epigenetic features at individual loci. By illuminating molecular-level structure-function relationships in eukaryotic chromatin, our findings establish organizational principles of nucleosome folding.

Dual Roles of Poly(dA:dT) Tracts in Replication Initiation and Fork Collapse.

Replication origins, fragile sites, and rDNA have been implicated as sources of chromosomal instability. However, the defining genomic features of replication origins and fragile sites are among the least understood elements of eukaryote genomes. Here, we map sites of replication initiation and breakage in primary cells at high resolution. We find that replication initiates between transcribed genes within nucleosome-depleted structures established by long asymmetrical poly(dA:dT) tracts flanking the initiation site. Paradoxically, long (>20 bp) (dA:dT) tracts are also preferential sites of polar replication fork stalling and collapse within early-replicating fragile sites (ERFSs) and late-replicating common fragile sites (CFSs) and at the rDNA replication fork barrier. Poly(dA:dT) sequences are fragile because long single-strand poly(dA) stretches at the replication fork are unprotected by the replication protein A (RPA). We propose that the evolutionary expansion of poly(dA:dT) tracts in eukaryotic genomes promotes replication initiation, but at the cost of chromosome fragility.

Transcriptome-wide Interrogation of the Functional Intronome by Spliceosome Profiling.

Full understanding of eukaryotic transcriptomes and how they respond to different conditions requires deep knowledge of all sites of intron excision. Although RNA sequencing (RNA-seq) provides much of this information, the low abundance of many spliced transcripts (often due to their rapid cytoplasmic decay) limits the ability of RNA-seq alone to reveal the full repertoire of spliced species. Here, we present "spliceosome profiling," a strategy based on deep sequencing of RNAs co-purifying with late-stage spliceosomes. Spliceosome profiling allows for unambiguous mapping of intron ends to single-nucleotide resolution and branchpoint identification at unprecedented depths. Our data reveal hundreds of new introns in S. pombe and numerous others that were previously misannotated. By providing a means to directly interrogate sites of spliceosome assembly and catalysis genome-wide, spliceosome profiling promises to transform our understanding of RNA processing in the nucleus, much as ribosome profiling has transformed our understanding mRNA translation in the cytoplasm.

Swi/Snf chromatin remodeling regulates transcriptional interference and gene repression.

Alternative transcription start sites can affect transcript isoform diversity and translation levels. In a recently described form of gene regulation, coordinated transcriptional and translational interference results in transcript isoform-dependent changes in protein expression. Specifically, a long undecoded transcript isoform (LUTI) is transcribed from a gene-distal promoter, interfering with expression of the gene-proximal promoter. Although transcriptional and chromatin features associated with LUTI expression have been described, the mechanism underlying LUTI-based transcriptional interference is not well understood. Using an unbiased genetic approach followed by functional genomics, we uncovered that the Swi/Snf chromatin remodeling complex is required for co-transcriptional nucleosome remodeling that leads to LUTI-based repression. We identified genes with tandem promoters that rely on Swi/Snf function for transcriptional interference during protein folding stress, including LUTI-regulated genes. This study provides clear evidence for Swi/Snf playing a direct role in gene repression via a cis transcriptional interference mechanism.

Insights into Nucleosome Organization in Mouse Embryonic Stem Cells through Chemical Mapping.

Nucleosome organization influences gene activity by controlling DNA accessibility to transcription machinery. Here, we develop a chemical biology approach to determine mammalian nucleosome positions genome-wide. We uncovered surprising features of nucleosome organization in mouse embryonic stem cells. In contrast to the prevailing model, we observe that for nearly all mouse genes, a class of fragile nucleosomes occupies previously designated nucleosome-depleted regions around transcription start sites and transcription termination sites. We show that nucleosomes occupy DNA targets for a subset of DNA-binding proteins, including CCCTC-binding factor (CTCF) and pluripotency factors. Furthermore, we provide evidence that promoter-proximal nucleosomes, with the +1 nucleosome in particular, contribute to the pausing of RNA polymerase II. Lastly, we find a characteristic preference for nucleosomes at exon-intron junctions. Taken together, we establish an accurate method for defining the nucleosome landscape and provide a valuable resource for studying nucleosome-mediated gene regulation in mammalian cells.

RNA Polymerase II Regulates Topoisomerase 1 Activity to Favor Efficient Transcription.

We report a mechanism through which the transcription machinery directly controls topoisomerase 1 (TOP1) activity to adjust DNA topology throughout the transcription cycle. By comparing TOP1 occupancy using chromatin immunoprecipitation sequencing (ChIP-seq) versus TOP1 activity using topoisomerase 1 sequencing (TOP1-seq), a method reported here to map catalytically engaged TOP1, TOP1 bound at promoters was discovered to become fully active only after pause-release. This transition coupled the phosphorylation of the carboxyl-terminal-domain (CTD) of RNA polymerase II (RNAPII) with stimulation of TOP1 above its basal rate, enhancing its processivity. TOP1 stimulation is strongly dependent on the kinase activity of BRD4, a protein that phosphorylates Ser2-CTD and regulates RNAPII pause-release. Thus the coordinated action of BRD4 and TOP1 overcame the torsional stress opposing transcription as RNAPII commenced elongation but preserved negative supercoiling that assists promoter melting at start sites. This nexus between transcription and DNA topology promises to elicit new strategies to intercept pathological gene expression.

Structural basis of transcription reduction by a promoter-proximal +1 nucleosome.

At active human genes, the +1 nucleosome is located downstream of the RNA polymerase II (RNA Pol II) pre-initiation complex (PIC). However, at inactive genes, the +1 nucleosome is found further upstream, at a promoter-proximal location. Here, we establish a model system to show that a promoter-proximal +1 nucleosome can reduce RNA synthesis in vivo and in vitro, and we analyze its structural basis. We find that the PIC assembles normally when the edge of the +1 nucleosome is located 18 base pairs (bp) downstream of the transcription start site (TSS). However, when the nucleosome edge is located further upstream, only 10 bp downstream of the TSS, the PIC adopts an inhibited state. The transcription factor IIH (TFIIH) shows a closed conformation and its subunit XPB contacts DNA with only one of its two ATPase lobes, inconsistent with DNA opening. These results provide a mechanism for nucleosome-dependent regulation of transcription initiation.

N6-methyldeoxyadenosine marks active transcription start sites in Chlamydomonas.

N(6)-methyldeoxyadenosine (6mA or m(6)A) is a DNA modification preserved in prokaryotes to eukaryotes. It is widespread in bacteria and functions in DNA mismatch repair, chromosome segregation, and virulence regulation. In contrast, the distribution and function of 6mA in eukaryotes have been unclear. Here, we present a comprehensive analysis of the 6mA landscape in the genome of Chlamydomonas using new sequencing approaches. We identified the 6mA modification in 84% of genes in Chlamydomonas. We found that 6mA mainly locates at ApT dinucleotides around transcription start sites (TSS) with a bimodal distribution and appears to mark active genes. A periodic pattern of 6mA deposition was also observed at base resolution, which is associated with nucleosome distribution near the TSS, suggesting a possible role in nucleosome positioning. The new genome-wide mapping of 6mA and its unique distribution in the Chlamydomonas genome suggest potential regulatory roles of 6mA in gene expression in eukaryotic organisms.

A One Precursor One siRNA Model for Pol IV-Dependent siRNA Biogenesis.

RNA-directed DNA methylation in Arabidopsis thaliana is driven by the plant-specific RNA Polymerase IV (Pol IV). It has been assumed that a Pol IV transcript can give rise to multiple 24-nt small interfering RNAs (siRNAs) that target DNA methylation. Here, we demonstrate that Pol IV-dependent RNAs (P4RNAs) from wild-type Arabidopsis are surprisingly short in length (30 to 40 nt) and mirror 24-nt siRNAs in distribution, abundance, strand bias, and 5'-adenine preference. P4RNAs exhibit transcription start sites similar to Pol II products and are featured with 5'-monophosphates and 3'-misincorporated nucleotides. The 3'-misincorporation preferentially occurs at methylated cytosines on the template DNA strand, suggesting a co-transcriptional feedback to siRNA biogenesis by DNA methylation to reinforce silencing locally. These results highlight an unusual mechanism of Pol IV transcription and suggest a "one precursor, one siRNA" model for the biogenesis of 24-nt siRNAs in Arabidopsis.

Position-dependent function of human sequence-specific transcription factors.

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.

Incomplete transcripts dominate the Mycobacterium tuberculosis transcriptome.

Mycobacterium tuberculosis (Mtb) is a bacterial pathogen that causes tuberculosis (TB), an infectious disease that is responsible for major health and economic costs worldwide1. Mtb encounters diverse environments during its life cycle and responds to these changes largely by reprogramming its transcriptional output2. However, the mechanisms of Mtb transcription and how they are regulated remain poorly understood. Here we use a sequencing method that simultaneously determines both termini of individual RNA molecules in bacterial cells3 to profile the Mtb transcriptome at high resolution. Unexpectedly, we find that most Mtb transcripts are incomplete, with their 5' ends aligned at transcription start sites and 3' ends located 200-500 nucleotides downstream. We show that these short RNAs are mainly associated with paused RNA polymerases (RNAPs) rather than being products of premature termination. We further show that the high propensity of Mtb RNAP to pause early in transcription relies on the binding of the σ-factor. Finally, we show that a translating ribosome promotes transcription elongation, revealing a potential role for transcription-translation coupling in controlling Mtb gene expression. In sum, our findings depict a mycobacterial transcriptome that prominently features incomplete transcripts resulting from RNAP pausing. We propose that the pausing phase constitutes an important transcriptional checkpoint in Mtb that allows the bacterium to adapt to environmental changes and could be exploited for TB therapeutics.

Genome-wide DNA-methylation landscape defines specialization of regulatory T cells in tissues.

Regulatory T cells (Treg cells) perform two distinct functions: they maintain self-tolerance, and they support organ homeostasis by differentiating into specialized tissue Treg cells. We found that epigenetic modifications defined the molecular characteristics of tissue Treg cells. Tagmentation-based whole-genome bisulfite sequencing revealed more than 11,000 regions that were methylated differentially in pairwise comparisons of tissue Treg cell populations and lymphoid T cells. Similarities in the epigenetic landscape led to the identification of a common tissue Treg cell population that was present in many organs and was characterized by gain and loss of DNA methylation that included many gene sites associated with the TH2 subset of helper T cells, such as the gene encoding cytokine IL-33 receptor ST2, as well as the production of tissue-regenerative factors. Furthermore, the ST2-expressing population was dependent on the transcriptional regulator BATF and could be expanded by IL-33. Thus, tissue Treg cells integrate multiple waves of epigenetic reprogramming that define their tissue-restricted specialization.

Convergent transcription at intragenic super-enhancers targets AID-initiated genomic instability.

Activation-induced cytidine deaminase (AID) initiates both somatic hypermutation (SHM) for antibody affinity maturation and DNA breakage for antibody class switch recombination (CSR) via transcription-dependent cytidine deamination of single-stranded DNA targets. Though largely specific for immunoglobulin genes, AID also acts on a limited set of off-targets, generating oncogenic translocations and mutations that contribute to B cell lymphoma. How AID is recruited to off-targets has been a long-standing mystery. Based on deep GRO-seq studies of mouse and human B lineage cells activated for CSR or SHM, we report that most robust AID off-target translocations occur within highly focal regions of target genes in which sense and antisense transcription converge. Moreover, we found that such AID-targeting "convergent" transcription arises from antisense transcription that emanates from super-enhancers within sense transcribed gene bodies. Our findings provide an explanation for AID off-targeting to a small subset of mostly lineage-specific genes in activated B cells.

Argonaute-bound small RNAs from promoter-proximal RNA polymerase II.

Argonaute (Ago) proteins mediate posttranscriptional gene repression by binding guide miRNAs to regulate targeted RNAs. To confidently assess Ago-bound small RNAs, we adapted a mouse embryonic stem cell system to express a single epitope-tagged Ago protein family member in an inducible manner. Here, we report the small RNA profile of Ago-deficient cells and show that Ago-dependent stability is a common feature of mammalian miRNAs. Using this criteria and immunopurification, we identified an Ago-dependent class of noncanonical miRNAs derived from protein-coding gene promoters, which we name transcriptional start site miRNAs (TSS-miRNAs). A subset of promoter-proximal RNA polymerase II (RNAPII) complexes produces hairpin RNAs that are processed in a DiGeorge syndrome critical region gene 8 (Dgcr8)/Drosha-independent but Dicer-dependent manner. TSS-miRNA activity is detectable from endogenous levels and following overexpression of mRNA constructs. Finally, we present evidence of differential expression and conservation in humans, suggesting important roles in gene regulation.

A variety of dicer substrates in human and C. elegans.

The endoribonuclease Dicer is known for its central role in the biogenesis of eukaryotic small RNAs/microRNAs. Despite its importance, Dicer target transcripts have not been directly mapped. Here, we apply biochemical methods to human cells and C. elegans and identify thousands of Dicer-binding sites. We find known and hundreds of additional miRNAs with high sensitivity and specificity. We also report structural RNAs, promoter RNAs, and mitochondrial transcripts as Dicer targets. Interestingly, most Dicer-binding sites reside on mRNAs/lncRNAs and are not significantly processed into small RNAs. These passive sites typically harbor small, Dicer-bound hairpins within intact transcripts and generally stabilize target expression. We show that passive sites can sequester Dicer and reduce microRNA expression. mRNAs with passive sites were in human and worm significantly associated with processing-body/granule function. Together, we provide the first transcriptome-wide map of Dicer targets and suggest conserved binding modes and functions outside of the miRNA pathway.

Acetyl-methyllysine marks chromatin at active transcription start sites.

Lysine residues in histones and other proteins can be modified by post-translational modifications that encode regulatory information1. Lysine acetylation and methylation are especially important for regulating chromatin and gene expression2-4. Pathways involving these post-translational modifications are targets for clinically approved therapeutics to treat human diseases. Lysine methylation and acetylation are generally assumed to be mutually exclusive at the same residue. Here we report cellular lysine residues that are both methylated and acetylated on the same side chain to form Nε-acetyl-Nε-methyllysine (Kacme). We show that Kacme is found on histone H4 (H4Kacme) across a range of species and across mammalian tissues. Kacme is associated with marks of active chromatin, increased transcriptional initiation and is regulated in response to biological signals. H4Kacme can be installed by enzymatic acetylation of monomethyllysine peptides and is resistant to deacetylation by some HDACs in vitro. Kacme can be bound by chromatin proteins that recognize modified lysine residues, as we demonstrate with the crystal structure of acetyllysine-binding protein BRD2 bound to a histone H4Kacme peptide. These results establish Kacme as a cellular post-translational modification with the potential to encode information distinct from methylation and acetylation alone and demonstrate that Kacme has all the hallmarks of a post-translational modification with fundamental importance to chromatin biology.

Genome-wide mapping of human DNA replication by optical replication mapping supports a stochastic model of eukaryotic replication.

The heterogeneous nature of eukaryotic replication kinetics and the low efficiency of individual initiation sites make mapping the location and timing of replication initiation in human cells difficult. To address this challenge, we have developed optical replication mapping (ORM), a high-throughput single-molecule approach, and used it to map early-initiation events in human cells. The single-molecule nature of our data and a total of >2,500-fold coverage of the human genome on 27 million fibers averaging ∼300 kb in length allow us to identify initiation sites and their firing probability with high confidence. We find that the distribution of human replication initiation is consistent with inefficient, stochastic activation of heterogeneously distributed potential initiation complexes enriched in accessible chromatin. These observations are consistent with stochastic models of initiation-timing regulation and suggest that stochastic regulation of replication kinetics is a fundamental feature of eukaryotic replication, conserved from yeast to humans.

Epigenome plasticity in plants.

Plant intra-individual and inter-individual variation can be determined by the epigenome, a set of covalent modifications of DNA and chromatin that can alter genome structure and activity without changes to the genome sequence. The epigenome of plant cells is plastic, that is, it can change in response to internal or external cues, such as during development or due to environmental changes, to create a memory of such events. Ongoing advances in technologies to read and write epigenomic patterns with increasing resolution, scale and precision are enabling the extent of plant epigenome variation to be more extensively characterized and functionally interrogated. In this Review, we discuss epigenome dynamics and variation within plants during development and in response to environmental changes, including stress, as well as between plants. We review known or potential functions of such plasticity and emphasize the importance of investigating the causality of epigenomic changes. Finally, we discuss emerging technologies that may underpin future research into plant epigenome plasticity.