Tissue-specific alternative splicing remodels protein-protein interaction networks

Origin of spliceosomal introns and alternative splicing

In this work we review the current knowledge on the prehistory, origins, and evolution of spliceosomal introns. First, we briefly outline the major features of the different types of introns, with particular emphasis on the nonspliceosomal self-splicing group II introns, which are widely thought to be the ancestors of spliceosomal introns. Next, we discuss the main scenarios proposed for the origin and proliferation of spliceosomal introns, an event intimately linked to eukaryogenesis. We then summarize the evidence that suggests that the last eukaryotic common ancestor (LECA) had remarkably high intron densities and many associated characteristics resembling modern intron-rich genomes. From this intron-rich LECA, the different eukaryotic lineages have taken very distinct evolutionary paths leading to profoundly diverged modern genome structures. Finally, we discuss the origins of alternative splicing and the qualitative differences in alternative splicing forms and functions across lineages.

Lineage-specific splicing of a brain-enriched alternative exon promotes glioblastoma progression

Tissue-specific alternative splicing is critical for the emergence of tissue identity during development, yet the role of this process in malignant transformation is undefined. Tissue-specific splicing involves evolutionarily conserved, alternative exons that represent only a minority of the total alternative exons identified. Many of these conserved exons have functional features that influence signaling pathways to profound biological effect. Here, we determined that lineage-specific splicing of a brain-enriched cassette exon in the membrane-binding tumor suppressor annexin A7 (ANXA7) diminishes endosomal targeting of the EGFR oncoprotein, consequently enhancing EGFR signaling during brain tumor progression. ANXA7 exon splicing was mediated by the ribonucleoprotein PTBP1, which is normally repressed during neuronal development. PTBP1 was highly expressed in glioblastomas due to loss of a brain-enriched microRNA (miR-124) and to PTBP1 amplification. The alternative ANXA7 splicing trait was present in precursor cells, suggesting that glioblastoma cells inherit the trait from a potential tumor-initiating ancestor and that these cells exploit this trait through accumulation of mutations that enhance EGFR signaling. Our data illustrate that lineage-specific splicing of a tissue-regulated alternative exon in a constituent of an oncogenic pathway eliminates tumor suppressor functions and promotes glioblastoma progression. This paradigm may offer a general model as to how tissue-specific regulatory mechanisms can reprogram normal developmental processes into oncogenic ones.

Distinct functions of alternatively spliced isoforms encoded by zebrafish mef2ca and mef2cb

In mammals, an array of MEF2C proteins is generated by alternative splicing (AS), yet specific functions have not been ascribed to each isoform. Teleost fish possess two MEF2C paralogues, mef2ca and mef2cb. In zebrafish, the Mef2cs function to promote cardiomyogenic differentiation and myofibrillogenesis in nascent skeletal myofibers. We found that zebrafish mef2ca and mef2cb are alternatively spliced in the coding exons 4–6 region and these splice variants differ in their biological activity. Of the two, mef2ca is more abundantly expressed in developing skeletal muscle, its activity is tuned through zebrafish development by AS. By 24 hpf, we found the prevalent expression of the highly active full length protein in differentiated muscle in the somites. The splicing isoform of mef2ca that lacks exon 5 (mef2ca 4–6), encodes a protein that has 50% lower transcriptional activity, and is found mainly earlier in development, before muscle differentiation. mef2ca transcripts including exon 5 (mef2ca 4–5–6) are present early in the embryo. Over-expression of this isoform alters the expression of genes involved in early dorso-ventral patterning of the embryo such as chordin, nodal related 1 and goosecoid, and induces severe developmental defects. AS of mef2cb generates a long splicing isoform in the exon 5 region (Mef2cbL) that predominates during somitogenesis. Mef2cbL contains an evolutionarily conserved domain derived from exonization of a fragment of intron 5, which confers the ability to induce ectopic muscle in mesoderm upon over-expression of the protein. Taken together, the data show that AS is a significant regulator of Mef2c activity.

•

mef2ca and mef2cb gene products are alternatively spliced in zebrafish.

•

Inclusion of exon 5 in mef2ca transcripts is regulated during zebrafish development.

•

Exon 5 confers on Mef2ca the ability to activate early patterning genes.

•

Mef2cb includes an extra octapeptide encoded by a region of intron 5.

•

Inclusion of the extra-octapeptide confers on Mef2cb pro-myogenic activity.

Activation-induced tumor necrosis factor receptor-associated factor 3 (Traf3) alternative splicing controls the noncanonical nuclear factor κB pathway and chemokine expression in human T cells

The noncanonical nuclear factor κB (ncNFκB) pathway regulates the expression of chemokines required for secondary lymphoid organ formation and thus plays a pivotal role in adaptive immunity. Whereas ncNFκB signaling has been well described in stromal cells and B cells, its role and regulation in T cells remain largely unexplored. ncNFκB activity critically depends on the upstream NFκB-inducing kinase (NIK). NIK expression is negatively regulated by the full-length isoform of TNF receptor-associated factor 3 (Traf3) as formation of a NIK-Traf3-Traf2 complex targets NIK for degradation. Here we show that T cell-specific and activation-dependent alternative splicing generates a Traf3 isoform lacking exon 8 (Traf3DE8) that, in contrast to the full-length protein, activates ncNFκB signaling. Traf3DE8 disrupts the NIK-Traf3-Traf2 complex and allows accumulation of NIK to initiate ncNFκB signaling in activated T cells. ncNFκB activity results in expression of several chemokines, among them B cell chemoattractant (CxCL13), both in a model T cell line and in primary human CD4(+) T cells. Because CxCL13 plays an important role in B cell migration and activation, our data suggest an involvement and provide a mechanistic basis for Traf3 alternative splicing and ncNFκB activation in contributing to T cell-dependent adaptive immunity.

Evolution at protein ends: major contribution of alternative transcription initiation and termination to the transcriptome and proteome diversity in mammals

Alternative splicing (AS), alternative transcription initiation (ATI) and alternative transcription termination (ATT) create the extraordinary complexity of transcriptomes and make key contributions to the structural and functional diversity of mammalian proteomes. Analysis of mammalian genomic and transcriptomic data shows that contrary to the traditional view, the joint contribution of ATI and ATT to the transcriptome and proteome diversity is quantitatively greater than the contribution of AS. Although the mean numbers of protein-coding constitutive and alternative nucleotides in gene loci are nearly identical, their distribution along the transcripts is highly non-uniform. On average, coding exons in the variable 5' and 3' transcript ends that are created by ATI and ATT contain approximately four times more alternative nucleotides than core protein-coding regions that diversify exclusively via AS. Short upstream exons that encompass alternative 5'-untranslated regions and N-termini of proteins evolve under strong nucleotide-level selection whereas in 3'-terminal exons that encode protein C-termini, protein-level selection is significantly stronger. The groups of genes that are subject to ATI and ATT show major differences in biological roles, expression and selection patterns.

Human proteins with target sites of multiple post-translational modification types are more prone to be involved in disease

Many proteins can be modified by multiple types of post-translational modifications (Mtp-proteins). Although some post-translational modifications (PTMs) have recently been found to be associated with life-threatening diseases like cancers and neurodegenerative disorders, the underlying mechanisms remain enigmatic to date. In this study, we examined the relationship of human Mtp-proteins and disease and systematically characterized features of these proteins. Our results indicated that Mtp-proteins are significantly more inclined to participate in disease than proteins carrying no known PTM sites. Mtp-proteins were found significantly enriched in protein complexes, having more protein partners and preferred to act as hubs/superhubs in protein-protein interaction (PPI) networks. They possess a distinct functional focus, such as chromatin assembly or disassembly, and reside in biased, multiple subcellular localizations. Moreover, most Mtp-proteins harbor more intrinsically disordered regions than the others. Mtp-proteins carrying PTM types biased toward locating in the ordered regions were mainly related to protein-DNA complex assembly. Examination of the energetic effects of PTMs on the stability of PPI revealed that only a small fraction of single PTM events influence the binding energy of >2 kcal/mol, whereas the binding energy can change dramatically by combinations of multiple PTM types. Our work not only expands the understanding of Mtp-proteins but also discloses the potential ability of Mtp-proteins to act as key elements in disease development.

Mutations in BIN1 associated with centronuclear myopathy disrupt membrane remodeling by affecting protein density and oligomerization

The regulation of membrane shapes is central to many cellular phenomena. Bin/Amphiphysin/Rvs (BAR) domain-containing proteins are key players for membrane remodeling during endocytosis, cell migration, and endosomal sorting. BIN1, which contains an N-BAR domain, is assumed to be essential for biogenesis of plasma membrane invaginations (T-tubules) in muscle tissues. Three mutations, K35N, D151N and R154Q, have been discovered so far in the BAR domain of BIN1 in patients with centronuclear myopathy (CNM), where impaired organization of T-tubules has been reported. However, molecular mechanisms behind this malfunction have remained elusive. None of the BIN1 disease mutants displayed a significantly compromised curvature sensing ability. However, two mutants showed impaired membrane tubulation both in vivo and in vitro, and displayed characteristically different behaviors. R154Q generated smaller membrane curvature compared to WT N-BAR. Quantification of protein density on membranes revealed a lower membrane-bound density for R154Q compared to WT and the other mutants, which appeared to be the primary reason for the observation of impaired deformation capacity. The D151N mutant was unable to tubulate liposomes under certain experimental conditions. At medium protein concentrations we found 'budding' structures on liposomes that we hypothesized to be intermediates during the tubulation process except for the D151N mutant. Chemical crosslinking assays suggested that the D151N mutation impaired protein oligomerization upon membrane binding. Although we found an insignificant difference between WT and K35N N-BAR in in vitro assays, depolymerizing actin in live cells allowed tubulation of plasma membranes through the K35N mutant. Our results provide insights into the membrane-involved pathophysiological mechanisms leading to human disease.

LSD1 Neurospecific Alternative Splicing Controls Neuronal Excitability in Mouse Models of Epilepsy

Alternative splicing in the brain is dynamic and instrumental to adaptive changes in response to stimuli. Lysine-specific demethylase 1 (LSD1/KDM1A) is a ubiquitously expressed histone H3Lys4 demethylase that acts as a transcriptional co-repressor in complex with its molecular partners CoREST and HDAC1/2. In mammalian brain, alternative splicing of LSD1 mini-exon E8a gives rise to neuroLSD1, a neurospecific isoform that, upon phosphorylation, acts as a dominant-negative causing disassembly of the co-repressor complex and de-repression of target genes. Here we show that the LSD1/neuroLSD1 ratio changes in response to neuronal activation and such effect is mediated by neurospecific splicing factors NOVA1 and nSR100/SRRM4 together with a novel cis-silencer. Indeed, we found that, in response to epileptogenic stimuli, downregulation of NOVA1 reduces exon E8a splicing and expression of neuroLSD1. Using behavioral and EEG analyses we observed that neuroLSD1-specific null mice are hypoexcitable and display decreased seizure susceptibility. Conversely, in a mouse model of Rett syndrome characterized by hyperexcitability, we measured higher levels of NOVA1 protein and upregulation of neuroLSD1. In conclusion, we propose that, in the brain, correct ratio between LSD1 and neuroLSD1 contributes to excitability and, when altered, could represent a pathogenic event associated with neurological disorders involving altered E/I.

Analysis of the robustness of network-based disease-gene prioritization methods reveals redundancy in the human interactome and functional diversity of disease-genes

Complex biological systems usually pose a trade-off between robustness and fragility where a small number of perturbations can substantially disrupt the system. Although biological systems are robust against changes in many external and internal conditions, even a single mutation can perturb the system substantially, giving rise to a pathophenotype. Recent advances in identifying and analyzing the sequential variations beneath human disorders help to comprehend a systemic view of the mechanisms underlying various disease phenotypes. Network-based disease-gene prioritization methods rank the relevance of genes in a disease under the hypothesis that genes whose proteins interact with each other tend to exhibit similar phenotypes. In this study, we have tested the robustness of several network-based disease-gene prioritization methods with respect to the perturbations of the system using various disease phenotypes from the Online Mendelian Inheritance in Man database. These perturbations have been introduced either in the protein-protein interaction network or in the set of known disease-gene associations. As the network-based disease-gene prioritization methods are based on the connectivity between known disease-gene associations, we have further used these methods to categorize the pathophenotypes with respect to the recoverability of hidden disease-genes. Our results have suggested that, in general, disease-genes are connected through multiple paths in the human interactome. Moreover, even when these paths are disturbed, network-based prioritization can reveal hidden disease-gene associations in some pathophenotypes such as breast cancer, cardiomyopathy, diabetes, leukemia, parkinson disease and obesity to a greater extend compared to the rest of the pathophenotypes tested in this study. Gene Ontology (GO) analysis highlighted the role of functional diversity for such diseases.

Phosphorylation of SRSF1 by SRPK1 regulates alternative splicing of tumor-related Rac1b in colorectal cells

The premessenger RNA of the majority of human genes can generate various transcripts through alternative splicing, and different tissues or disease states show specific patterns of splicing variants. These patterns depend on the relative concentrations of the splicing factors present in the cell nucleus, either as a consequence of their expression levels or of post-translational modifications, such as protein phosphorylation, which are determined by signal transduction pathways. Here, we analyzed the contribution of protein kinases to the regulation of alternative splicing variant Rac1b that is overexpressed in certain tumor types. In colorectal cells, we found that depletion of AKT2, AKT3, GSK3β, and SRPK1 significantly decreased endogenous Rac1b levels. Although knockdown of AKT2 and AKT3 affected only Rac1b protein levels suggesting a post-splicing effect, the depletion of GSK3β or SRPK1 decreased Rac1b alternative splicing, an effect mediated through changes in splicing factor SRSF1. In particular, the knockdown of SRPK1 or inhibition of its catalytic activity reduced phosphorylation and subsequent translocation of SRSF1 to the nucleus, limiting its availability to promote the inclusion of alternative exon 3b into the Rac1 pre-mRNA. Altogether, the data identify SRSF1 as a prime regulator of Rac1b expression in colorectal cells and provide further mechanistic insight into how the regulation of alternative splicing events by protein kinases can contribute to sustain tumor cell survival.

Illuminating the Transcriptome through the Genome

Sequencing the human genome was a huge milestone in genetic research that revealed almost the total DNA sequence required to create a human being. However, in order to function, the DNA genome needs to be expressed as an RNA transcriptome. This article reviews how knowledge of genome sequence information has led to fundamental discoveries in how the transcriptome is processed, with a focus on new system-wide insights into how pre-mRNAs that are encoded by split genes in the genome are rearranged by splicing into functional mRNAs. These advances have been made possible by the development of new post-genome technologies to probe splicing patterns. Transcriptome-wide approaches have characterised a "splicing code" that is embedded within and has a significant role in deciphering the genome, and is deciphered by RNA binding proteins. These analyses have also found that most human genes encode multiple mRNA isoforms, and in some cases proteins, leading in turn to a re-assessment of what exactly a gene is. Analysis of the transcriptome has given insights into how the genome is packaged and transcribed, and is helping to explain important aspects of genome evolution.

Intrinsically disordered proteins and intrinsically disordered protein regions

Intrinsically disordered proteins (IDPs) and IDP regions fail to form a stable structure, yet they exhibit biological activities. Their mobile flexibility and structural instability are encoded by their amino acid sequences. They recognize proteins, nucleic acids, and other types of partners; they accelerate interactions and chemical reactions between bound partners; and they help accommodate posttranslational modifications, alternative splicing, protein fusions, and insertions or deletions. Overall, IDP-associated biological activities complement those of structured proteins. Recently, there has been an explosion of studies on IDP regions and their functions, yet the discovery and investigation of these proteins have a long, mostly ignored history. Along with recent discoveries, we present several early examples and the mechanisms by which IDPs contribute to function, which we hope will encourage comprehensive discussion of IDPs and IDP regions in biochemistry textbooks. Finally, we propose future directions for IDP research.

CHIP-MYTH: a novel interactive proteomics method for the assessment of agonist-dependent interactions of the human β₂-adrenergic receptor

G-protein coupled receptors (GPCRs) are involved in a variety of disease processes and comprise major drug targets. However, the complexity of integral membrane proteins such as GPCRs makes the identification of their interacting partners and subsequent drug development challenging. A comprehensive understanding of GPCR protein interaction networks is needed to design effective therapeutic strategies to inhibit these drug targets. Here, we developed a novel split-ubiquitin membrane yeast two-hybrid (MYTH) technology called CHIP-MYTH, which allows the unbiased characterization of interaction partners of full-length GPCRs in a drug-dependent manner. This was achieved by coupling DNA microarray technology to the MYTH approach, which allows a quantitative evaluation of interacting partners of a given integral membrane protein in the presence or absence of drug. As a proof of principle, we applied the CHIP-MYTH approach to the human β2-adrenergic receptor (β2AR), a target of interest in the treatment of asthma, chronic obstructive pulmonary disease (COPD), neurological disease, cardiovascular disease, and obesity. A CHIP-MYTH screen was performed in the presence or absence of salmeterol, a long-acting β2AR-agonist. Our results suggest that β2AR activation with salmeterol can induce the dissociation of heterotrimeric G-proteins, Gαβγ, into Gα and Gβγ subunits, which in turn activates downstream signaling cascades. Using CHIP-MYTH, we confirmed previously known and identified novel β2AR interactors involved in GPCR-mediated signaling cascades. Several of these interactions were confirmed in mammalian cells using LUminescence-based Mammalian IntERactome (LUMIER) and co-immunoprecipitation assays. In summary, the CHIP-MYTH approach is ideal for conducting comprehensive protein-protein interactions (PPI) screenings of full-length GPCRs in the presence or absence of drugs, thus providing a valuable tool to further our understanding of GPCR-mediated signaling.

Prioritizing protein complexes implicated in human diseases by network optimization

Background

The detection of associations between protein complexes and human inherited diseases is of great importance in understanding mechanisms of diseases. Dysfunctions of a protein complex are usually defined by its member disturbance and consequently result in certain diseases. Although individual disease proteins have been widely predicted, computational methods are still absent for systematically investigating disease-related protein complexes.

Results

We propose a method, MAXCOM, for the prioritization of candidate protein complexes. MAXCOM performs a maximum information flow algorithm to optimize relationships between a query disease and candidate protein complexes through a heterogeneous network that is constructed by combining protein-protein interactions and disease phenotypic similarities. Cross-validation experiments on 539 protein complexes show that MAXCOM can rank 382 (70.87%) protein complexes at the top against protein complexes constructed at random. Permutation experiments further confirm that MAXCOM is robust to the network structure and parameters involved. We further analyze protein complexes ranked among top ten for breast cancer and demonstrate that the SWI/SNF complex is potentially associated with breast cancer.

Conclusions

MAXCOM is an effective method for the discovery of disease-related protein complexes based on network optimization. The high performance and robustness of this approach can facilitate not only pathologic studies of diseases, but also the design of drugs targeting on multiple proteins.

A survey of software for genome-wide discovery of differential splicing in RNA-Seq data

Alternative splicing is a major contributor to cellular diversity. Therefore the identification and quantification of differentially spliced transcripts in genome-wide transcript analysis is an important consideration. Here, I review the software available for analysis of RNA-Seq data for differential splicing and discuss intrinsic challenges for differential splicing analyses. Three approaches to differential splicing analysis are described, along with their associated software implementations, their strengths, limitations, and caveats. Suggestions for future work include more extensive experimental validation to assess accuracy of the software predictions and consensus formats for outputs that would facilitate visualizations, data exchange, and downstream analyses.

Alternative splicing in cancer: implications for biology and therapy

Alternative splicing has critical roles in normal development and can promote growth and survival in cancer. Aberrant splicing, the production of noncanonical and cancer-specific mRNA transcripts, can lead to loss-of-function in tumor suppressors or activation of oncogenes and cancer pathways. Emerging data suggest that aberrant splicing products and loss of canonically spliced variants correlate with stage and progression in malignancy. Here, we review the splicing landscape of TP53, BARD1 and AR to illuminate roles for alternative splicing in cancer. We also examine the intersection between alternative splicing pathways and novel therapeutic approaches.

Splicing in the human brain

It has become increasingly clear over the past decade that RNA has important functions in human cells beyond its role as an intermediate translator of DNA to protein. It is now known that RNA plays highly specific roles in pathways involved in regulatory, structural, and catalytic functions. The complexity of RNA production and regulation has become evident with the advent of high-throughput methods to study the transcriptome. Deep sequencing has revealed an enormous diversity of RNA types and transcript isoforms in human cells. The transcriptome of the human brain is particularly interesting as it contains more expressed genes than other tissues and also displays an extreme diversity of transcript isoforms, indicating that highly complex regulatory pathways are present in the brain. Several of these regulatory proteins are now identified, including RNA-binding proteins that are neuron specific. RNA-binding proteins also play important roles in regulating the splicing process and the temporal and spatial isoform production. While significant progress has been made in understanding the human transcriptome, many questions still remain regarding the basic mechanisms of splicing and subcellular localization of RNA. A long-standing question is to what extent the splicing of pre-mRNA is cotranscriptional and posttranscriptional, respectively. Recent data, including studies of the human brain, indicate that splicing is primarily cotranscriptional in human cells. This chapter describes the current understanding of splicing and splicing regulation in the human brain and discusses the recent global sequence-based analyses of transcription and splicing.

High-resolution functional annotation of human transcriptome: predicting isoform functions by a novel multiple instance-based label propagation method

Alternative transcript processing is an important mechanism for generating functional diversity in genes. However, little is known about the precise functions of individual isoforms. In fact, proteins (translated from transcript isoforms), not genes, are the function carriers. By integrating multiple human RNA-seq data sets, we carried out the first systematic prediction of isoform functions, enabling high-resolution functional annotation of human transcriptome. Unlike gene function prediction, isoform function prediction faces a unique challenge: the lack of the training data--all known functional annotations are at the gene level. To address this challenge, we modelled the gene-isoform relationships as multiple instance data and developed a novel label propagation method to predict functions. Our method achieved an average area under the receiver operating characteristic curve of 0.67 and assigned functions to 15 572 isoforms. Interestingly, we observed that different functions have different sensitivities to alternative isoform processing, and that the function diversity of isoforms from the same gene is positively correlated with their tissue expression diversity. Finally, we surveyed the literature to validate our predictions for a number of apoptotic genes. Strikingly, for the famous 'TP53' gene, we not only accurately identified the apoptosis regulation function of its five isoforms, but also correctly predicted the precise direction of the regulation.

Regulated aggregative multicellularity in a close unicellular relative of metazoa

The evolution of metazoans from their unicellular ancestors was one of the most important events in the history of life. However, the cellular and genetic changes that ultimately led to the evolution of multicellularity are not known. In this study, we describe an aggregative multicellular stage in the protist Capsaspora owczarzaki, a close unicellular relative of metazoans. Remarkably, transition to the aggregative stage is associated with significant upregulation of orthologs of genes known to establish multicellularity and tissue architecture in metazoans. We further observe transitions in regulated alternative splicing during the C. owczarzaki life cycle, including the deployment of an exon network associated with signaling, a feature of splicing regulation so far only observed in metazoans. Our results reveal the existence of a highly regulated aggregative stage in C. owczarzaki and further suggest that features of aggregative behavior in an ancestral protist may had been co-opted to develop some multicellular properties currently seen in metazoans.

DOI:http://dx.doi.org/10.7554/eLife.01287.001

When living things made from many cells evolved from single-celled ancestors, it was a breakthrough in the history of life—and one that has occurred more than once. In fact, multicellular life has evolved independently at least 25 times, in groups as diverse as animals, fungi, plants, slime molds and seaweeds. There are broadly two ways to become multicellular. The most complex multicellular species, such as animals, will replicate a single cell, over and over, without separating the resultant cells. However, in species that are only occasionally multicellular, free-living cells tend instead to join together in one mass of many cells.

Evolution is constrained by its raw materials; so looking at the living relatives of a given species, or group, can lead to a better understanding of its evolution because its relatives contain clues about its ancestors. To gain insights into how animal multicellular life might have began; Sebé-Pedrós et al. studied the life cycle of the amoeboid organism Capsaspora owczarzaki. Found within the bodies of freshwater snails, this single-celled amoeba is a close relative of multicellular animals and could resemble one of their earliest ancestors.

At certain stages of the life cycle Sebé-Pedrós et al. noticed that individual amoebae gathered together to form a multicellular mass—something that had not been seen before in such a close relative of the animals. Moreover, the genes that ‘switched on’ when the amoebae began to aggregate are also found in animals; where, together with other genes, they control development and the formation of tissues. Sebé-Pedrós et al. suggest that the first multicellular animals could have recycled the genes that control the aggregation of single-celled species: in other words, genes that once controlled the changes that happen at different times in a life cycle, now control the changes that develop between different tissues at the same time.

Sebé-Pedrós et al. also observed that alternative splicing—a process that allows different proteins to be made from a single gene—occurs via two different mechanisms during the life cycle of Capsaspora. Most of the time Capsaspora employs a form of alternative splicing that is often seen in plants and fungi, and only rarely in animals; for the rest of the time it uses a form of alternative splicing similar to that used by animal cells.

The evolution of complex alternative splicing mechanisms is a hallmark feature of multicellular animals. The exploitation of two major forms of alternative splicing in Capsaspora could thus reflect an important transition during evolution that resulted in an increased diversity of proteins and in more complex gene regulation. Such a transition may ultimately have paved the way for the increased specialization of cell types seen in animals.

This glimpse into the possible transitions in gene regulation that contributed to the birth of multicellular animals indicates that the single-celled ancestor of the animals was likely more complex than previously thought. Future analyses of the animals’ close relatives may further improve our understanding of how single-celled organisms became multicellular animals.

DOI:http://dx.doi.org/10.7554/eLife.01287.002

Integrative approaches for finding modular structure in biological networks

A central goal of systems biology is to elucidate the structural and functional architecture of the cell. To this end, large and complex networks of molecular interactions are being rapidly generated for humans and model organisms. A recent focus of bioinformatics research has been to integrate these networks with each other and with diverse molecular profiles to identify sets of molecules and interactions that participate in a common biological function - that is, 'modules'. Here, we classify such integrative approaches into four broad categories, describe their bioinformatic principles and review their applications.

Identification and characterization of a neuron-specific isoform of protrudin

Protrudin is a membrane protein that regulates polarized vesicular transport. Now, we have identified a novel isoform of protrudin (protrudin-L) that contains an additional seven amino acids between the FFAT motif and the coiled-coil domain compared with the conventional isoform (protrudin-S) as a result of alternative splicing of a microexon (exon L). Protrudin-L mRNA was found to be mostly restricted to the central nervous system in mice, whereas protrudin-S mRNA was detected in all tissues examined. With the use of a splicing reporter minigene that produces two distinct fluorescent proteins in a manner dependent on the splicing pattern of protrudin transcripts, we found that most neurons express protrudin-L, whereas astrocytes express both protrudin isoforms and oligodendrocytes express only protrudin-S. Protrudin-L associated to a greater extent with vesicle-associated membrane protein-associated protein (VAP) than protrudin-S. Expression of protrudin-L in hippocampal neurons of protrudin-deficient mice also promoted neurite outgrowth more efficiently than protrudin-S. Our results suggest that protrudin-L is a neuron-specific protrudin isoform that promotes axonal elongation and contributes to the establishment of neuronal polarity.

Promiscuity as a functional trait: intrinsically disordered regions as central players of interactomes

Because of their pervasiveness in eukaryotic genomes and their unique properties, understanding the role that ID (intrinsically disordered) regions in proteins play in the interactome is essential for gaining a better understanding of the network. Especially critical in determining this role is their ability to bind more than one partner using the same region. Studies have revealed that proteins containing ID regions tend to take a central role in protein interaction networks; specifically, they act as hubs, interacting with multiple different partners across time and space, allowing for the co-ordination of many cellular activities. There appear to be three different modules within ID regions responsible for their functionally promiscuous behaviour: MoRFs (molecular recognition features), SLiMs (small linear motifs) and LCRs (low complexity regions). These regions allow for functionality such as engaging in the formation of dynamic heteromeric structures which can serve to increase local activity of an enzyme or store a collection of functionally related molecules for later use. However, the use of promiscuity does not come without a cost: a number of diseases that have been associated with ID-containing proteins seem to be caused by undesirable interactions occurring upon altered expression of the ID-containing protein.

Postnatal isoform switch and protein localization of LEF1 and TCF7L2 transcription factors in cortical, thalamic, and mesencephalic regions of the adult mouse brain

β-Catenin signaling, leading to the activation of lymphoid enhancer-binding factor 1/T cell factor (LEF1/TCF) transcription factors, plays a well-established role in transcription regulation during development and tissue homeostasis. In the adult organism, the activity of this pathway has been found in stem cell niches and postmitotic thalamic neurons. Recently, studies show that mutations in components of β-catenin signaling networks have been associated with several psychiatric disorders, indicating the involvement of β-catenin and LEF1/TCF proteins in the proper functioning of the brain. Here, we report a comprehensive analysis of LEF1/TCF protein localization and the expression profile of their isoforms in cortical, thalamic, and midbrain regions in mice. We detected LEF1 and TCF7L2 proteins in neurons of the thalamus and dorsal midbrain, i.e., subcortical regions specialized in the integration of diverse sources of sensory information. These neurons also exhibited nuclear localization of β-catenin, suggesting the involvement of β-catenin/LEF1/TCF7L2 in the regulation of gene expression in these regions. Analysis of alternative splicing and promoter usage identified brain-specific TCF7L2 isoforms and revealed a developmentally coordinated transition in the composition of LEF1 and TCF7L2 isoforms. In the case of TCF7L2, the typical brain isoforms lack the so-called C clamp; in addition, the dominant-negative isoforms are predominant in the embryonic thalamus but disappear postnatally. The present study provides a necessary framework to understand the role of LEF1/TCF factors in thalamic and midbrain development until adulthood and predicts that the regulatory role of these proteins in the adult brain is significantly different from their role in the embryonic brain or other non-neural tissues.

An integrative framework identifies alternative splicing events in colorectal cancer development

Alternative splicing (AS) is a common mechanism which creates diverse RNA isoforms from a single gene, potentially increasing protein variety. Growing evidence suggests that this mechanism is closely related to cancer progression. In this study, whole transcriptome analysis was performed with GeneChip Human exon 1.0 ST Array from 80 samples comprising 23 normal colon mucosa, 30 primary colorectal cancer and 27 liver metastatic specimens from 46 patients, to identify AS events in colorectal cancer progression. Differentially expressed genes and exons were estimated and AS events were reconstructed by combining exon-level analyses with AltAnalyze algorithms and transcript-level estimations (MMBGX probabilistic method). The number of AS genes in the transition from normal colon mucosa to primary tumor was the most abundant, but fell considerably in the next transition to liver metastasis. 206 genes with probable AS events in colon cancer development and progression were identified, that are involved in processes and pathways relevant to tumor biology, as cell-cell and cell-matrix interactions. Several AS events in VCL, CALD1, B3GNT6 and CTHRC1 genes, differentially expressed during tumor development were validated, at RNA and at protein level. Taken together, these results demonstrate that cancer-specific AS is common in early phases of colorectal cancer natural history.

Complexity of the alternative splicing landscape in plants

Alternative splicing (AS) of precursor mRNAs (pre-mRNAs) from multiexon genes allows organisms to increase their coding potential and regulate gene expression through multiple mechanisms. Recent transcriptome-wide analysis of AS using RNA sequencing has revealed that AS is highly pervasive in plants. Pre-mRNAs from over 60% of intron-containing genes undergo AS to produce a vast repertoire of mRNA isoforms. The functions of most splice variants are unknown. However, emerging evidence indicates that splice variants increase the functional diversity of proteins. Furthermore, AS is coupled to transcript stability and translation through nonsense-mediated decay and microRNA-mediated gene regulation. Widespread changes in AS in response to developmental cues and stresses suggest a role for regulated splicing in plant development and stress responses. Here, we review recent progress in uncovering the extent and complexity of the AS landscape in plants, its regulation, and the roles of AS in gene regulation. The prevalence of AS in plants has raised many new questions that require additional studies. New tools based on recent technological advances are allowing genome-wide analysis of RNA elements in transcripts and of chromatin modifications that regulate AS. Application of these tools in plants will provide significant new insights into AS regulation and crosstalk between AS and other layers of gene regulation.

The relationship between gene isoform multiplicity, number of exons and protein divergence

At present we know that phenotypic differences between organisms arise from a variety of sources, like protein sequence divergence, regulatory sequence divergence, alternative splicing, etc. However, we do not have yet a complete view of how these sources are related. Here we address this problem, studying the relationship between protein divergence and the ability of genes to express multiple isoforms. We used three genome-wide datasets of human-mouse orthologs to study the relationship between isoform multiplicity co-occurrence between orthologs (the fact that two orthologs have more than one isoform) and protein divergence. In all cases our results showed that there was a monotonic dependence between these two properties. We could explain this relationship in terms of a more fundamental one, between exon number of the largest isoform and protein divergence. We found that this last relationship was present, although with variations, in other species (chimpanzee, cow, rat, chicken, zebrafish and fruit fly). In summary, we have identified a relationship between protein divergence and isoform multiplicity co-occurrence and explained its origin in terms of a simple gene-level property. Finally, we discuss the biological implications of these findings for our understanding of inter-species phenotypic differences.

Drift and conservation of differential exon usage across tissues in primate species

Alternative usage of exons provides genomes with plasticity to produce different transcripts from the same gene, modulating the function, localization, and life cycle of gene products. It affects most human genes. For a limited number of cases, alternative functions and tissue-specific roles are known. However, recent high-throughput sequencing studies have suggested that much alternative isoform usage across tissues is nonconserved, raising the question of the extent of its functional importance. We address this question in a genome-wide manner by analyzing the transcriptomes of five tissues for six primate species, focusing on exons that are 1:1 orthologous in all six species. Our results support a model in which differential usage of exons has two major modes: First, most of the exons show only weak differences, which are dominated by interspecies variability and may reflect neutral drift and noisy splicing. These cases dominate the genome-wide view and explain why conservation appears to be so limited. Second, however, a sizeable minority of exons show strong differences between tissues, which are mostly conserved. We identified a core set of 3,800 exons from 1,643 genes that show conservation of strongly tissue-dependent usage patterns from human at least to macaque. This set is enriched for exons encoding protein-disordered regions and untranslated regions. Our findings support the theory that isoform regulation is an important target of evolution in primates, and our method provides a powerful tool for discovering potentially functional tissue-dependent isoforms.

Extracting reaction networks from databases-opening Pandora's box

Large quantities of information describing the mechanisms of biological pathways continue to be collected in publicly available databases. At the same time, experiments have increased in scale, and biologists increasingly use pathways defined in online databases to interpret the results of experiments and generate hypotheses. Emerging computational techniques that exploit the rich biological information captured in reaction systems require formal standardized descriptions of pathways to extract these reaction networks and avoid the alternative: time-consuming and largely manual literature-based network reconstruction. Here, we systematically evaluate the effects of commonly used knowledge representations on the seemingly simple task of extracting a reaction network describing signal transduction from a pathway database. We show that this process is in fact surprisingly difficult, and the pathway representations adopted by various knowledge bases have dramatic consequences for reaction network extraction, connectivity, capture of pathway crosstalk and in the modelling of cell-cell interactions. Researchers constructing computational models built from automatically extracted reaction networks must therefore consider the issues we outline in this review to maximize the value of existing pathway knowledge.

Alternative pre-mRNA splicing in neurons: growing up and extending its reach

Alternative pre-mRNA splicing determines the protein output of most neuronally expressed genes. Many examples have been described of protein function being modulated by coding changes in different mRNA isoforms. Several recent studies demonstrate that, through the coupling of splicing to other processes of mRNA metabolism, alternative splicing can also act as an on/off switch for gene expression. Other regulated splicing events may determine how an mRNA is utilized in its later cytoplasmic life by changing its localization or translation. These studies make clear that the multiple steps of post-transcriptional gene regulation are strongly linked. Together, these regulatory process play key roles in all aspects of the cell biology of neurons, from their initial differentiation, to their choice of connections, and finally to their function with mature circuits.

Intragenic epigenetic changes modulate NCAM alternative splicing in neuronal differentiation

Alternative splicing contributes to cell type-specific transcriptomes. Here, we show that changes in intragenic chromatin marks affect NCAM (neural cell adhesion molecule) exon 18 (E18) alternative splicing during neuronal differentiation. An increase in the repressive marks H3K9me2 and H3K27me3 along the gene body correlated with inhibition of polymerase II elongation in the E18 region, but without significantly affecting total mRNA levels. Treatment with the general DNA methylation inhibitor 5-azacytidine and BIX 01294, a specific inhibitor of H3K9 dimethylation, inhibited the differentiation-induced E18 inclusion, pointing to a role for repressive marks in sustaining NCAM splicing patterns typical of mature neurons. We demonstrate that intragenic deployment of repressive chromatin marks, induced by intronic small interfering RNAs targeting NCAM intron 18, promotes E18 inclusion in undifferentiated N2a cells, confirming the chromatin changes observed upon differentiation to be sufficient to induce alternative splicing. Combined with previous evidence that neuronal depolarization causes H3K9 acetylation and subsequent E18 skipping, our results show how two alternative epigenetic marks regulate NCAM alternative splicing and E18 levels in different cellular contexts.

GLiMMPS: robust statistical model for regulatory variation of alternative splicing using RNA-seq data

To characterize the genetic variation of alternative splicing, we develop GLiMMPS, a robust statistical method for detecting splicing quantitative trait loci (sQTLs) from RNA-seq data. GLiMMPS takes into account the individual variation in sequencing coverage and the noise prevalent in RNA-seq data. Analyses of simulated and real RNA-seq datasets demonstrate that GLiMMPS outperforms competing statistical models. Quantitative RT-PCR tests of 26 randomly selected GLiMMPS sQTLs yielded a validation rate of 100%. As population-scale RNA-seq studies become increasingly affordable and popular, GLiMMPS provides a useful tool for elucidating the genetic variation of alternative splicing in humans and model organisms.

Transcriptome analysis of human tissues and cell lines reveals one dominant transcript per gene

BACKGROUND

RNA sequencing has opened new avenues for the study of transcriptome composition. Significant evidence has accumulated showing that the human transcriptome contains in excess of a hundred thousand different transcripts. However, it is still not clear to what extent this diversity prevails when considering the relative abundances of different transcripts from the same gene.

RESULTS

Here we show that, in a given condition, most protein coding genes have one major transcript expressed at significantly higher level than others, that in human tissues the major transcripts contribute almost 85 percent to the total mRNA from protein coding loci, and that often the same major transcript is expressed in many tissues. We detect a high degree of overlap between the set of major transcripts and a recently published set of alternatively spliced transcripts that are predicted to be translated utilizing proteomic data. Thus, we hypothesize that although some minor transcripts may play a functional role, the major ones are likely to be the main contributors to the proteome. However, we still detect a non-negligible fraction of protein coding genes for which the major transcript does not code a protein.

CONCLUSIONS

Overall, our findings suggest that the transcriptome from protein coding loci is dominated by one transcript per gene and that not all the transcripts that contribute to transcriptome diversity are equally likely to contribute to protein diversity. This observation can help to prioritize candidate targets in proteomics research and to predict the functional impact of the detected changes in variation studies.

Alternative splicing tends to avoid partial removals of protein-protein interaction sites

BACKGROUND

Anecdotal evidence of the involvement of alternative splicing (AS) in the regulation of protein-protein interactions has been reported by several studies. AS events have been shown to significantly occur in regions where a protein interaction domain or a short linear motif is present. Several AS variants show partial or complete loss of interface residues, suggesting that AS can play a major role in the interaction regulation by selectively targeting the protein binding sites. In the present study we performed a statistical analysis of the alternative splicing of a non-redundant dataset of human protein-protein interfaces known at molecular level to determine the importance of this way of modulation of protein-protein interactions through AS.

RESULTS

Using a Cochran-Mantel-Haenszel chi-square test we demonstrated that the alternative splicing-mediated partial removal of both heterodimeric and homodimeric binding sites occurs at lower frequencies than expected, and this holds true even if we consider only those isoforms whose sequence is less different from that of the canonical protein and which therefore allow to selectively regulate functional regions of the protein. On the other hand, large removals of the binding site are not significantly prevented, possibly because they are associated to drastic structural changes of the protein. The observed protection of the binding sites from AS is not preferentially directed towards putative hot spot interface residues, and is widespread to all protein functional classes.

CONCLUSIONS

Our findings indicate that protein-protein binding sites are generally protected from alternative splicing-mediated partial removals. However, some cases in which the binding site is selectively removed exist, and here we discuss one of them.

CERKL knockdown causes retinal degeneration in zebrafish

The human CERKL gene is responsible for common and severe forms of retinal dystrophies. Despite intense in vitro studies at the molecular and cellular level and in vivo analyses of the retina of murine knockout models, CERKL function remains unknown. In this study, we aimed to approach the developmental and functional features of cerkl in Danio rerio within an Evo-Devo framework. We show that gene expression increases from early developmental stages until the formation of the retina in the optic cup. Unlike the high mRNA-CERKL isoform multiplicity shown in mammals, the moderate transcriptional complexity in fish facilitates phenotypic studies derived from gene silencing. Moreover, of relevance to pathogenicity, teleost CERKL shares the two main human protein isoforms. Morpholino injection has been used to generate a cerkl knockdown zebrafish model. The morphant phenotype results in abnormal eye development with lamination defects, failure to develop photoreceptor outer segments, increased apoptosis of retinal cells and small eyes. Our data support that zebrafish Cerkl does not interfere with proliferation and neural differentiation during early developmental stages but is relevant for survival and protection of the retinal tissue. Overall, we propose that this zebrafish model is a powerful tool to unveil CERKL contribution to human retinal degeneration.

Dynamic integration of splicing within gene regulatory pathways

Precursor mRNA splicing is one of the most highly regulated processes in metazoan species. In addition to generating vast repertoires of RNAs and proteins, splicing has a profound impact on other gene regulatory layers, including mRNA transcription, turnover, transport, and translation. Conversely, factors regulating chromatin and transcription complexes impact the splicing process. This extensive crosstalk between gene regulatory layers takes advantage of dynamic spatial, physical, and temporal organizational properties of the cell nucleus, and further emphasizes the importance of developing a multidimensional understanding of splicing control.

Distinct types of disorder in the human proteome: functional implications for alternative splicing

Intrinsically disordered regions have been associated with various cellular processes and are implicated in several human diseases, but their exact roles remain unclear. We previously defined two classes of conserved disordered regions in budding yeast, referred to as "flexible" and "constrained" conserved disorder. In flexible disorder, the property of disorder has been positionally conserved during evolution, whereas in constrained disorder, both the amino acid sequence and the property of disorder have been conserved. Here, we show that flexible and constrained disorder are widespread in the human proteome, and are particularly common in proteins with regulatory functions. Both classes of disordered sequences are highly enriched in regions of proteins that undergo tissue-specific (TS) alternative splicing (AS), but not in regions of proteins that undergo general (i.e., not tissue-regulated) AS. Flexible disorder is more highly enriched in TS alternative exons, whereas constrained disorder is more highly enriched in exons that flank TS alternative exons. These latter regions are also significantly more enriched in potential phosphosites and other short linear motifs associated with cell signaling. We further show that cancer driver mutations are significantly enriched in regions of proteins associated with TS and general AS. Collectively, our results point to distinct roles for TS alternative exons and flanking exons in the dynamic regulation of protein interaction networks in response to signaling activity, and they further suggest that alternatively spliced regions of proteins are often functionally altered by mutations responsible for cancer.

Integrative avenues for exploring the dynamics and evolution of protein interaction networks

Over the past decade, the study of protein interaction networks (PINs) has shed light on the organizing principles of living cells. However, PINs have been mostly mapped in one single condition. We outline three of the most promising avenues of investigation in this field, namely the study of first, how PINs are rewired by mutations and environmental perturbations; secondly, how inter-species interactions affect PIN achitectures; thirdly, what mechanisms and forces drive PIN evolution. These investigations will unravel the dynamics and condition dependence of PINs and will thus lead to a better functional annotation of network architecture. One major challenge to reach these goals is the integration of PINs with other cellular regulatory networks in the context of complex cellular phenotypes.

The impact of splicing on protein domain architecture

Many proteins are composed of protein domains, functional units of common descent. Multidomain forms are common in all eukaryotes making up more than half of the proteome and the evolution of novel domain architecture has been accelerated in metazoans. It is also becoming increasingly clear that alternative splicing is prevalent among vertebrates. Given that protein domains are defined as structurally, functionally and evolutionarily distinct units, one may speculate that some alternative splicing events may lead to clean excisions of protein domains, thus generating a number of different domain architectures from one gene template. However, recent findings indicate that smaller alternative splicing events, in particular in disordered regions, might be more prominent than domain architectural changes. The problem of identifying protein isoforms is, however, still not resolved. Clearly, many splice forms identified through detection of mRNA sequences appear to produce 'nonfunctional' proteins, such as proteins with missing internal secondary structure elements. Here, we review the state of the art methods for identification of functional isoforms and present a summary of what is known, thus far, about alternative splicing with regard to protein domain architectures.

The switches.ELM resource: a compendium of conditional regulatory interaction interfaces

Short linear motifs (SLiMs) are protein interaction sites that play an important role in cell regulation by controlling protein activity, localization, and local abundance. The functionality of a SLiM can be modulated in a context-dependent manner to induce a gain, loss, or exchange of binding partners, which will affect the function of the SLiM-containing protein. As such, these conditional interactions underlie molecular decision-making in cell signaling. We identified multiple types of pre- and posttranslational switch mechanisms that can regulate the function of a SLiM and thereby control its interactions. The collected examples of experimentally characterized SLiM-based switch mechanisms were curated in the freely accessible switches.ELM resource (http://switches.elm.eu.org). On the basis of these examples, we defined and integrated rules to analyze SLiMs for putative regulatory switch mechanisms. We applied these rules to known validated SLiMs, providing evidence that more than half of these are likely to be pre- or posttranslationally regulated. In addition, we showed that posttranslationally modified sites are enriched around SLiMs, which enables cooperative and integrative regulation of protein interaction interfaces. We foresee switches.ELM complementing available resources to extend our knowledge of the molecular mechanisms underlying cell signaling.

Alternative splicing: a pivotal step between eukaryotic transcription and translation

Alternative splicing was discovered simultaneously with splicing over three decades ago. Since then, an enormous body of evidence has demonstrated the prevalence of alternative splicing in multicellular eukaryotes, its key roles in determining tissue- and species-specific differentiation patterns, the multiple post- and co-transcriptional regulatory mechanisms that control it, and its causal role in hereditary disease and cancer. The emerging evidence places alternative splicing in a central position in the flow of eukaryotic genetic information, between transcription and translation, in that it can respond not only to various signalling pathways that target the splicing machinery but also to transcription factors and chromatin structure.

Exploring the binding diversity of intrinsically disordered proteins involved in one-to-many binding

Molecular recognition features (MoRFs) are intrinsically disordered protein regions that bind to partners via disorder-to-order transitions. In one-to-many binding, a single MoRF binds to two or more different partners individually. MoRF-based one-to-many protein-protein interaction (PPI) examples were collected from the Protein Data Bank, yielding 23 MoRFs bound to 2-9 partners, with all pairs of same-MoRF partners having less than 25% sequence identity. Of these, 8 MoRFs were bound to 2-9 partners having completely different folds, whereas 15 MoRFs were bound to 2-5 partners having the same folds but with low sequence identities. For both types of partner variation, backbone and side chain torsion angle rotations were used to bring about the conformational changes needed to enable close fits between a single MoRF and distinct partners. Alternative splicing events (ASEs) and posttranslational modifications (PTMs) were also found to contribute to distinct partner binding. Because ASEs and PTMs both commonly occur in disordered regions, and because both ASEs and PTMs are often tissue-specific, these data suggest that MoRFs, ASEs, and PTMs may collaborate to alter PPI networks in different cell types. These data enlarge the set of carefully studied MoRFs that use inherent flexibility and that also use ASE-based and/or PTM-based surface modifications to enable the same disordered segment to selectively associate with two or more partners. The small number of residues involved in MoRFs and in their modifications by ASEs or PTMs may simplify the evolvability of signaling network diversity.

The evolutionary landscape of alternative splicing in vertebrate species

How species with similar repertoires of protein-coding genes differ so markedly at the phenotypic level is poorly understood. By comparing organ transcriptomes from vertebrate species spanning ~350 million years of evolution, we observed significant differences in alternative splicing complexity between vertebrate lineages, with the highest complexity in primates. Within 6 million years, the splicing profiles of physiologically equivalent organs diverged such that they are more strongly related to the identity of a species than they are to organ type. Most vertebrate species-specific splicing patterns are cis-directed. However, a subset of pronounced splicing changes are predicted to remodel protein interactions involving trans-acting regulators. These events likely further contributed to the diversification of splicing and other transcriptomic changes that underlie phenotypic differences among vertebrate species.

Plant protein interactomes

Protein-protein interactions are a critical element of biological systems, and the analysis of interaction partners can provide valuable hints about unknown functions of a protein. In recent years, several large-scale protein interaction studies have begun to unravel the complex networks through which plant proteins exert their functions. Two major classes of experimental approaches are used for protein interaction mapping: analysis of direct interactions using binary methods such as yeast two-hybrid or split ubiquitin, and analysis of protein complexes through affinity purification followed by mass spectrometry. In addition, bioinformatics predictions can suggest interactions that have evaded detection by other methods or those of proteins that have not been investigated. Here we review the major approaches to construct, analyze, use, and carry out quality control on plant protein interactome networks. We present experimental and computational approaches for large-scale mapping, methods for validation or smaller-scale functional studies, important bioinformatics resources, and findings from recently published large-scale plant interactome network maps.

Intrinsically disordered proteins and conformational noise: implications in cancer

Intrinsically disordered proteins, IDPs, are proteins that lack a rigid 3D structure under physiological conditions, at least in vitro. Despite the lack of structure, IDPs play important roles in biological processes and transition from disorder to order upon binding to their targets. With multiple conformational states and rapid conformational dynamics, they engage in myriad and often "promiscuous" interactions. These stochastic interactions between IDPs and their partners, defined here as conformational noise, is an inherent characteristic of IDP interactions. The collective effect of conformational noise is an ensemble of protein network configurations, from which the most suitable can be explored in response to perturbations, conferring protein networks with remarkable flexibility and resilience. Moreover, the ubiquitous presence of IDPs as transcriptional factors and, more generally, as hubs in protein networks, is indicative of their role in propagation of transcriptional (genetic) noise. As effectors of transcriptional and conformational noise, IDPs rewire protein networks and unmask latent interactions in response to perturbations. Thus, noise-driven activation of latent pathways could underlie state-switching events such as cellular transformation in cancer. To test this hypothesis, we created a model of a protein network with the topological characteristics of a cancer protein network and tested its response to a perturbation in presence of IDP hubs and conformational noise. Because numerous IDPs are found to be epigenetic modifiers and chromatin remodelers, we hypothesize that they could further channel noise into stable, heritable genotypic changes.