Alternative splicing during mammalian organ development

Alternative splicing coupled to nonsense-mediated decay coordinates downregulation of non-neuronal genes in developing mouse neurons

BACKGROUND

The functional coupling between alternative pre-mRNA splicing (AS) and the mRNA quality control mechanism called nonsense-mediated decay (NMD) can modulate transcript abundance. Previous studies have identified several examples of such a regulation in developing neurons. However, the systems-level effects of AS-NMD in this context are poorly understood.

RESULTS

We developed an R package, factR2, which offers a comprehensive suite of AS-NMD analysis functions. Using this tool, we conducted a longitudinal analysis of gene expression in pluripotent stem cells undergoing induced neuronal differentiation. Our analysis uncovers hundreds of AS-NMD events with significant potential to regulate gene expression. Notably, this regulation is significantly overrepresented in specific functional groups of developmentally downregulated genes. Particularly strong association with gene downregulation is detected for alternative cassette exons stimulating NMD upon their inclusion into mature mRNA. By combining bioinformatic analyses with CRISPR/Cas9 genome editing and other experimental approaches we show that NMD-stimulating cassette exons regulated by the RNA-binding protein PTBP1 dampen the expression of their genes in developing neurons. We also provided evidence that the inclusion of NMD-stimulating cassette exons into mature mRNAs is temporally coordinated with NMD-independent gene repression mechanisms.

CONCLUSIONS

Our study provides an accessible workflow for the discovery and prioritization of AS-NMD targets. It further argues that the AS-NMD pathway plays a widespread role in developing neurons by facilitating the downregulation of functionally related non-neuronal genes.

Lineage-specific splicing regulation of MAPT gene in the primate brain

Divergence of precursor messenger RNA (pre-mRNA) alternative splicing (AS) is widespread in mammals, including primates, but the underlying mechanisms and functional impact are poorly understood. Here, we modeled cassette exon inclusion in primate brains as a quantitative trait and identified 1,170 (∼3%) exons with lineage-specific splicing shifts under stabilizing selection. Among them, microtubule-associated protein tau (MAPT) exons 2 and 10 underwent anticorrelated, two-step evolutionary shifts in the catarrhine and hominoid lineages, leading to their present inclusion levels in humans. The developmental-stage-specific divergence of exon 10 splicing, whose dysregulation can cause frontotemporal lobar degeneration (FTLD), is mediated by divergent distal intronic MBNL-binding sites. Competitive binding of these sites by CRISPR-dCas13d/gRNAs effectively reduces exon 10 inclusion, potentially providing a therapeutically compatible approach to modulate tau isoform expression. Our data suggest adaptation of MAPT function and, more generally, a role for AS in the evolutionary expansion of the primate brain.

Single-cell long-read sequencing-based mapping reveals specialized splicing patterns in developing and adult mouse and human brain

RNA isoforms influence cell identity and function. However, a comprehensive brain isoform map was lacking. We analyze single-cell RNA isoforms across brain regions, cell subtypes, developmental time points and species. For 72% of genes, full-length isoform expression varies along one or more axes. Splicing, transcription start and polyadenylation sites vary strongly between cell types, influence protein architecture and associate with disease-linked variation. Additionally, neurotransmitter transport and synapse turnover genes harbor cell-type variability across anatomical regions. Regulation of cell-type-specific splicing is pronounced in the postnatal day 21-to-postnatal day 28 adolescent transition. Developmental isoform regulation is stronger than regional regulation for the same cell type. Cell-type-specific isoform regulation in mice is mostly maintained in the human hippocampus, allowing extrapolation to the human brain. Conversely, the human brain harbors additional cell-type specificity, suggesting gain-of-function isoforms. Together, this detailed single-cell atlas of full-length isoform regulation across development, anatomical regions and species reveals an unappreciated degree of isoform variability across multiple axes.

KIS counteracts PTBP2 and regulates alternative exon usage in neurons

Alternative RNA splicing is an essential and dynamic process in neuronal differentiation and synapse maturation, and dysregulation of this process has been associated with neurodegenerative diseases. Recent studies have revealed the importance of RNA-binding proteins in the regulation of neuronal splicing programs. However, the molecular mechanisms involved in the control of these splicing regulators are still unclear. Here, we show that KIS, a kinase upregulated in the developmental brain, imposes a genome-wide alteration in exon usage during neuronal differentiation in mice. KIS contains a protein-recognition domain common to spliceosomal components and phosphorylates PTBP2, counteracting the role of this splicing factor in exon exclusion. At the molecular level, phosphorylation of unstructured domains within PTBP2 causes its dissociation from two co-regulators, Matrin3 and hnRNPM, and hinders the RNA-binding capability of the complex. Furthermore, KIS and PTBP2 display strong and opposing functional interactions in synaptic spine emergence and maturation. Taken together, our data uncover a post-translational control of splicing regulators that link transcriptional and alternative exon usage programs in neuronal development.

GFAP-isoforms in the nervous system: Understanding the need for diversity

Glial fibrillary acidic protein (GFAP) is an intermediate filament (IF) protein expressed in specific types of glial cells in the nervous system. The expression of GFAP is highly regulated during brain development and in neurological diseases. The presence of distinct GFAP-isoforms in various cell types, developmental stages, and diseases indicates that GFAP (post-)transcriptional regulation has a role in glial cell physiology and pathology. GFAP-isoforms differ in sub-cellular localisation, IF-network assembly properties, and IF-dynamics which results in distinct molecular interactions and mechanical properties of the IF-network. Therefore, GFAP (post-)transcriptional regulation is likely a mechanism by which radial glia, astrocytes, and glioma cells can modulate cellular function.

Antisense oligonucleotide therapeutic approach for Timothy syndrome

Timothy syndrome (TS) is a severe, multisystem disorder characterized by autism, epilepsy, long-QT syndrome and other neuropsychiatric conditions1. TS type 1 (TS1) is caused by a gain-of-function variant in the alternatively spliced and developmentally enriched CACNA1C exon 8A, as opposed to its counterpart exon 8. We previously uncovered several phenotypes in neurons derived from patients with TS1, including delayed channel inactivation, prolonged depolarization-induced calcium rise, impaired interneuron migration, activity-dependent dendrite retraction and an unanticipated persistent expression of exon 8A2-6. We reasoned that switching CACNA1C exon utilization from 8A to 8 would represent a potential therapeutic strategy. Here we developed antisense oligonucleotides (ASOs) to effectively decrease the inclusion of exon 8A in human cells both in vitro and, following transplantation, in vivo. We discovered that the ASO-mediated switch from exon 8A to 8 robustly rescued defects in patient-derived cortical organoids and migration in forebrain assembloids. Leveraging a transplantation platform previously developed7, we found that a single intrathecal ASO administration rescued calcium changes and in vivo dendrite retraction of patient neurons, suggesting that suppression of CACNA1C exon 8A expression is a potential treatment for TS1. Broadly, these experiments illustrate how a multilevel, in vivo and in vitro stem cell model-based approach can identify strategies to reverse disease-relevant neural pathophysiology.

Alternative mRNA splicing events and regulators in epidermal differentiation

Alternative splicing (AS) of messenger RNAs occurs in ∼95% of multi-exon human genes and generates diverse RNA and protein isoforms. We investigated AS events associated with human epidermal differentiation, a process crucial for skin function. We identified 6,413 AS events, primarily involving cassette exons. We also predicted 34 RNA-binding proteins (RBPs) regulating epidermal AS, including 19 previously undescribed candidate regulators. From these results, we identified FUS as an RBP that regulates the balance between keratinocyte proliferation and differentiation. Additionally, we characterized the function of a cassette exon AS event in MAP3K7, which encodes a kinase involved in cell signaling. We found that a switch from the short to long isoform of MAP3K7, triggered during differentiation, enforces the demarcation between proliferating basal progenitors and overlying differentiated strata. Our findings indicate that AS occurs extensively in the human epidermis and has critical roles in skin homeostasis.

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.

Random genetic drift sets an upper limit on mRNA splicing accuracy in metazoans

Most eukaryotic genes undergo alternative splicing (AS), but the overall functional significance of this process remains a controversial issue. It has been noticed that the complexity of organisms (assayed by the number of distinct cell types) correlates positively with their genome-wide AS rate. This has been interpreted as evidence that AS plays an important role in adaptive evolution by increasing the functional repertoires of genomes. However, this observation also fits with a totally opposite interpretation: given that 'complex' organisms tend to have small effective population sizes (Ne), they are expected to be more affected by genetic drift, and hence more prone to accumulate deleterious mutations that decrease splicing accuracy. Thus, according to this 'drift barrier' theory, the elevated AS rate in complex organisms might simply result from a higher splicing error rate. To test this hypothesis, we analyzed 3496 transcriptome sequencing samples to quantify AS in 53 metazoan species spanning a wide range of Ne values. Our results show a negative correlation between Ne proxies and the genome-wide AS rates among species, consistent with the drift barrier hypothesis. This pattern is dominated by low abundance isoforms, which represent the vast majority of the splice variant repertoire. We show that these low abundance isoforms are depleted in functional AS events, and most likely correspond to errors. Conversely, the AS rate of abundant isoforms, which are relatively enriched in functional AS events, tends to be lower in more complex species. All these observations are consistent with the hypothesis that variation in AS rates across metazoans reflects the limits set by drift on the capacity of selection to prevent gene expression errors.

Comprehensive Atlas of Alternative Splicing Reveals NSRP1 Promoting Adipogenesis through CCDC18

Alternative splicing (AS) plays a crucial role in regulating gene expression, function, and diversity. However, limited reports exist on the identification and comparison of AS in Eastern and Western pigs. Here, we analyzed 243 transcriptome data from eight tissues, integrating information on transcription factors (TFs), selection signals, splicing factors (SFs), and quantitative trait loci (QTL) to comprehensively study alternative splicing events (ASEs) in pigs. Five ASE types were identified, with Mutually Exclusive Exon (MXE) and Skipped Exon (SE) ASEs being the most prevalent. A significant portion of genes with ASEs (ASGs) showed conservation across all eight tissues (63.21-76.13% per tissue). Differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) exhibited tissue specificity, with blood and adipose tissues having more DASGs. Functional enrichment analysis revealed coDASG_DEGs in adipose were enriched in pathways associated with adipose deposition and immune inflammation, while coDASG_DEGs in blood were enriched in pathways related to immune inflammation and metabolism. Adipose deposition in Eastern pigs might be linked to the down-regulation of immune-inflammation-related pathways and reduced insulin resistance. The TFs, selection signals, and SFs appeared to regulate ASEs. Notably, ARID4A (TF), NSRP1 (SF), ANKRD12, IFT74, KIAA2026, CCDC18, NEXN, PPIG, and ROCK1 genes in adipose tissue showed potential regulatory effects on adipose-deposition traits. NSRP1 could promote adipogenesis by regulating alternative splicing and expression of CCDC18. Conducting an in-depth investigation into AS, this study has successfully identified key marker genes essential for pig genetic breeding and the enhancement of meat quality, which will play important roles in promoting the diversity of pork quality and meeting market demand.

Metabolic regulation of mRNA splicing

Alternative mRNA splicing enables the diversification of the proteome from a static genome and confers plasticity and adaptiveness on cells. Although this is often explored in development, where hard-wired programs drive the differentiation and specialization, alternative mRNA splicing also offers a way for cells to react to sudden changes in outside stimuli such as small-molecule metabolites. Fluctuations in metabolite levels and availability in particular convey crucial information to which cells react and adapt. We summarize and highlight findings surrounding the metabolic regulation of mRNA splicing. We discuss the principles underlying the biochemistry and biophysical properties of mRNA splicing, and propose how these could intersect with metabolite levels. Further, we present examples in which metabolites directly influence RNA-binding proteins and splicing factors. We also discuss the interplay between alternative mRNA splicing and metabolite-responsive signaling pathways. We hope to inspire future research to obtain a holistic picture of alternative mRNA splicing in response to metabolic cues.

Splicing factor SRSF1 is essential for homing of precursor spermatogonial stem cells in mice

Spermatogonial stem cells (SSCs) are essential for continuous spermatogenesis and male fertility. The underlying mechanisms of alternative splicing (AS) in mouse SSCs are still largely unclear. We demonstrated that SRSF1 is essential for gene expression and splicing in mouse SSCs. Crosslinking immunoprecipitation and sequencing data revealed that spermatogonia-related genes (e.g. Plzf, Id4, Setdb1, Stra8, Tial1/Tiar, Bcas2, Ddx5, Srsf10, Uhrf1, and Bud31) were bound by SRSF1 in the mouse testes. Specific deletion of Srsf1 in mouse germ cells impairs homing of precursor SSCs leading to male infertility. Whole-mount staining data showed the absence of germ cells in the testes of adult conditional knockout (cKO) mice, which indicates Sertoli cell-only syndrome in cKO mice. The expression of spermatogonia-related genes (e.g. Gfra1, Pou5f1, Plzf, Dnd1, Stra8, and Taf4b) was significantly reduced in the testes of cKO mice. Moreover, multiomics analysis suggests that SRSF1 may affect survival of spermatogonia by directly binding and regulating Tial1/Tiar expression through AS. In addition, immunoprecipitation mass spectrometry and co-immunoprecipitation data showed that SRSF1 interacts with RNA splicing-related proteins (e.g. SART1, RBM15, and SRSF10). Collectively, our data reveal the critical role of SRSF1 in spermatogonia survival, which may provide a framework to elucidate the molecular mechanisms of the posttranscriptional network underlying homing of precursor SSCs.

Enhanced bovine genome annotation through integration of transcriptomics and epi-transcriptomics datasets facilitates genomic biology

BACKGROUND

The accurate identification of the functional elements in the bovine genome is a fundamental requirement for high-quality analysis of data informing both genome biology and genomic selection. Functional annotation of the bovine genome was performed to identify a more complete catalog of transcript isoforms across bovine tissues.

RESULTS

A total of 160,820 unique transcripts (50% protein coding) representing 34,882 unique genes (60% protein coding) were identified across tissues. Among them, 118,563 transcripts (73% of the total) were structurally validated by independent datasets (PacBio isoform sequencing data, Oxford Nanopore Technologies sequencing data, de novo assembled transcripts from RNA sequencing data) and comparison with Ensembl and NCBI gene sets. In addition, all transcripts were supported by extensive data from different technologies such as whole transcriptome termini site sequencing, RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression, chromatin immunoprecipitation sequencing, and assay for transposase-accessible chromatin using sequencing. A large proportion of identified transcripts (69%) were unannotated, of which 86% were produced by annotated genes and 14% by unannotated genes. A median of two 5' untranslated regions were expressed per gene. Around 50% of protein-coding genes in each tissue were bifunctional and transcribed both coding and noncoding isoforms. Furthermore, we identified 3,744 genes that functioned as noncoding genes in fetal tissues but as protein-coding genes in adult tissues. Our new bovine genome annotation extended more than 11,000 annotated gene borders compared to Ensembl or NCBI annotations. The resulting bovine transcriptome was integrated with publicly available quantitative trait loci data to study tissue-tissue interconnection involved in different traits and construct the first bovine trait similarity network.

CONCLUSIONS

These validated results show significant improvement over current bovine genome annotations.

Revisiting the development of cerebellar inhibitory interneurons in the light of single-cell genetic analyses

The present review aims to provide a short update of our understanding of the inhibitory interneurons of the cerebellum. While these cells constitute but a minority of all cerebellar neurons, their functional significance is increasingly being recognized. For one, inhibitory interneurons of the cerebellar cortex are now known to constitute a clearly more diverse group than their traditional grouping as stellate, basket, and Golgi cells suggests, and this diversity is now substantiated by single-cell genetic data. The past decade or so has also provided important information about interneurons in cerebellar nuclei. Significantly, developmental studies have revealed that the specification and formation of cerebellar inhibitory interneurons fundamentally differ from, say, the cortical interneurons, and define a mode of diversification critically dependent on spatiotemporally patterned external signals. Last, but not least, in the past years, dysfunction of cerebellar inhibitory interneurons could also be linked with clinically defined deficits. I hope that this review, however fragmentary, may stimulate interest and help focus research towards understanding the cerebellum.

PTBP1 as a potential regulator of disease

Polypyrimidine tract-binding protein 1 (PTBP1) is a member of the heterogeneous nuclear ribonucleoprotein (hnRNP) family, which plays a key role in alternative splicing of precursor mRNA and RNA metabolism. PTBP1 is universally expressed in various tissues and binds to multiple downstream transcripts to interfere with physiological and pathological processes such as the tumor growth, body metabolism, cardiovascular homeostasis, and central nervous system damage, showing great prospects in many fields. The function of PTBP1 involves the regulation and interaction of various upstream molecules, including circular RNAs (circRNAs), microRNAs (miRNAs) and long non-coding RNAs (lncRNAs). These regulatory systems are inseparable from the development and treatment of diseases. Here, we review the latest knowledge regarding the structure and molecular functions of PTBP1 and summarize its functions and mechanisms of PTBP1 in various diseases, including controversial studies. Furthermore, we recommend future studies on PTBP1 and discuss the prospects of targeting PTBP1 in new clinical therapeutic approaches.

Analysis of alternative splicing in chicken macrophages transfected with overexpression/knockdown of RIP2 gene

Receptor-interacting protein 2 (RIP2) plays a critical role in the transduction of many signaling pathways and is associated with many diseases. Alternative splicing (AS) is an essential and ubiquitous regulatory mechanism of gene expression that contributes to distinct transcript variants and many different kinds of proteins. In this present study, we characterized genome-wide AS events in wild-type chicken macrophages (WT) and RIP2 overexpression/knockdown chicken macrophages (oeRIP2/shRIP2) by high-throughput RNA sequencing technology. A total of 1901, 2061, and 817 differentially expressed (DE) AS genes were identified in the comparison of oeRIP2 vs. WT, oeRIP2 vs. shRIP2, and shRIP2 vs. WT, respectively. These DE AS genes participated in many important KEGG pathways, including regulation of autophagy, Wnt signaling pathway, Ubiquitin mediated proteolysis, MAPK signaling pathway, and Focal adhesion, etc. In conclusion, this research provided a broad atlas of the genome-wide scale of the AS event landscape in RIP2 overexpression/knockdown and wild-type chicken macrophages. This research also provides the theoretical basis of the gene network related to RIP2.

Transcript Isoform Diversity of Ampliconic Genes on the Y Chromosome of Great Apes

Y chromosomal ampliconic genes (YAGs) are important for male fertility, as they encode proteins functioning in spermatogenesis. The variation in copy number and expression levels of these multicopy gene families has been studied in great apes; however, the diversity of splicing variants remains unexplored. Here, we deciphered the sequences of polyadenylated transcripts of all nine YAG families (BPY2, CDY, DAZ, HSFY, PRY, RBMY, TSPY, VCY, and XKRY) from testis samples of six great ape species (human, chimpanzee, bonobo, gorilla, Bornean orangutan, and Sumatran orangutan). To achieve this, we enriched YAG transcripts with capture probe hybridization and sequenced them with long (Pacific Biosciences) reads. Our analysis of this data set resulted in several findings. First, we observed evolutionarily conserved alternative splicing patterns for most YAG families except for BPY2 and PRY. Second, our results suggest that BPY2 transcripts and proteins originate from separate genomic regions in bonobo versus human, which is possibly facilitated by acquiring new promoters. Third, our analysis indicates that the PRY gene family, having the highest representation of noncoding transcripts, has been undergoing pseudogenization. Fourth, we have not detected signatures of selection in the five YAG families shared among great apes, even though we identified many species-specific protein-coding transcripts. Fifth, we predicted consensus disorder regions across most gene families and species, which could be used for future investigations of male infertility. Overall, our work illuminates the YAG isoform landscape and provides a genomic resource for future functional studies focusing on infertility phenotypes in humans and critically endangered great apes.

eQTL mapping in fetal-like pancreatic progenitor cells reveals early developmental insights into diabetes risk

The impact of genetic regulatory variation active in early pancreatic development on adult pancreatic disease and traits is not well understood. Here, we generate a panel of 107 fetal-like iPSC-derived pancreatic progenitor cells (iPSC-PPCs) from whole genome-sequenced individuals and identify 4065 genes and 4016 isoforms whose expression and/or alternative splicing are affected by regulatory variation. We integrate eQTLs identified in adult islets and whole pancreas samples, which reveal 1805 eQTL associations that are unique to the fetal-like iPSC-PPCs and 1043 eQTLs that exhibit regulatory plasticity across the fetal-like and adult pancreas tissues. Colocalization with GWAS risk loci for pancreatic diseases and traits show that some putative causal regulatory variants are active only in the fetal-like iPSC-PPCs and likely influence disease by modulating expression of disease-associated genes in early development, while others with regulatory plasticity likely exert their effects in both the fetal and adult pancreas by modulating expression of different disease genes in the two developmental stages.

Isoform-resolved transcriptome of the human preimplantation embryo

Human preimplantation development involves extensive remodeling of RNA expression and splicing. However, its transcriptome has been compiled using short-read sequencing data, which fails to capture most full-length mRNAs. Here, we generate an isoform-resolved transcriptome of early human development by performing long- and short-read RNA sequencing on 73 embryos spanning the zygote to blastocyst stages. We identify 110,212 unannotated isoforms transcribed from known genes, including highly conserved protein-coding loci and key developmental regulators. We further identify 17,964 isoforms from 5,239 unannotated genes, which are largely non-coding, primate-specific, and highly associated with transposable elements. These isoforms are widely supported by the integration of published multi-omics datasets, including single-cell 8CLC and blastoid studies. Alternative splicing and gene co-expression network analyses further reveal that embryonic genome activation is associated with splicing disruption and transient upregulation of gene modules. Together, these findings show that the human embryo transcriptome is far more complex than currently known, and will act as a valuable resource to empower future studies exploring development.

The role of alternative pre-mRNA splicing in cancer progression

Alternative pre-mRNA splicing is a critical mechanism that generates multiple mRNA from a single gene, thereby increasing the diversity of the proteome. Recent research has highlighted the significance of specific splicing isoforms in cellular processes, particularly in regulating cell numbers. In this review, we examine the current understanding of the role of alternative splicing in controlling cancer cell growth and discuss specific splicing factors and isoforms and their molecular mechanisms in cancer progression. These isoforms have been found to intricately control signaling pathways crucial for cell cycle progression, proliferation, and apoptosis. Furthermore, studies have elucidated the characteristics and functional importance of splicing factors that influence cell numbers. Abnormal expression of oncogenic splicing isoforms and splicing factors, as well as disruptions in splicing caused by genetic mutations, have been implicated in the development and progression of tumors. Collectively, these findings provide valuable insights into the complex interplay between alternative splicing and cell proliferation, thereby suggesting the potential of alternative splicing as a therapeutic target for cancer.

Comparative analyses of dynamic transcriptome profiles highlight key response genes and dominant isoforms for muscle development and growth in chicken

BACKGROUND

Modern breeding strategies have resulted in significant differences in muscle mass between indigenous chicken and specialized broiler. However, the molecular regulatory mechanisms that underlie these differences remain elusive. The aim of this study was to identify key genes and regulatory mechanisms underlying differences in breast muscle development between indigenous chicken and specialized broiler.

RESULTS

Two time-series RNA-sequencing profiles of breast muscles were generated from commercial Arbor Acres (AA) broiler (fast-growing) and Chinese indigenous Lushi blue-shelled-egg (LS) chicken (slow-growing) at embryonic days 10, 14, and 18, and post-hatching day 1 and weeks 1, 3, and 5. Principal component analysis of the transcriptome profiles showed that the top four principal components accounted for more than 80% of the total variance in each breed. The developmental axes between the AA and LS chicken overlapped at the embryonic stages but gradually separated at the adult stages. Integrative investigation of differentially-expressed transcripts contained in the top four principal components identified 44 genes that formed a molecular network associated with differences in breast muscle mass between the two breeds. In addition, alternative splicing analysis revealed that genes with multiple isoforms always had one dominant transcript that exhibited a significantly higher expression level than the others. Among the 44 genes, the TNFRSF6B gene, a mediator of signal transduction pathways and cell proliferation, harbored two alternative splicing isoforms, TNFRSF6B-X1 and TNFRSF6B-X2. TNFRSF6B-X1 was the dominant isoform in both breeds before the age of one week. A switching event of the dominant isoform occurred at one week of age, resulting in TNFRSF6B-X2 being the dominant isoform in AA broiler, whereas TNFRSF6B-X1 remained the dominant isoform in LS chicken. Gain-of-function assays demonstrated that both isoforms promoted the proliferation of chicken primary myoblasts, but only TNFRSF6B-X2 augmented the differentiation and intracellular protein content of chicken primary myoblasts.

CONCLUSIONS

For the first time, we identified several key genes and dominant isoforms that may be responsible for differences in muscle mass between slow-growing indigenous chicken and fast-growing commercial broiler. These findings provide new insights into the regulatory mechanisms underlying breast muscle development in chicken.

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.

Cooperative and competitive regulation of the astrocytic transcriptome by neurons and endothelial cells: Impact on astrocyte maturation

Astrocytes have essential roles in central nervous system (CNS) health and disease. During development, immature astrocytes show complex interactions with neurons, endothelial cells, and other glial cell types. Our work and that of others have shown that these interactions are important for astrocytic maturation. However, whether and how these cells work together to control this process remains poorly understood. Here, we test the hypothesis that cooperative interactions of astrocytes with neurons and endothelial cells promote astrocytic maturation. Astrocytes were cultured alone, with neurons, endothelial cells, or a combination of both. This was followed by astrocyte sorting, RNA sequencing, and bioinformatic analysis to detect transcriptional changes. Across culture configurations, 7302 genes were differentially expressed by 4 or more fold and organized into 8 groups that demonstrate cooperative and antagonist effects of neurons and endothelia on astrocytes. We also discovered that neurons and endothelial cells caused splicing of 200 and 781 mRNAs, respectively. Changes in gene expression were validated using quantitative PCR, western blot (WB), and immunofluorescence analysis. We found that the transcriptomic data from the three-culture configurations correlated with protein expression of three representative targets (FAM107A, GAT3, and GLT1) in vivo. Alternative splicing results also correlated with cortical tissue isoform representation of a target (Fibronectin 1) at different developmental stages. By comparing our results to published transcriptomes of immature and mature astrocytes, we found that neurons or endothelia shift the astrocytic transcriptome toward a mature state and that the presence of both cell types has a greater effect on maturation than either cell alone. These results increase our understanding of cellular interactions/pathways that contribute to astrocytic maturation. They also provide insight into how alterations to neurons and/or endothelial cells may alter astrocytes with implications for astrocytic changes in CNS disorders and diseases.

Mechanisms of mRNA processing defects in inherited THOC6 intellectual disability syndrome

THOC6 is the genetic basis of autosomal recessive THOC6 Intellectual Disability Syndrome (TIDS). THOC6 facilitates the formation of the Transcription Export complex (TREX) tetramer, composed of four THO monomers. The TREX tetramer supports mammalian mRNA processing that is distinct from yeast TREX dimer functions. Human and mouse TIDS model systems allow novel THOC6-dependent TREX tetramer functions to be investigated. Biallelic loss-of-functon(LOF) THOC6 variants do not influence the expression and localization of TREX members in human cells, but our data suggests reduced binding affinity of ALYREF. Impairment of TREX nuclear export functions were not detected in cells with biallelic THOC6 LOF. Instead, mRNA mis-splicing was observed in human and mouse neural tissue, revealing novel insights into THOC6-mediated TREX coordination of mRNA processing. We demonstrate that THOC6 is required for regulation of key signaling pathways in human corticogenesis that dictate the transition from proliferative to neurogenic divisions that may inform TIDS neuropathology.

Comprehensive isoform-level analysis reveals the contribution of alternative isoforms to venom evolution and repertoire diversity

Animal venom systems have emerged as valuable models for investigating how novel polygenic phenotypes may arise from gene evolution by varying molecular mechanisms. However, a significant portion of venom genes produce alternative mRNA isoforms that have not been extensively characterized, hindering a comprehensive understanding of venom biology. In this study, we present a full-length isoform-level profiling workflow integrating multiple RNA sequencing technologies, allowing us to reconstruct a high-resolution transcriptome landscape of venom genes in the parasitoid wasp Pteromalus puparum Our findings demonstrate that more than half of the venom genes generate multiple isoforms within the venom gland. Through mass spectrometry analysis, we confirm that alternative splicing contributes to the diversity of venom proteins, acting as a mechanism for expanding the venom repertoire. Notably, we identified seven venom genes that exhibit distinct isoform usages between the venom gland and other tissues. Furthermore, evolutionary analyses of venom serpin3 and orcokinin further reveal that the co-option of an ancient isoform and a newly evolved isoform, respectively, contributes to venom recruitment, providing valuable insights into the genetic mechanisms driving venom evolution in parasitoid wasps. Overall, our study presents a comprehensive investigation of venom genes at the isoform level, significantly advancing our understanding of alternative isoforms in venom diversity and evolution and setting the stage for further in-depth research on venoms.

A brain specific alternatively spliced isoform of nonmuscle myosin IIA lacks its mechanoenzymatic activities

Recent genomic studies reported that 90 to 95% of human genes can undergo alternative splicing, by which multiple isoforms of proteins are synthesized. However, the functional consequences of most of the isoforms are largely unknown. Here, we report a novel alternatively spliced isoform of nonmuscle myosin IIA (NM IIA), called NM IIA2, which is generated by the inclusion of 21 amino acids near the actin-binding region (loop 2) of the head domain of heavy chains. Expression of NM IIA2 is found exclusively in the brain tissue, where it reaches a maximum level at 24 h during the circadian rhythm. The actin-dependent Mg2+-ATPase activity and in vitro motility assays reveal that NM IIA2 lacks its motor activities but localizes with actin filaments in cells. Interestingly, NM IIA2 can also make heterofilaments with NM IIA0 (noninserted isoform of NM IIA) and can retard the in vitro motility of NM IIA, when the two are mixed. Altogether, our findings provide the functional importance of a previously unknown alternatively spliced isoform, NM IIA2, and its potential physiological role in regulating NM IIA activity in the brain.

Integrative analysis of transcriptome dynamics during human craniofacial development identifies candidate disease genes

Craniofacial disorders arise in early pregnancy and are one of the most common congenital defects. To fully understand how craniofacial disorders arise, it is essential to characterize gene expression during the patterning of the craniofacial region. To address this, we performed bulk and single-cell RNA-seq on human craniofacial tissue from 4-8 weeks post conception. Comparisons to dozens of other human tissues revealed 239 genes most strongly expressed during craniofacial development. Craniofacial-biased developmental enhancers were enriched +/- 400 kb surrounding these craniofacial-biased genes. Gene co-expression analysis revealed that regulatory hubs are enriched for known disease causing genes and are resistant to mutation in the normal healthy population. Combining transcriptomic and epigenomic data we identified 539 genes likely to contribute to craniofacial disorders. While most have not been previously implicated in craniofacial disorders, we demonstrate this set of genes has increased levels of de novo mutations in orofacial clefting patients warranting further study.

Cardiac splicing as a diagnostic and therapeutic target

Despite advances in therapeutics for heart failure and arrhythmias, a substantial proportion of patients with cardiomyopathy do not respond to interventions, indicating a need to identify novel modifiable myocardial pathobiology. Human genetic variation associated with severe forms of cardiomyopathy and arrhythmias has highlighted the crucial role of alternative splicing in myocardial health and disease, given that it determines which mature RNA transcripts drive the mechanical, structural, signalling and metabolic properties of the heart. In this Review, we discuss how the analysis of cardiac isoform expression has been facilitated by technical advances in multiomics and long-read and single-cell sequencing technologies. The resulting insights into the regulation of alternative splicing - including the identification of cardiac splice regulators as therapeutic targets and the development of a translational pipeline to evaluate splice modulators in human engineered heart tissue, animal models and clinical trials - provide a basis for improved diagnosis and therapy. Finally, we consider how the medical and scientific communities can benefit from facilitated acquisition and interpretation of splicing data towards improved clinical decision-making and patient care.

Cytoplasmic aggregation of mutant FUS causes multistep RNA splicing perturbations in the course of motor neuron pathology

Dysfunction of the RNA-binding protein (RBP) FUS implicated in RNA metabolism can cause amyotrophic lateral sclerosis (ALS) and other neurodegenerative diseases. Mutations affecting FUS nuclear localization can drive RNA splicing defects and stimulate the formation of non-amyloid inclusions in affected neurons. However, the mechanism by which FUS mutations contribute to the development of ALS remains uncertain. Here we describe a pattern of RNA splicing changes in the dynamics of the continuous proteinopathy induced by mislocalized FUS. We show that the decrease in intron retention of FUS-associated transcripts represents the hallmark of the pathogenesis of ALS and is the earliest molecular event in the course of progression of the disease. As FUS aggregation increases, the pattern of RNA splicing changes, becoming more complex, including a decrease in the inclusion of neuron-specific microexons and induction of cryptic exon splicing due to the sequestration of additional RBPs into FUS aggregates. Crucially, the identified features of the pathological splicing pattern are also observed in ALS patients in both sporadic and familial cases. Our data provide evidence that both a loss of nuclear FUS function due to mislocalization and the subsequent cytoplasmic aggregation of mutant protein lead to the disruption of RNA splicing in a multistep fashion during FUS aggregation.

Evidence-Based Guide to Using Artificial Introns for Tissue-Specific Knockout in Mice

Up until recently, methods for generating floxed mice either conventionally or by CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)-Cas9 (CRISPR-associated protein 9) editing have been technically challenging, expensive and error-prone, or time-consuming. To circumvent these issues, several labs have started successfully using a small artificial intron to conditionally knockout (KO) a gene of interest in mice. However, many other labs are having difficulty getting the technique to work. The key problem appears to be either a failure in achieving correct splicing after the introduction of the artificial intron into the gene or, just as crucial, insufficient functional KO of the gene's protein after Cre-induced removal of the intron's branchpoint. Presented here is a guide on how to choose an appropriate exon and where to place the recombinase-regulated artificial intron (rAI) in that exon to prevent disrupting normal gene splicing while maximizing mRNA degradation after recombinase treatment. The reasoning behind each step in the guide is also discussed. Following these recommendations should increase the success rate of this easy, new, and alternative technique for producing tissue-specific KO mice.

TreeTerminus -creating transcript trees using inferential replicate counts

A certain degree of uncertainty is always associated with the transcript abundance estimates. The uncertainty may make many downstream analyses, such as differential testing, difficult for certain transcripts. Conversely, gene-level analysis, though less ambiguous, is often too coarse-grained. We introduce TreeTerminus, a data-driven approach for grouping transcripts into a tree structure where leaves represent individual transcripts and internal nodes represent an aggregation of a transcript set. TreeTerminus constructs trees such that, on average, the inferential uncertainty decreases as we ascend the tree topology. The tree provides the flexibility to analyze data at nodes that are at different levels of resolution in the tree and can be tuned depending on the analysis of interest. We evaluated TreeTerminus on two simulated and two experimental datasets and observed an improved performance compared to transcripts (leaves) and other methods under several different metrics.

The role of Limch1 alternative splicing in skeletal muscle function

Postnatal skeletal muscle development is a highly dynamic period associated with widespread alternative splicing changes required to adapt tissues to adult function. These splicing events have significant implications because the reversion of adult mRNA isoforms to fetal isoforms is observed in forms of muscular dystrophy. LIMCH1 is a stress fiber-associated protein that is alternatively spliced to generate uLIMCH1, a ubiquitously expressed isoform, and mLIMCH1, a skeletal muscle-specific isoform containing six additional exons simultaneously included after birth in the mouse. CRISPR/Cas9 was used to delete the six alternatively spliced exons of LIMCH1 in mice, thereby forcing the constitutive expression of the predominantly fetal isoform, uLIMCH1. mLIMCH1 knockout mice had significant grip strength weakness in vivo, and maximum force generated was decreased ex vivo. Calcium-handling deficits were observed during myofiber stimulation that could explain the mechanism by which mLIMCH1 knockout leads to muscle weakness. In addition, LIMCH1 is mis-spliced in myotonic dystrophy type 1, with the muscleblind-like (MBNL) family of proteins acting as the likely major regulator of Limch1 alternative splicing in skeletal muscle.

Splicing complexity as a pivotal feature of alternative exons in mammalian species

BACKGROUND

As a significant process of post-transcriptional gene expression regulation in eukaryotic cells, alternative splicing (AS) of exons greatly contributes to the complexity of the transcriptome and indirectly enriches the protein repertoires. A large number of studies have focused on the splicing inclusion of alternative exons and have revealed the roles of AS in organ development and maturation. Notably, AS takes place through a change in the relative abundance of the transcript isoforms produced by a single gene, meaning that exons can have complex splicing patterns. However, the commonly used percent spliced-in (Ψ) values only define the usage rate of exons, but lose information about the complexity of exons' linkage pattern. To date, the extent and functional consequence of splicing complexity of alternative exons in development and evolution is poorly understood.

RESULTS

By comparing splicing complexity of exons in six tissues (brain, cerebellum, heart, liver, kidney, and testis) from six mammalian species (human, chimpanzee, gorilla, macaque, mouse, opossum) and an outgroup species (chicken), we revealed that exons with high splicing complexity are prevalent in mammals and are closely related to features of genes. Using traditional machine learning and deep learning methods, we found that the splicing complexity of exons can be moderately predicted with features derived from exons, among which length of flanking exons and splicing strength of downstream/upstream splice sites are top predictors. Comparative analysis among human, chimpanzee, gorilla, macaque, and mouse revealed that, alternative exons tend to evolve to an increased level of splicing complexity and higher tissue specificity in splicing complexity. During organ development, not only developmentally regulated exons, but also 10-15% of non-developmentally regulated exons show dynamic splicing complexity.

CONCLUSIONS

Our analysis revealed that splicing complexity is an important metric to characterize the splicing dynamics of alternative exons during the development and evolution of mammals.

Splicing factors in the heart: Uncovering shared and unique targets

Alternative splicing generates specialized protein isoforms that allow the heart to adapt during development and disease. The recent discovery that mutations in the splicing factor RNA-binding protein 20 (RBM20) cause a severe form of familial dilated cardiomyopathy has sparked a great interest in alternative splicing in the field of cardiology. Since then, identification of splicing factors controlling alternative splicing in the heart has grown at a rapid pace. Despite the intriguing observation that a certain overlap exists between the targets of some splicing factors, an integrated and systematic analysis of their splicing networks is missing. Here, we compared the splicing networks of individual splicing factors by re-analyzing original RNA-sequencing data from eight previously published mouse models, in which a single splicing factor has been genetically deleted (i.e. HNRNPU, MBNL1/2, QKI, RBM20, RBM24, RBPMS, SRSF3, SRSF4). We show that key splicing events in Camk2d, Ryr2, Tpm1, Tpm2 and Pdlim5 require the combined action of the majority of these splicing factors. Additionally, we identified common targets and pathways among splicing factors, with the largest overlap between the splicing networks of MBNL, QKI and RBM24. We also re-analyzed a large-scale RNA-sequencing study on hearts of 128 heart failure patients. Here, we observed that MBNL1, QKI and RBM24 expression varied greatly. This variation in expression correlated with differential splicing of their downstream targets as found in mice, suggesting that aberrant splicing by MBNL1, QKI and RBM24 might contribute to the disease mechanism in heart failure.

Single-cell long-read mRNA isoform regulation is pervasive across mammalian brain regions, cell types, and development

RNA isoforms influence cell identity and function. Until recently, technological limitations prevented a genome-wide appraisal of isoform influence on cell identity in various parts of the brain. Using enhanced long-read single-cell isoform sequencing, we comprehensively analyze RNA isoforms in multiple mouse brain regions, cell subtypes, and developmental timepoints from postnatal day 14 (P14) to adult (P56). For 75% of genes, full-length isoform expression varies along one or more axes of phenotypic origin, underscoring the pervasiveness of isoform regulation across multiple scales. As expected, splicing varies strongly between cell types. However, certain gene classes including neurotransmitter release and reuptake as well as synapse turnover, harbor significant variability in the same cell type across anatomical regions, suggesting differences in network activity may influence cell-type identity. Glial brain-region specificity in isoform expression includes strong poly(A)-site regulation, whereas neurons have stronger TSS regulation. Furthermore, developmental patterns of cell-type specific splicing are especially pronounced in the murine adolescent transition from P21 to P28. The same cell type traced across development shows more isoform variability than across adult anatomical regions, indicating a coordinated modulation of functional programs dictating neural development. As most cell-type specific exons in P56 mouse hippocampus behave similarly in newly generated data from human hippocampi, these principles may be extrapolated to human brain. However, human brains have evolved additional cell-type specificity in splicing, suggesting gain-of-function isoforms. Taken together, we present a detailed single-cell atlas of full-length brain isoform regulation across development and anatomical regions, providing a previously unappreciated degree of isoform variability across multiple scales of the brain.

MARVEL: an integrated alternative splicing analysis platform for single-cell RNA sequencing data

Alternative splicing is an important source of heterogeneity underlying gene expression between individual cells but remains an understudied area due to the paucity of computational tools to analyze splicing dynamics at single-cell resolution. Here, we present MARVEL, a comprehensive R package for single-cell splicing analysis applicable to RNA sequencing generated from the plate- and droplet-based methods. We performed extensive benchmarking of MARVEL against available tools and demonstrated its utility by analyzing multiple publicly available datasets in diverse cell types, including in disease. MARVEL enables systematic and integrated splicing and gene expression analysis of single cells to characterize the splicing landscape and reveal biological insights.

Transcript Isoform Diversity of Ampliconic Genes on the Y Chromosome of Great Apes

Y-chromosomal Ampliconic Genes (YAGs) are important for male fertility, as they encode proteins functioning in spermatogenesis. The variation in copy number and expression levels of these multicopy gene families has been recently studied in great apes, however, the diversity of splicing variants remains unexplored. Here we deciphered the sequences of polyadenylated transcripts of all nine YAG families (BPY2, CDY, DAZ, HSFY, PRY, RBMY, TSPY, VCY, and XKRY) from testis samples of six great ape species (human, chimpanzee, bonobo, gorilla, Bornean orangutan, and Sumatran orangutan). To achieve this, we enriched YAG transcripts with capture-probe hybridization and sequenced them with long (Pacific Biosciences) reads. Our analysis of this dataset resulted in several findings. First, we uncovered a high diversity of YAG transcripts across great apes. Second, we observed evolutionarily conserved alternative splicing patterns for most YAG families except for BPY2 and PRY. Our results suggest that BPY2 transcripts and predicted proteins in several great ape species (bonobo and the two orangutans) have independent evolutionary origins and are not homologous to human reference transcripts and proteins. In contrast, our results suggest that the PRY gene family, having the highest representation of transcripts without open reading frames, has been undergoing pseudogenization. Third, even though we have identified many species-specific protein-coding YAG transcripts, we have not detected any signatures of positive selection. Overall, our work illuminates the YAG isoform landscape and its evolutionary history, and provides a genomic resource for future functional studies focusing on infertility phenotypes in humans and critically endangered great apes.

A novel splicing variant of DJ-1 in Parkinson's disease induces mitochondrial dysfunction

Several studies have identified mutations in neuroprotective genes in a few cases of Parkinson's disease (PD); however, the role of alternative splicing changes in PD remains unelucidated. Based on the transcriptome analysis of substantia nigra (SN) tissues obtained from PD cases and age-matched healthy controls, we identified a novel alternative splicing variant of DJ-1, lacking exon 6 (DJ-1 ^ΔE6), frequently detected in the SN of patients with PD. We found that the exon 6 skipping of DJ-1 induces mitochondrial dysfunction and impaired antioxidant capability. According to an in silico modeling study, the exon 6 skipping of DJ-1 disrupts the structural state suitable for the oxidation of the cysteine 106 residue that is a prerequisite for activating its neuroprotective roles. Our results suggest that change in DJ-1 alternative splicing may contribute to PD progression and provide an insight for studying PD etiology and its potential therapeutic targets.

A cell-type-specific alternative splicing regulator shapes synapse properties in a trans-synaptic manner

The specification of synaptic properties is fundamental for the function of neuronal circuits. "Terminal selector" transcription factors coordinate terminal gene batteries that specify cell-type-specific properties. Moreover, pan-neuronal splicing regulators have been implicated in directing neuronal differentiation. However, the cellular logic of how splicing regulators instruct specific synaptic properties remains poorly understood. Here, we combine genome-wide mapping of mRNA targets and cell-type-specific loss-of-function studies to uncover the contribution of the RNA-binding protein SLM2 to hippocampal synapse specification. Focusing on pyramidal cells and somatostatin (SST)-positive GABAergic interneurons, we find that SLM2 preferentially binds and regulates alternative splicing of transcripts encoding synaptic proteins. In the absence of SLM2, neuronal populations exhibit normal intrinsic properties, but there are non-cell-autonomous synaptic phenotypes and associated defects in a hippocampus-dependent memory task. Thus, alternative splicing provides a critical layer of gene regulation that instructs specification of neuronal connectivity in a trans-synaptic manner.

Elucidation of the Landscape of Alternatively Spliced Genes and Features in the Dorsal Striatum of Aggressive/Aggression-Deprived Mice in the Model of Chronic Social Conflicts

Both aggressive and aggression-deprived (AD) individuals represent pathological cases extensively studied in psychiatry and substance abuse disciplines. We employed the animal model of chronic social conflicts curated in our laboratory for over 30 years. In the study, we pursued the task of evaluation of the key events in the dorsal striatum transcriptomes of aggression-experienced mice and AD species, as compared with the controls, using RNA-seq profiling. We evaluated the alternative splicing-mediated transcriptome dynamics based on the RNA-seq data. We confined our attention to the exon skipping (ES) events as the major AS type for animals. We report the concurrent posttranscriptional and posttranslational regulation of the ES events observed in the phosphorylation cycles (in phosphoproteins and their targets) in the neuron-specific genes of the striatum. Strikingly, we found that major neurospecific splicing factors (Nova1, Ptbp1, 2, Mbnl1, 2, and Sam68) related to the alternative splicing regulation of cAMP genes (Darpp-32, Grin1, Ptpn5, Ppp3ca, Pde10a, Prkaca, Psd95, and Adora1) are upregulated specifically in aggressive individuals as compared with the controls and specifically AD animals, assuming intense switching between isoforms in the cAMP-mediated (de)phosphorylation signaling cascade. We found that the coding alternative splicing events were mostly attributed to synaptic plasticity and neural development-related proteins, while the nonsense-mediated decay-associated splicing events are mostly attributed to the mRNA processing of genes, including the spliceosome and splicing factors. In addition, considering the gene families, the transporter (Slc) gene family manifested most of the ES events. We found out that the major molecular systems employing AS for their plasticity are the ‘spliceosome’, ‘chromatin rearrangement complex’, ‘synapse’, and ‘neural development/axonogenesis’ GO categories. Finally, we state that approximately 35% of the exon skipping variants in gene coding regions manifest the noncoding variants subject to nonsense-mediated decay, employed as a homeostasis-mediated expression regulation layer and often associated with the corresponding gene expression alteration.

Post-Transcriptional Modification by Alternative Splicing and Pathogenic Splicing Variants in Cardiovascular Development and Congenital Heart Defects

Advancements in genomics, bioinformatics, and genome editing have uncovered new dimensions in gene regulation. Post-transcriptional modifications by the alternative splicing of mRNA transcripts are critical regulatory mechanisms of mammalian gene expression. In the heart, there is an expanding interest in elucidating the role of alternative splicing in transcriptome regulation. Substantial efforts were directed toward investigating this process in heart development and failure. However, few studies shed light on alternative splicing products and their dysregulation in congenital heart defects (CHDs). While elegant reports showed the crucial roles of RNA binding proteins (RBPs) in orchestrating splicing transitions during heart development and failure, the impact of RBPs dysregulation or genetic variation on CHDs has not been fully addressed. Herein, we review the current understanding of alternative splicing and RBPs' roles in heart development and CHDs. Wediscuss the impact of perinatal splicing transition and its dysregulation in CHDs. We further summarize the discoveries made of causal splicing variants in key transcription factors that are implicated in CHDs. An improved understanding of the roles of alternative splicing in heart development and CHDs may potentially inform novel preventive and therapeutic advancements for newborn infants with CHDs.

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.

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.

The molecular evolution of spermatogenesis across mammals

The testis produces gametes through spermatogenesis and evolves rapidly at both the morphological and molecular level in mammals^1-6, probably owing to the evolutionary pressure on males to be reproductively successful⁷. However, the molecular evolution of individual spermatogenic cell types across mammals remains largely uncharacterized. Here we report evolutionary analyses of single-nucleus transcriptome data for testes from 11 species that cover the three main mammalian lineages (eutherians, marsupials and monotremes) and birds (the evolutionary outgroup), and include seven primates. We find that the rapid evolution of the testis was driven by accelerated fixation rates of gene expression changes, amino acid substitutions and new genes in late spermatogenic stages, probably facilitated by reduced pleiotropic constraints, haploid selection and transcriptionally permissive chromatin. We identify temporal expression changes of individual genes across species and conserved expression programs controlling ancestral spermatogenic processes. Genes predominantly expressed in spermatogonia (germ cells fuelling spermatogenesis) and Sertoli (somatic support) cells accumulated on X chromosomes during evolution, presumably owing to male-beneficial selective forces. Further work identified transcriptomal differences between X- and Y-bearing spermatids and uncovered that meiotic sex-chromosome inactivation (MSCI) also occurs in monotremes and hence is common to mammalian sex-chromosome systems. Thus, the mechanism of meiotic silencing of unsynapsed chromatin, which underlies MSCI, is an ancestral mammalian feature. Our study illuminates the molecular evolution of spermatogenesis and associated selective forces, and provides a resource for investigating the biology of the testis across mammals.

AKT1 phosphorylates RBM17 to promote Sox2 transcription by modulating alternative splicing of FOXM1 to enhance cancer stem cell properties in colorectal cancer cells

Colorectal cancer (CRC) is one of the leading causes of cancer-related death worldwide. The existence of cancer stem cells (CSC) causes tumor relapses, metastasis, and resistance to conventional therapy. Alternative splicing has been shown to affect physiological and pathological processes. Accumulating evidence has confirmed that targeting alternative splicing could be an effective strategy to treat CRC. Currently, the role of alternative splicing in the regulation of CSC properties in CRC has not been elucidated. Here, we show that RBM17 displays oncogenic roles in CRC cells. RBM17 enhances cell proliferation and reduces chemotherapeutic-induced apoptosis in CRC cells. Besides, RBM17 increases CD133 positive and ALDEFLUOR positive populations and promotes sphere formation in CRC cells. In mechanism studies, we found that FOXM1 is critical for RBM17 enhanced CSC properties. Moreover, FOXM1 alternative splicing is essential for RBM17 enhanced CSC properties in CRC cells. Additionally, RBM17 enhances CSC characteristics by controlling FOXM1 expression to promote Sox2 expression. Furthermore, AKT1 works as an upstream kinase to control RBM17-mediated FOXM1 alternative splicing and enhancement of CSC properties in CRC cells. Our study reveals that AKT1-RBM17-FOXM1-Sox2 axis could be a potential target for modulating alternative splicing to reduce CSC properties in CRC cells.

The landscape of expression and alternative splicing variation across human traits

Understanding the consequences of individual transcriptome variation is fundamental to deciphering human biology and disease. We implement a statistical framework to quantify the contributions of 21 individual traits as drivers of gene expression and alternative splicing variation across 46 human tissues and 781 individuals from the Genotype-Tissue Expression project. We demonstrate that ancestry, sex, age, and BMI make additive and tissue-specific contributions to expression variability, whereas interactions are rare. Variation in splicing is dominated by ancestry and is under genetic control in most tissues, with ribosomal proteins showing a strong enrichment of tissue-shared splicing events. Our analyses reveal a systemic contribution of types 1 and 2 diabetes to tissue transcriptome variation with the strongest signal in the nerve, where histopathology image analysis identifies novel genes related to diabetic neuropathy. Our multi-tissue and multi-trait approach provides an extensive characterization of the main drivers of human transcriptome variation in health and disease.

Broad misappropriation of developmental splicing profile by cancer in multiple organs

Oncogenesis mimics key aspects of embryonic development. However, the underlying mechanisms are incompletely understood. Here, we demonstrate that the splicing events specifically active during human organogenesis, are broadly reactivated in the organ-specific tumor. Such events are associated with key oncogenic processes and predict proliferation rates in cancer cell lines as well as patient survival. Such events preferentially target nitrosylation and transmembrane-region domains, whose coordinated splicing in multiple genes respectively affect intracellular transport and N-linked glycosylation. We infer critical splicing factors potentially regulating embryonic splicing events and show that such factors are potential oncogenic drivers and are upregulated specifically in malignant cells. Multiple complementary analyses point to MYC and FOXM1 as potential transcriptional regulators of critical splicing factors in brain and liver. Our study provides a comprehensive demonstration of a splicing-mediated link between development and cancer, and suggest anti-cancer targets including splicing events, and their upstream splicing and transcriptional regulators.