Resurrection of an Urbilaterian U1A/U2B″/SNF protein

Integrating Phylogenetics With Intron Positions Illuminates the Origin of the Complex Spliceosome

Eukaryotic genes are characterized by the presence of introns that are removed from pre-mRNA by a spliceosome. This ribonucleoprotein complex is comprised of multiple RNA molecules and over a hundred proteins, which makes it one of the most complex molecular machines that originated during the prokaryote-to-eukaryote transition. Previous works have established that these introns and the spliceosomal core originated from self-splicing introns in prokaryotes. Yet, how the spliceosomal core expanded by recruiting many additional proteins remains largely elusive. In this study, we use phylogenetic analyses to infer the evolutionary history of 145 proteins that we could trace back to the spliceosome in the last eukaryotic common ancestor. We found that an overabundance of proteins derived from ribosome-related processes was added to the prokaryote-derived core. Extensive duplications of these proteins substantially increased the complexity of the emerging spliceosome. By comparing the intron positions between spliceosomal paralogs, we infer that most spliceosomal complexity postdates the spread of introns through the proto-eukaryotic genome. The reconstruction of early spliceosomal evolution provides insight into the driving forces behind the emergence of complexes with many proteins during eukaryogenesis.

Phylogenetic analysis and stress response of the plant U2 small nuclear ribonucleoprotein B″ gene family

Background

Alternative splicing (AS) is an important channel for gene expression regulation and protein diversification, in addition to a major reason for the considerable differences in the number of genes and proteins in eukaryotes. In plants, U2 small nuclear ribonucleoprotein B″ (U2B″), a component of splicing complex U2 snRNP, plays an important role in AS. Currently, few studies have investigated plant U2B″, and its mechanism remains unclear.

Result

Phylogenetic analysis, including gene and protein structures, revealed that U2B″ is highly conserved in plants and typically contains two RNA recognition motifs. Subcellular localisation showed that OsU2B″ is located in the nucleus and cytoplasm, indicating that it has broad functions throughout the cell. Elemental analysis of the promoter region showed that it responded to numerous external stimuli, including hormones, stress, and light. Subsequent qPCR experiments examining response to stress (cold, salt, drought, and heavy metal cadmium) corroborated the findings. The prediction results of protein–protein interactions showed that its function is largely through a single pathway, mainly through interaction with snRNP proteins.

Conclusion

U2B″ is highly conserved in the plant kingdom, functions in the nucleus and cytoplasm, and participates in a wide range of processes in plant growth and development.

Phylogeny and conservation of plant U2A/U2A', a core splicing component in U2 spliceosomal complex

This study systematically identifies 112 U2A genes from 80 plant species by combinatory bioinformatics analysis, which is important for understanding their phylogenetic history, expression profiles and for predicting specific functions. In eukaryotes, a pre-mRNA can generate multiple transcripts by removing certain introns and joining corresponding exons, thus greatly expanding the transcriptome and proteome diversity. The spliceosome is a mega-Dalton ribonucleoprotein (RNP) complex that is essential for the process of splicing. In spliceosome components, the U2 small nuclear ribonucleoprotein (U2 snRNP) forms the pre-spliceosome by association with the branch site. An essential component that promotes U2 snRNP assembly, named U2A, has been extensively identified in humans, yeast and nematodes. However, studies examining U2A genes in plants are scarce. In this study, we performed a comprehensive analysis and identified a total of 112 U2A genes from 80 plant species representing dicots, monocots, mosses and algae. Comparisons of the gene structures, protein domains, and expression patterns of 112 U2A genes indicated that the conserved functions were likely retained by plant U2A genes and important for responses to internal and external stimuli. In addition, analysis of alternative transcripts and splice sites of U2A genes indicated that the fifth intron contained a conserved alternative splicing event that might be important for its molecular function. Our work provides a general understanding of this splicing factor family in terms of genes and proteins, and it will serve as a fundamental resource that will contribute to further mechanistic characterization in plants.

Characterization of Pre-mRNA Splicing and Spliceosomal Machinery in Porphyridium purpureum and Evolutionary Implications for Red Algae

Pre-mRNA splicing is a highly conserved eukaryotic process, but our understanding of it is limited by a historical focus on well-studied organisms such as humans and yeast. There is considerable diversity in mechanisms and components of pre-mRNA splicing, especially in lineages that have evolved under the pressures of genome reduction. The ancestor of red algae is thought to have undergone genome reduction prior to the lineage's radiation, resulting in overall gene and intron loss in extant groups. Previous studies on the extremophilic red alga Cyanidioschyzon merolae revealed an intron-sparse genome with a highly reduced spliceosome. To determine whether these features applied to other red algae, we investigated multiple aspects of pre-mRNA splicing in the mesophilic red alga Porphyridium purpureum. Through strand-specific RNA-Seq, we observed high levels of intron retention across a large number of its introns, and nearly half of the transcripts for these genes are not spliced at all. We also discovered a relationship between variability of 5' splice site sequences and levels of splicing. To further investigate the connections between intron retention and splicing machinery, we bioinformatically assembled the P. purpureum spliceosome, and biochemically verified the presence of snRNAs. While most other core spliceosomal components are present, our results suggest highly divergent or missing U1 snRNP proteins, despite the presence of an uncharacteristically long U1 snRNA. These unusual aspects highlight the diverse nature of pre-mRNA splicing that can be seen in lesser-studied eukaryotes, raising the importance of investigating fundamental eukaryotic processes outside of model organisms.

Molecular principles underlying dual RNA specificity in the Drosophila SNF protein

The first RNA recognition motif of the Drosophila SNF protein is an example of an RNA binding protein with multi-specificity. It binds different RNA hairpin loops in spliceosomal U1 or U2 small nuclear RNAs, and only in the latter case requires the auxiliary U2A' protein. Here we investigate its functions by crystal structures of SNF alone and bound to U1 stem-loop II, U2A' or U2 stem-loop IV and U2A', SNF dynamics from NMR spectroscopy, and structure-guided mutagenesis in binding studies. We find that different loop-closing base pairs and a nucleotide exchange at the tips of the loops contribute to differential SNF affinity for the RNAs. U2A' immobilizes SNF and RNA residues to restore U2 stem-loop IV binding affinity, while U1 stem-loop II binding does not require such adjustments. Our findings show how U2A' can modulate RNA specificity of SNF without changing SNF conformation or relying on direct RNA contacts.

Ancestral Haloalkane Dehalogenases Show Robustness and Unique Substrate Specificity

Ancestral sequence reconstruction (ASR) represents a powerful approach for empirical testing structure-function relationships of diverse proteins. We employed ASR to predict sequences of five ancestral haloalkane dehalogenases (HLDs) from the HLD-II subfamily. Genes encoding the inferred ancestral sequences were synthesized and expressed in Escherichia coli, and the resurrected ancestral enzymes (AncHLD1-5) were experimentally characterized. Strikingly, the ancestral HLDs exhibited significantly enhanced thermodynamic stability compared to extant enzymes (ΔT_m up to 24 °C), as well as higher specific activities with preference for short multi-substituted halogenated substrates. Moreover, multivariate statistical analysis revealed a shift in the substrate specificity profiles of AncHLD1 and AncHLD2. This is extremely difficult to achieve by rational protein engineering. The study highlights that ASR is an efficient approach for the development of novel biocatalysts and robust templates for directed evolution.

Robustness of Reconstructed Ancestral Protein Functions to Statistical Uncertainty

Hypotheses about the functions of ancient proteins and the effects of historical mutations on them are often tested using ancestral protein reconstruction (APR)-phylogenetic inference of ancestral sequences followed by synthesis and experimental characterization. Usually, some sequence sites are ambiguously reconstructed, with two or more statistically plausible states. The extent to which the inferred functions and mutational effects are robust to uncertainty about the ancestral sequence has not been studied systematically. To address this issue, we reconstructed ancestral proteins in three domain families that have different functions, architectures, and degrees of uncertainty; we then experimentally characterized the functional robustness of these proteins when uncertainty was incorporated using several approaches, including sampling amino acid states from the posterior distribution at each site and incorporating the alternative amino acid state at every ambiguous site in the sequence into a single "worst plausible case" protein. In every case, qualitative conclusions about the ancestral proteins' functions and the effects of key historical mutations were robust to sequence uncertainty, with similar functions observed even when scores of alternate amino acids were incorporated. There was some variation in quantitative descriptors of function among plausible sequences, suggesting that experimentally characterizing robustness is particularly important when quantitative estimates of ancient biochemical parameters are desired. The worst plausible case method appears to provide an efficient strategy for characterizing the functional robustness of ancestral proteins to large amounts of sequence uncertainty. Sampling from the posterior distribution sometimes produced artifactually nonfunctional proteins for sequences reconstructed with substantial ambiguity.

Urb-RIP - An Adaptable and Efficient Approach for Immunoprecipitation of RNAs and Associated RNAs/Proteins

Post-transcriptional regulation of gene expression is an important process that is mediated by interactions between mRNAs and RNA binding proteins (RBP), non-coding RNAs (ncRNA) or ribonucleoproteins (RNP). Key to the study of post-transcriptional regulation of mRNAs and the function of ncRNAs such as long non-coding RNAs (lncRNAs) is an understanding of what factors are interacting with these transcripts. While several techniques exist for the enrichment of a transcript whether it is an mRNA or an ncRNA, many of these techniques are cumbersome or limited in their application. Here we present a novel method for the immunoprecipitation of mRNAs and ncRNAs, Urb—RNA immunoprecipitation (Urb-RIP). This method employs the RRM1 domain of the “resurrected” snRNA-binding protein Urb to enrich messages containing a stem-loop tag. Unlike techniques which employ the MS2 protein, which require large repeats of the MS2 binding element, Urb-RIP requires only one stem-loop. This method routinely provides over ~100-fold enrichment of tagged messages. Using this technique we have shown enrichment of tagged mRNAs and lncRNAs as well as miRNAs and RNA-binding proteins bound to those messages. We have confirmed, using Urb-RIP, interaction between RNA PolIII transcribed lncRNA BC200 and polyA binding protein.

Resurrecting ancestral structural dynamics of an antiviral immune receptor: adaptive binding pocket reorganization repeatedly shifts RNA preference

Background

Although resurrecting ancestral proteins is a powerful tool for understanding the molecular-functional evolution of gene families, nearly all studies have examined proteins functioning in relatively stable biological processes. The extent to which more dynamic systems obey the same ‘rules’ governing stable processes is unclear. Here we present the first detailed investigation of the functional evolution of the RIG-like receptors (RLRs), a family of innate immune receptors that detect viral RNA in the cytoplasm.

Results

Using kinetic binding assays and molecular dynamics simulations of ancestral proteins, we demonstrate how a small number of adaptive protein-coding changes repeatedly shifted the RNA preference of RLRs throughout animal evolution by reorganizing the shape and electrostatic distribution across the RNA binding pocket, altering the hydrogen bond network between the RLR and its RNA target. In contrast to observations of proteins involved in metabolism and development, we find that RLR-RNA preference ‘flip flopped’ between two functional states, and shifts in RNA preference were not always coupled to gene duplications or speciation events. We demonstrate at least one reversion of RLR-RNA preference from a derived to an ancestral function through a novel structural mechanism, indicating multiple structural implementations of similar functions.

Conclusions

Our results suggest a model in which frequent shifts in selection pressures imposed by an evolutionary arms race preclude the long-term functional optimization observed in stable biological systems. As a result, the evolutionary dynamics of immune receptors may be less constrained by structural epistasis and historical contingency.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0818-6) contains supplementary material, which is available to authorized users.

Co-evolution of SNF spliceosomal proteins with their RNA targets in trans-splicing nematodes

Although the mechanism of pre-mRNA splicing has been well characterized, the evolution of spliceosomal proteins is poorly understood. The U1A/U2B″/SNF family (hereafter referred to as the SNF family) of RNA binding spliceosomal proteins participates in both the U1 and U2 small interacting nuclear ribonucleoproteins (snRNPs). The highly constrained nature of this system has inhibited an analysis of co-evolutionary trends between the proteins and their RNA binding targets. Here we report accelerated sequence evolution in the SNF protein family in Phylum Nematoda, which has allowed an analysis of protein:RNA co-evolution. In a comparison of SNF genes from ecdysozoan species, we found a correlation between trans-splicing species (nematodes) and increased phylogenetic branch lengths of the SNF protein family, with respect to their sister clade Arthropoda. In particular, we found that nematodes (~70-80 % of pre-mRNAs are trans-spliced) have experienced higher rates of SNF sequence evolution than arthropods (predominantly cis-spliced) at both the nucleotide and amino acid levels. Interestingly, this increased evolutionary rate correlates with the reliance on trans-splicing by nematodes, which would alter the role of the SNF family of spliceosomal proteins. We mapped amino acid substitutions to functionally important regions of the SNF protein, specifically to sites that are predicted to disrupt protein:RNA and protein:protein interactions. Finally, we investigated SNF's RNA targets: the U1 and U2 snRNAs. Both are more divergent in nematodes than arthropods, suggesting the RNAs have co-evolved with SNF in order to maintain the necessarily high affinity interaction that has been characterized in other species.

A compare-and-contrast NMR dynamics study of two related RRMs: U1A and SNF

The U1A/U2B″/SNF family of small nuclear ribonucleoproteins uses a phylogenetically conserved RNA recognition motif (RRM1) to bind RNA stemloops in U1 and/or U2 small nuclear RNA (snRNA). RRMs are characterized by their α/β sandwich topology, and these RRMs use their β-sheet as the RNA binding surface. Unique to this RRM family is the tyrosine-glutamine-phenylalanine (YQF) triad of solvent-exposed residues that are displayed on the β-sheet surface; the aromatic residues form a platform for RNA nucleobases to stack. U1A, U2B″, and SNF have very different patterns of RNA binding affinity and specificity, however, so here we ask how YQF in Drosophila SNF RRM1 contributes to RNA binding, as well as to domain stability and dynamics. Thermodynamic double-mutant cycles using tyrosine and phenylalanine substitutions probe the communication between those two residues in the free and bound states of the RRM. NMR experiments follow corresponding changes in the glutamine side-chain amide in both U1A and SNF, providing a physical picture of the RRM1 β-sheet surface. NMR relaxation and dispersion experiments compare fast (picosecond to nanosecond) and intermediate (microsecond-to-millisecond) dynamics of U1A and SNF RRM1. We conclude that there is a network of amino acid interactions involving Tyr-Gln-Phe in both SNF and U1A RRM1, but whereas mutations of the Tyr-Gln-Phe triad result in small local responses in U1A, they produce extensive microsecond-to-millisecond global motions throughout SNF that alter the conformational states of the RRM.

Climbing the vertebrate branch of U1A/U2B″ protein evolution

In the vertebrate lineage of the U1A/U2B″/SNF protein family, the U1A and U2B″ proteins bind to RNA stem-loops in the U1 or U2 snRNPs, respectively. However, their specialization is fairly recent, as they evolved from a single ancestral protein. The progress of their specialization (subfunctionalization) can be monitored by the amino acid sequence changes that give rise to their modern RNA-binding specificity. Using ancestral sequence reconstruction to predict the intermediates on the evolutionary branch, a probable path of sequential changes is defined for U1A and U2B″. The RNA-binding affinity for U1A/U2B″ protein ancestors was measured using modern U1 and U2 snRNA stem-loops and RNA stem-loop variants to understand how the proteins' RNA specificities evolved.

Binding affinity and cooperativity control U2B″/snRNA/U2A' RNP formation

The U1A and U2B″ proteins are components of the U1 and U2 snRNPs, respectively, where they bind to snRNA stemloops. While localization of U1A and U2B″ to their respective snRNP is a well-known phenomenon, binding of U2B″ to U2 snRNA is typically thought to be accompanied by the U2A′ protein. The molecular mechanisms that lead to formation of the RNA/U2B″/U2A′ complex and its localization to the U2 snRNP are investigated here, using a combination of in vitro RNA–protein and protein–protein fluorescence and isothermal titration calorimetry binding experiments. We find that U2A′ protein binds to U2B″ with nanomolar affinity but binds to U1A with only micromolar affinity. In addition, there is RNA-dependent cooperativity (linkage) between protein–protein and protein–RNA binding. The unique combination of tight binding and cooperativity ensures that the U2A′/U2B″ complex is partitioned only to the U2 snRNP.

Linkage and allostery in snRNP protein/RNA complexes

Drosophila SNF is a member of the U1A/U2B″/SNF protein family that is found in U1 and U2 snRNPs, where it binds to Stemloop II and Stemloop IV of U1 and U2 snRNA, respectively. SNF also binds to the U2A' protein, but only in the U2 snRNP. Although previous reports have implicated U2A' as a necessary auxiliary protein for the binding of SNF to Stemloop IV, there are no mechanisms that explain the partitioning of U2A' to the U2 snRNP and its absence from the U1 snRNP. Using in vitro RNA binding isotherms and isothermal titration calorimetry, the thermodynamics of SNF/RNA/U2A' ternary complex formation have now been characterized. There is a very large binding cooperativity unique to Stemloop IV that favors formation of the SLIV/SNF/U2A' complex. The binding cooperativity, or heterotropic linkage, is interpreted with respect to linked conformational equilibria of both SNF and its RNA ligand and so represents an example of protein-RNA allostery.