Distinctive properties of signal sequences from bacterial lipoproteins

We have compared a number of attributes (hydrophobicity, amino acid size, charge and secondary structure propensities) of signal sequences from bacterial lipoproteins with the same attributes of signal peptides from other prokaryotic proteins (non-lipoproteins). Lipoprotein leader sequences tend to be shorter, more hydrophobic and bulky, and they have stronger conformational preferences, the most conspicuous being a predicted beta-turn comprising positions 2 or 3 of the mature protein. Another distinctive feature is a maximum in the local energy profile between positions -1 and +2. With one exception (beta-lactamase III), the lipoproteins do not have Pro in their signal peptides, and they tend to have fewer Ser and Thr but more Gly than non-lipoproteins. Lipoproteins also lack a net negative charge in the N-terminal regions of the mature proteins. The signal peptides of the bacteriocin plasmid-coded lysis proteins appear to be unique in that they have all the ascribed features of lipoprotein signals; these characteristics can be used to guide signal peptide mutagenesis experiments and to construct new secretion vehicles.     

Enhanced Synonymous Site Divergence in Positively Selected VertebrateAntimicrobial Peptide Genes

Nonrandom patterns associated with adaptively evolving genes can shed light on how selection and mutation produce rapid changes in sequences. I examine such patterns in two independent families of antimicrobial peptide genes: those in frogs, which are known to have evolved under positive selection, and those in flatfishes, which I show have also evolved under positive selection. I address two recently proposed hypotheses about the molecular evolution of antimicrobial peptide genes. The first is that the mature peptide region is replicated by an error-prone polymerase that increases the mutation rate and the transversion/transition ratio compared to the signal sequence of the same genes. The second is that mature peptides evolve in a coordinated fashion with their propieces, such that a change in net charge in one molecular region prompts an opposite change in charge in the other region. I test these hypotheses using alternative methods that minimize alignment errors, correct for phylogenetic nonindependence, reduce sequence saturation, and account for differing selection pressures on different regions of the gene. In both gene families I show that divergence at both synonymous and nonsynonymous sites within the mature peptide region is enhanced. However, in neither gene family is there evidence of an increased mutational transversion/transition ratio or coordinated evolution. My observations are consistent with either an elevated mutation rate in an adaptively evolving gene region or widespread selection on “silent” sites. These hypotheses challenge the assumption that mutations are random and can be measured by the synonymous substitution rate. [Reviewing Editor: Dr. Willie J. Swanson]     

Atlantic Cod Piscidin and Its Diversification through Positive Selection

Piscidins constitute a family of cationic antimicrobial peptides that are thought to play an important role in the innate immune response of teleosts. On the one hand they show a remarkable diversity, which indicates that they are shaped by positive selection, but on the other hand they are ancient and have specific targets, suggesting that they are constrained by purifying selection. Until now piscidins had only been found in fish species from the superorder Acanthopterygii but we have recently identified a piscidin gene in Atlantic cod (Gadus morhua), thus showing that these antimicrobial peptides are not restricted to evolutionarily modern teleosts. Nucleotide diversity was much higher in the regions of the piscidin gene that code for the mature peptide and its pro domain than in the signal peptide. Maximum likelihood analyses with different evolution models revealed that the piscidin gene is under positive selection. Charge or hydrophobicity-changing amino acid substitutions observed in positively selected sites within the mature peptide influence its amphipathic structure and can have a marked effect on its function. This diversification might be associated with adaptation to new habitats or rapidly evolving pathogens.     

Towards a more nuanced understanding of the relationship between sex-biased gene expression and rates of protein-coding sequence evolution

Genes that are differentially expressed between the sexes (sex-biased genes) are among the fastest evolving genes in animal genomes. The majority of sex-biased expression is attributable to genes that are primarily expressed in sex-limited reproductive tissues, and these reproductive genes are often rapidly evolving because of intra- and intersexual selection pressures. Additionally, studies of multiple taxa have revealed that genes with sex-biased expression are also expressed in a limited number of tissues. This is worth noting because narrowly expressed genes are known to evolve faster than broadly expressed genes. Therefore, it is not clear whether sex-biased genes are rapidly evolving because they have sexually dimorphic expression, because they are expressed in sex-limited reproductive tissues, or because they are narrowly expressed. To determine the extend to which other confounding variables can explain the rapid evolution of sex-biased genes, I analyzed the rates of evolution of sex-biased genes in Drosophila melanogaster and Mus musculus in light of tissue-specific measures of expression. I find that genes with sex-biased expression in somatic tissues shared by both sexes are often evolving faster than non-sex-biased genes, but this is best explained by the narrow expression profiles of sex-biased genes. Sex-biased genes in sex-limited tissues in D. melanogaster, however, evolve faster than other narrowly expressed genes. Therefore, the rapid evolution of sex-biased genes is limited only to those genes primarily expressed in sex-limited reproductive tissues.     

Molecular evolution of imprinted genes: no evidence for antagonistic coevolution.

Genomically imprinted genes are those for which expression is dependent on the sex of the parent from which they are derived. Numerous theories have been proposed for the evolution of genomic imprinting: one theory is that it is an intra-individual manifestation of classical parent -offspring conflict. This theory is unique in predicting that an arms race may develop between maternally and paternally derived genes for the control of foetal growth demands. Such antagonistic coevolution may be mediated through changes in the structure of the proteins concerned. Comparable coevolution is the most likely explanation for the rapid changes seen in antigenic components of parasites and antigen recognition components of immune systems. We have examined the evolution of insulin-like growth factor Igf2, and its antagonistic receptor Igf2r) and find that in contrast to immune genes, at the sites of mutual binding they are highly conserved. In addition, we have analysed the rate of molecular evolution of seven imprinted genes including Igf2 and Igf2r), sequenced in both mouse and rat, and had that this is the same as that of nonimprinted receptors and significantly lower than that of immune genes controlling for differences in mutation rates. Contrary to the expectations of the conflict hypothesis, we hence find no evidence for antagonistic coevolution of imprinted genes mediated by changes in sequence.     

Adaptive radiation of venomous marine snail lineages and the accelerated evolution of venom peptide genes

An impressive biodiversity (>10,000 species) of marine snails (suborder Toxoglossa or superfamily Conoidea) have complex venoms, each containing approximately 100 biologically active, disulfide-rich peptides. In the genus Conus, the most intensively investigated toxoglossan lineage (～500 species), a small set of venom gene superfamilies undergo rapid sequence hyperdiversification within their mature toxin regions. Each major lineage of Toxoglossa has its own distinct set of venom gene superfamilies. Two recently identified venom gene superfamilies are expressed in the large Turridae clade, but not in Conus. Thus, as major venomous molluscan clades expand, a small set of lineage-specific venom gene superfamilies undergo accelerated evolution. The juxtaposition of extremely conserved signal sequences with hypervariable mature peptide regions is unprecedented and raises the possibility that in these gene superfamilies, the signal sequences are conserved as a result of an essential role they play in enabling rapid sequence evolution of the region of the gene that encodes the active toxin.     

cDNA cloning of two novel T-superfamily conotoxins from Conus leopardus

The full-length cDNAs of two novel T-superfamily conotoxins,Lp5.1 and Lp5.2,were clonedfrom a vermivorous cone snail Conus leopardus using 3'/5'-rapid amplification of cDNA ends.The cDNA ofLp5.1 encodes a precursor of 65 residues,including a 22-residue signal peptide,a 28-residue propeptide anda 15-residue mature peptide.Lp5.1 is processed at the common signal site -X-Arg- immediately before themature peptide sequences.In the case of Lp5.2,the precursor includes a 25-residue signal peptide anda 43-residue sequence comprising the propeptide and mature peptide,which is probably cleaved to yield a29-residue propeptide and a 14-residue mature toxin.Although these two conotoxins share a similar signalsequence and a conserved disulfide pattern with the known T-superfamily,the pro-region and mature peptidesare of low identity,especially Lp5.2 with an identity as low as 10.7% compared with the reference Mr5.1a.The elucidated cDNAs of these two toxins will facilitate a better understanding of the species distribution,the sequence diversity of T-superfamily conotoxins,the special gene structure and the evolution of thesepeptides.     

Low-complexity sequences and single amino acid repeats: not just "junk" peptide sequences

For decades proteins were thought to interact in a "lock and key" system, which led to the definition of a paradigm linking stable three-dimensional structure to biological function. As a consequence, any non-structured peptide was considered to be nonfunctional and to evolve neutrally. Surprisingly, the most commonly shared peptides between eukaryotic proteomes are low-complexity sequences that in most conditions do not present a stable three-dimensional structure. However, because these sequences evolve rapidly and because the size variation of a few of them can have deleterious effects, low-complexity sequences have been suggested to be the target of selection. Here we review evidence that supports the idea that these simple sequences should not be considered just "junk" peptides and that selection drives the evolution of many of them.     

Molecular diversity and conformity of neurohormonal peptides: clues to an adaptive role in evolution

Peptide regulators are probably the most widely distributed signalling agents in the animal kingdom. They are found in both invertebrates and vertebrates where they are produced in endocrine, neuroendocrine and neuronal tissues. As more of these ubiquitous messengers have been characterized it has become evident that they may be grouped into families with clearly defined relationships. Amino acid sequences of the mature, final product indicate relationships between for example cholecystokinin (CCK) and gastrin. More detailed examination of peptide precursors can give further insights into family trees and in the case of the secretin-vasoactive intestinal polypeptide family result in the identification of a novel co-coded peptide. Such dual coding has led to the hypothesis of gene-duplication in peptide evolution, a phenomenon admirably exemplified by the glucagon family and the opioid family. A further example of peptide diversity is evident when mRNA processing is examined. Here a single gene encoding two (or more) structural sequences can produce multiple mRNAs, each encoding a unique peptide sequence. The mRNA produced varies according to the tissue site. The calcitonin and Tachykinin family are fine examples with different peptides produced in neurones and endocrine tissues. A remarkable conservatism of peptide sequences is found in the insulin family where distinct relationships are evident between insulin, insulin-related growth factors (IGF) and insect prothoracicotrophic hormone. Such relationships are supported by examination of the genes for insulin and IGF. Peptide regulators do not evolve in isolation and it is clear that their receptors are also exposed to adaptive pressures. Receptor classes for the Tachykinin family are well characterized, with receptors being identified as falling into two categories, SP-P type and SP-E type. Similar situations obtain for the opioids. Much of this information is based on mammalian studies, however recent work on gastrin/CCK receptors in a range of vertebrates indicates a phylogenetic diversity between brain and gut receptors.     

Sequence and structure space of RNA-binding peptides

Studies of RNA-binding peptides, and recent combinatorial library experiments in particular, have demonstrated that diverse peptide sequences and structures can be used to recognize specific RNA sites. The identification of large numbers of sequences capable of binding to a particular site has provided extensive phylogenetic information used to deduce basic principles of recognition. The high frequency at which RNA-binding peptides are found in large sequence libraries suggests plausible routes to evolve sequence-specific binders, facilitating the design of new binding molecules and perhaps reflecting characteristics of natural evolution.     

Creating Prebiotic Sanctuary: Self-Assembling Supramolecular Peptide Structures Bind and Stabilize RNA

Any attempt to uncover the origins of life must tackle the known ‘blind watchmaker problem’. That is to demonstrate the likelihood of the emergence of a prebiotic system simple enough to be formed spontaneously and yet complex enough to allow natural selection that will lead to Darwinistic evolution. Studies of short aromatic peptides revealed their ability to self-assemble into ordered and stable structures. The unique physical and chemical characteristics of these peptide assemblies point out to their possible role in the origins of life. We have explored mechanisms by which self-assembling short peptides and RNA fragments could interact together and go through a molecular co-evolution, using diphenylalanine supramolecular assemblies as a model system. The spontaneous formation of these self-assembling peptides under prebiotic conditions, through the salt-induced peptide formation (SIPF) pathway was demonstrated. These peptide assemblies possess the ability to bind and stabilize ribonucleotides in a sequence-depended manner, thus increase their relative fitness. The formation of these peptide assemblies is dependent on the homochirality of the peptide monomers: while homochiral peptides (L-Phe-L-Phe and D-Phe-D-Phe) self-assemble rapidly in aqueous environment, heterochiral diastereoisomers (L-Phe-D-Phe and D-Phe-L-Phe) do not tend to self-assemble. This characteristic consists with the homochirality of all living matter. Finally, based on these findings, we propose a model for the role of short self-assembling peptides in the prebiotic molecular evolution and the origin of life.     

Circular dichroism studies on synthetic peptides corresponding to the cleavage site region of precursor proteins

The conformations of synthetic peptides which span the region in which the precursor part of proteins (signal sequences) destined for export are cleaved by signal peptidases, were investigated by circular dichroism spectroscopy. Pentapeptides comprising amino acids only from the carboxy-terminus of signal sequences or the amino terminus of the mature protein do not have any preferred conformation in a variety of solvents. Octa- and nonapeptides containing amino acids from the carboxy-terminal protion of signal sequences and the amino-terminus of the mature portions of precursor proteins tend to adopt beta-turn conformations in trifluoroethanol and micelles of sodium dodecylsulphate. Hence, in addition to the distribution of amino acids with small side chains at the carboxy terminus of signal sequences, it is conceivable that signal peptidases also recognize a beta-turn conformation in the cleavage site region of precursor proteins.     

Rapid evolution of yeast centromeres in the absence of drive

To find the most rapidly evolving regions in the yeast genome we compared most of chromosome III from three closely related lineages of the wild yeast Saccharomyces paradoxus. Unexpectedly, the centromere appears to be the fastest-evolving part of the chromosome, evolving even faster than DNA sequences unlikely to be under selective constraint (i.e., synonymous sites after correcting for codon usage bias and remnant transposable elements). Centromeres on other chromosomes also show an elevated rate of nucleotide substitution. Rapid centromere evolution has also been reported for some plants and animals and has been attributed to selection for inclusion in the egg or the ovule at female meiosis. But Saccharomyces yeasts have symmetrical meioses with all four products surviving, thus providing no opportunity for meiotic drive. In addition, yeast centromeres show the high levels of polymorphism expected under a neutral model of molecular evolution. We suggest that yeast centromeres suffer an elevated rate of mutation relative to other chromosomal regions and they change through a process of "centromere drift," not drive.     

Alternatively and constitutively spliced exons are subject to different evolutionary forces

There has been a controversy on whether alternatively spliced exons (ASEs) evolve faster than constitutively spliced exons (CSEs). Although it has been noted that ASEs are subject to weaker selective constraints than CSEs, so they evolve faster, there have also been studies that indicated slower evolution in ASEs than in CSEs. In this study, we retrieve more than 5,000 human-mouse orthologous exons and calculate the synonymous (KS) and nonsynonymous (KA) substitution rates in these exons. Our results show that ASEs have higher KA values and higher KA/KS ratios than CSEs, indicating faster amino acid-level evolution in ASEs. The faster evolution may be in part due to weaker selective constraints. It is also possible that the faster rate is in part due to faster functional evolution in ASEs. On the other hand, the majority of ASEs have lower KS values than CSEs. With reference to the substitution rate in introns, we show that the KS values in ASEs are close to the neutral substitution rate, whereas the synonymous substitution rate in CSEs has likely been accelerated. The elevated synonymous rate in CSEs is not related to CpG dinucleotides or low-complexity regions of protein but may be weakly related to codon usage bias. The overall trends of higher KA and lower KS in ASEs than in CSEs are also observed in human-rat and mouse-rat comparisons. Therefore, our observations hold for mammals of different molecular clocks.     

Evidence of positively selected sites in mammalian alpha-defensins

Alpha-defensins are a family of mammalian antimicrobial peptides that exhibit variable activity against a panel of microbes, including bacteria, fungi, and enveloped viruses. We have employed a maximum-likelihood approach to detect evidence of positive selection (adaptive evolution) in the evolution of these important molecules of the innate immune response. We have identified 14 amino acid sites that are predicted to be subject to positive selection. Furthermore, we show that all these sites are located in the mature antimicrobial peptide and not in the prepropeptide region of the molecule, implying that they are of functional importance. These results suggest that mammalian alpha-defensins have been under selective pressure to evolve in response to potentially infectious challenges by fast-evolving microbes.     

Adaptation or biased gene conversion? Extending the null hypothesis of molecular evolution

The analysis of evolutionary rates is a popular approach to characterizing the effect of natural selection at the molecular level. Sequences contributing to species adaptation are expected to evolve faster than nonfunctional sequences because favourable mutations have a higher fixation probability than neutral ones. Such an accelerated rate of evolution might be due to factors other than natural selection, in particular GC-biased gene conversion. This is true of neutral sequences, but also of constrained sequences, which can be illustrated using the mouse Fxy gene. Several criteria can discriminate between the natural selection and biased gene conversion models. These criteria suggest that the recently reported human accelerated regions are most likely the result of biased gene conversion. We argue that these regions, far from contributing to human adaptation, might represent the Achilles' heel of our genome.     

Electrostatics in the ribosomal tunnel modulate chain elongation rates

Electrostatic potentials along the ribosomal exit tunnel are nonuniform and negative. The significance of electrostatics in the tunnel remains relatively uninvestigated, yet they are likely to play a role in translation and secondary folding of nascent peptides. To probe the role of nascent peptide charges in ribosome function, we used a molecular tape measure that was engineered to contain different numbers of charged amino acids localized to known regions of the tunnel and measured chain elongation rates. Positively charged arginine or lysine sequences produce transient arrest (pausing) before the nascent peptide is fully elongated. The rate of conversion from transiently arrested to full-length nascent peptide is faster for peptides containing neutral or negatively charged residues than for those containing positively charged residues. We provide experimental evidence that extraribosomal mechanisms do not account for this charge-specific pausing. We conclude that pausing is due to charge-specific interactions between the tunnel and the nascent peptide.     

Novel alpha-conotoxins identified by gene sequencing from cone snails native to Hainan, and their sequence diversity.

Conotoxins (CTX) from the venom of marine cone snails (genus Conus) represent large families of proteins, which show a similar precursor organization with surprisingly conserved signal sequence of the precursor peptides, but highly diverse pharmacological activities. By using the conserved sequences found within the genes that encode the alpha-conotoxin precursors, a technique based on RT-PCR was used to identify, respectively, two novel peptides (LiC22, LeD2) from the two worm-hunting Conus species Conus lividus, and Conus litteratus, and one novel peptide (TeA21) from the snail-hunting Conus species Conus textile, all native to Hainan in China. The three peptides share an alpha4/7 subfamily alpha-conotoxins common cysteine pattern (CCX(4)CX(7)C, two disulfide bonds), which are competitive antagonists of nicotinic acetylcholine receptor (nAChRs). The cDNA of LiC22N encodes a precursor of 40 residues, including a propeptide of 19 residues and a mature peptide of 21 residues. The cDNA of LeD2N encodes a precursor of 41 residues, including a propeptide of 21 residues and a mature peptide of 16 residues with three additional Gly residues. The cDNA of TeA21N encodes a precursor of 38 residues, including a propeptide of 20 residues and a mature peptide of 17 residues with an additional residue Gly. The additional residue Gly of LeD2N and TeA21N is a prerequisite for the amidation of the preceding C-terminal Cys. All three sequences are processed at the common signal site -X-Arg- immediately before the mature peptide sequences. The properties of the alpha4/7 conotoxins known so far were discussed in detail. Phylogenetic analysis of the new conotoxins in the present study and the published homologue of alpha4/7 conotoxins from the other Conus species were performed systematically. Patterns of sequence divergence for the three regions of signal, proregion, and mature peptides, both nucleotide acids and residue substitutions in DNA and peptide levels, as well as Cys codon usage were analyzed, which suggest how these separate branches originated. Percent identities of the DNA and amino acid sequences of the signal region exhibited high conservation, whereas the sequences of the mature peptides ranged from almost identical to highly divergent between inter- and intra-species. Notably, the diversity of the proregion was also high, with an intermediate percentage of divergence between that observed in the signal and in the toxin regions. The data presented are new and are of importance, and should attract the interest of researchers in this field. The elucidated cDNAs of these toxins will facilitate a better understanding of the relationship of their structure and function, as well as the process of their evolutionary relationships.     

Primary and secondary structure of the pore-forming peptide of pathogenic Entamoeba histolytica.

A pore-forming peptide is implicated in the potent cytolytic activity of pathogenic Entamoeba histolytica. Using NH2-terminal sequence information of this peptide, the corresponding cDNA was isolated. The cDNA-deduced amino acid sequence revealed a putative signal peptide and a mature peptide of 77 amino acids including six cysteine residues. Computer-aided secondary structure analysis predicted that the peptide would be composed of four adjacent alpha-helices, and CD spectroscopy indicated an all alpha-helical conformation. The tertiary structure appears to be stabilized by three disulfide bonds; the pore-forming activity was not sensitive to heat but was lost in the presence of reducing agents. Sequence homology was found to the saposins and to surfactant-associated protein B, both mammalian polypeptides of similar size and secondary structure but of non-lytic function. In particular, the six cysteine residues were found to be conserved, suggesting a common motif for stabilizing a favourable tertiary structure. Compared with previously characterized toxic peptides also containing three disulfide bonds, the amoeba peptide may represent a distinct class of biologically active peptides.