Genetic regulation of nascent RNA maturation revealed by direct RNA nanopore sequencing.

Quantitative trait loci analyses have revealed an important role for genetic variants in regulating alternative splicing (AS) and alternative cleavage and polyadenylation (APA) in humans. Yet, these studies are generally performed with mature mRNA, so they report on the outcome rather than the processes of RNA maturation and thus may overlook how variants directly modulate pre-mRNA processing. The order in which the many introns of a human gene are removed can substantially influence AS, while nascent RNA polyadenylation can affect RNA stability and decay. However, how splicing order and poly(A) tail length are regulated by genetic variation has never been explored. Here, we used direct RNA nanopore sequencing to investigate allele-specific pre-mRNA maturation in 12 human lymphoblastoid cell lines. We found frequent splicing order differences between alleles and uncovered significant single nucleotide polymorphism (SNP)-splicing order associations in 17 genes. This included SNPs located in or near splice sites as well as more distal intronic and exonic SNPs. Moreover, several genes showed allele-specific poly(A) tail lengths, many of which also had a skewed allelic abundance ratio. HLA class I transcripts, which encode proteins that play an essential role in antigen presentation, showed the most allele-specific splicing orders, which frequently co-occurred with allele-specific AS, APA or poly(A) tail length differences. Together, our results expose new layers of genetic regulation of pre-mRNA maturation and highlight the power of long-read RNA sequencing for allele-specific analyses.

Context-specific inhibition of mitochondrial ribosomes by phenicol and oxazolidinone antibiotics.

The antibiotics chloramphenicol (CHL) and oxazolidinones including linezolid (LZD) are known to inhibit mitochondrial translation. This can result in serious, potentially deadly, side effects when used therapeutically. Although the mechanism by which CHL and LZD inhibit bacterial ribosomes has been elucidated in detail, their mechanism of action against mitochondrial ribosomes has yet to be explored. CHL and oxazolidinones bind to the ribosomal peptidyl transfer center (PTC) of the bacterial ribosome and prevent incorporation of incoming amino acids under specific sequence contexts, causing ribosomes to stall only at certain sequences. Through mitoribosome profiling, we show that inhibition of mitochondrial ribosomes is similarly context-specific - CHL and LZD lead to mitoribosome stalling primarily when there is an alanine, serine, or threonine in the penultimate position of the nascent peptide chain. We further validate context-specific stalling through in vitro translation assays. A high resolution cryo-EM structure of LZD bound to the PTC of the human mitoribosome shows extensive similarity to the mode of bacterial inhibition and also suggests potential avenues for altering selectivity. Our findings could help inform the rational development of future, less mitotoxic, antibiotics, which are critically needed in the current era of increasing antimicrobial resistance.

Structural basis of LRPPRC-SLIRP-dependent translation by the mitoribosome.

In mammalian mitochondria, mRNAs are cotranscriptionally stabilized by the protein factor LRPPRC (leucine-rich pentatricopeptide repeat-containing protein). Here, we characterize LRPPRC as an mRNA delivery factor and report its cryo-electron microscopy structure in complex with SLIRP (SRA stem-loop-interacting RNA-binding protein), mRNA and the mitoribosome. The structure shows that LRPPRC associates with the mitoribosomal proteins mS39 and the N terminus of mS31 through recognition of the LRPPRC helical repeats. Together, the proteins form a corridor for handoff of the mRNA. The mRNA is directly bound to SLIRP, which also has a stabilizing function for LRPPRC. To delineate the effect of LRPPRC on individual mitochondrial transcripts, we used RNA sequencing, metabolic labeling and mitoribosome profiling, which showed a transcript-specific influence on mRNA translation efficiency, with cytochrome c oxidase subunit 1 and 2 translation being the most affected. Our data suggest that LRPPRC-SLIRP acts in recruitment of mitochondrial mRNAs to modulate their translation. Collectively, the data define LRPPRC-SLIRP as a regulator of the mitochondrial gene expression system.

RNA polymerases reshape chromatin architecture and couple transcription on individual fibers.

RNA polymerases must initiate and pause within a complex chromatin environment, surrounded by nucleosomes and other transcriptional machinery. This environment creates a spatial arrangement along individual chromatin fibers ripe for both competition and coordination, yet these relationships remain largely unknown owing to the inherent limitations of traditional structural and sequencing methodologies. To address this, we employed long-read chromatin fiber sequencing (Fiber-seq) in Drosophila to visualize RNA polymerase (Pol) within its native chromatin context with single-molecule precision along up to 30 kb fibers. We demonstrate that Fiber-seq enables the identification of individual Pol II, nucleosome, and transcription factor footprints, revealing Pol II pausing-driven destabilization of downstream nucleosomes. Furthermore, we demonstrate pervasive direct distance-dependent transcriptional coupling between nearby Pol II genes, Pol III genes, and transcribed enhancers, modulated by local chromatin architecture. Overall, transcription initiation reshapes surrounding nucleosome architecture and couples nearby transcriptional machinery along individual chromatin fibers.

Genome-wide quantification of RNA flow across subcellular compartments reveals determinants of the mammalian transcript life cycle.

Dissecting the regulatory mechanisms controlling mammalian transcripts from production to degradation requires quantitative measurements of mRNA flow across the cell. We developed subcellular TimeLapse-seq to measure the rates at which RNAs are released from chromatin, exported from the nucleus, loaded onto polysomes, and degraded within the nucleus and cytoplasm in human and mouse cells. These rates varied substantially, yet transcripts from genes with related functions or targeted by the same transcription factors and RNA-binding proteins flowed across subcellular compartments with similar kinetics. Verifying these associations uncovered a link between DDX3X and nuclear export. For hundreds of RNA metabolism genes, most transcripts with retained introns were degraded by the nuclear exosome, while the remaining molecules were exported with stable cytoplasmic lifespans. Transcripts residing on chromatin for longer had extended poly(A) tails, whereas the reverse was observed for cytoplasmic mRNAs. Finally, machine learning identified molecular features that predicted the diverse life cycles of mRNAs.

Central dogma rates in human mitochondria.

In human cells, the nuclear and mitochondrial genomes engage in a complex interplay to produce dual-encoded oxidative phosphorylation (OXPHOS) complexes. The coordination of these dynamic gene expression processes is essential for producing matched amounts of OXPHOS protein subunits. This review focuses on our current understanding of the mitochondrial central dogma rates, highlighting the striking differences in gene expression rates between mitochondrial and nuclear genes. We synthesize a coherent model of mitochondrial gene expression kinetics, highlighting the emerging principles and emphasizing where more precise measurements would be beneficial. Such an understanding is pivotal for grasping the unique aspects of mitochondrial function and its role in cellular energetics, and it has profound implications for aging, metabolic disorders, and neurodegenerative diseases.

Timing is everything: advances in quantifying splicing kinetics.

Splicing is a highly regulated process critical for proper pre-mRNA maturation and the maintenance of a healthy cellular environment. Splicing events are impacted by ongoing transcription, neighboring splicing events, and cis and trans regulatory factors on the respective pre-mRNA transcript. Within this complex regulatory environment, splicing kinetics have the potential to influence splicing outcomes but have historically been challenging to study in vivo. In this review, we highlight recent technological advancements that have enabled measurements of global splicing kinetics and of the variability of splicing kinetics at single introns. We demonstrate how identifying features that are correlated with splicing kinetics has increased our ability to form potential models for how splicing kinetics may be regulated in vivo.

Post-transcriptional splicing can occur in a slow-moving zone around the gene.

Splicing is the stepwise molecular process by which introns are removed from pre-mRNA and exons are joined together to form mature mRNA sequences. The ordering and spatial distribution of these steps remain controversial, with opposing models suggesting splicing occurs either during or after transcription. We used single-molecule RNA FISH, expansion microscopy, and live-cell imaging to reveal the spatiotemporal distribution of nascent transcripts in mammalian cells. At super-resolution levels, we found that pre-mRNA formed clouds around the transcription site. These clouds indicate the existence of a transcription-site-proximal zone through which RNA move more slowly than in the nucleoplasm. Full-length pre-mRNA undergo continuous splicing as they move through this zone following transcription, suggesting a model in which splicing can occur post-transcriptionally but still within the proximity of the transcription site, thus seeming co-transcriptional by most assays. These results may unify conflicting reports of co-transcriptional versus post-transcriptional splicing.

A kinetic dichotomy between mitochondrial and nuclear gene expression processes.

Oxidative phosphorylation (OXPHOS) complexes, encoded by both mitochondrial and nuclear DNA, are essential producers of cellular ATP, but how nuclear and mitochondrial gene expression steps are coordinated to achieve balanced OXPHOS subunit biogenesis remains unresolved. Here, we present a parallel quantitative analysis of the human nuclear and mitochondrial messenger RNA (mt-mRNA) life cycles, including transcript production, processing, ribosome association, and degradation. The kinetic rates of nearly every stage of gene expression differed starkly across compartments. Compared with nuclear mRNAs, mt-mRNAs were produced 1,100-fold more, degraded 7-fold faster, and accumulated to 160-fold higher levels. Quantitative modeling and depletion of mitochondrial factors LRPPRC and FASTKD5 identified critical points of mitochondrial regulatory control, revealing that the mitonuclear expression disparities intrinsically arise from the highly polycistronic nature of human mitochondrial pre-mRNA. We propose that resolving these differences requires a 100-fold slower mitochondrial translation rate, illuminating the mitoribosome as a nexus of mitonuclear co-regulation.

Proteome-scale movements and compartment connectivity during the eukaryotic cell cycle.

Cell cycle progression relies on coordinated changes in the composition and subcellular localization of the proteome. By applying two distinct convolutional neural networks on images of millions of live yeast cells, we resolved proteome-level dynamics in both concentration and localization during the cell cycle, with resolution of ∼20 subcellular localization classes. We show that a quarter of the proteome displays cell cycle periodicity, with proteins tending to be controlled either at the level of localization or concentration, but not both. Distinct levels of protein regulation are preferentially utilized for different aspects of the cell cycle, with changes in protein concentration being mostly involved in cell cycle control and changes in protein localization in the biophysical implementation of the cell cycle program. We present a resource for exploring global proteome dynamics during the cell cycle, which will aid in understanding a fundamental biological process at a systems level.

Single-nucleoid architecture reveals heterogeneous packaging of mitochondrial DNA.

Cellular metabolism relies on the regulation and maintenance of mitochondrial DNA (mtDNA). Hundreds to thousands of copies of mtDNA exist in each cell, yet because mitochondria lack histones or other machinery important for nuclear genome compaction, it remains unresolved how mtDNA is packaged into individual nucleoids. In this study, we used long-read single-molecule accessibility mapping to measure the compaction of individual full-length mtDNA molecules at near single-nucleotide resolution. We found that, unlike the nuclear genome, human mtDNA largely undergoes all-or-none global compaction, with most nucleoids existing in an inaccessible, inactive state. Highly accessible mitochondrial nucleoids are co-occupied by transcription and replication components and selectively form a triple-stranded displacement loop structure. In addition, we showed that the primary nucleoid-associated protein TFAM directly modulates the fraction of inaccessible nucleoids both in vivo and in vitro, acting consistently with a nucleation-and-spreading mechanism to coat and compact mitochondrial nucleoids. Together, these findings reveal the primary architecture of mtDNA packaging and regulation in human cells.

Regulators of mitonuclear balance link mitochondrial metabolism to mtDNA expression.

Mitochondrial oxidative phosphorylation (OXPHOS) complexes are assembled from proteins encoded by both nuclear and mitochondrial DNA. These dual-origin enzymes pose a complex gene regulatory challenge for cells requiring coordinated gene expression across organelles. To identify genes involved in dual-origin protein complex synthesis, we performed fluorescence-activated cell-sorting-based genome-wide screens analysing mutant cells with unbalanced levels of mitochondrial- and nuclear-encoded subunits of Complex IV. We identified genes involved in OXPHOS biogenesis, including two uncharacterized genes: PREPL and NME6. We found that PREPL specifically impacts Complex IV biogenesis by acting at the intersection of mitochondrial lipid metabolism and protein synthesis, whereas NME6, an uncharacterized nucleoside diphosphate kinase, controls OXPHOS biogenesis through multiple mechanisms reliant on its NDPK domain. Firstly, NME6 forms a complex with RCC1L, which together perform nucleoside diphosphate kinase activity to maintain local mitochondrial pyrimidine triphosphate levels essential for mitochondrial RNA abundance. Secondly, NME6 modulates the activity of mitoribosome regulatory complexes, altering mitoribosome assembly and mitochondrial RNA pseudouridylation. Taken together, we propose that NME6 acts as a link between compartmentalized mitochondrial metabolites and mitochondrial gene expression.

Elongation rate of RNA polymerase II affects pausing patterns across 3' UTRs.

Yeast mRNAs are polyadenylated at multiple sites in their 3' untranslated regions (3' UTRs), and poly(A) site usage is regulated by the rate of transcriptional elongation by RNA polymerase II (Pol II). Slow Pol II derivatives favor upstream poly(A) sites, and fast Pol II derivatives favor downstream poly(A) sites. Transcriptional elongation and polyadenylation are linked at the nucleotide level, presumably reflecting Pol II dwell time at each residue that influences the level of polyadenylation. Here, we investigate the effect of Pol II elongation rate on pausing patterns and the relationship between Pol II pause sites and poly(A) sites within 3' UTRs. Mutations that affect Pol II elongation rate alter sequence preferences at pause sites within 3' UTRs, and pausing preferences differ between 3' UTRs and coding regions. In addition, sequences immediately flanking the pause sites show preferences that are largely independent of Pol II speed. In wild-type cells, poly(A) sites are preferentially located < 50 nucleotides upstream from Pol II pause sites, but this spatial relationship is diminished in cells harboring Pol II speed mutants. Based on a random forest classifier, Pol II pause sites are modestly predicted by the distance to poly(A) sites but are better predicted by the chromatin landscape in Pol II speed derivatives. Transcriptional regulatory proteins can influence the relationship between Pol II pausing and polyadenylation but in a manner distinct from Pol II elongation rate derivatives. These results indicate a complex relationship between Pol II pausing and polyadenylation.

RNA polymerase II pausing temporally coordinates cell cycle progression and erythroid differentiation.

Controlled release of promoter-proximal paused RNA polymerase II (RNA Pol II) is crucial for gene regulation. However, studying RNA Pol II pausing is challenging, as pause-release factors are almost all essential. In this study, we identified heterozygous loss-of-function mutations in SUPT5H, which encodes SPT5, in individuals with ß-thalassemia. During erythropoiesis in healthy human cells, cell cycle genes were highly paused as cells transition from progenitors to precursors. When the pathogenic mutations were recapitulated by SUPT5H editing, RNA Pol II pause release was globally disrupted, and as cells began transitioning from progenitors to precursors, differentiation was delayed, accompanied by a transient lag in erythroid-specific gene expression and cell cycle kinetics. Despite this delay, cells terminally differentiate, and cell cycle phase distributions normalize. Therefore, hindering pause release perturbs proliferation and differentiation dynamics at a key transition during erythropoiesis, identifying a role for RNA Pol II pausing in temporally coordinating the cell cycle and erythroid differentiation.

Intronic small nucleolar RNAs regulate host gene splicing through base pairing with their adjacent intronic sequences.

BACKGROUND: Small nucleolar RNAs (snoRNAs) are abundant noncoding RNAs best known for their involvement in ribosomal RNA maturation. In mammals, most expressed snoRNAs are embedded in introns of longer genes and produced through transcription and splicing of their host. Intronic snoRNAs were long viewed as inert passengers with little effect on host expression. However, a recent study reported a snoRNA influencing the splicing and ultimate output of its host gene. Overall, the general contribution of intronic snoRNAs to host expression remains unclear. RESULTS: Computational analysis of large-scale human RNA-RNA interaction datasets indicates that 30% of detected snoRNAs interact with their host transcripts. Many snoRNA-host duplexes are located near alternatively spliced exons and display high sequence conservation suggesting a possible role in splicing regulation. The study of the model SNORD2-EIF4A2 duplex indicates that the snoRNA interaction with the host intronic sequence conceals the branch point leading to decreased inclusion of the adjacent alternative exon. Extended SNORD2 sequence containing the interacting intronic region accumulates in sequencing datasets in a cell-type-specific manner. Antisense oligonucleotides and mutations that disrupt the formation of the snoRNA-intron structure promote the splicing of the alternative exon, shifting the EIF4A2 transcript ratio away from nonsense-mediated decay. CONCLUSIONS: Many snoRNAs form RNA duplexes near alternative exons of their host transcripts, placing them in optimal positions to control host output as shown for the SNORD2-EIF4A2 model system. Overall, our study supports a more widespread role for intronic snoRNAs in the regulation of their host transcript maturation.

Pre-mRNA splicing order is predetermined and maintains splicing fidelity across multi-intronic transcripts.

Combinatorially, intron excision within a given nascent transcript could proceed down any of thousands of paths, each of which would expose different dynamic landscapes of cis-elements and contribute to alternative splicing. In this study, we found that post-transcriptional multi-intron splicing order in human cells is largely predetermined, with most genes spliced in one or a few predominant orders. Strikingly, these orders were conserved across cell types and stages of motor neuron differentiation. Introns flanking alternatively spliced exons were frequently excised last, after their neighboring introns. Perturbations to the spliceosomal U2 snRNA altered the preferred splicing order of many genes, and these alterations were associated with the retention of other introns in the same transcript. In one gene, early removal of specific introns was sufficient to induce delayed excision of three proximal introns, and this delay was caused by two distinct cis-regulatory mechanisms. Together, our results demonstrate that multi-intron splicing order in human cells is predetermined, is influenced by a component of the spliceosome and ensures splicing fidelity across long pre-mRNAs.

RNA Polymerase II pausing temporally coordinates cell cycle progression and erythroid differentiation.

The controlled release of promoter-proximal paused RNA polymerase II (Pol II) into productive elongation is a major step in gene regulation. However, functional analysis of Pol II pausing is difficult because factors that regulate pause release are almost all essential. In this study, we identified heterozygous loss-of-function mutations in SUPT5H , which encodes SPT5, in individuals with ß-thalassemia unlinked to HBB mutations. During erythropoiesis in healthy human cells, cell cycle genes were highly paused at the transition from progenitors to precursors. When the pathogenic mutations were recapitulated by SUPT5H editing, Pol II pause release was globally disrupted, and the transition from progenitors to precursors was delayed, marked by a transient lag in erythroid-specific gene expression and cell cycle kinetics. Despite this delay, cells terminally differentiate, and cell cycle phase distributions normalize. Therefore, hindering pause release perturbs proliferation and differentiation dynamics at a key transition during erythropoiesis, revealing a role for Pol II pausing in the temporal coordination between the cell cycle and differentiation.

A kinetic dichotomy between mitochondrial and nuclear gene expression drives OXPHOS biogenesis.

Oxidative phosphorylation (OXPHOS) complexes, encoded by both mitochondrial and nuclear DNA, are essential producers of cellular ATP, but how nuclear and mitochondrial gene expression steps are coordinated to achieve balanced OXPHOS biogenesis remains unresolved. Here, we present a parallel quantitative analysis of the human nuclear and mitochondrial messenger RNA (mt-mRNA) life cycles, including transcript production, processing, ribosome association, and degradation. The kinetic rates of nearly every stage of gene expression differed starkly across compartments. Compared to nuclear mRNAs, mt-mRNAs were produced 700-fold higher, degraded 5-fold faster, and accumulated to 170-fold higher levels. Quantitative modeling and depletion of mitochondrial factors, LRPPRC and FASTKD5, identified critical points of mitochondrial regulatory control, revealing that the mitonuclear expression disparities intrinsically arise from the highly polycistronic nature of human mitochondrial pre-mRNA. We propose that resolving these differences requires a 100-fold slower mitochondrial translation rate, illuminating the mitoribosome as a nexus of mitonuclear co-regulation.

Genome-wide screens for mitonuclear co-regulators uncover links between compartmentalized metabolism and mitochondrial gene expression.

Mitochondrial oxidative phosphorylation (OXPHOS) complexes are assembled from proteins encoded by both nuclear and mitochondrial DNA. These dual-origin enzymes pose a complex gene regulatory challenge for cells, in which gene expression must be coordinated across organelles using distinct pools of ribosomes. How cells produce and maintain the accurate subunit stoichiometries for these OXPHOS complexes remains largely unknown. To identify genes involved in dual-origin protein complex synthesis, we performed FACS-based genome-wide screens analyzing mutant cells with unbalanced levels of mitochondrial- and nuclear-encoded subunits of cytochrome c oxidase (Complex IV). We identified novel genes involved in OXPHOS biogenesis, including two uncharacterized genes: PREPL and NME6 . We found that PREPL specifically regulates Complex IV biogenesis by interacting with mitochondrial protein synthesis machinery, while NME6, an uncharacterized nucleoside diphosphate kinase (NDPK), controls OXPHOS complex biogenesis through multiple mechanisms reliant on its NDPK domain. First, NME6 maintains local mitochondrial pyrimidine triphosphate levels essential for mitochondrial RNA abundance. Second, through stabilizing interactions with RCC1L, NME6 modulates the activity of mitoribosome regulatory complexes, leading to disruptions in mitoribosome assembly and mitochondrial RNA pseudouridylation. Taken together, we propose that NME6 acts as a link between compartmentalized mitochondrial metabolites and mitochondrial gene expression. Finally, we present these screens as a resource, providing a catalog of genes involved in mitonuclear gene regulation and OXPHOS biogenesis.