Mass-spectrometry (MS) - based screening revealed several proteins bound to FMR1 RNA with expanded CGG repeats; some affect the yield of polyglycine-containing proteins, alleviating their toxicity.

(A) Scheme of two RNA molecules used to perform MS-based screening. The FMR1 RNA contains the entire length of the 5’ untranslated region (5’UTR) of FMR1 with expanded CGG (x99) repeats forming a hairpin structure (red). The open reading frame for the polyglycine-containing protein starts at a repeats-associated non-AUG initiated (RAN) translation-specific ACG codon. RAN translation can also be initiated from a GUG near-cognate start codon, which is not indicated in the scheme. The GC-rich RNA contains TMEM107 mRNA enriched with G and C nucleotide residues (GC content > 70%; similar to FMR1 RNA) with the open reading frame starting at the canonical AUG codon (blue). Both RNAs are tagged with three MS2 stem-loop aptamers (grey) interacting with an MS2 protein tagged with an in vivo biotinylating peptide used to pull down proteins interacting with the RNAs.

(B) The volcano plot representing proteins captured during MS-based screening showing the magnitude of enrichment (log2 fold change) and the statistical significance ( log P-value); red dots indicate proteins significantly enriched (P < 0.05) on FMR1 RNA compared to GC-rich RNA. Three proteins, RPS26, LUC7L3, and DHX15, tested in a subsequent validation experiment are marked. DDX3X is also indicated as this protein has been previously described in the context of interaction with FMR1 mRNA.

(C) The scheme of RNA used for the transient overexpression of FMR99xG (mutant, long, 99 polyglycine tract-containing protein). The construct contains the entire length of the 5’UTR of FMR1 with expanded CGG repeats forming a hairpin structure (red) tagged with enhanced green fluorescent protein (eGFP) (green). Western blot analysis of FMR99xG and Vinculin for HEK293T cells with insufficient DHX15, RPS26, and LUC7L3 induced by specific short interfering RNA (siRNA) treatment. To detect FMR99xG, the 9FM antibody was used. The upper bands were used for quantification. The graph presents the mean signal for FMR99xG normalized to Vinculin from N = 3 biologically independent samples with the standard deviation (SD). An unpaired Student’s t-test was used to calculate statistical significance: *, P < 0.05; ***, P < 0.001; ns, non-significant.

(D) Results of luciferase assay for cells with overexpression of FMR16xG fused with nanoluciferase pretreated with different concentrations of siCtrl or siRPS26 (see legend). The graph presents luminescence signals (means from N = 8 with SDs) normalized to siCtrl treated cells. Western blot represents the siRNA-concentration dependent effect on RPS26 downregulation. An unpaired Student’s t-test was used to calculate statistical significance: ****, P < 0.0001; ns, non-significant.

(E) Results of microscopic quantification of FMR99xG-positive aggregates (white arrows) in HeLa cells transfected with the FMR99xG construct on the background of normal (siCtrl) or reduced amount of RPS26 (siRPS26). Representative images were pseudo-colored and merged; green, GFP-positive inclusions; blue, nuclei stained with Hoechst 33342; scale bars, 50 μm. The histogram presents a number of FMRpolyG-GFP inclusions per nucleus in cells treated with siCtrl (blue) or siRPS26 (red) in N = 10 biologically independent samples with SDs. In average, 1100 cells were identified per single image and final value were normalized to 1 for control siRNA treated cells. An unpaired Student’s t-test was used to calculate statistical significance: *, P < 0.05.

(F) The influence of RPS26 silencing on apoptosis evoked in HeLa cells after 28h or 43h of FMR99xG overexpression or in mock treated cells. Apoptosis was measured as luminescence signals (relative luminescence units; RLU) derived from Annexin V fusion protein bound to phosphatidylserine (PS) exposed on the outer leaflet of cell membranes, an indicator for early apoptosis. The graph presents relative mean values from N = 6 biologically independent samples treated with either siCtrl (red) or siRPS26 (blue) with the SD normalized to mock controls (cells transfected only with the delivering reagent and cultured for 28h or 43h). An unpaired Student’s t-test was used to calculate statistical significance: *, P < 0.05; **, P < 0.01; ***, P < 0.001. Note, that significantly higher apoptosis was observed in cells expressing FMR99xG compare to control (mock), and that pretreatment with siRPS26 significantly reduced this phenotype.

RPS26 depletion affects proliferation of cells but does not inhibit global translation.

(A) The left graph represents the number of cells and right graph represents percentage (%) of dead cells positively stained with trypan blue after 48h post siRNA-induced silencing. Prior siRNA transfection equal amount of HEK293T cells were seeded (ca. 1 × 104). The graphs present means from N = 4 with SDs. FC, fold change. An unpaired Student’s t-test was used to calculate statistical significance: **, P < 0.01; ns, non-significant. Note, that the number of cells in samples treated with siRPS26 was lowered by 34% compared to siCtrl treated samples.

(B) Western blot analysis of the surface sensing of translation (SUnSET) assay 48h post RPS26 silencing (siRPS26) in HEK293T cells, with antibody against puromycin (PURO), to monitor efficiency of PURO incorporation into peptides during translation. Cells were incubated with siRPS26 or siCtrl and 10 min before protein isolation PURO was added. As a positive control for efficient translation inhibition cells were treated with cycloheximide (CHX), an inhibitor of translation elongation (CHX was added 5 min prior PURO treatment). As a negative control cells treated with DMSO only were used (negative for both PURO and CHX). The graph represents mean signal from whole lanes from N = 3 independent samples with SDs. An unpaired Student’s t-test was used to calculate statistical significance: ns, non-significant.

RPS26 insufficiency induces a lower production of polyglycine, but not FMRP, in multiple FXTAS cellular models.

(A) Representative microscopic images showing inducible expression of FMR95xG fused with eGFP in stable transgenic cell line containing a single copy of 5’UTR FMR1 99xCGG-eGFP transgene under the control of a doxycycline-inducible promoter: S-95xCGG model. Blue are nuclei stained with Hoechst; the green signal is derived from FMR95xG tagged with eGFP; scale bar 50 µm; +DOX and DOX indicate cells treated (or not) with doxycycline to induce transcription of the transgene from the doxycycline-dependent promoter.

(B & C) Results of western blot analyses of FMR95xG (with long, mutant polyglycine stretches) or FMR16xG (with short, normal polyglycine stretches), both fused with GFP, normalized to Vinculin and an RT-qPCR analysis of FMR1-GFP transgene expression normalized to GAPDH upon RPS26 silencing in S-95xCGG and S-16xCGG models, respectively. siRPS26_I and siRPS26_II indicate two different siRNAs used for RPS26 silencing. The graphs present means from N = 3 (B) or N = 5 biologically independent samples (C) with SDs. An unpaired Student’s t-test was used to calculate statistical significance: *, P < 0.05; **, P < 0.01; ns, non-significant.

(D) Representative microscopic images of cells stably expressing FMR99xCGG obtained after transduction with lentivirus containing the 5’UTR FMR1 99xCGG-eGFP transgene: L-99xCGG model. Nuclei are stained with Hoechst and green fluorescence signal is derived from FMR99xG tagged with GFP. Scale bar 50 µm.

(E) Results of western blot analyses of FMR99xG and FMRP normalized to Vinculin upon RPS26 silencing in the L-99xCGG model. The graphs present means from N = 3 biologically independent samples with SDs. An unpaired Student’s t-test was used to calculate statistical significance: *, P < 0.05; **, P < 0.01; ns, non-significant.

(F) Results of western blot analyses of FMRP levels normalized to Vinculin upon RPS26 silencing in fibroblasts derived from two males, either from healthy individual (CGGnorm – XY, 31xCGG) or FXTAS patient (CGGexp – XY, 81xCGG). The graphs present means from N = 3 biologically independent samples with SDs. An unpaired Student’s t-test was used to calculate statistical significance: ns, non-significant.

Changes in the proteome of RPS26-deficient cells are not robust.

(A) The volcano plot represents a stable isotope labeling using amino acids in cell culture (SILAC)-based quantitative proteomic analysis identifying proteins sensitive to RPS26 insufficiency. It shows the magnitude of protein-level changes (log2 fold change) vs. the statistical significance ( log P-value) 48h post siRPS26 treatment. Data were collected from three independent biological replicates for each group. Grey dots indicate proteins non-responding to RPS26 depletion (N = 1506). Red dots indicate proteins responding to RPS26 insufficiency (P < 0.05); EIF5 and PDCD4 are examples of Negative (N = 223; green) and Positive (N = 158; orange) responders, respectively. Protein groups were further analyzed in B, C, and D.

(B) Gene ontology (GO) analysis performed for positive responders of RPS26 insufficiency from the proteomic experiment shown in A. The graph presents significantly enriched GO terms (P < 0.05); statistical significance was calculated using Fisher’s Exact test with the Bonferroni correction. Note that for negative responders, no GO terms were significantly enriched.

(C) Data validation from the proteomic experiment described in A. Western blot analyses of PDCD4, ILF3, FMRP, and EIF5 proteins normalized to Vinculin for HEK293T cells treated with siRPS26. The graphs present means from N = 4 biologically independent samples with SDs. An unpaired Student’s t-test was used to calculate statistical significance: *, P < 0.05; **, P < 0.01; ***, P < 0.001; ns, non-significant.

(D) The percentage of GC content in 5’UTRs and coding sequences across extending sequence windows initiated from the start codon within three groups of transcripts: Negative responders (green), Positive responders (orange), and Background (BG) understood as the total transcriptome (red). To avoid biases in the GC-content mean, transcripts with 5’UTR sequences shorter than 20 nucleotides were excluded from the analysis, yielding the following sample sizes: Negative responders (N = 213), Positive responders (N = 147), and BG (N = 20,862). For example, position 6 in 5’ UTR sequences corresponds to a 6-nucleotide fragment (window from 6 to 1 positions upstream of ATG), while position 7 corresponds to a 7-nucleotide fragment (window from 7 to 1 positions). The solid line shows the mean GC content at a given position (i.e., within the window), and the shade indicates the standard error of the mean. P-values < 0.01 are denoted by green or yellow dots. These reflect pairwise comparisons of GC content between transcript groups (compared to BG) and were determined using a two-tailed paired t-test with Bonferroni correction.

(E) Box plot showing the results of a one-tailed Mann-Whitney U test comparing 5’UTR sequence lengths between positive and negative responders relative to the background transcriptome. To ensure a balanced comparison, the background transcriptome is represented by 1000 randomly selected transcripts from the complete human transcriptome. The plot highlights the median (indicated by the central line), the upper and lower quartiles (represented by the boxes), and the highest and lowest values, which are determined within a range of 1.5 times the interquartile range (IQR) (whiskers).

(F) Results of luciferase assay post RPS26 silencing showing differences in expression of transgene containing 5’UTR of FMR1 with 16xCGG fused with nanoluciferase from which RAN translation was initiated from either ACG (wild type) or GTG or CTG near cognate codons. The signals for activity of nanoluciferase fused with FMR16xG were normalized to firefly luciferase signals (means from N = 4 with SDs are presented on the graph). Note that firefly luciferase was translated from canonical ATG start codon. An unpaired Student’s t-test was used to calculate statistical significance: **, P < 0.01; ns, non-significant.

(G) Western blot analysis of FMR99xG post RPS26 silencing and treatment with control ASO (ASO_CTRL) or ASO targeting ACG codon, RAN translation initiation codon in 5’UTR of FMR1 (ASO_ACG). The graph represents the means of N = 3 with SDs. An unpaired Student’s t-test was used to calculate statistical significance: *, P < 0.05; **, P < 0.01; ***, P < 0.001; ns, non-significant.

Insufficiency of the TSR2 chaperone protein lowers the level of RPS26 and FMRpolyG but not FMRP and selected 40S components.

(A) Western blot analyses of RPS26, FMR95xG, FMRP but also Histone H3.3, a sensor of ribosomes depleted with RPS26, normalized to Vinculin in stable S-95xCGG cells treated with siTSR2. Graphs represent means from N = 5 biologically independent samples with SDs. An unpaired Student’s t-test was used to calculate statistical significance: **, P < 0.01; ***, P < 0.001; ns, non-significant.

(B) Western blot analyses of selected 40S ribosomal proteins RACK1, RPS6, and RPS15, normalized to Vinculin upon TSR2 silencing in the S-95xCGG model. The graphs present means from N = 5 biologically independent samples with SDs. An unpaired Student’s t-test was used to calculate statistical significance: *, P < 0.05; ns, non-significant.

(C) Western blot analysis of the results of SUnSET assay (upper) and the level of FMR99xG (lower) 48h post TSR2 or RPS6 silencing in HEK293T cells. CHX treatment served as a positive control for translation inhibition. Negative control (PURO negative and CHX negative) was protein extract from cells treated with DMSO only (see Figure 2 for more details). Graphs represent means form N = 3 independent samples with SDs. An unpaired Student’s t-test was used to calculate statistical significance: **, P < 0.01; ns, non-significant. Note that presented blots are cropped. Full size images are deposited on Zenodo server under the DOI: 10.5281/zenodo.13860369.

(D) Graphs represent the number of cells post 48h of TSR2 or RPS6 silencing. Prior siRNA transfection equal amount of HEK293T cells were seeded (∼1 × 104). The graphs present means from N = 4 with SDs. FC, fold change. An unpaired Student’s t-test was used to calculate statistical significance: *, P < 0.05; **, P < 0.01.

Silencing RPS25, the other component of 43S PIC, reduces the level of polyglycine produced from mutant FMR1 mRNA.

(A) A western blot analysis of FMR95xG normalized to Vinculin and an RT-qPCR analysis of the FMR1-GFP transgene expression normalized to GAPDH upon RPS25 silencing in the S-95xCGG model. The graphs present means from N = 3 biologically independent samples with SDs. An unpaired Student’s t-test was used to calculate statistical significance: *, P < 0.05; ns, non-significant.

(B) The graph represents the number of cells post 48h of RPS25 silencing. Prior siRNA transfection equal amount of HEK293T cells were seeded (∼1 × 104). The graph presents means from N = 4 with SDs. FC, fold change. An unpaired Student’s t-test was used to calculate statistical significance: ****, P < 0.0001.

(C) The graph represents viability of SH-SY5Y neuroblastoma cells measured by luminescence based on the reducing potential of viable cells. SH-SY5Y cells were treated with all siRNAs tested in this study. The graph presents means from N = 4 biologically independent samples with SDs. An unpaired Student’s t-test was used to calculate statistical significance: *, P < 0.05; **, P < 0.01.

(D) A western blot analysis of FMR99xG normalized to Vinculin from the human neuroblastoma cell line (SH-SY5Y) upon silencing of either RPS26 or RPS25. siRPS26_I and siRPS26_II indicate two different siRNAs used for RPS26 silencing. The upper bands were used for quantification. The graph presents means from N = 3 biologically independent samples with SDs. An unpaired Student’s t-test was used to calculate statistical significance: *, P < 0.05; **, P < 0.01; ns, non-significant.

MS-based screening revealed proteins interacting with FMR1 RNA.

(A) Sequences of FMR1 RNA and GC-rich RNA constructs used in MS2-based screening. Note that the sequences presented are a fragment of the plasmids used in the experiment.

(B) Results of western blot representing FMRpolyG derived from FRM1 RNA construct (lane 1&2) in comparison to FMRpolyG tagged with GFP (lane 3&4) derived from template construct (Addgene #63089). Tubulin was used as a loading control.

(C) Gene Ontology (GO) analysis performed on proteins detected in FMR1 RNA sample group (common among three technical replicates). Graph presents significantly enriched, selected GO terms.

(D) Results of western blot analysis of FMR99xG normalized to Vinculin upon NOP58, HSHP1, NAF1, MYBBP1A and DHX9 siRNA-based silencing in HEK293T cells. Graphs present means from N = 3 biologically independent samples with SDs. Unpaired Student’s t test was used to calculate statistical significance: ns, non-significant. Note that presented blots are cropped.

(E) RT-qPCR analysis of FMR1-GFP transgene expression derived from transient transfection system normalized to GAPDH upon RPS26 and DHX15 depletion in HEK293T cells. siRPS26_I and siRPS26_II indicate two independent siRNAs. Graphs present means from N = 3 biologically independent samples with SDs. Unpaired Student’s t test was used to calculate statistical significance: ns, non-significant

(F) Results of western blot demonstrating enrichment of RPS26 detected in 5’UTR 99xCGG and 23xCGG eluates compared to GC-rich RNA in biochemical assay. 50% of eluate fraction was loaded on the gel. All RNA molecules in the experiment were biotinylated (b) in order to perform RNA-protein pull down procedure from HEK293T cell extract. No RNA sample served as a negative control.

(G) Results of western blot demonstrating enrichment of SAM68 protein detected in FMR1 RNA sample eluate after MS2-based, in cellulo RNA-protein pull down procedure; the percentage of following fractions were loaded on the gel: Input, 5% of total lysate; eluate, 20% of immunoprecipitated fraction. Note that presented blot is cropped.

Characterization of S-95xCGG and S-16xCGG models and RPS26-responders.

(A) Results of western blot of FMR95xG and FMR16xCGG expression in S-95xCGG and S-16xCGG models, respectively, after induction of doxycycline-inducible promoter. Upper image represents the titration of doxycycline (in ng/ml). Following images demonstrate time-dependent course of FMR95xG/FMR16xCGG expression post promoter induction. The same amount of protein sample was loaded on gels. Vinculin serves as a loading control.

(B) The representative confocal microscopic images showing homogenous expression of FMR95xG fused with eGFP in S-95xCGG model after 35 days of doxycycline treatment; scale bar 10 µm. Note that no GFP-positive aggregates of polyglycine were detected.

(C) RT-qPCR analysis of endogenous FMR1 and FMR1-GFP transgene expression in S-95xCGG and S-16xCGG models with (+DOX) or without (-DOX) doxycycline induction. On the left, bars represent the mean cycle (Ct) for S-95xCGG model, and on the right, bars indicate relative FMR1-GFP mRNA level normalized to GAPDH for both models. Unpaired Student’s t test was used to calculate statistical significance: **, P < 0.01.

(D) Images of L-99xCGG model were taken with Nikon A1 confocal microscope, equipped with GaAsP Detectors. GFP was excited with 488 nm ion argon Laser and detected with 525/50 Detector. Arrows indicate FMRpolyG-GFP-positive aggregates (green); blue are nuclei; scale bar 10 µm.

(E) Results of western blot analysis of FMR95xG normalized to Vinculin upon RPS26 silencing in S-95xCGG. siRPS26_I and siRPS26_II indicate two different siRNAs used in the study. The graphs present means from N = 3 with SDs. An unpaired Student’s t-test was used to calculate statistical significance: *, P < 0.05.

(F) Results of western blot analysis of SUZ12 and Histone H3.3 proteins normalized to Vinculin upon RPS26 silencing in COS7 cells. Graphs present means from N = 3 biologically independent samples with SDs. Unpaired Student’s t test was used to calculate statistical significance: *, P < 0.05; **, P < 0.01.

Validation of SILAC-MS results and bioinformatic analysis of mRNAs encoding RPS26-sensitive proteins.

(A) Validation of SILAC-MS data. Results of western blot analysis of RPS6 & FUS (results are in line with proteomic analysis) and PCBP2 & EIF3J (results are not in line with proteomic analysis) normalized to Vinculin upon RPS26 silencing in HEK293T cells. Graphs present means from N = 4 biologically independent samples with SDs. Unpaired Student’s t test was used to calculate statistical significance: *, P < 0.05; ns, non-significant.

(B) Frequency of nucleotides across transcript sequence positions in the close vicinity of start codon within three groups of transcripts: negative responders (N = 218), positive responders (N = 158), and background transcripts (N = 22,160). The analyzed sequence positions span the last 20 nucleotides of the 5’ UTR sequence (from 20 to 1 downstream of the start codon) and the first 20 nucleotides of the coding sequence (from +1 to +20 upstream of the start codon). Each nucleotide symbol’s height on the graph is proportionate to its relative frequency at the corresponding sequence position. The cumulative height of the stack at each position sums to one, encompassing the frequencies of all four nucleotides.

(C) The percentage of GC content of 5’UTR across extending sequence windows initiated upstream from the start codon within groups of transcripts. The graph gathers analysis of transcripts determined with two P-value cut offs (P < 0.05 and P < 0.01) yielding different samples sizes for each group. Colors corresponding to given group as well as its size are indicated in the legend. The solid line shows the mean GC content at a given position (i.e., within the window), and the shade indicates the standard error of the mean. adj. P values (right panel) reflect pairwise comparisons of GC content between transcript groups and were determined using a two-tailed paired t-test with Bonferroni correction.

(D) Results of western blot analysis of FMRpolyG protein derived from wild type (WT) or mutated ( 4G<A) construct normalized to Vinculin upon RPS26 silencing in HEK293T cells. Graphs present means from N = 3 biologically independent samples with SDs. Unpaired Student’s t test was used to calculate statistical significance: ns, non-significant.

(E) Comparison of overrepresented hexamers in 5’UTR full-length sequences within two transcript groups: negative (N = 218) and positive responders (N = 158). The frequency of each hexamer is compared to the corresponding hexamer frequency in 5’UTR sequences of background transcripts (BG, N = 22,160) (indicated as log2 fold change). The P-value associated with each hexamer signifies the degree of overrepresentation (enrichment) and is calculated from the binomial distribution. Only hexamers with P-values < 10-5 are displayed.

Silencing of TSR2 reduces FMR99xG and FMR16xG level

(A) Results of western blot analysis of FMR99xG normalized to Vinculin upon TSR2 silencing in COS7 cell line. The graphs present means from N = 3 with SDs. An unpaired Student’s t-test was used to calculate statistical significance: *, P < 0.05.

(B) Results of western blot analysis of FMR16xG normalized to Vinculin upon TSR2 silencing in S-16xCGG model. Graphs present means from N = 5 biologically independent samples with SDs. Unpaired Student’s t test was used to calculate statistical significance: ***, P < 0.001.

RPS25 co-precipitates with FMR1 and affects FMR99xG level

(A) Results of western blot analysis of FMR99xG level upon RPS25 silencing in SH-SY5Y cell line.

(B) Results of western blot demonstrating enrichment of RPS25 detected in samples obtained with biotinylated 5’UTR 99xCGG RNA and to a lesser extend in 5’UTR (no repeats) & GC-rich RNA eluates (ca. 20% of eluate fraction was loaded on the gel). Input, 5% of total protein lysate was loaded. All RNA molecules in the experiment were biotinylated in order to perform RNA-protein pull down procedure. Beads only was a sample (with no RNA) which served as a negative control