Pore models of membrane fusion postulate that cylinders of integral membrane proteins can initiate a fusion pore after conformational rearrangement of pore subunits. In the fusion of yeast vacuoles, V-ATPase V0 sectors, which contain a central cylinder of membrane integral proteolipid subunits, associate to form a transcomplex that might resemble an intermediate postulated in some pore models. We tested the role of V0 sectors in vacuole fusion. V0 functions in fusion and proton translocation could be experimentally separated via the differential effects of mutations and inhibitory antibodies. Inactivation of the V0 subunit Vph1p blocked fusion in the terminal reaction stage that is independent of a proton gradient. Deltavph1 mutants were capable of docking and trans-SNARE pairing and of subsequent release of lumenal Ca2+, but they did not fuse. The Ca2+-releasing channel appears to be tightly coupled to V0 because inactivation of Vph1p by antibodies blocked Ca2+ release. Vph1 deletion on only one fusion partner sufficed to severely reduce fusion activity. The functional requirement for Vph1p correlates to V0 transcomplex formation in that both occur after docking and Ca2+ release. These observations establish V0 as a crucial factor in vacuole fusion acting downstream of trans-SNARE pairing.

Nuclei bind yeast vacuoles via nucleus-vacuole (NV) junctions. Under nutrient restriction, NV junctions invaginate and release vesicles filled with nuclear material into vacuoles, resulting in piecemeal microautophagy of the nucleus (PMN). We show that the electrochemical gradient across the vacuolar membrane promotes invagination of NV junctions. Existing invaginations persist independently of the gradient, but final release of PMN vesicles requires again V-ATPase activity. We find that NV junctions form a diffusion barrier on the vacuolar membrane that excludes V-ATPase but is enriched in the VTC complex and accessible to other membrane-integral proteins. V-ATPase exclusion depends on the NV junction proteins Nvj1p,Vac8p, and the electrochemical gradient. It also depends on factors of lipid metabolism, such as the oxysterol binding protein Osh1p and the enoyl-CoA reductase Tsc13p, which are enriched in NV junctions, and on Lag1p and Fen1p. Our observations suggest that NV junctions form in two separable steps: Nvj1p and Vac8p suffice to establish contact between the two membranes. The electrochemical potential and lipid-modifying enzymes are needed to establish the vacuolar diffusion barrier, invaginate NV junctions, and form PMN vesicles.

RNA polymerase (Pol) II transcribes protein-coding genes in the nucleus of eukaryotic cells and consists of 12 polypeptide subunits. It is unknown how Pol II is imported into the nucleus. Here we show that Pol II nuclear import requires the protein Iwr1 and provide evidence for cyclic Iwr1 function. Iwr1 binds Pol II in the active center cleft between the two largest subunits, maybe facilitating or sensing complete Pol II assembly in the cytoplasm. Iwr1 then uses an N-terminal bipartite nuclear localization signal that is recognized by karyopherin α to direct Pol II nuclear import. In the nucleus, Iwr1 is displaced from Pol II by transcription initiation factors and nucleic acids, enabling its export and recycling. Iwr1 function is Pol II specific, transcription independent, and apparently conserved from yeast to human.

Mediator is a modular multiprotein complex required for regulated transcription by RNA polymerase (Pol) II. Here, we show that the middle module of the Mediator core contains a submodule of unique structure and function that comprises the N-terminal part of subunit Med7 (Med7N) and the highly conserved subunit Med31 (Soh1). The Med7N/31 submodule shows a conserved novel fold, with two proline-rich stretches in Med7N wrapping around the right-handed four-helix bundle of Med31. In vitro, Med7N/31 is required for activated transcription and can act in trans when added exogenously. In vivo, Med7N/31 has a predominantly positive function on the expression of a specific subset of genes, including genes involved in methionine metabolism and iron transport. Comparative phenotyping and transcriptome profiling identify specific and overlapping functions of different Mediator submodules.

The rates of mRNA synthesis and degradation determine cellular mRNA levels and can be monitored by comparative dynamic transcriptome analysis (cDTA) that uses nonperturbing metabolic RNA labeling. Here we present cDTA data for 46 yeast strains lacking genes involved in mRNA degradation and metabolism. In these strains, changes in mRNA degradation rates are generally compensated by changes in mRNA synthesis rates, resulting in a buffering of mRNA levels. We show that buffering of mRNA levels requires the RNA exonuclease Xrn1. The buffering is rapidly established when mRNA synthesis is impaired, but is delayed when mRNA degradation is impaired, apparently due to Xrn1-dependent transcription repressor induction. Cluster analysis of the data defines the general mRNA degradation machinery, reveals different substrate preferences for the two mRNA deadenylase complexes Ccr4-Not and Pan2-Pan3, and unveils an interwoven cellular mRNA surveillance network.

Biogenesis of eukaryotic messenger ribonucleoprotein complexes (mRNPs) involves the synthesis, splicing, and 3' processing of pre-mRNA, and the assembly of mature mRNPs for nuclear export. We mapped 23 mRNP biogenesis factors onto the yeast transcriptome, providing 10(4)-10(6) high-confidence RNA interaction sites per factor. The data reveal how mRNP biogenesis factors recognize pre-mRNA elements in vivo. They define conserved interactions between splicing factors and pre-mRNA introns, including the recognition of intron-exon junctions and the branchpoint. They also identify a unified arrangement of 3' processing factors at pre-mRNA polyadenylation (pA) sites in yeast and human, which results from an A-U sequence bias at pA sites. Global data analysis indicates that 3' processing factors have roles in splicing and RNA surveillance, and that they couple mRNP biogenesis events to restrict nuclear export to mature mRNPs.

At the 3' ends of protein-coding genes, RNA polymerase (Pol) II is dephosphorylated at tyrosine residues (Tyr1) of its C-terminal domain (CTD). In addition, the associated cleavage-and-polyadenylation factor (CPF) cleaves the transcript and adds a poly(a) tail. Whether these events are coordinated and how they lead to transcription termination remains poorly understood. Here we show that CPF from Saccharomyces cerevisiae is a Pol II-CTD phosphatase and that the CPF subunit Glc7 dephosphorylates Tyr1 in vitro. In vivo, the activity of Glc7 is required for normal Tyr1 dephosphorylation at the polyadenylation site, for recruitment of termination factors Pcf11 and Rtt103 and for normal Pol II termination. These results show that transcription termination involves Tyr1 dephosphorylation of the CTD and indicate that pre-mRNA processing by CPF and transcription termination are coupled via Glc7-dependent Pol II-Tyr1 dephosphorylation.

We present One Hand Clapping (OHC), a method for the detection of condition-specific interactions between transcription factors (TFs) from genome-wide gene activity measurements. OHC is based on a mapping between transcription factors and their target genes. Given a single case-control experiment, it uses a linear regression model to assess whether the common targets of two arbitrary TFs behave differently than expected from the genes targeted by only one of the TFs. When applied to osmotic stress data in S. cerevisiae, OHC produces consistent results across three types of expression measurements: gene expression microarray data, RNA Polymerase II ChIP-chip binding data and messenger RNA synthesis rates. Among the eight novel, condition-specific TF pairs, we validate the interaction between Gcn4p and Arr1p experimentally. We apply OHC to a large gene activity dataset in S. cerevisiae and provide a compendium of condition-specific TF interactions.

To monitor eukaryotic mRNA metabolism, we developed comparative dynamic transcriptome analysis (cDTA). cDTA provides absolute rates of mRNA synthesis and decay in Saccharomyces cerevisiae (Sc) cells with the use of Schizosaccharomyces pombe (Sp) as an internal standard. cDTA uses nonperturbing metabolic labeling that supersedes conventional methods for mRNA turnover analysis. cDTA reveals that Sc and Sp transcripts that encode orthologous proteins have similar synthesis rates, whereas decay rates are fivefold lower in Sp, resulting in similar mRNA concentrations despite the larger Sp cell volume. cDTA of Sc mutants reveals that a eukaryote can buffer mRNA levels. Impairing transcription with a point mutation in RNA polymerase (Pol) II causes decreased mRNA synthesis rates as expected, but also decreased decay rates. Impairing mRNA degradation by deleting deadenylase subunits of the Ccr4-Not complex causes decreased decay rates as expected, but also decreased synthesis rates. Extended kinetic modeling reveals mutual feedback between mRNA synthesis and degradation that may be achieved by a factor that inhibits synthesis and enhances degradation.

The Mediator is a highly conserved, large multiprotein complex that is involved essentially in the regulation of eukaryotic mRNA transcription. It acts as a general transcription factor by integrating regulatory signals from gene-specific activators or repressors to the RNA Polymerase II. The internal network of interactions between Mediator subunits that conveys these signals is largely unknown. Here, we introduce MC EMiNEM, a novel method for the retrieval of functional dependencies between proteins that have pleiotropic effects on mRNA transcription. MC EMiNEM is based on Nested Effects Models (NEMs), a class of probabilistic graphical models that extends the idea of hierarchical clustering. It combines mode-hopping Monte Carlo (MC) sampling with an Expectation-Maximization (EM) algorithm for NEMs to increase sensitivity compared to existing methods. A meta-analysis of four Mediator perturbation studies in Saccharomyces cerevisiae, three of which are unpublished, provides new insight into the Mediator signaling network. In addition to the known modular organization of the Mediator subunits, MC EMiNEM reveals a hierarchical ordering of its internal information flow, which is putatively transmitted through structural changes within the complex. We identify the N-terminus of Med7 as a peripheral entity, entailing only local structural changes upon perturbation, while the C-terminus of Med7 and Med19 appear to play a central role. MC EMiNEM associates Mediator subunits to most directly affected genes, which, in conjunction with gene set enrichment analysis, allows us to construct an interaction map of Mediator subunits and transcription factors.

Termination of RNA polymerase II (Pol II) transcription is an important step in the transcription cycle, which involves the dislodgement of polymerase from DNA, leading to release of a functional transcript. Recent studies have identified the key players required for this process and showed that a common feature of these proteins is a conserved domain that interacts with the phosphorylated C-terminus of Pol II (CTD-interacting domain, CID). However, the mechanism by which transcription termination is achieved is not understood. Using genome-wide methods, here we show that the fission yeast CID-protein Seb1 is essential for termination of protein-coding and non-coding genes through interaction with S2-phosphorylated Pol II and nascent RNA. Furthermore, we present the crystal structures of the Seb1 CTD- and RNA-binding modules. Unexpectedly, the latter reveals an intertwined two-domain arrangement of a canonical RRM and second domain. These results provide important insights into the mechanism underlying eukaryotic transcription termination.

In the developing human neocortex, progenitor cells generate diverse cell types prenatally. Progenitor cells and newborn neurons respond to signaling cues, including neurotransmitters. While single-cell RNA sequencing has revealed cellular diversity, physiological heterogeneity has yet to be mapped onto these developing and diverse cell types. By combining measurements of intracellular Ca2+ elevations in response to neurotransmitter receptor agonists and RNA sequencing of the same single cells, we show that Ca2+ responses are cell-type-specific and change dynamically with lineage progression. Physiological response properties predict molecular cell identity and additionally reveal diversity not captured by single-cell transcriptomics. We find that the serotonin receptor HTR2A selectively activates radial glia cells in the developing human, but not mouse, neocortex, and inhibiting HTR2A receptors in human radial glia disrupts the radial glial scaffold. We show highly specific neurotransmitter signaling during neurogenesis in the developing human neocortex and highlight evolutionarily divergent mechanisms of physiological signaling.

The 3'-ends of eukaryotic pre-mRNAs are processed in the nucleus by a large multiprotein complex, the cleavage and polyadenylation factor (CPF). CPF cleaves RNA, adds a poly(A) tail and signals transcription termination. CPF harbors four enzymatic activities essential for these processes, but how these are coordinated remains poorly understood. Several subunits of CPF, including two protein phosphatases, are also found in the related 'associated with Pta1' (APT) complex, but the relationship between CPF and APT is unclear. Here, we show that the APT complex is physically distinct from CPF. The 21 kDa Syc1 protein is associated only with APT, and not with CPF, and is therefore the defining subunit of APT. Using ChIP-seq, PAR-CLIP and RNA-seq, we show that Syc1/APT has distinct, but possibly overlapping, functions from those of CPF. Syc1/APT plays a more important role in sn/snoRNA production whereas CPF processes the 3'-ends of protein-coding pre-mRNAs. These results define distinct protein machineries for synthesis of mature eukaryotic protein-coding and non-coding RNAs.

At the end of protein-coding genes, RNA polymerase (Pol) II undergoes a concerted transition that involves 3'-processing of the pre-mRNA and transcription termination. Here, we present a genome-wide analysis of the 3'-transition in budding yeast. We find that the 3'-transition globally requires the Pol II elongation factor Spt5 and factors involved in the recognition of the polyadenylation (pA) site and in endonucleolytic RNA cleavage. Pol II release from DNA occurs in a narrow termination window downstream of the pA site and requires the "torpedo" exonuclease Rat1 (XRN2 in human). The Rat1-interacting factor Rai1 contributes to RNA degradation downstream of the pA site. Defects in the 3'-transition can result in increased transcription at downstream genes.

RNA splicing is an essential part of eukaryotic gene expression. Although the mechanism of splicing has been extensively studied in vitro, in vivo kinetics for the two-step splicing reaction remain poorly understood. Here, we combine transient transcriptome sequencing (TT-seq) and mathematical modeling to quantify RNA metabolic rates at donor and acceptor splice sites across the human genome. Splicing occurs in the range of minutes and is limited by the speed of RNA polymerase elongation. Splicing kinetics strongly depends on the position and nature of nucleotides flanking splice sites, and on structural interactions between unspliced RNA and small nuclear RNAs in spliceosomal intermediates. Finally, we introduce the 'yield' of splicing as the efficiency of converting unspliced to spliced RNA and show that it is highest for mRNAs and independent of splicing kinetics. These results lead to quantitative models describing how splicing rates and yield are encoded in the human genome.

Eukaryotic genes often generate a variety of RNA isoforms that can lead to functionally distinct protein variants. The synthesis and stability of RNA isoforms is poorly characterized because current methods to quantify RNA metabolism use short-read sequencing and cannot detect RNA isoforms. Here we present nanopore sequencing-based isoform dynamics (nano-ID), a method that detects newly synthesized RNA isoforms and monitors isoform metabolism. Nano-ID combines metabolic RNA labeling, long-read nanopore sequencing of native RNA molecules, and machine learning. Nano-ID derives RNA stability estimates and evaluates stability determining factors such as RNA sequence, poly(A)-tail length, secondary structure, translation efficiency, and RNA-binding proteins. Application of nano-ID to the heat shock response in human cells reveals that many RNA isoforms change their stability. Nano-ID also shows that the metabolism of individual RNA isoforms differs strongly from that estimated for the combined RNA signal at a specific gene locus. Nano-ID enables studies of RNA metabolism at the level of single RNA molecules and isoforms in different cell states and conditions.

For transcription through chromatin, RNA polymerase (Pol) II associates with elongation factors (EFs). Here we show that many EFs crosslink to RNA emerging from transcribing Pol II in the yeast Saccharomyces cerevisiae. Most EFs crosslink preferentially to mRNAs, rather than unstable non-coding RNAs. RNA contributes to chromatin association of many EFs, including the Pol II serine 2 kinases Ctk1 and Bur1 and the histone H3 methyltransferases Set1 and Set2. The Ctk1 kinase complex binds RNA in vitro, consistent with direct EF-RNA interaction. Set1 recruitment to genes in vivo depends on its RNA recognition motifs (RRMs). These results strongly suggest that nascent RNA contributes to EF recruitment to transcribing Pol II. We propose that EF-RNA interactions facilitate assembly of the elongation complex on transcribed genes when RNA emerges from Pol II, and that loss of EF-RNA interactions upon RNA cleavage at the polyadenylation site triggers disassembly of the elongation complex.

Gene transcription can be activated by decreasing the duration of RNA polymerase II pausing in the promoter-proximal region, but how this is achieved remains unclear. Here we use a 'multi-omics' approach to demonstrate that the duration of polymerase pausing generally limits the productive frequency of transcription initiation in human cells ('pause-initiation limit'). We further engineer a human cell line to allow for specific and rapid inhibition of the P-TEFb kinase CDK9, which is implicated in polymerase pause release. CDK9 activity decreases the pause duration but also increases the productive initiation frequency. This shows that CDK9 stimulates release of paused polymerase and activates transcription by increasing the number of transcribing polymerases and thus the amount of mRNA synthesized per time. CDK9 activity is also associated with long-range chromatin interactions, suggesting that enhancers can influence the pause-initiation limit to regulate transcription.

Cryo-electron microscopy (cryo-EM) enables macromolecular structure determination in vitro and inside cells. In addition to aligning individual particles, accurate registration of sample motion and three-dimensional deformation during exposures are crucial for achieving high-resolution reconstructions. Here we describe M, a software tool that establishes a reference-based, multi-particle refinement framework for cryo-EM data and couples a comprehensive spatial deformation model to in silico correction of electron-optical aberrations. M provides a unified optimization framework for both frame-series and tomographic tilt-series data. We show that tilt-series data can provide the same resolution as frame-series data on a purified protein specimen, indicating that the alignment step no longer limits the resolution obtainable from tomographic data. In combination with Warp and RELION, M resolves to residue level a 70S ribosome bound to an antibiotic inside intact bacterial cells. Our work provides a computational tool that facilitates structural biology in cells.

Termination of RNA polymerase II (RNAPII) transcription in metazoans relies largely on the cleavage and polyadenylation (CPA) and integrator (INT) complexes originally found to act at the ends of protein-coding and small nuclear RNA (snRNA) genes, respectively. Here, we monitor CPA- and INT-dependent termination activities genome-wide, including at thousands of previously unannotated transcription units (TUs), producing unstable RNA. We verify the global activity of CPA occurring at pA sites indiscriminately of their positioning relative to the TU promoter. We also identify a global activity of INT, which is largely sequence-independent and restricted to a ~3-kb promoter-proximal region. Our analyses suggest two functions of genome-wide INT activity: it dampens transcriptional output from weak promoters, and it provides quality control of RNAPII complexes that are unfavorably configured for transcriptional elongation. We suggest that the function of INT in stable snRNA production is an exception from its general cellular role, the attenuation of non-productive transcription.