The molecular evolution of signal peptides

High-throughput optimisation of protein secretion in yeast via an engineered biosensor

Secretion of high-value proteins and enzymes is fundamental to the synthetic biology economy, allowing continuous fermentation during production and protein purification without cell lysis. Most eukaryotic protein secretion is encoded by an N-terminal signal peptide (SP); however, the strong impact of SP sequence variation on the secretion efficiency of a given protein is not well defined. Despite high natural SP sequence diversity, most recombinant protein secretion systems use only a few well-characterised SPs. Additionally, the selection of promoters and terminators can significantly affect secretion efficiency, yet screening numerous genetic constructs for optimal sequences remains inefficient. Here, we adapted a yeast G-protein-coupled receptor (GPCR) biosensor, to measure the concentration of a peptide tag that is co-secreted with any protein of interest (POI). Thus, protein secretion efficiency can be quantified via induction of a fluorescent reporter that is upregulated downstream of receptor activation. This enabled high-throughput screening of over 6000 combinations of promoters, SPs, and terminators, assembled using one-pot Combinatorial Golden Gate cloning. We demonstrate this biosensor can quickly identify best combinations for secretion and quantify secretion levels. Our results highlight the importance of SP optimisation as an initial step in designing heterologous protein expression strategies, demonstrating the value of high-throughput screening (HTS) approaches for maximising secretion efficiency.

Taken to extremes: Loss of plastid rpl32 in Streptophyta and Cuscuta's unconventional solution for its replacement

The evolution of plant genomes is riddled with exchanges of genetic material within one plant (endosymbiotic gene transfer/EGT) and between unrelated plants (horizontal gene transfer/HGT). These exchanges have left their marks on plant genomes. Parasitic plants with their special evolutionary niche provide ample examples for these processes because they are under a reduced evolutionary pressure to maintain autotrophy and thus to conserve their plastid genomes. On the other hand, the close physical connections with different hosts enabled them to acquire genetic material from other plants. Based on an analysis of an extensive dataset including the parasite Cuscuta campestris and other parasitic plant species, we identified a unique evolutionary history of rpl32 genes coding for an essential plastid ribosomal subunit in Cuscuta. Our analysis suggests that the gene was most likely sequestered by HGT from a member of the Oxalidales order serving as host to an ancestor of the Cuscuta subgenus Grammica. Oxalidales had suffered an ancestral EGT of rpl32 predating the evolution of the genus Cuscuta. The HGT subsequently relieved the plastid rpl32 from its evolutionary constraint and led to its loss from the plastid genome. The HGT-based acquisition in Cuscuta is supported by a high sequence similarity of the mature L32 protein between species of the subgenus Grammica and representatives of the Oxalidales, and by a surprisingly conserved transit peptide, whose functionality in Cuscuta was experimentally verified. The findings are discussed in view of an overall pattern of EGT events for plastid ribosomal subunits in Streptophyta.

A toxin-based approach to neuropeptide and peptide hormone discovery

Peptide hormones and neuropeptides form a diverse class of bioactive secreted molecules that control essential processes in animals. Despite breakthroughs in peptide discovery, many signaling peptides remain undiscovered. Recently, we demonstrated the use of somatostatin-mimicking toxins from cone snails to identify the invertebrate ortholog of somatostatin. Here, we show that this toxin-based approach can be systematically applied to discover other unknown secretory peptides that are likely to have signaling function. Using large sequencing datasets, we searched for homologies between cone snail toxins and secreted proteins from the snails' prey. We identified and confirmed expression of five toxin families that share strong similarities with unknown secretory peptides from mollusks and annelids and in one case also from ecdysozoans. Based on several lines of evidence we propose that these peptides likely act as signaling peptides that serve important physiological functions. Indeed, we confirmed that one of the identified peptides belongs to the family of crustacean hyperglycemic hormone, a peptide not previously observed in Spiralia. We propose that this discovery pipeline can be broadly applied to other systems in which one organism has evolved molecules to manipulate the physiology of another.

Hold the fold: how delayed folding aids protein secretion

In bacteria, N-terminal signal peptides mark proteins for transport across the plasma membrane. A recent study by Smets et al (2022) followed the folding of a pair of structural twins to shed light on how evolution has optimised the secretory process.

The Molecular Biodiversity of Protein Targeting and Protein Transport Related to the Endoplasmic Reticulum

Looking at the variety of the thousands of different polypeptides that have been focused on in the research on the endoplasmic reticulum from the last five decades taught us one humble lesson: no one size fits all. Cells use an impressive array of components to enable the safe transport of protein cargo from the cytosolic ribosomes to the endoplasmic reticulum. Safety during the transit is warranted by the interplay of cytosolic chaperones, membrane receptors, and protein translocases that together form functional networks and serve as protein targeting and translocation routes. While two targeting routes to the endoplasmic reticulum, SRP (signal recognition particle) and GET (guided entry of tail-anchored proteins), prefer targeting determinants at the N- and C-terminus of the cargo polypeptide, respectively, the recently discovered SND (SRP-independent) route seems to preferentially cater for cargos with non-generic targeting signals that are less hydrophobic or more distant from the termini. With an emphasis on targeting routes and protein translocases, we will discuss those functional networks that drive efficient protein topogenesis and shed light on their redundant and dynamic nature in health and disease.

Intramolecular quality control: HIV-1 envelope gp160 signal-peptide cleavage as a functional folding checkpoint

Removal of the membrane-tethering signal peptides that target secretory proteins to the endoplasmic reticulum is a prerequisite for proper folding. While generally thought to be removed co-translationally, we report two additional post-targeting functions for the HIV-1 gp120 signal peptide, which remains attached until gp120 folding triggers its removal. First, the signal peptide improves folding fidelity by enhancing conformational plasticity of gp120 by driving disulfide isomerization through a redox-active cysteine. Simultaneously, the signal peptide delays folding by tethering the N terminus to the membrane, until assembly with the C terminus. Second, its carefully timed cleavage represents intramolecular quality control and ensures release of (only) natively folded gp120. Postponed cleavage and the redox-active cysteine are both highly conserved and important for viral fitness. Considering the ∼15% proteins with signal peptides and the frequency of N-to-C contacts in protein structures, these regulatory roles of signal peptides are bound to be more common in secretory-protein biogenesis.

A modular two yeast species secretion system for the production and preparative application of unspecific peroxygenases

Fungal unspecific peroxygenases (UPOs) represent an enzyme class catalysing versatile oxyfunctionalisation reactions on a broad substrate scope. They are occurring as secreted, glycosylated proteins bearing a haem-thiolate active site and rely on hydrogen peroxide as the oxygen source. However, their heterologous production in a fast-growing organism suitable for high throughput screening has only succeeded once-enabled by an intensive directed evolution campaign. We developed and applied a modular Golden Gate-based secretion system, allowing the first production of four active UPOs in yeast, their one-step purification and application in an enantioselective conversion on a preparative scale. The Golden Gate setup was designed to be universally applicable and consists of the three module types: i) signal peptides for secretion, ii) UPO genes, and iii) protein tags for purification and split-GFP detection. The modular episomal system is suitable for use in Saccharomyces cerevisiae and was transferred to episomal and chromosomally integrated expression cassettes in Pichia pastoris. Shake flask productions in Pichia pastoris yielded up to 24 mg/L secreted UPO enzyme, which was employed for the preparative scale conversion of a phenethylamine derivative reaching 98.6 % ee. Our results demonstrate a rapid, modular yeast secretion workflow of UPOs yielding preparative scale enantioselective biotransformations.

Drosophila female reproductive tract gene expression reveals coordinated mating responses and rapidly evolving tissue-specific genes

Sexual reproduction in internally fertilizing species requires complex coordination between female and male reproductive systems and among the diverse tissues of the female reproductive tract (FRT). Here, we report a comprehensive, tissue-specific investigation of Drosophila melanogaster FRT gene expression before and after mating. We identified expression profiles that distinguished each tissue, including major differences between tissues with glandular or primarily nonglandular epithelium. All tissues were enriched for distinct sets of genes possessing secretion signals that exhibited accelerated evolution, as might be expected for genes participating in molecular interactions between the sexes within the FRT extracellular environment. Despite robust transcriptional differences between tissues, postmating responses were dominated by coordinated transient changes indicative of an integrated systems-level functional response. This comprehensive characterization of gene expression throughout the FRT identifies putative female contributions to postcopulatory events critical to reproduction and potentially reproductive isolation, as well as the putative targets of sexual selection and conflict.

Facets of Autophagy Based Unconventional Protein Secretion-The Road Less Traveled

Unconventional protein secretion (UCPS) of leaderless proteins bypasses the conventional endoplasmic reticulum (ER)-Golgi route. The proportion of UCPS in the secretome varies tremendously across eukaryotes. Interestingly, macroautophagy, an intracellular recycling process that is generally involved in cargo degradation, also participates in UCPS. This emerging field of secretory mode of autophagy is underexplored and has several unanswered questions regarding the composition of players, cargo, and the mechanisms that drive it. As secretomes vary considerably across cell types and physiological conditions, the contribution of secretory autophagy in healthy and pathophysiological states remain to be elucidated. Recent studies have begun to shed light on this enigmatic process.

Plastid transit peptides-where do they come from and where do they all belong? Multi-genome and pan-genomic assessment of chloroplast transit peptide evolution

Subcellular relocalization of proteins determines an organism's metabolic repertoire and thereby its survival in unique evolutionary niches. In plants, the plastid and its various morphotypes import a large and varied number of nuclear-encoded proteins to orchestrate vital biochemical reactions in a spatiotemporal context. Recent comparative genomics analysis and high-throughput shotgun proteomics data indicate that there are a large number of plastid-targeted proteins that are either semi-conserved or non-conserved across different lineages. This implies that homologs are differentially targeted across different species, which is feasible only if proteins have gained or lost plastid targeting peptides during evolution. In this study, a broad, multi-genome analysis of 15 phylogenetically diverse genera and in-depth analyses of pangenomes from Arabidopsis and Brachypodium were performed to address the question of how proteins acquire or lose plastid targeting peptides. The analysis revealed that random insertions or deletions were the dominant mechanism by which novel transit peptides are gained by proteins. While gene duplication was not a strict requirement for the acquisition of novel subcellular targeting, 40% of novel plastid-targeted genes were found to be most closely related to a sequence within the same genome, and of these, 30.5% resulted from alternative transcription or translation initiation sites. Interestingly, analysis of the distribution of amino acids in the transit peptides of known and predicted chloroplast-targeted proteins revealed monocot and eudicot-specific preferences in residue distribution.

Signal peptide recognition in Trypanosoma cruzi GP82 adhesin relies on its localization at protein N-terminus

Trypanosoma cruzi, the causative agent of Chagas disease, has a dense coat of GPI-anchored virulence factors. T. cruzi GPI-anchored adhesin GP82 is encoded by a repertoire of transcripts containing several in-frame initiation codons located up-stream from that adjacent to the predicted signal peptide (SP). Transfection of T. cruzi epimastigotes with constructs encoding GP82 starting at the SP or from the farthest up-stream methionine confirmed protein expression on the parasite cell surface, comparable to the native GP82. Proteins were fully functional, inducing parasite adhesion to HeLa cells and lysosome mobilization, events required for parasite invasion. Transgenic and native GP82 proteins showed indistinguishable electrophoretic mobility, suggesting similar processing of the SP. Deletion of SP generated a ~72 kDa protein devoid of N-linked oligosaccharides allowing irrefutable identification of GP82 precursor. SP transposition to an internal region of GP82 rendered the signal unrecognizable by the signal peptidase and incapable to direct the nascent protein for ER-membrane association. Altogether our data strongly suggests that GP82 SP fails to function as transmembrane domain and its recognition by the signal peptidase shows strict dependence on the signal localization at protein N-terminus. This report presents the first experimental characterization of the full-length GP82 and its signal peptide.

NommPred: Prediction of Mitochondrial and Mitochondrion-Related Organelle Proteins of Nonmodel Organisms

To estimate the functions of mitochondria of diverse eukaryotic nonmodel organisms in which the mitochondrial proteomes are not available, it is necessary to predict the protein sequence features of the mitochondrial proteins computationally. Various prediction methods that are trained using the proteins of model organisms belonging particularly to animals, plants, and fungi exist. However, such methods may not be suitable for predicting the proteins derived from nonmodel organisms because the sequence features of the mitochondrial proteins of diversified nonmodel organisms can differ from those of model organisms that are present only in restricted parts of the tree of eukaryotes. Here, we proposed NommPred, which predicts the mitochondrial proteins of nonmodel organisms that are widely distributed over eukaryotes. We used a gradient boosting machine to develop 2 predictors-one for predicting the proteins of mitochondria and the other for predicting the proteins of mitochondrion-related organelles that are highly reduced mitochondria. The performance of both predictors was found to be better than that of the best method available.

Evolutionary Interplay between Symbiotic Relationships and Patterns of Signal Peptide Gain and Loss

Can orthologous proteins differ in terms of their ability to be secreted? To answer this question, we investigated the distribution of signal peptides within the orthologous groups of Enterobacterales. Parsimony analysis and sequence comparisons revealed a large number of signal peptide gain and loss events, in which signal peptides emerge or disappear in the course of evolution. Signal peptide losses prevail over gains, an effect which is especially pronounced in the transition from the free-living or commensal to the endosymbiotic lifestyle. The disproportionate decline in the number of signal peptide-containing proteins in endosymbionts cannot be explained by the overall reduction of their genomes. Signal peptides can be gained and lost either by acquisition/elimination of the corresponding N-terminal regions or by gradual accumulation of mutations. The evolutionary dynamics of signal peptides in bacterial proteins represents a powerful mechanism of functional diversification.

Synonymous Mutations at the Beginning of the Influenza A Virus Hemagglutinin Gene Impact Experimental Fitness

The fitness effects of synonymous mutations can provide insights into biological and evolutionary mechanisms. We analyzed the experimental fitness effects of all single-nucleotide mutations, including synonymous substitutions, at the beginning of the influenza A virus hemagglutinin (HA) gene. Many synonymous substitutions were deleterious both in bulk competition and for individually isolated clones. Investigating protein and RNA levels of a subset of individually expressed HA variants revealed that multiple biochemical properties contribute to the observed experimental fitness effects. Our results indicate that a structural element in the HA segment viral RNA may influence fitness. Examination of naturally evolved sequences in human hosts indicates a preference for the unfolded state of this structural element compared to that found in swine hosts. Our overall results reveal that synonymous mutations may have greater fitness consequences than indicated by simple models of sequence conservation, and we discuss the implications of this finding for commonly used evolutionary tests and analyses.

Antagonistic Coevolution of MER Tyrosine Kinase Expression and Function

TYRO3, AXL, and MERTK (TAM) receptors are a family of receptor tyrosine kinases that maintain homeostasis through the clearance of apoptotic cells, and when defective, contribute to chronic inflammatory and autoimmune diseases such as atherosclerosis, multiple sclerosis, systemic lupus erythematosus, rheumatoid arthritis, and Crohn's disease. In addition, certain enveloped viruses utilize TAM receptors for immune evasion and entry into host cells, with several viruses preferentially hijacking MERTK for these purposes. Despite the biological importance of TAM receptors, little is understood of their recent evolution and its impact on their function. Using evolutionary analysis of primate TAM receptor sequences, we identified strong, recent positive selection in MERTK's signal peptide and transmembrane domain that was absent from TYRO3 and AXL. Reconstruction of hominid and primate ancestral MERTK sequences revealed three nonsynonymous single nucleotide polymorphisms in the human MERTK signal peptide, with a G14C mutation resulting in a predicted non-B DNA cruciform motif, producing a significant decrease in MERTK expression with no significant effect on MERTK trafficking or half-life. Reconstruction of MERTK's transmembrane domain identified three amino acid substitutions and four amino acid insertions in humans, which led to significantly higher levels of self-clustering through the creation of a new interaction motif. This clustering counteracted the effect of the signal peptide mutations through enhancing MERTK avidity, whereas the lower MERTK expression led to reduced binding of Ebola virus-like particles. The decreased MERTK expression counterbalanced by increased avidity is consistent with antagonistic coevolution to evade viral hijacking of MERTK.

Proteomic approaches to decipher cancer cell secretome

In this review, we give an overview of the actual proteomic approaches used in the study of cancer cells secretome. In particular, we describe the proteomic strategies to decipher cancer cell secretome initially focusing on the different aspects of sample preparation. We examine the issues related to the presence of low abundant proteins, the analysis of secreted proteins in the conditioned media with or without the removal of fetal bovine serum and strategies developed to reduce intracellular protein contamination. As regards the identification and quantification of secreted proteins, we described the different proteomic approaches used, i.e. gel-based, MS-based (label-based and label-free), and the antibody and array-based methods, together with some of the most recent applications in the field of cancer research. Moreover, we describe the bioinformatics tools developed for the in silico validation and characterization of cancer cells secretome. We also discuss the most important available tools for protein annotation and for prediction of classical and non-classical secreted proteins. In summary in this review advances, concerns and challenges in the field of cancer secretome analysis are discussed.

Distinct Structural Elements Govern the Folding, Stability, and Catalysis in the Outer Membrane Enzyme PagP

The outer membrane enzyme PagP is indispensable for lipid A palmitoylation in Gram-negative bacteria and has been implicated in resistance to host immune defenses. PagP possesses an unusual structure for an integral membrane protein, with a highly dynamic barrel domain that is tilted with respect to the membrane normal. In addition, it contains an N-terminal amphipathic helix. Recent functional and structural studies have shown that these molecular factors are critical for PagP to carry out its function in the challenging environment of the bacterial outer membrane. However, the precise contributions of the N-helix to folding and stability and residues that can influence catalytic rates remain to be addressed. Here, we identify a sequence-dependent stabilizing role for the N-terminal helix of PagP in the measured thermodynamic stability of the barrel. Using chimeric barrel sequences, we show that the Escherichia coli PagP N-terminal helix confers 2-fold greater stability to the Salmonella typhimurium barrel. Further, we find that the W78F substitution in S. typhimurium causes a nearly 20-fold increase in the specific activity in vitro for the phospholipase reaction, compared to that of E. coli PagP. Here, phenylalanine serves as a key regulator of catalysis, possibly by increasing the reaction rate. Through coevolution analysis, we detect an interaction network between seemingly unrelated segments of this membrane protein. Exchanging the structural and functional features between homologous PagP enzymes from E. coli and S. typhimurium has provided us with an understanding of the molecular factors governing PagP stability and function.

Prolonged Adaptive Evolution of a Defensive Gene in the Solanaceae

Although plants and their natural enemies may coevolve for prolonged periods, little is known about how long individual plant defensive genes are involved in the coevolutionary process. We address this issue by examining patterns of selection on the defensive gene threonine deaminase (TD). Tomato (Solanum lycopersicum) has two copies of this gene. One performs the canonical housekeeping function in amino acid metabolism of catalyzing the first reaction in the conversion of threonine to isoleucine. The second copy functions as an antinutritive defense against lepidopteran herbivores by depleting threonine in the insect gut. Wild tobacco (Nicotiana attenuata) also contains a defensive copy. We show that a single copy of TD underwent two or three duplications near the base of the Solanaceae. One copy retains the housekeeping function, whereas a second copy evolved defensive functions. Positive selection occurred on the branch of the TD2 gene tree subtending the common ancestor of the Nicotianoideae and Solanoideae. It also occurred within the Solanoideae clade but not within the Nicotianoideae clade. Finally, it occurred on most branches leading from the common ancestor to S. lycopersicum. Based on recent calibrations of the Solanaceae phylogeny, TD2 experienced adaptive substitutions for a period of 30-50 My. We suggest that the most likely explanation for this result is fluctuating herbivore abundances: When herbivores are rare, relaxed selection increases the likelihood that slightly disadvantageous mutations will be fixed by drift; when herbivores are common, increased selection causes the evolution of compensatory adaptive mutations. Alternative explanations are also discussed.

Plus ça change - evolutionary sequence divergence predicts protein subcellular localization signals

BACKGROUND

Protein subcellular localization is a central problem in understanding cell biology and has been the focus of intense research. In order to predict localization from amino acid sequence a myriad of features have been tried: including amino acid composition, sequence similarity, the presence of certain motifs or domains, and many others. Surprisingly, sequence conservation of sorting motifs has not yet been employed, despite its extensive use for tasks such as the prediction of transcription factor binding sites.

RESULTS

Here, we flip the problem around, and present a proof of concept for the idea that the lack of sequence conservation can be a novel feature for localization prediction. We show that for yeast, mammal and plant datasets, evolutionary sequence divergence alone has significant power to identify sequences with N-terminal sorting sequences. Moreover sequence divergence is nearly as effective when computed on automatically defined ortholog sets as on hand curated ones. Unfortunately, sequence divergence did not necessarily increase classification performance when combined with some traditional sequence features such as amino acid composition. However a post-hoc analysis of the proteins in which sequence divergence changes the prediction yielded some proteins with atypical (i.e. not MPP-cleaved) matrix targeting signals as well as a few misannotations.

CONCLUSION

We report the results of the first quantitative study of the effectiveness of evolutionary sequence divergence as a feature for protein subcellular localization prediction. We show that divergence is indeed useful for prediction, but it is not trivial to improve overall accuracy simply by adding this feature to classical sequence features. Nevertheless we argue that sequence divergence is a promising feature and show anecdotal examples in which it succeeds where other features fail.

NMR solution structure and SRP54M predicted interaction of the N-terminal sequence (1-30) of the ovine Doppel protein

Prion protein (PrP(C)) biosynthesis involves a multi-step process that includes translation and post-translational modifications. While PrP has been widely investigated, for the homolog Doppel (Dpl), limited knowledge is available. In this study, we focused on a vital step of eukaryotic protein biosynthesis: targeting by the signal recognition particle (SRP). Taking the ovine Dpl (OvDpl(1-30)) peptide as a template, we studied its behavior in two different hydrophobic environments using CD and NMR spectroscopy. In both trifluoroethanol (TFE) and dihexanoyl-sn-glycero-3-phosphatidylcholine (DHPC), the OvDpl(1-30) peptide revealed to fold in an alpha-helical conformation with a well-defined central region extending from residue Cys8 until Ser22. The NMR structure was subsequently included in a computational docking complex with the conserved M-domain of SRP54 protein (SRP54M), and further compared with the N-terminal structures of mouse Dpl and bovine PrP(C) proteins. This allowed the determination of (i) common predicted N-terminal/SRP54M polar contacts (Asp331, Gln335, Glu365 and Lys432) and (ii) different N-C orientations between prion and Dpl peptides at the SRP54M hydrophobic groove, that are in agreement with each peptide electrostatic potential. Together, these findings provide new insights into the biosynthesis of prion-like proteins. Besides they also show the role of protein conformational switches in signalization toward the endoplasmic membrane, a key event of major significance in the cell cycle. They are thus of general applicability to the study of the biological function of prion-like as well as other proteins.

All roads lead to Rome (but some may be harder to travel): SRP-independent translocation into the endoplasmic reticulum

Translocation into the endoplasmic reticulum (ER) is the first biogenesis step for hundreds of eukaryotic secretome proteins. Over the past 30 years, groundbreaking biochemical, structural and genetic studies have delineated one conserved pathway that enables ER translocation- the signal recognition particle (SRP) pathway. However, it is clear that this is not the only pathway which can mediate ER targeting and insertion. In fact, over the past decade, several SRP-independent pathways have been uncovered, which recognize proteins that cannot engage the SRP and ensure their subsequent translocation into the ER. These SRP-independent pathways face the same challenges that the SRP pathway overcomes: chaperoning the preinserted protein while in the cytosol, targeting it rapidly to the ER surface and generating vectorial movement that inserts the protein into the ER. This review strives to summarize the various mechanisms and machineries which mediate these stages of SRP-independent translocation, as well as examine why SRP-independent translocation is utilized by the cell. This emerging understanding of the various pathways utilized by secretory proteins to insert into the ER draws light to the complexity of the translocational task, and underlines that insertion into the ER might be more varied and tailored than previously appreciated.

Correlated evolution of positions within mammalian cis elements

Transcriptional regulation critically depends on proper interactions between transcription factors (TF) and their cognate DNA binding sites. The widely used model of TF-DNA binding – the Positional Weight Matrix (PWM) – presumes independence between positions within the binding site. However, there is evidence to show that the independence assumption may not always hold, and the extent of interposition dependence is not completely known. We hypothesize that the interposition dependence should partly be manifested as correlated evolution at the positions. We report a Maximum-Likelihood (ML) approach to infer correlated evolution at any two positions within a PWM, based on a multiple alignment of 5 mammalian genomes. Application to a genome-wide set of putative cis elements in human promoters reveals a prevalence of correlated evolution within cis elements. We found that the interdependence between two positions decreases with increasing distance between the positions. The interdependent positions tend to be evolutionarily more constrained and moreover, the dependence patterns are relatively similar across structurally related transcription factors. Although some of the detected mutational dependencies may be due to context-dependent genomic hyper-mutation, notably CG to TG, the majority is likely due to context-dependent preferences for specific nucleotide combinations within the cis elements. Patterns of evolution at individual nucleotide positions within mammalian TF binding sites are often significantly correlated, suggesting interposition dependence. The proposed methodology is also applicable to other classes of non-coding functional elements. A detailed investigation of mutational dependencies within specific motifs could reveal preferred nucleotide combinations that may help refine the DNA binding models.

A structural and functional comparison of nematode and crustacean PDH-like sequences

The elucidation of the whole genome of the nematode Caenorhabditis elegans allowed for the identification of ortholog genes belonging to the pigment dispersing hormone/factor (PDH/PDF) peptide family. Members of this peptide family are known from crustaceans, insects and nematodes and seem to exist exclusively in ecdysozoans where they play a role in different processes, ranging from the dispersion of integumental and eye (retinal) pigments in decapod crustaceans to circadian rhythms in insects and locomotion in C. elegans. Two pdf genes (pdf-1 and pdf-2) encoding three different peptides: PDF-1a, PDF-1b and PDF-2 have been identified in C. elegans. These three C. elegans PDH-like peptides are similar but not identical in primary structure to PDHs from decapod crustaceans. We investigate whether this divergence has an influence on the pigment dispersing function of the peptides in a decapod crustacean, namely the shrimp Palaemon pacificus. We show that C. elegans PDF-1a and b peptides display cross-functional activity by dispersing pigments in the epithelium of P. pacificus at physiological doses. Moreover, by means of a comparative amino acid sequence analysis of nematode and crustacean PDH-like peptides, we can pinpoint several potentially important residues for eliciting pigment dispersing activity in decapod crustaceans. Although there is no sequence information on a receptor for PDH in decapod crustaceans, we postulate that there is general conservation of the PDH/PDF signaling system based on structural similarities of precursor proteins and receptors (including those from a branchiopod crustacean and from C. elegans).

Nuclear export as a key arbiter of "mRNA identity" in eukaryotes

Over the past decade, various studies have indicated that most of the eukaryotic genome is transcribed at some level. The pervasiveness of transcription might seem surprising when one considers that only a quarter of the human genome comprises genes (including exons and introns) and less than 2% codes for protein. This conundrum is partially explained by the unique evolutionary pressures that are imposed on species with small population sizes, such as eukaryotes. These conditions promote the expansion of introns and non-functional intergenic DNA, and the accumulation of cryptic transcriptional start sites. As a result, the eukaryotic gene expression machinery must effectively evaluate whether or not a transcript has all the hallmarks of a protein-coding mRNA. If a transcript contains these features, then positive feedback loops are activated to further stimulate its transcription, processing, nuclear export and ultimately, translation. However if a transcript lacks features associated with "mRNA identity", then the RNA is degraded and/or used to inhibit further transcription and translation of the gene. Here we discuss how mRNA identity is assessed by the nuclear export machinery in order to extract meaningful information from the eukaryotic genome. In the process, we provide an explanation of why certain sequences that are enriched in protein-coding genes, such as the signal sequence coding region, promote mRNA nuclear export in vertebrates. This article is part of a Special Issue entitled: Nuclear Transport and RNA Processing.

Comparative studies of vertebrate Beta integrin genes and proteins: ancient genes in vertebrate evolution

Intregins are heterodimeric α- and β-subunit containing membrane receptor proteins which serve various cell adhesion roles in tissue repair, hemostasis, immune response, embryogenesis and metastasis. At least 18 α- (ITA or ITGA) and 8 β-integrin subunits (ITB or ITGB) are encoded on mammalian genomes. Comparative ITB amino acid sequences and protein structures and ITB gene locations were examined using data from several vertebrate genome projects. Vertebrate ITB genes usually contained 13-16 coding exons and encoded protein subunits with ~800 amino acids, whereas vertebrate ITB4 genes contained 36-39 coding exons and encoded larger proteins with ~1800 amino acids. The ITB sequences exhibited several conserved domains including signal peptide, extracellular β-integrin, β-tail domain and integrin β-cytoplasmic domains. Sequence alignments of the integrin β-cytoplasmic domains revealed highly conserved regions possibly for performing essential functions and its maintenance during vertebrate evolution. With the exception of the human ITB8 sequence, the other ITB sequences shared a predicted 19 residue α-helix for this region. Potential sites for regulating human ITB gene expression were identified which included CpG islands, transcription factor binding sites and microRNA binding sites within the 3'-UTR of human ITB genes. Phylogenetic analyses examined the relationships of vertebrate beta-integrin genes which were consistent with four major groups: 1: ITB1, ITB2, ITB7; 2: ITB3, ITB5, ITB6; 3: ITB4; and 4: ITB8 and a common evolutionary origin from an ancestral gene, prior to the appearance of fish during vertebrate evolution. The phylogenetic analyses revealed that ITB4 is the most likely primordial form of the vertebrate β integrin subunit encoding genes, that is the only β subunit expressed as a constituent of the sole integrin receptor 'α6β4' in the hemidesmosomes of unicellular organisms.

Impact of extracellularity on the evolutionary rate of mammalian proteins

It is of fundamental importance to understand the determinants of the rate of protein evolution. Eukaryotic extracellular proteins are known to evolve faster than intracellular proteins. Although this rate difference appears to be due to the lower essentiality of extracellular proteins than intracellular proteins in yeast, we here show that, in mammals, the impact of extracellularity is independent from the impact of gene essentiality. Our partial correlation analysis indicated that the impact of extracellularity on mammalian protein evolutionary rate is also independent from those of tissue-specificity, expression level, gene compactness, and the number of protein–protein interactions and, surprisingly, is the strongest among all the factors we examined. Similar results were also found from principal component regression analysis. Our findings suggest that different rules govern the pace of protein sequence evolution in mammals and yeasts.

Genomic imprinting - the story of the other half and the conflicts of silencing

Genomic imprinting is an epigenetic mechanism that produces functional differences between the paternal and maternal genomes and plays an essential role in mammalian development and growth. There are a number of genes in our genomes that are subject to genomic imprinting where one parent's copy of the gene is expressed while the other is silent. Silencing of one allele predetermines that any function ascribed to that gene are now dependant on the single active copy. Possession of only a single active allele can lead to deleterious health consequences in humans. If imprinted genes are crucial in mammalian development, one would also expect mutations in these genes to cause diseases. Since imprinting is an epigenetic mechanism, mistakes in maintaining epigenetic mark also cause imprinting disorders. Here we in this review focus on the current understanding of this unique genetic mechanism more than two decades after the first description of the imprinting phenomenon was given by McGrath and Solter. Although the possible molecular mechanisms by which imprinting is imposed and maintained are being identified, we have a long way to go in understanding the molecular mechanisms that regulate the expression of these oddly behaving genes, the function of imprinting and the evolution. Post genomic technologies might ultimately lead to a better understanding of the 'imprinting effects'.

Possible diversifying selection in the imprinted gene, MEDEA, in Arabidopsis

Coevolutionary conflict among imprinted genes that influence traits such as offspring growth may arise when maternal and paternal genomes have different evolutionary optima. This conflict is expected in outcrossing taxa with multiple paternity, but not self-fertilizing taxa. MEDEA (MEA) is an imprinted plant gene that influences seed growth. Disagreement exists regarding the type of selection acting on this gene. We present new data and analyses of sequence diversity of MEA in self-fertilizing and outcrossing Arabidopsis and its relatives, to help clarify the form of selection acting on this gene. Codon-based branch analysis among taxa (PAML) suggests that selection on the coding region is changing over time, and nonsynonymous substitution is elevated in at least one outcrossing branch. Codon-based analysis of diversity within outcrossing Arabidopsis lyrata ssp. petraea (OmegaMap) suggests that diversifying selection is acting on a portion of the gene, to cause elevated nonsynonymous polymorphism. Providing further support for balancing selection in A. lyrata, Hudson, Kreitman and Aguadé analysis indicates that diversity/divergence at silent sites in the MEA promoter and genic region is elevated relative to reference genes, and there are deviations from the neutral frequency spectrum. This combination of positive selection as well as balancing and diversifying selection in outcrossing lineages is consistent with other genes influence by evolutionary conflict, such as disease resistance genes. Consistent with predictions that conflict would be eliminated in self-fertilizing taxa, we found no evidence of positive, balancing, or diversifying selection in A. thaliana promoter or genic region.

Isolation and analysis of the genetic diversity of repertoires of VSG expression site containing telomeres from Trypanosoma brucei gambiense, T. b. brucei and T. equiperdum

Background

African trypanosomes (including Trypanosoma brucei) are unicellular parasites which multiply in the mammalian bloodstream. T. brucei has about twenty telomeric bloodstream form Variant Surface Glycoprotein (VSG) expression sites (BESs), of which one is expressed at a time in a mutually exclusive fashion. BESs are polycistronic transcription units, containing a variety of families of expression site associated genes (ESAGs) in addition to the telomeric VSG. These polymorphic ESAG families are thought to play a role in parasite-host adaptation, and it has been proposed that ESAG diversity might be related to host range. Analysis of the genetic diversity of these telomeric gene families has been confounded by the underrepresentation of telomeric sequences in standard libraries. We have previously developed a method to selectively isolate sets of trypanosome BES containing telomeres using Transformation associated recombination (TAR) cloning in yeast.

Results

Here we describe the isolation of repertoires of BES containing telomeres from three trypanosome subspecies: Trypanosoma brucei gambiense DAL 972 (causative agent of West-African trypanosomiasis), T. b. brucei EATRO 2340 (a nonhuman infective strain) and T. equiperdum STIB 818 (which causes a sexually transmitted disease in equines). We have sequenced and analysed the genetic diversity at four BES loci (BES promoter region, ESAG6, ESAG5 and ESAG2) from these three trypanosome BES repertoires.

Conclusion

With the exception of ESAG2, the BES sequence repertoires derived from T. b. gambiense are both less diverse than and nearly reciprocally monophyletic relative to those from T. b. brucei and T. equiperdum. Furthermore, although we find evidence for adaptive evolution in all three ESAG repertoires in T. b. brucei and T. equiperdum, only ESAG2 appears to be under diversifying selection in T. b. gambiense. This low level of variation in the T. b. gambiense BES sequence repertoires is consistent both with the relatively narrow host range of this subspecies and its apparent long-term clonality. However, our data does not show a clear correlation between size of trypanosome host range and either number of BESs or extent of ESAG genetic diversity.

The concept of translocational regulation

Biological processes are regulated to provide cells with exquisite adaptability to changing environmental conditions and cellular demands. The mechanisms regulating secretory and membrane protein translocation into the endoplasmic reticulum (ER) are unknown. A conceptual framework for translocational regulation is proposed based on our current mechanistic understanding of ER protein translocation and general principles of regulatory control.

The signal sequence coding region promotes nuclear export of mRNA

In eukaryotic cells, most mRNAs are exported from the nucleus by the transcription export (TREX) complex, which is loaded onto mRNAs after their splicing and capping. We have studied in mammalian cells the nuclear export of mRNAs that code for secretory proteins, which are targeted to the endoplasmic reticulum membrane by hydrophobic signal sequences. The mRNAs were injected into the nucleus or synthesized from injected or transfected DNA, and their export was followed by fluorescent in situ hybridization. We made the surprising observation that the signal sequence coding region (SSCR) can serve as a nuclear export signal of an mRNA that lacks an intron or functional cap. Even the export of an intron-containing natural mRNA was enhanced by its SSCR. Like conventional export, the SSCR-dependent pathway required the factor TAP, but depletion of the TREX components had only moderate effects. The SSCR export signal appears to be characterized in vertebrates by a low content of adenines, as demonstrated by genome-wide sequence analysis and by the inhibitory effect of silent adenine mutations in SSCRs. The discovery of an SSCR-mediated pathway explains the previously noted amino acid bias in signal sequences and suggests a link between nuclear export and membrane targeting of mRNAs.

In eukaryotic cells, precursors of messenger RNAs (mRNAs) are synthesized and processed in the nucleus. During processing, noncoding introns are spliced out, and a cap and poly-adenosine sequence are added to the beginning and end of the transcript, respectively. The resulting mature mRNA is exported from the nucleus to the cytoplasm by crossing the nuclear pore. Both the introns and the cap help to recruit factors that are necessary for nuclear export of an mRNA. Here we provide evidence for a novel mRNA export pathway that is specific for transcripts coding for secretory proteins. These proteins contain signal sequences that target them for translocation across the endoplasmic reticulum membrane. We made the surprising observation that the signal sequence coding region (SSCR) can serve as a nuclear export signal of an mRNA that lacks an intron or functional cap. Even the export of an intron-containing natural mRNA was enhanced by its SSCR. The SSCR export signal appears to be characterized in vertebrates by a low content of adenines. Our discovery of an SSCR-mediated pathway explains the previously noted amino acid bias in signal sequences, and suggests a link between nuclear export and membrane targeting of mRNAs.

Signal sequences, which target newly synthesized secretory proteins to the endoplasmic reticulum, are encoded by adenine-depleted nucleotide sequences that promote the nuclear export of the mRNAs.

Globin gene family evolution and functional diversification in annelids

Globins are the most common type of oxygen-binding protein in annelids. In this paper, we show that circulating intracellular globin (Alvinella pompejana and Glycera dibranchiata), noncirculating intracellular globin (Arenicola marina myoglobin) and extracellular globin from various annelids share a similar gene structure, with two conserved introns at canonical positions B12.2 and G7.0. Despite sequence divergence between intracellular and extracellular globins, these data strongly suggest that these three globin types are derived from a common ancestral globin-like gene and evolved by duplication events leading to diversification of globin types and derived functions. A phylogenetic analysis shows a distinct evolutionary history of annelid extracellular hemoglobins with respect to intracellular annelid hemoglobins and mollusc and arthropod extracellular hemoglobins. In addition, dehaloperoxidase (DHP) from the annelid, Amphitrite ornata, surprisingly exhibits close phylogenetic relationships to some annelid intracellular globins. We have characterized the gene structure of A. ornata DHP to confirm assumptions about its homology with globins. It appears that it has the same intron position as in globin genes, suggesting a common ancestry with globins. In A. ornata, DHP may be a derived globin with an unusual enzymatic function.

Substrate-specific translocational attenuation during ER stress defines a pre-emptive quality control pathway

Eukaryotic proteins entering the secretory pathway are translocated into the ER by signal sequences that vary widely in primary structure. We now provide a functional rationale for this long-observed sequence diversity by demonstrating that differences among signals facilitate substrate-selective modulation of protein translocation. We find that during acute ER stress, translocation of secretory and membrane proteins is rapidly and transiently attenuated in a signal sequence-selective manner. Their cotranslational rerouting to the cytosol for degradation reduces the burden of misfolded substrates entering the ER and represents a pathway for pre-emptive quality control (pQC). Bypassing the pQC pathway for the prion protein increases its rate of aggregation in the ER lumen during prolonged stress and renders cells less capable of viable recovery. Conversely, pharmacologically augmenting pQC during ER stress proved protective. Thus, protein translocation is a physiologically regulated process that is utilized for pQC as part of the ER stress response.

An overabundance of phase 0 introns immediately after the start codon in eukaryotic genes

BACKGROUND

A knowledge of the positions of introns in eukaryotic genes is important for understanding the evolution of introns. Despite this, there has been relatively little focus on the distribution of intron positions in genes.

RESULTS

In proteins with signal peptides, there is an overabundance of phase 1 introns around the region of the signal peptide cleavage site. This has been described before. But in proteins without signal peptides, a novel phenomenon is observed: There is a sharp peak of phase 0 intron positions immediately following the start codon, i.e. between codons 1 and 2. This effect is seen in a wide range of eukaryotes: Vertebrates, arthropods, fungi, and flowering plants. Proteins carrying this start codon intron are found to comprise a special class of relatively short, lysine-rich and conserved proteins with an overrepresentation of ribosomal proteins. In addition, there is a peak of phase 0 introns at position 5 in Drosophila genes with signal peptides, predominantly representing cuticle proteins.

CONCLUSION

There is an overabundance of phase 0 introns immediately after the start codon in eukaryotic genes, which has been described before only for human ribosomal proteins. We give a detailed description of these start codon introns and the proteins that contain them.

Molecular evolution of the avian growth hormone gene and comparison with its mammalian counterpart

The molecular evolution of all available avian growth hormone (GH) gene sequences was investigated using both maximum-likelihood and parsimony methods, and the patterns compared to those found in mammals. In contrast to the rapid bursts of evolution observed for mammalian GH, the evolutionary rate of the avian GH mature peptide appears to have been more constant. However several positively selected sites were identified at functionally important positions in the avian signal peptide by the site-specific likelihood method. This implies that sequence variation in the avian GH signal peptide may be adaptive, although more conservative parsimony methods failed to confirm this. Nevertheless, the differing patterns of avian and mammalian GH signal peptide molecular evolution are consistent with the apparently differing roles of GH in controlling growth in these taxonomic groups and support the hypothesis that signal peptide sequence variation may in fact be the basis for increased functional complexity.

A possible role of exon-shuffling in the evolution of signal peptides of human proteins

It was recently shown that there is a predominance of phase 1 introns near the cleavage site of signal peptides encoded by human genes. It was suggested that this biased distribution was due to intron insertion at AGmid R:G proto-splice sites. However, we found that there is no disproportional excess of AGmid R:G that would support insertion at proto-splice sites. In fact, all nGmid R:G sites are enriched in the vicinity of the cleavage site. Additional analyses support an alternative scenario in which exon-shuffling is largely responsible for such excess of phase 1 introns.

Varying signals of the effects of natural selection during teleost growth hormone gene evolution

The growth hormone (GH) gene of teleost fish exhibits a higher degree of variability compared with other vertebrate groups. However, the different selective constraints at the sequence level are not well understood. In this study, maximum-likelihood (ML) models of codon substitutions were used to investigate Darwinian adaptive evolution of the GH gene in teleost fishes. Complete GH gene sequences of 54 fish species were classified into 4 orders, and the variable nature of GH was examined by determining the dNand dSrate variation and the rates of molecular evolution for each teleost order. The results indicate that although the overall evolution rate for teleost GH is high ((1.15 ± 0.01) × 10–9substitutions/(aa site·y)) compared with the "slow phases" in mammals ((0.21 to 0.28 ± 0.05) × 10–9), the vital structure of this gene has been retained. While the majority of the amino acid changes appear to be due to relaxation of purifying selection, some positively selected sites were detected in regions with no specifically identified role in protein function. The positively selected regions observed in salmoniformes lineage suggests a possible role for positive selection driving functional divergence in paralogous forms of the GH gene after whole-genome duplication in this lineage.Key words: teleost fish, growth hormone, positive selection, synonymous substitution, non-synonymous substitution, molecular evolution.

The applicability of recurrent neural networks for biological sequence analysis

Selection of machine learning techniques requires a certain sensitivity to the requirements of the problem. In particular, the problem can be made more tractable by deliberately using algorithms that are biased toward solutions of the requisite kind. In this paper, we argue that recurrent neural networks have a natural bias toward a problem domain of which biological sequence analysis tasks are a subset. We use experiments with synthetic data to illustrate this bias. We then demonstrate that this bias can be exploitable using a data set of protein sequences containing several classes of subcellular localization targeting peptides. The results show that, compared with feed forward, recurrent neural networks will generally perform better on sequence analysis tasks. Furthermore, as the patterns within the sequence become more ambiguous, the choice of specific recurrent architecture becomes more critical.

Prediction of subcellular localization using sequence-biased recurrent networks

MOTIVATION

Targeting peptides direct nascent proteins to their specific subcellular compartment. Knowledge of targeting signals enables informed drug design and reliable annotation of gene products. However, due to the low similarity of such sequences and the dynamical nature of the sorting process, the computational prediction of subcellular localization of proteins is challenging.

RESULTS

We contrast the use of feed forward models as employed by the popular TargetP/SignalP predictors with a sequence-biased recurrent network model. The models are evaluated in terms of performance at the residue level and at the sequence level, and demonstrate that recurrent networks improve the overall prediction performance. Compared to the original results reported for TargetP, an ensemble of the tested models increases the accuracy by 6 and 5% on non-plant and plant data, respectively.

AVAILABILITY

The Protein Prowler incorporating the recurrent network predictor described in this paper is available online at http://pprowler.imb.uq.edu.au/

Widespread adaptive evolution in the human immunodeficiency virus type 1 genome

We investigated variable selective pressures among amino acid sites in HIV-1 genes. Selective pressure at the amino acid level was measured by using the nonsynonymous/synonymous substitution rate ratio (omega = dN/dS). To identify amino acid sites under positive selection with omega > 1, we applied maximum likelihood models that allow variable omega ratios among sites to analyze genomic sequences of 26 HIV-1 lineages including subtypes A, B, and C. Likelihood ratio tests detected sites under positive selection in each of the major genes in the genome: env, gag, pol, vif, and vpr. Positive selection was also detected in nef, tat, and vpu, although those genes are very small. The majority of positive selection sites is located in gp160. Positive selection was not detected if omega was estimated as an average across all sites, indicating the lack of power of the averaging approach. Candidate positive selection sites were mapped onto the available protein tertiary structures and immunogenic epitopes. We measured the physiochemical properties of amino acids and found that those at positive selection sites were more diverse than those at variable sites. Furthermore, amino acid residues at exposed positive selection sites were more physiochemically diverse than at buried positive selection sites. Our results demonstrate genomewide diversifying selection acting on the HIV-1.

Signal sequences control gating of the protein translocation channel in a substrate-specific manner

N-terminal signal sequences mediate targeting of nascent chains to the endoplasmic reticulum and facilitate opening of the protein translocation channel to the passage of substrate. We have assessed each of these steps for a diverse set of mammalian signals. While minimal differences were seen in their targeting function, signal sequences displayed a remarkable degree of variation in initiating nascent chain access to the lumenal environment. Such substrate-specific properties of signals were evolutionarily conserved, functionally matched to their respective mature domains, and important for the proper biogenesis of some proteins. Thus, the sequence variations of signals do not simply represent functional degeneracy, but instead encode critical differences in translocon gating that are coordinated with their respective passengers to facilitate efficient translocation.