Although transcriptional elongation by RNA polymerase II is coupled with many RNA-related processes, genomewide elongation rates remain unknown. We describe a method, called 4sUDRB-seq, based on reversible inhibition of transcription elongation coupled with tagging newly transcribed RNA with 4-thiouridine and high throughput sequencing to measure simultaneously with high confidence genome-wide transcription elongation rates in cells. We find that most genes are transcribed at about 3.5 Kb/min, with elongation rates varying between 2 Kb/min and 6 Kb/min. 4sUDRB-seq can facilitate genomewide exploration of the involvement of specific elongation factors in transcription and the contribution of deregulated transcription elongation to various pathologies.

Genome replication introduces a stepwise increase in the DNA template available for transcription. Genes replicated early in S phase experience this increase before late-replicating genes, raising the question of how expression levels are affected by DNA replication. We show that in budding yeast, messenger RNA (mRNA) synthesis rate is buffered against changes in gene dosage during S phase. This expression homeostasis depends on acetylation of H3 on its internal K56 site by Rtt109/Asf1. Deleting these factors, mutating H3K56 or up-regulating its deacetylation, increases gene expression in S phase in proportion to gene replication timing. Therefore, H3K56 acetylation on newly deposited histones reduces transcription efficiency from replicated DNA, complementing its role in guarding genome stability. Our study provides molecular insight into the mechanism maintaining expression homeostasis during DNA replication.

Structural variants and presence/absence polymorphisms are common in plant genomes, yet they are routinely overlooked in genome-wide association studies (GWAS). Here, we expand the type of genetic variants detected in GWAS to include major deletions, insertions, and rearrangements. We first use raw sequencing data directly to derive short sequences, k -mers, that mark a broad range of polymorphisms independently of a reference genome. We then link k -mers associated with phenotypes to specific genomic regions. Using this approach, we re-analyzed 2,000 traits in Arabidopsis thaliana, tomato, and maize populations. Associations identified with k -mers recapitulate those found with single-nucleotide polymorphisms (SNPs), but with stronger statistical support. Importantly, we discovered new associations with structural variants and with regions missing from reference genomes. Our results demonstrate the power of performing GWAS before linking sequence reads to specific genomic regions, which allows detection of a wider range of genetic variants responsible for phenotypic variation.

Gene expression shows a significant variation (noise) between genetically identical cells. Noise depends on the gene expression process regulated by the chromatin environment. We screened for chromatin factors that modulate noise in S. cerevisiae and analyzed the results using a theoretical model that infers regulatory mechanisms from the noise vs. mean relationship. Distinct activities of the Rpd3(L) and Set3 histone deacetylase complexes were predicted. Both HDACs repressed expression. Yet, Rpd3(L)C decreased the frequency of transcriptional bursts, while Set3C decreased the burst size, as did H2B mono-ubiquitination (ubH2B). We mapped the acetylation of H3 Lysine 9 (H3K9ac) upon deletion of multiple subunits of Set3C and Rpd3(L)C, and of ubH2B effectors. ubH2B and Set3C appear to function in the same pathway to reduce the probability that an elongating PolII produces a functional transcript (PolII processivity), while Rpd3(L)C likely represses gene expression at a step preceding elongation.

Various histone modifications decorate nucleosomes within transcribed genes. Among these, monoubiquitylation of histone H2B (H2Bub1) and methylation of histone H3 on lysines 36 (H3K36me2/3) and 79 (H3K79me2/3) correlate positively with gene expression. By measuring the progression of the transcriptional machinery along genes within live cells, we now report that H2B monoubiquitylation occurs cotranscriptionally and accurately reflects the advance of RNA polymerase II (Pol II). In contrast, H3K36me3 and H3K79me2 are less dynamic and represent Pol II movement less faithfully. High-resolution ChIP-seq reveals that H2Bub1 levels are selectively reduced at exons and decrease in an exon-dependent stepwise manner toward the 39 end of genes. Exonic depletion of H2Bub1 in gene bodies is highly correlated with Pol II pausing at exons, suggesting elongation rate changes associated with intron-exon structure. In support of this notion, H2Bub1 levels were found to be significantly correlated with transcription elongation rates measured in various cell lines. Overall, our data shed light on the organization of H2Bub1 within transcribed genes and single out H2Bub1 as a reliable marker for ongoing transcription elongation.