Widespread Expansion of Protein Interaction Capabilities by Alternative Splicing

Personal Mutanomes Meet Modern Oncology Drug Discovery and Precision Health

Recent remarkable advances in genome sequencing have enabled detailed maps of identified and interpreted genomic variation, dubbed "mutanomes." The availability of thousands of exome/genome sequencing data has prompted the emergence of new challenges in the identification of novel druggable targets and therapeutic strategies. Typically, mutanomes are viewed as one- or two-dimensional. The three-dimensional protein structural view of personal mutanomes sheds light on the functional consequences of clinically actionable mutations revealed in tumor diagnosis and followed up in personalized treatments, in a mutanome-informed manner. In this review, we describe the protein structural landscape of personal mutanomes and provide expert opinions on rational strategies for more streamlined oncological drug discovery and molecularly targeted therapies for each individual and each tumor. We provide the structural mechanism of orthosteric versus allosteric drugs at the atom-level via targeting specific somatic alterations for combating drug resistance and the "undruggable" challenges in solid and hematologic neoplasias. We discuss computational biophysics strategies for innovative mutanome-informed cancer immunotherapies and combination immunotherapies. Finally, we highlight a personal mutanome infrastructure for the emerging development of personalized cancer medicine using a breast cancer case study.

OTUB2 stabilizes U2AF2 to promote the Warburg effect and tumorigenesis via the AKT/mTOR signaling pathway in non-small cell lung cancer

Increasing evidence has confirmed that deubiquitinating enzymes play an important role in lung cancer progression. In the current study, we investigated the expression profile of deubiquitinating enzymes in non-small cell lung cancer (NSCLC) tissues and identified OTUB2 as an upregulated deubiquitinating enzyme. The role of OTUB2 in NSCLC is unknown. Methods: Quantitative, real-time PCR and Western blot were used to detect OTUB2 and U2AF2 expression in NSCLC tissues. The correlations between OTUB2 and U2AF2 expression and clinicopathologic features were then analyzed. We used In vitro Cell Counting Kit-8 (CCK-8) , colony formation , and trans-well invasion assays to investigate the function of OTUB2 and U2AF2 in tumorigenesis. The regulation of glycolysis by OTUB2 and U2AF2 was assessed by determining the extracellular acid ratio, glucose consumption, and lactate production. The mechanism of OTUB2 was explored through co-immunoprecipitation and mass spectrometry analyses. A xenograft model was also used to study the tumorigenesis role of OTUB2 In vivo. Results: OTUB2 expression was significantly upregulated in primary NSCLC tissues and greatly associated with metastasis, advanced tumor stages, poor survival, and recurrence. In NSCLC cell lines, OTUB2 promoted cell growth, colony formation, migration, and invasive activities. Mechanistic investigations showed that OTUB2 stimulated the Warburg effect and induced the activation of the serine/threonine kinase/mechanistic target of rapamycin kinase (AKT/mTOR) pathway in different NSCLC cells. More importantly, OTUB2 promoted NSCLC progression, which was largely dependent on the direct binding to and deubiquitination of U2AF2, at least in NSCLC cells. U2AF2 expression was also significantly upregulated in primary NSCLC tissues and dramatically associated with metastasis, advanced tumor stages, poor survival, and recurrence. Importantly, a positive correlation between the protein expression of OTUB2 and U2AF2 in NSCLC tissues was found. In vivo experiments indicated that OTUB2 promoted xenograft tumor growth of NSCLC cell. In addition, our results suggest that high expression of OTUB2, U2AF2 and PGK1 is significantly associated with worse prognosis in NSCLC patients. Conclusion: Taken together, the present study provides the first evidence that OTUB2 acts as a pivotal driver in NSCLC tumorigenesis by stabilizing U2AF2 and activating the AKT/mTOR pathway and the Warburg effect. It may serve as a new potential prognostic indicator and therapeutic target in NSCLC.

Genome-wide identification of alternative splicing events that regulate protein transport across the secretory pathway

Alternative splicing (AS) strongly increases proteome diversity and functionality in eukaryotic cells. Protein secretion is a tightly-controlled process, especially in a tissue-specific and differentiation-dependent manner. While previous work has focussed on transcriptional and post-translational regulatory mechanisms, the impact of AS on the secretory pathway remains largely unexplored. Here we integrate a published screen for modulators of protein transport and RNA-Seq analyses to identify over 200 AS events as secretion regulators. We confirm that splicing events along all stages of the secretory pathway regulate the efficiency of membrane trafficking using Morpholinos and CRISPR/Cas9. We furthermore show that these events are highly tissue-specific and adapt the secretory pathway during T-cell activation and adipocyte differentiation. Our data substantially advance the understanding of AS functionality, add a new regulatory layer to a fundamental cell biological process and provide a resource of alternative isoforms that control the secretory pathway.

Alternative Splicing and Protein Diversity: Plants Versus Animals

Plants, unlike animals, exhibit a very high degree of plasticity in their growth and development and employ diverse strategies to cope with the variations during diurnal cycles and stressful conditions. Plants and animals, despite their remarkable morphological and physiological differences, share many basic cellular processes and regulatory mechanisms. Alternative splicing (AS) is one such gene regulatory mechanism that modulates gene expression in multiple ways. It is now well established that AS is prevalent in all multicellular eukaryotes including plants and humans. Emerging evidence indicates that in plants, as in animals, transcription and splicing are coupled. Here, we reviewed recent evidence in support of co-transcriptional splicing in plants and highlighted similarities and differences between plants and humans. An unsettled question in the field of AS is the extent to which splice isoforms contribute to protein diversity. To take a critical look at this question, we presented a comprehensive summary of the current status of research in this area in both plants and humans, discussed limitations with the currently used approaches and suggested improvements to current methods and alternative approaches. We end with a discussion on the potential role of epigenetic modifications and chromatin state in splicing memory in plants primed with stresses.

Sequence-based structural analysis and evaluation of polymorphism in buffalo Nod-like receptor-1 gene

In this study, we have sequence characterized and analyzed the polymorphism in buffalo NOD1 (nucleotide-binding oligomerization domain 1) gene as well as its expression analysis. Full-length sequence analysis of NOD1 revealed this gene in buffalo being conserved with respect to the domain structures, similar to other species. Alternate splice variants having exon3 skipping also identified for the first time in the gene expressed in buffalo-purified peripheral blood mononuclear cells (PBMCs). Phylogenetically ruminant species were found to be clustering together and buffalo displaying maximum similarity with cattle. Sequencing of NOD1 across 12 Indian buffalo breeds identified 23 polymorphic sites within coding region, among which 16 were synonymous and 7 changes found to be non-synonymous. Four SNPs (single nucleotide polymorphisms) of them were genotyped in 393 animals belonging to 12 riverine, swamp and hybrid (riverine × swamp) buffalo populations of diverse phenotypes and utilities, showing variable allelic frequencies. Principal component analysis revealed, riverine and swamp buffaloes being distinctly placed with the distribution of breeds within the group based on the geographical isolation. Further, quantitative real-time PCR detected NOD1 expression in multiple tissues with PBMCs and lungs showing highest expression among the tissues examined. Structural analysis based on the translated amino acid sequence of buffalo NOD1 identified four protein interaction motifs LxxLL important for ligand binding. Molecular interaction analysis of iE-DAP and NOD1-LRR and their complex stability and binding-free energy studies indicated variable binding energies in buffalo and cattle NOD1. Overall, the study reveals unique structural features in buffalo NOD1, important for species-specific ligand interaction.

Advances in Analyzing Virus-Induced Alterations of Host Cell Splicing

Alteration of host cell splicing is a common feature of many viral infections which is underappreciated because of the complexity and technical difficulty of studying alternative splicing (AS) regulation. Recent advances in RNA sequencing technologies revealed that up to several hundreds of host genes can show altered mRNA splicing upon viral infection. The observed changes in AS events can be either a direct consequence of viral manipulation of the host splicing machinery or result indirectly from the virus-induced innate immune response or cellular damage. Analysis at a higher resolution with single-cell RNAseq, and at a higher scale with the integration of multiple omics data sets in a systems biology perspective, will be needed to further comprehend this complex facet of virus-host interactions.

HMGCLL1 is a predictive biomarker for deep molecular response to imatinib therapy in chronic myeloid leukemia

Achieving a deep molecular response (DMR) to tyrosine kinase inhibitor (TKI) therapy for chronic myeloid leukemia (CML) remains challenging and at present, there is no biomarker to predict DMR in this setting. Herein, we report that an HMGCLL1 genetic variant located in 6p12.1 can be used as a predictive genetic biomarker for intrinsic sensitivity to imatinib (IM) therapy. We measured DMR rate according to HMGCLL1 variant in a discovery set of CML patients (n = 201) and successfully replicated it in a validation set (n = 270). We also investigated the functional relevance of HMGCLL1 blockade with respect to response to TKI therapy and showed that small interfering RNA mediated blockade of HMGCLL1 isoform 3 results in significant decrease in viability of BCR-ABL1-positive cells including K562, CML-T1 or BaF3 cell lines with or without ABL1 kinase domain mutations such as T315I mutation. Decreased cell viability was also demonstrated in murine CML stem cells and human hematopoietic progenitor cells. RNA sequencing showed that blockade of HMGCLL1 was associated with G0/G1 arrest and the cell cycle. In summary, the HMGCLL1 gene polymorphism is a novel genetic biomarker for intrinsic sensitivity to IM therapy in CML patients that predicts DMR in this setting.

Integrating molecular networks with genetic variant interpretation for precision medicine

More reliable and cheaper sequencing technologies have revealed the vast mutational landscapes characteristic of many phenotypes. The analysis of such genetic variants has led to successful identification of altered proteins underlying many Mendelian disorders. Nevertheless the simple one‐variant one‐phenotype model valid for many monogenic diseases does not capture the complexity of polygenic traits and disorders. Although experimental and computational approaches have improved detection of functionally deleterious variants and important interactions between gene products, the development of comprehensive models relating genotype and phenotypes remains a challenge in the field of genomic medicine. In this context, a new view of the pathologic state as significant perturbation of the network of interactions between biomolecules is crucial for the identification of biochemical pathways associated with complex phenotypes. Seminal studies in systems biology combined the analysis of genetic variation with protein–protein interaction networks to demonstrate that even as biological systems evolve to be robust to genetic variation, their topologies create disease vulnerabilities. More recent analyses model the impact of genetic variants as changes to the “wiring” of the interactome to better capture heterogeneity in genotype–phenotype relationships. These studies lay the foundation for using networks to predict variant effects at scale using machine‐learning or algorithmic approaches. A wealth of databases and resources for the annotation of genotype–phenotype relationships have been developed to support developments in this area. This overview describes how study of the molecular interactome has generated insights linking the organization of biological systems to disease mechanism, and how this information can enable precision medicine.

This article is categorized under:

Translational, Genomic, and Systems Medicine > Translational Medicine

Biological Mechanisms > Cell Signaling

Models of Systems Properties and Processes > Mechanistic Models

Analytical and Computational Methods > Computational Methods

Protein Language: Post-Translational Modifications Talking to Each Other

Post-translational modifications (PTMs) are at the heart of many cellular signaling events. Apart from a single regulatory PTM, there are also PTMs that function in orchestrated manners. Such PTM crosstalk usually serves as a fine-tuning mechanism to adjust cellular responses to the slightest changes in the environment. While PTM crosstalk has been studied in depth in various species; in plants, this field is just emerging. In this review, we discuss recent studies on crosstalk between three of the most common protein PTMs in plant cells, being phosphorylation, ubiquitination, and sumoylation, and we highlight the diverse underlying mechanisms as well as signaling outputs of such crosstalk.

Identification and Quantification of Murine Mitochondrial Proteoforms Using an Integrated Top-Down and Intact-Mass Strategy

The development of effective strategies for the comprehensive identification and quantification of proteoforms in complex systems is a critical challenge in proteomics. Proteoforms, the specific molecular forms in which proteins are present in biological systems, are the key effectors of biological function. Thus, knowledge of proteoform identities and abundances is essential to unraveling the mechanisms that underlie protein function. We recently reported a strategy that integrates conventional top-down mass spectrometry with intact-mass determinations for enhanced proteoform identifications and the elucidation of proteoform families and applied it to the analysis of yeast cell lysate. In the present work, we extend this strategy to enable quantification of proteoforms, and we examine changes in the abundance of murine mitochondrial proteoforms upon differentiation of mouse myoblasts to myotubes. The integrated top-down and intact-mass strategy provided an increase of ∼37% in the number of identified proteoforms compared to top-down alone, which is in agreement with our previous work in yeast; 1779 unique proteoforms were identified using the integrated strategy compared to 1301 using top-down analysis alone. Quantitative comparison of proteoform differences between the myoblast and myotube cell types showed 129 observed proteoforms exhibiting statistically significant abundance changes (fold change >2 and false discovery rate <5%).

Alternative splicing rewires Hippo signaling pathway in hepatocytes to promote liver regeneration

During liver regeneration, most new hepatocytes arise via self-duplication; yet, the underlying mechanisms that drive hepatocyte proliferation following injury remain poorly defined. By combining high-resolution transcriptome- and polysome-profiling of hepatocytes purified from quiescent and toxin-injured mouse livers, we uncover pervasive alterations in the mRNA translation of metabolic and RNA processing factors, which modulate the protein levels of a set of splicing regulators. Specifically, downregulation of ESRP2 activates a neonatal alternative splicing program that rewires the Hippo signaling pathway in regenerating hepatocytes. We show that production of neonatal splice isoforms attenuates Hippo signaling, enables greater transcriptional activation of downstream target genes, and facilitates liver regeneration. We further demonstrate that ESRP2 deletion in mice causes excessive hepatocyte proliferation upon injury, whereas forced expression of ESRP2 inhibits proliferation by suppressing the expression of neonatal Hippo pathway isoforms. Thus, our findings reveal an ESRP2-Hippo pathway-alternative splicing axis that supports regeneration following chronic liver injury.

Profiles of alternative splicing in colorectal cancer and their clinical significance: A study based on large-scale sequencing data

Background

Alternative splicing (AS), as a potent and pervasive mechanism of transcriptional regulatory, expands the genome's coding capacity and involves in the initiation and progression of cancer. Systematic analysis of alternative splicing in colorectal cancer (CRC) is lacking and greatly needed.

Methods

RNA-Seq data and corresponding clinical information of CRC cohort were downloaded from the TCGA data portal. Then, a java application, known as SpliceSeq, was used to evaluate the RNA splicing patterns and calculate the Percent Spliced In (PSI) value. Differently expressed AS events (DEAS) were identified based on PSI value between paired CRC and adjacent tissues. DEAS and its splicing networks were further analyzed by bioinformatics methods. Kaplan-Meier, Cox proportional regression and unsupervised clustering analysis were used to evaluate the association between DEAS and patients' clinical features.

Results

After strict filtering, a total of 34,334 AS events were identified, among which 421 AS events were found expressed differently. Parent genes of these DEAS play a important role in regulating CRC-related processes such as protein kinase activity (FDR<0.0001), PI3K-Akt signaling pathway (FDR = 0.0024) and p53 signaling pathway (FDR = 0.0143). 37 DEAS events were found to be associated with OS, and 68 DEAS events were found to be associated with DFS. Stratifying patients according to the PSI value of AT in CXCL12 and RI in CSTF3 formed significant Kaplan-Meier curves in both OS and DFS survival analysis. Unsupervised clustering analysis using DEAS revealed four clusters with distinct survival patterns, and associated with consensus molecular subtypes.

Conclusions

Large differences of AS events in CRC appear to exist, and these differences are likely to be important determinants of both prognosis and biological regulation. Our identified CRC-related AS events and uncovered splicing networks are valuable in deciphering the underlying mechanisms of AS in CRC, and provide clues of therapeutic targets to further validations.

Origin of exon skipping-rich transcriptomes in animals driven by evolution of gene architecture

Background

Alternative splicing, particularly through intron retention and exon skipping, is a major layer of pre-translational regulation in eukaryotes. While intron retention is believed to be the most prevalent mode across non-animal eukaryotes, animals have unusually high rates of exon skipping. However, when and how this high prevalence of exon skipping evolved is unknown. Since exon skipping can greatly expand proteomes, answering these questions sheds light on the evolution of higher organismal complexity in metazoans.

Results

We used RNA-seq data to quantify exon skipping and intron retention frequencies across 65 eukaryotic species, with particular focus on early branching animals and unicellular holozoans. We found that only bilaterians have significantly increased their exon skipping frequencies compared to all other eukaryotic groups. Unlike in other eukaryotes, however, exon skipping in nearly all animals, including non-bilaterians, is strongly enriched for frame-preserving sequences, suggesting that exon skipping involvement in proteome expansion predated the increase in frequency. We also identified architectural features consistently associated with higher exon skipping rates within all studied eukaryotic genomes. Remarkably, these architectures became more prevalent during animal evolution, indicating co-evolution between genome architectures and exon skipping frequencies.

Conclusion

We suggest that the increase of exon skipping rates in animals followed a two-step process. First, exon skipping in early animals became enriched for frame-preserving events. Second, bilaterian ancestors dramatically increased their exon skipping frequencies, likely driven by the interplay between a shift in their genome architectures towards more exon definition and recruitment of frame-preserving exon skipping events to functionally diversify their cell-specific proteomes.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1499-9) contains supplementary material, which is available to authorized users.

Characterization of cis-acting elements that control oscillating alternative splicing

Alternative splicing (AS) in response to changing external conditions often requires alterations in the ability of sequence-specific RNA-binding proteins to bind to cis-acting sequences in their target pre-mRNA. While daily oscillations in AS events have been described in several organisms, cis-acting sequences that control time of the day-dependent AS remain largely elusive. Here we define cis-regulatory RNA elements that control body-temperature driven rhythmic AS using the mouse U2af26 gene as a model system. We identify a complex network of cis-regulatory sequences that regulate AS of U2af26, and show that the activity of two enhancer elements is necessary for oscillating AS. A minigene comprising these U2af26 regions recapitulates rhythmic splicing of the endogenous gene, which is controlled through temperature-regulated SR protein phosphorylation. Mutagenesis of the minigene delineates the cis-acting enhancer element for SRSF2 within exon 6 to single nucleotide resolution and reveals that the combined activity of SRSF2 and SRSF7 is required for oscillating U2af26 AS. By combining RNA-Seq with an siRNA screen and individual-nucleotide resolution cross-linking and immunoprecipitation (iCLIP), we identify a complex network of SR proteins that globally controls temperature-dependent rhythmic AS, with the direction of splicing depending on the position of the cis-acting elements. Together, we provide detailed insights into the sequence requirements that allow trans-acting factors to generate daily rhythms in AS.

Decoding a cancer-relevant splicing decision in the RON proto-oncogene using high-throughput mutagenesis

Mutations causing aberrant splicing are frequently implicated in human diseases including cancer. Here, we establish a high-throughput screen of randomly mutated minigenes to decode the cis-regulatory landscape that determines alternative splicing of exon 11 in the proto-oncogene MST1R (RON). Mathematical modelling of splicing kinetics enables us to identify more than 1000 mutations affecting RON exon 11 skipping, which corresponds to the pathological isoform RON∆165. Importantly, the effects correlate with RON alternative splicing in cancer patients bearing the same mutations. Moreover, we highlight heterogeneous nuclear ribonucleoprotein H (HNRNPH) as a key regulator of RON splicing in healthy tissues and cancer. Using iCLIP and synergy analysis, we pinpoint the functionally most relevant HNRNPH binding sites and demonstrate how cooperative HNRNPH binding facilitates a splicing switch of RON exon 11. Our results thereby offer insights into splicing regulation and the impact of mutations on alternative splicing in cancer.

Searching the overlap between network modules with specific betweeness (S2B) and its application to cross-disease analysis

Discovering disease-associated genes (DG) is strategic for understanding pathological mechanisms. DGs form modules in protein interaction networks and diseases with common phenotypes share more DGs or have more closely interacting DGs. This prompted the development of Specific Betweenness (S2B) to find genes associated with two related diseases. S2B prioritizes genes frequently and specifically present in shortest paths linking two disease modules. Top S2B scores identified genes in the overlap of artificial network modules more than 80% of the times, even with incomplete or noisy knowledge. Applied to Amyotrophic Lateral Sclerosis and Spinal Muscular Atrophy, S2B candidates were enriched in biological processes previously associated with motor neuron degeneration. Some S2B candidates closely interacted in network cliques, suggesting common molecular mechanisms for the two diseases. S2B is a valuable tool for DG prediction, bringing new insights into pathological mechanisms. More generally, S2B can be applied to infer the overlap between other types of network modules, such as functional modules or context-specific subnetworks. An R package implementing S2B is publicly available at https://github.com/frpinto/S2B .

Systematic Analysis of Survival-Associated Alternative Splicing Signatures in Gastrointestinal Pan-Adenocarcinomas

BACKGROUND

Gastrointestinal pan-adenocarcinomas, which mainly include adenocarcinomas of the esophagus, stomach, colon, and rectum, place a heavy burden on society owing to their poor prognoses. Since aberrant alternative splicing (AS) are starting to be considered as efficacious signatures for tumor prognosis predicting and therapeutic targets, systematic analysis of AS events is urgent.

METHODS

Prognosis-related AS events were selected by using univariate COX regression analysis. Gene functional enrichment analysis revealed the pathways enriched by prognosis-related AS. Then, prognostic signatures based on AS events were developed for prognosis prediction. Potential mechanism to regulate splicing events by splicing factors was analyzed via Pearson correlation and regulatory networks were constructed.

FINDINGS

A total of 967, 918, 674, and 406 AS events were identified as prognosis-related AS events in esophagus, stomach, colon, and rectum adenocarcinomas, respectively. Survival-associated AS events were distinguishing in the four subtypes of adenocarcinoma. Furthermore, computational algorithm results indicated that perturbation of ribosome and ubiquitin-mediated proteolysis pathways were the potential molecular mechanisms corresponding to inferior prognoses. Most notably, several prognostic signatures based on AS events displayed moderate performance in prognosis predicting. The area under curve values of the time-dependent receiver operating characteristic were 0.961, 0.871, 0.870, and 0.890 in esophagus, stomach, colon, and rectum adenocarcinomas. Survival-associated splicing factors were submitted to construct the AS regulatory network, which could be an underlying mechanism of AS events.

INTERPRETATION

AS may could be ideal indiactors in the prognosis of gastrointestinal pan-adenocarcinomas. Exploring interesting splicing regulatory networks is conducive to solve the puzzles of AS.

Alternative splicing in cardiomyopathy

Alternative splicing is an important mechanism used by the cell to generate greater transcriptomic and proteomic diversity from the genome. In the heart, alternative splicing is increasingly being recognised as an important layer of post-transcriptional gene regulation. Driven by rapidly evolving technologies in next-generation sequencing, alternative splicing has emerged as a crucial process governing complex biological processes during cardiac development and disease. The recent identification of several cardiac splice factors, such as RNA-binding motif protein 20 and 24, not only provided important insight into the mechanisms underlying alternative splicing but also revealed how these splicing factors impact functional properties of the heart. Here, we review our current knowledge of alternative splicing in the heart, with a particular focus on the factors controlling cardiac alternative splicing and their role in cardiomyopathies and subsequent heart failure.

Gene co-expression analysis for functional classification and gene-disease predictions

Gene co-expression networks can be used to associate genes of unknown function with biological processes, to prioritize candidate disease genes or to discern transcriptional regulatory programmes. With recent advances in transcriptomics and next-generation sequencing, co-expression networks constructed from RNA sequencing data also enable the inference of functions and disease associations for non-coding genes and splice variants. Although gene co-expression networks typically do not provide information about causality, emerging methods for differential co-expression analysis are enabling the identification of regulatory genes underlying various phenotypes. Here, we introduce and guide researchers through a (differential) co-expression analysis. We provide an overview of methods and tools used to create and analyse co-expression networks constructed from gene expression data, and we explain how these can be used to identify genes with a regulatory role in disease. Furthermore, we discuss the integration of other data types with co-expression networks and offer future perspectives of co-expression analysis.

Genome-Wide Association Analyses Reveal the Importance of Alternative Splicing in Diversifying Gene Function and Regulating Phenotypic Variation in Maize

Alternative splicing (AS) enhances transcriptome diversity and plays important roles in regulating plant processes. Although widespread natural variation in AS has been observed in plants, how AS is regulated and contribute to phenotypic variation is poorly understood. Here, we report a population-level transcriptome assembly and genome-wide association study to identify splicing quantitative trait loci (sQTLs) in developing maize (Zea mays) kernels from 368 inbred lines. We detected 19,554 unique sQTLs for 6570 genes. Most sQTLs showed small isoform usage changes without involving major isoform switching between genotypes. The sQTL-affected isoforms tend to display distinct protein functions. We demonstrate that nonsense-mediated mRNA decay, microRNA-mediated regulation, and small interfering peptide-mediated peptide interference are frequently involved in sQTL regulation. The natural variation in AS and overall mRNA level appears to be independently regulated with different cis-sequences preferentially used. We identified 214 putative trans-acting splicing regulators, among which ZmGRP1, encoding an hnRNP-like glycine-rich RNA binding protein, regulates the largest trans-cluster. Knockout of ZmGRP1 by CRISPR/Cas9 altered splicing of numerous downstream genes. We found that 739 sQTLs colocalized with previous marker-trait associations, most of which occurred without changes in overall mRNA level. Our findings uncover the importance of AS in diversifying gene function and regulating phenotypic variation.

Gene Regulatory Network Perturbation by Genetic and Epigenetic Variation

Gene regulatory networks underlie biological function and cellular physiology. Alternative splicing (AS) is a fundamental step in gene regulatory networks and plays a key role in development and disease. In addition to the identification of aberrant AS events, an increasing number of studies are focusing on molecular determinants of AS, including genetic and epigenetic regulators. We review here recent efforts to identify various deregulated AS events as well as their molecular determinants that alter biological functions, and discuss clinical features of AS and their druggable potential.

Nuclear PTEN safeguards pre-mRNA splicing to link Golgi apparatus for its tumor suppressive role

Dysregulation of pre-mRNA alternative splicing (AS) is closely associated with cancers. However, the relationships between the AS and classic oncogenes/tumor suppressors are largely unknown. Here we show that the deletion of tumor suppressor PTEN alters pre-mRNA splicing in a phosphatase-independent manner, and identify 262 PTEN-regulated AS events in 293T cells by RNA sequencing, which are associated with significant worse outcome of cancer patients. Based on these findings, we report that nuclear PTEN interacts with the splicing machinery, spliceosome, to regulate its assembly and pre-mRNA splicing. We also identify a new exon 2b in GOLGA2 transcript and the exon exclusion contributes to PTEN knockdown-induced tumorigenesis by promoting dramatic Golgi extension and secretion, and PTEN depletion significantly sensitizes cancer cells to secretion inhibitors brefeldin A and golgicide A. Our results suggest that Golgi secretion inhibitors alone or in combination with PI3K/Akt kinase inhibitors may be therapeutically useful for PTEN-deficient cancers.

How alternative splicing affects membrane-trafficking dynamics

The cell biology field has outstanding working knowledge of the fundamentals of membrane-trafficking pathways, which are of critical importance in health and disease. Current challenges include understanding how trafficking pathways are fine-tuned for specialized tissue functions in vivo and during development. In parallel, the ENCODE project and numerous genetic studies have revealed that alternative splicing regulates gene expression in tissues and throughout development at a post-transcriptional level. This Review summarizes recent discoveries demonstrating that alternative splicing affects tissue specialization and membrane-trafficking proteins during development, and examines how this regulation is altered in human disease. We first discuss how alternative splicing of clathrin, SNAREs and BAR-domain proteins influences endocytosis, secretion and membrane dynamics, respectively. We then focus on the role of RNA-binding proteins in the regulation of splicing of membrane-trafficking proteins in health and disease. Overall, our aim is to comprehensively summarize how trafficking is molecularly influenced by alternative splicing and identify future directions centered on its physiological relevance.

Expansion of DUB functionality generated by alternative isoforms - USP35, a case study

Protein ubiquitylation is a dynamic post-translational modification that can be reversed by deubiquitylating enzymes (DUBs). It is unclear how the small number (∼100) of DUBs present in mammalian cells regulate the thousands of different ubiquitylation events. Here, we analysed annotated transcripts of human DUBs and found ∼300 ribosome-associated transcripts annotated as protein coding, which thus increases the total number of DUBs. By using USP35, a poorly studied DUB, as a case study, we provide evidence that alternative isoforms contribute to the functional expansion of DUBs. We show that there are two different USP35 isoforms that localise to different intracellular compartments and have distinct functions. Our results reveal that isoform 1 is an anti-apoptotic factor that inhibits staurosporine- and TNF-related apoptosis-inducing ligand (TRAIL; also known as TNFSF10)-induced apoptosis. In contrast, USP35 isoform 2 is an integral membrane protein of the endoplasmic reticulum (ER) that is also present at lipid droplets. Manipulations of isoform 2 levels cause rapid ER stress, likely through deregulation of lipid homeostasis, and lead to cell death. Our work highlights how alternative isoforms provide functional expansion of DUBs and sets directions for future research.This article has an associated First Person interview with the first author of the paper.

Unexplored therapeutic opportunities in the human genome

A large proportion of biomedical research and the development of therapeutics is focused on a small fraction of the human genome. In a strategic effort to map the knowledge gaps around proteins encoded by the human genome and to promote the exploration of currently understudied, but potentially druggable, proteins, the US National Institutes of Health launched the Illuminating the Druggable Genome (IDG) initiative in 2014. In this article, we discuss how the systematic collection and processing of a wide array of genomic, proteomic, chemical and disease-related resource data by the IDG Knowledge Management Center have enabled the development of evidence-based criteria for tracking the target development level (TDL) of human proteins, which indicates a substantial knowledge deficit for approximately one out of three proteins in the human proteome. We then present spotlights on the TDL categories as well as key drug target classes, including G protein-coupled receptors, protein kinases and ion channels, which illustrate the nature of the unexplored opportunities for biomedical research and therapeutic development.

Comprehensive Analyses of Tissue-Specific Networks with Implications to Psychiatric Diseases

Recent advances in genome sequencing and "omics" technologies are opening new opportunities for improving diagnosis and treatment of human diseases. The precision medicine initiative in particular aims at developing individualized treatment options that take into account individual variability in genes and environment of each person. Systems biology approaches that group genes, transcripts and proteins into functionally meaningful networks will play crucial role in the future of personalized medicine. They will allow comparison of healthy and disease-affected tissues and organs from the same individual, as well as between healthy and disease-afflicted individuals. However, the field faces a multitude of challenges ranging from data integration to statistical and combinatorial issues in data analyses. This chapter describes computational approaches developed by us and the others to tackle challenges in tissue-specific network analyses, with the main focus on psychiatric diseases.

Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression

Major depressive disorder (MDD) is a common illness accompanied by considerable morbidity, mortality, costs, and heightened risk of suicide. We conducted a genome-wide association meta-analysis based in 135,458 cases and 344,901 controls and identified 44 independent and significant loci. The genetic findings were associated with clinical features of major depression and implicated brain regions exhibiting anatomical differences in cases. Targets of antidepressant medications and genes involved in gene splicing were enriched for smaller association signal. We found important relationships of genetic risk for major depression with educational attainment, body mass, and schizophrenia: lower educational attainment and higher body mass were putatively causal, whereas major depression and schizophrenia reflected a partly shared biological etiology. All humans carry lesser or greater numbers of genetic risk factors for major depression. These findings help refine the basis of major depression and imply that a continuous measure of risk underlies the clinical phenotype.

Network Analysis of UBE3A/E6AP-Associated Proteins Provides Connections to Several Distinct Cellular Processes

Perturbations in activity and dosage of the UBE3A ubiquitin-ligase have been linked to Angelman syndrome and autism spectrum disorders. UBE3A was initially identified as the cellular protein hijacked by the human papillomavirus E6 protein to mediate the ubiquitylation of p53, a function critical to the oncogenic potential of these viruses. Although a number of substrates have been identified, the normal cellular functions and pathways affected by UBE3A are largely unknown. Previously, we showed that UBE3A associates with HERC2, NEURL4, and MAPK6/ERK3 in a high-molecular-weight complex of unknown function that we refer to as the HUN complex (HERC2, UBE3A, and NEURL4). In this study, the combination of two complementary proteomic approaches with a rigorous network analysis revealed cellular functions and pathways in which UBE3A and the HUN complex are involved. In addition to finding new UBE3A-associated proteins, such as MCM6, SUGT1, EIF3C, and ASPP2, network analysis revealed that UBE3A-associated proteins are connected to several fundamental cellular processes including translation, DNA replication, intracellular trafficking, and centrosome regulation. Our analysis suggests that UBE3A could be involved in the control and/or integration of these cellular processes, in some cases as a component of the HUN complex, and also provides evidence for crosstalk between the HUN complex and CAMKII interaction networks. This study contributes to a deeper understanding of the cellular functions of UBE3A and its potential role in pathways that may be affected in Angelman syndrome, UBE3A-associated autism spectrum disorders, and human papillomavirus-associated cancers.

Full-length mRNA sequencing uncovers a widespread coupling between transcription initiation and mRNA processing

Background

The multifaceted control of gene expression requires tight coordination of regulatory mechanisms at transcriptional and post-transcriptional level. Here, we studied the interdependence of transcription initiation, splicing and polyadenylation events on single mRNA molecules by full-length mRNA sequencing.

Results

In MCF-7 breast cancer cells, we find 2700 genes with interdependent alternative transcription initiation, splicing and polyadenylation events, both in proximal and distant parts of mRNA molecules, including examples of coupling between transcription start sites and polyadenylation sites. The analysis of three human primary tissues (brain, heart and liver) reveals similar patterns of interdependency between transcription initiation and mRNA processing events. We predict thousands of novel open reading frames from full-length mRNA sequences and obtained evidence for their translation by shotgun proteomics. The mapping database rescues 358 previously unassigned peptides and improves the assignment of others. By recognizing sample-specific amino-acid changes and novel splicing patterns, full-length mRNA sequencing improves proteogenomics analysis of MCF-7 cells.

Conclusions

Our findings demonstrate that our understanding of transcriptome complexity is far from complete and provides a basis to reveal largely unresolved mechanisms that coordinate transcription initiation and mRNA processing.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1418-0) contains supplementary material, which is available to authorized users.

Species-conserved SYNGAP1 phenotypes associated with neurodevelopmental disorders

SYNGAP1 loss-of-function variants are causally associated with intellectual disability, severe epilepsy, autism spectrum disorder and schizophrenia. While there are hundreds of genetic risk factors for neurodevelopmental disorders (NDDs), this gene is somewhat unique because of the frequency and penetrance of loss-of-function variants found in patients combined with the range of brain disorders associated with SYNGAP1 pathogenicity. These clinical findings indicate that SYNGAP1 regulates fundamental neurodevelopmental processes that are necessary for brain development. Here, we describe four phenotypic domains that are controlled by Syngap1 expression across vertebrate species. Two domains, the maturation of cognitive functions and maintenance of excitatory-inhibitory balance, are defined exclusively through a review of the current literature. Two additional domains are defined by integrating the current literature with new data indicating that SYNGAP1/Syngap1 regulates innate survival behaviors and brain structure. These four phenotypic domains are commonly disrupted in NDDs, suggesting that a deeper understanding of developmental Syngap1 functions will be generalizable to other NDDs of known or unknown etiology. Therefore, we discuss the known molecular and cellular functions of Syngap1 and consider how these functions may contribute to the emergence of disease-relevant phenotypes. Finally, we identify major unexplored areas of Syngap1 neurobiology and discuss how a deeper understanding of this gene may uncover general principles of NDD pathobiology.

ProForma: A Standard Proteoform Notation

The Consortium for Top-Down Proteomics (CTDP) proposes a standardized notation, ProForma, for writing the sequence of fully characterized proteoforms. ProForma provides a means to communicate any proteoform by writing the amino acid sequence using standard one-letter notation and specifying modifications or unidentified mass shifts within brackets following certain amino acids. The notation is unambiguous, human-readable, and can easily be parsed and written by bioinformatic tools. This system uses seven rules and supports a wide range of possible use cases, ensuring compatibility and reproducibility of proteoform annotations. Standardizing proteoform sequences will simplify storage, comparison, and reanalysis of proteomic studies, and the Consortium welcomes input and contributions from the research community on the continued design and maintenance of this standard.

Alternative transcription of a shorter, non-anti-angiogenic thrombospondin-2 variant in cancer-associated blood vessels

Thrombospondin-2 (TSP2) is an anti-angiogenic matricellular protein that inhibits tumor growth and angiogenesis. Tumor-associated blood vascular endothelial cells (BECs) were isolated from human invasive bladder cancers and from matched normal bladder tissue by immuno-laser capture microdissection. Exon expression profiling analyses revealed a particularly high expression of a short TSP2 transcript containing only the last 9 (3′) exons of the full-length TSP2 transcript. Using 5′ and 3′ RACE (rapid amplification of cDNA ends) and Sanger sequencing, we confirmed the existence of the shorter transcript of TSP2 (sTSP2) and determined its sequence which completely lacked the anti-angiogenic thrombospondin type 1 repeats domain. The largest open reading frame predicted within the transcript comprises 209 amino acids and matches almost completely the C-terminal lectin domain of full-length TSP2. We produced recombinant sTSP2 and found that unlike the full-length TSP2, sTSP2 did not inhibit vascular endothelial growth factor-A-induced proliferation of cultured human BECs, but in contrast when combined with TSP2 blocked the inhibitory effects of TSP2 on BEC proliferation. In vivo studies with stably transfected A431 squamous cell carcinoma cells revealed that full-length TSP2, but not sTSP2, inhibited tumor growth and angiogenesis. This study reveals that the transcriptional program of tumor stromal cells can change to transcribe a new version of an endogenous angiogenesis inhibitor that has lost its anti-angiogenic activity.

Identification of new PDE9A isoforms and how their expression and subcellular compartmentalization in the brain change across the life span

3',5'-Cyclic nucleotide phosphodiesterases (PDEs) degrade 3',5' cyclic adenonosine monophosphate (cAMP) and 3',5' cyclic guanosine monophosphate (cGMP), with PDE9A having the highest affinity for cGMP. We show PDE9A6 and 3 novel PDE9 isoforms (PDE9X-100, PDE9X-120, and PDE9X-175) are reliably detected in the brain and lung of mice, whereas PDE9A2 and other isoforms are found elsewhere. PDE9A localizes to the membrane in all organs except the bladder, where it is cytosolic. Brain additionally shows PDE9 in the nuclear fraction. PDE9A mRNA expression/localization dramatically changes across neurodevelopment in a manner that is strikingly consistent between mice and humans (i.e., decreased expression in the hippocampus and cortex and inverted-U in the cerebellum). Study of the 4 PDE9 isoforms in the mouse brain from postnatal day 7 through 24 months similarly identifies dramatic effects of age on expression and subcellular compartmentalization that are isoform specific and brain region specific. Finally, PDE9A mRNA is elevated in the aged human hippocampus with dementia when there is a history of traumatic brain injury. Thus, brain PDE9 is localized to preferentially regulate nuclear- and membrane-proximal pools of cGMP, and its function likely changes across the life span.

Increased proteomic complexity in Drosophila hybrids during development

Cellular proteomes are thought to be optimized for function, leaving no room for proteome plasticity and, thus, evolution. However, hybrid animals that result from a viable cross of two different species harbor hybrid proteomes of unknown complexity. We charted the hybrid proteome of a viable cross between Drosophila melanogaster females and Drosophila simulans males with bottom-up proteomics. Developing hybrids harbored 20% novel proteins in addition to proteins that were also present in either parental species. In contrast, adult hybrids and developmentally failing embryos of the reciprocal cross showed less additional proteins (5 and 6%, respectively). High levels of heat shock proteins, proteasome-associated proteins, and proteasomal subunits indicated that proteostasis sustains the expanded complexity of the proteome in developing hybrids. We conclude that increased proteostasis gives way to proteomic plasticity and thus opens up additional space for rapid phenotypic variation during embryonic development.

Alternative start and termination sites of transcription drive most transcript isoform differences across human tissues

Most human genes generate multiple transcript isoforms. The differential expression of these isoforms can help specify cell types. Diverse transcript isoforms arise from the use of alternative transcription start sites, polyadenylation sites and splice sites; however, the relative contribution of these processes to isoform diversity in normal human physiology is unclear. To address this question, we investigated cell type-dependent differences in exon usage of over 18 000 protein-coding genes in 23 cell types from 798 samples of the Genotype-Tissue Expression Project. We found that about half of the expressed genes displayed tissue-dependent transcript isoforms. Alternative transcription start and termination sites, rather than alternative splicing, accounted for the majority of tissue-dependent exon usage. We confirmed the widespread tissue-dependent use of alternative transcription start sites in a second, independent dataset, Cap Analysis of Gene Expression data from the FANTOM consortium. Moreover, our results indicate that most tissue-dependent splicing involves untranslated exons and therefore may not increase proteome complexity. Thus, alternative transcription start and termination sites are the principal drivers of transcript isoform diversity across tissues, and may underlie the majority of cell type specific proteomes and functions.

Ovule identity mediated by pre-mRNA processing in Arabidopsis

Ovules are fundamental for plant reproduction and crop yield as they are the precursors of seeds. Therefore, ovule specification is a critical developmental program. In Arabidopsis thaliana, ovule identity is redundantly conferred by the homeotic D-class genes SHATTERPROOF1 (SHP1), SHP2 and SEEDSTICK (STK), phylogenetically related to the MADS-domain regulatory gene AGAMOUS (AG), essential in floral organ specification. Previous studies have shown that the HUA-PEP activity, comprised of a suite of RNA-binding protein (RBP) encoding genes, regulates AG pre-mRNA processing and thus flower patterning and organ identity. Here, we report that the HUA-PEP activity additionally governs ovule morphogenesis. Accordingly, in severe hua-pep backgrounds ovules transform into flower organ-like structures. These homeotic transformations are most likely due to the dramatic reduction in SHP1, SHP2 and STK activity. Our molecular and genome-wide profiling strategies revealed the accumulation of prematurely terminated transcripts of D-class genes in hua-pep mutants and reduced amounts of their respective functional messengers, which points to pre-mRNA processing misregulation as the origin of the ovule developmental defects in such backgrounds. RNA processing and transcription are coordinated by the RNA polymerase II (RNAPII) carboxyl-terminal domain (CTD). Our results show that HUA-PEP activity members can interact with the CTD regulator C-TERMINAL DOMAIN PHOSPHATASE-LIKE1 (CPL1), supporting a co-transcriptional mode of action for the HUA-PEP activity. Our findings expand the portfolio of reproductive developmental programs in which HUA-PEP activity participates, and further substantiates the importance of RNA regulatory mechanisms (pre-mRNA co-transcriptional regulation) for correct gene expression during plant morphogenesis.

The comparison of alternative splicing among the multiple tissues in cucumber

BACKGROUND

Alternative splicing (AS) is an important post-transcriptional process. It has been suggested that most AS events are subject to tissue-specific regulation. However, the global dynamics of AS in different tissues are poorly explored.

RESULTS

To analyse global changes in AS in multiple tissues, we identified the AS events and constructed a comprehensive catalogue of AS events within each tissue based on the genome-wide RNA-seq reads from ten tissues in cucumber. First, we found that 58% of the multi-exon genes underwent AS. We further obtained 565 genes with significantly more AS events compared with random genes. These genes were found significant enrichment in biological processes related to the regulation of actin filament length. Second, significantly different AS event profiles among ten tissues were found. The tissues with the same origin of development are more likely to have a relatively similar AS profile. Moreover, 7370 genes showed tissue-specific AS events and were highly enriched in biological processes related to the positive regulation of cellular component organization. Root-specificity AS genes were related to the cellular response to DNA damage stimulus. Third, the genes with different intron retention (IR) patterns among the ten tissues showed significant difference in GC percentages of the retained intron, and the number of exons and FPKM of the major transcripts.

CONCLUSIONS

Our study provided a comprehensive view of AS in multiple tissues. We revealed novel insights into the patterns of AS in multiple tissues and the tissue-specific AS in cucumber.

Alternative RNA splicing in the endothelium mediated in part by Rbfox2 regulates the arterial response to low flow

Low and disturbed blood flow drives the progression of arterial diseases including atherosclerosis and aneurysms. The endothelial response to flow and its interactions with recruited platelets and leukocytes determine disease progression. Here, we report widespread changes in alternative splicing of pre-mRNA in the flow-activated murine arterial endothelium in vivo. Alternative splicing was suppressed by depletion of platelets and macrophages recruited to the arterial endothelium under low and disturbed flow. Binding motifs for the Rbfox-family are enriched adjacent to many of the regulated exons. Endothelial deletion of Rbfox2, the only family member expressed in arterial endothelium, suppresses a subset of the changes in transcription and RNA splicing induced by low flow. Our data reveal an alternative splicing program activated by Rbfox2 in the endothelium on recruitment of platelets and macrophages and demonstrate its relevance in transcriptional responses during flow-driven vascular inflammation.

Expanding Proteoform Identifications in Top-Down Proteomic Analyses by Constructing Proteoform Families

In top-down proteomics, intact proteins are analyzed by tandem mass spectrometry and proteoforms, which are defined forms of a protein with specific sequences of amino acids and localized post-translational modifications, are identified using precursor mass and fragmentation data. Many proteoforms that are detected in the precursor scan (MS1) are not selected for fragmentation by the instrument and therefore remain unidentified in typical top-down proteomic workflows. Our laboratory has developed the open source software program Proteoform Suite to analyze MS1-only intact proteoform data. Here, we have adapted it to provide identifications of proteoform masses in precursor MS1 spectra of top-down data, supplementing the top-down identifications obtained using the MS2 fragmentation data. Proteoform Suite performs mass calibration using high-scoring top-down identifications and identifies additional proteoforms using calibrated, accurate intact masses. Proteoform families, the set of proteoforms from a given gene, are constructed and visualized from proteoforms identified by both top-down and intact-mass analyses. Using this strategy, we constructed proteoform families and identified 1861 proteoforms in yeast lysate, yielding an approximately 40% increase over the original 1291 proteoform identifications observed using traditional top-down analysis alone.

Proteoform Suite: Software for Constructing, Quantifying, and Visualizing Proteoform Families

We present an open-source, interactive program named Proteoform Suite that uses proteoform mass and intensity measurements from complex biological samples to identify and quantify proteoforms. It constructs families of proteoforms derived from the same gene, assesses proteoform function using gene ontology (GO) analysis, and enables visualization of quantified proteoform families and their changes. It is applied here to reveal systemic proteoform variations in the yeast response to salt stress.

A High-Resolution Genome-Wide CRISPR/Cas9 Viability Screen Reveals Structural Features and Contextual Diversity of the Human Cell-Essential Proteome

To interrogate genes essential for cell growth, proliferation and survival in human cells, we carried out a genome-wide clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9 screen in a B-cell lymphoma line using a custom extended-knockout (EKO) library of 278,754 single-guide RNAs (sgRNAs) that targeted 19,084 RefSeq genes, 20,852 alternatively spliced exons, and 3,872 hypothetical genes. A new statistical analysis tool called robust analytics and normalization for knockout screens (RANKS) identified 2,280 essential genes, 234 of which were unique. Individual essential genes were validated experimentally and linked to ribosome biogenesis and stress responses. Essential genes exhibited a bimodal distribution across 10 different cell lines, consistent with a continuous variation in essentiality as a function of cell type. Genes essential in more lines had more severe fitness defects and encoded the evolutionarily conserved structural cores of protein complexes, whereas genes essential in fewer lines formed context-specific modules and encoded subunits at the periphery of essential complexes. The essentiality of individual protein residues across the proteome correlated with evolutionary conservation, structural burial, modular domains, and protein interaction interfaces. Many alternatively spliced exons in essential genes were dispensable and were enriched for disordered regions. Fitness defects were observed for 44 newly evolved hypothetical reading frames. These results illuminate the contextual nature and evolution of essential gene functions in human cells.

Prediction of Protein-Protein Interactions

The authors provide an overview of physical protein-protein interaction prediction, covering the main strategies for predicting interactions, approaches for assessing predictions, and online resources for accessing predictions. This unit focuses on the main advancements in each of these areas over the last decade. The methods and resources that are presented here are not an exhaustive set, but characterize the current state of the field-highlighting key challenges and achievements. © 2017 by John Wiley & Sons, Inc.

Living Organisms Author Their Read-Write Genomes in Evolution

Evolutionary variations generating phenotypic adaptations and novel taxa resulted from complex cellular activities altering genome content and expression: (i) Symbiogenetic cell mergers producing the mitochondrion-bearing ancestor of eukaryotes and chloroplast-bearing ancestors of photosynthetic eukaryotes; (ii) interspecific hybridizations and genome doublings generating new species and adaptive radiations of higher plants and animals; and, (iii) interspecific horizontal DNA transfer encoding virtually all of the cellular functions between organisms and their viruses in all domains of life. Consequently, assuming that evolutionary processes occur in isolated genomes of individual species has become an unrealistic abstraction. Adaptive variations also involved natural genetic engineering of mobile DNA elements to rewire regulatory networks. In the most highly evolved organisms, biological complexity scales with "non-coding" DNA content more closely than with protein-coding capacity. Coincidentally, we have learned how so-called "non-coding" RNAs that are rich in repetitive mobile DNA sequences are key regulators of complex phenotypes. Both biotic and abiotic ecological challenges serve as triggers for episodes of elevated genome change. The intersections of cell activities, biosphere interactions, horizontal DNA transfers, and non-random Read-Write genome modifications by natural genetic engineering provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

Extended Linkers Improve the Detection of Protein-protein Interactions (PPIs) by Dihydrofolate Reductase Protein-fragment Complementation Assay (DHFR PCA) in Living Cells

Understanding the function of cellular systems requires describing how proteins assemble with each other into transient and stable complexes and to determine their spatial relationships. Among the tools available to perform these analyses on a large scale is Protein-fragment Complementation Assay based on the dihydrofolate reductase (DHFR PCA). Here we test how longer linkers between the fusion proteins and the reporter fragments affect the performance of this assay. We investigate the architecture of the RNA polymerases, the proteasome and the conserved oligomeric Golgi (COG) complexes in living cells and performed large-scale screens with these extended linkers. We show that longer linkers significantly improve the detection of protein-protein interactions and allow to measure interactions further in space than the standard ones. We identify new interactions, for instance between the retromer complex and proteins related to autophagy and endocytosis. Longer linkers thus contribute an enhanced additional tool to the existing toolsets for the detection and measurements of protein-protein interactions and protein proximity in living cells.

Microfluidic isoform sequencing shows widespread splicing coordination in the human transcriptome

Understanding transcriptome complexity is crucial for understanding human biology and disease. Technologies such as Synthetic long-read RNA sequencing (SLR-RNA-seq) delivered 5 million isoforms and allowed assessing splicing coordination. Pacific Biosciences and Oxford Nanopore increase throughput also but require high input amounts or amplification. Our new droplet-based method, sparse isoform sequencing (spISO-seq), sequences 100k–200k partitions of 10–200 molecules at a time, enabling analysis of 10–100 million RNA molecules. SpISO-seq requires less than 1 ng of input cDNA, limiting or removing the need for prior amplification with its associated biases. Adjusting the number of reads devoted to each molecule reduces sequencing lanes and cost, with little loss in detection power. The increased number of molecules expands our understanding of isoform complexity. In addition to confirming our previously published cases of splicing coordination (e.g., BIN1), the greater depth reveals many new cases, such as MAPT. Coordination of internal exons is found to be extensive among protein coding genes: 23.5%–59.3% (95% confidence interval) of highly expressed genes with distant alternative exons exhibit coordination, showcasing the need for long-read transcriptomics. However, coordination is less frequent for noncoding sequences, suggesting a larger role of splicing coordination in shaping proteins. Groups of genes with coordination are involved in protein–protein interactions with each other, raising the possibility that coordination facilitates complex formation and/or function. We also find new splicing coordination types, involving initial and terminal exons. Our results provide a more comprehensive understanding of the human transcriptome and a general, cost-effective method to analyze it.