Alternative splicing as a regulator of development and tissue identity

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.

Splicing factor YBX1 regulates bone marrow stromal cell fate during aging

Senescence and altered differentiation potential of bone marrow stromal cells (BMSCs) lead to age-related bone loss. As an important posttranscriptional regulatory pathway, alternative splicing (AS) regulates the diversity of gene expression and has been linked to induction of cellular senescence. However, the role of splicing factors in BMSCs during aging remains poorly defined. Herein, we found that the expression of the splicing factor Y-box binding protein 1 (YBX1) in BMSCs decreased with aging in mice and humans. YBX1 deficiency resulted in mis-splicing in genes linked to BMSC osteogenic differentiation and senescence, such as Fn1, Nrp2, Sirt2, Sp7, and Spp1, thus contributing to BMSC senescence and differentiation shift during aging. Deletion of Ybx1 in BMSCs accelerated bone loss in mice, while its overexpression stimulated bone formation. Finally, we identified a small compound, sciadopitysin, which attenuated the degradation of YBX1 and bone loss in old mice. Our study demonstrated that YBX1 governs cell fate of BMSCs via fine control of RNA splicing and provides a potential therapeutic target for age-related osteoporosis.

The nexus between RNA-binding proteins and their effectors

RNA-binding proteins (RBPs) regulate essentially every event in the lifetime of an RNA molecule, from its production to its destruction. Whereas much has been learned about RNA sequence specificity and general functions of individual RBPs, the ways in which numerous RBPs instruct a much smaller number of effector molecules, that is, the core engines of RNA processing, as to where, when and how to act remain largely speculative. Here, we survey the known modes of communication between RBPs and their effectors with a particular focus on converging RBP-effector interactions and their roles in reducing the complexity of RNA networks. We discern the emerging unifying principles and discuss their utility in our understanding of RBP function, regulation of biological processes and contribution to human disease.

Presence of immunogenic alternatively spliced insulin gene product in human pancreatic delta cells

AIMS/HYPOTHESIS

Transcriptome analyses revealed insulin-gene-derived transcripts in non-beta endocrine islet cells. We studied alternative splicing of human INS mRNA in pancreatic islets.

METHODS

Alternative splicing of insulin pre-mRNA was determined by PCR analysis performed on human islet RNA and single-cell RNA-seq analysis. Antisera were generated to detect insulin variants in human pancreatic tissue using immunohistochemistry, electron microscopy and single-cell western blot to confirm the expression of insulin variants. Cytotoxic T lymphocyte (CTL) activation was determined by MIP-1β release.

RESULTS

We identified an alternatively spliced INS product. This variant encodes the complete insulin signal peptide and B chain and an alternative C-terminus that largely overlaps with a previously identified defective ribosomal product of INS. Immunohistochemical analysis revealed that the translation product of this INS-derived splice transcript was detectable in somatostatin-producing delta cells but not in beta cells; this was confirmed by light and electron microscopy. Expression of this alternatively spliced INS product activated preproinsulin-specific CTLs in vitro. The exclusive presence of this alternatively spliced INS product in delta cells may be explained by its clearance from beta cells by insulin-degrading enzyme capturing its insulin B chain fragment and a lack of insulin-degrading enzyme expression in delta cells.

CONCLUSIONS/INTERPRETATION

Our data demonstrate that delta cells can express an INS product derived from alternative splicing, containing both the diabetogenic insulin signal peptide and B chain, in their secretory granules. We propose that this alternative INS product may play a role in islet autoimmunity and pathology, as well as endocrine or paracrine function or islet development and endocrine destiny, and transdifferentiation between endocrine cells. INS promoter activity is not confined to beta cells and should be used with care when assigning beta cell identity and selectivity.

DATA AVAILABILITY

The full EM dataset is available via www.nanotomy.org (for review: http://www.nanotomy.org/OA/Tienhoven2021SUB/6126-368/ ). Single-cell RNA-seq data was made available by Segerstolpe et al [13] and can be found at https://sandberglab.se/pancreas . The RNA and protein sequence of INS-splice was uploaded to GenBank (BankIt2546444 INS-splice OM489474).

Regulation of adult stem cell quiescence and its functions in the maintenance of tissue integrity

Adult stem cells are important for mammalian tissues, where they act as a cell reserve that supports normal tissue turnover and can mount a regenerative response following acute injuries. Quiescent stem cells are well established in certain tissues, such as skeletal muscle, brain, and bone marrow. The quiescent state is actively controlled and is essential for long-term maintenance of stem cell pools. In this Review, we discuss the importance of maintaining a functional pool of quiescent adult stem cells, including haematopoietic stem cells, skeletal muscle stem cells, neural stem cells, hair follicle stem cells, and mesenchymal stem cells such as fibro-adipogenic progenitors, to ensure tissue maintenance and repair. We discuss the molecular mechanisms that regulate the entry into, maintenance of, and exit from the quiescent state in mice. Recent studies revealed that quiescent stem cells have a discordance between RNA and protein levels, indicating the importance of post-transcriptional mechanisms, such as alternative polyadenylation, alternative splicing, and translation repression, in the control of stem cell quiescence. Understanding how these mechanisms guide stem cell function during homeostasis and regeneration has important implications for regenerative medicine.

Alternative Splicing Changes Promoted by NOVA2 Upregulation in Endothelial Cells and Relevance for Gastric Cancer

Angiogenesis is crucial for cancer progression. While several anti-angiogenic drugs are in use for cancer treatment, their clinical benefits are unsatisfactory. Thus, a deeper understanding of the mechanisms sustaining cancer vessel growth is fundamental to identify novel biomarkers and therapeutic targets. Alternative splicing (AS) is an essential modifier of human proteome diversity. Nevertheless, AS contribution to tumor vasculature development is poorly known. The Neuro-Oncological Ventral Antigen 2 (NOVA2) is a critical AS regulator of angiogenesis and vascular development. NOVA2 is upregulated in tumor endothelial cells (ECs) of different cancers, thus representing a potential driver of tumor blood vessel aberrancies. Here, we identified novel AS transcripts generated upon NOVA2 upregulation in ECs, suggesting a pervasive role of NOVA2 in vascular biology. In addition, we report that NOVA2 is also upregulated in ECs of gastric cancer (GC), and its expression correlates with poor overall survival of GC patients. Finally, we found that the AS of the Rap Guanine Nucleotide Exchange Factor 6 (RapGEF6), a newly identified NOVA2 target, is altered in GC patients and associated with NOVA2 expression, tumor angiogenesis, and poor patient outcome. Our findings provide a better understanding of GC biology and suggest that AS might be exploited to identify novel biomarkers and therapeutics for anti-angiogenic GC treatments.

Alternative splicing affects synapses in the hippocampus of offspring after maternal fructose exposure during gestation and lactation

Increased fructose over-intake is a global issue. Maternal fructose exposure during gestation and lactation can impair brain development in offspring. However, the effect on synapses is still unknown. For the diversification of RNA and biological functions, alternative splicing (AS) and alternative polyadenylation (APA) are essential. We constructed a maternal high-fructose diet model by administering 13% and 40% fructose water. The student's t-test analyzed the results of RT-qPCR. All other results were analyzed by one-way analysis of variance. The animal behavior experiment results revealed that conditioning and associative memory had been damaged. The proteins that form synapses were consistently low-expressed. In addition, compared with the control group, the Oxford Nanopore Technologies platform's full-length RNA-sequencing identified 298 different spliced genes (DSGs) and 51 differentially expressed alternative splicing (DEAS) genes in the 13% fructose group. 313 DSGs and 74 DEAS genes were in the 40% fructose group. Enrichment analysis based on these altered genes revealed some enlightening items and pathways. Our findings demonstrated the transcriptome mechanism underlying maternal fructose exposure during gestation and lactation and impaired synapse function during the transcripts' editing.

The interplay between non-coding RNAs and alternative splicing: from regulatory mechanism to therapeutic implications in cancer

Alternative splicing (AS) is a common and conserved process in eukaryotic gene regulation. It occurs in approximately 95% of multi-exon genes, greatly enriching the complexity and diversity of mRNAs and proteins. Recent studies have found that in addition to coding RNAs, non-coding RNAs (ncRNAs) are also inextricably linked with AS. Multiple different types of ncRNAs are generated by AS of precursor long non-coding (pre-lncRNAs) or precursor messenger RNAs (pre-mRNAs). Furthermore, ncRNAs, as a novel class of regulators, can participate in AS regulation by interacting with the cis-acting elements or trans-acting factors. Several studies have implicated abnormal expression of ncRNAs and ncRNA-related AS events in the initiation, progression, and therapy resistance in various types of cancers. Therefore, owing to their roles in mediating drug resistance, ncRNAs, AS-related factors and AS-related novel antigens may serve as promising therapeutic targets in cancer treatment. In this review, we summarize the interaction between ncRNAs and AS processes, emphasizing their great influences on cancer, especially on chemoresistance, and highlighting their potential values in clinical treatment.

Two Arabidopsis Splicing Factors, U2AF65a and U2AF65b, Differentially Control Flowering Time by Modulating the Expression or Alternative Splicing of a Subset of FLC Upstream Regulators

We investigated the transcriptomic changes in the shoot apices during floral transition in Arabidopsis mutants of two closely related splicing factors: AtU2AF65a (atu2af65a) and AtU2AF65b (atu2af65b). The atu2af65a mutants exhibited delayed flowering, while the atu2af65b mutants showed accelerated flowering. The underlying gene regulatory mechanism of these phenotypes was unclear. We performed RNA-seq analysis using shoot apices instead of whole seedlings and found that the atu2af65a mutants had more differentially expressed genes than the atu2af65b mutants when they were compared to wild type. The only flowering time gene that was significantly up- or down-regulated by more than two-fold in the mutants were FLOWERING LOCUS C (FLC), a major floral repressor. We also examined the expression and alternative splicing (AS) patterns of several FLC upstream regulators, such as COOLAIR, EDM2, FRIGIDA, and PP2A-b'ɤ, and found that those of COOLAIR, EDM2, and PP2A-b'ɤ were altered in the mutants. Furthermore, we demonstrated that AtU2AF65a and AtU2AF65b genes partially influenced FLC expression by analyzing these mutants in the flc-3 mutant background. Our findings indicate that AtU2AF65a and AtU2AF65b splicing factors modulate FLC expression by affecting the expression or AS patterns of a subset of FLC upstream regulators in the shoot apex, leading to different flowering phenotypes.

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.

Bi-allelic SNAPC4 variants dysregulate global alternative splicing and lead to neuroregression and progressive spastic paraparesis

The vast majority of human genes encode multiple isoforms through alternative splicing, and the temporal and spatial regulation of those isoforms is critical for organismal development and function. The spliceosome, which regulates and executes splicing reactions, is primarily composed of small nuclear ribonucleoproteins (snRNPs) that consist of small nuclear RNAs (snRNAs) and protein subunits. snRNA gene transcription is initiated by the snRNA-activating protein complex (SNAPc). Here, we report ten individuals, from eight families, with bi-allelic, deleterious SNAPC4 variants. SNAPC4 encoded one of the five SNAPc subunits that is critical for DNA binding. Most affected individuals presented with delayed motor development and developmental regression after the first year of life, followed by progressive spasticity that led to gait alterations, paraparesis, and oromotor dysfunction. Most individuals had cerebral, cerebellar, or basal ganglia volume loss by brain MRI. In the available cells from affected individuals, SNAPC4 abundance was decreased compared to unaffected controls, suggesting that the bi-allelic variants affect SNAPC4 accumulation. The depletion of SNAPC4 levels in HeLa cell lines via genomic editing led to decreased snRNA expression and global dysregulation of alternative splicing. Analysis of available fibroblasts from affected individuals showed decreased snRNA expression and global dysregulation of alternative splicing compared to unaffected cells. Altogether, these data suggest that these bi-allelic SNAPC4 variants result in loss of function and underlie the neuroregression and progressive spasticity in these affected individuals.

Etiology of oncogenic fusions in 5,190 childhood cancers and its clinical and therapeutic implication

Oncogenic fusions formed through chromosomal rearrangements are hallmarks of childhood cancer that define cancer subtype, predict outcome, persist through treatment, and can be ideal therapeutic targets. However, mechanistic understanding of the etiology of oncogenic fusions remains elusive. Here we report a comprehensive detection of 272 oncogenic fusion gene pairs by using tumor transcriptome sequencing data from 5190 childhood cancer patients. We identify diverse factors, including translation frame, protein domain, splicing, and gene length, that shape the formation of oncogenic fusions. Our mathematical modeling reveals a strong link between differential selection pressure and clinical outcome in CBFB-MYH11. We discover 4 oncogenic fusions, including RUNX1-RUNX1T1, TCF3-PBX1, CBFA2T3-GLIS2, and KMT2A-AFDN, with promoter-hijacking-like features that may offer alternative strategies for therapeutic targeting. We uncover extensive alternative splicing in oncogenic fusions including KMT2A-MLLT3, KMT2A-MLLT10, C11orf95-RELA, NUP98-NSD1, KMT2A-AFDN and ETV6-RUNX1. We discover neo splice sites in 18 oncogenic fusion gene pairs and demonstrate that such splice sites confer therapeutic vulnerability for etiology-based genome editing. Our study reveals general principles on the etiology of oncogenic fusions in childhood cancer and suggests profound clinical implications including etiology-based risk stratification and genome-editing-based therapeutics.

Preclinical Screening of Splice-Switching Antisense Oligonucleotides in PDAC Organoids

Aberrant alternative splicing is emerging as a cancer hallmark and a potential therapeutic target. It is the result of dysregulated splicing factors or genetic alterations in splicing-regulatory cis -elements. Targeting individual altered splicing events associated with cancer-cell dependencies is a potential therapeutic strategy, but several technical limitations need to be addressed. Patient-derived organoids (PDOs) are a promising platform to recapitulate key aspects of disease states and to facilitate drug development for precision medicine. Here, we report an efficient antisense-oligonucleotide (ASO) transfection method to systematically evaluate and screen individual splicing events as therapeutic targets in pancreatic ductal adenocarcinoma (PDAC) organoids. This optimized delivery method allows fast and efficient screening of ASOs that reverse oncogenic alternative splicing. In combination with advancements in chemical modifications and ASO-delivery strategies, this method has the potential to accelerate the discovery of anti-tumor ASO drugs that target pathological alternative splicing.

Regulation of pre-mRNA splicing: roles in physiology and disease, and therapeutic prospects

The removal of introns from mRNA precursors and its regulation by alternative splicing are key for eukaryotic gene expression and cellular function, as evidenced by the numerous pathologies induced or modified by splicing alterations. Major recent advances have been made in understanding the structures and functions of the splicing machinery, in the description and classification of physiological and pathological isoforms and in the development of the first therapies for genetic diseases based on modulation of splicing. Here, we review this progress and discuss important remaining challenges, including predicting splice sites from genomic sequences, understanding the variety of molecular mechanisms and logic of splicing regulation, and harnessing this knowledge for probing gene function and disease aetiology and for the design of novel therapeutic approaches.

Programmable macromolecule-based RNA-targeting therapies to treat human neurological disorders

Disruptions in RNA processing play critical roles in the pathogenesis of neurological diseases. In this Perspective, we discuss recent progress in the development of RNA-targeting therapeutic modalities. We focus on progress, limitations, and opportunities in a new generation of therapies engineered from RNA binding proteins and other endogenous RNA regulatory macromolecules to treat human neurological disorders.

Alternative splicing: a new breakthrough for understanding tumorigenesis and potential clinical applications

BACKGROUND

Alternative splicing (AS) is a post-transcriptional process that produces transcript variants, thus leading to transcriptome complexity. Recently, the scope of AS studies has been greatly expanded toward clinical applications owing to the abundance of RNA sequencing data.

OBJECTIVE

This review consists of two parts. We first summarize bioinformatic resources that are useful for large-scale cancer-related AS studies. We then highlight the research efforts to utilize AS events for predicting clinical outcomes and planning therapeutic strategies.

RESULTS

Computational approaches to interrogate AS events have been reviewed under three categories: (1) databases to provide functional and clinical annotation of AS events, (2) analytical tools to identify cancer-associated AS event, and (3) methods to identify splicing-related DNA variants and splicing-derived neoantigens. We also present the recent progress in exploring the clinical utility of AS under four categories: (1) identification of AS events for cancer prognosis, (2) utilization of AS events in molecular classification of various cancers, (3) regulatory mechanisms of AS underlying drug resistance, and (4) potential use of AS in cancer therapy.

CONCLUSION

This review will be helpful for understanding the biological implications of AS in cancer and facilitate the development of AS markers for cancer prognosis and treatment. We anticipate that future studies will lead to the application of genome-wide AS profiles in cancer precision medicine.

The physiology of alternative splicing

Alternative splicing is a substantial contributor to the high complexity of transcriptomes of multicellular eukaryotes. In this Review, we discuss the accumulated evidence that most of this complexity is reflected at the protein level and fundamentally shapes the physiology and pathology of organisms. This notion is supported not only by genome-wide analyses but, mainly, by detailed studies showing that global and gene-specific modulations of alternative splicing regulate highly diverse processes such as tissue-specific and species-specific cell differentiation, thermal regulation, neuron self-avoidance, infrared sensing, the Warburg effect, maintenance of telomere length, cancer and autism spectrum disorders (ASD). We also discuss how mastering the control of alternative splicing paved the way to clinically approved therapies for hereditary diseases.

Endothelial deletion of PTBP1 disrupts ventricular chamber development

The growth and maturation of the ventricular chamber require spatiotemporally precise synergy between diverse cell types. Alternative splicing deeply affects the processes. However, the functional properties of alternative splicing in cardiac development are largely unknown. Our study reveals that an alternative splicing factor polypyrimidine tract-binding protein 1 (PTBP1) plays a key role in ventricular chamber morphogenesis. During heart development, PTBP1 colocalizes with endothelial cells but is almost undetectable in cardiomyocytes. The endothelial-specific knockout of Ptbp1, in either endocardial cells or pan-endothelial cells, leads to a typical phenotype of left ventricular noncompaction (LVNC). Mechanistically, the deletion of Ptbp1 reduces the migration of endothelial cells, disrupting cardiomyocyte proliferation and ultimately leading to the LVNC. Further study shows that Ptbp1 deficiency changes the alternative splicing of β-arrestin-1 (Arrb1), which affects endothelial cell migration. In conclusion, as an alternative splicing factor, PTBP1 is essential during ventricular chamber development, and its deficiency can lead to congenital heart disease.

Depending on different apoptosis pathways, the effector Cscaspase-3 in Chilo suppressalis exposed to temperature and parasitic stress was induced

Apoptosis is a form of programmed cell death (PCD) that is largely triggered by caspases through both the mitochondria-dependent and mitochondria-independent pathways. The rice stem borer, Chilo suppressalis, serves as an economically important pest of rice, which is often suffered by temperature and parasitic stress under natural conditions. In the present study, effector Cscaspase-3 encoding caspase was obtained from the rice pest Chilo suppressalis. CsCaspase-3 possesses p20 and p10 subunits, two active sites, four substrate-binding sites, and two cleavage motifs. Real-time quantitative PCR showed that Cscaspase-3 was expressed at maximal levels in hemocytes; furthermore, transcription was most highly in female adults. Expression of Cscaspase-3 was induced by hot and cold temperatures, with the highest expression at 39 °C. Cscaspase-3 expression was also significantly induced at 10 h, 2 d, 5 d, and 7 d of parasitism. Flow cytometry results showed that both temperature and parasitism trigger apoptosis, but only parasitism induces apoptosis via the mitochondrial apoptosis pathway in C. suppressalis. RNAi-mediated silencing of Cscaspase-3 expression reduced C. suppressalis survival at -3 °C. This study provides a foundation for further studies of caspases in insects during biotic and abiotic stress.

Role of ADAM33 short isoform as a tumor suppressor in the pathogenesis of thyroid cancer via oncogenic function disruption of full-length ADAM33

Thyroid cancer is the most prevalent endocrine malignancy globally; however, its underlying pathogenesis remains unclarified. Reportedly, alternative splicing is involved in processes such as embryonic stem and precursor cell differentiation, cell lineage reprogramming, and epithelial-mesenchymal transitions. ADAM33-n, an alternative splicing isoform of ADAM33, encodes a small protein containing 138 amino acids of the N-terminal of full-length ADAM33, which constructs a chaperone-like domain that was previously reported to bind and block the proteolysis activity of ADAM33. In this study, we reported for the first time that ADAM33-n was downregulated in thyroid cancer. The results of cell counting kit-8 and colony formation assays showed that ectopic ADAM33-n in papillary thyroid cancer cell lines restricted cell proliferation and colony formation. Moreover, we demonstrated that ectopic ADAM33-n reversed the oncogenic function of full-length ADAM33 in cell growth and colony formation in the MDA-T32 and BCPAP cells. These findings indicate the tumor suppressor ability of ADAM33-n. Altogether, our study findings present a potential explanatory model of how the downregulation of the oncogenic gene ADAM33 promotes the pathogenesis of thyroid cancer.

Coordination of Alternative Splicing and Alternative Polyadenylation revealed by Targeted Long-Read Sequencing

Nervous system development is associated with extensive regulation of alternative splicing (AS) and alternative polyadenylation (APA). AS and APA have been extensively studied in isolation, but little is known about how these processes are coordinated. Here, the coordination of cassette exon (CE) splicing and APA in Drosophila was investigated using a targeted long-read sequencing approach we call Pull-a-Long-Seq (PL-Seq). This cost-effective method uses cDNA pulldown and Nanopore sequencing combined with an analysis pipeline to resolve the connectivity of alternative exons to alternative 3' ends. Using PL-Seq, we identified genes that exhibit significant differences in CE splicing depending on connectivity to short versus long 3'UTRs. Genomic long 3'UTR deletion was found to alter upstream CE splicing in short 3'UTR isoforms and ELAV loss differentially affected CE splicing depending on connectivity to alternative 3'UTRs. This work highlights the importance of considering connectivity to alternative 3'UTRs when monitoring AS events.

The short isoform of MS4A7 is a novel player in glioblastoma microenvironment, M2 macrophage polarization, and tumor progression

BACKGROUND

The unique intracranial tumor microenvironment (TME) contributes to the immunotherapy failure for glioblastoma (GBM), thus new functional protein targets are urgently needed. Alternative splicing is a widespread regulatory mechanism by which individual gene can express variant proteins with distinct functions. Moreover, proteins located in the cell plasma membrane facilitate targeted therapies. This study sought to obtain functional membrane protein isoforms from GBM TME.

METHODS

With combined single-cell RNA-seq and bulk RNA-seq analyses, novel candidate membrane proteins generated by prognostic splicing events were screened within GBM TME. The short isoform of MS4A7 (MS4A7-s) was selected for evaluation by RT-PCR and western blotting in clinical specimens. Its clinical relevance was evaluated in a GBM patient cohort. The function of MS4A7-s was identified by in vitro and in vivo experiments. MS4A7-s overexpression introduced transcriptome changes were analyzed to explore the potential molecular mechanism.

RESULTS

The main expression product, isoform MS4A7-s, generated by exon skipping, is an M2-specific plasma membrane protein playing a pro-oncogenic role in GBM TME. Higher expression of MS4A7-s correlates with poor prognosis in a GBM cohort. In vitro cell co-culture experiments, intracranial co-injection tumorigenesis assay, and RNA-seq suggest MS4A7-s promotes activation of glioma-associated macrophages' (GAMs) PI3K/AKT/GSK3β pathway, leading to M2 polarization, and drives malignant progression of GBM.

CONCLUSIONS

MS4A7-s, a novel splicing isoform of MS4A7 located on the surface of GAMs in GBM TME, is a predictor of patient outcome, which contributes to M2 polarization and the malignant phenotype of GBM. Targeting MS4A7-s may constitute a promising treatment for GBM.

Single-cell and long-read sequencing to enhance modelling of splicing and cell-fate determination

Single-cell sequencing technologies have revolutionised the life sciences and biomedical research. Single-cell sequencing provides high-resolution data on cell heterogeneity, allowing high-fidelity cell type identification, and lineage tracking. Computational algorithms and mathematical models have been developed to make sense of the data, compensate for errors and simulate the biological processes, which has led to breakthroughs in our understanding of cell differentiation, cell-fate determination and tissue cell composition. The development of long-read (a.k.a. third-generation) sequencing technologies has produced powerful tools for investigating alternative splicing, isoform expression (at the RNA level), genome assembly and the detection of complex structural variants (at the DNA level). In this review, we provide an overview of the recent advancements in single-cell and long-read sequencing technologies, with a particular focus on the computational algorithms that help in correcting, analysing, and interpreting the resulting data. Additionally, we review some mathematical models that use single-cell and long-read sequencing data to study cell-fate determination and alternative splicing, respectively. Moreover, we highlight the emerging opportunities in modelling cell-fate determination that result from the combination of single-cell and long-read sequencing technologies.

Capture RIC-seq reveals positional rules of PTBP1-associated RNA loops in splicing regulation

RNA-binding proteins (RBPs) bind at different positions of the pre-mRNA molecules to promote or reduce the usage of a particular exon. Seeking to understand the working principle of these positional effects, we develop a capture RIC-seq (CRIC-seq) method to enrich specific RBP-associated in situ proximal RNA-RNA fragments for deep sequencing. We determine hnRNPA1-, SRSF1-, and PTBP1-associated proximal RNA-RNA contacts and regulatory mechanisms in HeLa cells. Unexpectedly, the 3D RNA map analysis shows that PTBP1-associated loops in individual introns preferentially promote cassette exon splicing by accelerating asymmetric intron removal, whereas the loops spanning across cassette exon primarily repress splicing. These "positional rules" can faithfully predict PTBP1-regulated splicing outcomes. We further demonstrate that cancer-related splicing quantitative trait loci can disrupt RNA loops by reducing PTBP1 binding on pre-mRNAs to cause aberrant splicing in tumors. Our study presents a powerful method for exploring the functions of RBP-associated RNA-RNA proximal contacts in gene regulation and disease.

Transcriptomic analysis of testis and epididymis tissues from Banna mini-pig inbred line boars with single-molecule long-read sequencing†

In mammals, testis and epididymis are critical components of the male reproductive system for androgen production, spermatogenesis, sperm transportation, as well as sperm maturation. Here, we report single-molecule real-time sequencing data from the testis and epididymis of the Banna mini-pig inbred line (BMI), a promising laboratory animal for medical research. We obtained high-quality full-length transcriptomes and identified 9879 isoforms and 8761 isoforms in the BMI testis and epididymis, respectively. Most of the isoforms we identified have novel exon structures that will greatly improve the annotation of testis- and epididymis-expressed genes in pigs. We also found that 3055 genes (over 50%) were shared between BMI testis and epididymis, indicating widespread expression profiles of genes related to reproduction. We characterized extensive alternative splicing events in BMI testis and epididymis and showed that 96 testis-expressed genes and 79 epididymis-expressed genes have more than six isoforms, revealing the complexity of alternative splicing. We accurately defined the transcribed isoforms in BMI testis and epididymis by combining Pacific Biotechnology Isoform-sequencing (PacBio Iso-Seq) and Illumina RNA Sequencing (RNA-seq) techniques. The refined annotation of some key genes governing male reproduction will facilitate further understanding of the molecular mechanisms underlying BMI male sterility. In addition, the high-confident identification of 548 and 669 long noncoding RNAs (lncRNAs) in these two tissues has established a candidate gene set for future functional investigations. Overall, our study provides new insights into the role of the testis and epididymis during BMI reproduction, paving the path for further studies on BMI male infertility.

A Systematic Method for Detecting Abnormal mRNA Splicing and Assessing Its Clinical Impact in Individuals Undergoing Genetic Testing for Hereditary Cancer Syndromes

Nearly 14% of disease-causing germline variants result from the disruption of mRNA splicing. Most (67%) DNA variants predicted in silico to disrupt splicing are classified as variants of uncertain significance. An analytic workflow-splice effect event resolver (SPEER)-was developed and validated to use mRNA sequencing to reveal significant deviations in splicing, pinpoint the DNA variants potentially involved, and measure the deleterious effects of the altered splicing on mRNA transcripts, providing evidence for assessing the pathogenicity of the variant. SPEER was used to analyze leukocyte RNA encoding 63 hereditary cancer syndrome-related genes in 20,317 patients. Among 3563 patients (17.5%) with at least one DNA variant predicted to affect splicing, 971 (4.8%) had altered splicing with a deleterious effect on the transcript, and 40 had altered splicing due to a DNA variant located outside of the reportable range of the test. Integrating SPEER results into the interpretation of variants allowed variants of uncertain significance to be reclassified as pathogenic or likely pathogenic in 0.4%, and as benign or likely benign in 5.9%, of the 20,317 patients. SPEER-based evidence was associated with a significantly greater effect on classifications of pathogenic or likely pathogenic and benign or likely benign in nonwhite versus non-Hispanic white patients, illustrating that evidence derived from mRNA splicing analysis may help to reduce ethnic/ancestral disparities in genetic testing.

Cohesin regulates alternative splicing

Cohesin, a trimeric complex that establishes sister chromatid cohesion, has additional roles in chromatin organization and transcription. We report that among those roles is the regulation of alternative splicing through direct interactions and in situ colocalization with splicing factors. Degradation of cohesin results in marked changes in splicing, independent of its effects on transcription. Introduction of a single cohesin point mutation in embryonic stem cells alters splicing patterns, demonstrating causality. In primary human acute myeloid leukemia, mutations in cohesin are highly correlated with distinct patterns of alternative splicing. Cohesin also directly interacts with BRD4, another splicing regulator, to generate a pattern of splicing that is distinct from either factor alone, documenting their functional interaction. These findings identify a role for cohesin in regulating alternative splicing in both normal and leukemic cells and provide insights into the role of cohesin mutations in human disease.

Integrating transcription and splicing into cell fate: Transcription factors on the block

Transcription factors (TFs) are present in all life forms and conserved across great evolutionary distances in eukaryotes. From yeast to complex multicellular organisms, they are pivotal players of cell fate decision by orchestrating gene expression at diverse molecular layers. Notably, TFs fine-tune gene expression by coordinating RNA fate at both the expression and splicing levels. They regulate alternative splicing, an essential mechanism for cell plasticity, allowing the production of many mRNA and protein isoforms in precise cell and tissue contexts. Despite this apparent role in splicing, how TFs integrate transcription and splicing to ultimately orchestrate diverse cell functions and cell fate decisions remains puzzling. We depict substantial studies in various model organisms underlining the key role of TFs in alternative splicing for promoting tissue-specific functions and cell fate. Furthermore, we emphasize recent advances describing the molecular link between the transcriptional and splicing activities of TFs. As TFs can bind both DNA and/or RNA to regulate transcription and splicing, we further discuss their flexibility and compatibility for DNA and RNA substrates. Finally, we propose several models integrating transcription and splicing activities of TFs in the coordination and diversification of cell and tissue identities. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Processing > Splicing Mechanisms.

Identification of sex-specific splicing via comparative transcriptome analysis of gonads from sea cucumber Apostichopus japonicus

Alternative splicing (AS) is an essential post-transcriptional regulation mechanism for sex differentiation and gonadal development, which has rarely been reported in marine invertebrates. Sea cucumber (Apostichopus japonicus) is an economically important marine benthic echinoderm with a potential XX/XY sex determination mechanism, whose molecular mechanism in the gonadal differentiation has not been clearly understood. In this study, we analyzed available RNA-seq datasets of male and female gonads to explore if AS mechanism exerts an essential function in sex differentiation and gonadal development of A. japonicus. In our results, a total of 20,338 AS events from 7219 alternatively spliced genes, and 189 sexually differential alternative splicing (DAS) events from 156 genes were identified in gonadal transcriptome of sea cucumber. Gene Ontology analysis indicated that these DAS genes were significantly enriched in spermatogenesis-related GO terms. Maximal Clique Centrality (MCC) was then applied for protein-protein interaction (PPI) analysis to search for protein interactions and hub DAS gene. Among all DAS genes, we identified 10 DAS genes closely related to spermatogenesis and (or) sperm motility and a hub gene dnah1. Thus, this study revealed that alternative isoforms were generated from certain genes in female and male gonads through alternative splicing, which may provide direct evidence that alternative splicing mechanisms participate in female and male gonads. These results suggested a novel perspective for explaining the molecular mechanisms underlying gonadal differentiation between male and female sea cucumbers.

One locus, several functional RNAs-emerging roles of the mechanisms responsible for the sequence variability of microRNAs

With the development of modern molecular genetics, the original "one gene-one enzyme" hypothesis has been outdated. For protein coding genes, the discovery of alternative splicing and RNA editing provided the biochemical background for the RNA repertoire of a single locus, which also serves as an important pillar for the enormous protein variability of the genomes. Non-protein coding RNA genes were also revealed to produce several RNA species with distinct functions. The loci of microRNAs (miRNAs), encoding for small endogenous regulatory RNAs, were also found to produce a population of small RNAs, rather than a single defined product. This review aims to present the mechanisms contributing to the astonishing variability of miRNAs revealed by the new sequencing technologies. One important source is the careful balance of arm selection, producing sequentially different 5p- or 3p-miRNAs from the same pre-miRNA, thereby broadening the number of regulated target RNAs and the phenotypic response. In addition, the formation of 5', 3' and polymorphic isomiRs, with variable end and internal sequences also leads to a higher number of targeted sequences, and increases the regulatory output. These miRNA maturation processes, together with other known mechanisms such as RNA editing, further increase the potential outcome of this small RNA pathway. By discussing the subtle mechanisms behind the sequence diversity of miRNAs, this review intends to reveal this engaging aspect of the inherited "RNA world", how it contributes to the almost infinite molecular variability among living organisms, and how this variability can be exploited to treat human diseases.

Mendelian inheritance revisited: dominance and recessiveness in medical genetics

Understanding the consequences of genotype for phenotype (which ranges from molecule-level effects to whole-organism traits) is at the core of genetic diagnostics in medicine. Many measures of the deleteriousness of individual alleles exist, but these have limitations for predicting the clinical consequences. Various mechanisms can protect the organism from the adverse effects of functional variants, especially when the variant is paired with a wild type allele. Understanding why some alleles are harmful in the heterozygous state - representing dominant inheritance - but others only with the biallelic presence of pathogenic variants - representing recessive inheritance - is particularly important when faced with the deluge of rare genetic alterations identified by high throughput DNA sequencing. Both awareness of the specific quantitative and/or qualitative effects of individual variants and the elucidation of allelic and non-allelic interactions are essential to optimize genetic diagnosis and counselling.

hnRNPH1 establishes Sertoli-germ cell crosstalk through cooperation with PTBP1 and AR, and is essential for male fertility in mice

Spermatogenesis depends on the crosstalk of Sertoli cells (SCs) and germ cells. However, the gene regulatory network establishing the communications between SCs and germ cells remains unclear. Here, we report that heterogeneous nuclear ribonucleoprotein H1 (hnRNPH1) in SCs is essential for the establishment of crosstalk between SCs and germ cells. Conditional knockout of hnRNPH1 in mouse SCs leads to compromised blood-testis barrier function, delayed meiotic progression, increased germ cell apoptosis, sloughing of germ cells and, eventually, infertility of mice. Mechanistically, we discovered that hnRNPH1 could interact with the splicing regulator PTBP1 in SCs to regulate the pre-mRNA alternative splicing of the target genes functionally related to cell adhesion. Interestingly, we also found hnRNPH1 could cooperate with the androgen receptor, one of the SC-specific transcription factors, to modulate the transcription level of a group of genes associated with the cell-cell junction and EGFR pathway by directly binding to the gene promoters. Collectively, our findings reveal a crucial role for hnRNPH1 in SCs during spermatogenesis and uncover a potential molecular regulatory network involving hnRNPH1 in establishing Sertoli-germ cell crosstalk.

Identification of TSGA10 and GGNBP2 splicing variants in 5' untranslated region with distinct expression profiles in brain tumor samples

Introduction

Brain tumors (BTs) are perceived as one of the most common malignancies among children. The specific regulation of each gene can play a critical role in cancer progression. The present study aimed to determine the transcripts of the TSGA10 and GGNBP2 genes, considering the alternative 5'UTR region, and investigating the expression of these different transcripts in BTs.

Material and methods

Public data on brain tumor microarray datasets in GEO were analyzed with R software to evaluate the expression levels of TSGA10 and GGNBP2 genes (the Pheatmap package in R was also used to plot DEGs in a heat map). In addition, to validate our in-silico data analysis, RT-PCR was performed to determine the splicing variants of TSGA10 and GGNBP2 genes in testis and brain tumor samples. The expression levels of splice variants of these genes were analyzed in 30 brain tumor samples and two testicular tissue samples as a positive control.

Results

In silico results show that the differential expression levels of TSGA10 and GGNBP2 were significant in the GEO datasets of BTs compared to normal samples (with adjusted p-value<0.05 and log fold change > 1). This study's experimental results showed that the TSGA10 gene produces four different transcripts with two distinct promoter regions and splicing exon 4. The relative mRNA expression of transcripts without exon 4 was higher than transcripts with exon 4 in BT samples (p-value<001). In GGNBP2, exon 2 in the 5'UTR region and exon 6 in the coding sequence were spliced. The expression analysis results showed that the relative mRNA expression of transcript variants without exon 2 was higher than other transcript variants with exon 2 in BT samples (p-value<001).

Conclusion

The decreased expression levels of transcripts with longer 5'UTR in BT samples than in testicular or low-grade brain tumor samples may decrease their translation efficiency. Therefore, decreased amounts of TSGA10 and GGNBP2 as potential tumor suppressor proteins, especially in high-grade brain tumors, may cause cancer development by angiogenesis and metastasis.

The Pattern of RNA Editing Changes in Pleural Mesothelioma upon Epithelial-Mesenchymal Transition

Pleural mesothelioma (PM) is a cancer where epithelioid, biphasic and sarcomatoid histotypes are observed. Sarcomatoid PM is characterized by mesenchymal features. Multi-omics have been used to characterize the epithelial-to-mesenchymal (EMT) phenotype at the molecular level. We contribute to this effort by including the analysis of RNA editing. We extracted samples with the highest vs. lowest Epithelial score from two PM cohorts and observed increased RNA editing in introns and decreased RNA editing in 3'UTR upon EMT. The same was observed in primary PM primary cultures stratified by transcriptomics analysis into two groups, one of them enriched with mesenchymal features. Our data demonstrate that, as has been observed in other cancer types, RNA editing associates to EMT phenotype in PM.

Pancreatic microexons regulate islet function and glucose homeostasis

Pancreatic islets control glucose homeostasis by the balanced secretion of insulin and other hormones, and their abnormal function causes diabetes or hypoglycaemia. Here we uncover a conserved programme of alternative microexons included in mRNAs of islet cells, particularly in genes involved in vesicle transport and exocytosis. Islet microexons (IsletMICs) are regulated by the RNA binding protein SRRM3 and represent a subset of the larger neural programme that are particularly sensitive to SRRM3 levels. Both SRRM3 and IsletMICs are induced by elevated glucose levels, and depletion of SRRM3 in human and rat beta cell lines and mouse islets, or repression of particular IsletMICs using antisense oligonucleotides, leads to inappropriate insulin secretion. Consistently, mice harbouring mutations in Srrm3 display defects in islet cell identity and function, leading to hyperinsulinaemic hypoglycaemia. Importantly, human genetic variants that influence SRRM3 expression and IsletMIC inclusion in islets are associated with fasting glucose variation and type 2 diabetes risk. Taken together, our data identify a conserved microexon programme that regulates glucose homeostasis.

Transcriptome analyses in infertile men reveal germ cell-specific expression and splicing patterns

The process of spermatogenesis-when germ cells differentiate into sperm-is tightly regulated, and misregulation in gene expression is likely to be involved in the physiopathology of male infertility. The testis is one of the most transcriptionally rich tissues; nevertheless, the specific gene expression changes occurring during spermatogenesis are not fully understood. To better understand gene expression during spermatogenesis, we generated germ cell-specific whole transcriptome profiles by systematically comparing testicular transcriptomes from tissues in which spermatogenesis is arrested at successive steps of germ cell differentiation. In these comparisons, we found thousands of differentially expressed genes between successive germ cell types of infertility patients. We demonstrate our analyses' potential to identify novel highly germ cell-specific markers (TSPY4 and LUZP4 for spermatogonia; HMGB4 for round spermatids) and identified putatively misregulated genes in male infertility (RWDD2A, CCDC183, CNNM1, SERF1B). Apart from these, we found thousands of genes showing germ cell-specific isoforms (including SOX15, SPATA4, SYCP3, MKI67). Our approach and dataset can help elucidate genetic and transcriptional causes for male infertility.

RNA virus infections and their effect on host alternative splicing

It is evident that viral infection dramatically alters host gene expression, and these alterations have both pro- and anti-viral functions. While the effects of viral infection on transcription and translation have been comprehensively reviewed, less attention has been paid to the impact on alternative splicing of pre-messenger RNAs. Here we review salient examples of how viral infection leads to changes in alternative splicing and discuss how these changes impact infection.

Full-Length RNA Sequencing Provides Insights into Goldfish Evolution under Artificial Selection

Goldfish Carassius auratus is an ideal model for exploring fish morphology evolution. Although genes underlying several ornamental traits have been identified, little is known about the effects of artificial selection on embryo gene expression. In the present study, hybrid transcriptome sequencing was conducted to reveal gene expression profiles of Celestial-Eye (CE) and Ryukin (RK) goldfish embryos. Full-length transcriptome sequencing on the PacBio platform identified 54,218 and 54,106 transcript isoforms in CE and RK goldfish, respectively. Of particular note was that thousands of alternative splicing (AS) and alternative polyadenylation (APA) events were identified in both goldfish breeds, and most of them were inter-breed specific. RT-PCR and Sanger sequencing showed that most of the predicted AS and APA were correct. Moreover, abundant long non-coding RNA and fusion genes were detected, and again most of them were inter-breed specific. Through RNA-seq, we detected thousands of differentially expressed genes (DEGs) in each embryonic stage between the two goldfish breeds. KEGG enrichment analysis on DEGs showed extensive differences between CE and RK goldfish in gene expression. Taken together, our results demonstrated that artificial selection has led to far-reaching influences on goldfish gene expression, which probably laid the genetic basis for hundreds of goldfish variations.

All differential on the splicing front: Host alternative splicing alters the landscape of virus-host conflict

Alternative RNA splicing is a co-transcriptional process that richly increases proteome diversity, and is dynamically regulated based on cell species, lineage, and activation state. Virus infection in vertebrate hosts results in rapid host transcriptome-wide changes, and regulation of alternative splicing can direct a combinatorial effect on the host transcriptome. There has been a recent increase in genome-wide studies evaluating host alternative splicing during viral infection, which integrates well with prior knowledge on viral interactions with host splicing proteins. A critical challenge remains in linking how these individual events direct global changes, and whether alternative splicing is an overall favorable pathway for fending off or supporting viral infection. Here, we introduce the process of alternative splicing, discuss how to analyze splice regulation, and detail studies on genome-wide and splice factor changes during viral infection. We seek to highlight where the field can focus on moving forward, and how incorporation of a virus-host co-evolutionary perspective can benefit this burgeoning subject.

COL3A1: Potential prognostic predictor for head and neck cancer based on immune-microenvironment alternative splicing

We aimed to identify a novel prognostic biomarker for head and neck squamous cell carcinoma (HNSCC) based on tumor immunology-related alternative splicing (AS). Data for 502 HNSCC and 44 normal samples were obtained from the TCGA database and used to establish an AS-related risk model through univariate, least absolute shrinkage, and selection operator Cox regression analyses. Fresh HNSCC and normal oral tissues were surgically obtained from 44 HNSCC patients. Western blotting and quantitative reverse transcription-PCR were used to assess gene expression levels. Kaplan-Meier was performed to evaluate patients' overall survival (OS) rate. The CIBERSORT algorithm, single-sample gene set enrichment analysis, and immune checkpoint analyses were performed to compare immune activities between subgroups. The risk model was established using 10 pivotal AS events first. Collagen Type III Alpha 1 Chain (COL3A1) were screened based on |log2FC| ≥ 1 and FDR < 0.05 criteria. COL3A1 expression levels in HNSCC tissues were elevated relative to normal tissues (p < 0.001). Moreover, COL3A1 was a reliable biomarker for HNSCC patients' prognostic prediction in both cohorts (p < 0.001, p = 0.0085, respectively). COL3A1 protein (p = 0.0054) and mRNA (p < 0.0001) levels were correlated with HNSCC differentiation. Furthermore, the T stage was correlated with COL3A1 expression (p = 0.043), and COL3A1 expression was an independent prognostic predictor for HNSCC patients (p = 0.006). Compared with the risk model, COL3A1 was better at evaluating immune cell infiltrations, immune activities, and immune checkpoint gene expressions of HNSCC lesions.

Discovering riboswitches: the past and the future

Riboswitches are structured noncoding RNA domains used by many bacteria to monitor the concentrations of target ligands and regulate gene expression accordingly. In the past 20 years over 55 distinct classes of natural riboswitches have been discovered that selectively sense small molecules or elemental ions, and thousands more are predicted to exist. Evidence suggests that some riboswitches might be direct descendants of the RNA-based sensors and switches that were likely present in ancient organisms before the evolutionary emergence of proteins. We provide an overview of the current state of riboswitch research, focusing primarily on the discovery of riboswitches, and speculate on the major challenges facing researchers in the field.

Transcriptome-wide analysis reveals the coregulation of RNA-binding proteins and alternative splicing genes in the development of atherosclerosis

RNA-binding proteins (RBPs) are involved in the regulation of RNA splicing, stability, and localization. How RBPs control the development of atherosclerosis, is not fully understood. To explore the relevant RNA-binding proteins (RBPs) and alternative splicing events (ASEs) in atherosclerosis. We made a comprehensive work to integrate analyses of differentially expressed genes, including differential RBPs, and variable splicing characteristics related to different stages of atherosclerosis in dataset GSE104140. A total of 3712 differentially expressed genes (DEGs) were identified, including 2921 upregulated genes and 791 downregulated genes. Further analysis screened out 54 RBP genes, and 434 AS genes overlapped DEGs. We selected high expression ten RBP genes (SAMHD1, DDX60 L, TLR7, RBM47, MYEF2, RNASE6, PARP12, APOBEC3G, SMAD9, and RNASE1) for co-expression analysis. Meanwhile, we found seven regulated alternative splicing genes (RASGs) (ABI1, FXR1, CHID1, PLEC, PRKACB, BNIP2, PPP3CB) that could be regulated by RBPs. The co-expression network was used to further elucidate the regulatory and interaction relationship between RBPs and AS genes. Apoptotic process and innate immune response, revealed by the functional enrichment analysis of RASGs regulated by RBPs were closely related to atherosclerosis. In addition, 26 of the 344 alternative splicing genes regulated by the above 10 RBPs were transcription factors (TFs), We selected high expression nine TFs (TFDP1, RBBP7, STAT2, CREB5, ERG, ELF1, HMGN3, BCLAF1, and ZEB2) for co-expression analysis. The target genes of these TFs were mainly enriched in inflammatory and immune response pathways that were associated with atherosclerosis. indicating that AS abnormalities of these TFs may have a function in atherosclerosis. Furthermore, the expression of differentially expressed RBPs and the alternative splicing events of AS genes was validated by qRT-PCR in umbilical vein endothelial cells (HUVEC). The results showed that RBM47 were remarkedly difference in HUVEC treated with ox-LDL and the splicing ratio of AS in BCLAF1which is regulated by RBM47 significantly changed. In conclusion, the differentially expressed RBPs identified in our analysis may play important roles in the development of atherosclerosis by regulating the AS of these TF genes.

The internal ribosome entry site determines the neurotropic potential of enterovirus A71

The mechanisms underlying tissue-specific replication of enteroviruses are currently unclear. Although enterovirus A71 (EV-A71) and coxsackievirus A16 (CV-A16) are both common pathogens that cause hand-foot-mouth disease, they display quite different neurotropic properties. Herein, we characterized the role of the internal ribosomal entry site (IRES) in determining neurovirulence using an oral inoculation model of EV-A71. The receptor transgenic (hSCARB2-Tg) mice developed neurological symptoms after being infected with a mouse-adapted EV-A71 strain (MP4) via different administrative routes. Intragastric administration of the MP4 strain caused pathological changes such as neuronal loss and neuropil vacuolation, whereas replacing EV-A71 IRES with CV-A16 abolished the neuropathological phenotypes. The attenuated neurotropic potential of IRES-swapped EV-A71 was observed even in mice that received intraperitoneal and intracerebral inoculations. Fewer chimeric MP4 viral RNAs and proteins were detected in the mouse tissues, regardless of the inoculation route. Tissue-specific replication can be reflected in cell-based characterization. While chimeric MP4 virus replicated poorly in human intestinal C2BBe1 and neuroblastoma SH-SY5Y cells, its replication in susceptible rhabdomyosarcoma cells was not affected. Overall, our results demonstrated that the IRES determined the neurotropic potential of EV-A71 and CV-A16, emphasizing the importance of the IRES in tissue tropism, along with the host receptors.

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

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

PTBP1-activated co-transcriptional splicing controls epigenetic status of pluripotent stem cells

Many spliceosomal introns are excised from nascent transcripts emerging from RNA polymerase II (RNA Pol II). The extent of cell-type-specific regulation and possible functions of such co-transcriptional events remain poorly understood. We examined the role of the RNA-binding protein PTBP1 in this process using an acute depletion approach followed by the analysis of chromatin- and RNA Pol II-associated transcripts. We show that PTBP1 activates the co-transcriptional excision of hundreds of introns, a surprising effect given that this protein is known to promote intron retention. Importantly, some co-transcriptionally activated introns fail to complete their splicing without PTBP1. In a striking example, retention of a PTBP1-dependent intron triggers nonsense-mediated decay of transcripts encoding DNA methyltransferase DNMT3B. We provide evidence that this regulation facilitates the natural decline in DNMT3B levels in developing neurons and protects differentiation-specific genes from ectopic methylation. Thus, PTBP1-activated co-transcriptional splicing is a widespread phenomenon mediating epigenetic control of cellular identity.

Integrative analysis of Iso-Seq and RNA-seq reveals dynamic changes of alternative promoter, alternative splicing and alternative polyadenylation during Angiotensin II-induced senescence in rat primary aortic endothelial cells

In eukaryotes, alternative promoter (AP), alternative splicing (AS), and alternative polyadenylation (APA) are three crucial regulatory mechanisms that modulate message RNA (mRNA) diversity. Although AP, AS and APA are involved in diverse biological processess, whether they have dynamic changes in Angiotensin II (Ang II) induced senescence in rat primary aortic endothelial cells (RAECs), an important cellular model for studying cardiovascular disease, remains unclear. Here we integrated both PacBio single-molecule long-read isoform sequencing (Iso-Seq) and Illumina short-read RNA sequencing (RNA-seq) to analyze the changes of AP, AS and APA in Ang II-induced senescent RAECs. Iso-Seq generated 36,278 isoforms from 10,145 gene loci and 65.81% of these isoforms are novel, which were further cross-validated by public data obtained by other techonologies such as CAGE, PolyA-Seq and 3'READS. APA contributed most to novel isoforms, followed by AS and AP. Further investigation showed that AP, AS and APA could all contribute to the regulation of isoform, but AS has more dynamic changes compared to AP and APA upon Ang II stimulation. Genes undergoing AP, AS and APA in Ang II-treated cells are enriched in various pathways related to aging or senescence, suggesting that these molecular changes are involved in functional alterations during Ang II-induced senescence. Together, the present study largely improved the annotation of rat genome and revealed gene expression changes at isoform level, extending the understanding of the complexity of gene regulation in Ang II-treated RAECs, and also provided novel clues for discovering the regulatory mechanism undelying Ang II caused vascular senescence and diseases.

Comprehensive characteristics of pathological subtypes in testicular germ cell tumor: Gene expression, mutation and alternative splicing

Background

Testicular germ cell tumor (TGCT) is the most common tumor in young men, but molecular signatures, especially the alternative splicing (AS) between its subtypes have not yet been explored.

Methods

To investigate the differences between TGCT subtypes, we comprehensively analyzed the data of gene expression, alternative splicing (AS), and somatic mutation in TGCT patients from the TCGA database. The gene ontology (GO) enrichment analyses were used to explore the function of differentially expressed genes and spliced genes respectively, and Spearman correlation analysis was performed to explore the correlation between differential genes and AS events. In addition, the possible patterns in which AS regulates gene expression were elaborated by the ensemble database transcript atlas. And, we identified important transcription factors that regulate gene expression and AS and functionally validated them in TGCT cell lines.

Results

We found significant differences between expression and AS in embryonal carcinoma and seminoma, while mixed cell tumors were in between. GO enrichment analyses revealed that both differentially expressed and spliced genes were enriched in transcriptional regulatory pathways, and obvious correlation between expression and AS events was determined. By analyzing the transcript map and the sites where splicing occurs, we have demonstrated that AS regulates gene expression in a variety of ways. We further identified two pivot AS-related molecules (SOX2 and HDAC9) involved in AS regulation, which were validated in embryonal carcinoma and seminoma cell lines. Differences in somatic mutations between subtypes are also of concern, with our results suggesting that mutations in some genes (B3GNT8, CAPN7, FAT4, GRK1, TACC2, and TRAM1L1) occur only in embryonal carcinoma, while mutations in KIT, KARS, and NRAS are observed only in seminoma.

Conclusions

In conclusion, our analysis revealed the differences in gene expression, AS and somatic mutation among TGCT subtypes, providing a molecular basis for clinical diagnosis and precise therapy of TGCT patients.

Regionally restricted modulation of Sam68 expression and Arhgef9 alternative splicing in the hippocampus of a murine model of multiple sclerosis

Multiple sclerosis (MS) and its preclinical models are characterized by marked changes in neuroplasticity, including excitatory/inhibitory imbalance and synaptic dysfunction that are believed to underlie the progressive cognitive impairment (CI), which represents a significant clinical hallmark of the disease. In this study, we investigated several parameters of neuroplasticity in the hippocampus of the experimental autoimmune encephalomyelitis (EAE) SJL/J mouse model, characterized by rostral inflammatory and demyelinating lesions similar to Relapsing-Remitting MS. By combining morphological and molecular analyses, we found that the hippocampus undergoes extensive inflammation in EAE-mice, more pronounced in the CA3 and dentate gyrus (DG) subfields than in the CA1, associated with changes in GABAergic circuitry, as indicated by the increased expression of the interneuron marker Parvalbumin selectively in CA3. By laser-microdissection, we investigated the impact of EAE on the alternative splicing of Arhgef9, a gene encoding a post-synaptic protein playing an essential role in GABAergic synapses and whose mutations have been related to CI and epilepsy. Our results indicate that EAE induces a specific increase in inclusion of the alternative exon 11a only in the CA3 and DG subfields, in line with the higher local levels of inflammation. Consistently, we found a region-specific downregulation of Sam68, a splicing-factor that represses this splicing event. Collectively, our findings confirm a regionalized distribution of inflammation in the hippocampus of EAE-mice. Moreover, since neuronal circuit rearrangement and dynamic remodeling of structural components of the synapse are key processes that contribute to neuroplasticity, our study suggests potential new molecular players involved in EAE-induced hippocampal dysfunction.