An Empirical Codon Model for Protein Sequence Evolution

Mol. Biol. Evol. 2007 24:1464-1479 The first affiliation should have appeared as EMBL-EuropeanBioinformatics Institute, Hinxton, United Kingdom. On page     

Proteotyping：A New Approach Studying Influenza Virus Evolution at the Protein Level

Phylogenetic methods have been widely used to detect the evolution of influenza viruses.However,previous phylogenetic studies of influenza viruses do not make full use of the genetic information at the protein level and therefore cannot distinguish the subtle differences among viral genes.Proteotyping is a new approach to study influenza virus evolution.It aimed at mining the potential genetic information of the viral gene at the protein level by visualizing unique amino acid signatures(proteotypes).Neuraminidase gene fragments of some H5N1 avian influenza viruses were used as an example to illustrate how the proteotyping method worked.Bayesian analysis confirmed that the NA gene tree was mainly divided into three lineages.The NA proteotype analysis further suggested there might be multiple proteotypes within these three lineages and even within single genotypes.At the same time,some proteotypes might even involve more than one genotype.In particular,it also discovered some amino acids of viruses of some genotypes might co-reassort.All these results proved this approach could provide additional information in contrast to results from standard phylogenetic tree analysis.     

Heat Shock Protein 70 Family: Multiple Sequence Comparisons, Function, and Evolution

The heat shock protein 70 kDa sequences (HSP70) are of great importance as molecular chaperones in protein folding and transport. They are abundant under conditions of cellular stress. They are highly conserved in all domains of life: Archaea, eubacteria, eukaryotes, and organelles (mitochondria, chloroplasts). A multiple alignment of a large collection of these sequences was obtained employing our symmetric-iterative ITERALIGN program (Brocchieri and Karlin 1998). Assessments of conservation are interpreted in evolutionary terms and with respect to functional implications. Many archaeal sequences (methanogens and halophiles) tend to align best with the Gram-positive sequences. These two groups also miss a signature segment [about 25 amino acids (aa) long] present in all other HSP70 species (Gupta and Golding 1993). We observed a second signature sequence of about 4 aa absent from all eukaryotic homologues, significantly aligned in all prokaryotic sequences. Consensus sequences were developed for eight groups [Archaea, Gram-positive, proteobacterial Gram-negative, singular bacteria, mitochondria, plastids, eukaryotic endoplasmic reticulum (ER) isoforms, eukaryotic cytoplasmic isoforms]. All group consensus comparisons tend to summarize better the alignments than do the individual sequence comparisons. The global individual consensus ``matches' 87% with the consensus of consensuses sequence. A functional analysis of the global consensus identifies a (new) highly significant mixed charge cluster proximal to the carboxyl terminus of the sequence highlighting the hypercharge run EEDKKRRER (one-letter aa code used). The individual Archaea and Gram-positive sequences contain a corresponding significant mixed charge cluster in the location of the charge cluster of the consensus sequence. In contrast, the four Gram-negative proteobacterial sequences of the alignment do not have a charge cluster (even at the 5% significance level). All eukaryotic HSP70 sequences have the analogous charge cluster. Strikingly, several of the eukaryotic isoforms show multiple mixed charged clusters. These clusters were interpreted with supporting data related to HSP70 activity in facilitating chaperone, transport, and secretion function. We observed that the consensus contains only a single tryptophan residue and a single conserved cysteine. This is interpreted with respect to the target rule for disaggregating misfolded proteins. The mitochondrial HSP70 connections to bacterial HSP70 are analyzed, suggesting a polyphyletic split of Trypanosoma and Leishmania protist mitochondrial (Mt) homologues separated from Mt-animal/fungal/plant homologues. Moreover, the HSP70 sequences from the amitochondrial Entamoeba histolytica and Trichomonas vaginalis species were analyzed. The E. histolytica HSP70 is most similar to the higher eukaryotic cytoplasmic sequences, with significantly weaker alignments to ER sequences and much diminished matching to all eubacterial, mitochondrial, and chloroplast sequences. This appears to be at variance with the hypothesis that E. histolytica rather recently lost its mitochondrial organelle. T. vaginalis contains two HSP70 sequences, one Mt-like and the second similar to eukaryotic cytoplasmic sequences suggesting two diverse origins. Received: 29 January 1998 / Accepted: 14 May 1998     

Protein Sequence Modules

Abstract

Conserved protein sequence segments are commonly believed to correspond to functional sites in the protein sequence. A novel approach is proposed to profile the changing degree of conservation along the protein sequence, by evaluating the occurrence frequencies of all short oligopeptides of the given sequence in a large proteome database. Thus, a protein sequence conservation profile can be plotted for every protein. The profile indicates where along the sequences the potential functional (conserved) sites are located. The corresponding oligopeptides belonging to the sites are very frequent across many prokaryotic species. Analysis of a representative set of such profiles reveals a common feature of all examined proteins: they consist of sequence modules represented by the peaks of conservation. Typical size of the modules (peak-to-peak distance) is 25–30 amino acid residues.     

Messenger RNA Sequence Rather than Protein Sequence Determines the Level of Self-synthesis and Antigen Presentation of the EBV-encoded Antigen,EBNA1

Unique purine-rich mRNA sequences embedded in the coding sequences of a distinct group of gammaherpesvirus maintenance proteins underlie the ability of the latently infected cell to minimize immune recognition. The Epstein-Barr virus nuclear antigen, EBNA1, a well characterized lymphocryptovirus maintenance protein has been shown to inhibit in cis antigen presentation, due in part to a large internal repeat domain encoding glycine and alanine residues (GAr) encoded by a purine-rich mRNA sequence. Recent studies have suggested that it is the purine-rich mRNA sequence of this repeat region rather than the encoded GAr polypeptide that directly inhibits EBNA1 self-synthesis and contributes to immune evasion. To test this hypothesis, we generated a series of EBNA1 internal repeat frameshift constructs and assessed their effects on cis-translation and endogenous antigen presentation. Diverse peptide sequences resulting from alternative repeat reading frames did not alleviate the translational inhibition characteristic of EBNA1 self-synthesis or the ensuing reduced surface presentation of EBNA1-specific peptide-MHC class I complexes. Human cells expressing the EBNA1 frameshift variants were also poorly recognized by antigen-specific T-cells. Furthermore, a comparative analysis of the mRNA sequences of the corresponding repeat regions of different viral maintenance homologues highlights the high degree of identity between the nucleotide sequences despite very little homology in the encoded amino acid sequences. Based on these combined observations, we propose that the cis-translational inhibitory effect of the EBNA1 internal repeat sequence operates mechanistically at the nucleotide level, potentially through RNA secondary structural elements, and is unlikely to be mediated through the GAr polypeptide. The demonstration that the EBNA1 repeat mRNA sequence and not the encoded protein sequence underlies immune evasion in this class of virus suggests a novel approach to therapeutic development through the use of anti-sense strategies or small molecules targeting EBNA1 mRNA structure.     

Sequence Diversity of Bacillus thuringiensis Flagellin (H Antigen) Protein at the Intra-H Serotype Level

In Bacillus thuringiensis, the hag gene encodes flagellin, the protein responsible for eliciting the immunological reaction in H serotyping. Specific flagellin amino acid sequences have been correlated to specific B. thuringiensis H serotypes, H1 to H67. Ten H serotypes, however, contain three or more antigenic subfactors, labeled a, b, c, d, or e, and have been subdivided into 23 serovars. In the present study, we set out to analyze the sequence diversity of flagellins among serovars from the same H serotypes. We studied the hag genes in 39 B. thuringiensis strains representing the 23 serovars from the 10 H serotypes mentioned above. A serovar and a biovar from an 11th H serotype were also included. The hag genes were amplified and cloned and their nucleotide sequences were determined and translated into amino acid sequences, or the sequences were retrieved directly from GenBank when available. Strains of the H3 serotype contained two or three copies of the fla gene, an ortholog of the hag gene. Strains of the H6 serotype contained three copies. Strains of all other H serotypes each contained a single copy of the hag gene. Alignments of amino acid sequences from all copies in all strains of the H3 serotype revealed short signature sequences, GGAG and SGG, GPDPDDAVKNLT, and DITTTK, that appeared to be specific to the H3c, H3d, and H3e antigenic subfactors, respectively. Similar short signature sequences, GDIT, AFIK, TSAGKA, and SAPSKG, were revealed for H8b, H8c, H20b, and H20c, respectively. Amino acid sequences in the flagellin central variable region were highly conserved among serovars of the H3, H5, H11, and H20 serotypes and much more divergent among serovars of the H4, H10, H18, H24, and H28 serotypes. Two bootstrapped neighbor-joining trees were respectively generated from the alignments of the amino acid sequences translated from all copies of the hag genes in the B. thuringiensis strains of the H3 and H6 serotypes. Sequence identities and relationships were revealed. A third bootstrapped neighbor-joining tree was generated, this one from the alignment of the flagellin amino acid sequences from all the B. thuringiensis strains in the study. Eight clusters, I to VIII, were revealed. Although most clusters contained strains and serovars from the same H serotype, clusters VII and VIII contained serovars from different H serotypes.     

A Cluster of Cuticle Protein Genes of Drosophila Melanogaster at 65a: Sequence, Structure and Evolution

A 36-kb genomic DNA segment of the Drosophila melanogaster genome containing 12 clustered cuticle genes has been mapped and partially sequenced. The cluster maps at 65A 5-6 on the left arm of the third chromosome, in agreement with the previously determined location of a putative cluster encompassing the genes for the third instar larval cuticle proteins LCP5, LCP6 and LCP8. This cluster is the largest cuticle gene cluster discovered to date and shows a number of surprising features that explain in part the genetic complexity of the LCP5, LCP6 and LCP8 loci. The genes encoding LCP5 and LCP8 are multiple copy genes and the presence of extensive similarity in their coding regions gives the first evidence for gene conversion in cuticle genes. In addition, five genes in the cluster are intronless. Four of these five have arisen by retroposition. The other genes in the cluster have a single intron located at an unusual location for insect cuticle genes.     

Laboratory-Directed Protein Evolution

Systematic approaches to directed evolution of proteins have been documented since the 1970s. The ability to recruit new protein functions arises from the considerable substrate ambiguity of many proteins. The substrate ambiguity of a protein can be interpreted as the evolutionary potential that allows a protein to acquire new specificities through mutation or to regain function via mutations that differ from the original protein sequence. All organisms have evolutionarily exploited this substrate ambiguity. When exploited in a laboratory under controlled mutagenesis and selection, it enables a protein to “evolve” in desired directions. One of the most effective strategies in directed protein evolution is to gradually accumulate mutations, either sequentially or by recombination, while applying selective pressure. This is typically achieved by the generation of libraries of mutants followed by efficient screening of these libraries for targeted functions and subsequent repetition of the process using improved mutants from the previous screening. Here we review some of the successful strategies in creating protein diversity and the more recent progress in directed protein evolution in a wide range of scientific disciplines and its impacts in chemical, pharmaceutical, and agricultural sciences.     

Evolution of Serpin Specificity: Cooperative Interactions in the Reactive-Site Loop Sequence of Antithrombin Specifically Restrict the Inhibition of Activated Protein C

Protease cascades and their inhibitors are a common feature of many biological regulatory systems, and the various components of such cascades have been subjected to a long and concerted evolution. We present here evidence that in the coagulation cascade, the sequence of the protease-binding reactive-site loop of antithrombin has evolved such that the majority of its residues has been acquired not for the efficient inhibition of its target proteases, thrombin and factor Xa, but to avoid the inhibition of activated protein C (APC). We substituted residues of the reactive-site loop of antithrombin into α₁-antitrypsin and tested the chimeras against thrombin, factor Xa, and APC. With respect to factor Xa and thrombin, the difference in association rate between the fastest and the slowest inhibitors was 5.5- and 88-fold, respectively. However, with respect to APC the difference was 12,500-fold. While most of the variation in the inhibition rates of thrombin could be accounted for by P2 Gly-to-Pro substitutions, for APC almost every residue had an effect on inhibition. In 22 of 25 direct comparisons of antitrypsin residues with antithrombin residues, either singly or in blocs, the antithrombin residues caused a decrease in the rate of inhibition of APC. The antithrombin residue Asn393, at position P′3, emerged as particularly important for avoiding the inhibition of APC, however, its 190-fold effect was seen only when in conjunction with antithrombin P7 to P′2 residues. Cooperative effects among residues of the reactive-site loop thus emerged as critical for restricting the activity of this sequence against APC. Received: 15 November 1999 / Accepted: 2 August 2000     

Episodic Evolution of Protein Hormones: Molecular Evolution of Pituitary Prolactin

Previous studies have shown that pituitary growth hormone displays an episodic pattern of evolution, with a slow underlying evolutionary rate and occasional sustained bursts of rapid change. The present study establishes that pituitary prolactin shows a similar pattern. During much of tetrapod evolution the sequence of prolactin has been strongly conserved, showing a slow basal rate of change (approx 0.27 × 10⁹ substitutions/amino acid site/year). This rate has increased substantially (∼12- to 38-fold) on at least four occasions during eutherian evolution, during the evolution of primates, artiodactyls, rodents, and elephants. That these increases are real and not a consequence of inadvertant comparison of paralogous genes is shown (for at least the first three groups) by the fact that they are confined to mature protein coding sequence and not apparent in sequences coding for signal peptides or when synonymous substitutions are examined. Sequences of teleost prolactins differ markedly from those of tetrapods and lungfish, but during the course of teleost evolution the rate of change of prolactin has been less variable than that of growth hormone. It is concluded that the evolutionary pattern seen for prolactin shows long periods of near-stasis interrupted by occasional bursts of rapid change, resembling the pattern seen for growth hormone in general but not in detail. The most likely basis for these bursts appears to be adaptive evolution though the biological changes involved are relatively small. Received: 31 August 1999 / Accepted: 9 February 2000     

Man and Sea in the Mesolithic: Coastal Settlement above and below Present Sea Level

Man and Sea in the Mesolithic: Coastal Settlement above and below Present Sea Level. Anders Fischer. ed. Oxbow Monograph, 53. Oakville, CT. David Brown Book Company, 1995. 440 pp.     

Molecular Phylogeny of Grey Mullets Based on Mitochondrial DNA Sequence Analysis: Evidence of a Differential Rate of Evolution at the Intrafamily Level

The phylogeny of the grey mullets is considered problematic both at the intra- and interfamily level. Such a difficulty arises from the highly homogeneous morphology displayed by this group of fish and, consequently, from the paucity of the key morphological characters suitable to address their phylogeny and evolution. In the present work, we have approached the phylogenetic and evolutionary relationships of seven species of Mugilidae, six of which from the Mediterranean Sea, on the basis of the DNA sequences of two mitochondrial genes (cytochromeband 12S rRNA). Despite the morphological homogeneity exhibited by the taxa considered, the two species of the genusMugil(M. cephalusandM. curema) showed a remarkable genetic divergence compared to all the other members of the family. The relative rate test revealed a significantly higher rate of evolution along theMugillineage.     

Statistical Alignment with a Sequence Evolution Model Allowing Rate Heterogeneity along the Sequence

We present a stochastic sequence evolution model to obtain alignments and estimate mutation rates between two homologous sequences. The model allows two possible evolutionary behaviors along a DNA sequence in order to determine conserved regions and take its heterogeneity into account. In our model, the sequence is divided into slow and fast evolution regions. The boundaries between these sections are not known. It is our aim to detect them. The evolution model is based on a fragment insertion and deletion process working on fast regions only and on a substitution process working on fast and slow regions with different rates. This model induces a pair hidden Markov structure at the level of alignments, thus making efficient statistical alignment algorithms possible. We propose two complementary estimation methods, namely, a Gibbs sampler for Bayesian estimation and a stochastic version of the EM algorithm for maximum likelihood estimation. Both algorithms involve the sampling of alignments. We propose a partial alignment sampler, which is computationally less expensive than the typical whole alignment sampler. We show the convergence of the two estimation algorithms when used with this partial sampler. Our algorithms provide consistent estimates for the mutation rates and plausible alignments and sequence segmentations on both simulated and real data.     

cDNA Sequence, Protein Structure, and Evolution of the Single Hemocyanin from Aplysia californica, an Opisthobranch Gastropod

By protein immunobiochemistry and cDNA sequencing, we have found only a single hemocyanin polypeptide in an opisthobranch gastropod, the sea hare Aplysia californica, which contrasts with previously studied prosobranch gastropods, which express two distinct isoforms of this extracellular respiratory protein. We have cloned and sequenced the cDNA encoding the complete polypeptide of Aplysia californica hemocyanin (AcH). The cDNA comprises 11,433 bp, encompassing a 5UTR of 77 bp, a 3UTR of 1057 bp, and an open reading frame for a signal peptide of 20 amino acids plus a polypeptide of 3412 amino acids (M_r ca. 387 kDa). This polypeptide is the subunit of the cylindrical native hemocyanin (M_r ca. 8 MDa). It comprises eight different functional units (FUs: a, b, c, d, e, f, g, h) that have been identified immunobiochemically after limited proteolysis of AcH purified from the hemolymph. Each FU shows a highly conserved copper-A and copper-B site for reversible oxygen binding. FU AcH-h carries a specific C-terminal extension of ca. 100 amino acids that include two cysteines that may be utilized for disulfide bridge formation. Potential N-glycosylation sites are present in six FUs but lacking in AcH-b and AcH-c. On the basis of multiple sequence alignments, phylogenetic trees and a statistically firm molecular clock were calculated. The latter suggests that the last common ancestor of Haliotis and Aplysia lived 373±47 million years ago, in convincing agreement with fossil records from the early Devonian. However, the gene duplication yielding the two distinct hemocyanin isoforms found today in Haliotis tuberculata occurred 343±43 million years ago.[Reviewing Editor: Dr. Axel Meyer]The sequence reported in this paper has been deposited in the GenBank database under accession number AJ556169.