Noncoding mutations cause super-enhancer retargeting resulting in protein synthesis dysregulation during B cell lymphoma progression

Whole-genome sequencing of longitudinal tumor pairs representing transformation of follicular lymphoma to high-grade B cell lymphoma with MYC and BCL2 rearrangements (double-hit lymphoma) identified coding and noncoding genomic alterations acquired during lymphoma progression. Many of these transformation-associated alterations recurrently and focally occur at topologically associating domain resident regulatory DNA elements, including H3K4me3 promoter marks located within H3K27ac super-enhancer clusters in B cell non-Hodgkin lymphoma. One region found to undergo recurrent alteration upon transformation overlaps a super-enhancer affecting the expression of the PAX5/ZCCHC7 gene pair. ZCCHC7 encodes a subunit of the Trf4/5-Air1/2-Mtr4 polyadenylation-like complex and demonstrated copy number gain, chromosomal translocation and enhancer retargeting-mediated transcriptional upregulation upon lymphoma transformation. Consequently, lymphoma cells demonstrate nucleolar dysregulation via altered noncoding 5.8S ribosomal RNA processing. We find that a noncoding mutation acquired during lymphoma progression affects noncoding rRNA processing, thereby rewiring protein synthesis leading to oncogenic changes in the lymphoma proteome.

Systematic compositional analysis of exosomal extracellular vesicles produced by cells undergoing apoptosis, necroptosis and ferroptosis

Formation of extracellular vesicles (EVs) has emerged as a novel paradigm in cell-to-cell communication in health and disease. EVs are notably produced during cell death but it had remained unclear whether different modalities of regulated cell death (RCD) influence the biogenesis and composition of EVs. To this end, we performed a comparative analysis of steady-state (ssEVs) and cell death-associated EVs (cdEVs) following TNF-induced necroptosis (necEVs), anti-Fas-induced apoptosis (apoEVs), and ML162-induced ferroptosis (ferEVs) using the same cell line. For each RCD condition, we determined the biophysical and biochemical characteristics of the cell death-associated EVs (cdEVs), the protein cargo, and the presence of methylated ribosomal RNA. We found that the global protein content of all cdEVs was increased compared to steady-state EVs. Qualitatively, the isolated exosomal ssEVs and cdEVs, contained a largely overlapping protein cargo including some quantitative differences in particular proteins. All cdEVs were enriched for proteins involved in RNA splicing and nuclear export, and showed distinctive rRNA methylation patterns compared to ssEVs. Interestingly, necEVs and apoEVs, but strikingly not ferEVs, showed enrichment of proteins involved in ribosome biogenesis. Altogether, our work documents quantitative and qualitative differences between ssEVs and cdEVs.

Comprehensive identification of diverse ribosomal RNA modifications by targeted nanopore direct RNA sequencing and JACUSA2

Ribosomal RNAs are decorated by numerous post-transcriptional modifications whose exact roles in ribosome biogenesis, function, and human pathophysiology remain largely unknown. Here, we report a targeted direct rRNA sequencing approach involving a substrate selection step and demonstrate its suitability to identify differential modification sites in combination with the JACUSA2 software. We compared JACUSA2 to other tools designed for RNA modification detection and show that JACUSA2 outperforms other software with regard to detection of base modifications such as methylation, acetylation and aminocarboxypropylation. To illustrate its widespread usability, we applied our method to a collection of CRISPR-Cas9 engineered colon carcinoma cells lacking specific enzymatic activities responsible for particular rRNA modifications and systematically compared them to isogenic wild-type RNAs. Besides the numerous 2'-O methylated riboses and pseudouridylated residues, our approach was suitable to reliably identify differential base methylation and acetylation events. Importantly, our method does not require any prior knowledge of modification sites or the need to train complex models. We further report for the first time detection of human rRNA modifications by direct RNA-sequencing on Flongle flow cells, the smallest-scale nanopore flow cell available to date. The use of these smaller flow cells reduces RNA input requirements, making our workflow suitable for the analysis of samples with limited availability and clinical work.

Transcription-wide mapping of dihydrouridine reveals that mRNA dihydrouridylation is required for meiotic chromosome segregation

The epitranscriptome has emerged as a new fundamental layer of control of gene expression. Nevertheless, the determination of the transcriptome-wide occupancy and function of RNA modifications remains challenging. Here we have developed Rho-seq, an integrated pipeline detecting a range of modifications through differential modification-dependent rhodamine labeling. Using Rho-seq, we confirm that the reduction of uridine to dihydrouridine (D) by the Dus reductase enzymes targets tRNAs in E. coli and fission yeast. We find that the D modification is also present on fission yeast mRNAs, particularly those encoding cytoskeleton-related proteins, which is supported by large-scale proteome analyses and ribosome profiling. We show that the α-tubulin encoding mRNA nda2 undergoes Dus3-dependent dihydrouridylation, which affects its translation. The absence of the modification on nda2 mRNA strongly impacts meiotic chromosome segregation, resulting in low gamete viability. Applying Rho-seq to human cells revealed that tubulin mRNA dihydrouridylation is evolutionarily conserved.

Nucleolar stress controls mutant Huntington toxicity and monitors Huntington's disease progression

Transcriptional and cellular-stress surveillance deficits are hallmarks of Huntington's disease (HD), a fatal autosomal-dominant neurodegenerative disorder caused by a pathological expansion of CAG repeats in the Huntingtin (HTT) gene. The nucleolus, a dynamic nuclear biomolecular condensate and the site of ribosomal RNA (rRNA) transcription, is implicated in the cellular stress response and in protein quality control. While the exact pathomechanisms of HD are still unclear, the impact of nucleolar dysfunction on HD pathophysiology in vivo remains elusive. Here we identified aberrant maturation of rRNA and decreased translational rate in association with human mutant Huntingtin (mHTT) expression. The protein nucleophosmin 1 (NPM1), important for nucleolar integrity and rRNA maturation, loses its prominent nucleolar localization. Genetic disruption of nucleolar integrity in vulnerable striatal neurons of the R6/2 HD mouse model decreases the distribution of mHTT in a disperse state in the nucleus, exacerbating motor deficits. We confirmed NPM1 delocalization in the gradually progressing zQ175 knock-in HD mouse model: in the striatum at a presymptomatic stage and in the skeletal muscle at an early symptomatic stage. In Huntington's patient skeletal muscle biopsies, we found a selective redistribution of NPM1, similar to that in the zQ175 model. Taken together, our study demonstrates that nucleolar integrity regulates the formation of mHTT inclusions in vivo, and identifies NPM1 as a novel, readily detectable peripheral histopathological marker of HD progression.

A multi-scale map of cell structure fusing protein images and interactions

The cell is a multi-scale structure with modular organization across at least four orders of magnitude1. Two central approaches for mapping this structure-protein fluorescent imaging and protein biophysical association-each generate extensive datasets, but of distinct qualities and resolutions that are typically treated separately2,3. Here we integrate immunofluorescence images in the Human Protein Atlas4 with affinity purifications in BioPlex5 to create a unified hierarchical map of human cell architecture. Integration is achieved by configuring each approach as a general measure of protein distance, then calibrating the two measures using machine learning. The map, known as the multi-scale integrated cell (MuSIC 1.0), resolves 69 subcellular systems, of which approximately half are to our knowledge undocumented. Accordingly, we perform 134 additional affinity purifications and validate subunit associations for the majority of systems. The map reveals a pre-ribosomal RNA processing assembly and accessory factors, which we show govern rRNA maturation, and functional roles for SRRM1 and FAM120C in chromatin and RPS3A in splicing. By integration across scales, MuSIC increases the resolution of imaging while giving protein interactions a spatial dimension, paving the way to incorporate diverse types of data in proteome-wide cell maps.

Nopp140-chaperoned 2'-O-methylation of small nuclear RNAs in Cajal bodies ensures splicing fidelity

In this study, Bizarro et al. sought to understand the function and subcellular site of snRNA modification, and found that Cajal body (CB) localization of the protein Nopp140 is essential for concentration of small Cajal body-specific ribonucleoproteins (scaRNPs) in nuclear condensate and that phosphorylation by casein kinase 2 (CK2) at ∼80 serines targets Nopp140 to CBs. Nopp140 knockdown-mediated release of scaRNPs from CBs severely compromises 2′-O-methylation of spliceosomal snRNAs, identifying CBs as the site of scaRNP catalysis.

Spliceosomal small nuclear RNAs (snRNAs) are modified by small Cajal body (CB)-specific ribonucleoproteins (scaRNPs) to ensure snRNP biogenesis and pre-mRNA splicing. However, the function and subcellular site of snRNA modification are largely unknown. We show that CB localization of the protein Nopp140 is essential for concentration of scaRNPs in that nuclear condensate; and that phosphorylation by casein kinase 2 (CK2) at ∼80 serines targets Nopp140 to CBs. Transiting through CBs, snRNAs are apparently modified by scaRNPs. Indeed, Nopp140 knockdown-mediated release of scaRNPs from CBs severely compromises 2′-O-methylation of spliceosomal snRNAs, identifying CBs as the site of scaRNP catalysis. Additionally, alternative splicing patterns change indicating that these modifications in U1, U2, U5, and U12 snRNAs safeguard splicing fidelity. Given the importance of CK2 in this pathway, compromised splicing could underlie the mode of action of small molecule CK2 inhibitors currently considered for therapy in cholangiocarcinoma, hematological malignancies, and COVID-19.

Nopp140-chaperoned 2'-O-methylation of small nuclear RNAs in Cajal bodies ensures splicing fidelity

Spliceosomal small nuclear RNAs (snRNAs) are modified by small Cajal body (CB) specific ribonucleoproteins (scaRNPs) to ensure snRNP biogenesis and pre-mRNA splicing. However, the function and subcellular site of snRNA modification are largely unknown. We show that CB localization of the protein Nopp140 is essential for concentration of scaRNPs in that nuclear condensate; and that phosphorylation by casein kinase 2 (CK2) at some 80 serines targets Nopp140 to CBs. Transiting through CBs, snRNAs are apparently modified by scaRNPs. Indeed, Nopp140 knockdown-mediated release of scaRNPs from CBs severely compromises 2'-O-methylation of spliceosomal snRNAs, identifying CBs as the site of scaRNP catalysis. Additionally, alternative splicing patterns change indicating that these modifications in U1, U2, U5, and U12 snRNAs safeguard splicing fidelity. Given the importance of CK2 in this pathway, compromised splicing could underlie the mode of action of small molecule CK2 inhibitors currently considered for therapy in cholangiocarcinoma, hematological malignancies, and COVID-19.

The catalytic activity of the translation termination factor methyltransferase Mtq2-Trm112 complex is required for large ribosomal subunit biogenesis

The Mtq2-Trm112 methyltransferase modifies the eukaryotic translation termination factor eRF1 on the glutamine side chain of a universally conserved GGQ motif that is essential for release of newly synthesized peptides. Although this modification is found in the three domains of life, its exact role in eukaryotes remains unknown. As the deletion of MTQ2 leads to severe growth impairment in yeast, we have investigated its role further and tested its putative involvement in ribosome biogenesis. We found that Mtq2 is associated with nuclear 60S subunit precursors, and we demonstrate that its catalytic activity is required for nucleolar release of pre-60S and for efficient production of mature 5.8S and 25S rRNAs. Thus, we identify Mtq2 as a novel ribosome assembly factor important for large ribosomal subunit formation. We propose that Mtq2-Trm112 might modify eRF1 in the nucleus as part of a quality control mechanism aimed at proof-reading the peptidyl transferase center, where it will subsequently bind during translation termination.

The ribosomal RNA m5C methyltransferase NSUN-1 modulates healthspan and oogenesis in Caenorhabditis elegans

Our knowledge about the repertoire of ribosomal RNA modifications and the enzymes responsible for installing them is constantly expanding. Previously, we reported that NSUN-5 is responsible for depositing m⁵C at position C2381 on the 26S rRNA in Caenorhabditis elegans. Here, we show that NSUN-1 is writing the second known 26S rRNA m⁵C at position C2982. Depletion of nsun-1 or nsun-5 improved thermotolerance and slightly increased locomotion at midlife, however, only soma-specific knockdown of nsun-1 extended lifespan. Moreover, soma-specific knockdown of nsun-1 reduced body size and impaired fecundity, suggesting non-cell-autonomous effects. While ribosome biogenesis and global protein synthesis were unaffected by nsun-1 depletion, translation of specific mRNAs was remodeled leading to reduced production of collagens, loss of structural integrity of the cuticle, and impaired barrier function. We conclude that loss of a single enzyme required for rRNA methylation has profound and highly specific effects on organismal development and physiology.

HydraPsiSeq: a method for systematic and quantitative mapping of pseudouridines in RNA

Developing methods for accurate detection of RNA modifications remains a major challenge in epitranscriptomics. Next-generation sequencing-based mapping approaches have recently emerged but, often, they are not quantitative and lack specificity. Pseudouridine (ψ), produced by uridine isomerization, is one of the most abundant RNA modification. ψ mapping classically involves derivatization with soluble carbodiimide (CMCT), which is prone to variation making this approach only semi-quantitative. Here, we developed 'HydraPsiSeq', a novel quantitative ψ mapping technique relying on specific protection from hydrazine/aniline cleavage. HydraPsiSeq is quantitative because the obtained signal directly reflects pseudouridine level. Furthermore, normalization to natural unmodified RNA and/or to synthetic in vitro transcripts allows absolute measurements of modification levels. HydraPsiSeq requires minute amounts of RNA (as low as 10-50 ng), making it compatible with high-throughput profiling of diverse biological and clinical samples. Exploring the potential of HydraPsiSeq, we profiled human rRNAs, revealing strong variations in pseudouridylation levels at ∼20-25 positions out of total 104 sites. We also observed the dynamics of rRNA pseudouridylation throughout chondrogenic differentiation of human bone marrow stem cells. In conclusion, HydraPsiSeq is a robust approach for the systematic mapping and accurate quantification of pseudouridines in RNAs with applications in disease, aging, development, differentiation and/or stress response.

Synergistic defects in pre-rRNA processing from mutations in the U3-specific protein Rrp9 and U3 snoRNA

U3 snoRNA and the associated Rrp9/U3-55K protein are essential for 18S rRNA production by the SSU-processome complex. U3 and Rrp9 are required for early pre-rRNA cleavages at sites A0, A1 and A2, but the mechanism remains unclear. Substitution of Arg 289 in Rrp9 to Ala (R289A) specifically reduced cleavage at sites A1 and A2. Surprisingly, R289 is located on the surface of the Rrp9 β-propeller structure opposite to U3 snoRNA. To understand this, we first characterized the protein-protein interaction network of Rrp9 within the SSU-processome. This identified a direct interaction between the Rrp9 β-propeller domain and Rrp36, the strength of which was reduced by the R289A substitution, implicating this interaction in the observed processing phenotype. The Rrp9 R289A mutation also showed strong synergistic negative interactions with mutations in U3 that destabilize the U3/pre-rRNA base-pair interactions or reduce the length of their linking segments. We propose that the Rrp9 β-propeller and U3/pre-rRNA binding cooperate in the structure or stability of the SSU-processome. Additionally, our analysis of U3 variants gave insights into the function of individual segments of the 5′-terminal 72-nt sequence of U3. We interpret these data in the light of recently reported SSU-processome structures.

The 18S ribosomal RNA m6 A methyltransferase Mettl5 is required for normal walking behavior in Drosophila

RNA modifications have recently emerged as an important layer of gene regulation. N6-methyladenosine (m⁶ A) is the most prominent modification on eukaryotic messenger RNA and has also been found on noncoding RNA, including ribosomal and small nuclear RNA. Recently, several m⁶ A methyltransferases were identified, uncovering the specificity of m⁶ A deposition by structurally distinct enzymes. In order to discover additional m⁶ A enzymes, we performed an RNAi screen to deplete annotated orthologs of human methyltransferase-like proteins (METTLs) in Drosophila cells and identified CG9666, the ortholog of human METTL5. We show that CG9666 is required for specific deposition of m⁶ A on 18S ribosomal RNA via direct interaction with the Drosophila ortholog of human TRMT112, CG12975. Depletion of CG9666 yields a subsequent loss of the 18S rRNA m⁶ A modification, which lies in the vicinity of the ribosome decoding center; however, this does not compromise rRNA maturation. Instead, a loss of CG9666-mediated m⁶ A impacts fly behavior, providing an underlying molecular mechanism for the reported human phenotype in intellectual disability. Thus, our work expands the repertoire of m⁶ A methyltransferases, demonstrates the specialization of these enzymes, and further addresses the significance of ribosomal RNA modifications in gene expression and animal behavior.

The reverse transcription signature of N-1-methyladenosine in RNA-Seq is sequence dependent

The combination of Reverse Transcription (RT) and high-throughput sequencing has emerged as a powerful combination to detect modified nucleotides in RNA via analysis of either abortive RT-products or of the incorporation of mismatched dNTPs into cDNA. Here we simultaneously analyze both parameters in detail with respect to the occurrence of N-1-methyladenosine (m¹A) in the template RNA. This naturally occurring modification is associated with structural effects, but it is also known as a mediator of antibiotic resistance in ribosomal RNA. In structural probing experiments with dimethylsulfate, m¹A is routinely detected by RT-arrest. A specifically developed RNA-Seq protocol was tailored to the simultaneous analysis of RT-arrest and misincorporation patterns. By application to a variety of native and synthetic RNA preparations, we found a characteristic signature of m¹A, which, in addition to an arrest rate, features misincorporation as a significant component. Detailed analysis suggests that the signature depends on RNA structure and on the nature of the nucleotide 3′ of m¹A in the template RNA, meaning it is sequence dependent. The RT-signature of m¹A was used for inspection and confirmation of suspected modification sites and resulted in the identification of hitherto unknown m¹A residues in trypanosomal tRNA.

The human 18S rRNA base methyltransferases DIMT1L and WBSCR22-TRMT112 but not rRNA modification are required for ribosome biogenesis

An evolutionarily conserved quality control in ribosome biogenesis reveals that two human rRNA base methyltransferases associated with cell differentiation and cancer but, surprisingly, not their RNA-modifying activity are required for small ribosomal subunit biogenesis.

At the heart of the ribosome lie rRNAs, whose catalytic function in translation is subtly modulated by posttranscriptional modifications. In the small ribosomal subunit of budding yeast, on the 18S rRNA, two adjacent adenosines (A1781/A1782) are N⁶-dimethylated by Dim1 near the decoding site, and one guanosine (G1575) is N⁷-methylated by Bud23-Trm112 at a ridge between the P- and E-site tRNAs. Here we establish human DIMT1L and WBSCR22-TRMT112 as the functional homologues of yeast Dim1 and Bud23-Trm112. We report that these enzymes are required for distinct pre-rRNA processing reactions leading to synthesis of 18S rRNA, and we demonstrate that in human cells, as in budding yeast, ribosome biogenesis requires the presence of the modification enzyme rather than its RNA-modifying catalytic activity. We conclude that a quality control mechanism has been conserved from yeast to human by which binding of a methyltransferase to nascent pre-rRNAs is a prerequisite to processing, so that all cleaved RNAs are committed to faithful modification. We further report that 18S rRNA dimethylation is nuclear in human cells, in contrast to yeast, where it is cytoplasmic. Yeast and human ribosome biogenesis thus have both conserved and distinctive features.

Structural and functional studies of Bud23-Trm112 reveal 18S rRNA N7-G1575 methylation occurs on late 40S precursor ribosomes

The eukaryotic small ribosomal subunit carries only four ribosomal (r) RNA methylated bases, all close to important functional sites. N(7)-methylguanosine (m(7)G) introduced at position 1575 on 18S rRNA by Bud23-Trm112 is at a ridge forming a steric block between P- and E-site tRNAs. Here we report atomic resolution structures of Bud23-Trm112 in the apo and S-adenosyl-L-methionine (SAM)-bound forms. Bud23 and Trm112 interact through formation of a β-zipper involving main-chain atoms, burying an important hydrophobic surface and stabilizing the complex. The structures revealed that the coactivator Trm112 undergoes an induced fit to accommodate its methyltransferase (MTase) partner. We report important structural similarity between the active sites of Bud23 and Coffea canephora xanthosine MTase, leading us to propose and validate experimentally a model for G1575 coordination. We identify Bud23 residues important for Bud23-Trm112 complex formation and recruitment to pre-ribosomes. We report that though Bud23-Trm112 binds precursor ribosomes at an early nucleolar stage, m(7)G methylation occurs at a late step of small subunit biogenesis, implying specifically delayed catalytic activation. Finally, we show that Bud23-Trm112 interacts directly with the box C/D snoRNA U3-associated DEAH RNA helicase Dhr1 supposedly involved in central pseudoknot formation; this suggests that Bud23-Trm112 might also contribute to controlling formation of this irreversible and dramatic structural reorganization essential to overall folding of small subunit rRNA. Our study contributes important new elements to our understanding of key molecular aspects of human ribosomopathy syndromes associated with WBSCR22 (human Bud23) malfunction.

TOR regulates the subcellular distribution of DIM2, a KH domain protein required for cotranscriptional ribosome assembly and pre-40S ribosome export

Eukaryotic ribosome synthesis is a highly dynamic process that involves the transient association of scores of trans-acting factors to nascent pre-ribosomes. Many ribosome synthesis factors are nucleocytoplasmic shuttling proteins that engage the assembly pathway at early nucleolar stages and escort pre-ribosomes to the nucleoplasm and/or the cytoplasm. Here, we report that two 40S ribosome synthesis factors, the KH-domain protein DIM2 and the HEAT-repeats/Armadillo-domain and export factor RRP12, are nucleolar restricted upon nutritional, osmotic, and oxidative stress. Nucleolar entrapment of DIM2 and RRP12 was triggered by rapamycin treatment and was under the strict control of the target of rapamycin (TOR) signaling cascade. DIM2 binds pre-rRNAs directly through its KH domain at the 5'-end of ITS1 (D-A(2) segment) and, consistent with its requirements in early nucleolar pre-rRNA processing, is required for efficient cotranscriptional ribosome assembly. The substitution of a single and highly conserved amino acid (G207A) within the KH motif is sufficient to inhibit pre-rRNA processing in a fashion similar to genetic depletion of DIM2. DIM2 carries an evolutionarily conserved putative nuclear export sequence (NES) at its carboxyl-terminal end that is required for efficient pre-40S ribosome export. Strikingly, DIM2 and RRP12 are both involved in the nucleocytoplasmic translocation of pre-ribosomes, suggesting that this step in the ribosome assembly pathway has been selected as a regulatory target for the TOR pathway.

Characterization of the type III secretion locus of Bordetella pertussis

Multiple sequence comparisons of proteins of the LcrD/FlbF family allowed the design of primers that specifically amplify sequences coding for type III secretion components. Amplification of Bordetella pertussis DNA with these primers yielded a fragment that was further used as a probe for screening a genomic library. The nucleotide sequence of a positive clone revealed a 2100-bp gene, called bcrD, which specifies a 75-kDa polypeptide homologous to the Yersinia LcrD protein. Chromosome walking allowed the characterization of a 35-kb DNA segment that contains the entire locus and flanking housekeeping genes. The B. pertussis type III secretion locus consists of more than 30 open reading frames (ORFs), most of which are identical to annotated genes of Bordetella spp and share similarities with known type III secretion genes of related bacteria. In order to assess the function of this locus, we engineered a bcrD null mutant. However, none of the tested phenotypes, such as protein secretion, cellular invasion, cytotoxicity or mouse lung colonization, differentiated the mutant from its parental strain. Studies of bcrD and bscN expressions indicated that, under our experimental conditions, these genes are not expressed in vitro. Restriction analyses on pulsed-field gel electrophoresis allowed the type III locus mapping at coordinate position 1,590 kb on the Tohama I strain chromosome.