Comprehensive splice-site analysis using comparative genomics

Splicing analysis of STAT3 tandem donor suggests non-canonical binding registers for U1 and U6 snRNAs

Tandem donor splice sites (5'ss) are unique regions with at least two GU dinucleotides serving as splicing cleavage sites. The Δ3 tandem 5'ss are a specific subclass of 5'ss separated by 3 nucleotides which can affect protein function by inserting/deleting a single amino acid. One 5'ss is typically preferred, yet factors governing particular 5'ss choice are not fully understood. A highly conserved exon 21 of the STAT3 gene was chosen as a model to study Δ3 tandem 5'ss splicing mechanisms. Based on multiple lines of experimental evidence, endogenous U1 snRNA most likely binds only to the upstream 5'ss. However, the downstream 5'ss is used preferentially, and the splice site choice is not dependent on the exact U1 snRNA binding position. Downstream 5'ss usage was sensitive to exact nucleotide composition and dependent on the presence of downstream regulatory region. The downstream 5'ss usage could be best explained by two novel interactions with endogenous U6 snRNA. U6 snRNA enables the downstream 5'ss usage in STAT3 exon 21 by two mechanisms: (i) binding in a novel non-canonical register and (ii) establishing extended Watson-Crick base pairing with the downstream regulatory region. This study suggests that U6:5'ss interaction is more flexible than previously thought.

U6 snRNA m6A modification is required for accurate and efficient splicing of C. elegans and human pre-mRNAs

pre-mRNA splicing is a critical feature of eukaryotic gene expression. Both cis- and trans-splicing rely on accurately recognising splice site sequences by spliceosomal U snRNAs and associated proteins. Spliceosomal snRNAs carry multiple RNA modifications with the potential to affect different stages of pre-mRNA splicing. Here, we show that the conserved U6 snRNA m6A methyltransferase METT-10 is required for accurate and efficient cis- and trans-splicing of C. elegans pre-mRNAs. The absence of METT-10 in C. elegans and METTL16 in humans primarily leads to alternative splicing at 5' splice sites with an adenosine at +4 position. In addition, METT-10 is required for splicing of weak 3' cis- and trans-splice sites. We identified a significant overlap between METT-10 and the conserved splicing factor SNRNP27K in regulating 5' splice sites with +4A. Finally, we show that editing endogenous 5' splice site +4A positions to +4U restores splicing to wild-type positions in a mett-10 mutant background, supporting a direct role for U6 snRNA m6A modification in 5' splice site recognition. We conclude that the U6 snRNA m6A modification is important for accurate and efficient pre-mRNA splicing.

Alternative splicing in prostate cancer progression and therapeutic resistance

Prostate cancer (CaP) remains the second leading cause of cancer deaths in western men. CaP mortality results from diverse molecular mechanisms that mediate resistance to the standard of care treatments for metastatic disease. Recently, alternative splicing has been recognized as a hallmark of CaP aggressiveness. Alternative splicing events cause treatment resistance and aggressive CaP behavior and are determinants of the emergence of the two major types of late-stage treatment-resistant CaP, namely castration-resistant CaP (CRPC) and neuroendocrine CaP (NEPC). Here, we review recent multi-omics data that are uncovering the complicated landscape of alternative splicing events during CaP progression and the impact that different gene transcript isoforms can have on CaP cell biology and behavior. We discuss renewed insights in the molecular machinery by which alternative splicing occurs and contributes to the failure of systemic CaP therapies. The potential for alternative splicing events to serve as diagnostic markers and/or therapeutic targets is explored. We conclude by considering current challenges and promises associated with splicing-modulating therapies, and their potential for clinical translation into CaP patient care.

Comprehensive deciphering the alternative splicing patterns involved in leaf morphogenesis of Liriodendron chinense

Alternative splicing (AS), a pivotal post-transcriptional regulatory mechanism, profoundly amplifies diversity and complexity of transcriptome and proteome. Liriodendron chinense (Hemsl.) Sarg., an excellent ornamental tree species renowned for its distinctive leaf shape, which resembles the mandarin jacket. Despite the documented potential genes related to leaf development of L. chinense, the underlying post-transcriptional regulatory mechanisms remain veiled. Here, we conducted a comprehensive analysis of the transcriptome to clarify the genome-wide landscape of the AS pattern and the spectrum of spliced isoforms during leaf developmental stages in L. chinense. Our investigation unveiled 50,259 AS events, involving 10,685 genes (32.9%), with intron retention as the most prevalent events. Notably, the initial stage of leaf development witnessed the detection of 804 differentially AS events affiliated with 548 genes. Although both differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) were enriched into morphogenetic related pathways during the transition from fishhook (P2) to lobed (P7) leaves, there was only a modest degree of overlap between DASGs and DEGs. Furthermore, we conducted a comprehensively AS analysis on homologous genes involved in leaf morphogenesis, and most of which are subject to post-transcriptional regulation of AS. Among them, the AINTEGUMENTA-LIKE transcript factor LcAIL5 was characterization in detailed, which experiences skipping exon (SE), and two transcripts displayed disparate expression patterns across multiple stages. Overall, these findings yield a comprehensive understanding of leaf development regulation via AS, offering a novel perspective for further deciphering the mechanism of plant leaf morphogenesis.

Genomic survey and expression analysis of cellulose synthase superfamily and COBRA-like gene family in Zanthoxylum bungeanum stipule thorns

The Cellulose Synthase gene (CS) superfamily and COBRA-like (COBL) gene family are essential for synthesizing cellulose and hemicellulose, which play a crucial role in cell wall biosynthesis and the hardening of plant tissues. Our study identified 126 ZbCS and 31 ZbCOBL genes from the Zanthoxylum bungeanum (Zb) genome. Phylogenetic analysis and conservative domain analysis unfolded that ZbCS and ZbCOBL genes were divided into seven and two subfamilies, respectively. Gene duplication data suggested that more than 75% of these genes had tandem and fragment duplications. Codon usage patterns analysis indicated that the ZbCS and ZbCOBL genes prefer ending with A/T base, with weak codon preference. Furthermore, seven key ZbCS and five key ZbCOBL genes were identified based on the content of cellulose and hemicellulose and the expression characteristics of ZbCS and ZbCOBL genes in various stages of stipule thorns. Altogether, these results improve the understanding of CS and COBL genes and provide valuable reference data for cultivating Zb with soft thorns.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01432-x.

Exitrons: offering new roles to retained introns-the novel regulators of protein diversity and utility

Exitrons are exonic introns. This subclass of intron retention alternative splicing does not contain a Pre-Terminating stop Codon. Therefore, when retained, they are always a part of a protein. Intron retention is a frequent phenomenon predominantly found in plants, which results in either the degradation of the transcripts or can serve as a stable intermediate to be processed upon induction by specific signals or the cell status. Interestingly, exitrons have coding ability and may confer additional attributes to the proteins that retain them. Therefore, exitron-containing and exitron-spliced isoforms will be a driving force for creating protein diversity in the proteome of an organism. This review establishes a basic understanding of exitron, discussing its genesis, key features, identification methods and functions. We also try to depict its other potential roles. The present review also aims to provide a fundamental background to those who found such exitronic sequences in their gene(s) and to speculate the future course of studies.

Where the minor things are: a pan-eukaryotic survey suggests neutral processes may explain much of minor intron evolution

Spliceosomal introns are gene segments removed from RNA transcripts by ribonucleoprotein machineries called spliceosomes. In some eukaryotes a second 'minor' spliceosome is responsible for processing a tiny minority of introns. Despite its seemingly modest role, minor splicing has persisted for roughly 1.5 billion years of eukaryotic evolution. Identifying minor introns in over 3000 eukaryotic genomes, we report diverse evolutionary histories including surprisingly high numbers in some fungi and green algae, repeated loss, as well as general biases in their positional and genic distributions. We estimate that ancestral minor intron densities were comparable to those of vertebrates, suggesting a trend of long-term stasis. Finally, three findings suggest a major role for neutral processes in minor intron evolution. First, highly similar patterns of minor and major intron evolution contrast with both functionalist and deleterious model predictions. Second, observed functional biases among minor intron-containing genes are largely explained by these genes' greater ages. Third, no association of intron splicing with cell proliferation in a minor intron-rich fungus suggests that regulatory roles are lineage-specific and thus cannot offer a general explanation for minor splicing's persistence. These data constitute the most comprehensive view of minor introns and their evolutionary history to date, and provide a foundation for future studies of these remarkable genetic elements.

Generation of Mutants from the 57B Region of Drosophila melanogaster

The 57B region of Drosophila melanogaster includes a cluster of the three homeobox genes orthopedia (otp), Drosophila Retinal homeobox (DRx), and homeobrain (hbn). In an attempt to isolate mutants for these genes, we performed an EMS mutagenesis and isolated lethal mutants from the 57B region, among them mutants for otp, DRx, and hbn. With the help of two newly generated deletions from the 57B region, we mapped additional mutants to specific chromosomal intervals and identified several of these mutants from the 57B region molecularly. In addition, we generated mutants for CG15651 and RIC-3 by gene targeting and mutants for the genes CG9344, CG15649, CG15650, and ND-B14.7 using the CRISPR/Cas9 system. We determined the lethality period during development for most isolated mutants. In total, we analysed alleles from nine different genes from the 57B region of Drosophila, which could now be used to further explore the functions of the corresponding genes in the future.

Heterogeneous splicing patterns resulting from KIF5A variants associated with amyotrophic lateral sclerosis

Single-nucleotide variants (SNVs) in the gene encoding Kinesin Family Member 5A (KIF5A), a neuronal motor protein involved in anterograde transport along microtubules, have been associated with amyotrophic lateral sclerosis (ALS). ALS is a rapidly progressive and fatal neurodegenerative disease that primarily affects the motor neurons. Numerous ALS-associated KIF5A SNVs are clustered near the splice-site junctions of the penultimate exon 27 and are predicted to alter the carboxy-terminal (C-term) cargo-binding domain of KIF5A. Mis-splicing of exon 27, resulting in exon exclusion, is proposed to be the mechanism by which these SNVs cause ALS. Whether all SNVs proximal to exon 27 result in exon exclusion is unclear. To address this question, we designed an in vitro minigene splicing assay in human embryonic kidney 293 cells, which revealed heterogeneous site-specific effects on splicing: only 5' splice-site (5'ss) SNVs resulted in exon skipping. We also quantified splicing in select clustered, regularly interspaced, short palindromic repeats-edited human stem cells, differentiated to motor neurons, and in neuronal tissues from a 5'ss SNV knock-in mouse, which showed the same result. Moreover, the survival of representative 3' splice site, 5'ss, and truncated C-term variant KIF5A (v-KIF5A) motor neurons was severely reduced compared with wild-type motor neurons, and overt morphological changes were apparent. While the total KIF5A mRNA levels were comparable across the cell lines, the total KIF5A protein levels were decreased for v-KIF5A lines, suggesting an impairment of protein synthesis or stability. Thus, despite the heterogeneous effect on ribonucleic acid splicing, KIF5A SNVs similarly reduce the availability of the KIF5A protein, leading to axonal transport defects and motor neuron pathology.

Manual annotation of Drosophila genes: a Genomics Education Partnership protocol

Annotating the genomes of multiple species allows us to analyze the evolution of their genes. While many eukaryotic genome assemblies already include computational gene predictions, these predictions can benefit from review and refinement through manual gene annotation. The Genomics Education Partnership (GEP; https://thegep.org/) developed a structural annotation protocol for protein-coding genes that enables undergraduate student and faculty researchers to create high-quality gene annotations that can be utilized in subsequent scientific investigations. For example, this protocol has been utilized by the GEP faculty to engage undergraduate students in the comparative annotation of genes involved in the insulin signaling pathway in 27 Drosophila species, using D. melanogaster as the reference genome. Students construct gene models using multiple lines of computational and empirical evidence including expression data (e.g., RNA-Seq), sequence similarity (e.g., BLAST and multiple sequence alignment), and computational gene predictions. Quality control measures require each gene be annotated by at least two students working independently, followed by reconciliation of the submitted gene models by a more experienced student. This article provides an overview of the annotation protocol and describes how discrepancies in student submitted gene models are resolved to produce a final, high-quality gene set suitable for subsequent analyses. The protocol can be adapted to other scientific questions (e.g., expansion of the Drosophila Muller F element) and species (e.g., parasitoid wasps) to provide additional opportunities for undergraduate students to participate in genomics research. These student annotation efforts can substantially improve the quality of gene annotations in publicly available genomic databases.

U2 and U4 snDNA Comparative Chromosomal Mapping in the Neotropical Fish Genera Apareiodon and Parodon (Characiformes: Parodontidae)

Small nuclear DNA (snDNA) are valuable cytogenetic markers for comparative studies in chromosome evolution because different distribution patterns were found among species. Parodontidae, a Neotropical fish family, is known to have female heterogametic sex chromosome systems in some species. The U2 and U4 snDNA sites have been found to be involved in Z and W chromosome differentiation in Apareiodon sp., Apareiodon affinis, and Parodon hilarii. However, few studies have evaluated snDNA sites as propulsors of chromosome diversification among closely related fish species. In this study, we investigated the distribution of U2 and U4 snDNA clusters in the chromosomes of 10 populations/species belonging to Apareiodon and Parodon, aiming to identify chromosomal homeologies or diversification. In situ localization data revealed a submetacentric pair carrying the U2 snDNA site among the populations/species analyzed. Furthermore, all studied species demonstrated homeology in the location of U4 snDNA cluster in the proximal region of metacentric pair 1, besides an additional signal showing up with a divergence in Apareiodon. Comparative chromosomal mapping of U4 snDNA also helped to reinforce the proposal of the ZZ/ZW1W2 sex chromosome system origin in an A. affinis population. According to cytogenetic data, the study corroborates the diversification in Parodontidae paired species with uncertain taxonomy.

Sequence variations affect the 5' splice site selection of plant introns

Introns are noncoding sequences spliced out of pre-mRNAs by the spliceosome to produce mature mRNAs. The 5' ends of introns mostly begin with GU and have a conserved sequence motif of AG/GUAAGU that could base-pair with the core sequence of U1 snRNA of the spliceosome. Intriguingly, ∼ 1% of introns in various eukaryotic species begin with GC. This occurrence could cause misannotation of genes; however, the underlying splicing mechanism is unclear. We analyzed the sequences around the intron 5' splice site (ss) in Arabidopsis (Arabidopsis thaliana) and found sequences at the GC intron ss are much more stringent than those of GT introns. Mutational analysis at various positions of the intron 5' ss revealed that although mutations impair base pairing, different mutations at the same site can have different effects, suggesting that steric hindrance also affects splicing. Moreover, mutations of 5' ss often activate a hidden ss nearby. Our data suggest that the 5' ss is selected via a competition between the major ss and the nearby minor ss. This work not only provides insights into the splicing mechanism of intron 5' ss but also improves the accuracy of gene annotation and the study of the evolution of intron 5' ss.

Splicing defects in rare diseases: transcriptomics and machine learning strategies towards genetic diagnosis

Genomic variants affecting pre-messenger RNA splicing and its regulation are known to underlie many rare genetic diseases. However, common workflows for genetic diagnosis and clinical variant interpretation frequently overlook splice-altering variants. To better serve patient populations and advance biomedical knowledge, it has become increasingly important to develop and refine approaches for detecting and interpreting pathogenic splicing variants. In this review, we will summarize a few recent developments and challenges in using RNA sequencing technologies for rare disease investigation. Moreover, we will discuss how recent computational splicing prediction tools have emerged as complementary approaches for revealing disease-causing variants underlying splicing defects. We speculate that continuous improvements to sequencing technologies and predictive modeling will not only expand our understanding of splicing regulation but also bring us closer to filling the diagnostic gap for rare disease patients.

Introns: the "dark matter" of the eukaryotic genome

The emergence of introns was a significant evolutionary leap that is a major distinguishing feature between prokaryotic and eukaryotic genomes. While historically introns were regarded merely as the sequences that are removed to produce spliced transcripts encoding functional products, increasingly data suggests that introns play important roles in the regulation of gene expression. Here, we use an intron-centric lens to review the role of introns in eukaryotic gene expression. First, we focus on intron architecture and how it may influence mechanisms of splicing. Second, we focus on the implications of spliceosomal snRNAs and their variants on intron splicing. Finally, we discuss how the presence of introns and the need to splice them influences transcription regulation. Despite the abundance of introns in the eukaryotic genome and their emerging role regulating gene expression, a lot remains unexplored. Therefore, here we refer to introns as the "dark matter" of the eukaryotic genome and discuss some of the outstanding questions in the field.

Recent advances in the molecular understanding of medulloblastoma

Medulloblastoma is the most common pediatric malignant brain tumor composed of four molecular subgroups. Recent intensive genomics has greatly contributed to our understanding of medulloblastoma pathogenesis. Sequencing studies identified novel mutations involved in the cyclic AMP-dependent pathway or RNA processing in the Sonic Hedgehog (SHH) subgroup, and core-binding factor subunit alpha (CBFA) complex in the group 4 subgroup. Likewise, single-cell sequencing provided detailed insights into the cell of origin associated with brain development. In this review, we will summarize recent findings by sequencing analyses for medulloblastoma.

Pitfalls of using sequence databases for heterologous expression studies - a technical review

Synthesis of DNA fragments based on gene sequences that are available in public resources has become an efficient and affordable method that has gradually replaced traditional cloning efforts such as PCR cloning from cDNA. However, database entries based on genome sequencing results are prone to errors which can lead to false sequence information and, ultimately, errors in functional characterisation of proteins such as ion channels and transporters in heterologous expression systems. We have identified five common problems that repeatedly appear in public resources: (1) Not every gene has yet been annotated; (2) not all gene annotations are necessarily correct; (3) transcripts may contain automated corrections; (4) there are mismatches between gene, mRNA and protein sequences; and (5) splicing patterns often lack experimental validation. This technical review highlights and provides a strategy to bypass these issues in order to avoid critical mistakes that could impact future studies of any gene/protein of interest in heterologous expression systems.

ESPRESSO: Robust discovery and quantification of transcript isoforms from error-prone long-read RNA-seq data

Long-read RNA sequencing (RNA-seq) holds great potential for characterizing transcriptome variation and full-length transcript isoforms, but the relatively high error rate of current long-read sequencing platforms poses a major challenge. We present ESPRESSO, a computational tool for robust discovery and quantification of transcript isoforms from error-prone long reads. ESPRESSO jointly considers alignments of all long reads aligned to a gene and uses error profiles of individual reads to improve the identification of splice junctions and the discovery of their corresponding transcript isoforms. On both a synthetic spike-in RNA sample and human RNA samples, ESPRESSO outperforms multiple contemporary tools in not only transcript isoform discovery but also transcript isoform quantification. In total, we generated and analyzed ~1.1 billion nanopore RNA-seq reads covering 30 human tissue samples and three human cell lines. ESPRESSO and its companion dataset provide a useful resource for studying the RNA repertoire of eukaryotic transcriptomes.

Protein-to-genome alignment with miniprot

MOTIVATION

Protein-to-genome alignment is critical to annotating genes in non-model organisms. While there are a few tools for this purpose, all of them were developed over 10 years ago and did not incorporate the latest advances in alignment algorithms. They are inefficient and could not keep up with the rapid production of new genomes and quickly growing protein databases.

RESULTS

Here, we describe miniprot, a new aligner for mapping protein sequences to a complete genome. Miniprot integrates recent techniques such as k-mer sketch and vectorized dynamic programming. It is tens of times faster than existing tools while achieving comparable accuracy on real data.

AVAILABILITY AND IMPLEMENTATION

https://github.com/lh3/miniport.

Roles of minor spliceosome in intron recognition and the convergence with the better understood major spliceosome

Catalyzed by spliceosomes in the nucleus, RNA splicing removes intronic sequences from precursor RNAs in eukaryotes to generate mature RNA, which also significantly increases proteome complexity and fine-tunes gene expression. Most metazoans have two coexisting spliceosomes; the major spliceosome, which removes >99.5% of introns, and the minor spliceosome, which removes far fewer introns (only 770 at present have been predicted in the human genome). Both spliceosomes are large and dynamic machineries, each consisting of five small nuclear RNAs (snRNAs) and more than 100 proteins. However, the dynamic assembly, catalysis, and protein composition of the minor spliceosome are still poorly understood. With different splicing signals, minor introns are rare and usually distributed alone and flanked by major introns in genes, raising questions of how they are recognized by the minor spliceosome and how their processing deals with the splicing of neighboring major introns. Due to large numbers of introns and close similarities between the two machinery, cooperative, and competitive recognition by the two spliceosomes has been investigated. Functionally, many minor-intron-containing genes are evolutionarily conserved and essential. Mutations in the minor spliceosome exhibit a variety of developmental defects in plants and animals and are linked to numerous human diseases. Here, we review recent progress in the understanding of minor splicing, compare currently known components of the two spliceosomes, survey minor introns in a wide range of organisms, discuss cooperation and competition of the two spliceosomes in splicing of minor-intron-containing genes, and contributions of minor splicing mutations in development and diseases. This article is categorized under: RNA Processing > Processing of Small RNAs RNA Processing > Splicing Mechanisms RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry.

Biogenesis of telomerase RNA from a protein-coding mRNA precursor

Telomerase is a eukaryotic ribonucleoprotein (RNP) enzyme that adds DNA repeats onto chromosome ends to maintain genomic stability and confer cellular immortality in cancer and stem cells. The telomerase RNA (TER) component is essential for telomerase catalytic activity and provides the template for telomeric DNA synthesis. The biogenesis of TERs is extremely divergent across eukaryotic kingdoms, employing distinct types of transcription machinery and processing pathways. In ciliates and plants, TERs are transcribed by RNA polymerase III (Pol III), while animal and ascomycete fungal TERs are transcribed by RNA Pol II and share biogenesis pathways with small nucleolar RNA (snoRNA) and small nuclear RNA (snRNA), respectively. Here, we report an unprecedented messenger RNA (mRNA)-derived biogenesis pathway for the 1,291 nucleotide TER from the basidiomycete fungus Ustilago maydis. The U. maydis TER (UmTER) contains a 5'-monophosphate, distinct from the 5' 2,2,7-trimethylguanosine (TMG) cap common to animal and ascomycete fungal TERs. The mature UmTER is processed from the 3'-untranslated region (3'-UTR) of a larger RNA precursor that possesses characteristics of mRNA including a 5' 7-methyl-guanosine (m7G) cap, alternative splicing of introns, and a poly(A) tail. Moreover, this mRNA transcript encodes a protein called Early meiotic induction protein 1 (Emi1) that is conserved across dikaryotic fungi. A recombinant UmTER precursor expressed from an mRNA promoter is processed correctly to yield mature UmTER, confirming an mRNA-processing pathway for producing TER. Our findings expand the plethora of TER biogenesis mechanisms and demonstrate a pathway for producing a functional long noncoding RNA from a protein-coding mRNA precursor.

A high-resolution single-molecule sequencing-based Arabidopsis transcriptome using novel methods of Iso-seq analysis

BACKGROUND

Accurate and comprehensive annotation of transcript sequences is essential for transcript quantification and differential gene and transcript expression analysis. Single-molecule long-read sequencing technologies provide improved integrity of transcript structures including alternative splicing, and transcription start and polyadenylation sites. However, accuracy is significantly affected by sequencing errors, mRNA degradation, or incomplete cDNA synthesis.

RESULTS

We present a new and comprehensive Arabidopsis thaliana Reference Transcript Dataset 3 (AtRTD3). AtRTD3 contains over 169,000 transcripts-twice that of the best current Arabidopsis transcriptome and including over 1500 novel genes. Seventy-eight percent of transcripts are from Iso-seq with accurately defined splice junctions and transcription start and end sites. We develop novel methods to determine splice junctions and transcription start and end sites accurately. Mismatch profiles around splice junctions provide a powerful feature to distinguish correct splice junctions and remove false splice junctions. Stratified approaches identify high-confidence transcription start and end sites and remove fragmentary transcripts due to degradation. AtRTD3 is a major improvement over existing transcriptomes as demonstrated by analysis of an Arabidopsis cold response RNA-seq time-series. AtRTD3 provides higher resolution of transcript expression profiling and identifies cold-induced differential transcription start and polyadenylation site usage.

CONCLUSIONS

AtRTD3 is the most comprehensive Arabidopsis transcriptome currently. It improves the precision of differential gene and transcript expression, differential alternative splicing, and transcription start/end site usage analysis from RNA-seq data. The novel methods for identifying accurate splice junctions and transcription start/end sites are widely applicable and will improve single-molecule sequencing analysis from any species.

Major and minor U small nuclear RNAs genes characterization in a neotropical fish genome: Chromosomal remodeling and repeat units dispersion in Parodontidae

In association with many proteins, small nuclear RNAs (snRNAs) organize the spliceosomes that play a significant role in processing precursor mRNAs during gene expression. According to snRNAs genic arrangements, two kinds of spliceosomes (major and minor) can be organized into eukaryotic cells. Although in situ localization of U1 and U2 snDNAs have been performed in fish karyotypes, studies with genomic characterization and functionality of U snRNAs integrated into chromosomal changes on Teleostei are still scarce. This study aimed to achieve a genomic characterization of the U snRNAs genes in Apareiodon sp. (2n = 54, ZZ/ZW), apply these data to recognize functional/defective copies, and map chromosomal changes involving snDNAs in Parodontidae species karyotype diversification. Nine snRNA multigene families (U1, U2, U4, U5, U6, U11, U12, U4atac and U6atac) arranged in putatively functional copies in the genome were analyzed. Proximal Sequence Elements (PSE) and TATA-box promoters occurrence, besides an entire transcribed region and conserved secondary structures, qualify them for spliceosome activity. In addition, several defective copies or pseudogenes were identified for the snRNAs that make up the major spliceosome. In situ localization of snDNAs in five species of Parodontidae demonstrated that U1, U2, and U4 snDNAs were involved in chromosomal location changes or units dispersion. The U snRNAs defective/pseudogenes units dispersion could be favored by the probable occurrence of active retrotransposition enzymes in the Apareiodon genome. The U2 and U4 snDNAs sites were involved in independent events in the differentiation of sex chromosomes among Parodontidae lineages. The study characterized U snRNA genes that compose major and minor spliceosomes in the Apareiodon sp. genome and proposes that their defective copies trigger chromosome differentiation and diversification events in Parodontidae.

Alternative splicing modulation by G-quadruplexes

Alternative splicing is central to metazoan gene regulation, but the regulatory mechanisms are incompletely understood. Here, we show that G-quadruplex (G4) motifs are enriched ~3-fold near splice junctions. The importance of G4s in RNA is emphasised by a higher enrichment for the non-template strand. RNA-seq data from mouse and human neurons reveals an enrichment of G4s at exons that were skipped following depolarisation induced by potassium chloride. We validate the formation of stable RNA G4s for three candidate splice sites by circular dichroism spectroscopy, UV-melting and fluorescence measurements. Moreover, we find that sQTLs are enriched at G4s, and a minigene experiment provides further support for their role in promoting exon inclusion. Analysis of >1,800 high-throughput experiments reveals multiple RNA binding proteins associated with G4s. Finally, exploration of G4 motifs across eleven species shows strong enrichment at splice sites in mammals and birds, suggesting an evolutionary conserved splice regulatory mechanism.

Here the authors shows that G-quadruplexes, non-canonical DNA/RNA structures, can have a direct impact on alternative splicing and that binding of splicing regulators is affected by their presence.

U1RNP/lncRNA/Transcription Cycle Axis Promotes Tumorigenesis of Hepatocellular Carcinoma

As a component of the spliceosome, U1 small nuclear ribonucleoproteins (U1RNPs) play critical roles in RNA splicing, and recent studies have shown that U1RNPs could recruit long non-coding RNAs (lncRNAs) to chromatin which are involved in cancer development. However, the interplay of U1 snRNP, lncRNAs and downstream genes and signaling pathways are insufficiently understood in hepatocellular carcinoma (HCC). The expression of U1RNPs was found to be significantly higher in tumors than normal tissues in liver hepatocellular carcinomas of The Cancer Genome Atlas (TCGA-LIHC) dataset. LncRNAs with potential U1-binding sites (termed U1-lncRNAs) were found to be mostly located in the nucleus and their expression was higher in tumor than in normal tissues Bioinformatic analysis indicated that U1-lncRNAs worked with RNA-binding proteins and regulated the transcription cycle in HCC. A U1-lncRNA risk model was constructed using a TCGA dataset, and the AUCs of this risk model to predict 1-, 3- and 5-year overall survival were 0.82, 0.84 and 0.8, respectively. Furthermore, silencing of the small nuclear ribonucleoprotein D2 polypeptide (SNRPD2) resulted in impaired proliferation, G1/M cell cycle arrest and downregulation of transcription-cycle-related genes in HCC cell lines. Taken together, these results indicate that U1RNPs interact with lncRNAs and promote the transcription cycle process in HCC, which suggests that these could be novel biomarkers in the clinical management of HCC.

PRMT5 regulates ATF4 transcript splicing and oxidative stress response

Protein methyltransferase 5 (PRMT5) symmetrically dimethylates arginine residues leading to regulation of transcription and splicing programs. Although PRMT5 has emerged as an attractive oncology target, the molecular determinants of PRMT5 dependency in cancer remain incompletely understood. Our transcriptomic analysis identified PRMT5 regulation of the activating transcription factor 4 (ATF4) pathway in acute myelogenous leukemia (AML). PRMT5 inhibition resulted in the expression of unstable, intron-retaining ATF4 mRNA that is detained in the nucleus. Concurrently, the decrease in the spliced cytoplasmic transcript of ATF4 led to lower levels of ATF4 protein and downregulation of ATF4 target genes. Upon loss of functional PRMT5, cells with low ATF4 displayed increased oxidative stress, growth arrest, and cellular senescence. Interestingly, leukemia cells with EVI1 oncogene overexpression demonstrated dependence on PRMT5 function. EVI1 and ATF4 regulated gene signatures were inversely correlated. We show that EVI1-high AML cells have reduced ATF4 levels, elevated baseline reactive oxygen species and increased sensitivity to PRMT5 inhibition. Thus, EVI1-high cells demonstrate dependence on PRMT5 function and regulation of oxidative stress response. Overall, our findings identify the PRMT5-ATF4 axis to be safeguarding the cellular redox balance that is especially important in high oxidative stress states, such as those that occur with EVI1 overexpression.

Non-Canonical Splicing and Its Implications in Brain Physiology and Cancer

The advance of experimental and computational techniques has allowed us to highlight the existence of numerous different mechanisms of RNA maturation, which have been so far unknown. Besides canonical splicing, consisting of the removal of introns from pre-mRNA molecules, non-canonical splicing events may occur to further increase the regulatory and coding potential of the human genome. Among these, splicing of microexons, recursive splicing and biogenesis of circular and chimeric RNAs through back-splicing and trans-splicing processes, respectively, all contribute to expanding the repertoire of RNA transcripts with newly acquired regulatory functions. Interestingly, these non-canonical splicing events seem to occur more frequently in the central nervous system, affecting neuronal development and differentiation programs with important implications on brain physiology. Coherently, dysregulation of non-canonical RNA processing events is associated with brain disorders, including brain tumours. Herein, we summarize the current knowledge on molecular and regulatory mechanisms underlying canonical and non-canonical splicing events with particular emphasis on cis-acting elements and trans-acting factors that all together orchestrate splicing catalysis reactions and decisions. Lastly, we review the impact of non-canonical splicing on brain physiology and pathology and how unconventional splicing mechanisms may be targeted or exploited for novel therapeutic strategies in cancer.

Non-canonical splice junction processing increases the diversity of RBFOX2 splicing isoforms

The underlying mechanisms of splicing regulation through non-canonical splice junction processing remain largely unknown. Here, we identified two RBFOX2 splicing isoforms by alternative 3' splice site selection of exon 9; the non-canonical splice junction processed RBFOX2 transcript (RBFOX2-N.C.) was expressed by the selection of the 3' splice GG acceptor sequence. The cytoplasmic localization of RBFOX2-C., a canonical splice junction-processed RBFOX2 transcript, was different from that of RBFOX2-N.C., which showed nuclear localization. In addition, we confirmed that RBFOX2-C. showed a significantly stronger localization into stress granules than RBFOX2-N.C. upon sodium arsenite treatment. Next, we investigated the importance of non-canonical 3' splice GG sequence selection of specific cis-regulatory elements using minigene constructs of the RBFOX2 gene. We found that the non-canonical 3' splice GG sequence and suboptimal branch point site adjacent region were critical for RBFOX2-N.C. expression through a non-canonical 3' splice selection. Our results suggest a regulatory mechanism for the non-canonical 3' splice selection in the RBFOX2 gene, providing a basis for studies related to the regulation of alternative pre-mRNA splicing through non-canonical splice junction processing.

Identification and Quantification of Microexons Using Bulk and Single-Cell RNA-Seq Data

The analysis of RNA-seq has greatly improved the characterization and understanding of the transcriptome. In particular, RNA-seq experiments have extended catalogs of alternative splicing events. However, the analysis of RNAs-seq data for detection and quantification of microexons, extremely short exons of length up to 30 nt, require specialized computational workflows. Here, we describe MicroExonator, a reproducible computational workflow for microexon splicing analysis using bulk or single-cell RNA-seq data.

A novel role for nucleolin in splice site selection

Latent 5' splice sites, not normally used, are highly abundant in human introns, but are activated under stress and in cancer, generating thousands of nonsense mRNAs. A previously proposed mechanism to suppress latent splicing was shown to be independent of NMD, with a pivotal role for initiator-tRNA independent of protein translation. To further elucidate this mechanism, we searched for nuclear proteins directly bound to initiator-tRNA. Starting with UV-crosslinking, we identified nucleolin (NCL) interacting directly and specifically with initiator-tRNA in the nucleus, but not in the cytoplasm. Next, we show the association of ini-tRNA and NCL with pre-mRNA. We further show that recovery of suppression of latent splicing by initiator-tRNA complementation is NCL dependent. Finally, upon nucleolin knockdown we show activation of latent splicing in hundreds of coding transcripts having important cellular functions. We thus propose nucleolin, a component of the endogenous spliceosome, through its direct binding to initiator-tRNA and its effect on latent splicing, as the first protein of a nuclear quality control mechanism regulating splice site selection to protect cells from latent splicing that can generate defective mRNAs.

Alternative splicing and liver disease

Alternative splicing produces complex and dynamic changes in the protein isoforms that are necessary for the proper biological functioning of the metabolic pathways involved in liver development and hepatocyte homeostasis. Changes in the physiological state of alternatively spliced forms are increasingly linked to liver pathologies. This may occur when the expression or function of the set of proteins controlling the alternative splicing processes are altered by external effectors such as oxidative stress and other environmental variations. Studies addressing these modifications reveal a complex interplay between the expression levels of different proteins that regulate the alternative splicing process as well as the changes in alternative splicing. This interplay results in a cascade of different protein isoforms that correlate with the progression of non-alcoholic fatty liver disease, hepatocellular carcinoma, and alcoholic liver disease. However, research on the detailed molecular mechanism underlying the production of these isoforms is needed. It is imperative to identify the physiological processes affected by the differentially spliced isoforms and confirm their role on the onset and maintenance of the pathology. This is required to design potential therapeutic approaches targeting the key splicing changes to revert the pathological condition as well as identify prognostic markers. In this review, we describe the complexity of the splicing process through an example to encourage researchers to go down this path. Subsequently, rather than a catalog of splicing events we have hand-picked and discuss a few selected studies of specific liver pathologies and suggested ways to focus research on these areas.

Pseudoexon activation in disease by non-splice site deep intronic sequence variation - wild type pseudoexons constitute high-risk sites in the human genome

Accuracy of pre-messenger RNA (pre-mRNA) splicing is crucial for normal gene expression. Complex regulation supports the spliceosomal distinction between authentic exons and the many seemingly functional splice sites delimiting pseudoexons. Pseudoexons are nonfunctional intronic sequences that can be activated for aberrant inclusion in mRNA, which may cause disease. Pseudoexon activation is very challenging to predict, in particular when activation occurs by sequence variants that alter the splicing regulatory environment without directly affecting splice sites. As pseudoexon inclusion often evades detection due to activation of nonsense-mediated mRNA decay, and because conventional diagnostic procedures miss deep intronic sequence variation, pseudoexon activation is a heavily underreported disease mechanism. Pseudoexon characteristics have mainly been studied based on in silico predicted sequences. Moreover, because recognition of sequence variants that create or strengthen splice sites is possible by comparison with well-established consensus sequences, this type of pseudoexon activation is by far the most frequently reported. Here we review all known human disease-associated pseudoexons that carry functional splice sites and are activated by deep intronic sequence variants located outside splice site sequences. We delineate common characteristics that make this type of wild type pseudoexons distinct high-risk sites in the human genome.

Spliceator: multi-species splice site prediction using convolutional neural networks

Background

Ab initio prediction of splice sites is an essential step in eukaryotic genome annotation. Recent predictors have exploited Deep Learning algorithms and reliable gene structures from model organisms. However, Deep Learning methods for non-model organisms are lacking.

Results

We developed Spliceator to predict splice sites in a wide range of species, including model and non-model organisms. Spliceator uses a convolutional neural network and is trained on carefully validated data from over 100 organisms. We show that Spliceator achieves consistently high accuracy (89–92%) compared to existing methods on independent benchmarks from human, fish, fly, worm, plant and protist organisms.

Conclusions

Spliceator is a new Deep Learning method trained on high-quality data, which can be used to predict splice sites in diverse organisms, ranging from human to protists, with consistently high accuracy.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04471-3.

An ADAR1-dependent RNA editing event in the cyclin-dependent kinase CDK13 promotes thyroid cancer hallmarks

Background

Adenosine deaminases acting on RNA (ADARs) modify many cellular RNAs by catalyzing the conversion of adenosine to inosine (A-to-I), and their deregulation is associated with several cancers. We recently showed that A-to-I editing is elevated in thyroid tumors and that ADAR1 is functionally important for thyroid cancer cell progression. The downstream effectors regulated or edited by ADAR1 and the significance of ADAR1 deregulation in thyroid cancer remain, however, poorly defined.

Methods

We performed whole transcriptome sequencing to determine the consequences of ADAR1 deregulation for global gene expression, RNA splicing and editing. The effects of gene silencing or RNA editing were investigated by analyzing cell viability, proliferation, invasion and subnuclear localization, and by protein and gene expression analysis.

Results

We report an oncogenic function for CDK13 in thyroid cancer and identify a new ADAR1-dependent RNA editing event that occurs in the coding region of its transcript. CDK13 was significantly over-edited (c.308A > G) in tumor samples and functional analysis revealed that this editing event promoted cancer cell hallmarks. Finally, we show that CDK13 editing increases the nucleolar abundance of the protein, and that this event might explain, at least partly, the global change in splicing produced by ADAR1 deregulation.

Conclusions

Overall, our data support A-to-I editing as an important pathway in cancer progression and highlight novel mechanisms that might be used therapeutically in thyroid and other cancers.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12943-021-01401-y.

Splicing Outcomes of 5' Splice Site GT>GC Variants That Generate Wild-Type Transcripts Differ Significantly Between Full-Length and Minigene Splicing Assays

Combining data derived from a meta-analysis of human disease-associated 5' splice site GT>GC (i.e., +2T>C) variants and a cell culture-based full-length gene splicing assay (FLGSA) of forward engineered +2T>C substitutions, we recently estimated that ∼15-18% of +2T>C variants can generate up to 84% wild-type transcripts relative to their wild-type counterparts. Herein, we analyzed the splicing outcomes of 20 +2T>C variants that generate some wild-type transcripts in two minigene assays. We found a high discordance rate in terms of the generation of wild-type transcripts, not only between FLGSA and the minigene assays but also between the different minigene assays. In the pET01 context, all 20 wild-type minigene constructs generated the expected wild-type transcripts; of the 20 corresponding variant minigene constructs, 14 (70%) generated wild-type transcripts. In the pSPL3 context, only 18 of the 20 wild-type minigene constructs generated the expected wild-type transcripts whereas 8 of the 18 (44%) corresponding variant minigene constructs generated wild-type transcripts. Thus, in the context of a particular type of variant, we raise awareness of the limitations of minigene splicing assays and emphasize the importance of sequence context in regulating splicing. Whether or not our findings apply to other types of splice-altering variant remains to be investigated.

Chromosomal instability by mutations in the novel minor spliceosome component CENATAC

Aneuploidy is the leading cause of miscarriage and congenital birth defects, and a hallmark of cancer. Despite this strong association with human disease, the genetic causes of aneuploidy remain largely unknown. Through exome sequencing of patients with constitutional mosaic aneuploidy, we identified biallelic truncating mutations in CENATAC (CCDC84). We show that CENATAC is a novel component of the minor (U12-dependent) spliceosome that promotes splicing of a specific, rare minor intron subtype. This subtype is characterized by AT-AN splice sites and relatively high basal levels of intron retention. CENATAC depletion or expression of disease mutants resulted in excessive retention of AT-AN minor introns in ˜ 100 genes enriched for nucleocytoplasmic transport and cell cycle regulators, and caused chromosome segregation errors. Our findings reveal selectivity in minor intron splicing and suggest a link between minor spliceosome defects and constitutional aneuploidy in humans.

U5 snRNA Interactions With Exons Ensure Splicing Precision

Imperfect conservation of human pre-mRNA splice sites is necessary to produce alternative isoforms. This flexibility is combined with the precision of the message reading frame. Apart from intron-termini GU_AG and the branchpoint A, the most conserved are the exon-end guanine and +5G of the intron start. Association between these guanines cannot be explained solely by base-pairing with U1 snRNA in the early spliceosome complex. U6 succeeds U1 and pairs +5G in the pre-catalytic spliceosome, while U5 binds the exon end. Current U5 snRNA reconstructions by CryoEM cannot explain the conservation of the exon-end G. Conversely, human mutation analyses show that guanines of both exon termini can suppress splicing mutations. Our U5 hypothesis explains the mechanism of splicing precision and the role of these conserved guanines in the pre-catalytic spliceosome. We propose: (1) optimal binding register for human exons and U5-the exon junction positioned at U5Loop1 C₃₉|C₃₈; (2) common mechanism for base-pairing of human U5 snRNA with diverse exons and bacterial Ll.LtrB intron with new loci in retrotransposition-guided by base pair geometry; and (3) U5 plays a significant role in specific exon recognition in the pre-catalytic spliceosome. Statistical analyses showed increased U5 Watson-Crick pairs with the 5'exon in the absence of +5G at the intron start. In 5'exon positions -3 and -5, this effect is specific to U5 snRNA rather than U1 snRNA of the early spliceosome. Increased U5 Watson-Crick pairs with 3'exon position +1 coincide with substitutions of the conserved -3C at the intron 3'end. Based on mutation and X-ray evidence, we propose that -3C pairs with U2 G₃₁ juxtaposing the branchpoint and the 3'intron end. The intron-termini pair, formed in the pre-catalytic spliceosome to be ready for transition after branching, and the early involvement of the 3'intron end ensure that the 3'exon contacts U5 in the pre-catalytic complex. We suggest that splicing precision is safeguarded cooperatively by U5, U6, and U2 snRNAs that stabilize the pre-catalytic complex by Watson-Crick base pairing. In addition, our new U5 model explains the splicing effect of exon-start +1G mutations: U5 Watson-Crick pairs with exon +2C/+3G strongly promote exon inclusion. We discuss potential applications for snRNA therapeutics and gene repair by reverse splicing.

Benchmarking deep learning splice prediction tools using functional splice assays

Hereditary disorders are frequently caused by genetic variants that affect pre-messenger RNA splicing. Though genetic variants in the canonical splice motifs are almost always disrupting splicing, the pathogenicity of variants in the noncanonical splice sites (NCSS) and deep intronic (DI) regions are difficult to predict. Multiple splice prediction tools have been developed for this purpose, with the latest tools employing deep learning algorithms. We benchmarked established and deep learning splice prediction tools on published gold standard sets of 71 NCSS and 81 DI variants in the ABCA4 gene and 61 NCSS variants in the MYBPC3 gene with functional assessment in midigene and minigene splice assays. The selection of splice prediction tools included CADD, DSSP, GeneSplicer, MaxEntScan, MMSplice, NNSPLICE, SPIDEX, SpliceAI, SpliceRover, and SpliceSiteFinder-like. The best-performing splice prediction tool for the different variants was SpliceRover for ABCA4 NCSS variants, SpliceAI for ABCA4 DI variants, and the Alamut 3/4 consensus approach (GeneSplicer, MaxEntScacn, NNSPLICE and SpliceSiteFinder-like) for NCSS variants in MYBPC3 based on the area under the receiver operator curve. Overall, the performance in a real-time clinical setting is much more modest than reported by the developers of the tools.

Splicing regulation in hematopoiesis

PURPOSE OF REVIEW

Splicing mutations are among the most recurrent genetic perturbations in hematological malignancies, highlighting an important impact of splicing regulation in hematopoietic development. However, compared to our understanding of splicing factor mutations in hematological malignancies, studies of splicing components and alternative splicing in normal hematopoiesis have been less well investigated. Here, we outline the most recent findings on splicing regulation in normal hematopoiesis and discuss the important questions in the field.

RECENT FINDINGS

Recent studies have highlighted the critical role of splicing regulation in hematopoiesis, including characterization of splicing components in normal hematopoiesis, investigation of transcriptional alterations on splicing, and identification of stage-specific alternative splicing events during hematopoietic development.

SUMMARY

These interesting findings provide insights on hematopoietic regulation at a co-transcriptional level. More high-throughput RNA ribonucleic acid (RNA) sequencing and functional genomic screens are needed to advance our knowledge of critical alternative splicing patterns in shaping hematopoiesis.

Expansion and transformation of the minor spliceosomal system in the slime mold Physarum polycephalum

Spliceosomal introns interrupt nuclear genes and are removed from RNA transcripts ("spliced") by machinery called spliceosomes. Although the vast majority of spliceosomal introns are removed by the so-called major (or "U2") spliceosome, diverse eukaryotes also contain a rare second form, the minor ("U12") spliceosome, and associated ("U12-type") introns.^1-3 In all characterized species, U12-type introns are distinguished by several features, including being rare in the genome (∼0.5% of all introns),^4-6 containing extended evolutionarily conserved splicing motifs,^4,5,7,8 being generally ancient,^9,10 and being inefficiently spliced.^11-13 Here, we report a remarkable exception in the slime mold Physarum polycephalum. The P. polycephalum genome contains >20,000 U12-type introns-25 times more than any other species-enriched in a diversity of non-canonical splice boundaries as well as transformed splicing signals that appear to have co-evolved with the spliceosome due to massive gain of efficiently spliced U12-type introns. These results reveal an unappreciated dynamism of minor spliceosomal introns and spliceosomal introns in general.

Co-transcriptional splicing efficiencies differ within genes and between cell types

Pre-mRNA splicing is carried out by the spliceosome and involves splice site recognition, removal of introns, and ligation of exons. Components of the spliceosome have been shown to interact with the elongating RNA polymerase II (RNAPII) which is thought to allow splicing to occur concurrently with transcription. However, little is known about the regulation and efficiency of co-transcriptional splicing in human cells. In this study, we used Bru-seq and BruChase-seq to determine the co-transcriptional splicing efficiencies of 17,000 introns expressed across 6 human cell lines. We found that less than half of all introns across these 6 cell lines were co-transcriptionally spliced. Splicing efficiencies for individual introns showed variations across cell lines, suggesting that splicing may be regulated in a cell-type specific manner. Moreover, the splicing efficiency of introns varied within genes. The efficiency of co-transcriptional splicing did not correlate with gene length, intron position, splice site strengths, or the intron/neighboring exons GC content. However, we identified binding signals from multiple RNA binding proteins (RBPs) that correlated with splicing efficiency, including core spliceosomal machinery components-such as SF3B4, U2AF1 and U2AF2 showing higher binding signals in poorly spliced introns. In addition, multiple RBPs, such as BUD13, PUM1 and SND1, showed preferential binding in exons that flank introns with high splicing efficiencies. The nascent RNA splicing patterns presented here across multiple cell types add to our understanding of the complexity in RNA splicing, wherein RNA-binding proteins may play important roles in determining splicing outcomes in a cell type- and intron-specific manner.

Minor intron retention drives clonal hematopoietic disorders and diverse cancer predisposition

Most eukaryotes harbor two distinct pre-mRNA splicing machineries: the major spliceosome, which removes >99% of introns, and the minor spliceosome, which removes rare, evolutionarily conserved introns. Although hypothesized to serve important regulatory functions, physiologic roles of the minor spliceosome are not well understood. For example, the minor spliceosome component ZRSR2 is subject to recurrent, leukemia-associated mutations, yet functional connections among minor introns, hematopoiesis and cancers are unclear. Here, we identify that impaired minor intron excision via ZRSR2 loss enhances hematopoietic stem cell self-renewal. CRISPR screens mimicking nonsense-mediated decay of minor intron-containing mRNA species converged on LZTR1, a regulator of RAS-related GTPases. LZTR1 minor intron retention was also discovered in the RASopathy Noonan syndrome, due to intronic mutations disrupting splicing and diverse solid tumors. These data uncover minor intron recognition as a regulator of hematopoiesis, noncoding mutations within minor introns as potential cancer drivers and links among ZRSR2 mutations, LZTR1 regulation and leukemias.

BBS1 branchpoint variant is associated with non-syndromic retinitis pigmentosa

BACKGROUND

Inherited retinal diseases (IRDs) can be caused by variants in >270 genes. The Bardet-Biedl syndrome 1 (BBS1) gene is one of these genes and may be associated with syndromic and non-syndromic autosomal recessive retinitis pigmentosa (RP). Here, we identified a branchpoint variant in BBS1 and assessed its pathogenicity by in vitro functional analysis.

METHODS

Whole genome sequencing was performed for three unrelated monoallelic BBS1 cases with non-syndromic RP. A fourth case received MGCM 105 gene panel analysis. Functional analysis using a midigene splice assay was performed for the putative pathogenic branchpoint variant in BBS1. After confirmation of its pathogenicity, patients were clinically re-evaluated, including assessment of non-ocular features of Bardet-Biedl syndrome.

RESULTS

Clinical assessments of probands showed that all individuals displayed non-syndromic RP with macular involvement. Through detailed variant analysis and prioritisation, two pathogenic variants in BBS1, the most common missense variant, c.1169T>G (p.(Met390Arg)), and a branchpoint variant, c.592-21A>T, were identified. Segregation analysis confirmed that in all families, probands were compound heterozygous for c.1169T>G and c.592-21A>T. Functional analysis of the branchpoint variant revealed a complex splicing defect including exon 8 and exon 7/8 skipping, and partial in-frame deletion of exon 8.

CONCLUSION

A putative severe branchpoint variant in BBS1, together with a mild missense variant, underlies non-syndromic RP in four unrelated individuals. To our knowledge, this is the first report of a pathogenic branchpoint variant in IRDs that results in a complex splice defect. In addition, this research highlights the importance of the analysis of non-coding regions in order to provide a conclusive molecular diagnosis.

ZRSR1 cooperates with ZRSR2 in regulating splicing of U12-type introns in murine hematopoietic cells

Recurrent loss-of-function mutations of spliceosome gene, ZRSR2, occur in myelodysplastic syndromes (MDS). Mutation/loss of ZRSR2 in human myeloid cells primarily causes impaired splicing of the U12-type introns. In order to further investigate the role of this splice factor in RNA splicing and hematopoietic development, we generated mice lacking ZRSR2. Unexpectedly, Zrsr2-deficient mice developed normal hematopoiesis with no abnormalities in myeloid differentiation evident in either young or ≥1-year old knockout mice. Repopulation ability of Zrsr2-deficient hematopoietic stem cells was also unaffected in both competitive and non-competitive reconstitution assays. Myeloid progenitors lacking ZRSR2 exhibited mis-splicing of U12-type introns, however, this phenotype was moderate compared to the ZRSR2-deficient human cells. Our investigations revealed that a closely related homolog, Zrsr1, expressed in the murine hematopoietic cells, but not in human cells contributes to splicing of U12-type introns. Depletion of Zrsr1 in Zrsr2 KO myeloid cells exacerbated retention of the U12-type introns, thus highlighting a collective role of ZRSR1 and ZRSR2 in murine U12-spliceosome. We also demonstrate that aberrant retention of U12-type introns of MAPK9 and MAPK14 leads to their reduced protein expression. Overall, our findings highlight that both ZRSR1 and ZRSR2 are functional components of the murine U12-spliceosome, and depletion of both proteins is required to accurately model ZRSR2-mutant MDS in mice.

2passtools: two-pass alignment using machine-learning-filtered splice junctions increases the accuracy of intron detection in long-read RNA sequencing

Transcription of eukaryotic genomes involves complex alternative processing of RNAs. Sequencing of full-length RNAs using long reads reveals the true complexity of processing. However, the relatively high error rates of long-read sequencing technologies can reduce the accuracy of intron identification. Here we apply alignment metrics and machine-learning-derived sequence information to filter spurious splice junctions from long-read alignments and use the remaining junctions to guide realignment in a two-pass approach. This method, available in the software package 2passtools (https://github.com/bartongroup/2passtools), improves the accuracy of spliced alignment and transcriptome assembly for species both with and without existing high-quality annotations.

Early onset severe ATP1A2 epileptic encephalopathy: Clinical characteristics and underlying mutations

BACKGROUND

ATP1A2 mutations cause hemiplegic migraine with or without epilepsy or acute reversible encephalopathy. Typical onset is in adulthood or older childhood without subsequent severe long-term developmental impairments.

AIM

We aimed to describe the manifestations of early onset severe ATP1A2-related epileptic encephalopathy and its underlying mutations in a cohort of seven patients.

METHODS

A retrospective chart review of a cohort of seven patients was conducted. Response to open-label memantine therapy, used off-label due to its NMDA receptor antagonist effects, was assessed by the Global Rating Scale of Change (GRSC) and Clinical Global Impression Scale of Improvement (CGI-I) methodologies. Molecular modeling was performed using PyMol program.

RESULTS

Patients (age 2.5-20 years) had symptom onset at an early age (6 days-1 year). Seizures were either focal or generalized. Common features were: drug resistance, recurrent status epilepticus, etc., severe developmental delay with episodes of acute severe encephalopathy often with headaches, dystonias, hemiplegias, seizures, and developmental regression. All had variants predicted to be disease causing (p.Ile293Met, p.Glu1000Lys, c.1017+5G>A, p.Leu809Arg, and 3 patients with p.Met813Lys). Modeling revealed that mutations interfered with ATP1A2 ion binding and translocation sites. Memantine, given to five, was tolerated in all (mean treatment: 2.3 years, range 6 weeks-4.8 years) with some improvements reported in all five.

CONCLUSIONS

Our observations describe a distinctive clinical profile of seven unrelated probands with early onset severe ATP1A2-related epileptic encephalopathy, provide insights into structure-function relationships of ATP1A2 mutations, and support further studies of NMDAR antagonist therapy in ATP1A2-encephalopathy.

Structure of the activated human minor spliceosome

The minor spliceosome mediates splicing of the rare but essential U12-type precursor messenger RNA. Here, we report the atomic features of the activated human minor spliceosome determined by cryo-electron microscopy at 2.9-angstrom resolution. The 5' splice site and branch point sequence of the U12-type intron are recognized by the U6atac and U12 small nuclear RNAs (snRNAs), respectively. Five newly identified proteins stabilize the conformation of the catalytic center: The zinc finger protein SCNM1 functionally mimics the SF3a complex of the major spliceosome, the RBM48-ARMC7 complex binds the γ-monomethyl phosphate cap at the 5' end of U6atac snRNA, the U-box protein PPIL2 coordinates loop I of U5 snRNA and stabilizes U5 small nuclear ribonucleoprotein (snRNP), and CRIPT stabilizes U12 snRNP. Our study provides a framework for the mechanistic understanding of the function of the human minor spliceosome.

MicroExonator enables systematic discovery and quantification of microexons across mouse embryonic development

BACKGROUND

Microexons, exons that are ≤ 30 nucleotides, are a highly conserved and dynamically regulated set of cassette exons. They have key roles in nervous system development and function, as evidenced by recent results demonstrating the impact of microexons on behaviour and cognition. However, microexons are often overlooked due to the difficulty of detecting them using standard RNA-seq aligners.

RESULTS

Here, we present MicroExonator, a novel pipeline for reproducible de novo discovery and quantification of microexons. We process 289 RNA-seq datasets from eighteen mouse tissues corresponding to nine embryonic and postnatal stages, providing the most comprehensive survey of microexons available for mice. We detect 2984 microexons, 332 of which are differentially spliced throughout mouse embryonic brain development, including 29 that are not present in mouse transcript annotation databases. Unsupervised clustering of microexons based on their inclusion patterns segregates brain tissues by developmental time, and further analysis suggests a key function for microexons in axon growth and synapse formation. Finally, we analyse single-cell RNA-seq data from the mouse visual cortex, and for the first time, we report differential inclusion between neuronal subpopulations, suggesting that some microexons could be cell type-specific.

CONCLUSIONS

MicroExonator facilitates the investigation of microexons in transcriptome studies, particularly when analysing large volumes of data. As a proof of principle, we use MicroExonator to analyse a large collection of both mouse bulk and single-cell RNA-seq datasets. The analyses enabled the discovery of previously uncharacterized microexons, and our study provides a comprehensive microexon inclusion catalogue during mouse development.

Interpretation of BRCA2 Splicing Variants: A Case Series of Challenging Variant Interpretations and the Importance of Functional RNA Analysis

A substantial proportion of pathogenic variants associated with an increased risk of hereditary cancer are sequence variants affecting RNA splicing. The classification of these variants can be complex when both non-functional and functional transcripts are produced from the variant allele. We present four BRCA2 splice site variants with complex variant interpretations (BRCA2 c.68-3T>G, c.68-2A>G, c.425G>T, c.8331+2T>C). Evidence supporting a pathogenic classification is available for each variant, including in silico models, absence in population databases, and published functional data. However, comprehensive RNA analysis showed that some functional transcript may be produced by each variant. BRCA2 c.68-3T>G results in a partial splice defect. For BRCA2 c.68-2A>G and c.425G>T, aberrant splicing was shown to produce a potentially functional, in-frame transcript. BRCA2 c.8331+2T>C may utilize a functional GC donor in place of the wild-type GT donor. The severity of cancer history for carriers of these variants was also assessed using a history weighting algorithm and was not consistent with pathogenic controls (carriers of known pathogenic variants in BRCA2). Due to the conflicting evidence, our laboratory classifies these BRCA2 variants as variants of uncertain significance. This highlights the importance of evaluating new and existing evidence to ensure accurate variant classification and appropriate patient care.

Engineering rules that minimize germline silencing of transgenes in simple extrachromosomal arrays in C. elegans

Transgenes are prone to progressive silencing due to their structure, copy number, and genomic location. In C. elegans, repressive mechanisms are particularly strong in the germline with almost fully penetrant transgene silencing in simple extrachromosomal arrays and frequent silencing of single-copy transgene insertions. A class of non-coding DNA, Periodic An/Tn Clusters (PATCs) can prevent transgene-silencing in repressive chromatin or from small interfering RNAs (piRNAs). Here, we describe design rules (codon-optimization, intron and PATC inclusion, elevated temperature (25 °C), and vector backbone removal) for efficient germline expression from arrays in wildtype animals. We generate web-based tools to analyze PATCs and reagents for the convenient assembly of PATC-rich transgenes. An extensive collection of silencing resistant fluorescent proteins (e.g., gfp, mCherry, and tagBFP) can be used for dissecting germline regulatory elements and a set of enhanced enzymes (Mos1 transposase, Cas9, Cre, and Flp recombinases) enable efficient genetic engineering in C. elegans.