Widespread Expansion of Protein Interaction Capabilities by Alternative Splicing

Adhesion GPCR Latrophilin 3 regulates synaptic function of cone photoreceptors in a trans-synaptic manner

Cone photoreceptors mediate daylight vision in vertebrates. Changes in neurotransmitter release at cone synapses encode visual information and is subject to precise control by negative feedback from enigmatic horizontal cells. However, the mechanisms that orchestrate this modulation are poorly understood due to a virtually unknown landscape of molecular players. Here, we report a molecular player operating selectively at cone synapses that modulates effects of horizontal cells on synaptic release. Using an unbiased proteomic screen, we identified an adhesion GPCR Latrophilin3 (LPHN3) in horizontal cell dendrites that engages in transsynaptic control of cones. We detected and characterized a prominent splice isoform of LPHN3 that excludes a element with inhibitory influence on transsynaptic interactions. A gain-of-function mouse model specifically routing LPHN3 splicing to this isoform but not knockout of LPHN3 diminished Ca_V1.4 calcium channel activity profoundly disrupted synaptic release by cones and resulted in synaptic transmission deficits. These findings offer molecular insight into horizontal cell modulation on cone synaptic function and more broadly demonstrate the importance of alternative splicing in adhesion GPCRs for their physiological function.

Alternative Splicing-Based Differences Between Hepatocellular Carcinoma and Intrahepatic Cholangiocarcinoma: Genes, Immune Microenvironment, and Survival Prognosis

Alternative splicing (AS) event is a novel biomarker of tumor tumorigenesis and progression. However, the comprehensive analysis of hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (ICC) is lacking. Differentially expressed analysis was used to identify the differentially expressed alternative splicing (DEAS) events between HCC or ICC tissues and their normal tissues. The correlation between DEAS events and functional analyses or immune features was evaluated. The cluster analysis based on DEAS can accurately reflect the differences in the immune microenvironment between HCC and ICC. Forty-five immune checkpoints and 23 immune features were considered statistically significant in HCC, while only seven immune checkpoints and one immune feature in ICC. Then, the prognostic value of DEAS events was studied, and two transcripts with different basic cell functions (proliferation, cell cycle, invasion, and migration) were produced by ADHFE1 through alternative splicing. Moreover, four nomograms were established in conjunction with relevant clinicopathological factors. Finally, we found two most significant splicing factors and further showed their protein crystal structure. The joint analysis of the AS events in HCC and ICC revealed novel insights into immune features and clinical prognosis, which might provide positive implications in HCC and ICC treatment.

Methylation Modification, Alternative Splicing, and Noncoding RNA Play a Role in Cancer Metastasis through Epigenetic Regulation

Metastasis is the leading cause of cancer-related deaths. Understanding the pathogenesis of metastasis at the molecular levels is of great significance for cancer research. However, the molecular diagnosis or treatment of cancer metastasis is limited. Accumulating and growing evidence shows that epigenetic changes are present in all human cancers, and epigenetic regulation is an indispensable factor to promote tumor metastasis. With the deepening of research and the advancement of technology, the function and mechanism of epigenetic regulation, including DNA methylation, histone/RNA modification, and precursor messenger RNA alternative splicing and noncoding RNAs, has become more increasingly clear. At present, the application of epigenetic therapies in tumor treatment is becoming a feasible therapeutic route. In this review, we looked for the key molecules in epigenetic regulation and discuss their relative regulating mechanisms in cancer metastasis. Furthermore, we highlight promising therapeutic strategies, including monitoring serum DNA for diagnostic purposes and early phase clinical trial therapies that target DNA and histone methylation. This may also be beneficial in finding new targets for further prognosis and diagnosis of cancer metastasis.

The biological Maxwell's demons: exploring ideas about the information processing in biological systems

This work is based on ideas supported by some of the biologists who discovered foundational facts of twentieth-century biology and who argued that Maxwell's demons are physically implemented by biological devices. In particular, JBS Haldane first, and later J. Monod, A, Lwoff and F. Jacob argued that enzymes and molecular receptors implemented Maxwell's demons that operate in systems far removed from thermodynamic equilibrium and that were responsible for creating the biological order. Later, these ideas were extended to other biological processes. In this article, we argue that these biological Maxwell's demons (BMD) are systems that have information processing capabilities that allow them to select their inputs and direct their outputs toward targets. In this context, we propose the idea that these BMD are information catalysts in which the processed information has broad thermodynamic consequences.

OncoSplicing: an updated database for clinically relevant alternative splicing in 33 human cancers

Alternative splicing (AS) represents a crucial method in mRNA level to regulate gene expression and contributes to the protein complexity. Abnormal splicing has been reported to play roles in several diseases, including cancers. We developed the OncoSplicing database for visualization of survival-associated and differential alternative splicing in 2019. Here, we provide an updated version of OncoSplicing for an integrative view of clinically relevant alternative splicing based on 122 423 AS events across 33 cancers in the TCGA SpliceSeq project and 238 558 AS events across 32 cancers in the TCGA SplAdder project. The new version of the database contains several useful features, such as annotation of alternative splicing-associated transcripts, survival analysis based on median and optimal cut-offs, differential analysis between TCGA tumour samples and adjacent normal samples or GTEx normal samples, pan-cancer views of alternative splicing, splicing differences and results of Cox’PH regression, identification of clinical indicator-relevant and cancer-specific splicing events, and downloadable splicing data in the SplAdder project. Overall, the substantially updated version of OncoSplicing (www.oncosplicing.com) is a user-friendly and registration-free database for browsing and searching clinically relevant alternative splicing in human cancers.

Prioritizing de novo autism risk variants with calibrated gene- and variant-scoring models

Whole-exome and whole-genome sequencing studies in autism spectrum disorder (ASD) have identified hundreds of thousands of exonic variants. Only a handful of them, primarily loss-of-function variants, have been shown to increase the risk for ASD, while the contributory roles of other variants, including most missense variants, remain unknown. New approaches that combine tissue-specific molecular profiles with patients' genetic data can thus play an important role in elucidating the functional impact of exonic variation and improve understanding of ASD pathogenesis. Here, we integrate spatio-temporal gene co-expression networks from the developing human brain and protein-protein interaction networks to first reach accurate prioritization of ASD risk genes based on their connectivity patterns with previously known high-confidence ASD risk genes. We subsequently integrate these gene scores with variant pathogenicity predictions to further prioritize individual exonic variants based on the positive-unlabeled learning framework with gene- and variant-score calibration. We demonstrate that this approach discriminates among variants between cases and controls at the high end of the prediction range. Finally, we experimentally validate our top-scoring de novo mutation NP_001243143.1:p.Phe309Ser in the sodium/potassium-transporting ATPase ATP1A3 to disrupt protein binding with different partners.

RNA splicing programs define tissue compartments and cell types at single-cell resolution

The extent splicing is regulated at single-cell resolution has remained controversial due to both available data and methods to interpret it. We apply the SpliZ, a new statistical approach, to detect cell-type-specific splicing in >110K cells from 12 human tissues. Using 10X Chromium data for discovery, 9.1% of genes with computable SpliZ scores are cell-type-specifically spliced, including ubiquitously expressed genes MYL6 and RPS24. These results are validated with RNA FISH, single-cell PCR, and Smart-seq2. SpliZ analysis reveals 170 genes with regulated splicing during human spermatogenesis, including examples conserved in mouse and mouse lemur. The SpliZ allows model-based identification of subpopulations indistinguishable based on gene expression, illustrated by subpopulation-specific splicing of classical monocytes involving an ultraconserved exon in SAT1. Together, this analysis of differential splicing across multiple organs establishes that splicing is regulated cell-type-specifically.

Detecting Rewiring Events in Protein-Protein Interaction Networks Based on Transcriptomic Data

Proteins rarely carry out their cellular functions in isolation. Instead, eukaryotic proteins engage in about six interactions with other proteins on average. The aggregated protein interactome of an organism forms a “hairy ball”-type protein-protein interaction (PPI) network. Yet, in a typical human cell, only about half of all proteins are expressed at a particular time. Hence, it has become common practice to prune the full PPI network to the subset of expressed proteins. If RNAseq data is available, one can further resolve the specific protein isoforms present in a cell or tissue. Here, we review various approaches, software tools and webservices that enable users to construct context-specific or tissue-specific PPI networks and how these are rewired between two cellular conditions. We illustrate their different functionalities on the example of the interactions involving the human TNR6 protein. In an outlook, we describe how PPI networks may be integrated with epigenetic data or with data on the activity of splicing factors.

Network module-based drug repositioning for pulmonary arterial hypertension

Pulmonary arterial hypertension (PAH) is a progressive disorder characterized by pulmonary vascular remodeling leading to increased pulmonary vascular resistance and pulmonary arterial pressure. PAH is a highly morbid cardiopulmonary disease adversely affecting lifespan and quality of life. Despite increased awareness and advances of medical therapies in recent decades, long-term prognosis and survival remain poor for patients with PAH. Novel therapies that can target the underlying pathobiology of PAH and reverse pulmonary vascular remodeling are clearly needed. In this study, we develop a network module-based framework to examine potential drug repositioning for PAH. The rationale for this approach is that in order to have therapeutic effects, the targets of potential drugs must be significantly proximate to the disease module of interest in the human protein-protein interactome. Based on 15 existing drugs for treating PAH, our framework integrates drug-drug interactions, drug-drug chemical similarity, drug targets, and PAH disease proteins into the human interactome, and prioritizes candidate drugs for PAH. We identified 53 drugs that could potentially be repurposed for PAH. Many of these candidates have strong literature support. Compared to black-box-like machine learning models, network module-based drug repositioning can provide mechanistic insights into how repositioned drugs can target the underlying pathobiological mechanisms of PAH.

Selection and Evaluation of a Thornless and HLB-Tolerant Bud-Sport of Pummelo Citrus With an Emphasis on Molecular Mechanisms

The selection of elite bud-sports is an important breeding approach in horticulture. We discovered and evaluated a thornless pummelo bud-sport (TL) that grew more vigorously and was more tolerant to Huanglongbing (HLB) than the thorny wild type (W). To reveal the underlying molecular mechanisms, we carried out whole-genome sequencing of W, and transcriptome comparisons of W, TL, and partially recovered thorny "mutants" (T). The results showed W, TL, and T varied in gene expression, allelic expression, and alternative splicing. Most genes/pathways with significantly altered expression in TL compared to W remained similarly altered in T. Pathway and gene ontology enrichment analysis revealed that the expression of multiple pathways, including photosynthesis and cell wall biosynthesis, was altered among the three genotypes. Remarkably, two polar auxin transporter genes, PIN7 and LAX3, were expressed at a significantly lower level in TL than in both W and T, implying alternation of polar auxin transport in TL may be responsible for the vigorous growth and thornless phenotype. Furthermore, 131 and 68 plant defense-related genes were significantly upregulated and downregulated, respectively, in TL and T compared with W. These genes may be involved in enhanced salicylic acid (SA) dependent defense and repression of defense inducing callose deposition and programmed cell death. Overall, these results indicated that the phenotype changes of the TL bud-sport were associated with tremendous transcriptome alterations, providing new clues and targets for breeding and gene editing for citrus improvement.

Full-length isoform transcriptome of the developing human brain provides further insights into autism

Alternative splicing plays an important role in brain development, but its global contribution to human neurodevelopmental diseases (NDDs) requires further investigation. Here we examine the relationships between splicing isoform expression in the brain and de novo loss-of-function mutations from individuals with NDDs. We analyze the full-length isoform transcriptome of the developing human brain and observe differentially expressed isoforms and isoform co-expression modules undetectable by gene-level analyses. These isoforms are enriched in loss-of-function mutations and microexons, are co-expressed with a unique set of partners, and have higher prenatal expression. We experimentally test the effect of splice-site mutations and demonstrate exon skipping in five NDD risk genes, including SCN2A, DYRK1A, and BTRC. Our results suggest that the splice site mutation in BTRC reduces translational efficiency, likely affecting Wnt signaling through impaired degradation of β-catenin. We propose that functional effects of mutations should be investigated at the isoform- rather than gene-level resolution.

Mutation Edgotype Drives Fitness Effect in Human

Missense mutations are known to perturb protein-protein interaction networks (known as interactome networks) in different ways. However, it remains unknown how different interactome perturbation patterns (“edgotypes”) impact organismal fitness. Here, we estimate the fitness effect of missense mutations with different interactome perturbation patterns in human, by calculating the fractions of neutral and deleterious mutations that do not disrupt PPIs (“quasi-wild-type”), or disrupt PPIs either by disrupting the binding interface (“edgetic”) or by disrupting overall protein stability (“quasi-null”). We first map pathogenic mutations and common non-pathogenic mutations onto homology-based three-dimensional structural models of proteins and protein-protein interactions in human. Next, we perform structure-based calculations to classify each mutation as either quasi-wild-type, edgetic, or quasi-null. Using our predicted as well as experimentally determined interactome perturbation patterns, we estimate that >∼40% of quasi-wild-type mutations are effectively neutral and the remaining are mostly mildly deleterious, that >∼75% of edgetic mutations are only mildly deleterious, and that up to ∼75% of quasi-null mutations may be strongly detrimental. These estimates are the first such estimates of fitness effect for different network perturbation patterns in any interactome. Our results suggest that while mutations that do not disrupt the interactome tend to be effectively neutral, the majority of human PPIs are under strong purifying selection and the stability of most human proteins is essential to human life.

Assessing conservation of alternative splicing with evolutionary splicing graphs

Understanding how protein function has evolved and diversified is of great importance for human genetics and medicine. Here, we tackle the problem of describing the whole transcript variability observed in several species by generalizing the definition of splicing graph. We provide a practical solution to construct parsimonious evolutionary splicing graphs where each node is a minimal transcript building block defined across species. We show a clear link between the functional relevance, tissue regulation, and conservation of alternative transcripts on a set of 50 genes. By scaling up to the whole human protein-coding genome, we identify a few thousand genes where alternative splicing modulates the number and composition of pseudorepeats. We have implemented our approach in ThorAxe, an efficient, versatile, robust, and freely available computational tool.

Systemic characterization of alternative splicing related to prognosis and immune infiltration in malignant mesothelioma

BACKGROUND

Malignant mesothelioma (MM) is a relatively rare and highly lethal tumor with few treatment options. Thus, it is important to identify prognostic markers that can help clinicians diagnose mesothelioma earlier and assess disease activity more accurately. Alternative splicing (AS) events have been recognized as critical signatures for tumor diagnosis and treatment in multiple cancers, including MM.

METHODS

We systematically examined the AS events and clinical information of 83 MM samples from TCGA database. Univariate Cox regression analysis was used to identify AS events associated with overall survival. LASSO analyses followed by multivariate Cox regression analyses were conducted to construct the prognostic signatures and assess the accuracy of these prognostic signatures by receiver operating characteristic (ROC) curve and Kaplan-Meier survival analyses. The ImmuCellAI and ssGSEA algorithms were used to assess the degrees of immune cell infiltration in MM samples. The survival-related splicing regulatory network was established based on the correlation between survival-related AS events and splicing factors (SFs).

RESULTS

A total of 3976 AS events associated with overall survival were identified by univariate Cox regression analysis, and ES events accounted for the greatest proportion. We constructed prognostic signatures based on survival-related AS events. The prognostic signatures proved to be an efficient predictor with an area under the curve (AUC) greater than 0.9. Additionally, the risk score based on 6 key AS events proved to be an independent prognostic factor, and a nomogram composed of 6 key AS events was established. We found that the risk score was significantly decreased in patients with the epithelioid subtype. In addition, unsupervised clustering clearly showed that the risk score was associated with immune cell infiltration. The abundances of cytotoxic T (Tc) cells, natural killer (NK) cells and T-helper 17 (Th17) cells were higher in the high-risk group, whereas the abundances of induced regulatory T (iTreg) cells were lower in the high-risk group. Finally, we identified 3 SFs (HSPB1, INTS1 and LUC7L2) that were significantly associated with MM patient survival and then constructed a regulatory network between the 3 SFs and survival-related AS to reveal potential regulatory mechanisms in MM.

CONCLUSION

Our study provided a prognostic signature based on 6 key events, representing a better effective tumor-specific diagnostic and prognostic marker than the TNM staging system. AS events that are correlated with the immune system may be potential therapeutic targets for MM.

DMIL-IsoFun: predicting isoform function using deep multi-instance learning

MOTIVATION

Alternative splicing creates the considerable proteomic diversity and complexity on relatively limited genome. Proteoforms translated from alternatively spliced isoforms of a gene actually execute the biological functions of this gene, which reflect the functional knowledge of genes at a finer granular level. Recently, some computational approaches have been proposed to differentiate isoform functions using sequence and expression data. However, their performance is far from being desirable, mainly due to the imbalance and lack of annotations at isoform-level, and the difficulty of modeling gene-isoform relations.

RESULT

We propose a deep multi-instance learning based framework (DMIL-IsoFun) to differentiate the functions of isoforms. DMIL-IsoFun firstly introduces a multi-instance learning convolution neural network trained with isoform sequences and gene-level annotations to extract the feature vectors and initialize the annotations of isoforms, and then uses a class-imbalance Graph Convolution Network to refine the annotations of individual isoforms based on the isoform co-expression network and extracted features. Extensive experimental results show that DMIL-IsoFun improves the Smin and Fmax of state-of-the-art solutions by at least 29.6% and 40.8%. The effectiveness of DMIL-IsoFun is further confirmed on a testbed of human multiple-isoform genes, and Maize isoforms related with photosynthesis.

AVAILABILITY

The code and data are available at http://www.sdu-idea.cn/codes.php?name=DMIL-Isofun.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

DoChaP: the domain change presenter

Alternative splicing results in multiple transcripts of the same gene, possibly encoding for different protein isoforms with different domains. Whereas it is possible to manually determine the effect of alternative splicing on the domain composition for a single event, the process requires the tedious integration of several data sources; it is error prone and not feasible for genome-wide characterization of domains affected by differential splicing. To fulfill the need for an automated solution, we developed the Domain Change Presenter (DoChaP, https://dochap.bgu.ac.il/), a web server for the visualization of exon-domain associations. DoChaP visualizes all transcripts of a given gene, the encoded proteins and their domains, and enables a comparison between the transcripts and between their protein products. The colors and organization make the structural effect of alternative splicing events on protein structures easily identified. To enable the study of the conservation of exons structure, alternative splicing, and the effect of alternative splicing on protein domains, DoChaP also provides a two-species comparison of exon-domain associations. DoChaP thus provides a unique and easy-to-use visualization of the exon-domain association and conservation, and will facilitate the study of the structural effects of alternative splicing in health and disease.

Systematic Profiling of mRNA Splicing Reveals the Prognostic Predictor and Potential Therapeutic Target for Glioblastoma Multiforme

Despite many changes in alternative splicing events (ASEs) are frequently involved in various cancers, prognosis-related ASEs and drug treatment targets in glioblastoma multiforme (GBM) have not been well explored. ASEs participate in many biological behaviors in the initiation and progression of tumors, the aberrant ASE has been considered another hallmark of cancer, and the systematic study of alternative splicing may provide potential biomarkers for malignancies. In this study, we carried out a systematic analysis to characterize the ASE signatures in GBM cohort. Through comparing GBM tissues and nontumor tissues, a total of 48,191 differently expressed ASEs from 10,727 genes were obtained, and these aberrant ASEs play an important role in the oncogenic process. Then, we identified 514 ASEs independently associated with patient survival in GBM by univariate and multivariate Cox regression, including exon skip in CD3D, alternate acceptor site in POLD2, and exon skip in DCN. Those prognostic models built on ASEs of each splice type can accurately predict the outcome of GBM patients, and values for the area under curve were 0.97 in the predictive model based on alternate acceptor site. In addition, the splicing-regulatory network revealed an interesting correlation between survival-associated splicing factors and prognostic ASE corresponding genes. Moreover, these three hub splicing factors in splicing regulation network are the potential targets of some drugs. In conclusion, a systematic analysis of ASE signatures in GBM could serve as an indicator for identifying novel prognostic biomarkers and guiding clinical treatment.

Applications of western blot technique: From bench to bedside

Western blot (WB) or immunoblot is a workhorse method. It is commonly used by biologists for study of different aspects of protein biomolecules. In addition, it has been widely used in disease diagnosis. Despite some limitations such as long time, different applications of WB have not been limited. In the present review, we have summarized scientific and clinical applications of WB. In addition, we described some new generation of WB techniques.

Transcript isoforms of Reep6 have distinct functions in the retina

Much of the complexity of the eukaryotic cell transcriptome is due to the alternative splicing of mRNA. However, knowledge on how transcriptome complexity is translated into functional complexity remains limited. For example, although different isoforms of a gene may show distinct temporal and spatial expression patterns, it is largely unknown whether these isoforms encode proteins with distinct functions matching their expression pattern. In this report, we investigated the function and relationship of the two isoforms of Reep6, namely Reep6.1 and Reep6.2, in rod photoreceptor cells. These two isoforms result from the alternative splicing of exon 5 and show mutually exclusive expression patterns. Reep6.2 is the canonical isoform that is expressed in non-retinal tissues while Reep6.1 is the only expressed isoform in the adult retina. The Reep6.1 isoform-specific knockout mouse, Reep6E5/E5, is generated by deleting exon 5 and a homozygous deletion phenotypically displayed a rod degeneration phenotype comparable to a Reep6 full knockout mouse, indicating that the Reep6.1 isoform is essential for the rod photoreceptor cell survival. Consistent with the results obtained from a loss-of-function experiment, overexpression of Reep6.2 failed to rescue the rod degeneration phenotype of Reep6 knockout mice while overexpression of Reep6.1 does lead to rescue. These results demonstrate that, consistent with the expression pattern of the isoform, Reep6.1 has rod-specific functions that cannot be substituted by its canonical isoform. Our findings suggested that a strict regulation of splicing is required for the maintenance of photoreceptor cells.

Long noncoding RNA DGCR5 involves in tumorigenesis of esophageal squamous cell carcinoma via SRSF1-mediated alternative splicing of Mcl-1

Long noncoding RNAs (lncRNAs) emerge as essential roles in the regulation of alternative splicing (AS) in various malignancies. Serine- and arginine-rich splicing factor 1 (SRSF1)-mediated AS events are the most important molecular hallmarks in cancer. Nevertheless, the biological mechanism underlying tumorigenesis of lncRNAs correlated with SRSF1 in esophageal squamous cell carcinoma (ESCC) remains elusive. In this study, we found that lncRNA DiGeorge syndrome critical region gene 5 (DGCR5) was upregulated in ESCC clinical samples, which associated with poor prognosis. Through RNA interference and overexpression approaches, we confirmed that DGCR5 contributed to promote ESCC cell proliferation, migration, and invasion while inhibited apoptosis in vitro. Mechanistically, DGCR5 could directly bind with SRSF1 to increase its stability and thus stimulate alternative splicing events. Furthermore, we clarified that SRSF1 regulated the aberrant splicing of myeloid cell leukemia-1 (Mcl-1) and initiated a significant Mcl-1L (antiapoptotic) isoform switch, which contributed to the expression of the full length of Mcl-1. Moreover, the cell-derived xenograft (CDX) model was validated that DGCR5 could facilitate the tumorigenesis of ESCC in vivo. Collectively, our findings identified that the key biological role of lncRNA DGCR5 in alternative splicing regulation and emphasized DGCR5 as a potential biomarker and therapeutic target for ESCC.

Structural basis of the selective activation of enzyme isoforms: Allosteric response to activators of β1- and β2-containing AMPK complexes

AMP-activated protein kinase (AMPK) is a key energy sensor regulating the cell metabolism in response to energy supply and demand. The evolutionary adaptation of AMPK to different tissues is accomplished through the expression of distinct isoforms that can form up to 12 complexes, which exhibit notable differences in the sensitivity to allosteric activators. To shed light into the molecular determinants of the allosteric regulation of this energy sensor, we have examined the structural and dynamical properties of β1- and β2-containing AMPK complexes formed with small molecule activators A-769662 and SC4, and dissected the mechanical response leading to active-like enzyme conformations through the analysis of interaction networks between structural domains. The results reveal the mechanical sensitivity of the α2β1 complex, in contrast with a larger resilience of the α2β2 species, especially regarding modulation by A-769662. Furthermore, binding of activators to α2β1 consistently promotes the pre-organization of the ATP-binding site, favoring the adoption of activated states of the enzyme. These findings are discussed in light of the changes in the residue content of β-subunit isoforms, particularly regarding the β1Asn111 → β2Asp111 substitution as a key factor in modulating the mechanical sensitivity of β1- and β2-containing AMPK complexes. Our studies pave the way for the design of activators tailored for improving the therapeutic treatment of tissue-specific metabolic disorders.

IsoResolve: predicting splice isoform functions by integrating gene and isoform-level features with domain adaptation

MOTIVATION

High resolution annotation of gene functions is a central goal in functional genomics. A single gene may produce multiple isoforms with different functions through alternative splicing. Conventional approaches, however, consider a gene as a single entity without differentiating these functionally different isoforms. Towards understanding gene functions at higher resolution, recent efforts have focused on predicting the functions of isoforms. However, the performance of existing methods is far from satisfactory mainly because of the lack of isoform-level functional annotation.

RESULTS

We present IsoResolve, a novel approach for isoform function prediction, which leverages the information from gene function prediction models with domain adaptation (DA). IsoResolve treats gene-level and isoform-level features as source and target domains, respectively. It uses DA to project the two domains into a latent variable space in such a way that the latent variables from the two domains have similar distribution, which enables the gene domain information to be leveraged for isoform function prediction. We systematically evaluated the performance of IsoResolve in predicting functions. Compared with five state-of-the-art methods, IsoResolve achieved significantly better performance. IsoResolve was further validated by case studies of genes with isoform-level functional annotation.

AVAILABILITY AND IMPLEMENTATION

IsoResolve is freely available at https://github.com/genemine/IsoResolve.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Alternative splicing and gene co-expression network-based analysis of dizygotic twins with autism-spectrum disorder and their parents

Autism spectrum disorder (ASD) is a neurodevelopmental disorder with high heritability, however, understanding the complexity of the underlying genetic basis has proven to be a challenging task. We hypothesized that dissecting the aberrations in alternative splicing (AS) and their effects on expression networks might provide insight. Therefore, we performed AS and co-expression analyses of total RNA isolated from Peripheral Blood Mononuclear Cells (PBMCs) of two pairs of dizygotic (DZ) twins with non-syndromic autism and their parents. We identified 183 differential AS events in 146 genes, seven of them being Simons Foundation Autism Research Initiative (SFARI) Category 1-3 genes, three of which had previously been reported to be alternatively spliced in ASD post-mortem brains. Gene co-expression analysis identified 7 modules with 513 genes, 5 of which were SFARI Category 1 or Category 2 genes. Among differentially AS genes within the modules, ZNF322 and NR4A1 could be potentially interesting targets for further investigations.

HNRNPM controls circRNA biogenesis and splicing fidelity to sustain cancer cell fitness

High spliceosome activity is a dependency for cancer cells, making them more vulnerable to perturbation of the splicing machinery compared to normal cells. To identify splicing factors important for prostate cancer (PCa) fitness, we performed pooled shRNA screens in vitro and in vivo. Our screens identified heterogeneous nuclear ribonucleoprotein M (HNRNPM) as a regulator of PCa cell growth. RNA- and eCLIP-sequencing identified HNRNPM binding to transcripts of key homeostatic genes. HNRNPM binding to its targets prevents aberrant exon inclusion and backsplicing events. In both linear and circular mis-spliced transcripts, HNRNPM preferentially binds to GU-rich elements in long flanking proximal introns. Mimicry of HNRNPM-dependent linear-splicing events using splice-switching-antisense-oligonucleotides was sufficient to inhibit PCa cell growth. This suggests that PCa dependence on HNRNPM is likely a result of mis-splicing of key homeostatic coding and non-coding genes. Our results have further been confirmed in other solid tumors. Taken together, our data reveal a role for HNRNPM in supporting cancer cell fitness. Inhibition of HNRNPM activity is therefore a potential therapeutic strategy in suppressing growth of PCa and other solid tumors.

Capillary Zone Electrophoresis-Electron-Capture Collision-Induced Dissociation on a Quadrupole Time-of-Flight Mass Spectrometer for Top-Down Characterization of Intact Proteins

Mass spectrometry (MS)-based denaturing top-down proteomics (dTDP) requires high-capacity separation and extensive gas-phase fragmentation of proteoforms. Herein, we coupled capillary zone electrophoresis (CZE) to electron-capture collision-induced dissociation (ECciD) on an Agilent 6545 XT quadrupole time-of-flight (Q-TOF) mass spectrometer for dTDP for the first time. During ECciD, the protein ions were first fragmented using ECD, followed by further activation and fragmentation by applying a CID potential. In this pilot study, we optimized the CZE-ECciD method for small proteins (lower than 20 kDa) regarding the charge state of protein parent ions for fragmentation and the CID potential applied to maximize the protein backbone cleavage coverage and the number of sequence-informative fragment ions. The CZE-ECciD Q-TOF platform provided extensive backbone cleavage coverage for three standard proteins lower than 20 kDa from only single charge states in a single CZE-MS/MS run in the targeted MS/MS mode, including ubiquitin (97%, +7, 8.6 kDa), superoxide dismutase (SOD, 87%, +17, 16 kDa), and myoglobin (90%, +16, 17 kDa). The CZE-ECciD method produced comparable cleavage coverage of small proteins (i.e., myoglobin) with direct-infusion MS studies using electron transfer dissociation (ETD), activated ion-ETD, and combinations of ETD and collision-based fragmentation on high-end orbitrap mass spectrometers. The results render CZE-ECciD a new tool for dTDP to enhance both separation and gas-phase fragmentation of proteoforms.

LncRNA Ctcflos orchestrates transcription and alternative splicing in thermogenic adipogenesis

The recruitment of thermogenic brite adipocytes within white adipose tissue attenuates obesity and metabolic comorbidities, arousing interest in understanding the underlying regulatory mechanisms. The molecular network of brite adipogenesis, however, remains largely unresolved. In this light, long noncoding RNAs (lncRNAs) emerged as a versatile class of modulators that control many steps within the differentiation machinery. Leveraging the naturally varying propensities of different inbred mouse strains for white adipose tissue browning, we identify the nuclear lncRNA Ctcflos as a pivotal orchestrator of thermogenic gene expression during brite adipocyte differentiation. Mechanistically, Ctcflos acts as a pleiotropic regulator, being essential for the transcriptional recruitment of the early core thermogenic regulatory program and the modulation of alternative splicing to drive brite adipogenesis. This is showcased by Ctcflos‐regulated gene transcription and splicing of the key browning factor Prdm16 toward the isoform that is specific for the thermogenic gene program. Conclusively, our findings emphasize the mechanistic versatility of lncRNAs acting at several independent levels of gene expression for effective regulation of key differentiation factors to direct cell fate and function.

RBM10: Structure, functions, and associated diseases

RBM10 is a nuclear RNA-binding protein (RBP) that regulates the alternative splicing of primary transcripts. Recently, research on RBM10 has become increasingly active owing to its clinical importance, as indicated by studies on RBM0 mutations that cause TARP syndrome, an X-linked congenital pleiotropic developmental anomaly, and various cancers such as lung adenocarcinoma in adults. Herein, the molecular biology of RBM10 and its significance in medicine are reviewed, focusing on the gene and protein structures of RBM10, its cell biology, molecular functions and regulation, relationship with the paralogous protein RBM5, and the mutations of RBM10 and their associated diseases. Finally, the challenges in future studies of RBM10 are discussed in the concluding remarks.

Assessing the functional relevance of splice isoforms

Alternative splicing of messenger RNA can generate an array of mature transcripts, but it is not clear how many go on to produce functionally relevant protein isoforms. There is only limited evidence for alternative proteins in proteomics analyses and data from population genetic variation studies indicate that most alternative exons are evolving neutrally. Determining which transcripts produce biologically important isoforms is key to understanding isoform function and to interpreting the real impact of somatic mutations and germline variations. Here we have developed a method, TRIFID, to classify the functional importance of splice isoforms. TRIFID was trained on isoforms detected in large-scale proteomics analyses and distinguishes these biologically important splice isoforms with high confidence. Isoforms predicted as functionally important by the algorithm had measurable cross species conservation and significantly fewer broken functional domains. Additionally, exons that code for these functionally important protein isoforms are under purifying selection, while exons from low scoring transcripts largely appear to be evolving neutrally. TRIFID has been developed for the human genome, but it could in principle be applied to other well-annotated species. We believe that this method will generate valuable insights into the cellular importance of alternative splicing.

Characterization of alternative splicing events and prognostic signatures in breast cancer

Background

Breast cancer (BRCA) is one of the most common cancers worldwide. Abnormal alternative splicing (AS) frequently observed in cancers. This study aims to demonstrate AS events and signatures that might serve as prognostic indicators for BRCA.

Methods

Original data for all seven types of splice events were obtained from TCGA SpliceSeq database. RNA-seq and clinical data of BRCA cohorts were downloaded from TCGA database. Survival-associated AS events in BRCA were analyzed by univariate COX proportional hazards regression model. Prognostic signatures were constructed for prognosis prediction in patients with BRCA based on survival-associated AS events. Pearson correlation analysis was performed to measure the correlation between the expression of splicing factors (SFs) and the percent spliced in (PSI) values of AS events. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were conducted to demonstrate pathways in which survival-associated AS event is enriched.

Results

A total of 45,421 AS events in 21,232 genes were identified. Among them, 1121 AS events in 931 genes significantly correlated with survival for BRCA. The established AS prognostic signatures of seven types could accurately predict BRCA prognosis. The comprehensive AS signature could serve as independent prognostic factor for BRCA. A SF-AS regulatory network was therefore established based on the correlation between the expression levels of SFs and PSI values of AS events.

Conclusions

This study revealed survival-associated AS events and signatures that may help predict the survival outcomes of patients with BRCA. Additionally, the constructed SF-AS networks in BRCA can reveal the underlying regulatory mechanisms in BRCA.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08305-6.

A drug-repositioning screen using splicing-sensitive fluorescent reporters identifies novel modulators of VEGF-A splicing with anti-angiogenic properties

Alternative splicing of the vascular endothelial growth factor A (VEGF-A) terminal exon generates two protein families with differing functions. Pro-angiogenic VEGF-A_xxxa isoforms are produced via selection of the proximal 3' splice site of the terminal exon. Use of an alternative distal splice site generates the anti-angiogenic VEGF-A_xxxb proteins. A bichromatic splicing-sensitive reporter was designed to mimic VEGF-A alternative splicing and was used as a molecular tool to further investigate this alternative splicing event. Part of VEGF-A's terminal exon and preceding intron were inserted into a minigene construct followed by the coding sequences for two fluorescent proteins. A different fluorescent protein is expressed depending on which 3' splice site of the exon is used during splicing (dsRED denotes VEGF-A_xxxa and EGFP denotes VEGF-A_xxxb). The fluorescent output can be used to follow splicing decisions in vitro and in vivo. Following successful reporter validation in different cell lines and altering splicing using known modulators, a screen was performed using the LOPAC library of small molecules. Alterations to reporter splicing were measured using a fluorescent plate reader to detect dsRED and EGFP expression. Compounds of interest were further validated using flow cytometry and assessed for effects on endogenous VEGF-A alternative splicing at the mRNA and protein level. Ex vivo and in vitro angiogenesis assays were used to demonstrate the anti-angiogenic effect of the compounds. Furthermore, anti-angiogenic activity was investigated in a Matrigel in vivo model. To conclude, we have identified a set of compounds that have anti-angiogenic activity through modulation of VEGF-A terminal exon splicing.

FOXM1D potentiates PKM2-mediated tumor glycolysis and angiogenesis

Tumor growth, especially in the late stage, requires adequate nutrients and rich vasculature, in which PKM2 plays a convergent role. It has been reported that PKM2, together with FOXM1D, is upregulated in late-stage colorectal cancer and associated with metastasis; however, their underlying mechanism for promoting tumor progression remains elusive. Herein, we revealed that FOXM1D potentiates PKM2-mediated glycolysis and angiogenesis through multiple protein-protein interactions. In the presence of FBP, FOXM1D binds to tetrameric PKM2 and assembles a heterooctamer, restraining PKM2 metabolic activity by about a half and thereby promoting aerobic glycolysis. Furthermore, FOXM1D interacts with PKM2 and NF-κB and induces their nuclear translocation with the assistance of the nuclear transporter importin 4. Once in the nucleus, PKM2 and NF-κB complexes subsequently augment VEGFA transcription. The increased VEGFA is secreted extracellularly via exosomes, an event potentiated by the interaction of FOXM1 with VPS11, eventually promoting tumor angiogenesis. Based on these findings, our study provides another insight into the role of PKM2 in the regulation of glycolysis and angiogenesis.

Exons as units of phenotypic impact for truncating mutations in autism

Autism spectrum disorders (ASD) are a group of related neurodevelopmental diseases displaying significant genetic and phenotypic heterogeneity. Despite recent progress in understanding ASD genetics, the nature of phenotypic heterogeneity across probands remains unclear. Notably, likely gene-disrupting (LGD) de novo mutations affecting the same gene often result in substantially different ASD phenotypes. Nevertheless, we find that truncating mutations affecting the same exon frequently lead to strikingly similar intellectual phenotypes in unrelated ASD probands. Analogous patterns are observed for two independent proband cohorts and several other important ASD-associated phenotypes. We find that exons biased toward prenatal and postnatal expression preferentially contribute to ASD cases with lower and higher IQ phenotypes, respectively. These results suggest that exons, rather than genes, often represent a unit of effective phenotypic impact for truncating mutations in autism. The observed phenotypic patterns are likely mediated by nonsense-mediated decay (NMD) of splicing isoforms, with autism phenotypes usually triggered by relatively mild (15-30%) decreases in overall gene dosage. We find that each ASD gene with recurrent mutations can be characterized by a parameter, phenotype dosage sensitivity (PDS), which quantifies the relationship between changes in a gene's dosage and changes in a given disease phenotype. We further demonstrate analogous relationships between exon LGDs and gene expression changes in multiple human tissues. Therefore, similar phenotypic patterns may be also observed in other human genetic disorders.

Control Analysis of Protein-Protein Interaction Network Reveals Potential Regulatory Targets for MYCN

BACKGROUND

MYCN is an oncogenic transcription factor of the MYC family and plays an important role in the formation of tissues and organs during development before birth. Due to the difficulty in drugging MYCN directly, revealing the molecules in MYCN regulatory networks will help to identify effective therapeutic targets.

METHODS

We utilized network controllability theory, a recent developed powerful tool, to identify the potential drug target around MYCN based on Protein-Protein interaction network of MYCN. First, we constructed a Protein-Protein interaction network of MYCN based on public databases. Second, network control analysis was applied on network to identify driver genes and indispensable genes of the MYCN regulatory network. Finally, we developed a novel integrated approach to identify potential drug targets for regulating the function of the MYCN regulatory network.

RESULTS

We constructed an MYCN regulatory network that has 79 genes and 129 interactions. Based on network controllability theory, we analyzed driver genes which capable to fully control the network. We found 10 indispensable genes whose alternation will significantly change the regulatory pathways of the MYCN network. We evaluated the stability and correlation analysis of these genes and found EGFR may be the potential drug target which closely associated with MYCN.

CONCLUSION

Together, our findings indicate that EGFR plays an important role in the regulatory network and pathways of MYCN and therefore may represent an attractive therapeutic target for cancer treatment.

Regulatory roles and mechanisms of alternative RNA splicing in adipogenesis and human metabolic health

Alternative splicing (AS) regulates gene expression patterns at the post-transcriptional level and generates a striking expansion of coding capacities of genomes and cellular protein diversity. RNA splicing could undergo modulation and close interaction with genetic and epigenetic machinery. Notably, during the adipogenesis processes of white, brown and beige adipocytes, AS tightly interplays with the differentiation gene program networks. Here, we integrate the available findings on specific splicing events and distinct functions of different splicing regulators as examples to highlight the directive biological contribution of AS mechanism in adipogenesis and adipocyte biology. Furthermore, accumulating evidence has suggested that mutations and/or altered expression in splicing regulators and aberrant splicing alterations in the obesity-associated genes are often linked to humans’ diet-induced obesity and metabolic dysregulation phenotypes. Therefore, significant attempts have been finally made to overview novel detailed discussion on the prospects of splicing machinery with obesity and metabolic disorders to supply featured potential management mechanisms in clinical applicability for obesity treatment strategies.

Differential alternative RNA splicing and transcription events between tumors from African American and White patients in The Cancer Genome Atlas

Individuals of African ancestry suffer disproportionally from higher incidence, aggressiveness, and mortality for particular cancers. This disparity likely results from an interplay among differences in multiple determinants of health, including differences in tumor biology. We used The Cancer Genome Atlas (TCGA) SpliceSeq and TCGA aggregate expression datasets and identified differential alternative RNA splicing and transcription events (ARS/T) in cancers between self-identified African American (AA) and White (W) patients. We found that retained intron events were enriched among race-related ARS/T. In addition, on average, 12% of the most highly ranked race-related ARS/T overlapped between any two analyzed cancers. Moreover, the genes undergoing race-related ARS/T functioned in cancer-promoting pathways, and a number of race-related ARS/T were associated with patient survival. We built a web-application, CanSplice, to mine genomic datasets by self-identified race. The race-related targets have the potential to aid in the development of new biomarkers and therapeutics to mitigate cancer disparity.

Oncogenic translation directs spliceosome dynamics revealing an integral role for SF3A3 in breast cancer

Splicing is a central RNA-based process commonly altered in human cancers; however, how spliceosomal components are co-opted during tumorigenesis remains poorly defined. Here we unravel the core splice factor SF3A3 at the nexus of a translation-based program that rewires splicing during malignant transformation. Upon MYC hyperactivation, SF3A3 levels are modulated translationally through an RNA stem-loop in an eIF3D-dependent manner. This ensures accurate splicing of mRNAs enriched for mitochondrial regulators. Altered SF3A3 translation leads to metabolic reprogramming and stem-like properties that fuel MYC tumorigenic potential in vivo. Our analysis reveals that SF3A3 protein levels predict molecular and phenotypic features of aggressive human breast cancers. These findings unveil a post-transcriptional interplay between splicing and translation that governs critical facets of MYC-driven oncogenesis.

Biological impact of mutually exclusive exon switching

Alternative splicing can expand the diversity of proteomes. Homologous mutually exclusive exons (MXEs) originate from the same ancestral exon and result in polypeptides with similar structural properties but altered sequence. Why would some genes switch homologous exons and what are their biological impact? Here, we analyse the extent of sequence, structural and functional variability in MXEs and report the first large scale, structure-based analysis of the biological impact of MXE events from different genomes. MXE-specific residues tend to map to single domains, are highly enriched in surface exposed residues and cluster at or near protein functional sites. Thus, MXE events are likely to maintain the protein fold, but alter specificity and selectivity of protein function. This comprehensive resource of MXE events and their annotations is available at: http://gene3d.biochem.ucl.ac.uk/mxemod/. These findings highlight how small, but significant changes at critical positions on a protein surface are exploited in evolution to alter function.

Alternatively Splicing Interactomes Identify Novel Isoform-Specific Partners for NSD2

Nuclear receptor SET domain protein (NSD2) plays a fundamental role in the pathogenesis of Wolf-Hirschhorn Syndrome (WHS) and is overexpressed in multiple human myelomas, but its protein-protein interaction (PPI) patterns, particularly at the isoform/exon levels, are poorly understood. We explored the subcellular localizations of four representative NSD2 transcripts with immunofluorescence microscopy. Next, we used label-free quantification to perform immunoprecipitation mass spectrometry (IP-MS) analyses of the transcripts. Using the interaction partners for each transcript detected in the IP-MS results, we identified 890 isoform-specific PPI partners (83% are novel). These PPI networks were further divided into four categories of the exon-specific interactome. In these exon-specific PPI partners, two genes, RPL10 and HSPA8, were successfully confirmed by co-immunoprecipitation and Western blotting. RPL10 primarily interacted with Isoforms 1, 3, and 5, and HSPA8 interacted with all four isoforms, respectively. Using our extended NSD2 protein interactions, we constructed an isoform-level PPI landscape for NSD2 to serve as reference interactome data for NSD2 spliceosome-level studies. Furthermore, the RNA splicing processes supported by these isoform partners shed light on the diverse roles NSD2 plays in WHS and myeloma development. We also validated the interactions using Western blotting, RPL10, and the three NSD2 (Isoform 1, 3, and 5). Our results expand gene-level NSD2 PPI networks and provide a basis for the treatment of NSD2-related developmental diseases.

The functional analysis of ABCG transporters in the adaptation of pigeon pea (Cajanus cajan) to abiotic stresses

ATP-binding cassette (ABC) transporters are a class of proteins found in living organisms that mediate transmembrane transport by hydrolyzing ATP. They play a vital role in the physiological processes of growth and development in plants. The most numerous sub-type transporter in the ABC transporter family is the ABCG group and which have the most complex function in a plant’s response to abiotic stresses. Our study focused on the effect of ABCG transporters in the adaptation of the pigeon pea to adverse environments (such as drought, salt, temperature, etc.). We conducted a functional analysis of ABCG transporters in the pigeon pea and their role in response to abiotic stresses. A total of 51 ABCG genes (CcABCGs) were identified, and phylogenetic analysis was conducted. We also identified the physicochemical properties of the encoded proteins, predicted their subcellular localization, and identified of the conserved domains. Expression analysis showed that ABCG genes have different expression profiles with tissues and abiotic stresses. Our results showed that CcABCG28 was up-regulated at low temperatures, and CcABCG7 was up-regulated with drought and aluminum stress. The initial results revealed that ABCG transporters are more effective in the abiotic stress resistance of pigeon peas, which improves our understanding of their application in abiotic stress resistance.

e-MutPath: computational modeling reveals the functional landscape of genetic mutations rewiring interactome networks

Understanding the functional impact of cancer somatic mutations represents a critical knowledge gap for implementing precision oncology. It has been increasingly appreciated that the interaction profile mediated by a genomic mutation provides a fundamental link between genotype and phenotype. However, specific effects on biological signaling networks for the majority of mutations are largely unknown by experimental approaches. To resolve this challenge, we developed e-MutPath (edgetic Mutation-mediated Pathway perturbations), a network-based computational method to identify candidate ‘edgetic’ mutations that perturb functional pathways. e-MutPath identifies informative paths that could be used to distinguish disease risk factors from neutral elements and to stratify disease subtypes with clinical relevance. The predicted targets are enriched in cancer vulnerability genes, known drug targets but depleted for proteins associated with side effects, demonstrating the power of network-based strategies to investigate the functional impact and perturbation profiles of genomic mutations. Together, e-MutPath represents a robust computational tool to systematically assign functions to genetic mutations, especially in the context of their specific pathway perturbation effect.

Gene co-expression in the interactome: moving from correlation toward causation via an integrated approach to disease module discovery

In this study, we integrate the outcomes of co-expression network analysis with the human interactome network to predict novel putative disease genes and modules. We first apply the SWItch Miner (SWIM) methodology, which predicts important (switch) genes within the co-expression network that regulate disease state transitions, then map them to the human protein-protein interaction network (PPI, or interactome) to predict novel disease-disease relationships (i.e., a SWIM-informed diseasome). Although the relevance of switch genes to an observed phenotype has been recently assessed, their performance at the system or network level constitutes a new, potentially fascinating territory yet to be explored. Quantifying the interplay between switch genes and human diseases in the interactome network, we found that switch genes associated with specific disorders are closer to each other than to other nodes in the network, and tend to form localized connected subnetworks. These subnetworks overlap between similar diseases and are situated in different neighborhoods for pathologically distinct phenotypes, consistent with the well-known topological proximity property of disease genes. These findings allow us to demonstrate how SWIM-based correlation network analysis can serve as a useful tool for efficient screening of potentially new disease gene associations. When integrated with an interactome-based network analysis, it not only identifies novel candidate disease genes, but also may offer testable hypotheses by which to elucidate the molecular underpinnings of human disease and reveal commonalities between seemingly unrelated diseases.

DIGGER: exploring the functional role of alternative splicing in protein interactions

Alternative splicing plays a major role in regulating the functional repertoire of the proteome. However, isoform-specific effects to protein-protein interactions (PPIs) are usually overlooked, making it impossible to judge the functional role of individual exons on a systems biology level. We overcome this barrier by integrating protein-protein interactions, domain-domain interactions and residue-level interactions information to lift exon expression analysis to a network level. Our user-friendly database DIGGER is available at https://exbio.wzw.tum.de/digger and allows users to seamlessly switch between isoform and exon-centric views of the interactome and to extract sub-networks of relevant isoforms, making it an essential resource for studying mechanistic consequences of alternative splicing.

A splicing factor switch controls hematopoietic lineage specification of pluripotent stem cells

Alternative splicing (AS) leads to transcriptome diversity in eukaryotic cells and is one of the key regulators driving cellular differentiation. Although AS is of crucial importance for normal hematopoiesis and hematopoietic malignancies, its role in early hematopoietic development is still largely unknown. Here, by using high-throughput transcriptomic analyses, we show that pervasive and dynamic AS takes place during hematopoietic development of human pluripotent stem cells (hPSCs). We identify a splicing factor switch that occurs during the differentiation of mesodermal cells to endothelial progenitor cells (EPCs). Perturbation of this switch selectively impairs the emergence of EPCs and hemogenic endothelial progenitor cells (HEPs). Mechanistically, an EPC-induced alternative spliced isoform of NUMB dictates EPC specification by controlling NOTCH signaling. Furthermore, we demonstrate that the splicing factor SRSF2 regulates splicing of the EPC-induced NUMB isoform, and the SRSF2-NUMB-NOTCH splicing axis regulates EPC generation. The identification of this splicing factor switch provides a new molecular mechanism to control cell fate and lineage specification.

Noncoding RNAs regulate alternative splicing in Cancer

AS (alternative splicing) is a fundamental process by which a gene can generate multiple distinct mRNA transcripts to increase protein diversity. Defects in AS influence the occurrence and development of many diseases, including cancers, and are frequently found to participate in various aspects of cancer biology, such as promoting invasion, metastasis, apoptosis resistance and drug resistance. NcRNAs (noncoding RNAs) are an abundant class of RNAs that do not encode proteins. NcRNAs include miRNAs (microRNAs), lncRNAs (long noncoding RNAs), circRNAs (circular RNAs) and snRNAs (small nuclear RNAs) and have been proven to act as regulatory molecules that mediate cancer processes through AS. NcRNAs can directly or indirectly influence a plethora of molecular targets to regulate cis-acting elements, trans-acting factors, or pre-mRNA transcription at multiple levels, affecting the AS process and generating alternatively spliced isoforms. Consequently, ncRNA-mediated AS outcomes affect multiple cellular signaling pathways that promote or suppress cancer progression. In this review, we summarize the current mechanisms by which ncRNAs regulate AS in cancers and discuss their potential clinical applications as biomarkers and therapeutic targets.

Prognostic Alternative mRNA Splicing in Adrenocortical Carcinoma

BACKGROUND

This paper aims to identify alternative RNA splicing landscape and its prognostic value in adrenocortical carcinoma.

METHODS

The alternative splicing events data with corresponding clinical information data of 79 ACC patients were obtained from the Cancer Genome Atlas and SpliceSeq package. Prognosis-associated AS events by using univariate Cox regression analysis were selected. Gene functional enrichment analysis demonstrated the potential pathways enriched by survival-associated AS. Prognosis-related splicing events were submitted to develop moderate predictors using Lasso regression model.

RESULTS

One thousand five survival-associated alternative splicing events were identified. The prognostic genes included ATXN2L, MEIS1, IKBKB, COX4I1. Functional enrichment analysis suggested that prognostic splicing events are associated with Wnt signaling pathway. A prediction model including 12 alternative splicing events was constructed by Lasso regression using train set. ROC analysis showed good performance of the prediction model in test set. Then, a nomogram integrating the clinical-pathological factors and riskscore was constructed for predicting 1- and 3-year survival rate.

CONCLUSION

Our data provide a comprehensive bioinformatics analysis of AS events in ACC, providing biomarkers for disease progression and a potentially rich source of novel therapeutic targets.

Building genomes to understand biology

Genetic manipulation is one of the central strategies that biologists use to investigate the molecular underpinnings of life and its diversity. Thus, advances in genetic manipulation usually lead to a deeper understanding of biological systems. During the last decade, the construction of chromosomes, known as synthetic genomics, has emerged as a novel approach to genetic manipulation. By facilitating complex modifications to chromosome content and structure, synthetic genomics opens new opportunities for studying biology through genetic manipulation. Here, we discuss different classes of genetic manipulation that are enabled by synthetic genomics, as well as biological problems they each can help solve.

Survival-Associated Alternative Splicing Events and Prognostic Signatures in Pancreatic Cancer

Background

Alternative splicing (AS) is reported to be related to the biological process of multiple malignancies. This study is conducted to identify survival-associated AS events and prognostic signatures that may serve as prognostic indicators for pancreatic cancer (PC).

Methods

Univariate Cox analysis was used to determine the survival-associated AS events in PC. Prognostic signatures were constructed by LASSO Cox analysis based on seven types of survival-associated AS events. The correlation between the expression of splicing factors (SFs) and the percent spliced in values of AS events was analyzed by Pearson correlation analysis. Risk scores were calculated to determine high- or low-risk patients with different types of AS events. Gene functional annotation analysis was performed to reveal pathways in which prognostic AS is enriched.

Results

A total of 45,313 AS events in 10,624 genes were observed, and there were 1,565 AS events in 1,223 genes significantly correlated with overall survival for PC. Kaplan–Meier analysis, receiver-operator characteristic curve, univariate and multivariate Cox analyses showed that AS prognostic signatures could effectively predict prognosis of patients with PC. Splicing factors–AS regulatory networks were established to demonstrate the interaction between AS events and SFs.

Conclusion

The survival-associated AS events and prognostic signatures identified in this study can serve as useful tool for predicting prognosis of patients with PC. Moreover, the SF–AS regulatory networks may provide clues for the mechanisms underlying AS in PC.

Systematic Analyses of the Differentially Expressed Alternative Splicing Events in Gastric Cancer and Its Clinical Significance

Accumulation of evidence has indicated a close relationship between alternative splicing (AS) and gastric cancer (GC), whereas systematic analyses of the differentially expressed AS events (DEAS) between GC and normal tissues are lacking. RNA-Seq data and the corresponding clinical information were downloaded from TCGA GC cohort. The percent spliced-in (PSI) value calculated in the GC tissues and normal tissues was employed to quantify the DEAS. Further, survival-associated DEAS and DEAS signatures were identified by univariate and multivariate cox regression analyses. To evaluate the association between DEAS and patients' clinical features, Kaplan-Meier analysis, receiver operator characteristic (ROC) curve, Cox proportional regression and nomograms incorporating the DEAS signatures were performed. DEAS and their splicing networks were finally analyzed by bioinformatics methods. In addition, we use the method of random grouping to divide the samples into the training group and the test group. The final results of the two groups are consistent. After strict filtering, a total of 44,935 AS events were identified, among which 11,141 DEAS were preliminarily screened from 5032 genes. A total of 454 DEAS was associated with OS, and 872 DEAS were associated with DFS. The final prognostic signatures were constructed from the survival-associated DEAS with an area under the receiver operating characteristic (ROC) curve (AUC) greater than 0.6. Only ES in ABI1 was simultaneously associated with OS and DFS. Finally, we identified the splicing correlation network between the prognostic splicing factors (SF) and DEAS in GC. Our study provided a systematic portrait of survival-associated DEAS in GC and uncovered splicing networks that are valuable in deciphering the underlying mechanisms of AS in GC.