From structure to systems: high-resolution, quantitative genetic analysis of RNA polymerase II

Evolutionary rewiring of the dynamic network underpinning allosteric epistasis in NS1 of influenza A virus

Viral proteins frequently mutate to evade or antagonize host innate immune responses, yet the impact of these mutations on the molecular energy landscape remains unclear. Epistasis, the intramolecular communications between mutations, often renders the combined mutational effects unpredictable. Nonstructural protein 1 (NS1) is a major virulence factor of the influenza A virus (IAV) that activates host PI3K by binding to its p85β subunit. Here, we present the deep analysis for the impact of evolutionary mutations in NS1 that emerged between the 1918 pandemic IAV strain and its descendant PR8 strain. Our analysis reveal how the mutations rewired inter-residue communications which underlies long-range allosteric and epistatic networks in NS1. Our findings show that PR8 NS1 binds to p85β with approximately 10-fold greater affinity than 1918 NS1 due to allosteric mutational effects. Notably, these mutations also exhibited long-range epistatic effects. NMR chemical shift perturbation and methyl-axis order parameter analyses revealed that the mutations induced long-range structural and dynamic changes in PR8 NS1, enhancing its affinity to p85β. Complementary MD simulations and graph-based network analysis uncover how these mutations rewire dynamic residue interaction networks, which underlies the long-range epistasis and allosteric effects on p85β-binding affinity. Significantly, we find that conformational dynamics of residues with high betweenness centrality play a crucial role in communications between network communities and are highly conserved across influenza A virus evolution. These findings advance our mechanistic understanding of the allosteric and epistatic communications between distant residues and provides insight into their role in the molecular evolution of NS1.

Higher-order epistasis within Pol II trigger loop haplotypes

RNA polymerase II (Pol II) has a highly conserved domain, the trigger loop (TL), that controls transcription fidelity and speed. We previously probed pairwise genetic interactions between residues within and surrounding the TL and identified widespread incompatibility between TLs of different species when placed in the Saccharomyces cerevisiae Pol II context, indicating epistasis between the TL and its surrounding context. We sought to understand the nature of this incompatibility and probe higher order epistasis internal to the TL. We have employed deep mutational scanning with selected natural TL variants ("haplotypes"), and all possible intermediate substitution combinations between them and the yeast Pol II TL. We identified both positive and negative higher-order residue interactions within example TL haplotypes. Intricate higher-order epistasis formed by TL residues was sometimes only apparent from analysis of intermediate genotypes, emphasizing complexity of epistatic interactions. Furthermore, we distinguished TL substitutions with distinct classes of epistatic patterns, suggesting specific TL residues that potentially influence TL evolution. Our examples of complex residue interactions suggest possible pathways for epistasis to facilitate Pol II evolution.

mRNA initiation and termination are spatially coordinated

The expression of a precise mRNA transcriptome is crucial for establishing cell identity and function, with dozens of alternative isoforms produced for a single gene sequence. The regulation of mRNA isoform usage occurs by the coordination of co-transcriptional mRNA processing mechanisms across a gene. Decisions involved in mRNA initiation and termination underlie the largest extent of mRNA isoform diversity, but little is known about any relationships between decisions at both ends of mRNA molecules. Here, we systematically profile the joint usage of mRNA transcription start sites (TSSs) and polyadenylation sites (PASs) across tissues and species. Using both short and long read RNA-seq data, we observe that mRNAs preferentially using upstream TSSs also tend to use upstream PASs, and congruently, the usage of downstream sites is similarly paired. This observation suggests that mRNA 5' end choice may directly influence mRNA 3' ends. Our results suggest a novel "Positional Initiation-Termination Axis" (PITA), in which the usage of alternative terminal sites are coupled based on the order in which they appear in the genome. PITA isoforms are more likely to encode alternative protein domains and use conserved sites. PITA is strongly associated with the length of genomic features, such that PITA is enriched in longer genes with more area devoted to regions that regulate alternative 5' or 3' ends. Strikingly, we found that PITA genes are more likely than non-PITA genes to have multiple, overlapping chromatin structural domains related to pairing of ordinally coupled start and end sites. In turn, PITA coupling is also associated with fast RNA Polymerase II (RNAPII) trafficking across these long gene regions. Our findings indicate that a combination of spatial and kinetic mechanisms couple transcription initiation and mRNA 3' end decisions based on ordinal position to define the expression mRNA isoforms.

Quantitative analysis of transcription start site selection reveals control by DNA sequence, RNA polymerase II activity and NTP levels

Transcription start site (TSS) selection is a key step in gene expression and occurs at many promoter positions over a wide range of efficiencies. Here we develop a massively parallel reporter assay to quantitatively dissect contributions of promoter sequence, nucleoside triphosphate substrate levels and RNA polymerase II (Pol II) activity to TSS selection by 'promoter scanning' in Saccharomyces cerevisiae (Pol II MAssively Systematic Transcript End Readout, 'Pol II MASTER'). Using Pol II MASTER, we measure the efficiency of Pol II initiation at 1,000,000 individual TSS sequences in a defined promoter context. Pol II MASTER confirms proposed critical qualities of S. cerevisiae TSS -8, -1 and +1 positions, quantitatively, in a controlled promoter context. Pol II MASTER extends quantitative analysis to surrounding sequences and determines that they tune initiation over a wide range of efficiencies. These results enabled the development of a predictive model for initiation efficiency based on sequence. We show that genetic perturbation of Pol II catalytic activity alters initiation efficiency mostly independently of TSS sequence, but selectively modulates preference for the initiating nucleotide. Intriguingly, we find that Pol II initiation efficiency is directly sensitive to guanosine-5'-triphosphate levels at the first five transcript positions and to cytosine-5'-triphosphate and uridine-5'-triphosphate levels at the second position genome wide. These results suggest individual nucleoside triphosphate levels can have transcript-specific effects on initiation, representing a cryptic layer of potential regulation at the level of Pol II biochemical properties. The results establish Pol II MASTER as a method for quantitative dissection of transcription initiation in eukaryotes.

Structural basis of transcription: RNA Polymerase II substrate binding and metal coordination at 3.0 Å using a free-electron laser

Catalysis and translocation of multi-subunit DNA-directed RNA polymerases underlie all cellular mRNA synthesis. RNA polymerase II (Pol II) synthesizes eukaryotic pre-mRNAs from a DNA template strand buried in its active site. Structural details of catalysis at near atomic resolution and precise arrangement of key active site components have been elusive. Here we present the free electron laser (FEL) structure of a matched ATP-bound Pol II, revealing the full active site interaction network at the highest resolution to date, including the trigger loop (TL) in the closed conformation, bonafide occupancy of both site A and B Mg2+, and a putative third (site C) Mg2+ analogous to that described for some DNA polymerases but not observed previously for cellular RNA polymerases. Molecular dynamics (MD) simulations of the structure indicate that the third Mg2+ is coordinated and stabilized at its observed position. TL residues provide half of the substrate binding pocket while multiple TL/bridge helix (BH) interactions induce conformational changes that could propel translocation upon substrate hydrolysis. Consistent with TL/BH communication, a FEL structure and MD simulations of the hyperactive Rpb1 T834P bridge helix mutant reveals rearrangement of some active site interactions supporting potential plasticity in active site function and long-distance effects on both the width of the central channel and TL conformation, likely underlying its increased elongation rate at the expense of fidelity.

Pre-mRNA splicing and its cotranscriptional connections

Transcription of eukaryotic genes by RNA polymerase II (Pol II) yields RNA precursors containing introns that must be spliced out and the flanking exons ligated together. Splicing is catalyzed by a dynamic ribonucleoprotein complex called the spliceosome. Recent evidence has shown that a large fraction of splicing occurs cotranscriptionally as the RNA chain is extruded from Pol II at speeds of up to 5 kb/minute. Splicing is more efficient when it is tethered to the transcription elongation complex, and this linkage permits functional coupling of splicing with transcription. We discuss recent progress that has uncovered a network of connections that link splicing to transcript elongation and other cotranscriptional RNA processing events.

Cotranscriptional RNA processing and modification in plants

The activities of RNA polymerases shape the epigenetic landscape of genomes with profound consequences for genome integrity and gene expression. A fundamental event during the regulation of eukaryotic gene expression is the coordination between transcription and RNA processing. Most primary RNAs mature through various RNA processing and modification events to become fully functional. While pioneering results positioned RNA maturation steps after transcription ends, the coupling between the maturation of diverse RNA species and their transcription is becoming increasingly evident in plants. In this review, we discuss recent advances in our understanding of the crosstalk between RNA Polymerase II, IV, and V transcription and nascent RNA processing of both coding and noncoding RNAs.

Widespread epistasis shapes RNA Polymerase II active site function and evolution

Multi-subunit RNA Polymerases (msRNAPs) are responsible for transcription in all kingdoms of life. At the heart of these msRNAPs is an ultra-conserved active site domain, the trigger loop (TL), coordinating transcription speed and fidelity by critical conformational changes impacting multiple steps in substrate selection, catalysis, and translocation. Previous studies have observed several different types of genetic interactions between eukaryotic RNA polymerase II (Pol II) TL residues, suggesting that the TL's function is shaped by functional interactions of residues within and around the TL. The extent of these interaction networks and how they control msRNAP function and evolution remain to be determined. Here we have dissected the Pol II TL interaction landscape by deep mutational scanning in Saccharomyces cerevisiae Pol II. Through analysis of over 15000 alleles, representing all single mutants, a rationally designed subset of double mutants, and evolutionarily observed TL haplotypes, we identify interaction networks controlling TL function. Substituting residues creates allele-specific networks and propagates epistatic effects across the Pol II active site. Furthermore, the interaction landscape further distinguishes alleles with similar growth phenotypes, suggesting increased resolution over the previously reported single mutant phenotypic landscape. Finally, co-evolutionary analyses reveal groups of co-evolving residues across Pol II converge onto the active site, where evolutionary constraints interface with pervasive epistasis. Our studies provide a powerful system to understand the plasticity of RNA polymerase mechanism and evolution, and provide the first example of pervasive epistatic landscape in a highly conserved and constrained domain within an essential enzyme.

Characterization of RNA polymerase II trigger loop mutations using molecular dynamics simulations and machine learning

Catalysis and fidelity of multisubunit RNA polymerases rely on a highly conserved active site domain called the trigger loop (TL), which achieves roles in transcription through conformational changes and interaction with NTP substrates. The mutations of TL residues cause distinct effects on catalysis including hypo- and hyperactivity and altered fidelity. We applied molecular dynamics simulation (MD) and machine learning (ML) techniques to characterize TL mutations in the Saccharomyces cerevisiae RNA Polymerase II (Pol II) system. We did so to determine relationships between individual mutations and phenotypes and to associate phenotypes with MD simulated structural alterations. Using fitness values of mutants under various stress conditions, we modeled phenotypes along a spectrum of continual values. We found that ML could predict the phenotypes with 0.68 R2 correlation from amino acid sequences alone. It was more difficult to incorporate MD data to improve predictions from machine learning, presumably because MD data is too noisy and possibly incomplete to directly infer functional phenotypes. However, a variational auto-encoder model based on the MD data allowed the clustering of mutants with different phenotypes based on structural details. Overall, we found that a subset of loss-of-function (LOF) and lethal mutations tended to increase distances of TL residues to the NTP substrate, while another subset of LOF and lethal substitutions tended to confer an increase in distances between TL and bridge helix (BH). In contrast, some of the gain-of-function (GOF) mutants appear to cause disruption of hydrophobic contacts among TL and nearby helices.

Nucleotide addition and cleavage by RNA polymerase II: Coordination of two catalytic reactions using a single active site

RNA polymerase II (Pol II) incorporates complementary ribonucleotides into the growing RNA chain one at a time via the nucleotide addition cycle. The nucleotide addition cycle, however, is prone to misincorporation of noncomplementary nucleotides. Thus, to ensure transcriptional fidelity, Pol II backtracks and then cleaves the misincorporated nucleotides. These two reverse reactions, nucleotide addition and cleavage, are catalyzed in the same active site of Pol II, which is different from DNA polymerases or other endonucleases. Recently, substantial progress has been made to understand how Pol II effectively performs its dual role in the same active site. Our review highlights these recent studies and provides an overall model of the catalytic mechanisms of Pol II. In particular, RNA extension follows the two-metal-ion mechanism, and several Pol II residues play important roles to facilitate the catalysis. In sharp contrast, the cleavage reaction is independent of any Pol II residues. Interestingly, Pol II relies on its residues to recognize the misincorporated nucleotides during the backtracking process, prior to cleavage. In this way, Pol II efficiently compartmentalizes its two distinct catalytic functions using the same active site. Lastly, we also discuss a new perspective on the potential third Mg²⁺ in the nucleotide addition and intrinsic cleavage reactions.

Rate of transcription elongation and sequence-specific pausing by RNA polymerase I directly influence rRNA processing

One of the first steps in ribosome biogenesis is transcription of the ribosomal DNA by RNA polymerase I (Pol I). Processing of the resultant rRNA begins cotranscriptionally, and perturbation of Pol I transcription elongation results in defective rRNA processing. Mechanistic insight regarding the link between transcription elongation and ribosome assembly is lacking because of limited in vivo methods to assay Pol I transcription. Here, we use native elongating transcript sequencing (NET-Seq) with a strain of Saccharomyces cerevisiae containing a point mutation in Pol I, rpa190-F1205H, which results in impaired rRNA processing and ribosome assembly. We previously demonstrated that this mutation caused a mild reduction in the transcription elongation rate of Pol I in vitro; however, transcription elongation by the mutant has not been characterized in vivo. Here, our findings demonstrate that the mutant Pol I has an increased pause propensity during processive transcription elongation both in vitro and in vivo. NET-Seq reveals that rpa190-F1205H Pol I displays alternative pause site preferences in vivo. Specifically, the mutant is sensitized to A/G residues in the RNA:DNA hybrid and at the last incorporated nucleotide position. Furthermore, both NET-Seq and EM analysis of Miller chromatin spreads reveal pileups of rpa190-F1205H Pol I throughout the ribosomal DNA, particularly at the 5' end of the 35S gene. This combination of in vitro and in vivo analyses of a Pol I mutant provides novel insights into Pol I elongation properties and indicates how these properties are crucial for efficient cotranscriptional rRNA processing and ribosome assembly.

The Role of S. cerevisiae Sub1/PC4 in Transcription Elongation Depends on the C-Terminal Region and Is Independent of the ssDNA Binding Domain

Saccharomyces cerevisiae Sub1 (ScSub1) has been defined as a transcriptional stimulatory protein due to its homology to the ssDNA binding domain (ssDBD) of human PC4 (hPC4). Recently, PC4/Sub1 orthologues have been elucidated in eukaryotes, prokaryotes, and bacteriophages with functions related to DNA metabolism. Additionally, ScSub1 contains a unique carboxyl–terminal region (CT) of unknown function up to date. Specifically, it has been shown that Sub1 is required for transcription activation, as well as other processes, throughout the transcription cycle. Despite the progress that has been made in understanding the mechanism underlying Sub1′s functions, some questions remain unanswered. As a case in point: whether Sub1’s roles in initiation and elongation are differentially predicated on distinct regions of the protein or how Sub1′s functions are regulated. Here, we uncover some residues that are key for DNA–ScSub1 interaction in vivo, localized in the ssDBD, and required for Sub1 recruitment to promoters. Furthermore, using an array of genetic and molecular techniques, we demonstrate that the CT region is required for transcription elongation by RNA polymerase II (RNAPII). Altogether, our data indicate that Sub1 plays a dual role during transcription—in initiation through the ssDBD and in elongation through the CT region.

RNA helicase-dependent gene looping impacts messenger RNA processing

DDX5 and DDX17 are DEAD-box RNA helicase paralogs which regulate several aspects of gene expression, especially transcription and splicing, through incompletely understood mechanisms. A transcriptome analysis of DDX5/DDX17-depleted human cells confirmed the large impact of these RNA helicases on splicing and revealed a widespread deregulation of 3' end processing. In silico analyses and experiments in cultured cells showed the binding and functional contribution of the genome organizing factor CTCF to chromatin sites at or near a subset of DDX5/DDX17-dependent exons that are characterized by a high GC content and a high density of RNA Polymerase II. We propose the existence of an RNA helicase-dependent relationship between CTCF and the dynamics of transcription across DNA and/or RNA structured regions, that contributes to the processing of internal and terminal exons. Moreover, local DDX5/DDX17-dependent chromatin loops spatially connect RNA helicase-regulated exons with their cognate promoter, and we provide the first direct evidence that de novo gene looping modifies alternative splicing and polyadenylation. Overall our findings uncover the impact of DDX5/DDX17-dependent chromatin folding on pre-messenger RNA processing.

Genetic dissection of the RNA polymerase II transcription cycle

How DNA sequence affects the dynamics and position of RNA Polymerase II (Pol II) during transcription remains poorly understood. Here, we used naturally occurring genetic variation in F1 hybrid mice to explore how DNA sequence differences affect the genome-wide distribution of Pol II. We measured the position and orientation of Pol II in eight organs collected from heterozygous F1 hybrid mice using ChRO-seq. Our data revealed a strong genetic basis for the precise coordinates of transcription initiation and promoter proximal pause, allowing us to redefine molecular models of core transcriptional processes. Our results implicate DNA sequence, including both known and novel DNA sequence motifs, as key determinants of the position of Pol II initiation and pause. We report evidence that initiation site selection follows a stochastic process similar to Brownian motion along the DNA template. We found widespread differences in the position of transcription termination, which impact the primary structure and stability of mature mRNA. Finally, we report evidence that allelic changes in transcription often affect mRNA and ncRNA expression across broad genomic domains. Collectively, we reveal how DNA sequences shape core transcriptional processes at single nucleotide resolution in mammals.

From systems to structure - using genetic data to model protein structures

Understanding the effects of genetic variation is a fundamental problem in biology that requires methods to analyse both physical and functional consequences of sequence changes at systems-wide and mechanistic scales. To achieve a systems view, protein interaction networks map which proteins physically interact, while genetic interaction networks inform on the phenotypic consequences of perturbing these protein interactions. Until recently, understanding the molecular mechanisms that underlie these interactions often required biophysical methods to determine the structures of the proteins involved. The past decade has seen the emergence of new approaches based on coevolution, deep mutational scanning and genome-scale genetic or chemical-genetic interaction mapping that enable modelling of the structures of individual proteins or protein complexes. Here, we review the emerging use of large-scale genetic datasets and deep learning approaches to model protein structures and their interactions, and discuss the integration of structural data from different sources.

Dysregulation and therapeutic targeting of RNA splicing in cancer

High-throughput sequencing and functional characterization of the cancer transcriptome have uncovered cancer-specific dysregulation of RNA splicing across a variety of cancers. Alterations in the cancer genome and dysregulation of RNA splicing factors lead to missplicing, splicing alteration-dependent gene expression and, in some cases, generation of novel splicing-derived proteins. Here, we review recent advances in our understanding of aberrant splicing in cancer pathogenesis and present strategies to harness cancer-specific aberrant splicing for therapeutic intent.

Integrative structure determination of histones H3 and H4 using genetic interactions

Integrative structure modeling is increasingly used for determining the architectures of biological assemblies, especially those that are structurally heterogeneous. Recently, we reported on how to convert in vivo genetic interaction measurements into spatial restraints for structural modeling: first, phenotypic profiles are generated for each point mutation and thousands of gene deletions or environmental perturbations. Following, the phenotypic profile similarities are converted into distance restraints on the pairs of mutated residues. We illustrate the approach by determining the structure of the histone H3-H4 complex. The method is implemented in our open-source IMP program, expanding the structural biology toolbox by allowing structural characterization based on in vivo data without the need to purify the target system. We compare genetic interaction measurements to other sources of structural information, such as residue coevolution and deep-learning structure prediction of complex subunits. We also suggest that determining genetic interactions could benefit from new technologies, such as CRISPR-Cas9 approaches to gene editing, especially for mammalian cells. Finally, we highlight the opportunity for using genetic interactions to determine recalcitrant biomolecular structures, such as those of disordered proteins, transient protein assemblies, and host-pathogen protein complexes.

High-throughput functional characterization of protein phosphorylation sites in yeast

Phosphorylation is a critical post-translational modification involved in the regulation of almost all cellular processes. However, fewer than 5% of thousands of recently discovered phosphosites have been functionally annotated. In this study, we devised a chemical genetic approach to study the functional relevance of phosphosites in Saccharomyces cerevisiae. We generated 474 yeast strains with mutations in specific phosphosites that were screened for fitness in 102 conditions, along with a gene deletion library. Of these phosphosites, 42% exhibited growth phenotypes, suggesting that these are more likely functional. We inferred their function based on the similarity of their growth profiles with that of gene deletions and validated a subset by thermal proteome profiling and lipidomics. A high fraction exhibited phenotypes not seen in the corresponding gene deletion, suggestive of a gain-of-function effect. For phosphosites conserved in humans, the severity of the yeast phenotypes is indicative of their human functional relevance. This high-throughput approach allows for functionally characterizing individual phosphosites at scale.

Structural visualization of de novo transcription initiation by Saccharomyces cerevisiae RNA polymerase II

Previous structural studies of the initiation-elongation transition of RNA polymerase II (pol II) transcription have relied on the use of synthetic oligonucleotides, often artificially discontinuous to capture pol II in the initiating state. Here, we report multiple structures of initiation complexes converted de novo from a 33-subunit yeast pre-initiation complex (PIC) through catalytic activities and subsequently stalled at different template positions. We determine that PICs in the initially transcribing complex (ITC) can synthesize a transcript of ∼26 nucleotides before transitioning to an elongation complex (EC) as determined by the loss of general transcription factors (GTFs). Unexpectedly, transition to an EC was greatly accelerated when an ITC encountered a downstream EC stalled at promoter proximal regions and resulted in a collided head-to-end dimeric EC complex. Our structural analysis reveals a dynamic state of TFIIH, the largest of GTFs, in PIC/ITC with distinct functional consequences at multiple steps on the pathway to elongation.

Systems-level effects of allosteric perturbations to a model molecular switch

Molecular switch proteins whose cycling between states is controlled by opposing regulators^1,2 are central to biological signal transduction. As switch proteins function within highly connected interaction networks³, the fundamental question arises of how functional specificity is achieved when different processes share common regulators. Here we show that functional specificity of the small GTPase switch protein Gsp1 in Saccharomyces cerevisiae (the homologue of the human protein RAN)⁴ is linked to differential sensitivity of biological processes to different kinetics of the Gsp1 (RAN) switch cycle. We make 55 targeted point mutations to individual protein interaction interfaces of Gsp1 (RAN) and show through quantitative genetic⁵ and physical interaction mapping that Gsp1 (RAN) interface perturbations have widespread cellular consequences. Contrary to expectation, the cellular effects of the interface mutations group by their biophysical effects on kinetic parameters of the GTPase switch cycle and not by the targeted interfaces. Instead, we show that interface mutations allosterically tune the GTPase cycle kinetics. These results suggest a model in which protein partner binding, or post-translational modifications at distal sites, could act as allosteric regulators of GTPase switching. Similar mechanisms may underlie regulation by other GTPases, and other biological switches. Furthermore, our integrative platform to determine the quantitative consequences of molecular perturbations may help to explain the effects of disease mutations that target central molecular switches.

Ssl2/TFIIH function in transcription start site scanning by RNA polymerase II in Saccharomyces cerevisiae

In Saccharomyces cerevisiae, RNA polymerase II (Pol II) selects transcription start sites (TSSs) by a unidirectional scanning process. During scanning, a preinitiation complex (PIC) assembled at an upstream core promoter initiates at select positions within a window ~40-120 bp downstream. Several lines of evidence indicate that Ssl2, the yeast homolog of XPB and an essential and conserved subunit of the general transcription factor (GTF) TFIIH, drives scanning through its DNA-dependent ATPase activity, therefore potentially controlling both scanning rate and scanning extent (processivity). To address questions of how Ssl2 functions in promoter scanning and interacts with other initiation activities, we leveraged distinct initiation-sensitive reporters to identify novel ssl2 alleles. These ssl2 alleles, many of which alter residues conserved from yeast to human, confer either upstream or downstream TSS shifts at the model promoter ADH1 and genome-wide. Specifically, tested ssl2 alleles alter TSS selection by increasing or narrowing the distribution of TSSs used at individual promoters. Genetic interactions of ssl2 alleles with other initiation factors are consistent with ssl2 allele classes functioning through increasing or decreasing scanning processivity but not necessarily scanning rate. These alleles underpin a residue interaction network that likely modulates Ssl2 activity and TFIIH function in promoter scanning. We propose that the outcome of promoter scanning is determined by two functional networks, the first being Pol II activity and factors that modulate it to determine initiation efficiency within a scanning window, and the second being Ssl2/TFIIH and factors that modulate scanning processivity to determine the width of the scanning widow.

Characterization of an A3G-VifHIV-1-CRL5-CBFβ Structure Using a Cross-linking Mass Spectrometry Pipeline for Integrative Modeling of Host-Pathogen Complexes

Structural analysis of host-pathogen protein complexes remains challenging, largely due to their structural heterogeneity. Here, we describe a pipeline for the structural characterization of these complexes using integrative structure modeling based on chemical cross-links and residue-protein contacts inferred from mutagenesis studies. We used this approach on the HIV-1 Vif protein bound to restriction factor APOBEC3G (A3G), the Cullin-5 E3 ring ligase (CRL5), and the cellular transcription factor Core Binding Factor Beta (CBFβ) to determine the structure of the (A3G-Vif-CRL5-CBFβ) complex. Using the MS-cleavable DSSO cross-linker to obtain a set of 132 cross-links within this reconstituted complex along with the atomic structures of the subunits and mutagenesis data, we computed an integrative structure model of the heptameric A3G-Vif-CRL5-CBFβ complex. The structure, which was validated using a series of tests, reveals that A3G is bound to Vif mostly through its N-terminal domain. Moreover, the model ensemble quantifies the dynamic heterogeneity of the A3G C-terminal domain and Cul5 positions. Finally, the model was used to rationalize previous structural, mutagenesis and functional data not used for modeling, including information related to the A3G-bound and unbound structures as well as mapping functional mutations to the A3G-Vif interface. The experimental and computational approach described here is generally applicable to other challenging host-pathogen protein complexes.

RNA polymerase II speed: a key player in controlling and adapting transcriptome composition

RNA polymerase II (RNA Pol II) speed or elongation rate, i.e., the number of nucleotides synthesized per unit of time, is a major determinant of transcriptome composition. It controls co‐transcriptional processes such as splicing, polyadenylation, and transcription termination, thus regulating the production of alternative splice variants, circular RNAs, alternatively polyadenylated transcripts, or read‐through transcripts. RNA Pol II speed itself is regulated in response to intra‐ and extra‐cellular stimuli and can in turn affect the transcriptome composition in response to these stimuli. Evidence points to a potentially important role of transcriptome composition modification through RNA Pol II speed regulation for adaptation of cells to a changing environment, thus pointing to a function of RNA Pol II speed regulation in cellular physiology. Analyzing RNA Pol II speed dynamics may therefore be central to fully understand the regulation of physiological processes, such as the development of multicellular organisms. Recent findings also raise the possibility that RNA Pol II speed deregulation can be detrimental and participate in disease progression. Here, we review initial and current approaches to measure RNA Pol II speed, as well as providing an overview of the factors controlling speed and the co‐transcriptional processes which are affected. Finally, we discuss the role of RNA Pol II speed regulation in cell physiology.

A termination-independent role of Rat1 in cotranscriptional splicing

Rat1 is a 5′→3′ exoribonuclease in budding yeast. It is a highly conserved protein with homologs being present in fission yeast, flies, worms, mice and humans. Rat1 and its human homolog Xrn2 have been implicated in multiple nuclear processes. Here we report a novel role of Rat1 in mRNA splicing. We observed an increase in the level of unspliced transcripts in mutants of Rat1. Accumulation of unspliced transcripts was not due to the surveillance role of Rat1 in degrading unspliced mRNA, or an indirect effect of Rat1 function in termination of transcription or on the level of splicing factors in the cell, or due to an increased elongation rate in Rat1 mutants. ChIP-Seq analysis revealed Rat1 crosslinking to the introns of a subset of yeast genes. Mass spectrometry and coimmunoprecipitation revealed an interaction of Rat1 with the Clf1, Isy1, Yju2, Prp43 and Sub2 splicing factors. Furthermore, recruitment of splicing factors on the intron was compromised in the Rat1 mutant. Based on these findings we propose that Rat1 has a novel role in splicing of mRNA in budding yeast. Rat1, however, is not a general splicing factor as it crosslinks to only long introns with an average length of 400 nucleotides.

Chemical-genetic interrogation of RNA polymerase mutants reveals structure-function relationships and physiological tradeoffs

The multi-subunit bacterial RNA polymerase (RNAP) and its associated regulators carry out transcription and integrate myriad regulatory signals. Numerous studies have interrogated RNAP mechanism, and RNAP mutations drive Escherichia coli adaptation to many health- and industry-relevant environments, yet a paucity of systematic analyses hampers our understanding of the fitness trade-offs from altering RNAP function. Here, we conduct a chemical-genetic analysis of a library of RNAP mutants. We discover phenotypes for non-essential insertions, show that clustering mutant phenotypes increases their predictive power for drawing functional inferences, and demonstrate that some RNA polymerase mutants both decrease average cell length and prevent killing by cell-wall targeting antibiotics. Our findings demonstrate that RNAP chemical-genetic interactions provide a general platform for interrogating structure-function relationships in vivo and for identifying physiological trade-offs of mutations, including those relevant for disease and biotechnology. This strategy should have broad utility for illuminating the role of other important protein complexes.

Co-transcriptional splicing efficiencies differ within genes and between cell types

Pre-mRNA splicing is carried out by the spliceosome and involves splice site recognition, removal of introns, and ligation of exons. Components of the spliceosome have been shown to interact with the elongating RNA polymerase II (RNAPII) which is thought to allow splicing to occur concurrently with transcription. However, little is known about the regulation and efficiency of co-transcriptional splicing in human cells. In this study, we used Bru-seq and BruChase-seq to determine the co-transcriptional splicing efficiencies of 17,000 introns expressed across 6 human cell lines. We found that less than half of all introns across these 6 cell lines were co-transcriptionally spliced. Splicing efficiencies for individual introns showed variations across cell lines, suggesting that splicing may be regulated in a cell-type specific manner. Moreover, the splicing efficiency of introns varied within genes. The efficiency of co-transcriptional splicing did not correlate with gene length, intron position, splice site strengths, or the intron/neighboring exons GC content. However, we identified binding signals from multiple RNA binding proteins (RBPs) that correlated with splicing efficiency, including core spliceosomal machinery components-such as SF3B4, U2AF1 and U2AF2 showing higher binding signals in poorly spliced introns. In addition, multiple RBPs, such as BUD13, PUM1 and SND1, showed preferential binding in exons that flank introns with high splicing efficiencies. The nascent RNA splicing patterns presented here across multiple cell types add to our understanding of the complexity in RNA splicing, wherein RNA-binding proteins may play important roles in determining splicing outcomes in a cell type- and intron-specific manner.

Transcription Regulation Through Nascent RNA Folding

Folding of RNA into secondary structures through intramolecular base pairing determines an RNA's three-dimensional architecture and associated function. Simple RNA structures like stem loops can provide specialized functions independent of coding capacity, such as protein binding, regulation of RNA processing and stability, stimulation or inhibition of translation. RNA catalysis is dependent on tertiary structures found in the ribosome, tRNAs and group I and II introns. While the extent to which non-coding RNAs contribute to cellular maintenance is generally appreciated, the fact that both non-coding and coding RNA can assume relevant structural states has only recently gained attention. In particular, the co-transcriptional folding of nascent RNA of all classes has the potential to regulate co-transcriptional processing, RNP (ribonucleoprotein particle) formation, and transcription itself. Riboswitches are established examples of co-transcriptionally folded coding RNAs that directly regulate transcription, mainly in prokaryotes. Here we discuss recent studies in both prokaryotes and eukaryotes showing that structure formation may carry a more widespread regulatory logic during RNA synthesis. Local structures forming close to the catalytic center of RNA polymerases have the potential to regulate transcription by reducing backtracking. In addition, stem loops or more complex structures may alter co-transcriptional RNA processing or its efficiency. Several examples of functional structures have been identified to date, and this review provides an overview of physiologically distinct processes where co-transcriptionally folded RNA plays a role. Experimental approaches such as single-molecule FRET and in vivo structural probing to further advance our insight into the significance of co-transcriptional structure formation are discussed.

Efficient RNA polymerase II pause release requires U2 snRNP function

Transcription by RNA polymerase II (Pol II) is coupled to pre-mRNA splicing, but the underlying mechanisms remain poorly understood. Co-transcriptional splicing requires assembly of a functional spliceosome on nascent pre-mRNA, but whether and how this influences Pol II transcription remains unclear. Here we show that inhibition of pre-mRNA branch site recognition by the spliceosome component U2 snRNP leads to a widespread and strong decrease in new RNA synthesis from human genes. Multiomics analysis reveals that inhibition of U2 snRNP function increases the duration of Pol II pausing in the promoter-proximal region, impairs recruitment of the pause release factor P-TEFb, and reduces Pol II elongation velocity at the beginning of genes. Our results indicate that efficient release of paused Pol II into active transcription elongation requires the formation of functional spliceosomes and that eukaryotic mRNA biogenesis relies on positive feedback from the splicing machinery to the transcription machinery.

Role of promoters in regulating alternative splicing

Alternative splicing (AS) plays a critical role in enhancing proteome complexity in higher eukaryotes. Almost all the multi intron-containing genes undergo AS in humans. Splicing mainly occurs co-transcriptionally, where RNA polymerase II (RNA pol II) plays a crucial role in coordinating transcription and pre-mRNA splicing. Aberrant AS leads to non-functional proteins causative in various pathophysiological conditions such as cancers, neurodegenerative diseases, and muscular dystrophies. Transcription and pre-mRNA splicing are deeply interconnected and can influence each other's functions. Several studies evinced that specific promoters employed by RNA pol II dictate the RNA processing decisions. Promoter-specific recruitment of certain transcriptional factors or transcriptional coactivators influences splicing, and the extent to which these factors affect splicing has not been discussed in detail. Here, in this review, various DNA-binding proteins and their influence on promoter-specific AS are extensively discussed. Besides, this review highlights how the promoter-specific epigenetic changes might regulate AS.

Alternative RNA structures formed during transcription depend on elongation rate and modify RNA processing

Most RNA processing occurs co-transcriptionally. We interrogated nascent pol II transcripts by chemical and enzymatic probing and determined how the "nascent RNA structureome" relates to splicing, A-I editing and transcription speed. RNA folding within introns and steep structural transitions at splice sites are associated with efficient co-transcriptional splicing. A slow pol II mutant elicits extensive remodeling into more folded conformations with increased A-I editing. Introns that become more structured at their 3' splice sites get co-transcriptionally excised more efficiently. Slow pol II altered folding of intronic Alu elements where cryptic splicing and intron retention are stimulated, an outcome mimicked by UV, which decelerates transcription. Slow transcription also remodeled RNA folding around alternative exons in distinct ways that predict whether skipping or inclusion is favored, even though it occurs post-transcriptionally. Hence, co-transcriptional RNA folding modulates post-transcriptional alternative splicing. In summary, the plasticity of nascent transcripts has widespread effects on RNA processing.

Co-transcriptional splicing regulates 3' end cleavage during mammalian erythropoiesis

Pre-mRNA processing steps are tightly coordinated with transcription in many organisms. To determine how co-transcriptional splicing is integrated with transcription elongation and 3' end formation in mammalian cells, we performed long-read sequencing of individual nascent RNAs and precision run-on sequencing (PRO-seq) during mouse erythropoiesis. Splicing was not accompanied by transcriptional pausing and was detected when RNA polymerase II (Pol II) was within 75-300 nucleotides of 3' splice sites (3'SSs), often during transcription of the downstream exon. Interestingly, several hundred introns displayed abundant splicing intermediates, suggesting that splicing delays can take place between the two catalytic steps. Overall, splicing efficiencies were correlated among introns within the same transcript, and intron retention was associated with inefficient 3' end cleavage. Remarkably, a thalassemia patient-derived mutation introducing a cryptic 3'SS improved both splicing and 3' end cleavage of individual β-globin transcripts, demonstrating functional coupling between the two co-transcriptional processes as a determinant of productive gene output.

Genetic interaction mapping informs integrative structure determination of protein complexes

Determining structures of protein complexes is crucial for understanding cellular functions. Here, we describe an integrative structure determination approach that relies on in vivo measurements of genetic interactions. We construct phenotypic profiles for point mutations crossed against gene deletions or exposed to environmental perturbations, followed by converting similarities between two profiles into an upper bound on the distance between the mutated residues. We determine the structure of the yeast histone H3-H4 complex based on ~500,000 genetic interactions of 350 mutants. We then apply the method to subunits Rpb1-Rpb2 of yeast RNA polymerase II and subunits RpoB-RpoC of bacterial RNA polymerase. The accuracy is comparable to that based on chemical cross-links; using restraints from both genetic interactions and cross-links further improves model accuracy and precision. The approach provides an efficient means to augment integrative structure determination with in vivo observations.

Transcription and splicing: A two-way street

RNA synthesis by RNA polymerase II and RNA processing are closely coupled during the transcription cycle of protein-coding genes. This coupling affords opportunities for quality control and regulation of gene expression and the effects can go in both directions. For example, polymerase speed can affect splice site selection and splicing can increase transcription and affect the chromatin landscape. Here we review the many ways that transcription and splicing influence one another, including how splicing "talks back" to transcription. We will also place the connections between transcription and splicing in the context of other RNA processing events that define the exons that will make up the final mRNA. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing.

The transcriptional elongation rate regulates alternative polyadenylation in yeast

Yeast cells undergoing the diauxic response show a striking upstream shift in poly(A) site utilization, with increased use of ORF-proximal poly(A) sites resulting in shorter 3' mRNA isoforms for most genes. This altered poly(A) pattern is extremely similar to that observed in cells containing Pol II derivatives with slow elongation rates. Conversely, cells containing derivatives with fast elongation rates show a subtle downstream shift in poly(A) sites. Polyadenylation patterns of many genes are sensitive to both fast and slow elongation rates, and a global shift of poly(A) utilization is strongly linked to increased purine content of sequences flanking poly(A) sites. Pol II processivity is impaired in diauxic cells, but strains with reduced processivity and normal Pol II elongation rates have normal polyadenylation profiles. Thus, Pol II elongation speed is important for poly(A) site selection and for regulating poly(A) patterns in response to environmental conditions.

Universal promoter scanning by Pol II during transcription initiation in Saccharomyces cerevisiae

BACKGROUND

The majority of eukaryotic promoters utilize multiple transcription start sites (TSSs). How multiple TSSs are specified at individual promoters across eukaryotes is not understood for most species. In Saccharomyces cerevisiae, a pre-initiation complex (PIC) comprised of Pol II and conserved general transcription factors (GTFs) assembles and opens DNA upstream of TSSs. Evidence from model promoters indicates that the PIC scans from upstream to downstream to identify TSSs. Prior results suggest that TSS distributions at promoters where scanning occurs shift in a polar fashion upon alteration in Pol II catalytic activity or GTF function.

RESULTS

To determine the extent of promoter scanning across promoter classes in S. cerevisiae, we perturb Pol II catalytic activity and GTF function and analyze their effects on TSS usage genome-wide. We find that alterations to Pol II, TFIIB, or TFIIF function widely alter the initiation landscape consistent with promoter scanning operating at all yeast promoters, regardless of promoter class. Promoter architecture, however, can determine the extent of promoter sensitivity to altered Pol II activity in ways that are predicted by a scanning model.

CONCLUSIONS

Our observations coupled with previous data validate key predictions of the scanning model for Pol II initiation in yeast, which we term the shooting gallery. In this model, Pol II catalytic activity and the rate and processivity of Pol II scanning together with promoter sequence determine the distribution of TSSs and their usage.

Organismal benefits of transcription speed control at gene boundaries

RNA polymerase II (RNAPII) transcription is crucial for gene expression. RNAPII density peaks at gene boundaries, associating these key regions for gene expression control with limited RNAPII movement. The connections between RNAPII transcription speed and gene regulation in multicellular organisms are poorly understood. Here, we directly modulate RNAPII transcription speed by point mutations in the second largest subunit of RNAPII in Arabidopsis thaliana. A RNAPII mutation predicted to decelerate transcription is inviable, while accelerating RNAPII transcription confers phenotypes resembling auto-immunity. Nascent transcription profiling revealed that RNAPII complexes with accelerated transcription clear stalling sites at both gene ends, resulting in read-through transcription. The accelerated transcription mutant NRPB2-Y732F exhibits increased association with 5' splice site (5'SS) intermediates and enhanced splicing efficiency. Our findings highlight potential advantages of RNAPII stalling through local reduction in transcription speed to optimize gene expression for the development of multicellular organisms.

A Quantitative Genetic Interaction Map of HIV Infection

We have developed a platform for quantitative genetic interaction mapping using viral infectivity as a functional readout and constructed a viral host-dependency epistasis map (vE-MAP) of 356 human genes linked to HIV function, comprising >63,000 pairwise genetic perturbations. The vE-MAP provides an expansive view of the genetic dependencies underlying HIV infection and can be used to identify drug targets and study viral mutations. We found that the RNA deadenylase complex, CNOT, is a central player in the vE-MAP and show that knockout of CNOT1, 10, and 11 suppressed HIV infection in primary T cells by upregulating innate immunity pathways. This phenotype was rescued by deletion of IRF7, a transcription factor regulating interferon-stimulated genes, revealing a previously unrecognized host signaling pathway involved in HIV infection. The vE-MAP represents a generic platform that can be used to study the global effects of how different pathogens hijack and rewire the host during infection.

Spt5 modulates cotranscriptional spliceosome assembly in Saccharomyces cerevisiae

There is increasing evidence from yeast to humans that pre-mRNA splicing occurs mainly cotranscriptionally, such that splicing and transcription are functionally coupled. Currently, there is little insight into the contribution of the core transcription elongation machinery to cotranscriptional spliceosome assembly and pre-mRNA splicing. Spt5 is a member of the core transcription elongation machinery and an essential protein, whose absence in budding yeast causes defects in pre-mRNA splicing. To determine how Spt5 affects pre-mRNA splicing, we used the auxin-inducible degron system to conditionally deplete Spt5 in Saccharomyces cerevisiae and assayed effects on cotranscriptional spliceosome assembly and splicing. We show that Spt5 is needed for efficient splicing and for the accumulation of U5 snRNPs at intron-containing genes, and therefore for stable cotranscriptional assembly of spliceosomes. The defect in cotranscriptional spliceosome assembly can explain the relatively mild splicing defect as being a consequence of the failure of cotranscriptional splicing. Coimmunoprecipitation of Spt5 with core spliceosomal proteins and all spliceosomal snRNAs suggests a model whereby Spt5 promotes cotranscriptional pre-mRNA splicing by stabilizing the association of U5 snRNP with spliceosome complexes as they assemble on the nascent transcript. If this phenomenon is conserved in higher eukaryotes, it has the potential to be important for cotranscriptional regulation of alternative splicing.

Nascent RNA and the Coordination of Splicing with Transcription

At each active protein-encoding gene, nascent RNA is tethered to the DNA axis by elongating RNA polymerase II (Pol II) and is continuously altered by splicing and other processing events during its synthesis. This review discusses the development of three major methods that enable us to track the conversion of precursor messenger RNA (pre-mRNA) to messenger RNA (mRNA) products in vivo: live-cell imaging, metabolic labeling of RNA, and RNA-seq of purified nascent RNA. These approaches are complementary, addressing distinct issues of transcription rates and intron lifetimes alongside spatial information regarding the gene position of Pol II at which spliceosomes act. The findings will be placed in the context of active transcription units, each of which-because of the presence of nascent RNA, Pol II, and features of the chromatin environment-will recruit a potentially gene-specific constellation of RNA binding proteins and processing machineries.

CAGI 5 splicing challenge: Improved exon skipping and intron retention predictions with MMSplice

Pathogenic genetic variants often primarily affect splicing. However, it remains difficult to quantitatively predict whether and how genetic variants affect splicing. In 2018, the fifth edition of the Critical Assessment of Genome Interpretation proposed two splicing prediction challenges based on experimental perturbation assays: Vex-seq, assessing exon skipping, and MaPSy, assessing splicing efficiency. We developed a modular modeling framework, MMSplice, the performance of which was among the best on both challenges. Here we provide insights into the modeling assumptions of MMSplice and its individual modules. We furthermore illustrate how MMSplice can be applied in practice for individual genome interpretation, using the MMSplice VEP plugin and the Kipoi variant interpretation plugin, which are directly applicable to VCF files.

The telomeric Cdc13-Stn1-Ten1 complex regulates RNA polymerase II transcription

Specialized telomeric proteins have an essential role in maintaining genome stability through chromosome end protection and telomere length regulation. In the yeast Saccharomyces cerevisiae, the evolutionary conserved CST complex, composed of the Cdc13, Stn1 and Ten1 proteins, largely contributes to these functions. Here, we report genetic interactions between TEN1 and several genes coding for transcription regulators. Molecular assays confirmed this novel function of Ten1 and further established that it regulates the occupancies of RNA polymerase II and the Spt5 elongation factor within transcribed genes. Since Ten1, but also Cdc13 and Stn1, were found to physically associate with Spt5, we propose that Spt5 represents the target of CST in transcription regulation. Moreover, CST physically associates with Hmo1, previously shown to mediate the architecture of S-phase transcribed genes. The fact that, genome-wide, the promoters of genes down-regulated in the ten1-31 mutant are prefentially bound by Hmo1, leads us to propose a potential role for CST in synchronizing transcription with replication fork progression following head-on collisions.

Global donor and acceptor splicing site kinetics in human cells

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.

Spliceosome Profiling Visualizes Operations of a Dynamic RNP at Nucleotide Resolution

Tools to understand how the spliceosome functions in vivo have lagged behind advances in the structural biology of the spliceosome. Here, methods are described to globally profile spliceosome-bound pre-mRNA, intermediates, and spliced mRNA at nucleotide resolution. These tools are applied to three yeast species that span 600 million years of evolution. The sensitivity of the approach enables the detection of canonical and non-canonical events, including interrupted, recursive, and nested splicing. This application of statistical modeling uncovers independent roles for the size and position of the intron and the number of introns per transcript in substrate progression through the two catalytic stages. These include species-specific inputs suggestive of spliceosome-transcriptome coevolution. Further investigations reveal the ATP-dependent discard of numerous endogenous substrates after spliceosome assembly in vivo and connect this discard to intron retention, a form of splicing regulation. Spliceosome profiling is a quantitative, generalizable global technology used to investigate an RNP central to eukaryotic gene expression.

Perturbations of Transcription and Gene Expression-Associated Processes Alter Distribution of Cell Size Values in Saccharomyces cerevisiae

The question of what determines whether cells are big or small has been the focus of many studies because it is thought that such determinants underpin the coupling of cell growth with cell division. In contrast, what determines the overall pattern of how cell size is distributed within a population of wild type or mutant cells has received little attention. Knowing how cell size varies around a characteristic pattern could shed light on the processes that generate such a pattern and provide a criterion to identify its genetic basis. Here, we show that cell size values of wild type Saccharomyces cerevisiae cells fit a gamma distribution, in haploid and diploid cells, and under different growth conditions. To identify genes that influence this pattern, we analyzed the cell size distributions of all single-gene deletion strains in Saccharomyces cerevisiae. We found that yeast strains which deviate the most from the gamma distribution are enriched for those lacking gene products functioning in gene expression, especially those in transcription or transcription-linked processes. We also show that cell size is increased in mutants carrying altered activity substitutions in Rpo21p/Rpb1, the largest subunit of RNA polymerase II (Pol II). Lastly, the size distribution of cells carrying extreme altered activity Pol II substitutions deviated from the expected gamma distribution. Our results are consistent with the idea that genetic defects in widely acting transcription factors or Pol II itself compromise both cell size homeostasis and how the size of individual cells is distributed in a population.

Pre-mRNA Splicing in the Nuclear Landscape

Eukaryotic gene expression requires the cumulative activity of multiple molecular machines to synthesize and process newly transcribed pre-messenger RNA. Introns, the noncoding regions in pre-mRNA, must be removed by the spliceosome, which assembles on the pre-mRNA as it is transcribed by RNA polymerase II (Pol II). The assembly and activity of the spliceosome can be modulated by features including the speed of transcription elongation, chromatin, post-translational modifications of Pol II and histone tails, and other RNA processing events like 5'-end capping. Here, we review recent work that has revealed cooperation and coordination among co-transcriptional processing events and speculate on new avenues of research. We anticipate new mechanistic insights capable of unraveling the relative contribution of coupled processing to gene expression.

Genetic analysis reveals functions of atypical polyubiquitin chains

Although polyubiquitin chains linked through all lysines of ubiquitin exist, specific functions are well-established only for lysine-48 and lysine-63 linkages in Saccharomyces cerevisiae. To uncover pathways regulated by distinct linkages, genetic interactions between a gene deletion library and a panel of lysine-to-arginine ubiquitin mutants were systematically identified. The K11R mutant had strong genetic interactions with threonine biosynthetic genes. Consistently, we found that K11R mutants import threonine poorly. The K11R mutant also exhibited a strong genetic interaction with a subunit of the anaphase-promoting complex (APC), suggesting a role in cell cycle regulation. K11-linkages are important for vertebrate APC function, but this was not previously described in yeast. We show that the yeast APC also modifies substrates with K11-linkages in vitro, and that those chains contribute to normal APC-substrate turnover in vivo. This study reveals comprehensive genetic interactomes of polyubiquitin chains and characterizes the role of K11-chains in two biological pathways.

Use of genetic and chemical synthetic lethality as probes of complexity in bacterial cell systems

Different conditions and genomic contexts are known to have an impact on gene essentiality and interactions. Synthetic lethal interactions occur when a combination of perturbations, either genetic or chemical, result in a more profound fitness defect than expected based on the effect of each perturbation alone. Synthetic lethality in bacterial systems has long been studied; however, during the past decade, the emerging fields of genomics and chemical genomics have led to an increase in the scale and throughput of these studies. Here, we review the concepts of genomics and chemical genomics in the context of synthetic lethality and their revolutionary roles in uncovering novel biology such as the characterization of genes of unknown function and in antibacterial drug discovery. We provide an overview of the methodologies, examples and challenges of both genetic and chemical synthetic lethal screening platforms. Finally, we discuss how to apply genetic and chemical synthetic lethal approaches to rationalize the synergies of drugs, screen for new and improved antibacterial therapies and predict drug mechanism of action.

Adaptive transcription-splicing resynchronization upon losing an essential splicing factor

Essential genes form the core of a genome and are therefore thought to be indispensable for cellular viability. However, recent findings have challenged this notion in that cells may survive in the absence of some essential genes provided that relevant genetic modifiers are in existence. We therefore hypothesized that the loss of an essential gene may not always be fatefully detrimental; instead, it may pave the way towards genome evolution. We experimentally tested this hypothesis in the context of pre-messenger RNA splicing by evolving yeast cells harbouring a permanent loss of the essential splicing factor Prp28 in the presence of a genetic modifier. Here, we show that cellular fitness can be restored by compensatory mutations that alter either the splicing machinery per se or the Spt-Ada-Gcn5 acetyltransferase transcription co-activator complex in the cells with no Prp28. Biochemical and genetic analysis revealed that slowing down transcription compensates for splicing deficiency, which in turn boosts cellular fitness. In addition, we found that inefficient splicing also conversely decreases nascent RNA production. Taken together, our data suggest that transcription-splicing synchronization contributes to robustness in the gene-expression pathway and argue that the intrinsic interconnectivity within a biological system can be exploited for compensatory evolution and system re-optimization.