Variation in alternative splicing across human tissues

Disruption of nuclear speckle integrity dysregulates RNA splicing in C9ORF72-FTD/ALS

Expansion of an intronic (GGGGCC)n repeat within the C9ORF72 gene is the most common genetic cause of both frontotemporal dementia (FTD) and amyotrophic lateral sclerosis (ALS) (C9-FTD/ALS), characterized with aberrant repeat RNA foci and noncanonical translation-produced dipeptide repeat (DPR) protein inclusions. Here, we elucidate that the (GGGGCC)n repeat RNA co-localizes with nuclear speckles and alters their phase separation properties and granule dynamics. Moreover, the essential nuclear speckle scaffold protein SRRM2 is sequestered into the poly-GR cytoplasmic inclusions in the C9-FTD/ALS mouse model and patient postmortem tissues, exacerbating the nuclear speckle dysfunction. Impaired nuclear speckle integrity induces global exon skipping and intron retention in human iPSC-derived neurons and causes neuronal toxicity. Similar alternative splicing changes can be found in C9-FTD/ALS patient postmortem tissues. This work identified novel molecular mechanisms of global RNA splicing defects caused by impaired nuclear speckle function in C9-FTD/ALS and revealed novel potential biomarkers or therapeutic targets.

The role of HnrnpF/H as a driver of oligoteratozoospermia

Male subfertility or infertility is a common condition often characterized by men producing a low number of sperm with poor quality. To gain insight into this condition, we performed a quantitative proteomic analysis of semen samples obtained from infertile and fertile men. At least 6 proteins showed significant differences in regulation of alternatively spliced isoforms. To investigate this link between aberrant alternative splicing and production of poor-quality spermatozoa, we overexpressed the hnrnpH/F-orthologue Glorund (Glo) in Drosophila, which was also found to be abundant in poor quality human sperm. Transgenic animals produced low numbers of morphologically defective spermatozoa and aberrant formation of the "dense body," an organelle akin to the mammalian manchette. Furthermore, fertility trials demonstrated that transgenic flies were either completely infertile or highly subfertile. These findings suggest that dysregulation of hnrnpH/F is likely to result in the production of low-quality semen, leading to subfertility or infertility in men.

Regulation of DNA damage response by RNA/DNA-binding proteins: Implications for neurological disorders and aging

RNA-binding proteins (RBPs) are evolutionarily conserved across most forms of life, with an estimated 1500 RBPs in humans. Traditionally associated with post-transcriptional gene regulation, RBPs contribute to nearly every known aspect of RNA biology, including RNA splicing, transport, and decay. In recent years, an increasing subset of RBPs have been recognized for their DNA binding properties and involvement in DNA transactions. We refer to these RBPs with well-characterized DNA binding activity as RNA/DNA binding proteins (RDBPs), many of which are linked to neurological diseases. RDBPs are associated with both nuclear and mitochondrial DNA repair. Furthermore, the presence of intrinsically disordered domains in RDBPs appears to be critical for regulating their diverse interactions and plays a key role in controlling protein aggregation, which is implicated in neurodegeneration. In this review, we discuss the emerging roles of common RDBPs from the heterogeneous nuclear ribonucleoprotein (hnRNP) family, such as TAR DNA binding protein-43 (TDP43) and fused in sarcoma (FUS) in controlling DNA damage response (DDR). We also explore the implications of RDBP pathology in aging and neurodegenerative diseases and provide a prospective on the therapeutic potential of targeting RDBP pathology mediated DDR defects for motor neuron diseases and aging.

Transcriptomics : Approaches to Quantifying Gene Expression and Their Application to Studying the Human Brain

To date, the field of transcriptomics has been characterized by rapid methods development and technological advancement, with new technologies continuously rendering older ones obsolete.This chapter traces the evolution of approaches to quantifying gene expression and provides an overall view of the current state of the field of transcriptomics, its applications to the study of the human brain, and its place in the broader emerging multiomics landscape.

APOER2 splicing repertoire in Alzheimer's disease: Insights from long-read RNA sequencing

Disrupted alternative splicing plays a determinative role in neurological diseases, either as a direct cause or as a driver in disease susceptibility. Transcriptomic profiling of aged human postmortem brain samples has uncovered hundreds of aberrant mRNA splicing events in Alzheimer's disease (AD) brains, associating dysregulated RNA splicing with disease. We previously identified a complex array of alternative splicing combinations across apolipoprotein E receptor 2 (APOER2), a transmembrane receptor that interacts with both the neuroprotective ligand Reelin and the AD-associated risk factor, APOE. Many of the human APOER2 isoforms, predominantly featuring cassette splicing events within functionally important domains, are critical for the receptor's function and ligand interaction. However, a comprehensive repertoire and the functional implications of APOER2 isoforms under both physiological and AD conditions are not fully understood. Here, we present an in-depth analysis of the splicing landscape of human APOER2 isoforms in normal and AD states. Using single-molecule, long-read sequencing, we profiled the entire APOER2 transcript from the parietal cortex and hippocampus of Braak stage IV AD brain tissues along with age-matched controls and investigated several functional properties of APOER2 isoforms. Our findings reveal diverse patterns of cassette exon skipping for APOER2 isoforms, with some showing region-specific expression and others unique to AD-affected brains. Notably, exon 15 of APOER2, which encodes the glycosylation domain, showed less inclusion in AD compared to control in the parietal cortex of females with an APOE ɛ3/ɛ3 genotype. Also, some of these APOER2 isoforms demonstrated changes in cell surface expression, APOE-mediated receptor processing, and synaptic number. These variations are likely critical in inducing synaptic alterations and may contribute to the neuronal dysfunction underlying AD pathogenesis.

Reverse engineering neuron type-specific and type-orthogonal splicing-regulatory networks using single-cell transcriptomes

Cell type-specific alternative splicing (AS) enables differential gene isoform expression between diverse neuron types with distinct identities and functions. Current studies linking individual RNA-binding proteins (RBPs) to AS in a few neuron types underscore the need for holistic modeling. Here, we use network reverse engineering to derive a map of the neuron type-specific AS regulatory landscape from 133 mouse neocortical cell types defined by single-cell transcriptomes. This approach reliably inferred the regulons of 350 RBPs and their cell type-specific activities. Our analysis revealed driving factors delineating neuronal identities, among which we validated Elavl2 as a key RBP for MGE-specific splicing in GABAergic interneurons using an in vitro ESC differentiation system. We also identified a module of exons and candidate regulators specific for long- and short-projection neurons across multiple neuronal classes. This study provides a resource for elucidating splicing regulatory programs that drive neuronal molecular diversity, including those that do not align with gene expression-based classifications.

The Roles of hnRNP Family in the Brain and Brain-Related Disorders

Heterogeneous nuclear ribonucleoproteins (hnRNPs) belong to a complex family of RNA-binding proteins that are essential to control alternative splicing, mRNA trafficking, synaptic plasticity, stress granule formation, cell cycle regulation, and axonal transport. Over the past decade, hnRNPs have been associated with different brain disorders such as Alzheimer's disease, multiple sclerosis, and schizophrenia. Given their essential role in maintaining cell function and integrity, it is not surprising that dysregulated hnRNP levels lead to neurological implications. This review aims to explore the primary functions of hnRNPs in neurons, oligodendrocytes, microglia, and astrocytes, and their roles in brain disorders. We also discuss proteomics and other technologies and their potential for studying and evaluating hnRNPs in brain disorders, including the discovery of new therapeutic targets and possible pharmacological interventions.

Complete characterization of the yak testicular development using accurate full-length transcriptome sequencing

Alternative splicing is a prevalent phenomenon in testicular tissues. Due to the low assembly accuracy of short-read RNA sequencing technology in analyzing post-transcriptional regulatory events, full-length (FL) transcript sequencing is highly demanded to accurately determine FL splicing variants. In this study, we performed FL transcriptome sequencing of testicular tissues from 0.5, 1.5, 2.5, and 4-year-old yaks and 4-year-old cattle-yaks using Oxford Nanopore Technologies. The obtained sequencing data were predicted to have 47,185 open reading frames (ORFs), including 26,630 complete ORFs, detected 7645 fusion transcripts, 15,355 alternative splicing events, 25,798 simple sequence repeats, 7628 transcription factors, and 35,503 long non-coding RNAs. A total of 40,038 novel transcripts were obtained from the sequencing data, and the proportion was almost close to the number of known transcripts identified. Structural analysis and functional annotation of these novel transcripts resulted in the successful annotation of 9568 transcripts, with the highest and lowest annotation numbers in the Nr and KOG databases, respectively. Weighted gene co-expression network analysis revealed the key regulatory pathways and hub genes at various stages of yak testicular development. Our findings enhance our comprehension of transcriptome complexity, contribute to genome annotation refinement, and provide foundational data for further investigations into male sterility in cattle-yaks.

Long-read sequencing for 29 immune cell subsets reveals disease-linked isoforms

Alternative splicing events are a major causal mechanism for complex traits, but they have been understudied due to the limitation of short-read sequencing. Here, we generate a full-length isoform annotation of human immune cells from an individual by long-read sequencing for 29 cell subsets. This contains a number of unannotated transcripts and isoforms such as a read-through transcript of TOMM40-APOE in the Alzheimer's disease locus. We profile characteristics of isoforms and show that repetitive elements significantly explain the diversity of unannotated isoforms, providing insight into the human genome evolution. In addition, some of the isoforms are expressed in a cell-type specific manner, whose alternative 3'-UTRs usage contributes to their specificity. Further, we identify disease-associated isoforms by isoform switch analysis and by integration of several quantitative trait loci analyses with genome-wide association study data. Our findings will promote the elucidation of the mechanism of complex diseases via alternative splicing.

Uncovering a multitude of stage-specific splice variants and putative protein isoforms generated along mouse spermatogenesis

BACKGROUND

Mammalian testis is a highly complex and heterogeneous tissue. This complexity, which mostly derives from spermatogenic cells, is reflected at the transcriptional level, with the largest number of tissue-specific genes and long noncoding RNAs (lncRNAs) compared to other tissues, and one of the highest rates of alternative splicing. Although it is known that adequate alternative-splicing patterns and stage-specific isoforms are critical for successful spermatogenesis, so far only a very limited number of reports have addressed a detailed study of alternative splicing and isoforms along the different spermatogenic stages.

RESULTS

In the present work, using highly purified stage-specific testicular cell populations, we detected 33,002 transcripts expressed throughout mouse spermatogenesis not annotated so far. These include both splice variants of already annotated genes, and of hitherto unannotated genes. Using conservative criteria, we uncovered 13,471 spermatogenic lncRNAs, which reflects the still incomplete annotation of lncRNAs. A distinctive feature of lncRNAs was their lower number of splice variants compared to protein-coding ones, adding to the conclusion that lncRNAs are, in general, less complex than mRNAs. Besides, we identified 2,794 unannotated transcripts with high coding potential (including some arising from yet unannotated genes), many of which encode unnoticed putative testis-specific proteins. Some of the most interesting coding splice variants were chosen, and validated through RT-PCR. Remarkably, the largest number of stage-specific unannotated transcripts are expressed during early meiotic prophase stages, whose study has been scarcely addressed in former transcriptomic analyses.

CONCLUSIONS

We detected a high number of yet unannotated genes and alternatively spliced transcripts along mouse spermatogenesis, hence showing that the transcriptomic diversity of the testis is considerably higher than previously reported. This is especially prominent for specific, underrepresented stages such as those of early meiotic prophase, and its unveiling may constitute a step towards the understanding of their key events.

TREX tetramer disruption alters RNA processing necessary for corticogenesis in THOC6 Intellectual Disability Syndrome

THOC6 variants are the genetic basis of autosomal recessive THOC6 Intellectual Disability Syndrome (TIDS). THOC6 is critical for mammalian Transcription Export complex (TREX) tetramer formation, which is composed of four six-subunit THO monomers. The TREX tetramer facilitates mammalian RNA processing, in addition to the nuclear mRNA export functions of the TREX dimer conserved through yeast. Human and mouse TIDS model systems revealed novel THOC6-dependent, species-specific TREX tetramer functions. Germline biallelic Thoc6 loss-of-function (LOF) variants result in mouse embryonic lethality. Biallelic THOC6 LOF variants reduce the binding affinity of ALYREF to THOC5 without affecting the protein expression of TREX members, implicating impaired TREX tetramer formation. Defects in RNA nuclear export functions were not detected in biallelic THOC6 LOF human neural cells. Instead, mis-splicing was detected in human and mouse neural tissue, revealing novel THOC6-mediated TREX coordination of mRNA processing. We demonstrate that THOC6 is required for key signaling pathways known to regulate the transition from proliferative to neurogenic divisions during human corticogenesis. Together, these findings implicate altered RNA processing in the developmental biology of TIDS neuropathology.

Multiple Sgip1 splice variants inhibit cannabinoid receptor 1 internalization

Alternative splicing can often result in the expression of distinct protein isoforms from a single gene, with specific composition and properties. SH3-containing GRB2-like protein 3-interacting protein 1 (Sgip1) is a brain-enriched protein that regulates clathrin-mediated endocytosis and interferes with the internalization of cannabinoid receptor 1. Several research groups have studied the physiological importance of Sgip1, and four Sgip1 protein isoforms have been described to date, while the NCBI Gene database predicts the expression of 20 splice variants from the Sgip1 gene in mice. In this work, we cloned 15 Sgip1 splice variants from the mouse brain, including 11 novel splice variants. The cloned splice variants differed in exon composition within two Sgip1 regions: the membrane phospholipid-binding domain and the proline-rich region. All the Sgip1 splice isoforms had similar stability and comparable ability to inhibit the internalization of cannabinoid receptor 1. None of the isoforms influenced the internalization of the µ-opioid receptor. We confirm the expression of Sgip1 splice variants described in previous studies or predicted in silico. Our data provide a basis for further studies exploring the significance of Sgip1 splicing, and we suggest a new classification of Sgip1 splice variants to unify their nomenclature.

Decoding the role of aberrant RNA alternative splicing in hepatocellular carcinoma: a comprehensive review

During eukaryotic gene expression, alternative splicing of messenger RNA precursors is critical in increasing protein diversity and regulatory complexity. Multiple transcript isoforms could be produced by alternative splicing from a single gene; they could eventually be translated into protein isoforms with deleted, added, or altered domains or produce transcripts containing premature termination codons that could be targeted by nonsense-mediated mRNA decay. Alternative splicing can generate proteins with similar, different, or even opposite functions. Increasingly strong evidence indicates that abnormal RNA splicing is a prevalent and crucial occurrence in cellular differentiation, tissue advancement, and the development and progression of cancer. Aberrant alternative splicing could affect cancer cell activities such as growth, apoptosis, invasiveness, drug resistance, angiogenesis, and metabolism. This systematic review provides a comprehensive overview of the impact of abnormal RNA alternative splicing on the development and progression of hepatocellular carcinoma.

RBFOX2 regulated EYA3 isoforms partner with SIX4 or ZBTB1 to control transcription during myogenesis

Alternative splicing is a prevalent gene-regulatory mechanism, with over 95% of multi-exon human genes estimated to be alternatively spliced. Here, we describe a tissue-specific, developmentally regulated, highly conserved, and disease-associated alternative splicing event in exon 7 of the eyes absent homolog 3 (Eya3) gene. We discovered that EYA3 expression is vital to the proliferation and differentiation of myoblasts. Genome-wide transcriptomic analysis and mass spectrometry-based proteomic studies identified SIX homeobox 4 (SIX4) and zinc finger and BTB-domain containing 1 (ZBTB1), as major transcription factors that interact with EYA3 to dictate gene expression. EYA3 isoforms differentially regulate transcription, indicating that splicing aids in temporal control of gene expression during muscle cell differentiation. Finally, we identified RNA-binding fox-1 homolog 2 (RBFOX2) as the main regulator of EYA3 splicing. Together, our findings illustrate the interplay between alternative splicing and transcription during myogenesis.

A neuron-specific microexon ablates the novel DNA-binding function of a histone H3K4me0 reader PHF21A

How cell-type-specific chromatin landscapes emerge and progress during metazoan ontogenesis remains an important question. Transcription factors are expressed in a cell-type-specific manner and recruit chromatin-regulatory machinery to specific genomic loci. In contrast, chromatin-regulatory proteins are expressed broadly and are assumed to exert the same intrinsic function across cell types. However, human genetics studies have revealed an unexpected vulnerability of neurodevelopment to chromatin factor mutations with unknown mechanisms. Here, we report that 14 chromatin regulators undergo evolutionary-conserved neuron-specific splicing events involving microexons. Of the 14 chromatin regulators, two are integral components of a histone H3K4 demethylase complex; the catalytic subunit LSD1 and an H3K4me0-reader protein PHF21A adopt neuron-specific forms. We found that canonical PHF21A (PHF21A-c) binds to DNA by AT-hook motif, and the neuronal counterpart PHF21A-n lacks this DNA-binding function yet maintains H3K4me0 recognition intact. In-vitro reconstitution of the canonical and neuronal PHF21A-LSD1 complexes identified the neuronal complex as a hypomorphic H3K4 demethylating machinery with reduced nucleosome engagement. Furthermore, an autism-associated PHF21A missense mutation, 1285 G>A, at the last nucleotide of the common exon immediately upstream of the neuronal microexon led to impaired splicing of PHF21A -n. Thus, ubiquitous chromatin regulatory complexes exert unique intrinsic functions in neurons via alternative splicing of their subunits and potentially contribute to faithful human brain development.

Deciphering RNA splicing logic with interpretable machine learning

Machine learning methods, particularly neural networks trained on large datasets, are transforming how scientists approach scientific discovery and experimental design. However, current state-of-the-art neural networks are limited by their uninterpretability: Despite their excellent accuracy, they cannot describe how they arrived at their predictions. Here, using an "interpretable-by-design" approach, we present a neural network model that provides insights into RNA splicing, a fundamental process in the transfer of genomic information into functional biochemical products. Although we designed our model to emphasize interpretability, its predictive accuracy is on par with state-of-the-art models. To demonstrate the model's interpretability, we introduce a visualization that, for any given exon, allows us to trace and quantify the entire decision process from input sequence to output splicing prediction. Importantly, the model revealed uncharacterized components of the splicing logic, which we experimentally validated. This study highlights how interpretable machine learning can advance scientific discovery.

Alternative RNA Splicing in the Retina: Insights and Perspectives

Alternative splicing is a fundamental and highly regulated post-transcriptional process that enhances transcriptome and proteome diversity. This process is particularly important in neuronal tissues, such as the retina, which exhibit some of the highest levels of differentially spliced genes in the body. Alternative splicing is regulated both temporally and spatially during neuronal development, can be cell-type-specific, and when altered can cause a number of pathologies, including retinal degeneration. Advancements in high-throughput sequencing technologies have facilitated investigations of the alternative splicing landscape of the retina in both healthy and disease states. Additionally, innovations in human stem cell engineering, specifically in the generation of 3D retinal organoids, which recapitulate many aspects of the in vivo retinal microenvironment, have aided studies of the role of alternative splicing in human retinal development and degeneration. Here we review these advances and discuss the ongoing development of strategies for the treatment of alternative splicing-related retinal disease.

Alternative splicing impacts microRNA regulation within coding regions

MicroRNAs (miRNAs) are small non-coding RNA molecules that bind to target sites in different gene regions and regulate post-transcriptional gene expression. Approximately 95% of human multi-exon genes can be spliced alternatively, which enables the production of functionally diverse transcripts and proteins from a single gene. Through alternative splicing, transcripts might lose the exon with the miRNA target site and become unresponsive to miRNA regulation. To check this hypothesis, we studied the role of miRNA target sites in both coding and non-coding regions using six cancer data sets from The Cancer Genome Atlas (TCGA) and Parkinson's disease data from PPMI. First, we predicted miRNA target sites on mRNAs from their sequence using TarPmiR. To check whether alternative splicing interferes with this regulation, we trained linear regression models to predict miRNA expression from transcript expression. Using nested models, we compared the predictive power of transcripts with miRNA target sites in the coding regions to that of transcripts without target sites. Models containing transcripts with target sites perform significantly better. We conclude that alternative splicing does interfere with miRNA regulation by skipping exons with miRNA target sites within the coding region.

Alternative splicing of CEACAM1 by hypoxia-inducible factor-1α enhances tolerance to hepatic ischemia in mice and humans

Although alternative splicing (AS) drives transcriptional responses and cellular adaptation to environmental stresses, its contributions in organ transplantation have not been appreciated. We have shown that carcinoembryonic antigen-related cell adhesion molecule (Ceacam1; CD66a), a transmembrane biliary glycoprotein expressed in epithelial, endothelial, and immune cells, determines donor liver transplant quality. Here, we studied how AS of Ceacam1 affects ischemia-reperfusion injury (IRI) in mouse and human livers. We found that the short cytoplasmic isoform Ceacam1-S increased during early acute and late resolution phases of warm IRI injury in mice. Transfection of Ceacam1-deficient mouse hepatocytes with adenoviral Ceacam1-S mitigated hypoxia-induced loss of cellular adhesion by repressing the Ask1/p-p38 cell death pathway. Nucleic acid-blocking morpholinos, designed to selectively induce Ceacam1-S, protected hepatocyte cultures against temperature-induced stress in vitro. Luciferase and chromatin immunoprecipitation assays identified direct binding of hypoxia-inducible factor-1α (Hif-1α) to the mouse polypyrimidine tract binding protein 1 (Ptbp1) promoter region. Dimethyloxalylglycine protected mouse livers from warm IR stress and hepatocellular damage by inhibiting prolyl hydroxylase domain-containing protein 1 and promoting AS of Ceacam1-S. Last, analysis of 46 human donor liver grafts revealed that CEACAM1-S positively correlated with pretransplant HIF1A expression. This also correlated with better transplant outcomes, including reduced TIMP1, total bilirubin, proinflammatory MCP1, CXCL10 cytokines, immune activation markers IL17A, and incidence of delayed complications from biliary anastomosis. This translational study identified mouse Hif-1α-controlled AS of Ceacam1, through transcriptional regulation of Ptbp1 promoter region, as a functional underpinning of hepatoprotection against IR stress and tissue damage in liver transplantation.

Alternative splicing in neurodegenerative disease and the promise of RNA therapies

Alternative splicing generates a myriad of RNA products and protein isoforms of different functions from a single gene. Dysregulated alternative splicing has emerged as a new mechanism broadly implicated in the pathogenesis of neurodegenerative diseases such as Alzheimer disease, amyotrophic lateral sclerosis, frontotemporal dementia, Parkinson disease and repeat expansion diseases. Understanding the mechanisms and functional outcomes of abnormal splicing in neurological disorders is vital in developing effective therapies to treat mis-splicing pathology. In this Review, we discuss emerging research and evidence of the roles of alternative splicing defects in major neurodegenerative diseases and summarize the latest advances in RNA-based therapeutic strategies to target these disorders.

Characterising the RNA-binding protein atlas of the mammalian brain uncovers RBM5 misregulation in mouse models of Huntington's disease

RNA-binding proteins (RBPs) are key players regulating RNA processing and are associated with disorders ranging from cancer to neurodegeneration. Here, we present a proteomics workflow for large-scale identification of RBPs and their RNA-binding regions in the mammalian brain identifying 526 RBPs. Analysing brain tissue from males of the Huntington's disease (HD) R6/2 mouse model uncovered differential RNA-binding of the alternative splicing regulator RBM5. Combining several omics workflows, we show that RBM5 binds differentially to transcripts enriched in pathways of neurodegeneration in R6/2 brain tissue. We further find these transcripts to undergo changes in splicing and demonstrate that RBM5 directly regulates these changes in human neurons derived from embryonic stem cells. Finally, we reveal that RBM5 interacts differently with several known huntingtin interactors and components of huntingtin aggregates. Collectively, we demonstrate the applicability of our method for capturing RNA interactor dynamics in the contexts of tissue and disease.

RNA alternative splicing impacts the risk for alcohol use disorder

Alcohol use disorder (AUD) is a complex genetic disorder characterized by problems arising from excessive alcohol consumption. Identifying functional genetic variations that contribute to risk for AUD is a major goal. Alternative splicing of RNA mediates the flow of genetic information from DNA to gene expression and expands proteome diversity. We asked whether alternative splicing could be a risk factor for AUD. Herein, we used a Mendelian randomization (MR)-based approach to identify skipped exons (the predominant splicing event in brain) that contribute to AUD risk. Genotypes and RNA-seq data from the CommonMind Consortium were used as the training dataset to develop predictive models linking individual genotypes to exon skipping in the prefrontal cortex. We applied these models to data from the Collaborative Studies on Genetics of Alcoholism to examine the association between the imputed cis-regulated splicing outcome and the AUD-related traits. We identified 27 exon skipping events that were predicted to affect AUD risk; six of these were replicated in the Australian Twin-family Study of Alcohol Use Disorder. Their host genes are DRC1, ELOVL7, LINC00665, NSUN4, SRRM2 and TBC1D5. The genes downstream of these splicing events are enriched in neuroimmune pathways. The MR-inferred impacts of the ELOVL7 skipped exon on AUD risk was further supported in four additional large-scale genome-wide association studies. Additionally, this exon contributed to changes of gray matter volumes in multiple brain regions, including the visual cortex known to be involved in AUD. In conclusion, this study provides strong evidence that RNA alternative splicing impacts the susceptibility to AUD and adds new information on AUD-relevant genes and pathways. Our framework is also applicable to other types of splicing events and to other complex genetic disorders.

Schizophrenia-associated NRXN1 deletions induce developmental-timing- and cell-type-specific vulnerabilities in human brain organoids

De novo mutations and copy number deletions in NRXN1 (2p16.3) pose a significant risk for schizophrenia (SCZ). It is unclear how NRXN1 deletions impact cortical development in a cell type-specific manner and disease background modulates these phenotypes. Here, we leveraged human pluripotent stem cell-derived forebrain organoid models carrying NRXN1 heterozygous deletions in isogenic and SCZ patient genetic backgrounds and conducted single-cell transcriptomic analysis over the course of brain organoid development from 3 weeks to 3.5 months. Intriguingly, while both deletions similarly impacted molecular pathways associated with ubiquitin-proteasome system, alternative splicing, and synaptic signaling in maturing glutamatergic and GABAergic neurons, SCZ-NRXN1 deletions specifically perturbed developmental trajectories of early neural progenitors and accumulated disease-specific transcriptomic signatures. Using calcium imaging, we found that both deletions led to long-lasting changes in spontaneous and synchronous neuronal networks, implicating synaptic dysfunction. Our study reveals developmental-timing- and cell-type-dependent actions of NRXN1 deletions in unique genetic contexts.

Variants in SART3 cause a spliceosomopathy characterised by failure of testis development and neuronal defects

Squamous cell carcinoma antigen recognized by T cells 3 (SART3) is an RNA-binding protein with numerous biological functions including recycling small nuclear RNAs to the spliceosome. Here, we identify recessive variants in SART3 in nine individuals presenting with intellectual disability, global developmental delay and a subset of brain anomalies, together with gonadal dysgenesis in 46,XY individuals. Knockdown of the Drosophila orthologue of SART3 reveals a conserved role in testicular and neuronal development. Human induced pluripotent stem cells carrying patient variants in SART3 show disruption to multiple signalling pathways, upregulation of spliceosome components and demonstrate aberrant gonadal and neuronal differentiation in vitro. Collectively, these findings suggest that bi-allelic SART3 variants underlie a spliceosomopathy which we tentatively propose be termed INDYGON syndrome (Intellectual disability, Neurodevelopmental defects and Developmental delay with 46,XY GONadal dysgenesis). Our findings will enable additional diagnoses and improved outcomes for individuals born with this condition.

Hidden Genetic Regulation of Human Complex Traits via Brain Isoforms

Alternative splicing exists in most multi-exonic genes, and exploring these complex alternative splicing events and their resultant isoform expressions is essential. However, it has become conventional that RNA sequencing results have often been summarized into gene-level expression counts mainly due to the multiple ambiguous mapping of reads at highly similar regions. Transcript-level quantification and interpretation are often overlooked, and biological interpretations are often deduced based on combined transcript information at the gene level. Here, for the most variable tissue of alternative splicing, the brain, we estimate isoform expressions in 1,191 samples collected by the Genotype-Tissue Expression (GTEx) Consortium using a powerful method that we previously developed. We perform genome-wide association scans on the isoform ratios per gene and identify isoform-ratio quantitative trait loci (irQTL), which could not be detected by studying gene-level expressions alone. By analyzing the genetic architecture of the irQTL, we show that isoform ratios regulate educational attainment via multiple tissues including the frontal cortex (BA9), cortex, cervical spinal cord, and hippocampus. These tissues are also associated with different neuro-related traits, including Alzheimer's or dementia, mood swings, sleep duration, alcohol intake, intelligence, anxiety or depression, etc. Mendelian randomization (MR) analysis revealed 1,139 pairs of isoforms and neuro-related traits with plausible causal relationships, showing much stronger causal effects than on general diseases measured in the UK Biobank (UKB). Our results highlight essential transcript-level biomarkers in the human brain for neuro-related complex traits and diseases, which could be missed by merely investigating overall gene expressions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s43657-023-00100-6.

Functional Characterisation of the Circular RNA, circHTT(2-6), in Huntington's Disease

Trinucleotide repeat disorders comprise ~20 severe, inherited, human neuromuscular and neurodegenerative disorders, which result from an abnormal expansion of repetitive sequences in the DNA. The most common of these, Huntington's disease (HD), results from expansion of the CAG repeat region in exon 1 of the HTT gene via an unknown mechanism. Since non-coding RNAs have been implicated in the initiation and progression of many diseases, herein we focused on a circular RNA (circRNA) molecule arising from non-canonical splicing (backsplicing) of HTT pre-mRNA. The most abundant circRNA from HTT, circHTT(2-6), was found to be more highly expressed in the frontal cortex of HD patients, compared with healthy controls, and positively correlated with CAG repeat tract length. Furthermore, the mouse orthologue (mmu_circHTT(2-6)) was found to be enriched within the brain and specifically the striatum, a region enriched for medium spiny neurons that are preferentially lost in HD. Transgenic overexpression of circHTT(2-6) in two human cell lines-SH-SY5Y and HEK293-reduced cell proliferation and nuclear size without affecting cell cycle progression or cellular size, or altering the CAG repeat region length within HTT. CircHTT(2-6) overexpression did not alter total HTT protein levels, but reduced its nuclear localisation. As these phenotypic and genotypic changes resemble those observed in HD patients, our results suggest that circHTT(2-6) may play a functional role in the pathophysiology of this disease.

Alternative splicing shapes the transcriptome complexity in blackgram [Vigna mungo (L.) Hepper]

Vigna mungo, a highly consumed crop in the pan-Asian countries, is vulnerable to several biotic and abiotic stresses. Understanding the post-transcriptional gene regulatory cascades, especially alternative splicing (AS), may underpin large-scale genetic improvements to develop stress-resilient varieties. Herein, a transcriptome based approach was undertaken to decipher the genome-wide AS landscape and splicing dynamics in order to establish the intricacies of their functional interactions in various tissues and stresses. RNA sequencing followed by high-throughput computational analyses identified 54,526 AS events involving 15,506 AS genes that generated 57,405 transcripts isoforms. Enrichment analysis revealed their involvement in diverse regulatory functions and demonstrated that transcription factors are splicing-intensive, splice variants of which are expressed differentially across tissues and environmental cues. Increased expression of a splicing regulator NHP2L1/SNU13 was found to co-occur with lower intron retention events. The host transcriptome is significantly impacted by differential isoform expression of 1172 and 765 AS genes that resulted in 1227 (46.8% up and 53.2% downregulated) and 831 (47.5% up and 52.5% downregulated) transcript isoforms under viral pathogenesis and Fe²⁺ stressed condition, respectively. However, genes experiencing AS operate differently from the differentially expressed genes, suggesting AS is a unique and independent mode of regulatory mechanism. Therefore, it can be inferred that AS mediates a crucial regulatory role across tissues and stressful situations and the results would provide an invaluable resource for future endeavours in V. mungo genomics.

Cell-type-specific alternative splicing in the Arabidopsis germline

During sexual reproduction in flowering plants, the two haploid sperm cells (SCs) embedded within the cytoplasm of a growing pollen tube are carried to the embryo sac for double fertilization. Pollen development in flowering plants is a dynamic process that encompasses changes at transcriptome and epigenome levels. While the transcriptome of pollen and SCs in Arabidopsis (Arabidopsis thaliana) is well documented, previous analyses have mostly been based on gene-level expression. In-depth transcriptome analysis, particularly the extent of alternative splicing (AS) at the resolution of SC and vegetative nucleus (VN), is still lacking. Therefore, we performed RNA-seq analysis to generate a spliceome map of Arabidopsis SCs and VN isolated from mature pollen grains. Based on our de novo transcriptome assembly, we identified 58,039 transcripts, including 9,681 novel transcripts, of which 2,091 were expressed in SCs and 3,600 in VN. Four hundred and sixty-eight genes were regulated both at gene and splicing levels, with many having functions in mRNA splicing, chromatin modification, and protein localization. Moreover, a comparison with egg cell RNA-seq data uncovered sex-specific regulation of transcription and splicing factors. Our study provides insights into a gamete-specific AS landscape at unprecedented resolution.

ACTOR: a latent Dirichlet model to compare expressed isoform proportions to a reference panel

The relative proportion of RNA isoforms expressed for a given gene has been associated with disease states in cancer, retinal diseases, and neurological disorders. Examination of relative isoform proportions can help determine biological mechanisms, but such analyses often require a per-gene investigation of splicing patterns. Leveraging large public data sets produced by genomic consortia as a reference, one can compare splicing patterns in a data set of interest with those of a reference panel in which samples are divided into distinct groups, such as tissue of origin, or disease status. We propose A latent Dirichlet model to Compare expressed isoform proportions TO a Reference panel (ACTOR), a latent Dirichlet model with Dirichlet Multinomial observations to compare expressed isoform proportions in a data set to an independent reference panel. We use a variational Bayes procedure to estimate posterior distributions for the group membership of one or more samples. Using the Genotype-Tissue Expression project as a reference data set, we evaluate ACTOR on simulated and real RNA-seq data sets to determine tissue-type classifications of genes. ACTOR is publicly available as an R package at https://github.com/mccabes292/actor.

Alternative splicing is not a key source of chemerin isoforms diversity

BACKGROUND

Chemerin is a chemoattractant protein with adipokine and antimicrobial properties encoded by the retinoic acid receptor responder 2 (RARRES2) gene. Chemerin bioactivity largely depends on carboxyl-terminal proteolytic processing that generates chemerin isoforms with different chemotactic, regulatory, and antimicrobial potentials. While these mechanisms are relatively well known, the role of alternative splicing in generating isoform diversity remains obscure.

METHODS AND RESULTS

Using rapid amplification of cDNA ends (RACE) PCR, we determined RARRES2 transcript variants present in mouse and human tissues and identified novel transcript variant 4 of mouse Rarres2 encoding mChem153K. Moreover, analyses of real-time quantitative PCR (RT-qPCR) and publicly-available next-generation RNA sequencing (RNA-seq) datasets showed that different alternatively spliced variants of mouse Rarres2 are present in mouse tissues and their expression patterns were unaffected by inflammatory and infectious stimuli except brown adipose tissue. However, only one transcript variant of human RARRES2 was present in liver and adipose tissue.

CONCLUSION

Our findings indicate a limited role for alternative splicing in generating chemerin isoform diversity under all tested conditions.

Tau-induced deficits in nonsense-mediated mRNA decay contribute to neurodegeneration

INTRODUCTION

While brains of patients with Alzheimer's disease and related tauopathies have evidence of altered RNA processing, we lack a mechanistic understanding of how altered RNA processing arises in these disorders and if such changes are causally linked to neurodegeneration.

METHODS

Using Drosophila melanogaster models of tauopathy, we find that overall activity of nonsense-mediated mRNA decay (NMD), a key RNA quality-control mechanism, is reduced. Genetic manipulation of NMD machinery significantly modifies tau-induced neurotoxicity, suggesting that deficits in NMD are causally linked to neurodegeneration. Mechanistically, we find that deficits in NMD are a consequence of aberrant RNA export and RNA accumulation within nuclear envelope invaginations in tauopathy. We identify a pharmacological activator of NMD that suppresses neurodegeneration in tau transgenic Drosophila, indicating that tau-induced deficits in RNA quality control are druggable.

DISCUSSION

Our studies suggest that NMD activators should be explored for their potential therapeutic value to patients with tauopathies.

Identification of Alternative Splicing in Proteomes of Human Melanoma Cell Lines without RNA Sequencing Data

Alternative splicing is one of the main regulation pathways in living cells beyond simple changes in the level of protein expression. Most of the approaches proposed in proteomics for the identification of specific splicing isoforms require a preliminary deep transcriptomic analysis of the sample under study, which is not always available, especially in the case of the re-analysis of previously acquired data. Herein, we developed new algorithms for the identification and validation of protein splice isoforms in proteomic data in the absence of RNA sequencing of the samples under study. The bioinformatic approaches were tested on the results of proteome analysis of human melanoma cell lines, obtained earlier by high-resolution liquid chromatography and mass spectrometry (LC-MS). A search for alternative splicing events for each of the cell lines studied was performed against the database generated from all known transcripts (RefSeq) and the one composed of peptide sequences, which included all biologically possible combinations of exons. The identifications were filtered using the prediction of both retention times and relative intensities of fragment ions in the corresponding mass spectra. The fragmentation mass spectra corresponding to the discovered alternative splicing events were additionally examined for artifacts. Selected splicing events were further validated at the mRNA level by quantitative PCR.

Specific expression of alternatively spliced genes in the turkey (Meleagris gallopavo) reproductive tract revealed their function in spermatogenesis and post-testicular sperm maturation

The tissue-specific profile of alternatively spliced genes (ASGs) and their involvement in reproduction processes characteristic of turkey testis, epididymis, and ductus deferens were investigated for the first time in birds. Deep sequencing of male turkey reproductive tissue RNA samples (n = 6) was performed using Illumina RNA-Seq with 2 independent methods, rMATs and SUPPA2, for differential alternative splicing (DAS) event prediction. The expression of selected ASGs was validated using quantitative real-time reverse transcriptase-polymerase chain reaction. The testis was found to be the site of the highest number of posttranscriptional splicing events within the reproductive tract, and skipping exons were the most frequently occurring class of alternative splicing (AS) among the reproductive tract. Statistical analysis revealed 86, 229, and 6 DAS events in the testis/epididymis, testis/ductus deferens, and epididymis/ductus deferens comparison, respectively. Alternative splicing was found to be a mechanism of gene expression regulation within the turkey reproduction tract. In testis, modification was observed for spermatogenesis specific genes; the changes in 5' UTR could act as regulator of MEIG1 expression (a player during spermatocytes meiosis), and modification of 3' UTR led to diversification of CREM mRNA (modulator of gene expression related to the structuring of mature spermatozoa). Sperm tail formation can be regulated by changes in the 5' UTR of testicular SLC9A3R1 and gene silencing by producing dysfunctional variants of ODF2 in the testis and ATP1B3 in the epididymis. Predicted differentially ASGs in the turkey reproductive tract seem to be involved in the regulation of spermatogenesis, including acrosome formation and sperm tail formation and binding of sperm to the zona pellucida. Several ASGs were classified as cilia by actin and microtubule cytoskeleton organization. Such genes may play a role in the organization of sperm flagellum and post-testicular motility development. To our knowledge, this is the first functional investigation of alternatively spliced genes associated with tissue-specific processes in the turkey reproductive tract.

MAEL gene contributes to bovine testicular development through the m5C-mediated splicing

Knowledge of RNA molecules regulating testicular development and spermatogenesis in bulls is essential for elite bull selection and an ideal breeding program. Herein, we performed direct RNA sequencing (DRS) to explore the functional characterization of RNA molecules produced in the testicles of 9 healthy Simmental bulls at three testicular development stages (prepuberty, puberty, and postpuberty). We identified 5,043 differentially expressed genes associated with testicular weight. These genes exhibited more alternative splicing at sexual maturity, particularly alternative 3' (A3) and 5' (A5) splice sites usage and exon skipping (SE). The expression of hub genes in testicular developmental stages was also affected by both m6A and m5C RNA modifications. We found m5C-mediated splicing events significantly (p < 0.05) increased MAEL gene expression at the isoform level, likely promoting spermatogenesis. Our findings highlight the complexity of RNA processing and expression as well as the regulation of transcripts involved in testicular development and spermatogenesis.

LAFITE Reveals the Complexity of Transcript Isoforms in Subcellular Fractions

Characterization of the subcellular distribution of RNA is essential for understanding the molecular basis of biological processes. Here, the subcellular nanopore direct RNA-sequencing (DRS) of four lung cancer cell lines (A549, H1975, H358, and HCC4006) is performed, coupled with a computational pipeline, Low-abundance Aware Full-length Isoform clusTEr (LAFITE), to comprehensively analyze the full-length cytoplasmic and nuclear transcriptome. Using additional DRS and orthogonal data sets, it is shown that LAFITE outperforms current methods for detecting full-length transcripts, particularly for low-abundance isoforms that are usually overlooked due to poor read coverage. Experimental validation of six novel isoforms exclusively identified by LAFITE further confirms the reliability of this pipeline. By applying LAFITE to subcellular DRS data, the complexity of the nuclear transcriptome is revealed in terms of isoform diversity, 3'-UTR usage, m6A modification patterns, and intron retention. Overall, LAFITE provides enhanced full-length isoform identification and enables a high-resolution view of the RNA landscape at the isoform level.

Posttranscriptional regulation of neurofilament proteins and tau in health and disease

Neurofilament and tau proteins are neuron-specific cytoskeletal proteins that are enriched in axons, regulated by many of the same protein kinases, interact physically, and are the principal constituents of neurofibrillary lesions in major adult-onset dementias. Both proteins share functions related to the modulation of stability and functions of the microtubule network in axons, axonal transport and scaffolding of organelles, long-term synaptic potentiation, and learning and memory. Expression of these proteins is regulated not only at the transcriptional level but also through posttranscriptional control of pre-mRNA splicing, mRNA stability, transport, localization, local translation and degradation. Current evidence suggests that posttranscriptional determinants of their levels are usually regulated by RNA-binding proteins and microRNAs primarily through 3'-untranslated regions of neurofilament and tau mRNAs. Dysregulations of neurofilament and tau expression caused by mutations or pathologies of RNA-binding proteins such as TDP43, FUS and microRNAs are increasingly recognized in association with varied neurological disorders. In this review, we summarize the current understanding of posttranscriptional control of neurofilament and tau by examining the posttranscriptional regulation of neurofilament and tau by RNA-binding proteins and microRNAs implicated in health and diseases.

The alternative proteome in neurobiology

Translation involves the biosynthesis of a protein sequence following the decoding of the genetic information embedded in a messenger RNA (mRNA). Typically, the eukaryotic mRNA was considered to be inherently monocistronic, but this paradigm is not in agreement with the translational landscape of cells, tissues, and organs. Recent ribosome sequencing (Ribo-seq) and proteomics studies show that, in addition to currently annotated reference proteins (RefProt), other proteins termed alternative proteins (AltProts), and microproteins are encoded in regions of mRNAs thought to be untranslated or in transcripts annotated as non-coding. This experimental evidence expands the repertoire of functional proteins within a cell and potentially provides important information on biological processes. This review explores the hitherto overlooked alternative proteome in neurobiology and considers the role of AltProts in pathological and healthy neuromolecular processes.

SRSF10 is essential for progenitor spermatogonia expansion by regulating alternative splicing

Alternative splicing expands the transcriptome and proteome complexity and plays essential roles in tissue development and human diseases. However, how alternative splicing regulates spermatogenesis remains largely unknown. Here, using a germ cell-specific knockout mouse model, we demonstrated that the splicing factor Srsf10 is essential for spermatogenesis and male fertility. In the absence of SRSF10, spermatogonial stem cells can be formed, but the expansion of Promyelocytic Leukemia Zinc Finger (PLZF)-positive undifferentiated progenitors was impaired, followed by the failure of spermatogonia differentiation (marked by KIT expression) and meiosis initiation. This was further evidenced by the decreased expression of progenitor cell markers in bulk RNA-seq, and much less progenitor and differentiating spermatogonia in single-cell RNA-seq data. Notably, SRSF10 directly binds thousands of genes in isolated THY+ spermatogonia, and Srsf10 depletion disturbed the alternative splicing of genes that are preferentially associated with germ cell development, cell cycle, and chromosome segregation, including Nasp, Bclaf1, Rif1, Dazl, Kit, Ret, and Sycp1. These data suggest that SRSF10 is critical for the expansion of undifferentiated progenitors by regulating alternative splicing, expanding our understanding of the mechanism underlying spermatogenesis.

Bioinformatics detection of modulators controlling splicing factor-dependent intron retention in the human brain

Alternative RNA splicing is an important means of genetic control and transcriptome diversity. However, when alternative splicing events are studied independently, coordinated splicing modulated by common factors is often not recognized. As a result, the molecular mechanisms of how splicing regulators promote or repress splice site recognition in a context-dependent manner are not well understood. The functional coupling between multiple gene regulatory layers suggests that splicing is modulated by additional genetic or epigenetic components. Here, we developed a bioinformatics approach to identify causal modulators of splicing activity based on the variation of gene expression in large RNA sequencing datasets. We applied this approach in a neurological context with hundreds of dorsolateral prefrontal cortex samples. Our model is strengthened with the incorporation of genetic variants to impute gene expression in a Mendelian randomization-based approach. We identified novel modulators of the splicing factor SRSF1, including UIMC1 and the long noncoding RNA CBR3-AS1, that function over dozens of SRSF1 intron retention splicing targets. This strategy can be widely used to identify modulators of RNA-binding proteins involved in tissue-specific alternative splicing.

Prediction of transcript isoforms in 19 chicken tissues by Oxford Nanopore long-read sequencing

To identify and annotate transcript isoforms in the chicken genome, we generated Nanopore long-read sequencing data from 68 samples that encompassed 19 diverse tissues collected from experimental adult male and female White Leghorn chickens. More than 23.8 million reads with mean read length of 790 bases and average quality of 18.2 were generated. The annotation and subsequent filtering resulted in the identification of 55,382 transcripts at 40,547 loci with mean length of 1,700 bases. We predicted 30,967 coding transcripts at 19,461 loci, and 16,495 lncRNA transcripts at 15,512 loci. Compared to existing reference annotations, we found ∼52% of annotated transcripts could be partially or fully matched while ∼47% were novel. Seventy percent of novel transcripts were potentially transcribed from lncRNA loci. Based on our annotation, we quantified transcript expression across tissues and found two brain tissues (i.e., cerebellum and cortex) expressed the highest number of transcripts and loci. Furthermore, ∼22% of the transcripts displayed tissue specificity with the reproductive tissues (i.e., testis and ovary) exhibiting the most tissue-specific transcripts. Despite our wide sampling, ∼20% of Ensembl reference loci were not detected. This suggests that deeper sequencing and additional samples that include different breeds, cell types, developmental stages, and physiological conditions, are needed to fully annotate the chicken genome. The application of Nanopore sequencing in this study demonstrates the usefulness of long-read data in discovering additional novel loci (e.g., lncRNA loci) and resolving complex transcripts (e.g., the longest transcript for the TTN locus).

Proteomic Analysis of Adult Human Hippocampal Subfields Demonstrates Regional Heterogeneity in the Protein Expression

Background: Distinct hippocampal subfields are known to get affected during aging, psychiatric disorders, and various neurological and neurodegenerative conditions. To understand the biological processes associated with each subfield, it is important to understand its heterogeneity at the molecular level. To address this lacuna, we investigated the proteomic analysis of hippocampal subfields─the cornu ammonis sectors (CA1, CA2, CA3, CA4) and dentate gyrus (DG) from healthy adult human cohorts. Findings: Microdissection of hippocampal subfields from archived formalin-fixed paraffin-embedded tissue sections followed by TMT-based multiplexed proteomic analysis resulted in the identification of 5,593 proteins. Out of these, 890 proteins were found to be differentially abundant among the subfields. Further bioinformatics analysis suggested proteins related to gene splicing, transportation, myelination, structural activity, and learning processes to be differentially abundant in DG, CA4, CA3, CA2, and CA1, respectively. A subset of proteins was selected for immunohistochemistry-based validation in an independent set of hippocampal samples. Conclusions: We believe that our findings will effectively pave the way for further analysis of the hippocampal subdivisions and provide awareness of its subfield-specific association to various neurofunctional anomalies in the future. The current mass spectrometry data is deposited and publicly made available through ProteomeXchange Consortium via the PRIDE partner repository with the data set identifier PXD029697.

Transcriptome variation in human tissues revealed by long-read sequencing

Regulation of transcript structure generates transcript diversity and plays an important role in human disease1-7. The advent of long-read sequencing technologies offers the opportunity to study the role of genetic variation in transcript structure8-16. In this Article, we present a large human long-read RNA-seq dataset using the Oxford Nanopore Technologies platform from 88 samples from Genotype-Tissue Expression (GTEx) tissues and cell lines, complementing the GTEx resource. We identified just over 70,000 novel transcripts for annotated genes, and validated the protein expression of 10% of novel transcripts. We developed a new computational package, LORALS, to analyse the genetic effects of rare and common variants on the transcriptome by allele-specific analysis of long reads. We characterized allele-specific expression and transcript structure events, providing new insights into the specific transcript alterations caused by common and rare genetic variants and highlighting the resolution gained from long-read data. We were able to perturb the transcript structure upon knockdown of PTBP1, an RNA binding protein that mediates splicing, thereby finding genetic regulatory effects that are modified by the cellular environment. Finally, we used this dataset to enhance variant interpretation and study rare variants leading to aberrant splicing patterns.

Selection of Olduvai Domains during Evolution: A Role for Primate-Specific Splicing Super-Enhancer and RNA Guanine Quadruplex in Bipartite NBPF Exons

Olduvai protein domains (also known as DUF1220 or NBPF) have undergone the greatest human-specific increase in the copy number of any coding region in the genome. Their repeat number was strongly associated with the evolutionary expansion of brain volumes, neuron counts and cognitive abilities, as well as with disorders of the autistic spectrum. Nevertheless, the domain function and cellular mechanisms underlying the positive selection of Olduvai DNA sequences in higher primates remain obscure. Here, I show that the inclusion of Olduvai exon doublets in mature transcripts is facilitated by a potent splicing enhancer that was created through duplication within the first exon. The enhancer is the strongest among the NBPF transcripts and further promotes the already high splicing activity of the unexpanded first exons of the two-exon domains, safeguarding the expanded Olduvai exon doublets in the mature transcriptome. The duplication also creates a predicted RNA guanine quadruplex that may regulate the access to spliceosomal components of the super-enhancer and influence the splicing of adjacent exons. Thus, positive Olduvai selection during primate evolution is likely to result from a combination of multiple targets in gene expression pathways, including RNA splicing.

Derivation of splice junction-specific antibodies using a unique hapten targeting strategy and directed evolution

Alternative splicing of RNA occurs frequently in eukaryotic cells and can result in multiple protein isoforms that are nearly identical in amino acid sequence, but have unique biological roles. Moreover, the relative abundance of these unique isoforms can be correlative with diseased states and potentially used as biomarkers or therapeutic targets. However, due to high sequence similarities among isoforms, current proteomic methods are incapable of differentiating native protein isoforms derived from most alternative splicing events. Herein, a strategy employing a nonsynonymous, non-native amino acid (nnAA) pseudo-hapten (i.e. an amino acid or amino acid derivative that is different from the native amino acid at a particular position) as a targeting epitope in splice junction-spanning peptides was successful in directed antibody derivation. After isolating nnAA-specific antibodies, directed evolution reduced the antibody's binding dependence on the nnAA pseudo-hapten and improved binding to the native splice junction epitope. The resulting antibodies demonstrated codependent binding affinity to each exon of the splice junction and thus are splice junction- and isoform-specific. Furthermore, epitope scanning demonstrated that positioning of the nnAA pseudo-hapten within a peptide antigen can be exploited to predetermine the isolated antibody's specificity at, or near, amino acid resolution. Thus, this nnAA targeting strategy has the potential to robustly derive splice junction- and site-specific antibodies that can be used in a wide variety of research endeavors to unambiguously differentiate native protein isoforms.

The Long Non-Coding RNA GOMAFU in Schizophrenia: Function, Disease Risk, and Beyond

Neuropsychiatric diseases are among the most common brain developmental disorders, represented by schizophrenia (SZ). The complex multifactorial etiology of SZ remains poorly understood, which reflects genetic vulnerabilities and environmental risks that affect numerous genes and biological pathways. Besides the dysregulation of protein-coding genes, recent discoveries demonstrate that abnormalities associated with non-coding RNAs, including microRNAs and long non-coding RNAs (lncRNAs), also contribute to the pathogenesis of SZ. lncRNAs are an actively evolving family of non-coding RNAs that harbor greater than 200 nucleotides but do not encode for proteins. In general, lncRNA genes are poorly conserved. The large number of lncRNAs specifically expressed in the human brain, together with the genetic alterations and dysregulation of lncRNA genes in the SZ brain, suggests a critical role in normal cognitive function and the pathogenesis of neuropsychiatric diseases. A particular lncRNA of interest is GOMAFU, also known as MIAT and RNCR2. Growing evidence suggests the function of GOMAFU in governing neuronal development and its potential roles as a risk factor and biomarker for SZ, which will be reviewed in this article. Moreover, we discuss the potential mechanisms through which GOMAFU regulates molecular pathways, including its subcellular localization and interaction with RNA-binding proteins, and how interruption to GOMAFU pathways may contribute to the pathogenesis of SZ.

Exonic splicing code and protein binding sites for calcium

Auxilliary splicing sequences in exons, known as enhancers (ESEs) and silencers (ESSs), have been subject to strong selection pressures at the RNA and protein level. The protein component of this splicing code is substantial, recently estimated at ∼50% of the total information within ESEs, but remains poorly understood. The ESE/ESS profiles were previously associated with the Irving-Williams (I-W) stability series for divalent metals, suggesting that the ESE/ESS evolution was shaped by metal binding sites. Here, we have examined splicing activities of exonic sequences that encode protein binding sites for Ca2+, a weak binder in the I-W affinity order. We found that predicted exon inclusion levels for the EF-hand motifs and for Ca2+-binding residues in nonEF-hand proteins were higher than for average exons. For canonical EF-hands, the increase was centred on the EF-hand chelation loop and, in particular, on Ca2+-coordinating residues, with a 1>12>3∼5>9 hierarchy in the 12-codon loop consensus and usage bias at codons 1 and 12. The same hierarchy but a lower increase was observed for noncanonical EF-hands, except for S100 proteins. EF-hand loops preferentially accumulated exon splits in two clusters, one located in their N-terminal halves and the other around codon 12. Using splicing assays and published crosslinking and immunoprecipitation data, we identify candidate trans-acting factors that preferentially bind conserved GA-rich motifs encoding negatively charged amino acids in the loops. Together, these data provide evidence for the high capacity of codons for Ca2+-coordinating residues to be retained in mature transcripts, facilitating their exon-level expansion during eukaryotic evolution.

Hepatitis, testicular degeneration, and ataxia in DIDO3-deficient mice with altered mRNA processing

Background

mRNA processing is an essential step of gene expression; its malfunction can lead to different degrees of physiological disorder from subclinical disease to death. We previously identified Dido1 as a stemness marker and a gene involved in embryonic stem cell differentiation. DIDO3, the largest protein encoded by the Dido1 gene, is necessary for accurate mRNA splicing and correct transcription termination. The deletion of Dido1 exon16, which encodes the carboxy-terminal half of DIDO3, results in early embryonic lethality in mouse.

Results

We obtained mice bearing a Cre-LoxP conditional version of that deletion and studied the effects of inducing it ubiquitously in adult stages. DIDO3-deficient mice survive the deletion but suffer mild hepatitis, testicular degeneration, and progressive ataxia, in association with systemic alterations in mRNA splicing and transcriptional readthrough.

Conclusions

These results offer insight into the distinct vulnerabilities in mouse organs following impairment of the mRNA processing machinery, and could aid understanding of human health dependence on accurate mRNA metabolism.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-022-00804-8.

Alternative splicing diversifies the skeletal muscle transcriptome during prolonged spaceflight

BACKGROUND

As the interest in manned spaceflight increases, so does the requirement to understand the transcriptomic mechanisms that underlay the detrimental physiological adaptations of skeletal muscle to microgravity. While microgravity-induced differential gene expression (DGE) has been extensively investigated, the contribution of differential alternative splicing (DAS) to the plasticity and functional status of the skeletal muscle transcriptome has not been studied in an animal model. Therefore, by evaluating both DGE and DAS across spaceflight, we set out to provide the first comprehensive characterization of the transcriptomic landscape of skeletal muscle during exposure to microgravity.

METHODS

RNA-sequencing, immunohistochemistry, and morphological analyses were conducted utilizing total RNA and tissue sections isolated from the gastrocnemius and quadriceps muscles of 30-week-old female BALB/c mice exposed to microgravity or ground control conditions for 9 weeks.

RESULTS

In response to microgravity, the skeletal muscle transcriptome was remodeled via both DGE and DAS. Importantly, while DGE showed variable gene network enrichment, DAS was enriched in structural and functional gene networks of skeletal muscle, resulting in the expression of alternatively spliced transcript isoforms that have been associated with the physiological changes to skeletal muscle in microgravity, including muscle atrophy and altered fiber type function. Finally, RNA-binding proteins, which are required for regulation of pre-mRNA splicing, were themselves differentially spliced but not differentially expressed, an upstream event that is speculated to account for the downstream splicing changes identified in target skeletal muscle genes.

CONCLUSIONS

Our work serves as the first investigation of coordinate changes in DGE and DAS in large limb muscles across spaceflight. It opens up a new opportunity to understand (i) the molecular mechanisms by which splice variants of skeletal muscle genes regulate the physiological adaptations of skeletal muscle to microgravity and (ii) how small molecule splicing regulator therapies might thwart muscle atrophy and alterations to fiber type function during prolonged spaceflight.

Human APOER2 Isoforms Have Differential Cleavage Events and Synaptic Properties

Human apolipoprotein E receptor 2 (APOER2) is a type I transmembrane protein with a large extracellular domain (ECD) and a short cytoplasmic tail. APOER2-ECD contains several ligand-binding domains (LBDs) that are organized into exons with aligning phase junctions, which allows for in-frame exon cassette splicing events. We have identified 25 human APOER2 isoforms from cerebral cortex using gene-specific APOER2 primers, where the majority are exon-skipping events within the N-terminal LBD regions compared with six identified in the heart. APOER2 undergoes proteolytic cleavage in response to ligand binding that releases a C-terminal fragment (CTF) and transcriptionally active intracellular domain (ICD). We tested whether the diversity of human brain-specific APOER2 variants affects APOER2 cleavage. We found isoforms with differing numbers of ligand-binding repeats generated different amounts of CTFs compared with full-length APOER2 (APOER2-FL). Specifically, APOER2 isoforms lacking exons 5-8 (Δex5-8) and lacking exons 4-6 (Δex4-6) generated the highest and lowest amounts of CTF generation, respectively, in response to APOE peptide compared with APOER2-FL. The differential CTF generation of Δex5-8 and Δex4-6 coincides with the proteolytic release of the ICD, which mediates transcriptional activation facilitated by the Mint1 adaptor protein. Functionally, we demonstrated loss of mouse Apoer2 decreased miniature event frequency in excitatory synapses, which may be because of a decrease in the total number of synapses and/or VAMP2 positive neurons. Lentiviral infection with human APOER2-FL or Δex4-6 isoform in Apoer2 knockout neurons restored the miniature event frequency but not Δex5-8 isoform. These results suggest that human APOER2 isoforms have differential cleavage events and synaptic properties.SIGNIFICANCE STATEMENT Humans and mice share virtually the same number of protein-coding genes. However, humans have greater complexity of any higher eukaryotic organisms by encoding multiple protein forms through alternative splicing modifications. Alternative splicing allows pre-mRNAs transcribed from genes to be spliced in different arrangements, producing structurally and functionally distinct protein variants that increase proteomic diversity and are particularly prevalent in the human brain. Here, we identified 25 distinct human APOER2 splice variants from the cerebral cortex using gene-specific APOER2 primers, where the majority are exon-skipping events that exclude N-terminal ligand-binding regions of APOER2. We show that some of the APOER2 variants have differential proteolytic properties in response to APOE ligand and exhibit distinct synaptic properties.