Function of alternative splicing

Precision Transcriptome Editing

Manipulating RNA species in mammalian cells has emerged as an important strategy for precise gene expression control. Here we review recent advances in precision transcriptome editing with a focus on tools that engineer specific transcripts for abundance, translation, base editing, alternative isoforms, and chemical modifications. While some of these methods have demonstrated efficiency in therapeutically relevant cellular or in vivo models, most require further study on their clinical safety and efficacy. Precision transcriptome engineering holds great potential for both mechanistic study of RNA biology and future gene and cell-based therapeutic applications.

saseR: Juggling offsets unlocks RNA-seq tools for fast and Scalable differential usage, Aberrant Splicing and Expression Retrieval

RNA-seq data analysis relies on many different tools, each tailored to specific applications and coming with unique assumptions and restrictions. Indeed, tools for differential transcript usage, or diagnosing patients with rare diseases through splicing and expression outliers, either lack in performance, discard information, or do not scale to massive data compendia. Here, we show that replacing the normalisation offsets unlocks bulk RNA-seq workflows for scalable differential usage, aberrant splicing and expression analyses. Our method, saseR, is much faster than state-of-the-art methods, dramatically outperforms these to detect aberrant splicing, and provides a single workflow for various short- and long-read RNA-seq applications.

Alternative splicing perspective to prey preference of environmentally friendly biological agent Cryptolaemus montrouzieri

BACKGROUND

Cryptolaemus montrouzieri (Coccinellidae) is widely utilized as biological control agents in modern agriculture. A comprehensive understanding of its food preference can help guide mass rearing and safety management during field application of pest control. Although some studies have paid attentions to the impacts of prey shift on C. montrouzieri, little is known regarding the role of post-transcriptional regulations in its acclimation to unnatural preys.

RESULTS

We performed a genome-wide investigation on alternative splicing dynamics in C. montrouzieri in response to the predation transition from natural prey to unnatural ones. When feeding on undesired diets, 402-764 genes were differentially alternative spliced in C. montrouzieri. It is noteworthy that the majority of these genes (> 87%) were not differentially expressed, and these differentially spliced genes regulated distinct biological processes from differentially expressed genes, such as organ development and morphogenesis, locomotory behavior, and homeostasis processes. These suggested the functionally nonredendant role of alternative splicing in modulating physiological and metabolic responses of C. montrouzieri to the shift to undesired preys. In addition, the individuals feeding on aphids were subject to a lower level of changes in splicing than other alternative diets, which might be because of the similar chemical and microbial compositions. Our study further suggested a putative coupling of alternative splicing and nonsense-mediated decay (AS-NMD), which may play an important role in fine-tuning the protein repertoire of C. montrouzieri, and promoting its acclimation to predation changes.

CONCLUSION

These findings highlight the key role of alternative splicing in modulating the acclimation of ladybirds to prey shift and provide new genetic clues for the future application of ladybirds in biocontrol.

GWAS of multiple neuropathology endophenotypes identifies new risk loci and provides insights into the genetic risk of dementia

Genome-wide association studies (GWAS) have identified >80 Alzheimer's disease and related dementias (ADRD)-associated genetic loci. However, the clinical outcomes used in most previous studies belie the complex nature of underlying neuropathologies. Here we performed GWAS on 11 ADRD-related neuropathology endophenotypes with participants drawn from the following three sources: the National Alzheimer's Coordinating Center, the Religious Orders Study and Rush Memory and Aging Project, and the Adult Changes in Thought study (n = 7,804 total autopsied participants). We identified eight independent significantly associated loci, of which four were new (COL4A1, PIK3R5, LZTS1 and APOC2). Separately testing known ADRD loci, 19 loci were significantly associated with at least one neuropathology after false-discovery rate adjustment. Genetic colocalization analyses identified pleiotropic effects and quantitative trait loci. Methylation in the cerebral cortex at two sites near APOC2 was associated with cerebral amyloid angiopathy. Studies that include neuropathology endophenotypes are an important step in understanding the mechanisms underlying genetic ADRD risk.

PAIP1 binds to pre-mRNA and regulates alternative splicing of cancer pathway genes including VEGFA

BACKGROUND

Poly (A) binding protein interacting protein 1 (PAIP1) has been shown to causally contribute to the development and progression of cancer. However, the mechanisms of the PAIP1 regulation in tumor cells remain poorly understood.

RESULTS

Here, we used a recently developed UV cross-linking and RNA immunoprecipitation method (iRIP-seq) to map the direct and indirect interaction sites between PAIP1 and RNA on a transcriptome-wide level in HeLa cells. We found that PAIP1 not only binds to 3'UTRs, but also to pre-mRNAs/mRNAs with a strong bias towards the coding region and intron. PAIP1 binding sites are enriched in splicing enhancer consensus GA-rich motifs. RNA-seq analysis revealed that PAIP1 selectively modulates the alternative splicing of genes in some cancer hallmarks including cell migration, the mTOR signaling pathway and the HIF-1 signaling pathway. PAIP1-regulated alternative splicing events were strongly associated with PAIP1 binding, demonstrating that the binding may promote selection of the nearby splice sites. Deletion of a PAIP1 binding site containing seven repeats of GA motif reduced the PAIP1-mediated suppression of the exon 6 inclusion in a VEGFA mRNA isoform. Proteomic analysis of the PAIP1-interacted proteins revealed the enrichment of the spliceosome components and splicing factors.

CONCLUSIONS

These findings suggest that PAIP1 is both a polyadenylation and alternative splicing regulator, that may play a large role in RNA processing via its role in alternative splicing regulation.

Circular RNAs from bovine blastocysts can interact with miRNAs/tsRNAs from embryonic extracellular vesicles and regulate hatching

Circular RNAs (circRNAs) are endogenous biological macromolecules that regulate various biological processes including embryo development. However, little is known about which circRNAs are present in bovine preimplantation embryos and their respective roles. Here, we characterized the expression profile of circRNAs in bovine blastocysts for the first time. We detected 25,700 circRNAs in total, with 12,630 circRNAs uniquely expressed in blastocysts compared to degenerated embryos. CircRNA alternative splicing (AS) events were also found more frequently in blastocysts than in degenerated embryos (299 vs 258). Additionally, 410 circRNAs, among which 11 circRNAs with a high potential to encode polypeptides, were found differentially expressed between blastocysts and degenerated embryos. We further predicted and constructed a circRNA-miRNA-mRNA network, wherein differentially expressed circRNAs were shown to bind to bovine preimplantation embryo development-related miRNAs. Employing bioinformatic algorithms we found that differentially expressed circRNAs are associated with differentially expressed miRNAs and transfer RNA-derived small RNAs (tsRNAs) enclosed in embryonic extracellular vesicles (EVs). Furthermore, functional analysis revealed that knockdown of the evolutionarily conserved circAGO2 can inhibit blastocyst hatching. Overall, our study provides the first landscape of circRNAs in bovine preimplantation embryos and highlights the novel role of circRNAs as tsRNA binding partners influencing small RNA sorting and loading into EVs, with circAGO2 playing a regulatory role in bovine blastocyst hatching.

Role of long noncoding RNAs in the regulation of alternative splicing in glioblastoma

Glioblastoma multiforme (GBM) is a highly severe primary brain tumor. Despite extensive research, effective treatments remain elusive. Long noncoding RNAs (lncRNAs) play a significant role in both cancer and normal biology. They influence alternative splicing (AS), which is crucial in cancer. Advances in lncRNA-specific microarrays and next-generation sequencing have enhanced understanding of AS. Abnormal AS contributes to cancer invasion, metastasis, apoptosis, therapeutic resistance, and tumor development, including glioma. lncRNA-mediated AS affects several cellular signaling pathways, promoting or suppressing cancer malignancy. This review discusses the lncRNAs regulating AS in glioblastoma and their mechanisms.

The Liver-Protective Effects of the Essential Oil from Amomum villosum in Tilapia (Oreochromis niloticus): Antioxidant, Transcriptomic, and Metabolomic Modulations

This study investigates the effects of the essential oil from Amomum villosum (EOA) on liver-protective effects in Nile tilapia (Oreochromis niloticus), utilizing a multidisciplinary approach that integrates physiological assessments and transcriptomic and metabolomic analyses. Fish were fed diets containing 2 g/kg of EOA over a 56-day trial, with a no-EOA diet serving as the control. The results demonstrate that EOA supplementation improves liver histology, enhances antioxidant capacities, and reduces inflammation in tilapia. The transcriptomic analysis revealed significant alterations in gene expression profiles related to RNA splicing, metabolism, and disease pathways. The identification of differential genes and disease databases identified key target genes associated with the primary component of EOA for its anti-hepatobiliary disease effects. Furthermore, a molecular docking analysis of EOA major components with core differentially expressed genes in the hepatobiliary syndrome indicated that α-pinene is a potential Hsp90 inhibitor, which may prevent inflammation. A metabolomic analysis further demonstrated that EOA supplementation leads to notable changes in liver phospholipids, fatty acids, and carbohydrate metabolism. These findings underscore the potential of EOA as a natural additive for improving liver health in tilapia, offering valuable insights to the aquaculture industry for enhancing fish health and welfare in intensive farming systems.

Importance of Transcript Variants in Transcriptome Analyses

RNA sequencing (RNA-Seq) has become a widely adopted technique for studying gene expression. However, conventional RNA-Seq analyses rely on gene expression (GE) values that aggregate all the transcripts produced under a single gene identifier, overlooking the complexity of transcript variants arising from different transcription start sites or alternative splicing. Transcript variants may encode proteins with diverse functional domains, or noncoding RNAs. This study explored the implications of neglecting transcript variants in RNA-Seq analyses. Among the 1334 transcription factor (TF) genes expressed in mouse embryonic stem (ES) or trophoblast stem (TS) cells, 652 were differentially expressed in TS cells based on GE values (365 upregulated and 287 downregulated, ≥absolute 2-fold changes, false discovery rate (FDR) p-value ≤ 0.05). The 365 upregulated genes expressed 883 transcript variants. Further transcript expression (TE) based analyses identified only 174 (<20%) of the 883 transcripts to be upregulated. The remaining 709 transcripts were either downregulated or showed no significant changes. Meanwhile, the 287 downregulated genes expressed 856 transcript variants and only 153 (<20%) of the 856 transcripts were downregulated. The other 703 transcripts were either upregulated or showed no significant change. Additionally, the 682 insignificant TF genes (GE values < absolute 2-fold changes and/or FDR p-values > 0.05) between ES and TS cells expressed 2215 transcript variants. These included 477 (>21%) differentially expressed transcripts (276 upregulated and 201 downregulated, ≥absolute 2-fold changes, FDR p-value ≤ 0.05). Hence, GE based RNA-Seq analyses do not represent accurate expression levels due to divergent transcripts expression from the same gene. Our findings show that by including transcript variants in RNA-Seq analyses, we can generate a precise understanding of a gene's functional and regulatory landscape; ignoring the variants may result in an erroneous interpretation.

Predicting splicing patterns from the transcription factor binding sites in the promoter with deep learning

BACKGROUND

Alternative splicing is a pivotal mechanism of post-transcriptional modification that contributes to the transcriptome plasticity and proteome diversity in metazoan cells. Although many splicing regulations around the exon/intron regions are known, the relationship between promoter-bound transcription factors and the downstream alternative splicing largely remains unexplored.

RESULTS

In this study, we present computational approaches to unravel the regulatory relationship between promoter-bound transcription factor binding sites (TFBSs) and the splicing patterns. We curated a fine dataset that includes DNase I hypersensitive site sequencing and transcriptomes across fifteen human tissues from ENCODE. Specifically, we proposed different representations of TF binding context and splicing patterns to examine the associations between the promoter and downstream splicing events. While machine learning models demonstrated potential in predicting splicing patterns based on TFBS occupancies, the limitations in the generalization of predicting the splicing forms of singleton genes across diverse tissues was observed with carefully examination using different cross-validation methods. We further investigated the association between alterations in individual TFBS at promoters and shifts in exon splicing efficiency. Our results demonstrate that the convolutional neural network (CNN) models, trained on TF binding changes in the promoters, can predict the changes in splicing patterns. Furthermore, a systemic in silico substitutions analysis on the CNN models highlighted several potential splicing regulators. Notably, using empirical validation using K562 CTCFL shRNA knock-down data, we showed the significant role of CTCFL in splicing regulation.

CONCLUSION

In conclusion, our finding highlights the potential role of promoter-bound TFBSs in influencing the regulation of downstream splicing patterns and provides insights for discovering alternative splicing regulations.

Targeting alternative splicing of fibronectin in human renal proximal tubule epithelial cells with antisense oligonucleotides to reduce EDA+ fibronectin production and block an autocrine loop that drives renal fibrosis

TGFβ1 is a powerful regulator of fibrosis; secreted in a latent form, it becomes active after release from the latent complex. During tissue fibrosis, the EDA + isoform of cellular fibronectin is overexpressed. In pulmonary fibrosis it has been proposed that the fibronectin splice variant including an EDA domain (FN EDA+) activates latent TGFβ. Our work investigates the potential of blocking the 'splicing in' of EDA with antisense oligonucleotides to inhibit TGFβ1-induced EDA + fibronectin and to prevent the cascade of events initiated by TGFβ1 in human renal proximal tubule cells (PTEC). Human primary PTEC were treated with TGFβ1 for 48 h, medium removed and the cells transfected with RNase H-independent antisense oligonucleotides (ASO) designed to block EDA exon inclusion (ASO5). The efficacy of ASO to block EDA exon inclusion was assessed by EDA + fibronectin RNA and protein expression; the expression of TGFβ, αSMA (α smooth muscle actin), MMP2 (matrix metalloproteinse-2), MMP9 (matrix metalloproteinse-9), Collagen I, K Cadherin and connexin 43 was analysed. Targeting antisense oligonucleotides designed to block EDA exon inclusion in fibronectin pre mRNA were effective in reducing the amount of TGFβ1 -induced cellular EDA + fibronectin RNA and secreted EDA + fibronectin protein (assessed by western immunoblotting and immunocytochemistry) in human proximal tubule cells in an in vitro cell culture model. The effect was selective for EDA + exon with no effect on EDB + fibronectin RNA and total fibronectin mRNA. Exogenous TGFβ1 induced endogenous TGFβ, αSMA, MMP2, MMP9 and Col I mRNA. TGFβ1 treatment for 48h reduced the expression of K-Cadherin and increased the expression of connexin-43. These TGFβ1-induced pro-fibrotic changes were attenuated by ASO5 treatment. 48 h after the removal of exogenous TGFβ, further increases in αSMA, MMP2, MMP9 was observed; ASO5 significantly inhibited this subsequent increase. ASO5 treatment also significantly inhibited ability of the cell culture medium harvested at the end of the experiment (96h) to stimulate SMAD3 reporter cells. The role of endogenous TGFβ1 was confirmed by the use of a TGFβ receptor inhibitor. Our results demonstrate a critical role of FN EDA+ in a cycle of TGFβ driven pro-fibrotic responses in human PTEC and blocking its production with ASO technology offers a potential therapy to interrupt this vicious circle and hence limit the progression of renal fibrosis.

In silico and in cellulo approaches for functional annotation of human protein splice variants

The elegance of pre-mRNA splicing mechanisms continues to interest scientists even after over a half century, since the discovery of the fact that coding regions in genes are interrupted by non-coding sequences. The vast majority of human genes have several mRNA variants, coding structurally and functionally different protein isoforms in a tissue-specific manner and with a linkage to specific developmental stages of the organism. Alteration of splicing patterns shifts the balance of functionally distinct proteins in living systems, distorts normal molecular pathways, and may trigger the onset and progression of various pathologies. Over the past two decades, numerous studies have been conducted in various life sciences disciplines to deepen our understanding of splicing mechanisms and the extent of their impact on the functioning of living systems. This review aims to summarize experimental and computational approaches used to elucidate the functions of splice variants of a single gene based on our experience accumulated in the laboratory of interactomics of proteoforms at the Institute of Biomedical Chemistry (IBMC) and best global practices.

Genome-Wide Identification of a Maize Chitinase Gene Family and the Induction of Its Expression by Fusarium verticillioides (Sacc.) Nirenberg (1976) Infection

Maize chitinases are involved in chitin hydrolysis. Chitinases are distributed across various organisms including animals, plants, and fungi and are grouped into different glycosyl hydrolase families and classes, depending on protein structure. However, many chitinase functions and their interactions with other plant proteins remain unknown. The economic importance of maize (Zea mays L.) makes it relevant for studying the function of plant chitinases and their biological roles. This work aims to identify chitinase genes in the maize genome to study their gene structure, family/class classification, cis-related elements, and gene expression under biotic stress, such as Fusarium verticillioides infection. Thirty-nine chitinase genes were identified and found to be distributed in three glycosyl hydrolase (GH) families (18, 19 and 20). Likewise, the conserved domains and motifs were identified in each GH family member. The identified cis-regulatory elements are involved in plant development, hormone response, defense, and abiotic stress response. Chitinase protein-interaction network analysis predicted that they interact mainly with cell wall proteins. qRT-PCR analysis confirmed in silico data showing that ten different maize chitinase genes are induced in the presence of F. verticillioides, and that they could have several roles in pathogen infection depending on chitinase structure and cell wall localization.

Identification of Alternatively Spliced Novel Isoforms of Human HSPB8 Gene

HSPB8 is a heat shock protein belonging to a family of ATP-independent stress proteins called HSPB which are present far and wide in the cells of various organisms. They are committed to protein quality control (PQC) and strive to avert protein aggregation and to procreate a pool of non-native proteins that can be swiftly folded. Their fundamental expression or stress inducibility is regulated by various cis-elements localized in the HSPB regulatory regions. In the current study we have predicted and confirmed two alternatively spliced novel transcripts of HSPB8 gene in liver, brain, and heart. These spliced variants have smaller sizes owing to smaller N terminal regions and showed remarkable changes in their cellular localization. Novel isoform (HSPB8-N1) was predicted to be majorly localized to nuclear region while the reported isoform (HSPB8) and one of the novel isoforms (HSPB8-N2) were predicted to be cytoplasmic in nature. There were many changes observed in the phosphorylation sites of the novel isoforms as well. The newly reported isoforms lack several structural motifs that are essential for various functional endeavors of the HSPB8 protein. In silico analysis of the conceptually translated protein was carried out using various bioinformatics tools to gain an understanding of their properties in order to explore their possible potential in therapeutics.

Decoding RNA Metabolism by RNA-linked CRISPR Screening in Human Cells

RNAs undergo a complex choreography of metabolic processes in human cells that are regulated by thousands of RNA-associated proteins. While the effects of individual RNA-associated proteins on RNA metabolism have been extensively characterized, the full complement of regulators for most RNA metabolic events remain unknown. Here we present a massively parallel RNA-linked CRISPR (ReLiC) screening approach to measure the responses of diverse RNA metabolic events to knockout of 2,092 human genes encoding all known RNA-associated proteins. ReLiC screens highlight modular interactions between gene networks regulating splicing, translation, and decay of mRNAs. When combined with biochemical fractionation of polysomes, ReLiC reveals striking pathway-specific coupling between growth fitness and mRNA translation. Perturbing different components of the translation and proteostasis machineries have distinct effects on ribosome occupancy, while perturbing mRNA transcription leaves ribosome occupancy largely intact. Isoform-selective ReLiC screens capture differential regulation of intron retention and exon skipping by SF3b complex subunits. Chemogenomic screens using ReLiC decipher translational regulators upstream of mRNA decay and uncover a role for the ribosome collision sensor GCN1 during treatment with the anti-leukemic drug homoharringtonine. Our work demonstrates ReLiC as a versatile platform for discovering and dissecting regulatory principles of human RNA metabolism.

Characterization of alternative splicing during mammalian brain development reveals the extent of isoform diversity and potential effects on protein structural changes

Regulation of gene expression is critical for fate commitment of stem and progenitor cells during tissue formation. In the context of mammalian brain development, a plethora of studies have described how changes in the expression of individual genes characterize cell types across ontogeny and phylogeny. However, little attention has been paid to the fact that different transcripts can arise from any given gene through alternative splicing (AS). Considered a key mechanism expanding transcriptome diversity during evolution, assessing the full potential of AS on isoform diversity and protein function has been notoriously difficult. Here, we capitalize on the use of a validated reporter mouse line to isolate neural stem cells, neurogenic progenitors and neurons during corticogenesis and combine the use of short- and long-read sequencing to reconstruct the full transcriptome diversity characterizing neurogenic commitment. Extending available transcriptional profiles of the mammalian brain by nearly 50,000 new isoforms, we found that neurogenic commitment is characterized by a progressive increase in exon inclusion resulting in the profound remodeling of the transcriptional profile of specific cortical cell types. Most importantly, we computationally infer the biological significance of AS on protein structure by using AlphaFold2, revealing how radical protein conformational changes can arise from subtle changes in isoforms sequence. Together, our study reveals that AS has a greater potential to impact protein diversity and function than previously thought, independently from changes in gene expression.

Aging atlas reveals cell-type-specific effects of pro-longevity strategies

Organismal aging involves functional declines in both somatic and reproductive tissues. Multiple strategies have been discovered to extend lifespan across species. However, how age-related molecular changes differ among various tissues and how those lifespan-extending strategies slow tissue aging in distinct manners remain unclear. Here we generated the transcriptomic Cell Atlas of Worm Aging (CAWA, http://mengwanglab.org/atlas ) of wild-type and long-lived strains. We discovered cell-specific, age-related molecular and functional signatures across all somatic and germ cell types. We developed transcriptomic aging clocks for different tissues and quantitatively determined how three different pro-longevity strategies slow tissue aging distinctively. Furthermore, through genome-wide profiling of alternative polyadenylation (APA) events in different tissues, we discovered cell-type-specific APA changes during aging and revealed how these changes are differentially affected by the pro-longevity strategies. Together, this study offers fundamental molecular insights into both somatic and reproductive aging and provides a valuable resource for in-depth understanding of the diversity of pro-longevity mechanisms.

Surveying the landscape of RNA isoform diversity and expression across 9 GTEx tissues using long-read sequencing data

Even though alternative RNA splicing was discovered nearly 50 years ago (1977), we still understand very little about most isoforms arising from a single gene, including in which tissues they are expressed and if their functions differ. Human gene annotations suggest remarkable transcriptional complexity, with approximately 252,798 distinct RNA isoform annotations from 62,710 gene bodies (Ensembl v109; 2023), emphasizing the need to understand their biological effects. For example, 256 gene bodies have ≥50 annotated isoforms and 30 have ≥100, where one protein-coding gene (MAPK10) even has 192 distinct RNA isoform annotations. Whether such isoform diversity results from biological redundancy or spurious alternative splicing (i.e., noise), or whether individual isoforms have specialized functions (even if subtle) remains a mystery for most genes. Recent studies by Aguzzoli-Heberle et al., Leung et al., and Glinos et al. demonstrated long-read RNAseq enables improved RNA isoform quantification for essentially any tissue, cell type, or biological condition (e.g., disease, development, aging, etc.), making it possible to better assess individual isoform expression and function. While each study provided important discoveries related to RNA isoform diversity, deeper exploration is needed. We sought to quantify and characterize real isoform usage across tissues (compared to annotations). We used long-read RNAseq data from 58 GTEx samples across nine tissues (three brain, two heart, muscle, lung, liver, and cultured fibroblasts) generated by Glinos et al. and found considerable isoform diversity within and across tissues. Cerebellar hemisphere was the most transcriptionally complex tissue (22,522 distinct isoforms; 3,726 unique); liver was least diverse (12,435 distinct isoforms; 1,039 unique). We highlight gene clusters exhibiting high tissue-specific isoform diversity per tissue (e.g., TPM1 expresses 19 in heart's atrial appendage). We also validated 447 of the 700 new isoforms discovered by Aguzzoli-Heberle et al. and found that 88 were expressed in all nine tissues, while 58 were specific to a single tissue. This study represents a broad survey of the RNA isoform landscape, demonstrating isoform diversity across nine tissues and emphasizes the need to better understand how individual isoforms from a single gene body contribute to human health and disease.

Factors affecting the expression and stability of full-length and truncated SRSF3 proteins in human cancer cells

Alternative splicing plays a crucial role in increasing the diversity of mRNAs expressed in the genome. Serine/arginine-rich splicing factor 3 (SRSF3) is responsible for regulating the alternative splicing of its own mRNA and ensuring that its expression is balanced to maintain homeostasis. Moreover, the exon skipping of SRSF3 leads to the production of a truncated protein instead of a frameshift mutation that generates a premature termination codon (PTC). However, the precise regulatory mechanism involved in the splicing of SRSF3 remains unclear. In this study, we first established a platform for coexpressing full-length SRSF3 (SRSF3-FL) and SRSF3-PTC and further identified a specific antibody against the SRSF3-FL and truncated SRSF3 (SRSF3-TR) proteins. Next, we found that exogenously overexpressing SRSF3-FL or SRSF3-PTC failed to reverse the effects of digoxin, caffeine, or both in combination on this molecule and its targets. Endoplasmic reticulum-related pathways, transcription factors, and chemicals such as palmitic acid and phosphate were found to be involved in the regulation of SRSF3 expression. The downregulation of SRSF3-FL by palmitic acid and phosphate was mediated via different regulatory mechanisms in HeLa cells. In summary, we provide new insights into the altered expression of the SRSF3-FL and SRSF3-TR proteins for the identification of the functions of SRSF3 in cells.

RNA-seq analysis reveals transcriptome reprogramming and alternative splicing during early response to salt stress in tomato root

Salt stress is one of the dominant abiotic stress conditions that cause severe damage to plant growth and, in turn, limiting crop productivity. It is therefore crucial to understand the molecular mechanism underlying plant root responses to high salinity as such knowledge will aid in efforts to develop salt-tolerant crops. Alternative splicing (AS) of precursor RNA is one of the important RNA processing steps that regulate gene expression and proteome diversity, and, consequently, many physiological and biochemical processes in plants, including responses to abiotic stresses like salt stress. In the current study, we utilized high-throughput RNA-sequencing to analyze the changes in the transcriptome and characterize AS landscape during the early response of tomato root to salt stress. Under salt stress conditions, 10,588 genes were found to be differentially expressed, including those involved in hormone signaling transduction, amino acid metabolism, and cell cycle regulation. More than 700 transcription factors (TFs), including members of the MYB, bHLH, and WRKY families, potentially regulated tomato root response to salt stress. AS events were found to be greatly enhanced under salt stress, where exon skipping was the most prevalent event. There were 3709 genes identified as differentially alternatively spliced (DAS), the most prominent of which were serine/threonine protein kinase, pentatricopeptide repeat (PPR)-containing protein, E3 ubiquitin-protein ligase. More than 100 DEGs were implicated in splicing and spliceosome assembly, which may regulate salt-responsive AS events in tomato roots. This study uncovers the stimulation of AS during tomato root response to salt stress and provides a valuable resource of salt-responsive genes for future studies to improve tomato salt tolerance.

Landscape analysis of alternative splicing in kidney renal clear cell carcinoma and their clinical significance

A growing number of studies reveal that alternative splicing (AS) is associated with tumorigenesis, progression, and metastasis. Systematic analysis of alternative splicing signatures in renal cancer is lacking. In our study, we investigated the AS landscape of kidney renal clear cell carcinoma (KIRC) and identified AS predictive model to improve the prognostic prediction of KIRC. We obtained clinical data and gene expression profiles of KIRC patients from the TCGA database to evaluate AS events. The calculation results for seven types of AS events indicated that 46276 AS events from 10577 genes were identified. Next, we applied Cox regression analysis to identify 5864 prognostic-associated AS events. We used the Metascape database to verify the potential pathways of prognostic-associated AS. Moreover, we constructed KIRC prediction systems with prognostic-associated AS events by the LASSO Cox regression model. AUCs demonstrated that these prediction systems had excellent prognostic accuracy simultaneously. We identified 34 prognostic associated splicing factors (SFs) and constructed homologous regulatory networks. Furthermore, in vitro experiments were performed to validate the favorable effect of SFs FMR1 in KIRC. In conclusion, we overviewed AS events in KIRC and identified AS-based prognostic models to assist the survival prediction of KIRC patients. Our study may provide a novel predictive signature to improve the prognostic prediction of KIRC, which might facilitate KIRC patient counseling and individualized management.

Potential of Cell-Penetrating Peptide-Conjugated Antisense Oligonucleotides for the Treatment of SMA

Spinal muscular atrophy (SMA) is a severe neuromuscular disorder that is caused by mutations in the survival motor neuron 1 (SMN1) gene, hindering the production of functional survival motor neuron (SMN) proteins. Antisense oligonucleotides (ASOs), a versatile DNA-like drug, are adept at binding to target RNA to prevent translation or promote alternative splicing. Nusinersen is an FDA-approved ASO for the treatment of SMA. It effectively promotes alternative splicing in pre-mRNA transcribed from the SMN2 gene, an analog of the SMN1 gene, to produce a greater amount of full-length SMN protein, to compensate for the loss of functional protein translated from SMN1. Despite its efficacy in ameliorating SMA symptoms, the cellular uptake of these ASOs is suboptimal, and their inability to penetrate the CNS necessitates invasive lumbar punctures. Cell-penetrating peptides (CPPs), which can be conjugated to ASOs, represent a promising approach to improve the efficiency of these treatments for SMA and have the potential to transverse the blood-brain barrier to circumvent the need for intrusive intrathecal injections and their associated adverse effects. This review provides a comprehensive analysis of ASO therapies, their application for the treatment of SMA, and the encouraging potential of CPPs as delivery systems to improve ASO uptake and overall efficiency.

E3 SUMO ligase SIZ1 splicing variants localize and function according to external conditions

SIZ1 (SAP and MIZ1) is a member of the Siz/PIAS-type RING family of E3 SUMO (small ubiquitin-related modifier) ligases that play key roles in growth, development, and stress responses in plant and animal systems. Nevertheless, splicing variants of SIZ1 have not yet been characterized. Here, we identified four splicing variants of Arabidopsis (Arabidopsis thaliana) SIZ1, which encode three different protein isoforms. The SIZ1 gene encodes an 873-amino acid (aa) protein. Among the four SIZ1 splicing variants (SSVs), SSV1 and SSV4 encode identical 885 aa proteins; SSV2 encodes an 832 aa protein; and SSV3 encodes an 884 aa protein. SSV2 mainly localized to the plasma membrane, whereas SIZ1, SSV1/SSV4, and SSV3 localized to the nucleus. Interestingly, SIZ1 and all SSVs exhibited similar E3 SUMO ligase activities and preferred SUMO1 and SUMO2 for their E3 ligase activity. Transcript levels of SSV2 were substantially increased by heat treatment, while those of SSV1, SSV3, and SSV4 transcripts were unaffected by various abiotic stresses. SSV2 directly interacted with and sumoylated cyclic nucleotide-gated ion channel 6 (CNGC6), a positive thermotolerance regulator, enhancing the stability of CNGC6. Notably, transgenic siz1-2 mutants expressing SSV2 exhibited greater heat stress tolerance than wild-type plants, whereas those expressing SIZ1 were sensitive to heat stress. Furthermore, transgenic cngc6 plants overaccumulating a mutated mCNGC6 protein (K347R, a mutation at the sumoylation site) were sensitive to heat stress, similar to the cngc6 mutants, while transgenic cngc6 plants overaccumulating CNGC6 exhibited restored heat tolerance. Together, we propose that alternative splicing is an important mechanism that regulates the function of SSVs during development or under adverse conditions, including heat stress.

Deep proteomics reveals incorporation of unedited proteins into mitochondrial protein complexes in Arabidopsis

The mitochondrial proteome consists of numerous types of proteins which either are encoded and synthesized in the mitochondria, or encoded in the cell nucleus, synthesized in the cytoplasm and imported into the mitochondria. Their synthesis in the mitochondria, but not in the nucleus, relies on the editing of the primary transcripts of their genes at defined sites. Here, we present an in-depth investigation of the mitochondrial proteome of Arabidopsis (Arabidopsis thaliana) and a public online platform for the exploration of the data. For the analysis of our shotgun proteomic data, an Arabidopsis sequence database was created comprising all available protein sequences from the TAIR10 and Araport11 databases, supplemented with sequences of proteins translated from edited and nonedited transcripts of mitochondria. Amino acid sequences derived from partially edited transcripts were also added to analyze proteins encoded by the mitochondrial genome. Proteins were digested in parallel with six different endoproteases to obtain maximum proteome coverage. The resulting peptide fractions were finally analyzed using liquid chromatography coupled to ion mobility spectrometry and tandem mass spectrometry. We generated a "deep mitochondrial proteome" of 4,692 proteins. 1,339 proteins assigned to mitochondria by the SUBA5 database (https://suba.live) accounted for >80% of the total protein mass of our fractions. The coverage of proteins by identified peptides was particularly high compared to single-protease digests, allowing the exploration of differential splicing and RNA editing events at the protein level. We show that proteins translated from nonedited transcripts can be incorporated into native mitoribosomes and the ATP synthase complex. We present a portal for the use of our data, based on "proteomaps" with directly linked protein data. The portal is available at www.proteomeexplorer.de.

RNA binding proteins as mediators of pathological cardiac remodeling

RNA binding proteins (RBPs) play a central in the post-transcriptional regulation of gene expression, which can account for up to 50% of all variations in protein expression within a cell. Following their binding to target RNAs, RBPs most typically confer changes in gene expression through modulation of alternative spicing, RNA stabilization/degradation, or ribosome loading/translation rate. All of these post-transcriptional regulatory processes have been shown to play a functional role in pathological cardiac remodeling, and a growing body of evidence is beginning to identify the mechanistic contribution of individual RBPs and their cardiac RNA targets. This review highlights the mechanisms of RBP-dependent post-transcriptional gene regulation in cardiomyocytes and fibroblasts and our current understanding of how RNA binding proteins functionally contribute to pathological cardiac remodeling.

IFDlong: an isoform and fusion detector for accurate annotation and quantification of long-read RNA-seq data

Advancements in long-read transcriptome sequencing (long-RNA-seq) technology have revolutionized the study of isoform diversity. These full-length transcripts enhance the detection of various transcriptome structural variations, including novel isoforms, alternative splicing events, and fusion transcripts. By shifting the open reading frame or altering gene expressions, studies have proved that these transcript alterations can serve as crucial biomarkers for disease diagnosis and therapeutic targets. In this project, we proposed IFDlong, a bioinformatics and biostatistics tool to detect isoform and fusion transcripts using bulk or single-cell long-RNA-seq data. Specifically, the software performed gene and isoform annotation for each long-read, defined novel isoforms, quantified isoform expression by a novel expectation-maximization algorithm, and profiled the fusion transcripts. For evaluation, IFDlong pipeline achieved overall the best performance when compared with several existing tools in large-scale simulation studies. In both isoform and fusion transcript quantification, IFDlong is able to reach more than 0.8 Spearman's correlation with the truth, and more than 0.9 cosine similarity when distinguishing multiple alternative splicing events. In novel isoform simulation, IFDlong can successfully balance the sensitivity (higher than 90%) and specificity (higher than 90%). Furthermore, IFDlong has proved its accuracy and robustness in diverse in-house and public datasets on healthy tissues, cell lines and multiple types of diseases. Besides bulk long-RNA-seq, IFDlong pipeline has proved its compatibility to single-cell long-RNA-seq data. This new software may hold promise for significant impact on long-read transcriptome analysis. The IFDlong software is available at https://github.com/wenjiaking/IFDlong.

Diviner uncovers hundreds of novel human (and other) exons though comparative analysis of proteins

Background

Eukaryotic genes are often composed of multiple exons that are stitched together by splicing out the intervening introns. These exons may be conditionally joined in different combinations to produce a collection of related, but distinct, mRNA transcripts. For protein-coding genes, these products of alternative splicing lead to production of related protein variants (isoforms) of a gene. Complete labeling of the protein-coding content of a eukaryotic genome requires discovery of mRNA encoding all isoforms, but it is impractical to enumerate all possible combinations of tissue, developmental stage, and environmental context; as a result, many true exons go unlabeled in genome annotations.

Results

One way to address the combinatoric challenge of finding all isoforms in a single organism A is to leverage sequencing efforts for other organisms - each time a new organism is sequenced, it may be under a new combination of conditions, so that a previously unobserved isoform may be sequenced. We present Diviner, a software tool that identifies previously undocumented exons in organisms by comparing isoforms across species. We demonstrate Diviner's utility by locating hundreds of novel exons in the genomes of human, mouse, and rat, as well as in the ferret genome. Further, we provide analyses supporting the notion that most of the new exons reported by Diviner are likely to be part of a true (but unobserved) isoform of the containing species.

Integrative metabolomic and transcriptomic analyses reveals the accumulation patterns of key metabolites associated with flavonoids and terpenoids of Gynostemma pentaphyllum (Thunb.) Makino

Gynostemma pentaphyllum (Thunb.) Makino (G. pentaphyllum) is a medicinal and edible plant with multiple functions of liver protection, anti-tumor, anti-inflammation, balancing blood sugar and blood lipids. The nutritional value of the G. pentaphyllum plant is mainly due to its rich variety of biologically active substances, such as flavonoids, terpenes and polysaccharides. In this study, we performed a comprehensive analysis combining metabolomics and root, stem and leaf transcriptomic data of G. pentaphyllum. We used transcriptomics and metabolomics data to construct a dynamic regulatory network diagram of G. pentaphyllum flavonoids and terpenoids, and screened the transcription factors involved in flavonoids and terpenoids, including basic helix-loop-helix (bHLH), myb-related, WRKY, AP2/ERF. Transcriptome analysis results showed that among the DEGs related to the synthesis of flavonoids and terpenoids, dihydroflavonol 4-reductase (DFR) and geranylgeranyl diphosphate synthases (GGPPS) were core genes. This study presents a dynamic image of gene expression in different tissues of G. pentaphyllum, elucidating the key genes and metabolites of flavonoids and terpenoids. This study is beneficial to a deeper understanding of the medicinal plants of G. pentaphyllum, and also provides a scientific basis for further regulatory mechanisms of plant natural product synthesis pathways and drug development.

Differential transcript usage analysis incorporating quantification uncertainty via compositional measurement error regression modeling

Differential transcript usage (DTU) occurs when the relative expression of multiple transcripts arising from the same gene changes between different conditions. Existing approaches to detect DTU often rely on computational procedures that can have speed and scalability issues as the number of samples increases. Here we propose a new method, CompDTU, that uses compositional regression to model the relative abundance proportions of each transcript that are of interest in DTU analyses. This procedure leverages fast matrix-based computations that make it ideally suited for DTU analysis with larger sample sizes. This method also allows for the testing of and adjustment for multiple categorical or continuous covariates. Additionally, many existing approaches for DTU ignore quantification uncertainty in the expression estimates for each transcript in RNA-seq data. We extend our CompDTU method to incorporate quantification uncertainty leveraging common output from RNA-seq expression quantification tool in a novel method CompDTUme. Through several power analyses, we show that CompDTU has excellent sensitivity and reduces false positive results relative to existing methods. Additionally, CompDTUme results in further improvements in performance over CompDTU with sufficient sample size for genes with high levels of quantification uncertainty, while also maintaining favorable speed and scalability. We motivate our methods using data from the Cancer Genome Atlas Breast Invasive Carcinoma data set, specifically using RNA-seq data from primary tumors for 740 patients with breast cancer. We show greatly reduced computation time from our new methods as well as the ability to detect several novel genes with significant DTU across different breast cancer subtypes.

Widespread variation in molecular interactions and regulatory properties among transcription factor isoforms

Most human Transcription factors (TFs) genes encode multiple protein isoforms differing in DNA binding domains, effector domains, or other protein regions. The global extent to which this results in functional differences between isoforms remains unknown. Here, we systematically compared 693 isoforms of 246 TF genes, assessing DNA binding, protein binding, transcriptional activation, subcellular localization, and condensate formation. Relative to reference isoforms, two-thirds of alternative TF isoforms exhibit differences in one or more molecular activities, which often could not be predicted from sequence. We observed two primary categories of alternative TF isoforms: "rewirers" and "negative regulators", both of which were associated with differentiation and cancer. Our results support a model wherein the relative expression levels of, and interactions involving, TF isoforms add an understudied layer of complexity to gene regulatory networks, demonstrating the importance of isoform-aware characterization of TF functions and providing a rich resource for further studies.

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

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

Transcript Isoform Diversity of Y Chromosome Ampliconic Genes of Great Apes Uncovered Using Long Reads and Telomere-to-Telomere Reference Genome Assemblies

Y chromosomes of great apes harbor Ampliconic Genes (YAGs)-multi-copy gene families (BPY2, CDY, DAZ, HSFY, PRY, RBMY, TSPY, VCY, and XKRY) that encode proteins important for spermatogenesis. Previous work assembled YAG transcripts based on their targeted sequencing but not using reference genome assemblies, potentially resulting in an incomplete transcript repertoire. Here we used the recently produced gapless telomere-to-telomere (T2T) Y chromosome assemblies of great ape species (bonobo, chimpanzee, human, gorilla, Bornean orangutan, and Sumatran orangutan) and analyzed RNA data from whole-testis samples for the same species. We generated hybrid transcriptome assemblies by combining targeted long reads (Pacific Biosciences), untargeted long reads (Pacific Biosciences) and untargeted short reads (Illumina)and mapping them to the T2T reference genomes. Compared to the results from the reference-free approach, average transcript length was more than two times higher, and the total number of transcripts decreased three times, improving the quality of the assembled transcriptome. The reference-based transcriptome assemblies allowed us to differentiate transcripts originating from different Y chromosome gene copies and from their non-Y chromosome homologs. We identified two sources of transcriptome diversity-alternative splicing and gene duplication with subsequent diversification of gene copies. For each gene family, we detected transcribed pseudogenes along with protein-coding gene copies. We revealed previously unannotated gene copies of YAGs as compared to currently available NCBI annotations, as well as novel isoforms for annotated gene copies. This analysis paves the way for better understanding Y chromosome gene functions, which is important given their role in spermatogenesis.

Comprehensive gene profiling of the metabolic landscape of humanized livers in mice

BACKGROUND & AIMS

The human liver transcriptome is complex and highly dynamic, e.g. one gene may produce multiple distinct transcripts, each with distinct posttranscriptional modifications. Direct knowledge of transcriptome dynamics, however, is largely obscured by the inaccessibility of the human liver to treatments and the insufficient annotation of the human liver transcriptome at transcript and RNA modification levels.

METHODS

We generated mice that carry humanized livers of identical genetic background and subjected them to representative metabolic treatments. We then analyzed the humanized livers with nanopore single-molecule direct RNA sequencing to determine the expression level, m6A modification and poly(A) tail length of all RNA transcript isoforms. Our system allows for the de novo annotation of human liver transcriptomes to reflect metabolic responses and for the study of transcriptome dynamics in parallel.

RESULTS

Our analysis uncovered a vast number of novel genes and transcripts. Our transcript-level analysis of human liver transcriptomes also identified a multitude of regulated metabolic pathways that were otherwise invisible using conventional short-read RNA sequencing. We revealed for the first time the dynamic changes in m6A and poly(A) tail length of human liver transcripts, many of which are transcribed from key metabolic genes. Furthermore, we performed comparative analyses of gene regulation between humans and mice, and between two individuals using the liver-specific humanized mice, revealing that transcriptome dynamics are highly species- and genetic background-dependent.

CONCLUSION

Our work revealed a complex metabolic response landscape of the human liver transcriptome and provides a novel resource to understand transcriptome dynamics of the human liver in response to physiologically relevant metabolic stimuli (https://caolab.shinyapps.io/human_hepatocyte_landscape/).

IMPACT AND IMPLICATIONS

Direct knowledge of the human liver transcriptome is currently very limited, hindering the overall understanding of human liver pathophysiology. We combined a liver-specific humanized mouse model and long-read direct RNA sequencing technology to establish a de novo annotation of the human liver transcriptome and identified a multitude of regulated metabolic pathways that were otherwise invisible using conventional technologies. The extensive regulatory information on human genes we provided could enable basic scientists to infer the pathological relevance of their genes of interest and physician scientists to better pinpoint the changes in metabolic networks underlying a specific pathophysiology.

Unique transcriptomes of sensory and non-sensory neurons: insights from Splicing Regulatory States

Alternative Splicing (AS) programs serve as instructive signals of cell type specificity, particularly within the brain, which comprises dozens of molecularly and functionally distinct cell types. Among them, retinal photoreceptors stand out due to their unique transcriptome, making them a particularly well-suited system for studying how AS shapes cell type-specific molecular functions. Here, we use the Splicing Regulatory State (SRS) as a novel framework to discuss the splicing factors governing the unique AS pattern of photoreceptors, and how this pattern may aid in the specification of their highly specialized sensory cilia. In addition, we discuss how other sensory cells with ciliated structures, for which data is much scarcer, also rely on specific SRSs to implement a proteome specialized in the detection of sensory stimuli. By reviewing the general rules of cell type- and tissue-specific AS programs, firstly in the brain and subsequently in specialized sensory neurons, we propose a novel paradigm on how SRSs are established and how they can diversify. Finally, we illustrate how SRSs shape the outcome of mutations in splicing factors to produce cell type-specific phenotypes that can lead to various human diseases.

Environmentally relevant concentrations of 2,3,7,8-TCDD induced inhibition of multicellular alternative splicing and transcriptional dysregulation

Alternative splicing (AS), found in approximately 95 % of human genes, significantly amplifies protein diversity and is implicated in disease pathogenesis when dysregulated. However, the precise involvement of AS in the toxic mechanisms induced by TCDD (2,3,7,8-tetrachlorodibenzo-p-dioxin) remains incompletely elucidated. This study conducted a thorough global AS analysis in six human cell lines following TCDD exposure. Our findings revealed that environmentally relevant concentration (0.1 nM) of TCDD significantly suppressed AS events in all cell types, notably inhibiting diverse splicing events and reducing transcript diversity, potentially attributed to modifications in the splicing patterns of the inhibitory factor family, particularly hnRNP. And we identified 151 genes with substantial AS alterations shared among these cell types, particularly enriched in immune and metabolic pathways. Moreover, TCDD induced cell-specific changes in splicing patterns and transcript levels, with increased sensitivity notably in THP-1 monocyte, potentially linked to aberrant expression of pivotal genes within the spliceosome pathway (DDX5, EFTUD2, PUF60, RBM25, SRSF1, and CRNKL1). This study extends our understanding of disrupted alternative splicing and its relation to the multisystem toxicity of TCDD. It sheds light on how environmental toxins affect post-transcriptional regulatory processes, offering a fresh perspective for toxicology and disease etiology investigations.

Biosurfer for systematic tracking of regulatory mechanisms leading to protein isoform diversity

Long-read RNA sequencing has shed light on transcriptomic complexity, but questions remain about the functionality of downstream protein products. We introduce Biosurfer, a computational approach for comparing protein isoforms, while systematically tracking the transcriptional, splicing, and translational variations that underlie differences in the sequences of the protein products. Using Biosurfer, we analyzed the differences in 32,799 pairs of GENCODE annotated protein isoforms, finding a majority (70%) of variable N-termini are due to the alternative transcription start sites, while only 9% arise from 5' UTR alternative splicing. Biosurfer's detailed tracking of nucleotide-to-residue relationships helped reveal an uncommonly tracked source of single amino acid residue changes arising from the codon splits at junctions. For 17% of internal sequence changes, such split codon patterns lead to single residue differences, termed "ragged codons". Of variable C-termini, 72% involve splice- or intron retention-induced reading frameshifts. We found an unusual pattern of reading frame changes, in which the first frameshift is closely followed by a distinct second frameshift that restores the original frame, which we term a "snapback" frameshift. We analyzed long read RNA-seq-predicted proteome of a human cell line and found similar trends as compared to our GENCODE analysis, with the exception of a higher proportion of isoforms predicted to undergo nonsense-mediated decay. Biosurfer's comprehensive characterization of long-read RNA-seq datasets should accelerate insights of the functional role of protein isoforms, providing mechanistic explanation of the origins of the proteomic diversity driven by the alternative splicing. Biosurfer is available as a Python package at https://github.com/sheynkman-lab/biosurfer.

An array of signal-specific MoYpd1 isoforms determines full virulence in the pathogenic fungus Magnaporthe oryzae

Magnaporthe oryzae is placed first on a list of the world's top ten plant pathogens with the highest scientific and economic importance. The locus MGG_07173 occurs only once in the genome of M. oryzae and encodes the phosphotransfer protein MoYpd1p, which plays an important role in the high osmolarity glycerol (HOG) signaling pathway for osmoregulation. Originating from this locus, at least three MoYPD1 isoforms are produced in a signal-specific manner. The transcript levels of these MoYPD1-isoforms were individually affected by external stress. Salt (KCI) stress raised MoYPD1_T0 abundance, whereas osmotic stress by sorbitol elevates MoYPD1_T1 levels. In line with this, signal-specific nuclear translocation of green fluorescent protein-fused MoYpd1p isoforms in response to stress was observed. Mutant strains that produce only one of the MoYpd1p isoforms are less virulent, suggesting a combination thereof is required to invade the host successfully. In summary, we demonstrate signal-specific production of MoYpd1p isoforms that individually increase signal diversity and orchestrate virulence in M. oryzae.

Hemizygous splicing variant in CNKSR2 results in X-linked intellectual developmental disorder

BACKGROUND

Intellectual disability (ID) refers to a childhood-onset neurodevelopmental disorder with a prevalence of approximately 1%-3%.

METHODS

We performed whole exome sequencing for the patient with ID. And the splicing variant we found was validated by minigene assay.

RESULTS

Here, we report a boy with ID caused by a variant of CNKSR2. His neurological examination revealed hypsarrhythmia via electroencephalography and a right temporal polar arachnoid cyst via brain magnetic resonance imaging. A novel splicing variant in the CNKSR2 gene (NM_014927.5, c.1657+1G>A) was discovered by exome sequencing. The variant caused a 166 bp intron retention between exons 14 and 15, which was validated by a minigene assay. The variant was not reported in public databases such as gnomAD and the Exome Aggregation Consortium.

CONCLUSIONS

The variant was predicted to be damaging to correct the translation of the CNKRS2 protein and was classified as likely pathogenic according to the ACMG guidelines.

Alternative Splicing Reveals Acute Stress Response of Litopenaeus vannamei at High Alkalinity

Alkalinity is regarded as one of the primary stressors for aquatic animals in saline-alkaline water. Alternative splicing (AS) can significantly increase the diversity of transcripts and play key roles in stress response; however, the studies on AS under alkalinity stress of crustaceans are still limited. In the present study, we devoted ourselves to the study of AS under acute alkalinity stress at control (50 mg/L) and treatment groups (350 mg/L) by RNA-seq in pacific white shrimp (Litopenaeus vannamei). We identified a total of 10,556 AS events from 4865 genes and 619 differential AS (DAS) events from 519 DAS genes in pacific white shrimp. Functional annotation showed that the DAS genes primarily involved in spliceosome. Five splicing factors (SFs), U2AF1, PUF60, CHERP, SR140 and SRSF2 were significantly up-regulated and promoted AS. Furthermore, alkalinity activated the Leukocyte transendothelial migration, mTOR signaling pathway and AMPK signaling pathway, which regulated MAPK1, EIF3B and IGFP-RP1 associated with these pathways. We also studied three SFs (HSFP1, SRSF2 and NHE-RF1), which underwent AS to form different transcript isoforms. The above results demonstrated that AS was a regulatory mechanism in pacific white shrimp in response to acute alkalinity stress. SFs played vital roles in AS of pacific white shrimp, such as HSFP1, SRSF2 and NHE-RF1. DAS genes were significantly modified in immunity of pacific white shrimp to cope with alkalinity stress. This is the first study on the response of AS to acute alkalinity stress, which provided scientific basis for AS mechanism of crustaceans response to alkalinity stress.

Accounting for isoform expression increases power to identify genetic regulation of gene expression

A core problem in genetics is molecular quantitative trait locus (QTL) mapping, in which genetic variants associated with changes in the molecular phenotypes are identified. One of the most-studied molecular QTL mapping problems is expression QTL (eQTL) mapping, in which the molecular phenotype is gene expression. It is common in eQTL mapping to compute gene expression by aggregating the expression levels of individual isoforms from the same gene and then performing linear regression between SNPs and this aggregated gene expression level. However, SNPs may regulate isoforms from the same gene in different directions due to alternative splicing, or only regulate the expression level of one isoform, causing this approach to lose power. Here, we examine a broader question: which genes have at least one isoform whose expression level is regulated by genetic variants? In this study, we propose and evaluate several approaches to answering this question, demonstrating that "isoform-aware" methods-those that account for the expression levels of individual isoforms-have substantially greater power to answer this question than standard "gene-level" eQTL mapping methods. We identify settings in which different approaches yield an inflated number of false discoveries or lose power. In particular, we show that calling an eGene if there is a significant association between a SNP and any isoform fails to control False Discovery Rate, even when applying standard False Discovery Rate correction. We show that similar trends are observed in real data from the GEUVADIS and GTEx studies, suggesting the possibility that similar effects are present in these consortia.

Diosgenin biosynthesis investigation in medicinal herb (Tribulus terrestris) by transcriptome analysis

Next-generation sequencing (NGS) has revolutionized the analysis of specific genes, pathways, and their regulation in various species. Tribulus terrestris L., an annual medicinal herb of Zygophyllaceae family, has gained significant attention due to its diverse medicinal properties, including anti-inflammatory, antimicrobial, and anti-cancer effects. Diosgenin, a steroidal saponin, is the major bioactive compound responsible for the medicinal importance of T. terrestris. However, there is a paucity of information regarding the genes involved in the diosgenin biosynthetic pathway in T. terrestris. To address this gap, this study aimed to identify candidate genes associated with diosgenin biosynthesis through whole transcriptome profiling. A total of ∼7.9 GB of data, comprising 482 million reads, was obtained and assembled into 148,871 unigenes. Subsequently, functional annotations were assigned to 50 % of the unigenes using sequence similarity searches against the NCBI non-redundant (NR), Uniprot, KEGG, Pfam, GO, and COG databases, primarily based on Gene Ontology and KEGG-KAAS pathways. The majority of unigenes associated with the biosynthesis of the steroidal diosgenin backbone exhibited up-regulation in the fruit, leaf, and root tissues, except the SQE gene in root. The differential expression of selected genes was further validated through quantitative real-time polymerase chain reaction (qRT-PCR). Additionally, the study identified 21,026 unigenes related to transcription factors and 15,551 unigenes containing simple sequence repeats (SSR). Notably, di-nucleotide SSR motifs exhibited a high repeat frequency. These findings greatly enhance our understanding of the diosgenin biosynthesis pathway and provide a basis for future research in molecular investigation and metabolic engineering, specifically for boosting diosgenin content.

AI-Assisted Processing Pipeline to Boost Protein Isoform Detection

Proteomics, the study of proteins within biological systems, has seen remarkable advancements in recent years, with protein isoform detection emerging as one of the next major frontiers. One of the primary challenges is achieving the necessary peptide and protein coverage to confidently differentiate isoforms as a result of the protein inference problem and protein false discovery rate estimation challenge in large data. In this chapter, we describe the application of artificial intelligence-assisted peptide property prediction for database search engine rescoring by Oktoberfest, an approach that has proven effective, particularly for complex samples and extensive search spaces, which can greatly increase peptide coverage. Further, it illustrates a method for increasing isoform coverage by the PickedGroupFDR approach that is designed to excel when applied on large data. Real-world examples are provided to illustrate the utility of the tools in the context of rescoring, protein grouping, and false discovery rate estimation. By implementing these cutting-edge techniques, researchers can achieve a substantial increase in both peptide and isoform coverage, thus unlocking the potential of protein isoform detection in their studies and shedding light on their roles and functions in biological processes.

Identification of Key Prognostic Alternative Splicing Events of Costimulatory Molecule-Related Genes in Colon Cancer

OBJECTIVE

This study aimed to explore the key alternative splicing events in costimulatory molecule-related genes in colon cancer and to determine their correlation with prognosis.

METHODS

Gene expression RNA-sequencing data, clinical data, and SpliceSeq data of colon cancer were obtained from The Cancer Genome Atlas. Differentially expressed alternative splicing events in genes were identified, Followed by correlation analysis of genes corresponding to differentially expressed alternative splicing events with costimulatory molecule-related genes. Survival analysis was conducted using differentially expressed alternative splicing events in these genes and a prognostic model was constructed. Functional enrichment, proteinprotein interaction network, and splicing factor analyses were performed.

RESULTS

In total, 6504 differentially expressed alternative splicing events in 3949 genes were identified between tumor and normal tissues. Correlation analysis revealed 3499 differentially expressed alternative splicing events in 2168 costimulatory molecule-related genes. Moreover, 328 differentially expressed alternative splicing events in 288 costimulatory molecule-related genes were associated with overall survival. The prognostic models constructed using these showed considerable power in predicting survival. The ubiquitin A-52 residue ribosomal protein fusion product 1 and ribosomal protein S9 were the hub nodes in the protein-protein interaction network. Furthermore, one splicing factor, splicing factor proline and glutamine-rich, was significantly associated with patient prognosis. Four splicing factor-alternative splicing pairs were obtained from four alternative splicing events in three genes: TBC1 domain family member 8 B, complement factor H, and mitochondrial fission 1.

CONCLUSION

The identified differentially expressed alternative splicing events of costimulatory molecule-related genes may be used to predict patient prognosis and immunotherapy responses in colon cancer.

Effects of miR-722 on gene expression and alternative splicing in the liver of half-smooth tongue sole after infection with Vibrio anguillarum

MicroRNAs play crucial roles in various biological processes, including but not limited to differentiation, development, disease, and immunity. However, their immunoregulatory roles in half-smooth tongue sole are lacking. Our previous studies indicated that miR-722 could target C5aR1 to modulate the complement pathway to alleviate inflammatory response and even affect the mortality after the bacterial infection with Vibrio anguillarum. Driven by the purpose of revealing the underlying mechanisms, in this study, we investigated the effects of miR-722 on the gene expression and alternative splicing (AS) in the liver of half-smooth tongue sole after Vibrio anguillarum infection, with the approach of miR-722 overexpression/silencing and subsequent RNA-seq. Among the different comparisons, the I group (miR-722 inhibitor and V. anguillarum) versus blank control (PBS) exhibited the highest number of differentially expressed genes (DEGs), suggesting that the immune response was overactivated after inhibiting the miR-722. In addition, enrichment analyses were performed to reveal the functions of DEGs and differential AS (DAS) genes, reflecting the enrichment of RNA splicing and immune-related pathways including NF-κB and T cell receptor signaling pathway. Comparing the M group (miR-722 mimic and V. anguillarum) with the negative control (random sequence and V. anguillarum), two immune-related genes, cd48 and mapk8, were differentially expressed, of which mapk8 was also differentially spliced, indicating their importance in the immune response. Furthermore, representative gene analysis was performed, suggesting their corresponding functional changes due to AS. To verify the RNA-seq data, quantitative real-time PCR was employed with twenty pairs of primers for DEGs and DAS events. Overall, our results demonstrated that miR-722 could mediate the transcriptome-wide changes of gene expression and AS in half-smooth tongue sole, and provided insights into the regulatory role of miR-722 in immune responses, laying the foundation for further functional analyses and practical applications in aquaculture.

Improved detection of aberrant splicing with FRASER 2.0 and the intron Jaccard index

Detection of aberrantly spliced genes is an important step in RNA-seq-based rare-disease diagnostics. We recently developed FRASER, a denoising autoencoder-based method that outperformed alternative methods of detecting aberrant splicing. However, because FRASER's three splice metrics are partially redundant and tend to be sensitive to sequencing depth, we introduce here a more robust intron-excision metric, the intron Jaccard index, that combines the alternative donor, alternative acceptor, and intron-retention signal into a single value. Moreover, we optimized model parameters and filter cutoffs by using candidate rare-splice-disrupting variants as independent evidence. On 16,213 GTEx samples, our improved algorithm, FRASER 2.0, called typically 10 times fewer splicing outliers while increasing the proportion of candidate rare-splice-disrupting variants by 10-fold and substantially decreasing the effect of sequencing depth on the number of reported outliers. To lower the multiple-testing correction burden, we introduce an option to select the genes to be tested for each sample instead of a transcriptome-wide approach. This option can be particularly useful when prior information, such as candidate variants or genes, is available. Application on 303 rare-disease samples confirmed the relative reduction in the number of outlier calls for a slight loss of sensitivity; FRASER 2.0 recovered 22 out of 26 previously identified pathogenic splicing cases with default cutoffs and 24 when multiple-testing correction was limited to OMIM genes containing rare variants. Altogether, these methodological improvements contribute to more effective RNA-seq-based rare diagnostics by drastically reducing the amount of splicing outlier calls per sample at minimal loss of sensitivity.

A novel splice variant of goat CPT1a gene and their diverse mRNA expression profiles

The Alternative splicing (AS) of Carnitine palmitoyltransferase 1a (CPT1a) and their expression profiles had never been illuminated in goats until now. Herein, a novel splice transcript in the CPT1a gene that is predicted to result in the skipping of exons 6-19 (CPT1a-sv1) has been isolated in addition to the full-length transcript in goats. The result of RT-PCR showed that CPT1a-sv1 is 606 bp in length and consists of 6 exons. A novel exon 6 was consisted of partial exon 5 and partial exon 19, compared to that in CPT1a. RT-qPCR analysis showed that the expression patterns of CPT1a and CPT1a-sv1 are spatially different. In both kid and adult goats, the CPT1a transcript is strongly expressed in the liver, spleen, lung, kidney, and brain tissues. However, CPT1a-sv1 has a strong tissue-specific expression pattern, with moderate RNA levels in the liver and brain of kids, while highly expressed in the liver and minimally expressed in the brain of adults. We observed two transcripts to be involved in brain development. These findings improve our understanding of the function of the CPT1a gene in goats and provide information on the molecular mechanism of AS events.

Cryo-EM of a heterogeneous biochemical fraction elucidates multiple protein complexes from a multicellular thermophilic eukaryote

Biomolecular complexes and their interactions govern cellular structure and function. Understanding their architecture is a prerequisite for dissecting the cell's inner workings, but their higher-order assembly is often transient and challenging for structural analysis. Here, we performed cryo-EM on a single, highly heterogeneous biochemical fraction derived from Chaetomium thermophilum cell extracts to visualize the biomolecular content of the multicellular eukaryote. After cryo-EM single-particle image processing, results showed that a simultaneous three-dimensional structural characterization of multiple chemically diverse biomacromolecules is feasible. Namely, the thermophilic, eukaryotic complexes of (a) ATP citrate-lyase, (b) Hsp90, (c) 20S proteasome, (d) Hsp60 and (e) UDP-glucose pyrophosphorylase were characterized. In total, all five complexes have been structurally dissected in a thermophilic eukaryote in a total imaged sample area of 190.64 μm2, and two, in particular, 20S proteasome and Hsp60, exhibit side-chain resolution features. The C. thermophilum Hsp60 near-atomic model was resolved at 3.46 Å (FSC = 0.143) and shows a hinge-like conformational change of its equatorial domain, highly similar to the one previously shown for its bacterial orthologue, GroEL. This work demonstrates that cryo-EM of cell extracts will greatly accelerate the structural analysis of cellular complexes and provide unprecedented opportunities to annotate architectures of biomolecules in a holistic approach.

Comparative analysis of alternative splicing events in skeletal muscle of different sheep

This paper aims to investigate the relationship between genes with alternative splicing (AS) events and breed-specific differences in muscle development in two breeds of sheep. RNA-seq was utilized to identify genes with AS between Small-tailed Han sheep and Dorset sheep. The gene lists of differentially spliced genes were identified, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were conducted on these genes. In this study, 299 genes with 356 AS indicated significant differences between two diffrerent breeds. There are differences in 31 genes with 35 AS. Cassette, alt5' and alt3' exhibited the highest levels of enrichment across various significant levels. GO and KEGG enrichment analysis demonstrated a significant correlation between Wnt, TGF-beta, Notch and MAPK signaling pathways and the development of muscle in sheep. These findings indicate that genes with AS are linked to variations in muscle development in sheep. These results offer significant scientific and practical implications for improving the quality of sheep meat.

Alternative splicing in Atlantic salmon head kidney and SHK-1 cell line during the Piscirickettsia salmonis infection: A comparative transcriptome survey

Piscirickettsia salmonis, an intracellular bacterium in salmon aquaculture, is a big challenge because it is responsible for 54.2% of Atlantic salmon mortalities. In recent years, the high relevance of Alternative Splicing (AS) as a molecular mechanism associated with infectious conditions and host-pathogen interaction processes, especially in host immune activation, has been observed. Several studies have highlighted the role of AS in the host's immune response during viral, bacterial, and endoparasite infection. In the present study, we evaluated AS transcriptome profiles during P. salmonis infection in the two most used study models, SHK-1 cell line and salmon head kidney tissue. First, the SHK-1 cell line was exposed to P. salmonis infection at 0-, 7-, and 14-days post-infection (dpi). Following, total RNA was extracted for Illumina sequencing. On the other hand, RNA-Seq datasets of Atlantic salmon head kidney infected with the same P. salmonis strayingwase used. For both study models, the highest number of differentially alternative splicing (DAS) events was observed at 7 dpi, 16,830 DAS events derived from 9213 DAS genes in SHK-1 cells, and 13,820 DAS events from 7684 DAS genes in salmon HK. Alternative first exon (AF) was the most abundant AS type in the three infection times analyzed, representing 31% in SHK-1 cells and 228.6 in salmon HK; meanwhile, mutually exclusive exon (MX) was the least abundant. Notably, functional annotation of DAS genes in SHK-1 cells infected with P. salmonis showed a high presence of genes related to nucleotide metabolism. In contrast, the salmon head kidney exhibited many GO terms associated with immune response. Our findings reported the role of AS during P. salmonis infection in Atlantic salmon. These studies would contribute to a better understanding of the molecular bases that support the pathogen-host interaction, evidencing the contribution of AS regulating the transcriptional host response.

Hibernation-Induced microRNA Expression Promotes Signaling Pathways and Cell Cycle Dysregulation in Ictidomys tridecemlineatus Cardiac Tissue

The thirteen-lined ground squirrel Ictidomys tridecemlineatus is a rodent that lives throughout the United States and Canada and uses metabolic rate depression to facilitate circannual hibernation which helps it survive the winter. Metabolic rate depression is the reorganization of cellular physiology and molecular biology to facilitate a global downregulation of nonessential genes and processes, which conserves endogenous fuel resources and prevents the buildup of waste byproducts. Facilitating metabolic rate depression requires a complex interplay of regulatory approaches, including post-transcriptional modes such as microRNA. MicroRNA are short, single-stranded RNA species that bind to mRNA transcripts and target them for degradation or translational suppression. Using next-generation sequencing, we analyzed euthermic vs. hibernating cardiac tissue in I. tridecemlineatus to predict seven miRNAs (let-7e-5p, miR-122-5p, miR-2355-3p, miR-6715b-3p, miR-378i, miR-9851-3p, and miR-454-3p) that may be differentially regulated during hibernation. Gene ontology and KEGG pathway analysis suggested that these miRNAs cause a strong activation of ErbB2 signaling which causes downstream effects, including the activation of MAPK and PI3K/Akt signaling and concurrent decreases in p53 signaling and cell cycle-related processes. Taken together, these results predict critical miRNAs that may change during hibernation in the hearts of I. tridecemlineatus and identify key signaling pathways that warrant further study in this species.