Next Article in Journal
Regulation of Notch1 Signalling by Long Non-Coding RNAs in Cancers and Other Health Disorders
Next Article in Special Issue
Comprehensive Pan-Cancer Mutation Density Patterns in Enhancer RNA
Previous Article in Journal
Catalyst Design through Grafting of Diazonium Salts—A Critical Review on Catalyst Stability
Previous Article in Special Issue
Functional Relationships between Long Non-Coding RNAs and Estrogen Receptor Alpha: A New Frontier in Hormone-Responsive Breast Cancer Management
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Depletion of SNORA33 Abolishes ψ of 28S-U4966 and Affects the Ribosome Translational Apparatus

by
Alzbeta Chabronova
1,
Guus van den Akker
1,
Bas A. C. Housmans
1,
Marjolein M. J. Caron
1,
Andy Cremers
1,
Don A. M. Surtel
1,
Mandy J. Peffers
2,
Lodewijk W. van Rhijn
1,
Virginie Marchand
3,
Yuri Motorin
3,4 and
Tim J. M. Welting
1,5,*
1
Laboratory for Experimental Orthopedics, Department of Orthopedic Surgery, Maastricht University, 6229 HX Maastricht, The Netherlands
2
Institute of Life Course and Medical Sciences, University of Liverpool, Liverpool L8 7TX, UK
3
UAR2008 IBSLor CNRS-INSERM-Université de Lorraine, F54000 Nancy, France
4
UMR7365 IMOPA, CNRS-Université de Lorraine, F54000 Nancy, France
5
Laboratory for Experimental Orthopedics, Department of Orthopedic Surgery, Maastricht University Medical Center+ (MUMC+), 6229 HX Maastricht, The Netherlands
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(16), 12578; https://doi.org/10.3390/ijms241612578
Submission received: 5 July 2023 / Revised: 4 August 2023 / Accepted: 7 August 2023 / Published: 8 August 2023
(This article belongs to the Special Issue Non-coding RNA in Physiology and Pathophysiology)

Abstract

:
Eukaryotic ribosomes are complex molecular nanomachines translating genetic information from mRNAs into proteins. There is natural heterogeneity in ribosome composition. The pseudouridylation (ψ) of ribosomal RNAs (rRNAs) is one of the key sources of ribosome heterogeneity. Nevertheless, the functional consequences of ψ-based ribosome heterogeneity and its relevance for human disease are yet to be understood. Using HydraPsiSeq and a chronic disease model of non-osteoarthritic primary human articular chondrocytes exposed to osteoarthritic synovial fluid, we demonstrated that the disease microenvironment is capable of instigating site-specific changes in rRNA ψ profiles. To investigate one of the identified differential rRNA ψ sites (28S-ψ4966), we generated SNORA22 and SNORA33 KO SW1353 cell pools using LentiCRISPRv2/Cas9 and evaluated the ribosome translational capacity by 35S-Met/Cys incorporation, assessed the mode of translation initiation and ribosomal fidelity using dual luciferase reporters, and assessed cellular and ribosomal proteomes by LC-MS/MS. We uncovered that the depletion of SNORA33, but not SNORA22, reduced 28S-ψ4966 levels. The resulting loss of 28S-ψ4966 affected ribosomal protein composition and function and led to specific changes in the cellular proteome. Overall, our pioneering findings demonstrate that cells dynamically respond to disease-relevant changes in their environment by altering their rRNA pseudouridylation profiles, with consequences for ribosome function and the cellular proteome relevant to human disease.

1. Introduction

The eukaryotic ribosome represents the core cellular translation machinery. This complex cellular nanomachine consists of ribosomal proteins and ribosomal RNAs (rRNAs) [1]. Ribosomal RNAs are heavily post-transcriptionally modified [2]. Up to 228 modified nucleotides (post-transcriptional modifications, PTMs) have been reported in human 18S, 5.8S, and 28S rRNAs and this number keeps growing [2,3,4]. Pseudouridylation, with up to 104 already reported rRNA sites, is one of the most abundant rRNA modifications (exceeded only by 2′-O-methylation with ~110 sites) [3,4]. The pseudouridylation (ψ) of rRNA nucleotides is carried out by snoRNA–protein complexes (snoRNPs) composed of H/ACA snoRNAs and four core protein co-factors, NHP2, NOP10, GAR1, and Dyskerin (DKC1) [5,6,7]. The box H/ACA snoRNAs (SNORAs) site-directionally guide the pseudouridylation of specific nucleotides by base-pairing with the target rRNA sequence, allowing the pseudouridine synthase DKC1 to catalyze the isomerization of the target uridine to pseudouridine [5,6,7]. Many rRNA ψ sites are conserved among prokaryotic and eukaryotic species and are clustered in functionally important ribosomal regions that interact with tRNAs and mRNAs, including the peptidyl transferase center (PTC), the polypeptide exit tunnel, or the intersubunit bridge [2,8]. During evolution and with growing organismal complexity, the rRNA ψ content increased, presumably to reinforce the structural and functional stability of ribosomes [2,8]. Higher eukaryotes, especially mammals, have several additional rRNA ψ sites, usually located in the exterior regions of the ribosome [2,3].
PTMs of rRNAs are a major source of ribosome heterogeneity, with implications for ribosome functional specialization [9,10,11,12,13,14,15,16]. The first link between the ψ of rRNAs and human disease was made when mutations in the DKC1 gene were found to be implicated in the pathogenesis of X-linked dyskeratosis congenita (X-DC) [17,18]. X-DC is a rare, inherited, multisystemic syndrome characterized by bone marrow failure, reticular skin pigmentation, and oral leukoplakia, and it is associated with a high risk of cancer [19]. Studies using cells from X-DC patients, hypomorphic Dkc1 mutant mice, and a yeast strain with mutated Cbf5p (the yeast homologue of DKC1) reported deregulated translation, especially in connection to internal ribosome entry site (IRES)-mediated translation initiation [12,13]. A loss of Ψ in helix 69 of 28S rRNA, a highly conserved hotspot for rRNA PTMs that interacts with A- and P-site tRNAs and forms the intersubunit ribosome bridge, resulted in impaired growth, reduced amino acid incorporation, and defective ribosome subunits association, as well as faulty reading frame maintenance and stop codon recognition in yeast [20,21,22,23]. DKC1 is highly expressed in a variety of cancer tissues and high DKC1 expression correlates with poor prognosis [24]. Furthermore, the expression of several H/ACA snoRNAs is altered in cancer [25] and oncogene RAS was shown to be involved in the regulation of their expression (e.g., SNORA23, SNORA24, SNORA26, SNORA48, and SNORA6) [26]. A high-throughput, deep-sequencing-based method (HydraPsiSeq) for systematic rRNA ψ mapping and quantification was developed recently and allowed the first comprehensive analysis of ψ-based ribosome heterogeneity in human cells and the identification of a stable and variable group of ψ modifications [3]. Nevertheless, the roles of rRNA ψ-based ribosome heterogeneity in human disease are largely unknown.
Osteoarthritis (OA) is the most common degenerative joint disease and a major cause of pain and disability in adults [27]. Osteoarthritis represents an active disease process that affects the entire joint and all its tissues (e.g., articular cartilage, synovium, or subchondral bone) [28,29]. Disturbances in the homeostasis of joint tissues create a pathological, degenerative joint microenvironment. Synovial fluid located in the joint cavity is a plasma ultrafiltrate enriched in molecules (e.g., proteins, lipids, or metabolites) secreted by joint tissues, and as such it reflects systemic and local pathological changes that occur during osteoarthritis [30,31,32]. Chondrocytes, cells residing in the cartilage extracellular matrix, rely on synovial fluid for their nourishment and survival [33]. The major hallmark of osteoarthritis is the degeneration of the articular cartilage, which is a result of the “activation” of quiescent chondrocytes residing in the cartilage extracellular matrix [28]. Changes in chondrocyte activity and cellular phenotype are fueled by alterations in their protein expression programs [34]. Aberrations of ribosome biogenesis and activity have been reported in osteoarthritic chondrocytes [35]. However, very little is known about the role of ribosome heterogeneity and specialization in the pathobiology of this chronic disease. Our group recently published data demonstrating that osteoarthritic synovial fluid instigates site-specific changes in the 2′-O-me rRNA profiles of human chondrocytes and that these changes have consequences for ribosome function and translation [36]. In this work, we used primary human articular chondrocytes exposed to the osteoarthritic synovial fluid to investigate the effect of a chronic disease microenvironment on rRNA ψ profiles.

2. Results

2.1. A Chronic Disease Microenvironment Provokes Site-Specific Changes in rRNA ψ Profiles

To explore the effect of a chronic disease microenvironment on rRNA ψ profiles, we employed an in vitro osteoarthritis disease model in which we exposed human primary chondrocytes to osteoarthritic synovial fluid (OA-SF) for 14 days and performed rRNA ψ HydraPsiSeq profiling (Figure 1A) [3]. The 14-day-long exposure of chondrocytes to OA_SF was chosen based on our previous data. Our time course experiment showed that the treatment with OA-SF triggered OA-related changes in chondrocytes’ gene expression and signaling from day 3 onwards, and these changes appeared to be rather stable until day 13 (the latest time point that we measured) [37]. In this study, we chose a time point of 14 days due to expected slow turnover of ribosomes when compared to gene expression changes. Ribosomal RNA ψ is considered irreversible, and the average half-life of rRNA in the cytoplasm of mouse articular cartilage is 5–7 days [38]. There are no conclusive data on the rRNA half-life in primary human chondrocytes; however, it could be expected that 14 days of OA-SF treatment would allow for the complete turnover of rRNAs (and ribosomes) and provide sufficient time for changes in rRNA ψ profiles to take place. We detected 102 pseudouridylated rRNA nucleotides (two sites not detected). Of these, we observed 41 in 18S, 2 in 5.8S, and 59 in 28S (Figure S1 and Table S1). We identified seven differentially pseudouridylated sites in the OA-SF-treated cells, three sites in 18S (18S-ψ36, 18S-ψ210, 18S-ψ918), and four sites in 28S (28S-ψ1766, 28S-ψ3801, 28S-ψ4606, and 28S-ψ4966) (Figure 1B). In particular, 18S-ψ36 and 18S-ψ210 are both located in domain 5’ of the small human ribosomal subunit, specifically within helix 4 (18S-ψ36) and helix 9es3a (18S-ψ210), while 18S-ψ918 is located in helix 21es6d/22 of domain C. The remaining four differentially modified sites are located within the large human ribosomal subunit, domain II H38a (28S-ψ1766), domain 0 H64/70 (28S-ψ3801), and domains VI H69 (28S-ψ4606) and H101 (28S-ψ4966). All seven differentially ψ rRNA sites exhibited decreased modification levels in response to OA-SF. These data are the first evidence that a chronic disease microenvironment can induce changes in rRNA ψ profiles. There was no apparent clustering of these differential osteoarthritis-sensitive sites in terms of ribosomal helices or domains. Interestingly, decreased modification levels of 28S-ψ4966 were previously also reported in the cells of patients with familial X-DC [2]. Nevertheless, the relevance of 28S-U4966 and its conversion to ψ for ribosome function has not been investigated. Therefore, we undertook functional validation for this modification in our further work.

2.2. Depletion of SNORA33 but Not SNORA22 Reduces 28S-ψ4966 Levels

To functionally investigate 28S-ψ4966 and its role(s) in ribosome function and osteoarthritis pathobiology, we depleted chondrocytic SW1353 cells of snoRNAs guiding this specific pseudouridylation. Two intron-encoded H/ACA snoRNAs, SNORA22 and SNORA33, were previously predicted to guide the ψ of 28S-U4966 (Figure 2A,E) [39]. We used CRISPR/Cas9 gene editing and a double sgRNA targeting approach (Figure 2B,F) to generate SNORA22 and SNORA33 knockout (KO) cell pools along with a GFP-targeting CRISPR control cell pool. The efficiency of the gene editing was verified using the surveyor nuclease assay (Figure S2), and genomic DNA sequencing confirmed a 32 basepair (bp) deletion in the SNORA22 gene and 103 bp deletion in the SNORA33 gene (Figure S2C). Despite the complete loss of SNORA22 expression, we did not measure any decrease in 28S-ψ4966 levels (Figure 2C,D). We also did not measure any effects of SNORA22 deletion on its host gene (CCT6P1) expression or splicing. On the other hand, a 55% decrease in SNORA33 expression resulted in the ablation of ψ at 28S-U4966 (Figure 2G,H). These results imply that the ψ of 28S-U4966 is in fact guided only by SNORA33. Looking at the gene expression and splicing of the SNORA33 host gene RPS12, we noted a slight increase in RPS12 coding sequence (CDS) levels, but no effect on the splicing of exons 5 and 6 enclosing the SNORA33-hosting intron region (Figure 2G). In conclusion, our data demonstrate that the ψ of 28S-U4966 is guided by SNORA33. Besides the ablation of 28S-ψ4966 in SNORA33-depleted cells, we also measured a modest, but statistically significant, increase in ψ levels at three other sites, 28S-ψ4975, 28S-ψ3801, and 28S-ψ3741 (Figure 3A and Figure S4, Table S2). By mapping the sites with differential ψ status in SNORA33-depleted cells within the 3D human ribosome, we observed that these sites clustered in domains IV and VI of the large ribosomal subunit (Figure 3B). While 28S-U3801 and 28S-U3741 are located in domain IV, 28S-U4975 and 28S-U4966 are situated within domain VI. Interestingly, the ψ levels of 28S-U3801 were also regulated (decreased) in our original in vitro chronic disease model on which rRNA ψ HydraPsiSeq profiling was applied (Figure 1A). To corroborate the relevance of SNORA33, and by extension the 28S-U4966 pseudouridylation levels in osteoarthritis, we measured the SNORA33 expression levels in freshly isolated healthy and osteoarthritic human primary chondrocytes. Results showed that SNORA33 expression was indeed significantly decreased in osteoarthritic chondrocytes compared to healthy controls (Figure S3).

2.3. SNORA33 Depletion Affects Ribosome Composition, Function, and the Cellular Proteome

After establishing the SNORA33 KO cell pools, we proceeded with the functional validation of the relevance of 28S-ψ4966, by evaluating the ribosome composition and function in SNORA33-depleted cells. First, we probed whether the depletion of SNORA33 affected overall cellular fitness, by assessing the proliferation rate, total protein translation, and polysomal distribution. We did not detect any differences in these characteristics in the SNORA33-depleted cells compared to CRISPR controls (Figure 4). Next, we examined the protein compositions of ribosomes devoid of 28S-ψ4966. We used low salt procedures [40] for the ultracentrifugation of cytoplasmic extracts through sucrose cushions to isolate ribosomes as well as ribosome-associated factors and accessory proteins (Figure 5A). Analysis of the ribosomal proteome identified 1102 proteins, including 79 core ribosomal proteins, many ribosome biogenesis factors, and a large number of translation initiation, elongation, and termination factors. We also identified proteins involved in mRNA processing, stability, and decay and factors assisting protein folding and modification (Figure 5B,C; analyzed using the Reactome [41]). Seventy-four proteins showed differential abundance in isolated ribosomes lacking 28S-ψ4966 (Figure 5D, Table S3). Among the differentially abundant proteins were several ribosome biogenesis factors (e.g., NOB1, DNAJA1, PDS5A), translation initiation factors EIF3F and EIF3G, elongation factor EEF1D, eukaryotic peptide chain release factor 1 (ETF1), and proteins involved in ribosome quality control, LTN1 and NEMF (Table S4, Figure 5E). A single core ribosomal protein, RPS28, showed increased association with ribosomes depleted of 28S-ψ4966 (Figure 5E).
To examine whether ribosome heterogeneity based on the identified differential ribosomal protein composition and 28S ψ affects ribosome function, we probed the mode of translation initiation and translation fidelity in the SNORA33 KO cell pool using dual-luciferase reporter assays. We demonstrated the increased capacity of 28S-ψ4966-depleted ribosomes to initiate translation from CrPV and HCV IRES elements (Figure 6A). Translation accuracy was also affected, as evidenced by a greater incidence of -1 frameshift activity. Subsequently, we analyzed the cellular proteome of SNORA33-depleted cells. We detected 2377 cellular proteins, of which 28 proteins were differentially expressed (DE) in SNORA33-depleted cells compared to CRISPR controls (6 upregulated and 22 downregulated, Figure 6B, Table S5). The most upregulated protein was a subunit of the RNA polymerase I complex, POLR1A, which catalyzes the synthesis of 18S, 5.8S, and 28S rRNAs (47S transcript) [42]. Several inflammation-related factors were among the DE proteins, including programmed cell death protein 11 (PDCD11, also called NFBP) and α-2-macroglobulin (A2M) (Figure 6C). Importantly, we did not detect differential expression of the SNORA33 host gene RPS12. Of the proteins that were found to be differentially associated with the ribosomes of SNORA33-depleted cells, none was found to be DE in the cellular proteome. This implies that the ribosomal proteome data are reflective of distinct associations of these proteins with ribosomes, rather than their expression. Combined, our results demonstrate that SNORA33 depletion and the loss of 28S-ψ4966 affect ribosome composition and function, with consequences for the cellular proteome.

3. Discussion

The recent development of a high-throughput method for systematic ψ mapping, HydraPsiSeq, enabled a comprehensive analysis of ψ-based ribosome heterogeneity in human cells [3]. Nevertheless, the driving force behind the differential ψ of rRNA nucleotides and their functional consequences for human disease are largely unknown. We used an in vitro disease model to demonstrate that the rRNA pseudouridylation profiles of human primary cells are sensitive to disease-related changes in their microenvironment. We identified seven rRNA sites whose ψ-modification levels significantly decreased in response to a chronic osteoarthritis disease microenvironment. We selected one of them, 28S-ψ4966, for more detailed functional validation. Out of two snoRNAs (SNORA33 and SNORA22) that have been predicted to guide the ψ conversion of 28S-U4966, we determined that only SNORA33 is in fact responsible for this PTM in SW1353 cells. Depleting chondrocytic cells of SNORA33 using a CRISPR/Cas approach, we generated a cell pool devoid of 28S-ψ4966 and proceeded with the ribosome- and translation-focused experiments. The ablation of SNORA33 and thus 28S-ψ4966 affected the ribosomal protein composition, as evidenced by the differential association of core ribosomal protein RPS28 and many ribosome-associated factors related to translation initiation and elongation. Consequently, we detected alterations in ribosome translational modus and fidelity, as well as the differential expression of several inflammation-related proteins in the cellular proteome of the SNORA33-depleted SW1353 cells.
Osteoarthritis is a degenerative chronic joint disease. Synovial fluid is a major determinant of the cartilage (and chondrocytes) microenvironment that well reflects the pathological processes occurring in the joint during osteoarthritis [30,31,32]. Recently, we performed a direct comparison of non-osteoarthritic and osteoarthritic synovial fluid and we identified critical differences in their compositions and effects on chondrocytes [43]. We showed that OA-SF is strongly enriched in inflammatory mediators, growth factors, and disease-associated molecular patterns (DAMPs) and provokes differential chondrocyte signaling, inflammatory responses, and chondrocyte phenotype alterations relevant to osteoarthritis. This demonstrates that chondrocytes exposed to osteoarthritic synovial fluid represent a relevant model to study how the local joint disease microenvironment affects the pathomolecular characteristics of chondrocytes [36,37,43]. In the current study, we did not use non-OA-SF as a control. Instead, we utilized a culture medium supplemented with 0.9% NaCl in the same v/v concentration as OA-SF. This approach was employed due to the limited availability of human non-OA-SF. Importantly, the NaCl control was previously shown to be suitable for this type of comparison [36,37].
Recently published data from our group demonstrated that osteoarthritic synovial fluid instigates site-specific changes (18S-Gm1447, 5.8S-Um14, 28S-Am3739, 28S-Am3846, and 28S-Um4590) in the 2′-O-me rRNA profiles of human primary chondrocytes, with consequences for their translation [36]. Here, we extended the dataset of osteoarthritis-sensitive rRNA PTMs for seven rRNA ψ sites (18S-ψ36, 18S-ψ210, 18S-ψ918, 28S-ψ1766, 28S-ψ3801, 28S-ψ4606, and 28S-ψ4966). These data are valuable not only for the osteoarthritis field but also for research on ribosome heterogeneity in general. Thus far, studies on rRNA PTM-based ribosome heterogeneity in humans have focused mainly on measuring differences between different cell lines, cell types, and healthy and diseased tissues [3,11], while only a handful of studies investigated the regulatory mechanisms behind the differential rRNA PTMs [16]. Data on rRNA ψ are particularly limited as a high-throughput method for systematic ψ quantification was only recently developed in 2020 [3]. Our data represent a pioneering study that analyzed the full ψ signature of human primary cells in response to a chronic disease microenvironment. The magnitude of changes in rRNA ψ profiles in response to chronic disease microenvironment was relatively small, ranging from 5.8 to 12.1%. However, these changes were significant and consistent across all five individual donors. Although further research is needed to evaluate the relevance of this scale of ribosome rRNA heterogeneity for human diseases, we believe that it might contribute to pathological changes in articular cartilage in osteoarthritis [35].
Lower modification levels of one of the osteoarthritis-sensitive rRNA ψ sites were previously measured in fibroblasts and B-lymphoblastoid cell lines derived from X-DC patients when compared to profiles of their unaffected relatives [2]. In fact, only two sites exhibited decreased ψ levels in the rRNAs of X-DC cells, 28S-ψ4331 and the here identified osteoarthritis-sensitive 28S-ψ4966. For this reason, we decided to follow up on this modification in further translation-focused experiments.
Proteomic analysis of isolated ribosomes depleted of 28S-ψ4966 uncovered differential stoichiometry of several translation factors and a modest difference in the association of core ribosomal protein RPS28, which is located relatively far from 28S-ψ4966. RPS28 is part of the eukaryotic translation initiation complex and interacts with the 5’ UTRs of mRNAs [44]. Furthermore, it has been suggested that ribosomes containing RPS28 variants more efficiently translate mRNAs with certain 5’ or 3’ UTR features [45]. Interestingly, mutations in RPS28 have been implicated in Diamond–Blackfan anemia, a ribosomopathy characterized by red blood cell aplasia, skeletal anomalies, and short stature [46]. We speculate that the change in the association of RPS28 with the ribosomes of SNORA33-depleted cells might have contributed to the changes in ribosome function and cellular proteome that we observed. Nevertheless, a connection between SNORA33 depletion and the stoichiometry of RPS28 remains elusive. Besides RPS28, we also found a differential association of multiple translation factors regulating various steps of translation, including translation initiation and elongation but also termination, as well as ribosome stalling and ribosome recycling. These unique data accentuate the mutually dependent relationships between rRNA PTMs and protein-based features of ribosome heterogeneity.
IRES-mediated translation represents an alternative route of translation initiation that bypasses the cap-dependent mRNA scanning and directly binds the 40S ribosomal subunit to the vicinity of the start codon to initiate the translation in a cap-independent manner [47]. Although first discovered in viruses, the IRES-mediated translation of cellular mRNAs is becoming increasingly recognized as a mechanism activated under cellular stress [48]. The regulation behind the IRES-mediated translation in cellular systems is mostly incompletely understood. However, rRNA ψ is one of the few factors previously shown to influence the IRES-dependent translation of cellular mRNAs [12,15]. Studies of Dkc1 hypomorphic mice and the cells of X-DC patients showed that not only do they have distinct ψ profiles but they also exhibit altered translation of specific mRNAs harboring IRES elements, such as Bcl-xL, Xiap, p27, p53, or Vegf [2,13,49,50,51]. Generally, IRES elements require the help of certain canonical translation factors and/or specific cellular IRES-transacting factors (ITAFs) [52]. A heterogeneous ribonucleoprotein D (HNRNPD, also known as AU-rich element RNA-binding protein 1 (AUF1)) is an ITAF that was found to interact with the HCV IRES and thus promote translation from the HCV IRES element [52]. We found the increased association of HNRNPD with the ribosomes of SNORA33-depleted cells, which agrees with our IRES reporter data and increased translation from the HCV IRES. As the 28S-ψ4966 and ribosome IRES-binding site are located relatively far from each other (Figure S5), we suspect that the observed changes in translation initiation are mediated by ITAFs differentially associating with ribosomes. Besides the translational modus, we also evaluated the frequency of the errors occurring during the translation. Knockout of RAS-regulated SNORA24 and the decreased ψ of its rRNA targets 18S-U609 and 18S-U863 led to perturbations in aminoacyl-transfer RNA (aa-tRNA) selection, altered dynamics of the pre-translocation ribosome complex, and increased frequency of translational miscoding and stop codon readthrough [26]. Using a reporter construct based on the HIV-1 PRF sequence, we measured a higher frequency of −1 frameshifting in the SNORA33-depleted cells. Interestingly, ETF1 (also called eRF1), which was shown to decrease the incidence of the frameshift and thus inhibit HIV replication [53,54,55], was more abundantly associated with ribosomes lacking 28S-ψ4966. This may indicate the engagement of some compensatory mechanisms. Altogether, our data suggest that the modification level of 28S-ψ4966 affects ribosome protein stoichiometry, with consequences for translation regulation. In line with this, we also measured specific changes in the cellular proteome of the SNORA33-depleted cells, including decreased expression of several inflammation-regulating factors. One of them was PDCD11, also known as NFBP, which binds NF-κB subunits p50 and p65 and suppresses the expression of inflammatory cytokines [56,57]. A2M is an extracellular macromolecule that functions as a broad-spectrum protease inhibitor. A2M is expressed in articular cartilage and it is also present in synovial fluid [58]. It works as a negative regulator of cartilage catabolic enzymes and therefore inhibits cartilage degeneration [58]. Intraarticular injections of A2M were shown to attenuate the pathogenesis of post-traumatic osteoarthritis in the rat OA model [58,59]. Furthermore, a recent paper demonstrated that A2M binds and neutralizes IL-1β and thus inhibits the downstream NF-κB inflammatory signaling [60]. In this way, A2M inhibits the expression of cartilage-degenerating matrix metalloproteinases (MMPs) and TNFα, while it increases the expression of cartilage protective genes including collagen type II or aggrecan [60]. Inflammation plays a crucial role in osteoarthritis pathobiology [61]. Chondrocytes respond to pro-inflammatory factors in their microenvironment by downregulating their anabolic activities and, in contrast, upregulating the catabolic processes that eventually lead to cartilage degeneration. The downregulation of factors attenuating cellular inflammatory responses in SNORA33-depleted cells suggests that ribosome heterogeneity and the loss of 28S-ψ4966 promote chondrocyte inflammatory responses in osteoarthritis. In the future, more advanced translation measurements, such as Ribo-seq or SILAC mass-spec proteomics, would provide additional details about preferential translation regulation in chondrocytes in health and disease.
A close inspection of the full ψ profile of SNORA33-depleted KO cell pools revealed that besides the ablation of 28S-ψ4966, this also led to a minor increase in ψ levels at three other sites on 28S (28S-U4975, 28S-U3801, 28S-U3741). All four sites were affected by the SNORA33 depletion cluster within two regions of the large ribosomal subunit, domains IV (28S-U3801 and 28S-U3741) and VI (28S-U4975 and 28S-U4966). The concurrent regulation of ψ located close to each other could suggest steric effects. For example, the loss of 28S-ψ4966 might result in more space for the H/ACA snoRNP machinery to catalyze the isomerization of uridine at nearby 28S-U4975. Compensatory or competitive relationships cannot be excluded either. However, at this point, we know very little about the interrelationships and mutual regulation of individual rRNA PTMs to draw any definitive conclusion.
The manipulation of snoRNA levels is a commonly used approach to control the levels of specific rRNA PTMs and to investigate the roles of rRNA-based ribosome heterogeneity in cellular phenotypes and ribosome functions. However, it has certain limitations, mostly due to the potential non-canonical functions of certain snoRNAs. For this reason, the conclusions of these studies need to be interpreted carefully, considering both the direct effects of the loss of PTMs guided by a given snoRNA and its possible indirect (non-canonical) functions. To the best of our knowledge, no non-canonical functions have been reported for SNORA33. Furthermore, we presented data connecting the 28S-ψ4966 modification levels with changes in ribosomal protein composition and thus linking it with the observed ribosome (dys)function and cellular disease phenotype.
Taken together, we present the first data showing that a chronic disease microenvironment instigates site-specific changes in the rRNA ψ profiles of human primary cells. The ablation of a single, disease-sensitive ψ site, 28S-ψ4966, affected ribosome composition and function, with consequences for the cellular proteome, relevant for human disease. Our data broaden the knowledge of ribosome heterogeneity, specialization, and their relevance for human pathologies.

4. Materials and Methods

4.1. Human Articular Chondrocytes Isolation and Culture

Non-OA primary human articular chondrocytes (HACs) were isolated from the cartilage of patients undergoing arthroscopy, anterior cruciate ligament repair, or osteochondritis dissecans surgery according to the previously described procedure [62]. Material collection and use were approved by the METC from the Maastricht University Medical Center (approval number 2017-0183) and informed consent was acquired from all subjects. Non-OA HACs (n = 5, individual donors) were plated at 30.000 cells/cm2 in DMEM/F-12 low glucose with GlutaMAX (Gibco Life Technologies, Waltham, MA, USA, 31331-093), 10% fetal calf serum (FCS; Sigma-Aldrich, Zwijndrecht, The Netherlands, F7524), 1% antibiotic–antimycotic (A-A, Gibco Life Technologies, 15240-062), and 1% nonessential amino acids (NEAA; Gibco Life Technologies, 11140-035). The day after seeding, cells were treated with 20% (v/v) OA-SF (equal volume ratios of OA-SF from 14 OA patients to minimize interpatient variability) or 20% (v/v) 0.9% NaCl. The treatment lasted for 14 days and the medium was refreshed every other day.

4.2. Synovial Fluid Collection

Osteoarthritic synovial fluid was collected from patients undergoing total knee replacement surgery (n = 14, average age 67.1 ± 5.5 years) and stored at −80 °C until further use. Ethical approval was obtained from the Medical Ethics Committee (METC) at the Maastricht University Medical Center (approval number 2017-0183).

4.3. RNA Isolation and HydraPsiSeq

Total RNA was isolated by the RNeasy Plus Mini Kit (Qiagen, Manchester, UK, 74134) and its quality was determined with an RNA 6000 Nano Kit (Agilent, Abcoude, The Netherlands, 5067-1511) and 2100 Bioanalyzer (Agilent). All samples had an RNA integrity (RIN) number greater than 9. HydraPsiSeq of rRNAs was performed as previously described [3]. Total RNA was randomly fragmented at uridine residues by hydrazine (pseudouridines are resistant) and then treated with aniline for RNA scission at abasic sites. Fragments were converted into a library (NEBNext Small RNA Library kit), multiplexed, and sequenced by the Illumina HiSeq 1000 instrument (Illumina, Evry, France). The presence of pseudouridines was detected as the protection of modified U residues against hydrazine cleavage. The depth of the “gap” in the U-protection profile is proportional to the modification level and allows the calculation of a quantitative PsiScore. The value of PsiScore 1.0 represents a fully modified ψ position, while unmodified U shows a PsiScore close to 0 or even a negative value if more intensively cleaved compared to the neighboring U residues.

4.4. Visualization of the Human Ribosome and PTM Mapping

PTM sites were visualized using the PyMOL software 2.5 (Warren L. DeLano, sourced in Maastricht, The Netherlands) and the 3D structure of the human 80S ribosome from the Protein Data Bank (PDB ID: 6QZP) [63,64].

4.5. Lentiviral CRISPR/Cas9 Targeting snoRNAs

LentiCRISPR v2 (Addgene plasmid #52961) was a gift from Dr. Feng Zhang; the pCMVR8.74 (Addgene plasmid #22036) and pMD2.G (Addgene plasmid #12259) plasmids were gifts from Dr. Didier Trono [65]. SgRNA sequences (Table S6) targeting snoRNAs were based on a previous publication [66]. They were annealed and cloned into LentiCRISPRv2 using BsmBI overhangs and sequence-verified. LentiCRISPRv2 with sgRNAs targeting GFP was used as a CRISPR control [67,68]. HEK293T cells were transfected with transfer and production plasmid DNAs (12 μg DNA per 10 cm dish, second-generation lentiviral production system) in equimolar ratios using polyethylenimine (PEI, PEI/DNA ratio 2.5:1; Polysciences, 07923966-2) and lentiviruses were harvested and concentrated using the Lenti-X Concentrator (Takara, San Jose, CA, USA, 631232). Viral titers were determined by p24 ELISA (Fujirebio, Gent, Belgium, 80563) and the multiplicity of infection (MOI) was assessed by serial dilution and transduction. SW1353 cells were transduced at MOI = 1 using 8 μg/mL polybrene (Sigma-Aldrich, Zwijndrecht, The Netherlands, H9268). After 24 h, cells were selected in a medium supplemented with 2 μg/mL puromycin for 3 days (Sigma-Aldrich, P8833). SW1353 (ATCC, HTB-94) were cultured in DMEM/F-12 low glucose with GlutaMAX (Gibco Life Technologies, 31331-093) supplemented with 10% FCS (Sigma-Aldrich, F7524) and 1% A-A (Gibco Life Technologies, 15240-062). HEK293T (ATCC, CRL-3216) were cultured in DMEM high-glucose medium (Gibco Life Technologies, 41966029) supplemented with 10% FCS (Sigma-Aldrich, F7524), 1% A-A (Gibco Life Technologies, 15240-062), and 1% sodium pyruvate (Gibco Life Technologies, 11360070). Cells were cultured at 37 °C in a humidified atmosphere with 5% CO2.

4.6. DNA Isolation and Surveyor Assay

To assess the targeting of the DNA with CRISPR/Cas9, genomic DNA (gDNA) was isolated with the DNeasy Blood & Tissue Kit (Qiagen, Manchester, UK, 69504). High-fidelity Phusion polymerase in HF buffer (Thermo Scientific, Waltham, MA, USA, F530L) and 200 ng gDNA were used for the PCR amplification of the genomic regions spanning the sgRNA target sites (primer sequences are listed in Table S6). The MinElute PCR Purification Kit (Qiagen, 28004) was used to purify and concentrate PCR amplicons. Amplified DNA (300 ng) was hybridized and digested by Surveyor Nuclease S at 42 °C (Surveyor Mutation Detection Assay IDT, 706020); fragments were separated on agarose gel (1.5%) and detected with ethidium bromide and the Chemidoc MP imaging system (Bio-Rad, Hercules, CA, USA).

4.7. RT-qPCR

Total RNA (600 ng) was reversely transcribed using random hexamers (Promega, Madison, WI, USA, C1181). A qPCR was performed using Takyon™ No Rox SYBR Master Mix blue dTTP (Eurogentec, Liege, Belgium, UF-NSMT-B0710), cDNA (6 ng), forward and reverse primers (300 nM), and the following protocol: 50 °C 2 min, denaturation at 95 °C 10 min, 40 cycles of amplification (15 s 95 °C and 1 min 60 °C) (Bio-Rad CFX96 Real-Time PCR Detection System). Data were analyzed using the standard curve method (Bio-Rad CFX Manager Software version 1.1). The relative quantification of target gene expression was normalized to a reference gene and gene expression data were log-transformed for the analysis. Primer sequences are listed in Table S6.

4.8. Proliferation Assay

Cells were seeded subconfluently (10,000 cells/cm2) and cultured in a proliferation medium for 6 days. On the day of sampling (Days 0, 2, 4, and 6) cells were washed with PBS, fixated with 1% glutaraldehyde (VWR, L150739.1000), and stained with 0.1% crystal violet (Sigma-Aldrich, C-3886) in 200 mM boric acid (pH 9.0) (Sigma-Aldrich, B7901). Plates were washed with deionized water, the bound dye was solubilized using 10% acetic acid (VWR, Leicestershire, UK, 20102292), and the optical density was measured at 590 nm using a Multiskan FC Microplate Photometer (Thermo Scientific).

4.9. Total Protein Synthesis

SNORA33 KO and CRISPR control cell pools were cultured for 30 min in methionine and cysteine-free medium supplemented with EasyTag™ EXPRESS 35S Protein Labeling Mix (25 μCi/mL; PerkinElmer Waltham, MA, USA, NEG772002MC) and harvested by scraping in RIPA buffer. A Tri-Carb 2910 TR scintillation counter (PerkinElmer) was used to measure the radioactive signal and data were normalized to the total protein content of the well, as determined by a standard BCA assay (Sigma-Aldrich, 71285-3).

4.10. Polysome Profiling

SNORA33 KO and CRISPR control cell pools at 70–80% confluency were pre-incubated with cycloheximide (CHX, 100 µg/mL, Sigma-Aldrich, C1988) for 20 min at 37 °C, harvested by scraping in cold NaCl (0.9%) supplemented with CHX (100 µg/mL), centrifuged at 1.100 rpm (Hettich Rotanta 460, Tuttlingen, Germany) for 10 min, and gently resuspended in polysome extraction buffer (20 mM Tris–HCl (pH 7.5), 100 mM KCl, 5 mM MgCl2, 0.5% NP-40, CHX 100 μg/mL, protease inhibitor cocktail (Sigma-Aldrich, 11836170001), RNasin (40 U/mL, Promega, N2515)) [69] and incubated on ice for 10 min. To remove nuclei and cellular debris, cytoplasmic extracts were centrifuged at 13.200 rpm (Hettich Micro 200R, Tuttlingen, Germany) at 4 °C for 10 min. Linear sucrose gradients (10–50%) were prepared in open-top Polyclear tubes (Senton, 7030) by layering and mixing an equal amount of 50% and 10% sucrose solutions (0.5 M NaCl, 100 mM Tris–HCl (pH7.5), 50 mM MgCl2) using a Gradient Master 108 (BioComp). Equal amounts of cytoplasmic extracts were loaded onto the sucrose gradients and centrifuged at 39.000 rpm at 4 °C for 1.5 h (SW41 Ti Swinging-Bucket Rotor, Beckman-Coulter Inc, Indianapolis, IN, USA, 331362) [69]. A piston gradient fractionator (BioComp) coupled to a fraction collector (Gilson) and continuous A260 monitoring (Triax Flow Cell, Science Services, Munich, Germany) were used to separate sucrose gradient fractions (24 fractions, 500 μL each).

4.11. Translational Modus and Fidelity Assays

To assess ribosome translation characteristics, IRES-mediated translation initiation and translation fidelity (stop codon skipping and frameshift) were measured using dual-luciferase reporters (DLRs; gifts from Dr. S.R. Thompson and Dr. J. Dinman, respectively; sequences in Table S6). Plasmid DNA was transfected with Fugene6 Transfection Reagent (Promega, E2691; Fugene:DNA 4:1, 0.5 μg DNA/well of 24-well plate). Cells were co-transfected with lacZ plasmid (10% of transfected DNA). After 4 h of incubation, the medium was supplemented with 2% A-A culture medium and incubated for an additional 20 h. Firefly and Renilla luciferase activity were measured 48 h post-transfection by the Dual-Luciferase® Reporter (DLR™) Assay System (Promega, E1910) and a Tristar2 LB942 (Berthold Technologies). The β-galactosidase activity was measured with the β-Gal Assay Kit (Invitrogen, K1455-01) and a Multiskan FC Microplate Photometer (Thermo Scientific). Raw RLU values are listed in Table S7 [70]. The Fluc/Rluc ratio was corrected for the β-galactosidase activity of the well.

4.12. Label-Free Liquid Chromatography–Tandem Mass Spectrometry (LC-MS/MS) Analysis of Cellular and Ribosomal Proteome

For the analysis of the cellular proteomes, SNORA33 KO and CRISPR control cell pools were harvested at 70–80% confluency by scraping in a 25 mM ammonium bicarbonate buffer (Sigma-Aldrich, 09830) supplemented with 7.5U of Benzonase Nuclease (Merck Millipore, Darmstadt, Germany, 70664-3) and a protease inhibitor (Sigma-Aldrich, 11836170001). Samples were sonicated on ice and centrifuged, and supernatants were flash-frozen in liquid nitrogen. For the analyses of the ribosomal proteome, SNORA33 KO and CRISPR control cell pools were pre-incubated with CHX (100 µg/mL, Sigma-Aldrich, C1988) for 5 min at 37 °C and harvested by scraping in cold NaCl (0.9%) supplemented with CHX (100 µg/mL). Samples were spun at 1.100 rpm (Hettich Rotanta 460) for 5 min and gently resuspended in cytoplasmic extraction buffer (20 mM Tris pH 7.5, 5 mM MgCl2, 10 mM NaCl, 0.15% NP-40), freshly supplemented with CHX (100 µg/mL), a protease inhibitor cocktail (Sigma-Aldrich, 11836170001), and recombinant RNasin (Promega, N2515). Cells were lysed using a Dounce homogenizer (50 strokes/sample). To remove nuclei and mitochondria, cytoplasmic extracts were centrifuged first at 3.900 rpm (Hettich Micro 200R) at 4 °C for 10 min and afterwards at 14.000 rpm (Hettich Micro 200R) at 4 °C for 10 min. Equal amounts of cytoplasmic extracts were loaded on low-salt sucrose cushions for the isolation of ribosomes and associated proteins (1 M and 0.7 M sucrose cushion; 20 mM Tris pH 7.5, 5 mM MgCl2, 25 mM NaCl) [15] and ultracentrifuged at 32.200 rpm at 4 °C for 17 h (SW 41 Ti rotor, Beckman-Coulter, 331362). The ribosomal pellet was resuspended in 25 mM ammonium bicarbonate buffer (Sigma-Aldrich, 09830) supplemented with 7.5U of Benzonase Nuclease (Merck Millipore, 70664-3) and a protease inhibitor (Sigma-Aldrich, 11836170001) and flash-frozen in liquid nitrogen. In-solution tryptic digestion of the cell lysates was performed as previously described [71]. In short, 25 mM ammonium bicarbonate (Fluka Chemicals Ltd., Gillingham, UK) containing 0.05% (w/v) RapiGest (Waters, Elstree, UK) was added to 100 μg protein of each sample to produce a final volume of 160 µL and heated at 80 °C for 10 min. DL dithiothreitol (Sigma-Aldrich) was added (3 mM final concentration) and incubated at 60 °C for 10 min. Iodoacetamide (Sigma-Aldrich) was added (9 mM final concentration) and incubated at room temperature for 30 min in the dark. Tryptic digestion was undertaken using Lys-C endopeptidase (FUJIFILM Wako Pure Chemical, Osaka, Japan) at a ratio of 1:25 177 (enzyme:protein) and incubated at 37 °C for 16 h, followed by a second trypsin supplementation for 2 h. Digestion was halted by the addition of trifluoroacetic acid (Sigma-Aldrich) to 0.5% (v/v final concentration) and rotation for 30 min at 37 °C. Samples were centrifuged at 4 °C at 13,000 rpm for 10 min and the supernatant was removed. Five hundred ng of each tryptic digest was analyzed using LC-MS/MS, on a 2 h gradient for the cellular proteome and a 30 min gradient for the ribosomal proteome. Data analyses were performed as previously described [71]. The Progenesis QI software (V4, Waters) was used for protein quantification (only unique peptides were considered) [72]. Data are available via ProteomeXchange with the identifier PXD042545. We used Reactome to categorize all identified proteins in the ribosomal proteome [41].

4.13. Statistical Testing

All statistical tests were performed using GraphPad Prism 5.0.1. (San Diego, CA, USA). A normal distribution was tested or assumed (indicated in figure legends). Differences between the two groups were determined by a two-tailed unpaired t-test or Mann–Whitney U test. Two paired groups were analyzed by a two-tailed paired t-test. The LC-MS/MS cellular and ribosomal proteome data were log-transformed and analyzed by one-way analysis of variance (ANOVA). Additional details regarding the statistical analyses are indicated in the figure legends. Data are presented as mean ± SD. Ns—not significant, * p < 0.05, ** p < 0.01, *** p < 0.001.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms241612578/s1.

Author Contributions

Conceptualization, A.C. (Alzbeta Chabronova), G.v.d.A. and T.J.M.W.; Resources, B.A.C.H., A.C. (Andy Cremers), D.A.M.S. and T.J.M.W.; Methodology, A.C. (Alzbeta Chabronova), G.v.d.A., T.J.M.W., V.M. and Y.M.; Investigation, A.C. (Alzbeta Chabronova), G.v.d.A., A.C. (Andy Cremers), D.A.M.S., M.J.P., V.M., Y.M. and T.J.M.W.; Formal analysis, A.C. (Alzbeta Chabronova), G.v.d.A., M.J.P., V.M. and Y.M.; Visualization, A.C. (Alzbeta Chabronova) and G.v.d.A.; Writing—Original Draft, A.C. (Alzbeta Chabronova), G.v.d.A. and T.J.M.W.; Writing—Review and Editing, A.C. (Alzbeta Chabronova), GG.v.d.A., B.A.C.H., M.M.J.C., A.C. (Andy Cremers), D.A.M.S., M.J.P., L.W.v.R., V.M., Y.M. and T.J.M.W.; Funding Acquisition, L.W.v.R. and T.J.M.W.; Supervision, G.v.d.A. and T.J.M.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by a grant from Stichting de Weijerhorst (Bewegen zonder Pijn) and grants from the ReumaNederland (grants 17-2-401 and LLP14). Mandy Peffers is funded through a Wellcome Trust Clinical Intermediate Fellowship (grant 107471/Z/15/Z) and supported by Versus Arthritis as part of the MRC Versus Arthritis Centre for Integrated Research into Musculoskeletal Ageing (CIMA).

Institutional Review Board Statement

Material collection and use were approved by the METC from the Maastricht University Medical Center (approval number 2017-0183).

Informed Consent Statement

Informed consent was acquired from all subjects involved in this study.

Data Availability Statement

Proteomic data are available via ProteomeXchange with the identifier PXD042545.

Acknowledgments

The visualization of the human ribosome and PTM mapping was kindly performed by Kanin Wichapong from the Cardiovascular Research Institute Maastricht (CARIM), Department of Biochemistry, Maastricht University, The Netherlands. Panel A in Figure 1 and Figure 5 were created with BioRender.com (accessed on 14 October 2022).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

2′-O-me2′-O-methylation
A-Aantibiotic–antimycotic
ANOVAanalysis of variance
ANSanisomycin
AP1G1adaptor related protein complex 1 subunit gamma 1
BCAbicinchoninic acid
CHXcycloheximide
CPMcounts per minute
DEdifferentially expressed
DLRdual-luciferase
DMEM/F-12Dulbecco’s modified eagle medium/nutrient mixture F-12
dTTPdeoxythymidine triphosphate
ELISAenzyme-linked immunoassay
EMTemetine
FAfucidic acid
FBLfibrillarin
FCfold change
FCSfetal calf serum
gDNAgenomic DNA
GFPgreen fluorescent protein
sgRNAsingle guide RNA
HAChuman articular chondrocyte
KClpotassium chloride
KOknockout
Mmale
MgCl2magnesium chloride
NaClsodium chloride
NPM1nucleophosmin
nmnanometer
OAosteoarthritis
OA-SFosteoarthritic synovial fluid
PBSphosphate-buffered saline
RLU relative light unit
RINRNA integrity number
RIPA radioimmunoprecipitation assay
rRNAribosomal RNA
RT-qPCRquantitative reverse transcription polymerase chain reaction
snoRNAsmall nucleolar RNA

References

  1. Khatter, H.; Myasnikov, A.G.; Natchiar, S.K.; Klaholz, B.P. Structure of the human 80S ribosome. Nature 2015, 520, 640–645. [Google Scholar] [CrossRef]
  2. Taoka, M.; Nobe, Y.; Yamaki, Y.; Sato, K.; Ishikawa, H.; Izumikawa, K.; Yamauchi, Y.; Hirota, K.; Nakayama, H.; Takahashi, N.; et al. Landscape of the complete RNA chemical modifications in the human 80S ribosome. Nucleic Acids Res. 2018, 46, 9289–9298. [Google Scholar] [CrossRef]
  3. Marchand, V.; Pichot, F.; Neybecker, P.; Ayadi, L.; Bourguignon-Igel, V.; Wacheul, L.; Lafontaine, D.L.J.; Pinzano, A.; Helm, M.; Motorin, Y. HydraPsiSeq: A method for systematic and quantitative mapping of pseudouridines in RNA. Nucleic Acids Res. 2020, 48, e110. [Google Scholar] [CrossRef]
  4. Motorin, Y.; Quinternet, M.; Rhalloussi, W.; Marchand, V. Constitutive and variable 2′-O-methylation (Nm) in human ribosomal RNA. RNA Biol. 2021, 18 (Suppl. 1), 88–97. [Google Scholar] [CrossRef] [PubMed]
  5. Ojha, S.; Malla, S.; Lyons, S.M. snoRNPs: Functions in Ribosome Biogenesis. Biomolecules 2020, 10, 783. [Google Scholar] [CrossRef] [PubMed]
  6. Kiss, T. Small nucleolar RNAs: An abundant group of noncoding RNAs with diverse cellular functions. Cell 2002, 109, 145–148. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Ganot, P.; Bortolin, M.L.; Kiss, T. Site-specific pseudouridine formation in preribosomal RNA is guided by small nucleolar RNAs. Cell 1997, 89, 799–809. [Google Scholar] [CrossRef] [Green Version]
  8. Decatur, W.A.; Fournier, M.J. rRNA modifications and ribosome function. Trends Biochem. Sci. 2002, 27, 344–351. [Google Scholar] [CrossRef]
  9. Gay, D.M.; Lund, A.H.; Jansson, M.D. Translational control through ribosome heterogeneity and functional specialization. Trends Biochem. Sci. 2022, 47, 66–81. [Google Scholar] [CrossRef]
  10. Genuth, N.R.; Barna, M. The Discovery of Ribosome Heterogeneity and Its Implications for Gene Regulation and Organismal Life. Mol. Cell 2018, 71, 364–374. [Google Scholar] [CrossRef] [Green Version]
  11. Jaafar, M.; Paraqindes, H.; Gabut, M.; Diaz, J.J.; Marcel, V.; Durand, S. 2’O-Ribose Methylation of Ribosomal RNAs: Natural Diversity in Living Organisms, Biological Processes, and Diseases. Cells 2021, 10, 1948. [Google Scholar] [CrossRef] [PubMed]
  12. Jack, K.; Bellodi, C.; Landry, D.M.; Niederer, R.O.; Meskauskas, A.; Musalgaonkar, S.; Kopmar, N.; Krasnykh, O.; Dean, A.M.; Thompson, S.R.; et al. rRNA pseudouridylation defects affect ribosomal ligand binding and translational fidelity from yeast to human cells. Mol. Cell 2011, 44, 660–666. [Google Scholar] [CrossRef] [Green Version]
  13. Penzo, M.; Montanaro, L. Turning Uridines around: Role of rRNA Pseudouridylation in Ribosome Biogenesis and Ribosomal Function. Biomolecules 2018, 8, 38. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Rocchi, L.; Pacilli, A.; Sethi, R.; Penzo, M.; Schneider, R.J.; Trere, D.; Brigotti, M.; Montanaro, L. Dyskerin depletion increases VEGF mRNA internal ribosome entry site-mediated translation. Nucleic Acids Res. 2013, 41, 8308–8318. [Google Scholar] [CrossRef] [Green Version]
  15. Penzo, M.; Rocchi, L.; Brugiere, S.; Carnicelli, D.; Onofrillo, C.; Coute, Y.; Brigotti, M.; Montanaro, L. Human ribosomes from cells with reduced dyskerin levels are intrinsically altered in translation. FASEB J. 2015, 29, 3472–3482. [Google Scholar] [CrossRef] [Green Version]
  16. Jansson, M.D.; Hafner, S.J.; Altinel, K.; Tehler, D.; Krogh, N.; Jakobsen, E.; Andersen, J.V.; Andersen, K.L.; Schoof, E.M.; Menard, P.; et al. Regulation of translation by site-specific ribosomal RNA methylation. Nat. Struct. Mol. Biol. 2021, 28, 889–899. [Google Scholar] [CrossRef] [PubMed]
  17. Heiss, N.S.; Knight, S.W.; Vulliamy, T.J.; Klauck, S.M.; Wiemann, S.; Mason, P.J.; Poustka, A.; Dokal, I. X-linked dyskeratosis congenita is caused by mutations in a highly conserved gene with putative nucleolar functions. Nat. Genet. 1998, 19, 32–38. [Google Scholar] [CrossRef]
  18. Knight, S.W.; Heiss, N.S.; Vulliamy, T.J.; Greschner, S.; Stavrides, G.; Pai, G.S.; Lestringant, G.; Varma, N.; Mason, P.J.; Dokal, I.; et al. X-linked dyskeratosis congenita is predominantly caused by missense mutations in the DKC1 gene. Am. J. Hum. Genet. 1999, 65, 50–58. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Savage, S.A.; Alter, B.P. Dyskeratosis congenita. Hematol. Oncol. Clin. N. Am. 2009, 23, 215–231. [Google Scholar] [CrossRef]
  20. Ejby, M.; Sorensen, M.A.; Pedersen, S. Pseudouridylation of helix 69 of 23S rRNA is necessary for an effective translation termination. Proc. Natl. Acad. Sci. USA 2007, 104, 19410–19415. [Google Scholar] [CrossRef]
  21. Gutgsell, N.S.; Deutscher, M.P.; Ofengand, J. The pseudouridine synthase RluD is required for normal ribosome assembly and function in Escherichia coli. RNA 2005, 11, 1141–1152. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Baudin-Baillieu, A.; Fabret, C.; Liang, X.H.; Piekna-Przybylska, D.; Fournier, M.J.; Rousset, J.P. Nucleotide modifications in three functionally important regions of the Saccharomyces cerevisiae ribosome affect translation accuracy. Nucleic Acids Res. 2009, 37, 7665–7677. [Google Scholar] [CrossRef] [Green Version]
  23. Liang, X.H.; Liu, Q.; Fournier, M.J. rRNA modifications in an intersubunit bridge of the ribosome strongly affect both ribosome biogenesis and activity. Mol. Cell 2007, 28, 965–977. [Google Scholar] [CrossRef]
  24. Soung, Y.H.; Lee, J.W.; Kim, S.Y.; Nam, S.W.; Park, W.S.; Lee, J.Y.; Yoo, N.J.; Lee, S.H. Absence of DKC1 exon 3 mutation in common human cancers. Acta Oncol. 2006, 45, 342–343. [Google Scholar] [CrossRef] [PubMed]
  25. Liang, J.; Wen, J.; Huang, Z.; Chen, X.P.; Zhang, B.X.; Chu, L. Small Nucleolar RNAs: Insight Into Their Function in Cancer. Front. Oncol. 2019, 9, 587. [Google Scholar] [CrossRef] [Green Version]
  26. McMahon, M.; Contreras, A.; Holm, M.; Uechi, T.; Forester, C.M.; Pang, X.; Jackson, C.; Calvert, M.E.; Chen, B.; Quigley, D.A.; et al. A single H/ACA small nucleolar RNA mediates tumor suppression downstream of oncogenic RAS. Elife 2019, 8, e48847. [Google Scholar] [CrossRef]
  27. Cui, A.; Li, H.; Wang, D.; Zhong, J.; Chen, Y.; Lu, H. Global, regional prevalence, incidence and risk factors of knee osteoarthritis in population-based studies. EClinicalMedicine 2020, 29–30, 100587. [Google Scholar] [CrossRef] [PubMed]
  28. Loeser, R.F.; Goldring, S.R.; Scanzello, C.R.; Goldring, M.B. Osteoarthritis: A disease of the joint as an organ. Arthritis Rheum. 2012, 64, 1697–1707. [Google Scholar] [CrossRef] [Green Version]
  29. Thijssen, E.; van Caam, A.; van der Kraan, P.M. Obesity and osteoarthritis, more than just wear and tear: Pivotal roles for inflamed adipose tissue and dyslipidaemia in obesity-induced osteoarthritis. Rheumatology 2015, 54, 588–600. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Balakrishnan, L.; Nirujogi, R.S.; Ahmad, S.; Bhattacharjee, M.; Manda, S.S.; Renuse, S.; Kelkar, D.S.; Subbannayya, Y.; Raju, R.; Goel, R.; et al. Proteomic analysis of human osteoarthritis synovial fluid. Clin. Proteom. 2014, 11, 6. [Google Scholar] [CrossRef] [Green Version]
  31. Carlson, A.K.; Rawle, R.A.; Adams, E.; Greenwood, M.C.; Bothner, B.; June, R.K. Application of global metabolomic profiling of synovial fluid for osteoarthritis biomarkers. Biochem. Biophys. Res. Commun. 2018, 499, 182–188. [Google Scholar] [CrossRef]
  32. Kosinska, M.K.; Liebisch, G.; Lochnit, G.; Wilhelm, J.; Klein, H.; Kaesser, U.; Lasczkowski, G.; Rickert, M.; Schmitz, G.; Steinmeyer, J. A lipidomic study of phospholipid classes and species in human synovial fluid. Arthritis Rheum. 2013, 65, 2323–2333. [Google Scholar] [CrossRef] [PubMed]
  33. Ingelmark, B.E. The nutritive supply and nutritional value of synovial fluid. Acta Orthop. Scand. 1950, 20, 144–155. [Google Scholar] [CrossRef] [PubMed]
  34. Sandell, L.J.; Aigner, T. Articular cartilage and changes in arthritis. An introduction: Cell biology of osteoarthritis. Arthritis Res. 2001, 3, 107–113. [Google Scholar] [CrossRef] [Green Version]
  35. van den Akker, G.G.H.; Caron, M.M.J.; Peffers, M.J.; Welting, T.J.M. Ribosome dysfunction in osteoarthritis. Curr. Opin. Rheumatol. 2022, 34, 61–67. [Google Scholar] [CrossRef]
  36. Chabronova, A.; van den Akker, G.G.H.; Housmans, B.A.C.; Caron, M.M.J.; Cremers, A.; Surtel, D.A.M.; Wichapong, K.; Peffers, M.M.J.; van Rhijn, L.W.; Marchand, V.; et al. Ribosomal RNA-based epitranscriptomic regulation of chondrocyte translation and proteome in osteoarthritis. Osteoarthr. Cartil. 2023, 31, 374–385. [Google Scholar] [CrossRef] [PubMed]
  37. Housmans, B.A.C.; Neefjes, M.; Surtel, D.A.M.; Vitik, M.; Cremers, A.; van Rhijn, L.W.; van der Kraan, P.M.; van den Akker, G.G.H.; Welting, T.J.M. Synovial fluid from end-stage osteoarthritis induces proliferation and fibrosis of articular chondrocytes via MAPK and RhoGTPase signaling. Osteoarthr. Cartil. 2022, 30, 862–874. [Google Scholar] [CrossRef]
  38. Tonna, E.A.; Singh, I.J. An autoradiographic investigation of 3H-uridine utilization by aging mouse cartilage cells. Exp. Gerontol. 1976, 11, 231–234. [Google Scholar] [CrossRef]
  39. Kiss, A.M.; Jady, B.E.; Bertrand, E.; Kiss, T. Human box H/ACA pseudouridylation guide RNA machinery. Mol. Cell. Biol. 2004, 24, 5797–5807. [Google Scholar] [CrossRef] [Green Version]
  40. Belin, S.; Hacot, S.; Daudignon, L.; Therizols, G.; Pourpe, S.; Mertani, H.C.; Rosa-Calatrava, M.; Diaz, J.J. Purification of ribosomes from human cell lines. Curr. Protoc. Cell Biol. 2010, 49, 3–40. [Google Scholar] [CrossRef]
  41. Gillespie, M.; Jassal, B.; Stephan, R.; Milacic, M.; Rothfels, K.; Senff-Ribeiro, A.; Griss, J.; Sevilla, C.; Matthews, L.; Gong, C.; et al. The reactome pathway knowledgebase 2022. Nucleic Acids Res. 2022, 50, D687–D692. [Google Scholar] [CrossRef] [PubMed]
  42. Seither, P.; Coy, J.F.; Pouska, A.; Grummt, I. Molecular cloning and characterization of the cDNA encoding the largest subunit of mouse RNA polymerase I. Mol. Gen. Genet. 1997, 255, 180–186. [Google Scholar] [CrossRef] [PubMed]
  43. Housmans, B.A.C.; van den Akker, G.G.H.; Neefjes, M.; Timur, U.T.; Cremers, A.; Peffers, M.J.; Caron, M.M.J.; van Rhijn, L.W.; Emans, P.J.; Boymans, T.A.E.J.; et al. Direct comparison of non-osteoarthritic and osteoarthritic synovial fluid-induced intracellular chondrocyte signaling and phenotype changes. Osteoarthr. Cartil. 2023, 31, 60–71. [Google Scholar] [CrossRef] [PubMed]
  44. Pisarev, A.V.; Kolupaeva, V.G.; Yusupov, M.M.; Hellen, C.U.; Pestova, T.V. Ribosomal position and contacts of mRNA in eukaryotic translation initiation complexes. EMBO J. 2008, 27, 1609–1621. [Google Scholar] [CrossRef] [Green Version]
  45. Jiao, J.; Kavdia, K.; Pagala, V.; Palmer, L.; Finkelstein, D.; Fan, Y.; Peng, J.; Demontis, F. An age-downregulated ribosomal RpS28 protein variant regulates the muscle proteome. G3 (Bethesda) 2021, 11, jkab165. [Google Scholar] [CrossRef]
  46. Gripp, K.W.; Curry, C.; Olney, A.H.; Sandoval, C.; Fisher, J.; Chong, J.X.; Genomics, U.W.C.f.M.; Pilchman, L.; Sahraoui, R.; Stabley, D.L.; et al. Diamond-Blackfan anemia with mandibulofacial dystostosis is heterogeneous, including the novel DBA genes TSR2 and RPS28. Am. J. Med. Genet. A 2014, 164A, 2240–2249. [Google Scholar] [CrossRef] [Green Version]
  47. Pelletier, J.; Sonenberg, N. Internal initiation of translation of eukaryotic mRNA directed by a sequence derived from poliovirus RNA. Nature 1988, 334, 320–325. [Google Scholar] [CrossRef]
  48. Godet, A.C.; David, F.; Hantelys, F.; Tatin, F.; Lacazette, E.; Garmy-Susini, B.; Prats, A.C. IRES Trans-Acting Factors, Key Actors of the Stress Response. Int. J. Mol. Sci. 2019, 20, 924. [Google Scholar] [CrossRef] [Green Version]
  49. Bellodi, C.; Krasnykh, O.; Haynes, N.; Theodoropoulou, M.; Peng, G.; Montanaro, L.; Ruggero, D. Loss of function of the tumor suppressor DKC1 perturbs p27 translation control and contributes to pituitary tumorigenesis. Cancer Res. 2010, 70, 6026–6035. [Google Scholar] [CrossRef] [Green Version]
  50. Ruggero, D.; Grisendi, S.; Piazza, F.; Rego, E.; Mari, F.; Rao, P.H.; Cordon-Cardo, C.; Pandolfi, P.P. Dyskeratosis congenita and cancer in mice deficient in ribosomal RNA modification. Science 2003, 299, 259–262. [Google Scholar] [CrossRef]
  51. Yoon, A.; Peng, G.; Brandenburger, Y.; Zollo, O.; Xu, W.; Rego, E.; Ruggero, D. Impaired control of IRES-mediated translation in X-linked dyskeratosis congenita. Science 2006, 312, 902–906. [Google Scholar] [CrossRef]
  52. Paek, K.Y.; Kim, C.S.; Park, S.M.; Kim, J.H.; Jang, S.K. RNA-binding protein hnRNP D modulates internal ribosome entry site-dependent translation of hepatitis C virus RNA. J. Virol. 2008, 82, 12082–12093. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Kobayashi, Y.; Zhuang, J.; Peltz, S.; Dougherty, J. Identification of a cellular factor that modulates HIV-1 programmed ribosomal frameshifting. J. Biol. Chem. 2010, 285, 19776–19784. [Google Scholar] [CrossRef] [Green Version]
  54. Mathew, S.F.; Crowe-McAuliffe, C.; Graves, R.; Cardno, T.S.; McKinney, C.; Poole, E.S.; Tate, W.P. The highly conserved codon following the slippery sequence supports -1 frameshift efficiency at the HIV-1 frameshift site. PLoS ONE 2015, 10, e0122176. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Brass, A.L.; Dykxhoorn, D.M.; Benita, Y.; Yan, N.; Engelman, A.; Xavier, R.J.; Lieberman, J.; Elledge, S.J. Identification of host proteins required for HIV infection through a functional genomic screen. Science 2008, 319, 921–926. [Google Scholar] [CrossRef]
  56. Sweet, T.; Khalili, K.; Sawaya, B.E.; Amini, S. Identification of a novel protein from glial cells based on its ability to interact with NF-kappaB subunits. J. Cell. Biochem. 2003, 90, 884–891. [Google Scholar] [CrossRef]
  57. Yang, R.; Zhan, M.; Guo, M.; Yuan, H.; Wang, Y.; Zhang, Y.; Zhang, W.; Chen, S.; de The, H.; Chen, Z.; et al. Yolk sac-derived Pdcd11-positive cells modulate zebrafish microglia differentiation through the NF-kappaB-Tgfbeta1 pathway. Cell Death Differ. 2021, 28, 170–183. [Google Scholar] [CrossRef] [PubMed]
  58. Wang, S.; Wei, X.; Zhou, J.; Zhang, J.; Li, K.; Chen, Q.; Terek, R.; Fleming, B.C.; Goldring, M.B.; Ehrlich, M.G.; et al. Identification of alpha2-macroglobulin as a master inhibitor of cartilage-degrading factors that attenuates the progression of posttraumatic osteoarthritis. Arthritis Rheumatol. 2014, 66, 1843–1853. [Google Scholar] [CrossRef] [Green Version]
  59. Zhang, Y.; Wei, X.; Browning, S.; Scuderi, G.; Hanna, L.S.; Wei, L. Targeted designed variants of alpha-2-macroglobulin (A2M) attenuate cartilage degeneration in a rat model of osteoarthritis induced by anterior cruciate ligament transection. Arthritis Res. Ther. 2017, 19, 175. [Google Scholar] [CrossRef] [Green Version]
  60. Sun, C.; Cao, C.; Zhao, T.; Guo, H.; Fleming, B.C.; Owens, B.; Beveridge, J.; McAllister, S.; Wei, L. A2M inhibits inflammatory mediators of chondrocytes by blocking IL-1beta/NF-kappaB pathway. J. Orthop. Res. 2022, 41, 241–248. [Google Scholar] [CrossRef]
  61. Goldring, M.B.; Otero, M. Inflammation in osteoarthritis. Curr. Opin. Rheumatol. 2011, 23, 471–478. [Google Scholar] [CrossRef]
  62. Caron, M.M.; Emans, P.J.; Coolsen, M.M.; Voss, L.; Surtel, D.A.; Cremers, A.; van Rhijn, L.W.; Welting, T.J. Redifferentiation of dedifferentiated human articular chondrocytes: Comparison of 2D and 3D cultures. Osteoarthr. Cartil. 2012, 20, 1170–1178. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Berman, H.M.; Westbrook, J.; Feng, Z.; Gilliland, G.; Bhat, T.N.; Weissig, H.; Shindyalov, I.N.; Bourne, P.E. The Protein Data Bank. Nucleic Acids Res. 2000, 28, 235–242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Natchiar, S.K.; Myasnikov, A.G.; Kratzat, H.; Hazemann, I.; Klaholz, B.P. Visualization of chemical modifications in the human 80S ribosome structure. Nature 2017, 551, 472–477. [Google Scholar] [CrossRef] [PubMed]
  65. Sanjana, N.E.; Shalem, O.; Zhang, F. Improved vectors and genome-wide libraries for CRISPR screening. Nat. Methods 2014, 11, 783–784. [Google Scholar] [CrossRef] [Green Version]
  66. Pauli, C.; Liu, Y.; Rohde, C.; Cui, C.; Fijalkowska, D.; Gerloff, D.; Walter, C.; Krijgsveld, J.; Dugas, M.; Edemir, B.; et al. Site-specific methylation of 18S ribosomal RNA by SNORD42A is required for acute myeloid leukemia cell proliferation. Blood 2020, 135, 2059–2070. [Google Scholar] [CrossRef]
  67. Di Ceglie, I.; van den Akker, G.G.; Ascone, G.; Ten Harkel, B.; Hacker, H.; van de Loo, F.A.; Koenders, M.I.; van der Kraan, P.M.; de Vries, T.J.; Vogl, T.; et al. Genetic modification of ER-Hoxb8 osteoclast precursors using CRISPR/Cas9 as a novel way to allow studies on osteoclast biology. J. Leukoc. Biol. 2017, 101, 957–966. [Google Scholar] [CrossRef]
  68. Shalem, O.; Sanjana, N.E.; Hartenian, E.; Shi, X.; Scott, D.A.; Mikkelson, T.; Heckl, D.; Ebert, B.L.; Root, D.E.; Doench, J.G.; et al. Genome-scale CRISPR-Cas9 knockout screening in human cells. Science 2014, 343, 84–87. [Google Scholar] [CrossRef] [Green Version]
  69. Panda, A.C.; Martindale, J.L.; Gorospe, M. Polysome Fractionation to Analyze mRNA Distribution Profiles. Bio. Protoc. 2017, 7, e2126. [Google Scholar] [CrossRef] [Green Version]
  70. van den Akker, G.G.H.; Zacchini, F.; Housmans, B.A.C.; van der Vloet, L.; Caron, M.M.J.; Montanaro, L.; Welting, T.J.M. Current Practice in Bicistronic IRES Reporter Use: A Systematic Review. Int. J. Mol. Sci. 2021, 22, 5193. [Google Scholar] [CrossRef]
  71. Ripmeester, E.G.J.; Caron, M.M.J.; van den Akker, G.G.H.; Surtel, D.A.M.; Cremers, A.; Balaskas, P.; Dyer, P.; Housmans, B.A.C.; Chabronova, A.; Smagul, A.; et al. Impaired chondrocyte U3 snoRNA expression in osteoarthritis impacts the chondrocyte protein translation apparatus. Sci. Rep. 2020, 10, 13426. [Google Scholar] [CrossRef] [PubMed]
  72. Peffers, M.J.; Thorpe, C.T.; Collins, J.A.; Eong, R.; Wei, T.K.; Screen, H.R.; Clegg, P.D. Proteomic analysis reveals age-related changes in tendon matrix composition, with age- and injury-specific matrix fragmentation. J. Biol. Chem. 2014, 289, 25867–25878. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. The chronic disease microenvironment induces site-specific changes in the pseudouridylation levels of ribosomal RNAs in primary human chondrocytes. (A) A schematic of the experimental design. Non-OA human articular chondrocytes of 5 individual donors were exposed to a chronic disease microenvironment represented by the synovial fluid of end-stage OA patients (OA-SF, pool of 14 donors, 20% (v/v)) for 14 days. The culture medium was refreshed every other day. After 14 days of culture, total RNA was isolated and used for the pseudouridylation profiling of rRNAs by HydraPsiSeq. (B) Differentially pseudouridylated (ψ) rRNA nucleotides (n = 5). Statistical significance was assessed by paired t-test with the assumption of a normal distribution of the data. Full ψ rRNA profiles are shown in Figure S1 and PsiScore values are listed in Table S1. (C) Location of 28S-U4966 within domain VI of the large ribosomal subunit of the human ribosome. * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 1. The chronic disease microenvironment induces site-specific changes in the pseudouridylation levels of ribosomal RNAs in primary human chondrocytes. (A) A schematic of the experimental design. Non-OA human articular chondrocytes of 5 individual donors were exposed to a chronic disease microenvironment represented by the synovial fluid of end-stage OA patients (OA-SF, pool of 14 donors, 20% (v/v)) for 14 days. The culture medium was refreshed every other day. After 14 days of culture, total RNA was isolated and used for the pseudouridylation profiling of rRNAs by HydraPsiSeq. (B) Differentially pseudouridylated (ψ) rRNA nucleotides (n = 5). Statistical significance was assessed by paired t-test with the assumption of a normal distribution of the data. Full ψ rRNA profiles are shown in Figure S1 and PsiScore values are listed in Table S1. (C) Location of 28S-U4966 within domain VI of the large ribosomal subunit of the human ribosome. * p < 0.05, ** p < 0.01, *** p < 0.001.
Ijms 24 12578 g001
Figure 2. CRISPR/Cas9-mediated SNORA33 depletion leads to the loss of pseudouridylation at 28S-U4966. (A,E) Two box H/ACA snoRNAs, SNORA22 and SNORA33, are predicted to guide the ψ of 28S-U4966. (B,F) A scheme of the double sgRNA CRISPR/Cas9-mediated approach to generate SNORA22 and SNORA33 SW1353 KO cell pools. Top: Exons of host genes CCT6P1 and RPS12 flanking the intron-encoded SNORA22 and SNORA33, respectively. Bottom: Scheme of SNORA22 and SNORA33 sequences with indicated sgRNA1 and sgRNA2 cleavage sites and rRNA complementary sequences. The cell pool of GFP targeting CRISPR/Cas9-treated cells was used as a control (CRISPR Ctrl). (C) Expression levels of SNORA22, expression and splicing (exons 2 and 3) of its host gene CCT6P1 measured by RT-qPCR (n = 3). (G) Expression levels of SNORA33, expression and splicing (exons 5 and 6) of its host gene RPS12 measured by RT-qPCR (n = 3). Expression levels were normalized to the reference gene (PPIA) expression and plotted as fold change to CRISPR Ctrl. Data were analyzed by unpaired t-test with the assumption of a normal distribution. (D,H) Pseudouridylation levels of 28S-U4966 (n = 3) measured by HydraPsiSeq and analyzed by unpaired t-test with the assumption of a normal distribution of data. Negative PsiScore for 28S-ψ4966 in the case of SNORA33 KO cell pool indicates even more intensive hydrazine cleavage at this position compared to other neighboring U residues. Ns – not significant, * p < 0.05, *** p < 0.001.
Figure 2. CRISPR/Cas9-mediated SNORA33 depletion leads to the loss of pseudouridylation at 28S-U4966. (A,E) Two box H/ACA snoRNAs, SNORA22 and SNORA33, are predicted to guide the ψ of 28S-U4966. (B,F) A scheme of the double sgRNA CRISPR/Cas9-mediated approach to generate SNORA22 and SNORA33 SW1353 KO cell pools. Top: Exons of host genes CCT6P1 and RPS12 flanking the intron-encoded SNORA22 and SNORA33, respectively. Bottom: Scheme of SNORA22 and SNORA33 sequences with indicated sgRNA1 and sgRNA2 cleavage sites and rRNA complementary sequences. The cell pool of GFP targeting CRISPR/Cas9-treated cells was used as a control (CRISPR Ctrl). (C) Expression levels of SNORA22, expression and splicing (exons 2 and 3) of its host gene CCT6P1 measured by RT-qPCR (n = 3). (G) Expression levels of SNORA33, expression and splicing (exons 5 and 6) of its host gene RPS12 measured by RT-qPCR (n = 3). Expression levels were normalized to the reference gene (PPIA) expression and plotted as fold change to CRISPR Ctrl. Data were analyzed by unpaired t-test with the assumption of a normal distribution. (D,H) Pseudouridylation levels of 28S-U4966 (n = 3) measured by HydraPsiSeq and analyzed by unpaired t-test with the assumption of a normal distribution of data. Negative PsiScore for 28S-ψ4966 in the case of SNORA33 KO cell pool indicates even more intensive hydrazine cleavage at this position compared to other neighboring U residues. Ns – not significant, * p < 0.05, *** p < 0.001.
Ijms 24 12578 g002
Figure 3. The effect of SNORA33 depletion on rRNA ψ profile. (A) Differentially ψ rRNA nucleotides in SNORA33-depleted cells (n = 3). Statistical significance was assessed by unpaired t-test with the assumption of a normal distribution of data. Full ψ rRNA profiles are shown in Figure S4 and PsiScore values are listed in Table S2. (B) Negative PsiScore for 28S-ψ4966 in the case of SNORA33 KO cell pool indicates even more intensive hydrazine cleavage at this position compared to other neighboring U residues. The locations of differentially ψ sites identified in SNORA33-depleted cell pools mapped within a human ribosome. 28S-U3801 and 28S-U3741 are located in domain IV, and 28S-U4975 and 28S-U4966 are situated within domain VI of the large ribosomal subunit. * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 3. The effect of SNORA33 depletion on rRNA ψ profile. (A) Differentially ψ rRNA nucleotides in SNORA33-depleted cells (n = 3). Statistical significance was assessed by unpaired t-test with the assumption of a normal distribution of data. Full ψ rRNA profiles are shown in Figure S4 and PsiScore values are listed in Table S2. (B) Negative PsiScore for 28S-ψ4966 in the case of SNORA33 KO cell pool indicates even more intensive hydrazine cleavage at this position compared to other neighboring U residues. The locations of differentially ψ sites identified in SNORA33-depleted cell pools mapped within a human ribosome. 28S-U3801 and 28S-U3741 are located in domain IV, and 28S-U4975 and 28S-U4966 are situated within domain VI of the large ribosomal subunit. * p < 0.05, ** p < 0.01, *** p < 0.001.
Ijms 24 12578 g003
Figure 4. SNORA33 depletion does not affect proliferation or translation capacity. (A) Growth rate of SNORA33-depleted cells measured by crystal violet staining. Data (n = 4) are presented as fold change to Day 0 and were analyzed by two-way ANOVA. (B) The total rate of protein synthesis of SNORA33-depleted cells (n = 3) measured by [35S]methionine/cysteine incorporation. Data are presented as fold changes to CRISPR control. Statistical significance was assessed by unpaired t-test with the assumption of a normal distribution of the data. (C) Polysome profiles of CRISPR control and SNORA33-depleted cells (n = 1).
Figure 4. SNORA33 depletion does not affect proliferation or translation capacity. (A) Growth rate of SNORA33-depleted cells measured by crystal violet staining. Data (n = 4) are presented as fold change to Day 0 and were analyzed by two-way ANOVA. (B) The total rate of protein synthesis of SNORA33-depleted cells (n = 3) measured by [35S]methionine/cysteine incorporation. Data are presented as fold changes to CRISPR control. Statistical significance was assessed by unpaired t-test with the assumption of a normal distribution of the data. (C) Polysome profiles of CRISPR control and SNORA33-depleted cells (n = 1).
Ijms 24 12578 g004
Figure 5. Depletion of SNORA33 affects ribosomal protein composition. (A) The scheme of the experimental set-up. Cytoplasmic extracts of SNORA33-depleted and CRIRP Ctrl cells (n = 3) were depleted of mitochondria and nuclei and ultracentrifuged through the sucrose cushions. Ribosomal pellets were resuspended and analyzed by label-free LC-MS/MS. (B) A pie chart of all identified ribosomal proteins and their functional allocation. (C) Translation factors identified in the ribosomal proteome. (D) A volcano plot of all identified ribosomal proteins. Normal distribution of the data and equal variance of populations was assumed. Statistical significance was assessed by one-way ANOVA. The dotted lines represent cut-off values (FC ≥ 1.1; p < 0.05). Ribosomal proteins significantly less or more abundant in ribosomes of SNORA33-depleted cells are in red and blue, respectively. (E) Protein abundance of selection of differentially abundant core and accessory ribosomal proteins in SNORA33 KO cell pool. Data are presented as fold changes to CRISPR control. * p < 0.05, ** p < 0.01.
Figure 5. Depletion of SNORA33 affects ribosomal protein composition. (A) The scheme of the experimental set-up. Cytoplasmic extracts of SNORA33-depleted and CRIRP Ctrl cells (n = 3) were depleted of mitochondria and nuclei and ultracentrifuged through the sucrose cushions. Ribosomal pellets were resuspended and analyzed by label-free LC-MS/MS. (B) A pie chart of all identified ribosomal proteins and their functional allocation. (C) Translation factors identified in the ribosomal proteome. (D) A volcano plot of all identified ribosomal proteins. Normal distribution of the data and equal variance of populations was assumed. Statistical significance was assessed by one-way ANOVA. The dotted lines represent cut-off values (FC ≥ 1.1; p < 0.05). Ribosomal proteins significantly less or more abundant in ribosomes of SNORA33-depleted cells are in red and blue, respectively. (E) Protein abundance of selection of differentially abundant core and accessory ribosomal proteins in SNORA33 KO cell pool. Data are presented as fold changes to CRISPR control. * p < 0.05, ** p < 0.01.
Ijms 24 12578 g005
Figure 6. SNORA33 depletion affects ribosome function and cellular proteome. (A) The effect of SNORA33 depletion on translational modus and fidelity. Dual-luciferase reporters were used to measure the mode of translation initiation from IRES elements (CrPV, HCV, and P53) and translation fidelity (UAA—stop codon skipping, and −1 PRF—frameshift). Data (n = 9) are corrected for the β-galactosidase activity of the co-transfected lacZ gene and are plotted as fold changes to CRISPR Ctrl. The D’Agostino–Pearson test was used to assess normality and data were analyzed by a nonparametric Mann–Whitney U test. (B) LC-MS/MS analysis of the cellular proteome of SNORA33-depleted cells and CRISPR controls. A volcano plot of all proteins identified in the cellular proteome (n = 4). Normal distribution of the data and equal variance of populations were assumed. Statistical significance was assessed by one-way ANOVA. The dotted lines represent cut-off values (FC ≥ 1.5; p < 0.05). Cellular proteins significantly downregulated and upregulated in SNORA33-depleted cells are in red and blue, respectively. (C) Protein abundance of significantly DE inflammatory factors in SNORA33 KO cell pool. Data are presented as fold changes to CRISPR control. Ns – not significant, * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 6. SNORA33 depletion affects ribosome function and cellular proteome. (A) The effect of SNORA33 depletion on translational modus and fidelity. Dual-luciferase reporters were used to measure the mode of translation initiation from IRES elements (CrPV, HCV, and P53) and translation fidelity (UAA—stop codon skipping, and −1 PRF—frameshift). Data (n = 9) are corrected for the β-galactosidase activity of the co-transfected lacZ gene and are plotted as fold changes to CRISPR Ctrl. The D’Agostino–Pearson test was used to assess normality and data were analyzed by a nonparametric Mann–Whitney U test. (B) LC-MS/MS analysis of the cellular proteome of SNORA33-depleted cells and CRISPR controls. A volcano plot of all proteins identified in the cellular proteome (n = 4). Normal distribution of the data and equal variance of populations were assumed. Statistical significance was assessed by one-way ANOVA. The dotted lines represent cut-off values (FC ≥ 1.5; p < 0.05). Cellular proteins significantly downregulated and upregulated in SNORA33-depleted cells are in red and blue, respectively. (C) Protein abundance of significantly DE inflammatory factors in SNORA33 KO cell pool. Data are presented as fold changes to CRISPR control. Ns – not significant, * p < 0.05, ** p < 0.01, *** p < 0.001.
Ijms 24 12578 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chabronova, A.; van den Akker, G.; Housmans, B.A.C.; Caron, M.M.J.; Cremers, A.; Surtel, D.A.M.; Peffers, M.J.; van Rhijn, L.W.; Marchand, V.; Motorin, Y.; et al. Depletion of SNORA33 Abolishes ψ of 28S-U4966 and Affects the Ribosome Translational Apparatus. Int. J. Mol. Sci. 2023, 24, 12578. https://doi.org/10.3390/ijms241612578

AMA Style

Chabronova A, van den Akker G, Housmans BAC, Caron MMJ, Cremers A, Surtel DAM, Peffers MJ, van Rhijn LW, Marchand V, Motorin Y, et al. Depletion of SNORA33 Abolishes ψ of 28S-U4966 and Affects the Ribosome Translational Apparatus. International Journal of Molecular Sciences. 2023; 24(16):12578. https://doi.org/10.3390/ijms241612578

Chicago/Turabian Style

Chabronova, Alzbeta, Guus van den Akker, Bas A. C. Housmans, Marjolein M. J. Caron, Andy Cremers, Don A. M. Surtel, Mandy J. Peffers, Lodewijk W. van Rhijn, Virginie Marchand, Yuri Motorin, and et al. 2023. "Depletion of SNORA33 Abolishes ψ of 28S-U4966 and Affects the Ribosome Translational Apparatus" International Journal of Molecular Sciences 24, no. 16: 12578. https://doi.org/10.3390/ijms241612578

APA Style

Chabronova, A., van den Akker, G., Housmans, B. A. C., Caron, M. M. J., Cremers, A., Surtel, D. A. M., Peffers, M. J., van Rhijn, L. W., Marchand, V., Motorin, Y., & Welting, T. J. M. (2023). Depletion of SNORA33 Abolishes ψ of 28S-U4966 and Affects the Ribosome Translational Apparatus. International Journal of Molecular Sciences, 24(16), 12578. https://doi.org/10.3390/ijms241612578

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop