Widespread Expansion of Protein Interaction Capabilities by Alternative Splicing

A comprehensive overview of microbiome data in the light of machine learning applications: categorization, accessibility, and future directions

Metagenomics, Metabolomics, and Metaproteomics have significantly advanced our knowledge of microbial communities by providing culture-independent insights into their composition and functional potential. However, a critical challenge in this field is the lack of standard and comprehensive metadata associated with raw data, hindering the ability to perform robust data stratifications and consider confounding factors. In this comprehensive review, we categorize publicly available microbiome data into five types: shotgun sequencing, amplicon sequencing, metatranscriptomic, metabolomic, and metaproteomic data. We explore the importance of metadata for data reuse and address the challenges in collecting standardized metadata. We also, assess the limitations in metadata collection of existing public repositories collecting metagenomic data. This review emphasizes the vital role of metadata in interpreting and comparing datasets and highlights the need for standardized metadata protocols to fully leverage metagenomic data's potential. Furthermore, we explore future directions of implementation of Machine Learning (ML) in metadata retrieval, offering promising avenues for a deeper understanding of microbial communities and their ecological roles. Leveraging these tools will enhance our insights into microbial functional capabilities and ecological dynamics in diverse ecosystems. Finally, we emphasize the crucial metadata role in ML models development.

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.

Standigm ASK™: knowledge graph and artificial intelligence platform applied to target discovery in idiopathic pulmonary fibrosis

Standigm ASK™ revolutionizes healthcare by addressing the critical challenge of identifying pivotal target genes in disease mechanisms-a fundamental aspect of drug development success. Standigm ASK™ integrates a unique combination of a heterogeneous knowledge graph (KG) database and an attention-based neural network model, providing interpretable subgraph evidence. Empowering users through an interactive interface, Standigm ASK™ facilitates the exploration of predicted results. Applying Standigm ASK™ to idiopathic pulmonary fibrosis (IPF), a complex lung disease, we focused on genes (AMFR, MDFIC and NR5A2) identified through KG evidence. In vitro experiments demonstrated their relevance, as TGFβ treatment induced gene expression changes associated with epithelial-mesenchymal transition characteristics. Gene knockdown reversed these changes, identifying AMFR, MDFIC and NR5A2 as potential therapeutic targets for IPF. In summary, Standigm ASK™ emerges as an innovative KG and artificial intelligence platform driving insights in drug target discovery, exemplified by the identification and validation of therapeutic targets for IPF.

Predicting the Structural Impact of Human Alternative Splicing

Protein structure prediction with neural networks is a powerful new method for linking protein sequence, structure, and function, but structures have generally been predicted for only a single isoform of each gene, neglecting splice variants. To investigate the structural implications of alternative splicing, we used AlphaFold2 to predict the structures of more than 11,000 human isoforms. We employed multiple metrics to identify splicing-induced structural alterations, including template matching score, secondary structure composition, surface charge distribution, radius of gyration, accessibility of post-translational modification sites, and structure-based function prediction. We identified examples of how alternative splicing induced clear changes in each of these properties. Structural similarity between isoforms largely correlated with degree of sequence identity, but we identified a subset of isoforms with low structural similarity despite high sequence similarity. Exon skipping and alternative last exons tended to increase the surface charge and radius of gyration. Splicing also buried or exposed numerous post-translational modification sites, most notably among the isoforms of BAX. Functional prediction nominated numerous functional differences among isoforms of the same gene, with loss of function compared to the reference predominating. Finally, we used single-cell RNA-seq data from the Tabula Sapiens to determine the cell types in which each structure is expressed. Our work represents an important resource for studying the structure and function of splice isoforms across the cell types of the human body.

Improving the RNA velocity approach with single-cell RNA lifecycle (nascent, mature and degrading RNAs) sequencing technologies

We presented an experimental method called FLOUR-seq, which combines BD Rhapsody and nanopore sequencing to detect the RNA lifecycle (including nascent, mature, and degrading RNAs) in cells. Additionally, we updated our HIT-scISOseq V2 to discover a more accurate RNA lifecycle using 10x Chromium and Pacbio sequencing. Most importantly, to explore how single-cell full-length RNA sequencing technologies could help improve the RNA velocity approach, we introduced a new algorithm called 'Region Velocity' to more accurately configure cellular RNA velocity. We applied this algorithm to study spermiogenesis and compared the performance of FLOUR-seq with Pacbio-based HIT-scISOseq V2. Our findings demonstrated that 'Region Velocity' is more suitable for analyzing single-cell full-length RNA data than traditional RNA velocity approaches. These novel methods could be useful for researchers looking to discover full-length RNAs in single cells and comprehensively monitor RNA lifecycle in cells.

Systematic analysis of alternative exon-dependent interactome remodeling reveals multitasking functions of gene regulatory factors

Alternative splicing significantly expands biological complexity, particularly in the vertebrate nervous system. Increasing evidence indicates that developmental and tissue-dependent alternative exons often control protein-protein interactions; yet, only a minor fraction of these events have been characterized. Using affinity purification-mass spectrometry (AP-MS), we show that approximately 60% of analyzed neural-differential exons in proteins previously implicated in transcriptional regulation result in the gain or loss of interaction partners, which in some cases form unexpected links with coupled processes. Notably, a neural exon in Chtop regulates its interaction with the Prmt1 methyltransferase and DExD-Box helicases Ddx39b/a, affecting its methylation and activity in promoting RNA export. Additionally, a neural exon in Sap30bp affects interactions with RNA processing factors, modulating a critical function of Sap30bp in promoting the splicing of <100 nt "mini-introns" that control nuclear RNA levels. AP-MS is thus a powerful approach for elucidating the multifaceted functions of proteins imparted by context-dependent alternative exons.

Regulation of alternative splicing and polyadenylation in neurons

Cell-type-specific gene expression is a fundamental feature of multicellular organisms and is achieved by combinations of regulatory strategies. Although cell-restricted transcription is perhaps the most widely studied mechanism, co-transcriptional and post-transcriptional processes are also central to the spatiotemporal control of gene functions. One general category of expression control involves the generation of multiple transcript isoforms from an individual gene, whose balance and cell specificity are frequently tightly regulated via diverse strategies. The nervous system makes particularly extensive use of cell-specific isoforms, specializing the neural function of genes that are expressed more broadly. Here, we review regulatory strategies and RNA-binding proteins that direct neural-specific isoform processing. These include various classes of alternative splicing and alternative polyadenylation events, both of which broadly diversify the neural transcriptome. Importantly, global alterations of splicing and alternative polyadenylation are characteristic of many neural pathologies, and recent genetic studies demonstrate how misregulation of individual neural isoforms can directly cause mutant phenotypes.

ITPR1-AS1 promotes small cell lung cancer metastasis by facilitating P21HRAS splicing and stabilizing DDX3X to activate the cRaf-MEK-ERK cascade

The mechanisms underlying the involvement of long non-coding RNAs (lncRNAs) in the metastasis of small cell lung cancer (SCLC) remain largely unknown. Here, we identified that the lncRNA ITPR1-AS1 was upregulated in SCLC and lymph node metastasis tissues and positively correlated with SCLC malignant features. The overexpression of ITPR1-AS1 in SCLC was an independent risk factor for the overall survival of patients with SCLC. Our data confirmed that ITPR1-AS1 induces SCLC cell metastasis both in vitro and in vivo. Mechanistically, ITPR1-AS1 acts as a scaffold to enhance the interaction between SRC-associated in mitosis 68 kDa and heterogeneous nuclear ribonucleoprotein A1, which facilitates the alternative splicing of the H-Ras proto-oncogene (HRAS) pre-mRNA (P21HRAS). Moreover, we observed that ITPR1-AS1 could associate in a complex with and maintain the stability of DEAD-box polypeptide 3 (DDX3X), which inhibited the latter's ubiquitination and degradation. Our data provide evidence that ITPR1-AS1 activates the cRaf-MEK-ERK cascade by upregulating P21HRAS production and stabilizing DDX3X, to promote SCLC metastasis.

Extensive Alternative Splicing Patterns in Systemic Lupus Erythematosus Highlight Sexual Differences

Substantial evidence highlights divergences in immune responses between men and women. Women are more susceptible to autoimmunity, whereas men suffer from the more severe presentation of autoimmune disorders. The molecular mechanism of this sexual dimorphism remains elusive. Herein, we conducted a comprehensive analysis of sex differences in whole-blood gene expression focusing on alternative splicing (AS) events in systemic lupus erythematosus (SLE), which is a prototype sex-biased disease. This study included 79 SLE patients with active disease and 58 matched healthy controls who underwent whole-blood RNA sequencing. Sex differences in splicing events were widespread, existent in both SLE and a healthy state. However, we observed distinct gene sets and molecular pathways targeted by sex-dependent AS in SLE patients as compared to healthy subjects, as well as a notable sex dissimilarity in intron retention events. Sexually differential spliced genes specific to SLE patients were enriched for dynamic cellular processes including chromatin remodeling, stress and inflammatory responses. Remarkably, the extent of sexual differences in AS in the SLE patients and healthy individuals exceeded those in gene expression. Overall, this study reveals an unprecedent variation in sex-dependent splicing events in SLE and the healthy state, with potential implications for understanding the molecular basis of sexual dimorphism in autoimmunity.

Metabolic reprogramming of cancer cells by JMJD6-mediated pre-mRNA splicing is associated with therapeutic response to splicing inhibitor

Dysregulated pre-mRNA splicing and metabolism are two hallmarks of MYC-driven cancers. Pharmacological inhibition of both processes has been extensively investigated as potential therapeutic avenues in preclinical and clinical studies. However, how pre-mRNA splicing and metabolism are orchestrated in response to oncogenic stress and therapies is poorly understood. Here, we demonstrate that Jumonji Domain Containing 6, Arginine Demethylase and Lysine Hydroxylase, JMJD6, acts as a hub connecting splicing and metabolism in MYC-driven neuroblastoma. JMJD6 cooperates with MYC in cellular transformation by physically interacting with RNA binding proteins involved in pre-mRNA splicing and protein homeostasis. Notably, JMJD6 controls the alternative splicing of two isoforms of glutaminase (GLS), namely kidney-type glutaminase (KGA) and glutaminase C (GAC), which are rate-limiting enzymes of glutaminolysis in the central carbon metabolism in neuroblastoma. Further, we show that JMJD6 is correlated with the anti-cancer activity of indisulam, a "molecular glue" that degrades splicing factor RBM39, which complexes with JMJD6. The indisulam-mediated cancer cell killing is at least partly dependent on the glutamine-related metabolic pathway mediated by JMJD6. Our findings reveal a cancer-promoting metabolic program is associated with alternative pre-mRNA splicing through JMJD6, providing a rationale to target JMJD6 as a therapeutic avenue for treating MYC-driven cancers.

DASES: a database of alternative splicing for esophageal squamous cell carcinoma

Esophageal carcinoma ranks as the sixth leading cause of cancer-related mortality globally, with esophageal squamous cell carcinoma (ESCC) being particularly prevalent among Asian populations. Alternative splicing (AS) plays a pivotal role in ESCC development and progression by generating diverse transcript isoforms. However, the current landscape lacks a specialized database focusing on alternative splicing events (ASEs) derived from a large number of ESCC cases. Additionally, most existing AS databases overlook the contribution of long non-coding RNAs (lncRNAs) in ESCC molecular mechanisms, predominantly focusing on mRNA-based ASE identification. To address these limitations, we deployed DASES (http://www.hxdsjzx.cn/DASES). Employing a combination of publicly available and in-house ESCC RNA-seq datasets, our extensive analysis of 346 samples, with 93% being paired tumor and adjacent non-tumor tissues, led to the identification of 257 novel lncRNAs in esophageal squamous cell carcinoma. Leveraging a paired comparison of tumor and adjacent normal tissues, DASES identified 59,094 ASEs that may be associated with ESCC. DASES fills a critical gap by providing comprehensive insights into ASEs in ESCC, encompassing lncRNAs and mRNA, thus facilitating a deeper understanding of ESCC molecular mechanisms and serving as a valuable resource for ESCC research communities.

HnRNPR strongly represses splicing of a critical exon associated with spinal muscular atrophy through binding to an exonic AU-rich element

BACKGROUND

Spinal muscular atrophy (SMA) is a motor neuron disease caused by mutations of survival of motor neuron 1 (SMN1) gene, which encodes the SMN protein. SMN2, a nearly identical copy of SMN1, with several single-nucleotide substitutions leading to predominant skipping of its exon 7, is insufficient to compensate for loss of SMN1. Heterogeneous nuclear ribonucleoprotein R (hnRNPR) has been previously shown to interact with SMN in the 7SK complex in motoneuron axons and is implicated in the pathogenesis of SMA. Here, we show that hnRNPR also interacts with SMN1/2 pre-mRNAs and potently inhibits exon 7 inclusion.

METHODS

In this study, to examine the mechanism that hnRNPR regulates SMN1/2 splicing, deletion analysis in an SMN2 minigene system, RNA-affinity chromatography, co-overexpression analysis and tethering assay were performed. We screened antisense oligonucleotides (ASOs) in a minigene system and identified a few that markedly promoted SMN2 exon 7 splicing.

RESULTS

We pinpointed an AU-rich element located towards the 3' end of the exon that mediates splicing repression by hnRNPR. We uncovered that both hnRNPR and Sam68 bind to the element in a competitive manner, and the inhibitory effect of hnRNPR is much stronger than Sam68. Moreover, we found that, among the four hnRNPR splicing isoforms, the exon 5-skipped one has the minimal inhibitory effect, and ASOs inducing hnRNPR exon 5 skipping also promote SMN2 exon 7 inclusion.

CONCLUSION

We identified a novel mechanism that contributes to mis-splicing of SMN2 exon 7.

Bioinformatic analysis of smoothelin family members supports tissue-specific functions of unique C-terminal calponin homology domains

Smoothelins are cytoskeletal proteins with a single C-terminal calponin homology domain type 2 (CHD2). Little is known about the significance of variation in SMTN CHD2 domains, addressed here through analysis of public databases. A conserved 152 nt penultimate constitutive exon present in all SMTNs encodes helices II-IV of CHD2 with high identity (nt/aa 63/65%). Variable CHD2s of SMTN (helices IV-VI) are generated by alternative splicing of 165 nt exon E20. E20 and the CHD2 it encodes have high homology with the terminal constitutive exon of SMTNL1 (E8; nt/aa 72/75% identity). Unique to these CHD2 variants are a conserved extended nine amino acid C-terminal tail containing KTKK ubiquitination motifs. When E20 of SMTN is skipped (SMTN E20-), constitutive terminal E21 codes for helices IV-VI of CHD2. SMTN E21 has high identity with the terminal exon of SMTNL2 (E8; nt/aa 75/81% identity of aligned sequences) except for coding for a unique extended C-terminus (24 nt; 8aa) conserved only in mammals. SMTN isoform expression is tissue-specific: SMTNE20- and SMTNE20+ are highly expressed in SMC and non-muscle cells, respectively, while SMTNL1 + 2 are highly expressed in skeletal muscle cells. Tissue-specific expression of SMTN CHD2s with unique helices IV-VI suggest tissue-specific functions that require further study.

Expression and splicing mediate distinct biological signals

BACKGROUND

Through alternative splicing, most human genes produce multiple isoforms in a cell-, tissue-, and disease-specific manner. Numerous studies show that alternative splicing is essential for development, diseases, and their treatments. Despite these important examples, the extent and biological relevance of splicing are currently unknown.

RESULTS

To solve this problem, we developed pairedGSEA and used it to profile transcriptional changes in 100 representative RNA-seq datasets. Our systematic analysis demonstrates that changes in splicing, on average, contribute to 48.1% of the biological signal in expression analyses. Gene-set enrichment analysis furthermore indicates that expression and splicing both convey shared and distinct biological signals.

CONCLUSIONS

These findings establish alternative splicing as a major regulator of the human condition and suggest that most contemporary RNA-seq studies likely miss out on critical biological insights. We anticipate our results will contribute to the transition from a gene-centric to an isoform-centric research paradigm.

Genome-wide analysis revealed the dysregulation of RNA binding protein-correlated alternative splicing events in myocardial ischemia reperfusion injury

BACKGROUND

Myocardial ischemia reperfusion injury (MIRI), the tissue damage which is caused by the returning of blood supply to tissue after a period of ischemia, greatly reduces the therapeutic effect of treatment of myocardial infarction. But the underlying functional mechanisms of MIRI are still unclear.

METHODS

We constructed mouse models of MIRI, extracted injured and healthy myocardial tissues, and performed transcriptome sequencing experiments (RNA-seq) to systematically investigate the dysregulated transcriptome of MIRI, especially the alternative splicing (AS) regulation and RNA binding proteins (RBPs). Selected RBPs and MIRI-associated AS events were then validated by RT-qPCR experiments.

RESULTS

The differentially expressed gene (DEG) analyses indicated that transcriptome profiles were changed by MIRI and that DEGs' enriched functions were consistent with MIRI's dysregulated pathways. Furthermore, the AS profile was synergistically regulated and showed clear differences between the mouse model and the healthy samples. The exon skipping events significantly increased in MIRI model samples, while the opposite cassette exon events significantly decreased. According to the functional analysis, regulated alternative splicing genes (RASGs) were enriched in protein transport, cell division /cell cycle, RNA splicing, and endocytosis pathways, which were associated with the development of MIRI. Meanwhile, 493 differentially expressed RBPs (DE RBPs) were detected, most of which were correlated with the changed ratios of AS events. In addition, nine DE RBP genes were validated, including Eif5, Pdia6, Tagln2, Vasp, Zfp36l2, Grsf1, Idh2, Ndrg2, and Uqcrc1. These nine DE RBPs were correlated with RASGs enriched in translation process, cell growth and division, and endocytosis pathways, highly consistent with the functions of all RASGs. Finally, we validated the AS ratio changes of five regulated alternative splicing events (RASEs) derived from important regulatory genes, including Mtmr3, Cdc42, Cd47, Fbln2, Vegfa, and Fhl2.

CONCLUSION

Our study emphasized the critical roles of the dysregulated AS profiles in MIRI development, investigated the potential functions of MIRI-associated RASGs, and identified regulatory RBPs involved in AS regulation. We propose that the identified RASEs and RBPs could serve as important regulators and potential therapeutic targets in MIRI treatment in the future.

Cysteine Counting via Isotopic Chemical Labeling for Intact Mass Proteoform Identifications in Tissue

Top-down proteomics, the tandem mass spectrometric analysis of intact proteoforms, is the dominant method for proteoform characterization in complex mixtures. While this strategy produces detailed molecular information, it also requires extensive instrument time per mass spectrum obtained and thus compromises the depth of proteoform coverage that is accessible on liquid chromatography time scales. Such a top-down analysis is necessary for making original proteoform identifications, but once a proteoform has been confidently identified, the extensive characterization it provides may no longer be required for a subsequent identification of the same proteoform. We present a strategy to identify proteoforms in tissue samples on the basis of the combination of an intact mass determination with a measured count of the number of cysteine residues present in each proteoform. We developed and characterized a cysteine tagging chemistry suitable for the efficient and specific labeling of cysteine residues within intact proteoforms and for providing a count of the cysteine amino acids present. On simple protein mixtures, the tagging chemistry yields greater than 98% labeling of all cysteine residues, with a labeling specificity of greater than 95%. Similar results are observed on more complex samples. In a proof-of-principle study, proteoforms present in a human prostate tumor biopsy were characterized. Observed proteoforms, each characterized by an intact mass and a cysteine count, were grouped into proteoform families (groups of proteoforms originating from the same gene). We observed 2190 unique experimental proteoforms, 703 of which were grouped into 275 proteoform families.

Isoform switching leads to downregulation of cytokine producing genes in estrogen receptor positive breast cancer

Objective: Estrogen receptor breast cancer (BC) is characterized by the expression of estrogen receptors. It is the most common cancer among women, with an incidence rate of 2.26 million cases worldwide. The aim of this study was to identify differentially expressed genes and isoform switching between estrogen receptor positive and triple negative BC samples. Methods: The data were collected from ArrayExpress, followed by preprocessing and subsequent mapping from HISAT2. Read quantification was performed by StringTie, and then R package ballgown was used to perform differential expression analysis. Functional enrichment analysis was conducted using Enrichr, and then immune genes were shortlisted based on the ScType marker database. Isoform switch analysis was also performed using the IsoformSwitchAnalyzeR package. Results: A total of 9,771 differentially expressed genes were identified, of which 86 were upregulated and 117 were downregulated. Six genes were identified as mainly associated with estrogen receptor positive BC, while a novel set of ten genes were found which have not previously been reported in estrogen receptor positive BC. Furthermore, alternative splicing and subsequent isoform usage in the immune system related genes were determined. Conclusion: This study identified the differential usage of isoforms in the immune system related genes in cancer cells that suggest immunosuppression due to the dysregulation of CXCR chemokine receptor binding, iron ion binding, and cytokine activity.

Tissue Usage Preference and Intrinsically Disordered Region Remodeling of Alternative Splicing Derived Proteoforms in the Heart

A computational analysis of mass spectrometry data was performed to uncover alternative splicing derived protein variants across chambers of the human heart. Evidence for 216 non-canonical isoforms was apparent in the atrium and the ventricle, including 52 isoforms not documented on SwissProt and recovered using an RNA sequencing derived database. Among non-canonical isoforms, 29 show signs of regulation based on statistically significant preferences in tissue usage, including a ventricular enriched protein isoform of tensin-1 (TNS1) and an atrium-enriched PDZ and LIM Domain 3 (PDLIM3) isoform 2 (PDLIM3-2/ALP-H). Examined variant regions that differ between alternative and canonical isoforms are highly enriched in intrinsically disordered regions, and over two-thirds of such regions are predicted to function in protein binding and/or RNA binding. The analysis here lends further credence to the notion that alternative splicing diversifies the proteome by rewiring intrinsically disordered regions, which are increasingly recognized to play important roles in the generation of biological function from protein sequences.

Top-Down Protein Analysis by Tandem-Trapped Ion Mobility Spectrometry/Mass Spectrometry (Tandem-TIMS/MS) Coupled with Ultraviolet Photodissociation (UVPD) and Parallel Accumulation/Serial Fragmentation (PASEF) MS/MS Analysis

"Top-down" proteomics analyzes intact proteins and identifies proteoforms by their intact mass as well as the observed fragmentation pattern in tandem mass spectrometry (MS/MS) experiments. Recently, hybrid ion mobility spectrometry-mass spectrometry (IM/MS) methods have gained traction for top-down experiments, either by allowing top-down analysis of individual isomers or alternatively by improving signal/noise and dynamic range for fragment ion assignment. We recently described the construction of a tandem-trapped ion mobility spectrometer/mass spectrometer (tandem-TIMS/MS) coupled with an ultraviolet (UV) laser and demonstrated a proof-of-principle for top-down analysis by UV photodissociation (UVPD) at 2-3 mbar. The present work builds on this with an exploration of a top-down method that couples tandem-TIMS/MS with UVPD and parallel-accumulation serial fragmentation (PASEF) MS/MS analysis. We first survey types and structures of UVPD-specific fragment ions generated in the 2-3 mbar pressure regime of our instrument. Notably, we observe UVPD-induced fragment ions with multiple conformations that differ from those produced in the absence of UV irradiation. Subsequently, we discuss how MS/MS spectra of top-down fragment ions lend themselves ideally for probability-based scoring methods developed in the bottom-up proteomics field and how the ability to record automated PASEF-MS/MS spectra resolves ambiguities in the assignment of top-down fragment ions. Finally, we describe the coupling of tandem-TIMS/MS workflows with UVPD and PASEF-MS/MS analysis for native top-down protein analysis.

Extending protein interaction networks using proteoforms and small molecules

MOTIVATION

Biological network analysis for high-throughput biomedical data interpretation relies heavily on topological characteristics. Networks are commonly composed of nodes representing genes or proteins that are connected by edges when interacting. In this study, we use the rich information available in the Reactome pathway database to build biological networks accounting for small molecules and proteoforms modeled using protein isoforms and post-translational modifications to study the topological changes induced by this refinement of the network representation.

RESULTS

We find that improving the interactome modeling increases the number of nodes and interactions, but that isoform and post-translational modification annotation is still limited compared to what can be expected biologically. We also note that small molecule information can distort the topology of the network due to the high connectedness of these molecules, which does not necessarily represent the reality of biology. However, by restricting the connections of small molecules to the context of biochemical reactions, we find that these improve the overall connectedness of the network and reduce the prevalence of isolated components and nodes. Overall, changing the representation of the network alters the prevalence of articulation points and bridges globally but also within and across pathways. Hence, some molecules can gain or lose in biological importance depending on the level of detail of the representation of the biological system, which might in turn impact network-based studies of diseases or druggability.

AVAILABILITY AND IMPLEMENTATION

Networks are constructed based on data publicly available in the Reactome Pathway knowledgebase: reactome.org.

SNORA73B promotes endometrial cancer progression through targeting MIB1 and regulating host gene RCC1 alternative splicing

Endometrial cancer (EC) is a common gynaecological malignant tumour with unclear pathogenesis. Small nucleolar RNA (snoRNA) is involved in many biological processes, including those of cancers. Using the Cancer Genome Atlas (TCGA) database, the expression pattern of a snoRNA, SNORA73B, was analysed. The biological functions of SNORA73B were assessed by in vitro proliferation, apoptosis, migration, and invasion assays and in vivo by the xenograft model. RNA sequencing (RNA-seq) and RNA immunoprecipitation assays were performed to determine the relationship between SNORA73B and its target genes. High-performance liquid chromatography (HPLC) was performed to detect the pseudouridine content of the mindbomb E3 ubiquitin protein ligase 1 gene (MIB1). The stability of MIB1 mRNA was evaluated using a transcription inhibitor, actinomycin D. By performing co-immunoprecipitation assays, the change in the ubiquitin levels of the Jagged canonical Notch ligand 1 (Jag 1), caused by SNORA73B and MIB1, was identified. RNA-seq and qRT-PCR were performed to detect the alternative splicing of the regulator of the chromosome condensation 1 gene (RCC1). The TCGA database analysis showed that SNORA73B was highly expressed in EC. SNORA73B promoted cell proliferation, migration, and invasion and inhibited apoptosis. SNORA73B modified the pseudouridine content in MIB1 and increased the stability of MIB1 mRNA and protein; thus, it affected Jag 1 ubiquitination and further activated the Notch pathway. SNORA73B also affected the alternative splicing of RCC1, increasing the number of transcripts, RCC1-T2 and RCC1-T3, which promoted cell proliferation, migration, and invasion. SNORA73B can be a potential target for EC.

Long-read sequencing reveals the landscape of aberrant alternative splicing and novel therapeutic target in colorectal cancer

BACKGROUND

Alternative splicing complexity plays a vital role in carcinogenesis and cancer progression. Improved understanding of novel splicing events and the underlying regulatory mechanisms may contribute new insights into developing new therapeutic strategies for colorectal cancer (CRC).

METHODS

Here, we combined long-read sequencing technology with short-read RNA-seq methods to investigate the transcriptome complexity in CRC. By using experiment assays, we explored the function of newly identified splicing isoform TIMP1 Δ4-5. Moreover, a CRISPR/dCasRx-based strategy to induce the TIMP1 exon 4-5 exclusion was introduced to inhibit neoplasm growth.

RESULTS

A total of 90,703 transcripts were identified, of which > 62% were novel compared with current transcriptome annotations. These novel transcripts were more likely to be sample specific, expressed at relatively lower levels with more exons, and oncogenes displayed a characteristic to generate more transcripts in CRC. Clinical outcome data analysis showed that 1472 differentially expressed alternative splicing events (DEAS) were tightly associated with CRC patients' prognosis, and many novel isoforms were likely to be important determinants for patient survival. Among these, newly identified splicing isoform TIMP1 Δ4-5 was significantly downregulated in CRC. Further in vitro and in vivo assays demonstrated that ectopic expression of TIMP1 Δ4-5 significantly suppresses tumor cell growth and metastasis. Serine/arginine-rich splicing factor 1 (SRSF1) acts as a onco-splicing regulator through sustaining the inclusion of TIMP1 exon 4-5. Furthermore, CRISPR/dCasRx-based strategies designed to induce TIMP1 exon 4-5 exclusion have the potential to restrain the CRC growth.

CONCLUSIONS

This data provides a rich resource for deeper studies of gastrointestinal malignancies. Newly identified splicing isoform TIMP1 Δ4-5 plays an important role in mediating CRC progression and may be a potential therapy target in CRC.

Building alternative splicing and evolution-aware sequence-structure maps for protein repeats

Alternative splicing of repeats in proteins provides a mechanism for rewiring and fine-tuning protein interaction networks. In this work, we developed a robust and versatile method, ASPRING, to identify alternatively spliced protein repeats from gene annotations. ASPRING leverages evolutionary meaningful alternative splicing-aware hierarchical graphs to provide maps between protein repeats sequences and 3D structures. We re-think the definition of repeats by explicitly accounting for transcript diversity across several genes/species. Using a stringent sequence-based similarity criterion, we detected over 5,000 evolutionary conserved repeats by screening virtually all human protein-coding genes and their orthologs across a dozen species. Through a joint analysis of their sequences and structures, we extracted specificity-determining sequence signatures and assessed their implication in experimentally resolved and modelled protein interactions. Our findings demonstrate the widespread alternative usage of protein repeats in modulating protein interactions and open avenues for targeting repeat-mediated interactions.

Functional insight into a neurodevelopmental disorder caused by missense variants in an RNA-binding protein, RBM10

The RNA-binding motif protein 10, RBM10, is an RNA splicing regulator essential for development. Loss-of-function RBM10 variants are associated with TARP syndrome, a severe X-linked recessive condition in males. We report a 3-year-old male with a mild phenotype, consisting of cleft palate, hypotonia, developmental delay, and minor dysmorphisms, associated with a missense RBM10 variant, c.943T>C, p.Ser315Pro, affecting the RRM2 RNA-binding domain. His clinical features were similar to a previously reported case associated with a missense variant. The p.Ser315Pro mutant protein was expressed normally in the nucleus, but its expression level and protein stability were slightly reduced. Nuclear magnetic resonance spectroscopy showed that the structure and the RNA-binding ability of the RRM2 domain with the p.Ser315Pro were unaffected. However, it affects the alternative splicing regulations of downstream genes, NUMB and TNRC6A, and its splicing alteration patterns were variable depending on target transcripts. In summary, a novel germline missense RBM10 p.Ser315Pro variant that causes functional changes in the expression of its downstream genes results in a non-lethal phenotype associated with developmental delays. The functional alteration effects depend on the residues affected by missense variants. Our findings are expected to bring broader insights into the RBM10-associated genotype-phenotype relationships by delineating the molecular mechanism of RBM10 functions.

Alternative splicing regulates adaptor protein binding, trafficking, and activity of the Vps10p domain receptor SorCS2 in neuronal development

The Vps10p domain receptor SorCS2 is crucial for the development and function of the nervous system and essential for brain-derived neurotrophic factor (BDNF)-induced changes in neuronal morphology and plasticity. SorCS2 regulates the subcellular trafficking of the BDNF signaling receptor TrkB as well as selected neurotransmitter receptors in a manner that is dependent on the SorCS2 intracellular domain (ICD). However, the cellular machinery and adaptor protein (AP) interactions that regulate receptor trafficking via the SorCS2 ICD are unknown. We here identify four splice variants of human SorCS2 differing in the insertion of an acidic cluster motif and/or a serine residue within the ICD. We show that each variant undergoes posttranslational proteolytic processing into a one- or two-chain receptor, giving rise to eight protein isoforms, the expression of which differs between neuronal and nonneuronal tissues and is affected by cellular stressors. We found that the only variants without the serine were able to rescue BDNF-induced branching of SorCS2 knockout hippocampal neurons, while variants without the acidic cluster showed increased interactions with clathrin-associated APs AP-1, AP-2, and AP-3. Using yeast two-hybrid screens, we further discovered that all variants bound dynein light chain Tctex-type 3; however, only variants with an acidic cluster motif bound kinesin light chain 1. Accordingly, splice variants showed markedly different trafficking properties and localized to different subcellular compartments. Taken together, our findings demonstrate the existence of eight functional SorCS2 isoforms with differential capacity for interactions with cytosolic ligands dynein light chain Tctex-type 3 and kinesin light chain 1, which potentially allows cell-type specific SorCS2 trafficking and BDNF signaling.

Capillary Zone Electrophoresis-Tandem Mass Spectrometry for Top-Down Proteomics of Mouse Brain Integral Membrane Proteins

Mass spectrometry (MS)-based top-down characterization of integral membrane proteins (IMPs) is crucial for understanding their functions in biological processes. However, it is technically challenging due to their low solubility in typical MS-compatible buffers. In this work, for the first time, we developed an efficient capillary zone electrophoresis (CZE)-tandem MS (MS/MS) method for the top-down proteomics (TDP) of IMPs enriched from mouse brains. Our technique employs a sample buffer containing 30% (v/v) formic acid and 60% (v/v) methanol for solubilizing IMPs and utilizes a separation buffer of 30% (v/v) acetic acid and 30% (v/v) methanol for maintaining the solubility of IMPs during CZE separation. Single-shot CZE-MS/MS identified 51 IMP proteoforms from the mouse brain sample. Coupling size exclusion chromatography (SEC) to CZE-MS/MS enabled the identification of 276 IMP proteoforms from the mouse brain sample containing 1-4 transmembrane domains. This proof-of-concept work demonstrates the high potential of CZE-MS/MS for the large-scale TDP of IMPs.

Increased NMDARs in neurons and glutamine synthetase in astrocytes underlying autistic-like behaviors of Gabrb1-/- mice

Mutations of the GABA-A receptor subunit β1 (GABRB1) gene are found in autism patients. However, it remains unclear how mutations in Gabrb1 may lead to autism. We generated Gabrb1-/- mouse model, which showed autistic-like behaviors. We carried out RNA-seq on the hippocampus and found glutamatergic pathway may be involved. We further carried out single-cell RNA sequencing on the whole brain followed by qRT-PCR, immunofluorescence, electrophysiology, and metabolite detection on specific cell types. We identified the up-regulated Glul/Slc38a3 in astrocytes, Grin1/Grin2b in neurons, glutamate, and the ratio of Glu/GABA in the hippocampus. Consistent with these results, increased NMDAR-currents and reduced GABAAR-currents in the CA1 neurons were detected in Gabrb1-/- mice. NMDAR antagonist memantine or Glul inhibitor methionine sulfoximine could rescue the abnormal behaviors in Gabrb1-/- mice. Our data reveal that upregulation of the glutamatergic synapse pathway, including NMDARs at neuronal synapses and glutamine exported by astrocytes, may lead to autistic-like behaviors.

N-terminal proteoforms may engage in different protein complexes

Alternative translation initiation and alternative splicing may give rise to N-terminal proteoforms, proteins that differ at their N-terminus compared with their canonical counterparts. Such proteoforms can have altered localizations, stabilities, and functions. Although proteoforms generated from splice variants can be engaged in different protein complexes, it remained to be studied to what extent this applies to N-terminal proteoforms. To address this, we mapped the interactomes of several pairs of N-terminal proteoforms and their canonical counterparts. First, we generated a catalogue of N-terminal proteoforms found in the HEK293T cellular cytosol from which 22 pairs were selected for interactome profiling. In addition, we provide evidence for the expression of several N-terminal proteoforms, identified in our catalogue, across different human tissues, as well as tissue-specific expression, highlighting their biological relevance. Protein-protein interaction profiling revealed that the overlap of the interactomes for both proteoforms is generally high, showing their functional relation. We also showed that N-terminal proteoforms can be engaged in new interactions and/or lose several interactions compared with their canonical counterparts, thus further expanding the functional diversity of proteomes.

Patterns of transcription factor binding and epigenome at promoters allow interpretable predictability of multiple functions of non-coding and coding genes

Understanding the biological roles of all genes only through experimental methods is challenging. A computational approach with reliable interpretability is needed to infer the function of genes, particularly for non-coding RNAs. We have analyzed genomic features that are present across both coding and non-coding genes like transcription factor (TF) and cofactor ChIP-seq (823), histone modifications ChIP-seq (n = 621), cap analysis gene expression (CAGE) tags (n = 255), and DNase hypersensitivity profiles (n = 255) to predict ontology-based functions of genes. Our approach for gene function prediction was reliable (>90% balanced accuracy) for 486 gene-sets. PubMed abstract mining and CRISPR screens supported the inferred association of genes with biological functions, for which our method had high accuracy. Further analysis revealed that TF-binding patterns at promoters have high predictive strength for multiple functions. TF-binding patterns at the promoter add an unexplored dimension of explainable regulatory aspects of genes and their functions. Therefore, we performed a comprehensive analysis for the functional-specificity of TF-binding patterns at promoters and used them for clustering functions to reveal many latent groups of gene-sets involved in common major cellular processes. We also showed how our approach could be used to infer the functions of non-coding genes using the CRISPR screens of coding genes, which were validated using a long non-coding RNA CRISPR screen. Thus our results demonstrated the generality of our approach by using gene-sets from CRISPR screens. Overall, our approach opens an avenue for predicting the involvement of non-coding genes in various functions.

Proteins with amino acid repeats constitute a rapidly evolvable and human-specific essentialome

Protein products of essential genes, indispensable for organismal survival, are highly conserved and bring about fundamental functions. Interestingly, proteins that contain amino acid homorepeats that tend to evolve rapidly are enriched in eukaryotic essentialomes. Why are proteins with hypermutable homorepeats enriched in conserved and functionally vital essential proteins? We solve this functional versus evolutionary paradox by demonstrating that human essential proteins with homorepeats bring about crosstalk across biological processes through high interactability and have distinct regulatory functions affecting expansive global regulation. Importantly, essential proteins with homorepeats rapidly diverge with the amino acid substitutions frequently affecting functional sites, likely facilitating rapid adaptability. Strikingly, essential proteins with homorepeats influence human-specific embryonic and brain development, implying that the presence of homorepeats could contribute to the emergence of human-specific processes. Thus, we propose that homorepeat-containing essential proteins affecting species-specific traits can be potential intervention targets across pathologies, including cancers and neurological disorders.

Pilot investigation of magnetic nanoparticle-based immobilized metal affinity chromatography for efficient enrichment of phosphoproteoforms for mass spectrometry-based top-down proteomics

Protein phosphorylation is a vital and common post-translational modification (PTM) in cells, modulating various biological processes and diseases. Comprehensive top-down proteomics of phosphorylated proteoforms (phosphoproteoforms) in cells and tissues is essential for a better understanding of the roles of protein phosphorylation in fundamental biological processes and diseases. Mass spectrometry (MS)-based top-down proteomics of phosphoproteoforms remains challenging due to their relatively low abundance. Herein, we investigated magnetic nanoparticle-based immobilized metal affinity chromatography (IMAC, Ti4+, and Fe3+) for selective enrichment of phosphoproteoforms for MS-based top-down proteomics. The IMAC method achieved reproducible and highly efficient enrichment of phosphoproteoforms from simple and complex protein mixtures. It outperformed one commercial phosphoprotein enrichment kit regarding the capture efficiency and recovery of phosphoproteins. Reversed-phase liquid chromatography (RPLC)-tandem mass spectrometry (MS/MS) analyses of yeast cell lysates after IMAC (Ti4+ or Fe3+) enrichment produced roughly 100% more phosphoproteoform identifications compared to without IMAC enrichment. Importantly, the phosphoproteoforms identified after Ti4+-IMAC or Fe3+-IMAC enrichment correspond to proteins with much lower overall abundance compared to that identified without the IMAC treatment. We also revealed that Ti4+-IMAC and Fe3+-IMAC could enrich different pools of phosphoproteoforms from complex proteomes and the combination of those two methods will be useful for further improving the phosphoproteoform coverage from complex samples. The results clearly demonstrate the value of our magnetic nanoparticle-based Ti4+-IMAC and Fe3+-IMAC for advancing top-down MS characterization of phosphoproteoforms in complex biological systems.

Quantitative Analysis of Isoform Switching in Cancer

Over the past 8 years, multiple studies examined the phenomenon of isoform switching in human cancers and discovered that isoform switching is widespread, with hundreds to thousands of such events per cancer type. Although all of these studies used slightly different definitions of isoform switching, which in part led to a rather poor overlap of their results, they all leveraged transcript usage, a proportion of the transcript's expression in the total expression level of the parent gene, to detect isoform switching. However, how changes in transcript usage correlate with changes in transcript expression is not sufficiently explored. In this article, we adopt the most common definition of isoform switching and use a state-of-the-art tool for the analysis of differential transcript usage, SatuRn, to detect isoform switching events in 12 cancer types. We analyze the detected events in terms of changes in transcript usage and the relationship between transcript usage and transcript expression on a global scale. The results of our analysis suggest that the relationship between changes in transcript usage and changes in transcript expression is far from straightforward, and that such quantitative information can be effectively used for prioritizing isoform switching events for downstream analyses.

Transcription factors and splice factors - interconnected regulators of stem cell differentiation

Purpose of review

The underlying molecular mechanisms that direct stem cell differentiation into fully functional, mature cells remain an area of ongoing investigation. Cell state is the product of the combinatorial effect of individual factors operating within a coordinated regulatory network. Here, we discuss the contribution of both gene regulatory and splicing regulatory networks in defining stem cell fate during differentiation and the critical role of protein isoforms in this process.

Recent findings

We review recent experimental and computational approaches that characterize gene regulatory networks, splice regulatory networks, and the resulting transcriptome and proteome they mediate during differentiation. Such approaches include long-read RNA sequencing, which has demonstrated high-resolution profiling of mRNA isoforms, and Cas13-based CRISPR, which could make possible high-throughput isoform screening. Collectively, these developments enable systems-level profiling of factors contributing to cell state.

Summary

Overall, gene and splice regulatory networks are important in defining cell state. The emerging high-throughput systems-level approaches will characterize the gene regulatory network components necessary in driving stem cell differentiation.

Auxin-inducible degron system: an efficient protein degradation tool to study protein function

Targeted protein degradation, with its rapid protein depletion kinetics, allows the measurement of acute changes in the cell. The auxin-inducible degron (AID) system, rapidly degrades AID-tagged proteins only in the presence of auxin. The AID system being inducible makes the study of essential genes and dynamic processes like cell differentiation, cell cycle and genome organization feasible. The AID degradation system has been adapted to yeast, protozoans, C. elegans, Drosophila, zebrafish, mouse and mammalian cell lines. Using the AID system, researchers have unveiled novel functions for essential proteins at developmental stages that were previously difficult to investigate due to early lethality. This comprehensive review discusses the development, advancements, applications and drawbacks of the AID system and compares it with other available protein degradation systems.

SRSF5 Regulates the Expression of BQ323636.1 to Modulate Tamoxifen Resistance in ER-Positive Breast Cancer

About 70% of breast cancer patients are oestrogen receptor-positive (ER +ve). Adjuvant endocrine therapy using tamoxifen (TAM) is an effective approach for preventing local recurrence and metastasis. However, around half of the patients will eventually develop resistance. Overexpression of BQ323636.1 (BQ) is one of the mechanisms that confer TAM resistance. BQ is an alternative splice variant of NCOR2. The inclusion of exon 11 generates mRNA for NCOR2, while the exclusion of exon 11 produces mRNA for BQ. The expression of SRSF5 is low in TAM-resistant breast cancer cells. Modulation of SRSF5 can affect the alternative splicing of NCOR2 to produce BQ. In vitro and in vivo studies confirmed that the knockdown of SRSF5 enhanced BQ expression, and conferred TAM resistance; in contrast, SRSF5 overexpression reduced BQ expression and, thus, reversed TAM resistance. Clinical investigation using a tissue microarray confirmed the inverse correlation of SRSF5 and BQ. Low SRSF5 expression was associated with TAM resistance, local recurrence and metastasis. Survival analyses showed that low SRSF5 expression was associated with poorer prognosis. We showed that SRPK1 can interact with SRSF5 to phosphorylate it. Inhibition of SRPK1 by a small inhibitor, SRPKIN-1, suppressed the phosphorylation of SRSF5. This enhanced the proportion of SRSF5 interacting with exon 11 of NCOR2, reducing the production of BQ mRNA. As expected, SRPKIN-1 reduced TAM resistance. Our study confirms that SRSF5 is essential for BQ expression. Modulating the activity of SRSF5 in ER +ve breast cancer will be a potential approach to combating TAM resistance.

To Explore the Binding Affinity of Human γ-Secretase Activating Protein (GSAP) Isoform 4 with APP-C99 Peptides

γ-Secretase activating protein (GSAP) is known to play an important role in the β-amyloid pathway. It acts as a modulator and accentuates the truncation of the amyloid precursor protein C-99 fragment through the γ-secretase complex. GSAP has four isoforms, out of which canonical isoform 1, a 16 kDa C-terminal portion, has been extensively studied, whereas the function of other three isoforms remains unknown. Here, we explore the GSAP isoform 4 (GSAP_I4) expression and purification from inclusion bodies followed by the refolding of the protein. The secondary structure of GSAP_I4 is predicted using circular dichroism. The protein is further characterized by western blotting and mass spectroscopy analysis. Additionally, biochemical assays and in silico molecular docking and molecular simulation are performed to investigate the binding of GSAP_I4 and APP-C99 peptide fragments. The results reflect that although GSAP_I1 and GSAP_I4 share high sequence similarity, the isoform 4 does not show any affinity toward APP-C99 peptide fragments. This hints toward the fact that GSAP_I4 might have a different role in the living system that is yet unexplored.

Evaluation of the clinical significance of global mRNA alternative splicing in patients with acute myeloid leukemia

Aberrant alternative splicing (AS) is involved in leukemogenesis. This study explored the clinical impact of alterations in global AS patterns in 341 patients with acute myeloid leukemia (AML) newly diagnosed at the National Taiwan University Hospital and validated it using The Cancer Genome Atlas (TCGA) cohort. While studying normal cord blood CD34⁺ /CD38^- cells, we found that AML cells exhibited significantly different global splicing patterns. AML with mutated TP53 had a particularly high degree of genome-wide aberrations in the splicing patterns. Aberrance in the global splicing pattern was an independent unfavorable prognostic factor affecting the overall survival of patients with AML receiving standard intensive chemotherapy. The integration of global splicing patterns into the 2022 European LeukemiaNet risk classification could stratify AML patients into four groups with distinct prognoses in both our experimental and TCGA cohorts. We further identified four genes with AS alterations that harbored prognostic significance in both of these cohorts. Moreover, these survival-associated AS events are involved in several important cellular processes that might be associated with poor response to intensive chemotherapy. In summary, our study demonstrated the clinical and biological implications of differential global splicing patterns in AML patients. Further studies with larger prospective cohorts are required to confirm these findings.

Comparison of alternative splicing (AS) events in adipose tissue of polled dorset versus small tail han sheep

Background

During the alternative splicing (AS), the exons of primary transcripts are spliced in various arrangements, resulting in structurally and functionally distinct mRNAs and proteins. This study aimed to examine genes with AS events from Small Tail Han sheep and Dorset sheep to explore the mechanism of adipose developments.

Methods

This study identified the genes with AS events in adipose tissues of two different sheep with next-generation sequencing. In this paper, genes with significantly different AS events were performed gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses.

Results

364 genes with 411 A S events showed significant differences in adipose tissues between the two breeds; 108 genes with 120 A S events were extremely significant differences between the two breeds. We identified several novel genes that are related with adipose growth and development. The results of KEGG and GO analysis indicated that oocyte meiosis, mitogen-activated protein kinase (Wnt), mitogen-activated protein kinase (MAPK) signaling pathway, etc. Were closely related to the adipose tissue developments.

Conclusions

This paper revealed that the genes with AS events are important for adipose tissues in sheep, exploring the mechanisms of AS events associated with adipose tissue developments in sheep of different breeds.

A fluorescent probe to simultaneously detect both O-GlcNAcase and phosphatase

O-GlcNAc modification of proteins often has crosstalk with protein phosphorylation. These posttranslational modifications are highly dynamic events that modulate a wide range of cellular processes. Owing to the physiological and pathological significance of protein O-GlcNAcylation and phosphorylation, we designed the fluorescent probe, βGlcNAc-CM-Rhod-P, to differentially detect activities of O-GlcNAcase (OGA) and phosphatase, enzymes that are responsible for these modifications. βGlcNAc-CM-Rhod-P was comprised of a βGlcNAc-conjugated coumarin (βGlcNAc-CM) acting as an OGA substrate, a phosphorylated rhodol (Rhod-P) as a phosphatase substrate and a piperazine bridge. Because the emission wavelength maxima of CM and Rhod liberated from the probe are greatly different (100 nm), spectral interference is avoided. The results of this study revealed that treatment of βGlcNAc-CM-Rhod-P with OGA promotes formation of the GlcNAc-cleaved probe, CM-Rhod-P, and a consequent increase in the intensity of fluorescence associated with free CM. Also, it was found that exposure of the probe to phosphatase produces a dephosphorylated probe, βGlcNAc-CM-Rhod, which displays strong fluorescence arising from free Rhod. On the other hand, when incubated with both OGA and phosphatase, βGlcNAc-CM-Rhod-P was converted to CM-Rhod which lacked both βGlcNAc and phosphoryl groups, in conjunction with increases in the intensities of fluorescence arising from both free CM and Rhod. This probe was employed to detect activities of OGA and phosphatase in cell lysates and to fluorescently image both enzymes in cells. Collectively, the findings indicate that βGlcNAc-CM-Rhod-P can be utilized as a chemical tool to simultaneously determine activities of OGA and phosphatase.

Malignant DFFB isoform switching promotes leukemia survival in relapse pediatric T-cell acute lymphoblastic leukemia

With modern treatment most children with acute lymphoblastic leukemia (ALL) survive without relapse. However, for children who relapse the prognosis is still poor, especially in children with T-cell phenotype (T-ALL) and remains the major cause of death. The exact mechanism of relapse is currently not known. While contribution of RNA processing alteration has been linked to other hematological malignancies, its contribution in pediatric T-ALL may provide new insights. Almost all human genes express more than one alternative splice isoform. Thus, gene modulation producing a diverse repertoire of the transcriptome and proteome have become a significant molecular marker of cancer and a potential therapeutic vulnerability. To study this, we performed RNA-sequencing analysis on patient-derived samples followed by splice isoform-specific PCR. We uncovered a distinct RNA splice isoform expression pattern characteristic for relapse samples compared to the leukemia samples from the time of diagnosis. We also identified deregulated splicing and apoptosis pathways specific for relapse T-ALL. Moreover, patients with T-ALL displayed pro-survival splice isoform switching favoring pro-survival isoforms compared to normal healthy stem cells. Cumulatively, pro-survival isoform switching and DFFB isoform regulation of SOX2 and MYCN may play a role in T-ALL proliferation and survival, thus serving as a potential therapeutic option.

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.

Perspective on the potential of tandem-ion mobility /mass spectrometry methods for structural proteomics applications

Cellular processes are usually carried out collectively by the entirety of all proteins present in a biological cell, i.e. the proteome. Mass spectrometry-based methods have proven particularly successful in identifying and quantifying the constituent proteins of proteomes, including different molecular forms of a protein. Nevertheless, protein sequences alone do not reveal the function or dysfunction of the identified proteins. A straightforward way to assign function or dysfunction to proteins is characterization of their structures and dynamics. However, a method capable to characterize detailed structures of proteins and protein complexes in a large-scale, systematic manner within the context of cellular processes does not yet exist. Here, we discuss the potential of tandem-ion mobility / mass spectrometry (tandem-IM/MS) methods to provide such ability. We highlight the capability of these methods using two case studies on the protein systems ubiquitin and avidin using the tandem-TIMS/MS technology developed in our laboratory and discuss these results in the context of other developments in the broader field of tandem-IM/MS.

Identification of EMT-related alternative splicing event of TMC7 to promote invasion and migration of pancreatic cancer

Objective

Epithelial-to-mesenchymal transition (EMT) is tightly associated with the invasion and metastasis of pancreatic cancer with rapid progression and poor prognosis. Notably, gene alternative splicing (AS) event plays a critical role in regulating the progression of pancreatic cancer. Therefore, this study aims to identify the EMT-related AS event in pancreatic cancer.

Methods

The EMT-related gene sets, transcriptomes, and matched clinical data were obtained from the MSigDB, The Cancer Genome Atlas (TCGA), International Cancer Genome Consortium (ICGC), and Gene Expression Omnibus (GEO) databases. Key gene AS events associated with liver metastasis were identified by prognostic analysis, gene set variation analysis (GSVA), and correlation analysis in pancreatic cancer. The cell line and organoid model was constructed to evaluate these key gene AS events in regulating pancreatic cancer in vitro. Furthermore, we established an EMT-related gene set consisting of 13 genes by prognostic analysis, the role of which was validated in two other databases. Finally, the human pancreatic cancer tissue and organoid model was used to evaluate the correlation between the enrichment of this gene set and liver metastasis.

Results

Prognostic analysis and correlation analysis revealed that eight AS events were closely associated with the prognosis of pancreatic cancer. Furthermore, the expression of TMC7 and CHECK1 AS events was increased in the metastatic lesions of the human tissue and organoid model. Additionally, the knockdown of exon 17 of TMC7 significantly inhibited the proliferation, invasion, and migration of pancreatic cancer cells in 2D and 3D cell experiments. Finally, the expression of exon 17 of TMC17 exhibited a significant correlation with the poor prognosis in pancreatic ductal adenocarcinoma (PDAC).

Conclusion

The AS events of TMC7 and CHECK1 were associated with liver metastasis in pancreatic cancer. Moreover, exon 17 of TMC7 could be a potential therapeutic target in pancreatic cancer.

Identifying transcript-level differential expression in primary human immune cells

BACKGROUND

Multipotential hematopoietic stem cells differentiate into a wide variety of immune cells with a diversity of functions, including the ability to respond to a variety of stimuli. Importantly, numerous studies have demonstrated the importance of gene transcription in defining cell identity and functions. While these studies have primarily been performed at the level of the gene, it is known that key immune genes such as CD44 and CD45 generate multiple different transcripts that are differentially expressed across different immune cells, and that encode proteins with different sequences and functions. Prior genomic surveys have shown that the mechanisms for generating diversity in expressed transcripts (alternate splicing, alternate transcription start sites, etc.) are very active in immune cells, but have been lacking in terms of identifying genes with multiple transcripts, that are differentially expressed, and likely to affect cell functions.

METHODS

We first identified the set of genes that had at least two transcripts expressed in our RNA sequencing dataset generated from purified populations of neutrophils, monocytes and five lymphocyte populations (B, NK, γδ T, CD4 + T and CD8 + T) from twelve healthy donors. Next, we developed a heuristic approach to identify genes where two or more transcripts have distinct expression patterns across lymphoid and/or myeloid populations. We then focused our annotation and interpretation on differentially expressed transcripts that affect the coding sequence. This process was repeated to identify transcripts that were differentially expressed between monocytes and populations of macrophages and LPS-stimulated macrophages derived from these monocytes in vitro.

RESULTS

We found that over 55 % of genes had two or more expressed transcripts, with an average ∼3 transcripts per gene, and that 70 % of these had at least two of the transcripts that encoded proteins with different sequences. As expected, we identified a complex pattern of differential expression for multiple transcripts encoding the CD45 transmembrane protein, but we also found similar evidence for ten other genes (CD300A, FYB1, GPI, LITAF, PSMA1, PTMA, RPL32, SEPTIN9, SH3BP2, SH3KBP1) when comparing the expression patterns of transcripts within myeloid and lymphoid cells. We also identified five genes with differentially expressed transcripts associated with the transition from monocytes to macrophages (FNBP1, KLF6, and SEPTIN9) or between macrophages and LPS-stimulated macrophages (CD44, OAZ2, and SEPTIN9). For the most part, we found that the different transcripts of these genes are expected to impact specific biological functions, for example the different transcripts of SEPTIN9 likely regulate the cytoskeleton in immune cells via their interactions with actins filaments and microtubules.

CONCLUSIONS

This analytic approach successfully identified multi-transcript genes that are differentially expressed across immune cells and could be applied to other transcriptomic data.

DATA AVAILABILITY STATEMENT

Researchers can request access to the individual-level data from the current study by contacting the Montreal Heart Institute ethics committee at the following institutional email address: cer.icm@icm-mhi.org.

The Role of Alternative Splicing Factors, DDB2-Related Ageing and DNA Damage Repair in the Progression and Prognosis of Stomach Adenocarcinoma Patients

DNA damage response is a key signal transduction pathway in triggering ageing and tumor progression. Abnormal alternative splicing (AS) is associated with tumors and ageing. However, the role of AS factors associated with DNA damage repair and ageing in stomach adenocarcinoma (STAD) remains unclear. We downloaded the percentage of splicing (PSI) values for AS in STAD from the TCGA SpliceSeq database. The PSI values of DNA repair gene AS events were integrated with STAD patient survival data for Cox regression analysis. The prediction model for the overall survival (OS) was constructed by the clinical traits. The tumor immune microenvironment was analyzed by CIBERSORT and ESTIMATE. We detected 824 AS events originating from 166 DNA repair genes. Cox regression analysis provided 21 prognostic AS events connected with OS statistically, and a prognostic prediction model was constructed. The expression of these AS factors was higher in STAD tumors. DDB2 high senescence levels were associated with active immune responses and better survival in STAD patients. We built a novel prognostic model founded on DNA repair genes with AS events and identified that DDB2 may be a potential biomarker to apply in clinics.

Thermal proteome profiling: Insights into protein modifications, associations, and functions

Tracking proteins' biophysical characteristics on a proteome-wide scale can provide valuable information on their functions and interactions. Thermal proteome profiling (TPP) is a multiplexed quantitative proteomics approach that measures changes in protein thermal stability-a key biophysical property-across different cellular states. Developed in 2014, as a target-deconvolution assay for drugs and other small molecules, TPP has since evolved to a system-level biochemical omics technique providing insights into context-dependent changes in protein states. In this review, we summarise key advances in the experimental and data analysis pipeline that have aided this transformation and discuss the recent developments and applications of TPP.

Functional Impact and Regulation of Alternative Splicing in Mouse Heart Development and Disease

Alternative splicing (AS) plays a major role in the generation of transcript diversity. In the heart, roles have been described for some AS variants, but the global impact and regulation of AS patterns are poorly understood. Here, we studied the AS profiles in heart disease, their relationship with heart development, and the regulatory mechanisms controlling AS dynamics in the mouse heart. We found that AS profiles characterized the different groups and that AS and gene expression changes affected independent genes and biological functions. Moreover, AS changes, specifically in heart disease, were associated with potential protein-protein interaction changes. While developmental transitions were mainly driven by the upregulation of MBNL1, AS changes in disease were driven by a complex regulatory network, where PTBP1 played a central role. Indeed, PTBP1 over-expression was sufficient to induce cardiac hypertrophy and diastolic dysfunction, potentially by perturbing AS patterns.

Liver RBFOX2 regulates cholesterol homeostasis via Scarb1 alternative splicing in mice

RNA alternative splicing (AS) expands the regulatory potential of eukaryotic genomes. The mechanisms regulating liver-specific AS profiles and their contribution to liver function are poorly understood. Here, we identify a key role for the splicing factor RNA-binding Fox protein 2 (RBFOX2) in maintaining cholesterol homeostasis in a lipogenic environment in the liver. Using enhanced individual-nucleotide-resolution ultra-violet cross-linking and immunoprecipitation, we identify physiologically relevant targets of RBFOX2 in mouse liver, including the scavenger receptor class B type I (Scarb1). RBFOX2 function is decreased in the liver in diet-induced obesity, causing a Scarb1 isoform switch and alteration of hepatocyte lipid homeostasis. Our findings demonstrate that specific AS programmes actively maintain liver physiology, and underlie the lipotoxic effects of obesogenic diets when dysregulated. Splice-switching oligonucleotides targeting this network alleviate obesity-induced inflammation in the liver and promote an anti-atherogenic lipoprotein profile in the blood, underscoring the potential of isoform-specific RNA therapeutics for treating metabolism-associated diseases.

Where protein structure and cell diversity meet

Protein-protein interaction networks - interactomes - are charted with the hope to understand how phenotypes emerge and how they are altered in disease states. Early efforts to map interactomes have focused on the assembly of context agnostic, reference networks. However, recent studies have mapped interactomes across different cell lines and tissues, finding highly variable interactomes due to the rewiring of protein-protein interactions in different contexts. Increasing evidence points to significant links between protein structure and interactome diversity seen across cell types and tissues. We discuss how recent findings support the key role of alternative splicing and phosphorylation, two well-established regulators of protein structural and functional diversity, in defining cell type- and tissue-specific interactomes. Moreover, we show that intrinsically disordered protein regions are most favorably equipped to support interactome rewiring by acting as hubs of protein structure and function regulation.