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.

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.

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.

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.

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.

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.

Magnetic Resonance Acoustic Radiation Force Imaging (MR-ARFI) for the monitoring of High Intensity Focused Ultrasound (HIFU) ablation in anisotropic tissue.

OBJECTIVE: We introduce a non-invasive MR-Acoustic Radiation Force Imaging (ARFI)-based elastography method that provides both the local shear modulus and temperature maps for the monitoring of High Intensity Focused Ultrasound (HIFU) therapy. MATERIALS AND METHODS: To take tissue anisotropy into account, the local shear modulus µ is determined in selected radial directions around the focal spot by fitting the phase profiles to a linear viscoelastic model, including tissue-specific mechanical relaxation time τ. MR-ARFI was evaluated on a calibrated phantom, then applied to the monitoring of HIFU in a gel phantom, ex vivo and in vivo porcine muscle tissue, in parallel with MR-thermometry. RESULTS: As expected, the shear modulus polar maps reflected the isotropy of phantoms and the anisotropy of muscle. In the HIFU monitoring experiments, both the shear modulus polar map and the thermometry map were updated with every pair of MR-ARFI phase images acquired with opposite MR-ARFI-encoding. The shear modulus was found to decrease (phantom and ex vivo) or increase (in vivo) during heating, before remaining steady during the cooling phase. The mechanical relaxation time, estimated pre- and post-HIFU, was found to vary in muscle tissue. DISCUSSION: MR-ARFI allowed for monitoring of viscoelasticity changes around the HIFU focal spot even in anisotropic muscle tissue.

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.

Monitoring MR-guided high intensity focused ultrasound therapy using transient supersonic shear wave MR-elastography.

Objective.The aim of the paper is to propose an all-in-one method based on magnetic resonance-supersonic shear wave imaging (MR-SSI) and proton resonance frequency shift (PRFS) to monitor high intensity focused ultrasound (HIFU) thermal ablations.Approach.Mechanical properties have been shown to be related to tissue damage induced by thermal ablations. Monitoring elasticity in addition to temperature changes may help in ensuring the efficacy and the accuracy of HIFU therapies. For this purpose, an MR-SSI method has been developed where the ultrasonic transducer is used for both mechanical wave generation and thermal ablation. Transient quasi-planar shear waves are generated using the acoustic radiation force, and their propagation is monitored in motion-sensitized phase MR images. Using a single-shot gradient-echo echo-planar-imaging sequence, MR images can be acquired at a sufficiently high temporal resolution to provide an update of PRFS thermometry and MR-SSI elastography maps in real time.Main results.The proposed method was first validated on a calibrated elasticity phantom, in which both the possibility to detect inclusions with different stiffness and repeatability were demonstrated. The standard deviation between the 8 performed measurements was 2% on the background of the phantom and 11%, at most, on the inclusions. A second experiment consisted in performing a HIFU heating in a gelatin phantom. The temperature increase was estimated to be 9 °C and the shear modulus was found to decrease from 2.9 to 1.8 kPa, reflecting the gel softening around the HIFU focus, whereas it remained steady in non-heated areas.Significance.The proposed MR-SSI technique allows monitoring HIFU ablations using thermometry and elastography simultaneously, without the need for an additional external mechanical exciter such as those used in MR elastography.

Deep Intronic FGF14 GAA Repeat Expansion in Late-Onset Cerebellar Ataxia.

BACKGROUND: The late-onset cerebellar ataxias (LOCAs) have largely resisted molecular diagnosis. METHODS: We sequenced the genomes of six persons with autosomal dominant LOCA who were members of three French Canadian families and identified a candidate pathogenic repeat expansion. We then tested for association between the repeat expansion and disease in two independent case-control series - one French Canadian (66 patients and 209 controls) and the other German (228 patients and 199 controls). We also genotyped the repeat in 20 Australian and 31 Indian index patients. We assayed gene and protein expression in two postmortem cerebellum specimens and two induced pluripotent stem-cell (iPSC)-derived motor-neuron cell lines. RESULTS: In the six French Canadian patients, we identified a GAA repeat expansion deep in the first intron of FGF14, which encodes fibroblast growth factor 14. Cosegregation of the repeat expansion with disease in the families supported a pathogenic threshold of at least 250 GAA repeats ([GAA]≥250). There was significant association between FGF14 (GAA)≥250 expansions and LOCA in the French Canadian series (odds ratio, 105.60; 95% confidence interval [CI], 31.09 to 334.20; P<0.001) and in the German series (odds ratio, 8.76; 95% CI, 3.45 to 20.84; P<0.001). The repeat expansion was present in 61%, 18%, 15%, and 10% of French Canadian, German, Australian, and Indian index patients, respectively. In total, we identified 128 patients with LOCA who carried an FGF14 (GAA)≥250 expansion. Postmortem cerebellum specimens and iPSC-derived motor neurons from patients showed reduced expression of FGF14 RNA and protein. CONCLUSIONS: A dominantly inherited deep intronic GAA repeat expansion in FGF14 was found to be associated with LOCA. (Funded by Fondation Groupe Monaco and others.).

RNA Polymerase III Subunit Mutations in Genetic Diseases.

RNA polymerase (Pol) III transcribes small untranslated RNAs such as 5S ribosomal RNA, transfer RNAs, and U6 small nuclear RNA. Because of the functions of these RNAs, Pol III transcription is best known for its essential contribution to RNA maturation and translation. Surprisingly, it was discovered in the last decade that various inherited mutations in genes encoding nine distinct subunits of Pol III cause tissue-specific diseases rather than a general failure of all vital functions. Mutations in the POLR3A, POLR3C, POLR3E and POLR3F subunits are associated with susceptibility to varicella zoster virus-induced encephalitis and pneumonitis. In addition, an ever-increasing number of distinct mutations in the POLR3A, POLR3B, POLR1C and POLR3K subunits cause a spectrum of neurodegenerative diseases, which includes most notably hypomyelinating leukodystrophy. Furthermore, other rare diseases are also associated with mutations in genes encoding subunits of Pol III (POLR3H, POLR3GL) and the BRF1 component of the TFIIIB transcription initiation factor. Although the causal relationship between these mutations and disease development is widely accepted, the exact molecular mechanisms underlying disease pathogenesis remain enigmatic. Here, we review the current knowledge on the functional impact of specific mutations, possible Pol III-related disease-causing mechanisms, and animal models that may help to better understand the links between Pol III mutations and disease.

Revealing nascent RNA processing dynamics with nano-COP.

During maturation, eukaryotic precursor RNAs undergo processing events including intron splicing, 3'-end cleavage, and polyadenylation. Here we describe nanopore analysis of co-transcriptional processing (nano-COP), a method for probing the timing and patterns of RNA processing. An extension of native elongating transcript sequencing, which quantifies transcription genome-wide through short-read sequencing of nascent RNA 3' ends, nano-COP uses long-read nascent RNA sequencing to observe global patterns of RNA processing. First, nascent RNA is stringently purified through a combination of 4-thiouridine metabolic labeling and cellular fractionation. In contrast to cDNA or short-read-based approaches relying on reverse transcription or amplification, the sample is sequenced directly through nanopores to reveal the native context of nascent RNA. nano-COP identifies both active transcription sites and splice isoforms of single RNA molecules during synthesis, providing insight into patterns of intron removal and the physical coupling between transcription and splicing. The nano-COP protocol yields data within 3 d.

Channelopathies Are a Frequent Cause of Genetic Ataxias Associated with Cerebellar Atrophy.

BACKGROUND: Cerebellar atrophy is a nonspecific imaging finding observed in a number of neurological disorders. Genetic ataxias associated with cerebellar atrophy are a heterogeneous group of conditions, rendering the approach to diagnosis challenging. OBJECTIVES: To define the spectrum of genetic ataxias associated with cerebellar atrophy in a Canadian cohort and the diagnostic yield of exome sequencing for this group of conditions. METHODS: A total of 92 participants from 66 families with cerebellar atrophy were recruited for this multicenter prospective cohort study. Exome sequencing was performed for all participants between 2011 and 2017 as part of 1 of 2 national research programs, Finding of Rare Genetic Disease Genes or Enhanced Care for Rare Genetic Diseases in Canada. RESULTS: A genetic diagnosis was established in 53% of families (35/66). Pathogenic variants were found in 21 known genes, providing a diagnosis for 31/35 families (89%), and in 4 novel genes, accounting for 4/35 families (11%). Of the families, 31/66 (47%) remained without a genetic diagnosis. The most common diagnoses were channelopathies, which were established in 9/35 families (26%). Additional clinical findings provided useful clues to specific diagnoses. CONCLUSIONS: We report on the high frequency of channelopathies as a cause of genetic ataxias associated with cerebellar atrophy and the utility of exome sequencing for this group of conditions.

Splicing Kinetics and Coordination Revealed by Direct Nascent RNA Sequencing through Nanopores.

Understanding how splicing events are coordinated across numerous introns in metazoan RNA transcripts requires quantitative analyses of transient RNA processing events in living cells. We developed nanopore analysis of co-transcriptional processing (nano-COP), in which nascent RNAs are directly sequenced through nanopores, exposing the dynamics and patterns of RNA splicing without biases introduced by amplification. Long nano-COP reads reveal that, in human and Drosophila cells, splicing occurs after RNA polymerase II transcribes several kilobases of pre-mRNA, suggesting that metazoan splicing transpires distally from the transcription machinery. Inhibition of the branch-site recognition complex SF3B rapidly diminished global co-transcriptional splicing. We found that splicing order does not strictly follow the order of transcription and is associated with cis-acting elements, alternative splicing, and RNA-binding factors. Further, neighboring introns in human cells tend to be spliced concurrently, implying that splicing of these introns occurs cooperatively. Thus, nano-COP unveils the organizational complexity of RNA processing.

The leukodystrophy mutation Polr3b R103H causes homozygote mouse embryonic lethality and impairs RNA polymerase III biogenesis.

Recessive mutations in the ubiquitously expressed POLR3A and POLR3B genes are the most common cause of POLR3-related hypomyelinating leukodystrophy (POLR3-HLD), a rare childhood-onset disorder characterized by deficient cerebral myelin formation and cerebellar atrophy. POLR3A and POLR3B encode the two catalytic subunits of RNA Polymerase III (Pol III), which synthesizes numerous small non-coding RNAs. We recently reported that mice homozygous for the Polr3a mutation c.2015G > A (p.Gly672Glu) have no neurological abnormalities and thus do not recapitulate the human POLR3-HLD phenotype. To determine if other POLR3-HLD mutations can cause a leukodystrophy phenotype in mouse, we characterized mice carrying the Polr3b mutation c.308G > A (p.Arg103His). Surprisingly, homozygosity for this mutation was embryonically lethal with only wild-type and heterozygous animals detected at embryonic day 9.5. Using proteomics in a human cell line, we found that the POLR3B R103H mutation severely impairs assembly of the Pol III complex. We next generated Polr3aG672E/G672E/Polr3b+/R103Hdouble mutant mice but observed that this additional mutation was insufficient to elicit a neurological or transcriptional phenotype. Taken together with our previous study on Polr3a G672E mice, our results indicate that missense mutations in Polr3a and Polr3b can variably impair mouse development and Pol III function. Developing a proper model of POLR3-HLD is crucial to gain insights into the pathophysiological mechanisms involved in this devastating neurodegenerative disease.