RNA

Small RNAs (sRNAs) regulate the expression of many genes including those involved in antibiotic resistance and bacterial virulence, making them potential therapeutic targets. A molecule that binds an sRNA could interfere with its ability to bind its target mRNA and disrupt the regulation mechanism. Randomization and screening of natural arginine rich motif (ARM) peptides led to peptides capable of interfering with the sRNA MicF's ability to regulate ompF in Escherichia coli. Molecular dynamics simulations suggested that this effect was not a result of a direct disruption of the MicF-ompF interaction. Instead, the peptides interfere with binding of the chaperone Hfq, which is required for MicF-mediated regulation. Subsequent testing demonstrated peptide specificity for MicF over two other Hfq scaffolds and the ability to disrupt regulation of two additional MicF targets. Together, these findings support the use of synthetic ARMs as a potential tool for modulating sRNA function in bacteria.

Satellite RNAs (satRNAs) are parasitic subviral agents whose biological activities are mediated largely by specific sequence determinants and structured RNA elements. Establishing biologically supported global RNA structures is therefore essential for understanding how RNA architecture underlies satRNA functions. Cucumber mosaic virus (CMV) satRNA is one of the best-characterized models for investigating satRNA structure-function relationships; however, a biologically supported global RNA architecture of CMV satRNA has yet to be established. Here, we applied AlphaFold3 modeling to predict the global structure of CMV satRNA T1 (sat-T1). Initial full-length structure modeling generated multiple long-distance interactions that lacked biological support. We therefore used fragment-based modeling, combined with chemical probing, evolutionary covariation, and compensatory mutagenesis, to derive a biologically supported global secondary structure. To determine whether structurally distant regions could interact in the context of the full-length RNA, we engineered a structure-guided T1-ZD mutant that preserved the supported secondary structure while reducing alternative base-pairing possibilities. Full-length AlphaFold3 modeling of T1-ZD largely recapitulated the proposed architecture, while one predicted model revealed a long-distance interaction that was subsequently supported by compensatory mutagenesis analysis. These findings suggest that the 3' terminus of sat-T1 may undergo conformational switching between alternative structural states. Together, our work establishes a biologically supported global RNA architecture for CMV sat-T1 and provides a structural framework for investigating the molecular basis of satRNA function.

Intrinsically disordered regions are widespread in RNA-processing machines. In Escherichia coli, RNase E uses its intrinsically disordered C-terminal domain (CTD) to recruit RNAs to the N-terminal catalytic domain, including mRNAs targeted by regulatory small RNAs (sRNAs), but the basis of substrate recognition and specificity is unclear. We engineered a protease-cleavable RNase E and used split-CRAC to isolate RNAs crosslinked to the AR2 sub-domain of the intrinsically disordered CTD fragment. AR2 preferentially engaged mRNAs and was depleted of sRNAs and sRNA-containing hybrids, supporting recognition of the mRNA. AR2 contacts concentrated on accessible A-rich motifs surrounding ribosome-binding sites and start codons, and purified AR2 recognised this motif in vitro. AR2 also contacted an AUAA motif in the rne translation-initiation region, and AR2 deletion increased RNase E abundance implicating this interaction in autoregulation. These findings define a relatively short AR2-binding motif and are consistent with CTD interactions with the 30S subunit that may provide additional specificity for a subset of mRNA translation initiation regions. SIGNIFICANCE STATEMENTMost RNA turnover in bacterial cells is carried out by the RNA degradosome, yet how this molecular machine checks and selects RNAs for degradation remains incompletely understood. We show that an intrinsically disordered region of the degradosome enzyme RNase E, termed AR2, preferentially binds A-rich sequences near sites of translation initiation. Through recognition of this shared sequence feature in a common functional context, AR2 may help the degradosome recognise messenger RNAs as a functional class. AR2 also contributes to feedback control of RNase E expression by recognising its own messenger RNA.

Background: RNA molecules form thermodynamic ensembles, but interpretation often requires a single representative structure. Existing base-pair centroid estimators assess agreement at the level of individual base pairs and do not directly target nesting depth along the sequence. Methods: We introduce Mountain Centroid, which minimizes expected squared mountain-profile distance, and derive dynamic programming algorithms with and without RNA pairing constraints. We also combine the Mountain Centroid objective with the base-pair centroid gain. Results: Across 21,254 RNAStrAlign sequences, Mountain Centroid had lower median normalized mean squared mountain distance (NMSMD) than minimum-free-energy (MFE) and base-pair centroid ({gamma} = 1) structures, whereas its median base-pair F1 was lower. Imposing RNA pairing constraints improved base-pair F1 for 59.35% of sequences and reduced it for 3.58%. At an illustrative weight, the combined objective had median base-pair F1 similar to MFE while retaining lower median NMSMD than MFE and all tested {gamma}-centroid settings. Conclusions: Mountain Centroid represents an RNA structural ensemble with a single secondary structure that reflects how nesting depth varies across nucleotide positions. Combining mountain-profile and individual-base-pair criteria allows their relative contributions to be varied.

Previous studies revealed a crosstalk between intronic polyadenylation (IPA) and the DNA damage response (DDR). Indeed, genotoxic agents, including radiations and anticancer drugs (e.g., cisplatin that crosslinks DNA), regulate the ratio of IPA to last-exon transcripts in many genes. Conversely, multiple genes involved in the DDR, especially homologous recombination, are regulated at the IPA level. The U1 small nuclear RNA (snRNA) widely represses IPA, thereby enhancing full-length gene transcription. However, besides its implication in IPA regulation by ultraviolet-C radiation, little is known about U1 snRNA effects on the DDR and on cell sensitivity to genotoxic agents. Here, we show that U1 snRNA blockade using an antisense oligonucleotide (U1-AMO) in lung cancer cell lines enhances cell growth inhibition by cisplatin, through an increase in cisplatin-induced DNA damage. 3-seq analysis indicates that U1 snRNA blockade represses full-length mRNA expression of multiple genes of the nucleotide-excision repair and Fanconi anemia pathways, which are involved in the repair of cisplatin-DNA crosslinks. Our 3-seq analyses also reveal that moderate doses of U1-AMO and cisplatin upregulate the IPA:LE isoform ratio in overlapping but distinct sets of genes, and that U1-AMO prevents cisplatin effects on the IPA:LE ratio in a large subset of genes. Altogether, these data extend the crosstalk between IPA and the DDR and suggest that U1 snRNA targeting may be used to sensitize cancer cells to genotoxic agents.

RNA recognition motif (RRM) proteins frequently contain multiple RNA-binding domains connected by flexible linkers, yet the contribution of transient interdomain interactions to RNA recognition remains incompletely understood. Here, we investigated the structural organization of the tandem RRMs of the Drosophila melanogaster splicing regulator Sex-lethal (Sxl) using solution NMR spectroscopy in combination with rational protein engineering, restrained docking and RNA-binding studies. Progressive extension of the native interdomain linker resulted in a gradual decrease in rotational coupling between the two RRMs and continuous chemical shift changes, demonstrating that the RNA-free protein samples a dynamic conformational ensemble rather than behaving as two independently tumbling domains. NMR-guided docking identified a compact arrangement compatible with the experimental data and suggested a transient interface partially overlapping the RNA-binding surfaces. Surprisingly, a mutant designed to weaken this interface produced the opposite effect: instead of increasing interdomain mobility, it exhibited enhanced rotational coupling while remaining natively folded, indicating a redistribution of the conformational ensemble rather than disruption of the domain architecture. Both linker extension and the mutant reduced RNA-binding affinity, and the mutant additionally diminished sequence discrimination, demonstrating that perturbations shifting the conformational equilibrium in either direction compromise RNA recognition. Together, our results demonstrate that RNA recognition by Sxl is governed not by a single apo structure but by a finely balanced conformational ensemble, and that perturbing this equilibrium in either direction compromises high-affinity and sequence-selective RNA binding.

N2-methylguanosine (m2G) is widely found at multiple positions in tRNAs across the three domains of life. Tryptophan tRNA from Thermococcus kodakarensis contains m2G at position 67. We previously proposed that the tRNA m2G methyltransferase Trm14 is responsible for m2G67 formation in tRNATrp from T. kodakarensis, although Trm14 was originally identified as the enzyme catalyzing m2G6 formation in tRNACys in Methanocaldococcus jannaschii. Thus, it remained unclear whether Trm14 could also methylate G67. Here, we characterized archaeal Trm14. Biochemical analyses using recombinant T. kodakarensis Trm14 revealed that the enzyme catalyzes m2G formation at positions 6 and 67 in T. kodakarensis tRNACys and tRNATrp transcripts, respectively. Mass spectrometric analyses demonstrated the loss of m2G6 and m2G67 in native tRNACys and tRNATrp, respectively, from a T. kodakarensis trm14 gene disruptant strain, providing direct evidence for the dual-site specificity of T. kodakarensis Trm14. The growth phenotype of the trm14 gene disruptant strain was comparable to that of the wild-type strain. In contrast, a trm14/trm11 double disruptant, in which trm11 encodes the tRNA m2G10/m22G10 methyltransferase, exhibited severe growth retardation at 95 {degrees}C. This suggests that m2G6/m2G67 and m2G10/m22G10 cooperatively contribute to cellular fitness at high temperatures. Biochemical analyses revealed that Trm14 methylates all 46 T. kodakarensis tRNA transcripts. Furthermore, we found that recombinant M. jannaschii Trm14 methylated both positions. In contrast, the bacterial ortholog TrmN modified only position 6 in tRNA. Overall, this study expands our understanding of archaeal Trm14 by demonstrating its broader substrate specificity and the physiological significance of these modifications under hyperthermophilic conditions.

The majority of eukaryotes encode some intron-containing pre-tRNAs. Splicing of these pre-tRNAs requires a dedicated tRNA splicing machinery. The fungal and trypanosome tRNA ligase, Trl1, and the human RNA ligase, RTCB, catalyze an essential step in tRNA splicing. However, Trl1 and RTCB are nonhomologous and biochemically and structurally distinct from each other. Therefore, Trl1 could serve as a broad-spectrum antifungal and anti-trypanosomal target. While the functions and requirements of the three catalytic Trl1 domains have been extensively characterized in the model yeast Saccharomyces cerevisiae, the roles of Trl1 orthologs in pathogenic fungi remain unexplored. Here, we validate Trl1 as one of the few promising novel drug targets for the development of antifungal therapeutics. Functional analyses of the three Trl1 domains show that only the "sealing" domain is essential for growth and viability in Candida albicans and Aspergillus fumigatus. In contrast, the two "healing" domains are dispensable in these pathogenic fungi, suggesting the presence of redundant healing enzymes, unlike in S. cerevisiae. These findings indicate that only the sealing domain is a good drug target. Our analysis also shows that the Mucor enzyme, which only contains the sealing domain, is essential. Using a Caenorhabditis elegans infection model of C. albicans, we further demonstrated that inhibiting Trl1 expression protects worms during an established infection. In contrast to these fungal pathogens, we show that all three domains of Trl1 are essential in Trypanosoma brucei. Our findings show that the essentiality of the Trl1 sealing is conserved in important human pathogens and provides an impetus for future drug development. SIGNIFICANCEFungal infections are an important cause of human disease and death and difficult to treat and there is an urgent need to develop additional drugs. Based on studies in yeast, one promising target for antifungal drug development is the tRNA splicing pathway. Human tRNA ligase is fundamentally distinct from the fungal one. To investigate the possibility of developing tRNA ligase-targeting drugs, we investigated the function of the catalytic domains of fungal tRNA ligase in different fungal pathogens. Surprisingly, only the first domain is essential in these pathogens and yeast is not a good model fungus. In contrast, all three domains of Trypanosome tRNA ligase are essential. These findings provide an impetus for future drug development.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=79 SRC="FIGDIR/small/745202v1_ufig1.gif" ALT="Figure 1"> View larger version (21K): org.highwire.dtl.DTLVardef@52e3edorg.highwire.dtl.DTLVardef@1f21637org.highwire.dtl.DTLVardef@909ee8org.highwire.dtl.DTLVardef@b0c907_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOGRAPHICAL ABSTRACTC_FLOATNO C_FIG Crosslinking and immunoprecipitation followed by sequencing (CLIP-seq) is widely used to identify the RNA targets of RNA-binding proteins (RBPs). However, its application to non-canonical RBPs lacking canonical RNA-binding domains and frequently displaying low or transient RNA occupancy, is limited by low signal-to-noise ratios, high input requirements and error-prone ligation steps during library preparation. To overcome these limitations, we developed soniCLIP, a streamlined CLIP-seq workflow that replaces RNase-mediated RNA fragmentation with sonication and uses a ligation-free strategy for library construction. soniCLIP is optimized for reproducible identification of enriched RBP-associated RNA regions from limited starting material. We benchmarked soniCLIP against the widespread eCLIP approach and observed reproducible recovery of known RBP-associated regions and target recovery comparable to ENCODE eCLIP, while requiring only 10% (500 {micro}g) of protein input. We further applied soniCLIP to the glycolytic enzyme and non-canonical RBP pyruvate kinase M2 (PKM2). We identified 197 significantly enriched RNA regions and validated selected targets by RIP-qRT-PCR and in vitro binding assays. By combining reduced input requirements, high reproducibility, a shortened 3.5-day workflow and the elimination of gel-based purification, soniCLIP provides an efficient and robust approach for the identification of RNA targets of canonical and non-canonical RBPs.

We present a simple and efficient way to add cis dialdehydes, phosphoryl groups, or cis diols to the 3 prime end of oligonucleotides using periodate oxidized nucleotides (oNTPs) and terminal deoxynucleotidyl transferase (TdT). The 3 prime end cis dialdehyde-modified oligos are generated by incubating TdT with an oligo for several minutes followed by addition of a oNTP and incubated at 30 degrees C for 30 minutes to an hour. After allowing the addition of the cis dialdehydes, heating the reaction mixture at 90 to 95 degrees C for 10 minutes yields oligonucleotides with 3 prime phosphoryl groups. The 3 prime cis diol modified oligos are synthesized by starting with 3 prime cis diol nucleotides (HO-NTPs). The cis dialdehyde modified oligonucleotides and cis diols may then be used for a variety of investigations such as studying the interaction of proteins with the 3 prime end of DNA, and possibly RNA. As an example, we demonstrate the efficacy of using an oligonucleotide modified with oGMP at the 3 prime end as an affinity label for TdT and identified a peptide fragment that has been shown to contain two of three aspartate residues found to be in the TdT active site.

Chemically induced proximity has not been systematically applied to control RNA fate. Here, a programmable platform was developed to identify RNA-binding proteins (RBPs) that can be recruited by small molecules to destabilize RNA. Using microtubule-associated protein Tau (MAPT) pre-mRNA as a model target, a heterobifunctional molecule was designed to bind both a ligandable structure in MAPT pre-mRNA and FKBP12F36V-tagged RBPs. Screening of a library of tagged RBPs identified several proteins that reduced MAPT RNA levels, including zinc finger protein 36 (ZFP36) and nanos C2HC-type zinc finger 3 (NANOS3). The approach was then extended from engineered proteins to an endogenous RBP. Using small molecule ligandability maps, a cysteine-reactive ligand for ZFP36 was identified. When this ligand was linked to the MAPT-binding small molecule, endogenous ZFP36 was recruited to MAPT mRNA, reducing its abundance in cells. Genetic and chemical controls demonstrated that activity was dependent on both RNA binding and ZFP36 recruitment, supporting an induced-proximity mechanism. These studies establish a general strategy for identifying new recruitable RBP effectors and should advance ribonuclease-targeting chimera (RiboTAC) technology by expanding the repertoire of effector proteins that can be harnessed for RNA degradation. More broadly, new effectors can be discovered through model reporter-based screens and translated to endogenous systems by mining known protein binders and ligandability maps, providing a systematic path to develop small molecules that control RNA stability, including RNAs targeted through structured regions of nuclear pre-mRNAs.

RNA surveillance and decay is carried out in part by helicase containing complexes that identify, capture and sometimes modify RNA before delivering it to the RNA exosome complex for processing or degradation. The MTREC core complex includes a Mtr4-like protein (Mtl1) helicase and Red1 that works with other cofactors and the RNA exosome in Schizosaccharomyces pombe to degrade nuclear transcripts in processes that can result in formation of facultative heterochromatin. The activities of Mtl1 remain uncharacterized as do contributions of Red1 to Mtl1 within MTREC. Here, we reconstitute the MTREC core complex, resolve structures by cryo-electron microscopy, and compare MTREC activities to S. pombe Mtr4 and Mtl1. We show that Mtl1 is more active relative to MTREC and Mtr4, that MTREC binds RNA better than Mtl1, and that Red1 includes an autoinhibitory coiled-coil domain that dimerizes MTREC and contacts the Mtl1 RecA domains to disrupt its ATPase active site. Together, these data suggest that Red1 may endow MTREC to bind RNA while slowing translocation so that it remains associated with RNA long enough to chaperone it to the RNA exosome for processing or degradation.

This study proposes a novel therapeutic strategy to suppress cancer growth by modulating the MDM2-p53 axis via Alternative Polyadenylation (APA). MDM2 normally promotes tumorigenesis by ubiquitinating and degrading the tumor suppressor p53. In cancer cells, preferential use of proximal polyadenylation signals (PAS) results in shortened 3'UTRs, allowing oncogenic transcripts like MDM2 to evade nuclear sequestration mediated by Inverted Alu (IRAlu) double-stranded RNA structures. We hypothesized that forcing distal PAS usage would elongate the MDM2 mRNA, promoting its nuclear retention and reducing protein translation, thereby restoring p53 activity. Using CRISPR-Cas9, we targeted and deleted the most frequent proximal PAS in the MDM2 3'UTR of A549 cells. Successful genome editing was confirmed via PCR. As expected, Western blot analysis showed a significant reduction in MDM2 expression in PAS-edited cells. However, experimental outcomes contradicted our initial hypothesis: edited cells exhibited higher viability under doxorubicin treatment compared to wild-type cells. Furthermore, despite decreased MDM2 levels, a concurrent reduction in phosphorylated p53 (p-p53) was observed. These unexpected results suggest that MDM2 3'UTR elongation may trigger a non-canonical regulatory mechanism that bypasses the traditional MDM2-p53 interaction. This study highlights the complexity of post-transcriptional regulation and suggests that APA-mediated gene modulation can induce unforeseen compensatory survival pathways in cancer cells, necessitating further investigation into the broader functional landscape of elongated 3'UTRs.

The QxxR motif is evolutionarily conserved within DEAD-box RNA helicases, including Drosophila Me31B and human DDX6, which post-transcriptionally regulate gene expression during animal development. A pathogenic H372R substitution (QxHR to QxRR) in the QxxR motif of human DDX6 has been associated with various developmental defects, but how this motif contributes to DDX6-family protein function remains unclear. Here, we used Drosophila Me31B as an in vivo model to investigate the QxxR motifs developmental role. We generated a Drosophila strain carrying the corresponding H333R missense mutation in Me31B and characterized its effects on female fertility, embryonic viability, germline development, and Me31B-associated molecular pathways. The me31BH333R mutation reduced female fertility in a gene dose-dependent manner, with homozygous mutant females being sterile. Embryos from the mutant females also exhibited primordial germ cell defects. Despite these developmental phenotypes, the me31BH333R mutation did not significantly alter Me31B protein abundance, global ovarian transcriptome or proteome profiles, or representative germ plasm mRNA and protein localization. In contrast, bait-normalized IP-MS analysis revealed altered enrichment of selected Me31B-associated proteins, including increased association of known Me31B interactors Trailer hitch (Tral) and Ypsilon Schachtel (Yps). These findings establish Me31BH333R as an in vivo model for investigating the conserved QxxR motif and suggest that disruption of this motif compromises development not through broad changes in gene expression, but potentially through altered composition or regulation of Me31B-containing ribonucleoprotein complexes.

A crucial aspect of the bioinformatics workflow in small RNA-sequencing is the alignment of reads to a database of reference ncRNAs. Alignment algorithms such as Bowtie, Burrows-Wheeler Aligner (BWA), and Spliced Transcripts Alignment to a Reference (STAR) - which are designed for aligning reads to a reference genome - are typically used. Aligning short RNA-sequenced reads to a database of non-coding RNAs (ncRNAs) is fundamentally a different task than aligning longer reads to a genome due to ncRNAs (i) having roughly the same number of nucleotides as the reads being aligned and (ii) being subsequences of other ncRNAs. To account for these differences, we developed the novel alignment algorithm LevenMap. Of all reads which exactly matched a reference ncRNA in a publicly available dataset, LevenMap aligned 100.0% of them to their respective ncRNA while all other aligners mapped less than 40% of these reads to their corresponding ncRNA. Furthermore, the mean ratio (length of read) / (length of corresponding reference ncRNA) of all aligned reads was 1.0 and 0.998 for LevenMap with at most zero and one mismatch(es) allowed, respectively; this ratio was no more than 0.51 for all other aligners. Overall, LevenMap is designed to account for the nuances of aligning small RNA-sequencing data to a database of reference ncRNAs and yields more biologically relevant counts compared to traditional aligners in this context. LevenMap is free and publicly available on GitHub: https://github.com/hdlugas/LevenMap.

The 5' untranslated region (5' UTR) shapes translation initiation, so its design is central to mRNA therapeutics and to improving protein-production cell lines. Deep-learning models that predict translation efficiency, measured as mean ribosome load (MRL), from the 5' UTR sequence have been combined with genetic algorithms (GAs) for sequence optimization. However, optimizing against a model trained on offline data risks reward hacking that exploits the models estimation error outside the training distribution, yielding sequences that score highly in prediction yet fail to perform in the wet lab. Yet for 5' UTR design, few studies have systematically examined which region should be treated as untrustworthy (the definition of out-of-distribution, OOD) or which constraints keep the search away from it. We present a constrained optimization that keeps candidates within a trust region where the predictors validated accuracy holds; here "reliable" denotes keeping candidates within the training distribution over which prediction has been validated, not a guarantee of measured performance. As the OOD score, we compare the k-nearest-neighbor (KNN) distance in the predictors embedding space against a pseudo-perplexity (PPPL) from the encoder and LM head, and show that for nucleotide sequences--whose vocabulary is small--PPPL fails to separate in- vs out-of-distribution, whereas the KNN distance is an effective OOD score that can define a trust region even from unlabeled native UTR sequences. Using the KNN distance as a hard GA constraint keeps all candidates inside the trust region while maintaining predicted MRL: under unconstrained optimization most final-generation candidates (72-96% across seeds) left the trust region (self-KNN p95), whereas the hard constraint holds predicted MRL at the unconstrained level and yields about 4.3x more selectable low-risk candidates than post-hoc filtering of the unconstrained output. Comparing an output extrapolation guard, reference-sequence similarity and structural accessibility (RNAplfold), we find that the guard and the similarity constraint also suppress OOD as a side effect, whereas making accessibility a secondary objective broadens the search without suppressing OOD.

Bacterial small RNAs (sRNAs) regulate gene expression by base pairing with target mRNAs, yet transcriptome-wide interactome mapping has shown that many sRNA-mRNA interactions detected in vivo have modest or no regulatory effect using orthogonal reporter assays. The features that determine functional outcome remain poorly defined. Here, we integrated Hfq-CLASH interactome mapping with matched transcriptomic and proteomic profiling in Escherichia coli and developed an interpretable machine-learning framework to identify the determinants that distinguish functional from non-functional interactions. Using sequence, structural, thermodynamic, duplex and protein-occupancy features, transcriptomic and proteomic responses were predicted with above-chance performance, achieving AUCs of 0.78 and 0.74, respectively. Feature attribution revealed that physical pairing alone is insufficient for regulation; instead, regulatory outcome is shaped by a coordinated interplay between RNA secondary structure, thermodynamic accessibility and local protein-binding context. Target-side Hfq occupancy emerged as a positive predictor of functional regulation, whereas AR2-domain occupancy on the sRNA was associated with non-responsive interactions, suggesting that distinct ribonucleoprotein states may separate productive regulation from non-productive binding. These findings indicate that the regulatory fate of an sRNA-mRNA interaction is an emergent property of its biophysical context and protein-binding environment, rather than a direct consequence of physical pairing alone.

Haploinsufficiency in SYNGAP1 causes a severe neurodevelopmental syndrome. SYNGAP1 protein is mainly detected in neuronal synapses. However, SYNGAP1 RNA is more widely expressed and strongly regulated via alternative splicing: alternative 3' splice site (A3SS) inclusion leads to non-productive transcripts that are degraded through nonsense-mediated decay. Recently, splice-switching oligonucleotides (SSOs) that redirect SYNGAP1 splicing to increase SYNGAP1 protein levels were developed. However, we hypothesized that during neuronal maturation, non-productive splicing may decrease to enhance functional transcripts in mature neurons. This would reduce the abundance of the SSO target transcript, limiting the potential for SSO treatment to increase neuronal SYNGAP1 expression. Using neural differentiation of human induced pluripotent stem cells, we show that the A3SS transcript is abundant in neural progenitors, astrocytes, microglia and immature neurons, with minimal presence in mature neurons. These data imply that SSOs targeting A3SS might lack therapeutic efficacy to rescue the neuronal phenotypes associated with SYNGAP1 haploinsufficiency.

MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression primarily in the cytoplasm. However, emerging evidence highlights their additional roles in the nucleus. In particular, spliceosomal miRNAs have been implicated in novel regulatory functions, including the modulation of gene expression. Here, we investigate the nuclear role of spliceosomal miR-99b in breast cancer cells, focusing on its interaction with the long non-coding RNA (lncRNA) SPACA6-AS1. Using non-tumorigenic (MCF-10A) and breast cancer cell lines (MCF-7 and MDA-MB-231), we demonstrate that spliceosomal miR-99b expression increases with malignancy and correlates with elevated SPACA6-AS1 pre-mRNA levels. Notably, miR-99b exhibits full complementarity to the 5-prime splice junction of SPACA6-AS1, suggesting a direct role in splicing regulation. Functional assays reveal that inhibition of miR-99b reduces SPACA6-AS1 pre-mRNA levels, whereas its overexpression enhances pre-mRNA accumulation, indicating that miR-99b promotes the formation or stabilization of the unspliced transcript. Furthermore, increased miR-99b expression is associated with altered ratios of SPACA6 isoforms, supporting a broader role in RNA-level regulation of gene expression. Phenotypically, miR-99b enhances breast cancer cell migration and is required for efficient invasion, particularly in highly aggressive cancerous cells. Our findings uncover a novel nuclear function of miR-99b in modulating lncRNA splicing and gene expression. This spliceosomal miR-99b-SPACA6-AS1 axis represents a previously unrecognized regulatory pathway that contributes to breast cancer progression and may provide a potential target for diagnostic and therapeutic strategies.

Transfer RNA (tRNA) modifications are essential for translational accuracy and cellular adaptation to environmental changes. While short-term modification dynamics are well documented, the impact of prolonged stress exposure on the global tRNA landscape remains largely unexplored. This study provides the first systematic profiling of tRNA modifications in Saccharomyces cerevisiae following long-term exposure to distinct stress types: heat, suboptimal pH, oxidative stress (paraquat and diamide), osmotic stress (NaCl and KCl), and genotoxic stress (MMS). To this end, we used a broad-range UPLC-MS protocol to quantify global changes in tRNA modification and identify stress-specific signatures. The results revealed that long-term stress triggers a global reprogramming of the tRNA epitranscriptome in a stress-specific and time-dependent manner. Importantly, while our findings confirm the previously reported temperature-sensitivity of wobble uridine thiolation, we also identified a complete or partial loss of 5-methoxycarbonylmethyl-2-thiouridine (mcm5s2U34) modification upon exposure to paraquat and pH stress. Furthermore, this loss of thiolation is accompanied by an accumulation of the non-thiolated 5-methoxycarbonylmethyl (mcm5U) precursor, suggesting a stress-dependent impairment of the thiolation pathway. Next, we applied MarathonRT-based tRNA sequencing and showed that these modification changes occur independently of tRNA isoacceptor abundance. To further integrate these results, we devised a modification deviation (MDm) index, which indicates that the observed reprogramming is primarily linked to events that are independent from changes in tRNA abundance. Together, this study provides a comprehensive atlas of tRNA modification dynamics under prolonged stress, addressing a critical gap in our understanding of RNA-based translational control and establishes the MDm index as a robust quantitative framework to decouple the influence of tRNA abundance from global modification signals, providing a necessary metric for the field to interpret epitranscriptomic reprogramming. TABLE OF CONTENTS GRAPHIC O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=108 SRC="FIGDIR/small/745200v1_ufig1.gif" ALT="Figure 1"> View larger version (41K): org.highwire.dtl.DTLVardef@2cf571org.highwire.dtl.DTLVardef@1a5f3c2org.highwire.dtl.DTLVardef@266d61org.highwire.dtl.DTLVardef@d7214f_HPS_FORMAT_FIGEXP M_FIG C_FIG