Manuel Ares

Splicing deficiency may represent a critical yet underexplored form of genomic erosion in non-recombining regions. Across four phytoplankton species diverged ∼333–639 million years ago, genes within U (female) and V (male) “UV” mating-type regions—non-recombining chromosomal regions that determine mating compatibility—show strikingly elevated intron retention relative to genes in other genomic regions. Long-read data reveal abundant aberrant, likely non-functional mRNA isoforms despite preserved coding potential. This preservation suggests that splicing defects arose early in UV evolution and have persisted over deep time. We propose that these defects arise from evolutionary changes in sequence composition and chromatin organization that accompany recombination suppression, such as reduced GC content, altered nucleosome occupancy, and disrupted methylation, that collectively compromise splicing fidelity. Unlike sex chromosomes, which often degenerate through gene loss, splicing-deficient UV regions in green algae retain hundreds of genes, indicating that transcript-level dysfunction provides an alternative route to functional decay. Our results identify chromatin-mediated splicing deficiency as a novel axis of genomic erosion and position algal UV systems as models for studying how recombination suppression reshapes RNA processing fidelity in essential, non-recombining genomes.

Pre-mRNA secondary structures are hypothesized to regulate RNA processing pathways, but such structures have been difficult to visualize in vivo. Here, we characterize Saccharomyces cerevisiae pre-mRNA structures through transcriptome-wide dimethyl sulfate probing, enriching for low-abundance pre-mRNA through splicing inhibition. We cross-validate structures found from phylogenetic and mutational studies and identify structures within the majority of measured introns (79 of 88). We find widespread formation of ‘zipper stems’ between the 5′ splice site and branch point, ‘downstream stems’ between the branch point and the 3′ splice site, and previously uncharacterized long stems that distinguish pre-mRNA from spliced mRNA. Multi-dimensional chemical mapping reveals intron structures that independently form in vitro without the presence of binding partners, and structure ensemble prediction suggests that such structures appear in introns across the Saccharomyces genus. We further develop a high-throughput functional assay to characterize variants of RNA structure (VARS-seq), applying it to 135 sets of stems across 7 introns, identifying structured elements that alter retained intron levels at a distance from canonical splice sites. This transcriptome-wide inference of intron RNA structures introduces alternative paradigms and model systems for understanding how pre-mRNA folding influences gene expression.

A monthly journal publishing high-quality, peer-reviewed research on all topics related to RNA and its metabolism in all organisms

A previous study using Thermostable Group II Intron Reverse Transcriptase sequencing (TGIRT-seq) found human plasma contains short (≤300 nt) structured full-length excised linear intron (FLEXI) RNAs with potential to serve as blood-based biomarkers. Here, TGIRT-seq identified >9,000 different FLEXI RNAs in human cell lines, including relatively abundant FLEXIs with cell-type-specific expression patterns. Analysis of public CLIP-seq datasets identified 126 RNA-binding proteins (RBPs) that have binding sites within the region corresponding to the FLEXI or overlapping FLEXI splice sites in pre-mRNAs, including 53 RBPs with binding sites for ≥30 different FLEXIs. These included splicing factors, transcription factors, a chromatin remodeling protein, cellular growth regulators, and proteins with cytoplasmic functions. Analysis of ENCODE datasets identified subsets of these RBPs whose knockdown impacted FLEXI host gene mRNA levels or proximate alternative splicing, indicating functional interactions. Hierarchical clustering identified six subsets of RBPs whose FLEXI binding sites were co-enriched in six subsets of functionally related host genes: AGO1-4 and DICER, including but not limited to agotrons or mirtron pre-miRNAs; DKC1, NOLC1, SMNDC1, and AATF (Apoptosis Antagonizing Transcription Factor), including but not limited to snoRNA-encoding FLEXIs; two subsets of alternative splicing factors; and two subsets that included RBPs with cytoplasmic functions (e.g., LARP4, PABPC4, METAP2, and ZNF622) together with regulatory proteins. Cell fractionation experiments showed cytoplasmic enrichment of FLEXI RNAs with binding sites for RBPs with cytoplasmic functions. The subsets of host genes encoding FLEXIs with binding sites for different subsets of RBPs were co-enriched with non-FLEXI other short and long introns with binding sites for the same RBPs, suggesting overarching mechanisms for coordinately regulating expression of functionally related genes. Our findings identify FLEXIs as a previously unrecognized large class of cellular RNAs and provide a comprehensive roadmap for further analyzing their biological functions and the relationship of their RBPs to cellular regulatory mechanisms.

Rare, full-length circular intron RNAs distinct from lariats have been reported in several species, but their biogenesis is not understood. We envisioned and tested a hypothesis for their formation usingSaccharomyces cerevisiae, documenting full-length and novel processed circular RNAs from multiple introns. Evidence implicates a previously undescribed catalytic activity of the intron lariat spliceosome (ILS) in which the 3′-OH of the lariat tail (with optional trimming and adenylation by the nuclear 3′ processing machinery) attacks the branch, joining the intron 3′ end to the 5′ splice site in a 3′–5′ linked circle. Human U2 and U12 spliceosomes produce analogous full-length and processed circles. Postsplicing catalytic activity of the spliceosome may promote intron transposition during eukaryotic genome evolution.

Rare, full length circular intron RNAs distinct from lariats have been reported in several species, but their biogenesis is not understood. We envision and test a hypothesis for their formation using Saccharomyces cerevisiae , documenting full length and novel processed circular RNAs from multiple introns. Evidence implicates a previously undescribed catalytic activity of the intron-lariat spliceosome (ILS) in which the 3’-OH of the lariat tail (with optional trimming and adenylation by the nuclear 3’ processing machinery) attacks the branch, joining the intron 3’ end to the 5’ splice site in a 3’-5’ linked circle. Human U2 and U12 spliceosomes produce analogous full length and processed circles. Post-splicing catalytic activity of the spliceosome may promote intron transposition during eukaryotic genome evolution.

Intron branchpoint (BP) recognition by the U2 snRNP is a critical step of splicing, vulnerable to recurrent cancer mutations and bacterial natural product inhibitors. The BP binds a conserved pocket in the SF3B1 (human) or Hsh155 (yeast) U2 snRNP protein. Amino acids that line this pocket affect the binding of splicing inhibitors like Pladienolide-B (Plad-B), such that organisms differ in their sensitivity. To study the mechanism of splicing inhibitor action in a simplified system, we modified the naturally Plad-B resistant yeastSaccharomyces cerevisiaeby changing 14 amino acids in the Hsh155 BP pocket to those from human. This humanized yeast grows normally, and splicing is largely unaffected by the mutation. Splicing is inhibited within minutes after the addition of Plad-B, and different introns appear inhibited to different extents. Intron-specific inhibition differences are also observed during cotranscriptional splicing in Plad-B using single-molecule intron tracking to minimize gene-specific transcription and decay rates that cloud estimates of inhibition by standard RNA-seq. Comparison of Plad-B intron sensitivities to those of the structurally distinct inhibitor Thailanstatin-A reveals intron-specific differences in sensitivity to different compounds. This work exposes a complex relationship between the binding of different members of this class of inhibitors to the spliceosome and intron-specific rates of BP recognition and catalysis. Introns with variant BP sequences seem particularly sensitive, echoing observations from mammalian cells, where monitoring individual introns is complicated by multi-intron gene architecture and alternative splicing. The compact yeast system may hasten the characterization of splicing inhibitors, accelerating improvements in selectivity and therapeutic efficacy.

Intron branch point (BP) recognition by the U2 snRNP is a critical step of splicing, vulnerable to recurrent cancer mutations and bacterial natural product inhibitors. The BP binds a conserved pocket in the SF3B1 (human) or Hsh155 (yeast) U2 snRNP protein. Amino acids that line this pocket affect binding of splicing inhibitors like Pladienolide-B (Plad-B), such that organisms differ in their sensitivity. To study the mechanism of splicing inhibitor action in a simplified system, we modified the naturally Plad-B resistant yeast Saccharomyces cerevisiae by changing 14 amino acids in the Hsh155 BP pocket to those from human. This humanized yeast grows normally, and splicing is largely unaffected by the mutation. Splicing is inhibited within minutes after addition of Plad-B, and different introns appear inhibited to different extents. Intron-specific inhibition differences are also observed during co-transcriptional splicing in Plad-B using single-molecule intron tracking (SMIT) to minimize gene-specific transcription and decay rates that cloud estimates of inhibition by standard RNA-seq. Comparison of Plad-B intron sensitivities to those of the structurally distinct inhibitor Thailanstatin-A reveals intron-specific differences in sensitivity to different compounds. This work exposes a complex relationship between binding of different members of this class of inhibitors to the spliceosome and intron-specific rates of BP recognition and catalysis. Introns with variant BP sequences seem particularly sensitive, echoing observations from mammalian cells, where monitoring individual introns is complicated by multi-intron gene architecture and alternative splicing. The compact yeast system may hasten characterization of splicing inhibitors, accelerating improvements in selectivity and therapeutic efficacy.

There is massive variation in intron numbers across eukaryotic genomes, yet the major drivers of intron content during evolution remain elusive. Rapid intron loss and gain in some lineages contrast with long-term evolutionary stasis in others. Episodic intron gain could be explained by recently discovered specialized transposons called Introners, but so far Introners are only known from a handful of species. Here, we performed a systematic search across 3,325 eukaryotic genomes and identified 27,563 Introner-derived introns in 175 genomes (5.2%). Species with Introners span remarkable phylogenetic diversity, from animals to basal protists, representing lineages whose last common ancestor dates to over 1.7 billion years ago. Aquatic organisms were 6.5 times more likely to contain Introners than terrestrial organisms. Introners exhibit mechanistic diversity but most are consistent with DNA transposition, indicating that Introners have evolved convergently hundreds of times from nonautonomous transposable elements. Transposable elements and aquatic taxa are associated with high rates of horizontal gene transfer, suggesting that this combination of factors may explain the punctuated and biased diversity of species containing Introners. More generally, our data suggest that Introners may explain the episodic nature of intron gain across the eukaryotic tree of life. These results illuminate the major source of ongoing intron creation in eukaryotic genomes.

Pre-mRNA secondary structures are hypothesized to play widespread roles in regulating RNA processing pathways, but these structures have been difficult to visualize in vivo . Here, we characterize S. cerevisiae pre-mRNA structures through transcriptome-wide dimethyl sulfate (DMS) probing, enriching for low-abundance pre-mRNA through splicing inhibition. We cross-validate structures found from phylogenetic and mutational studies and identify new structures within the majority of probed introns (102 of 161). We find widespread formation of “zipper stems” between the 5’ splice site and branch point, “downstream stems” between the branch point and the 3’ splice site, and previously uncharacterized long stems that distinguish pre-mRNA from spliced mRNA. Multi-dimensional chemical mapping reveals examples where intron structures can form in vitro without the presence of binding partners, and structure ensemble prediction suggests that such structures appear in introns across the Saccharomyces genus. We develop a high-throughput functional assay to characterize variants of RNA structure (VARS-seq) and we apply the method on 135 sets of stems across 7 introns, identifying structured elements that alter retained intron levels at a distance from canonical splice sites. This transcriptome-wide inference of intron RNA structures suggests new ideas and model systems for understanding how pre-mRNA folding influences gene expression.