Amino Acid Substitution Rates Research Articles

Dissection of Functional Residues in Receptor Activity-Modifying Proteins through Phylogenetic and Statistical Analyses

Type I and type-II functional divergences have been stated to highlight specific residues carrying out differential functions in evolutionary-divergent protein clusters from a single common ancestor. Briefly, type I analysis is based on residue constraints reflecting a gain of function just in one cluster of an entire family of proteins; while the type-II approach is based on residue constraints showing a different chemical nature in every cluster of a protein family. This last evidence is understood as differential functionality among clusters. The Receptor Activity-Modifying Proteins constitute a family characterized by its paralogous distribution in vertebrates. They are known as G-Protein Coupled Receptor modulators. Although several studies have determined their involvement in ligand binding, specificity, and enhancement of signal transduction, the responsible residues supporting those functions are unclear. Using different bioinformatic approaches, we predicted residues involved in different RAMP functional tasks. Many residues localized in an extracellular coil of RAMP proteins were predicted to be under functional divergence suggesting a gain of function in their respective proteins. Interestingly, the transmembrane region also showed important results for residues playing relevant roles where most of them showed a biased distribution on the structure. A relevant role was conferred by the enrichment of type-II residues observed in their sequences. We show a collection of residues explaining possible gain of function and differential functionality in RAMP proteins. These residues are still experimentally unexplored with regards to functionality. Finally, an evolutionary history could be discerned. Mainly, the RAMP2 cluster has evolved in a higher manner than other RAMP clusters. However, a deacceleration in the aminoacid substitution rate of RAMP2 was observed in mammals. Such effect could be caused by the co-evolution of ligands and receptors interacting with RAMP2 through evolution and/or the specialization of this cluster in GPCR modulation.

Evaluating Evolutionary Constraint on the Rapidly Evolving Gene matK Using Protein Composition

The rapidly evolving chloroplast matK gene has nucleotide and amino acid substitution rates suggestive of progression toward a pseudogene state. However, molecular evidence has demonstrated that matK is expressed and functional. We explore in this paper the underlying factors behind the mode and tempo of matK evolution that allow this protein coding gene to accommodate such elevated rates of substitution and yet maintain functionality. Conservative amino acid replacement may reconcile the fast evolutionary rate in matK with conservation in protein function. Based on this premise, we have examined putative amino acid sequences for MATK from across green plants to determine constraint on this protein as indicated by variation in composition of amino acid side chain category. Amino acids in the MATK ORF were divided into six categories based on chemical properties of their side chains: nonpolar, uncharged (pH 7), basic, acidic, aromatic, and "special" (amino acids that specifically affect protein structure, i.e., proline, glycine, and cysteine). The amount of standard deviation (SD) in side chain composition was used as a measure of variation and constraint, where a low SD implied high evolutionary constraint and a high SD implied low constraint. Further, we used secondary structure prediction to evaluate if conservation observed in side chain composition was reflected in stable predicted structure. The results of this study demonstrate evolutionary constraint on MATK, identify three regions of functional importance, show highly conserved secondary structure, and support the putative function of MATK as a group II intron maturase.

Parameters of proteome evolution from histograms of amino-acid sequence identities of paralogous proteins

BackgroundThe evolution of the full repertoire of proteins encoded in a given genome is mostly driven by gene duplications, deletions, and sequence modifications of existing proteins. Indirect information about relative rates and other intrinsic parameters of these three basic processes is contained in the proteome-wide distribution of sequence identities of pairs of paralogous proteins.ResultsWe introduce a simple mathematical framework based on a stochastic birth-and-death model that allows one to extract some of this information and apply it to the set of all pairs of paralogous proteins in H. pylori, E. coli, S. cerevisiae, C. elegans, D. melanogaster, and H. sapiens. It was found that the histogram of sequence identities p generated by an all-to-all alignment of all protein sequences encoded in a genome is well fitted with a power-law form ~ p-γ with the value of the exponent γ around 4 for the majority of organisms used in this study. This implies that the intra-protein variability of substitution rates is best described by the Gamma-distribution with the exponent α ≈ 0.33. Different features of the shape of such histograms allow us to quantify the ratio between the genome-wide average deletion/duplication rates and the amino-acid substitution rate.ConclusionWe separately measure the short-term ("raw") duplication and deletion rates , which include gene copies that will be removed soon after the duplication event and their dramatically reduced long-term counterparts rdup, rdel. High deletion rate among recently duplicated proteins is consistent with a scenario in which they didn't have enough time to significantly change their functional roles and thus are to a large degree disposable. Systematic trends of each of the four duplication/deletion rates with the total number of genes in the genome were analyzed. All but the deletion rate of recent duplicates were shown to systematically increase with Ngenes. Abnormally flat shapes of sequence identity histograms observed for yeast and human are consistent with lineages leading to these organisms undergoing one or more whole-genome duplications. This interpretation is corroborated by our analysis of the genome of Paramecium tetraurelia where the p-4 profile of the histogram is gradually restored by the successive removal of paralogs generated in its four known whole-genome duplication events.

Estimation of DNA Sequence Context-dependent Mutation Rates Using Primate Genomic Sequences

It is understood that DNA and amino acid substitution rates are highly sequence context-dependent, e.g., C --> T substitutions in vertebrates may occur much more frequently at CpG sites and that cysteine substitution rates may depend on support of the context for participation in a disulfide bond. Furthermore, many applications rely on quantitative models of nucleotide or amino acid substitution, including phylogenetic inference and identification of amino acid sequence positions involved in functional specificity. We describe quantification of the context dependence of nucleotide substitution rates using baboon, chimpanzee, and human genomic sequence data generated by the NISC Comparative Sequencing Program. Relative mutation rates are reported for the 96 classes of mutations of the form 5' alphabetagamma 3' --> 5' alphadeltagamma 3', where alpha, beta, gamma, and delta are nucleotides and beta not equal delta, based on maximum likelihood calculations. Our results confirm that C --> T substitutions are enhanced at CpG sites compared with other transitions, relatively independent of the identity of the preceding nucleotide. While, as expected, transitions generally occur more frequently than transversions, we find that the most frequent transversions involve the C at CpG sites (CpG transversions) and that their rate is comparable to the rate of transitions at non-CpG sites. A four-class model of the rates of context-dependent evolution of primate DNA sequences, CpG transitions > non-CpG transitions approximately CpG transversions > non-CpG transversions, captures qualitative features of the mutation spectrum. We find that despite qualitative similarity of mutation rates among different genomic regions, there are statistically significant differences.

Evidence for positive Darwinian selection on the hepcidin gene of Perciform and Pleuronectiform fishes

Hepcidin is a small cysteine-rich peptide that plays an important role in antimicrobial activity and in maintaining iron homeostasis in vertebrates. Here we report on the underlying mechanism that maintains high sequence diversities among the hepcidin-like variants of perciform and pleuronectiform fishes. In contrast to mammals, maximum likelihood-based codon substitution analyses revealed that positive Darwinian selection (nonsynonymous to synonymous substitution, omega > 1) is the likely cause of accelerated rate of amino acid substitutions in the hepcidin mature peptide region of these fishes. Comparison of models incorporating positive selection (omega > 1) at certain sites with models not incorporating positive selection (omega < 1) failed to reject (p = 0) the evidence of positive selection among the codon sites of percifom and pleuronectiform hepcidin. The adaptive evolution of this peptide in perciform and pleuronectiform fishes might be directed by pathogens when the host is exposed to new habitats/environments.

Gene Duplication and Adaptive Evolution of Digestive Proteases in Drosophila arizonae Female Reproductive Tracts

It frequently has been postulated that intersexual coevolution between the male ejaculate and the female reproductive tract is a driving force in the rapid evolution of reproductive proteins. The dearth of research on female tracts, however, presents a major obstacle to empirical tests of this hypothesis. Here, we employ a comparative EST approach to identify 241 candidate female reproductive proteins in Drosophila arizonae, a repleta group species in which physiological ejaculate–female coevolution has been documented. Thirty-one of these proteins exhibit elevated amino acid substitution rates, making them candidates for molecular coevolution with the male ejaculate. Strikingly, we also discovered 12 unique digestive proteases whose expression is specific to the D. arizonae lower female reproductive tract. These enzymes belong to classes most commonly found in the gastrointestinal tracts of a diverse array of organisms. We show that these proteases are associated with recent, lineage-specific gene duplications in the Drosophila repleta species group, and exhibit strong signatures of positive selection. Observation of adaptive evolution in several female reproductive tract proteins indicates they are active players in the evolution of reproductive tract interactions. Additionally, pervasive gene duplication, adaptive evolution, and rapid acquisition of a novel digestive function by the female reproductive tract points to a novel coevolutionary mechanism of ejaculate–female interaction.

Differential affinity of mammalian histone H1 somatic subtypes for DNA and chromatin

BackgroundHistone H1 is involved in the formation and maintenance of chromatin higher order structure. H1 has multiple isoforms; the subtypes differ in timing of expression, extent of phosphorylation and turnover rate. In vertebrates, the amino acid substitution rates differ among subtypes by almost one order of magnitude, suggesting that each subtype might have acquired a unique function. We have devised a competitive assay to estimate the relative binding affinities of histone H1 mammalian somatic subtypes H1a-e and H1° for long chromatin fragments (30–35 nucleosomes) in physiological salt (0.14 M NaCl) at constant stoichiometry.ResultsThe H1 complement of native chromatin was perturbed by adding an additional amount of one of the subtypes. A certain amount of SAR (scaffold-associated region) DNA was present in the mixture to avoid precipitation of chromatin by excess H1. SAR DNA also provided a set of reference relative affinities, which were needed to estimate the relative affinities of the subtypes for chromatin from the distribution of the subtypes between the SAR and the chromatin. The amounts of chromatin, SAR and additional H1 were adjusted so as to keep the stoichiometry of perturbed chromatin similar to that of native chromatin. H1 molecules freely exchanged between the chromatin and SAR binding sites. In conditions of free exchange, H1a was the subtype of lowest affinity, H1b and H1c had intermediate affinities and H1d, H1e and H1° the highest affinities. Subtype affinities for chromatin differed by up to 19-fold. The relative affinities of the subtypes for chromatin were equivalent to those estimated for a SAR DNA fragment and a pUC19 fragment of similar length. Avian H5 had an affinity ~12-fold higher than H1e for both DNA and chromatin.ConclusionH1 subtypes freely exchange in vitro between chromatin binding sites in physiological salt (0.14 M NaCl). The large differences in relative affinity of the H1 subtypes for chromatin suggest that differential affinity could be functionally relevant and thus contribute to the functional differentiation of the subtypes. The conservation of the relative affinities for SAR and non-SAR DNA, in spite of a strong preference for SAR sequences, indicates that differential affinity alone cannot be responsible for the heterogeneous distribution of some subtypes in cell nuclei.

Molecular clock: an anti-neo-Darwinian legacy.

S early as the time of the rediscovery of Mendel’slaws, George H. F. Nuttal measured the amount ofprecipitateofnormalserafromgreatapes,monkeys,andsome other mammals. His crude method, using rabbitantiserum directed against whole human serum, in-dicated that except for flying lemurs, the amount ofprecipitatedeclinedwiththeevolutionarydistancefromhumans (Nuttal 1904; Klein 1995). The immunolog-ical method was later refined and played an importantrole in reconstructing primate phylogenies (Goodman1962). In 1962, almost in parallel, Emile Zuckerkandlcollaborated with Linus Pauling at Caltech on hemoglo-bin evolution and expressed the idea of ‘‘molecular an-thropology’’asanewdiscipline(Z uckerkandl1963).Theidea was optimistic and ahead of the times, but MorrisGoodman shared it. In the same year, Zuckerkandland Pauling (1962) calibrated the amino acid substi-tution rate in mammalian hemoglobins and estimatedthe divergence times of orthologous and paralogoushemoglobins. Clearly, the immunological and proteinsequence data had already provided the germ for im-munological and protein clocks (see also Margoliash1963 for cytochrome c; Doolittle and Blomba¨ck1964for fibrinopeptides).The time was ripe for Zuckerkandl and Pauling(1965, p. 138) to advocate the concept of a molecularevolutionary clock: ‘‘Anyone who recognizes the valueof the immunological approach for estimating phyleticdistance with certain limits should find it impossible todeny that the comparison of amino acid sequences ispotentially an even better tool. It is only potentially lessequivocal, more accurate, suited for absolute instead ofonly relative evaluations, and able to extrapolate fromthe present to the past.’’ The stochastic nature of themolecular clock was well recognized and described by aPoisson process for the first time. It was also pointed outthat, for a molecular clock to exist, amino acid changesmust be limited almost exclusively to functionally nearlyneutral changes, although not only random genetic driftbut also Darwinian selection was invoked for fixation ofsuch changes. Thus, the discovery of a molecular clocksupportedtheconceptofnearneutrality at the molecularlevel. From a historical point of view, it is of interest to askwho was responsible for the monumental proposition of amolecular clock, Zuckerkandl or Pauling. Thirty yearslater, Pauling recalled, ‘‘I think it was my idea, but I amnot sure. We were just collaborating on these studies.Perhaps it was Emile’s idea’’ (M organ 1998, p. 166).This recollection might be a Freudian memory lapse (E.Zuckerkandl, personal communication). The idea mustbe Emile Zuckerkandl’s, since it is clear in Z uckerkandland Pauling (1962), an article almost entirely writtenby Zuckerkandl (see also Zuckerkandl 1987).The high evolutionary rate estimated from hemoglo-bin and other proteins was a key to development of theneutral theory by Motoo Kimura, rather than, as some-times asserted, the large extent of electrophoreticallyobserved polymorphism (Dietrich 1994; Sua´rez andBarahona 1996). Indeed, Kimura’s 1968 article beginswith a discussion of the amino acid substitution rate ob-tained from Zuckerkandl and Pauling (1965). Fromthis, Kimura estimated the nucleotide substitution rateon the basis of codon degeneracy and extrapolated thisratetotheentiregenomeonthebasisofthetotalnumberof base pairs estimated by Muller (1958). This rate wastoo high to be accounted for by natural selection, ac-cording to Haldane’s (1957) cost of natural selection.Kimura often said that Muller’s estimate (4 3 10

Analysis of the complete hepatitis B virus genome in patients with genotype C chronic hepatitis in relation to HBeAg and anti‐HBe

To investigate the relationship between viral factors and the development of chronic hepatitis B, the entire hepatitis B virus (HBV) genome of chronic carriers at different disease stages were analyzed. Eighty genotype C HBV carriers including 12 hepatitis B e antigen (HBeAg) positive asymptomatic carriers (Group A), 49 HBeAg positive patients with chronic liver diseases (Group B) and 19 anti-HBe positive patients with chronic liver diseases (Group C) were studied. HBV nucleic acid from serum samples was sequenced directly and compared with GenBank reference sequences HBV X01587 and M12906. On phylogenetic analysis, 76 cases were genotype C2. Of the 76 genotype C2 cases, the nucleotide and amino acid substitution rates in the precore/core region were significantly higher in Groups B and C than in Group A, also in Group C than in Group B. The nucleotide substitution rates in the full genome and the core promoter region were significantly higher in Group C than in Group A, also in group C than in Group B. The nucleotide and amino acid substitution rates in the X region were significantly higher in Group C than in Group A. The amino acid substitution rate in the pre-S2 region was significantly higher in Group C than in Group B. Deletion mutations were found mainly in Groups B and C. This whole genome analysis of HBV chronic carriers suggested that the nucleotide substitutions and deletions in HBV were closely associated with the pathogenesis of chronic HBV infection.

RATES OF MOLECULAR EVOLUTION IN BACTERIA ARE RELATIVELY CONSTANT DESPITE SPORE DORMANCY

Rates of molecular evolution are known to vary considerably among lineages, partially due to differences in life-history traits such as generation time. The generation-time effect has been well documented in some eukaryotes, but its prevalence in prokaryotes is unknown. "Because many species of Firmicute bacteria spend long periods of time as metabolically dormant spores, which could result in fewer DNA substitutions per unit time, they present an excellent system for testing predictions of the molecular clock hypothesis." To test whether spore-forming bacteria evolve more slowly than their non-spore-forming relatives, I used phylogenetic methods to determine if there were differences in rates of amino acid substitution between spore-forming and non-spore-forming lineages of Firmicute bacteria. Although rates of evolution do vary among lineages, I find no evidence for an effect of spore-formation on evolutionary rate and, furthermore, evolutionary rates are similar to those calculated for enteric bacteria. These results support the notion that variation in generation time does not affect evolutionary rates in bacterial lineages.

Adaptive evolution after duplication of penaeidin antimicrobial peptides

Penaeidin antimicrobial peptides in penaeid shrimps are an important component of their innate immune system that provides immunity against infection caused by several gram-positive bacteria and filamentous fungal species. Despite the knowledge on the identification and characterization of these peptides in penaeid shrimps, little is known about the evolutionary pattern of these peptides and the underlying genetic mechanisms that maintain high sequence diversities in the penaeidin gene family. Based on the phylogenetic analyses and maximum likelihood-based codon substitution analyses, here we present the convincing evidence that multiple copies of penaeidins have evolved by gene duplication, and positive Darwinian selection (adaptive evolution) is the likely cause of accelerated rate of amino acid substitutions among these duplicated genes. While the average ratio of non-synonymous to synonymous substitutions (ω) for the entire coding region of both active domains is 0.9805, few codon sites showed significantly higher ω (3.73). The likelihood ratio tests that compare models incorporating positive selection (ω > 1) at certain codon sites with models not incorporating positive selection (ω < 1), failed to reject ( p = 0) the evidence of positive Darwinian selection. The rapid adaptive evolution of this gene family might be directed by the pathogens and the faster rate of amino acid substitutions in the N-terminal proline-rich and C-terminal cysteine-rich domains could be due to their direct involvement in the protection against pathogens. When the host expose to different habitats/environment an accelerated rate of amino acid substitutions in both the active domains may also be expected.

Adaptive evolution of the symbiotic gene NORK is not correlated with shifts of rhizobial specificity in the genus Medicago

BackgroundThe NODULATION RECEPTOR KINASE (NORK) gene encodes a Leucine-Rich Repeat (LRR)-containing receptor-like protein and controls the infection by symbiotic rhizobia and endomycorrhizal fungi in Legumes. The occurrence of numerous amino acid changes driven by directional selection has been reported in this gene, using a limited number of messenger RNA sequences, but the functional reason of these changes remains obscure. The Medicago genus, where changes in rhizobial associations have been previously examined, is a good model to test whether the evolution of NORK is influenced by rhizobial interactions.ResultsWe sequenced a region of 3610 nucleotides (encoding a 392 amino acid-long region of the NORK protein) in 32 Medicago species. We confirm that positive selection in NORK has occurred within the Medicago genus and find that the amino acid positions targeted by selection occur in sites outside of solvent-exposed regions in LRRs, and other sites in the N-terminal region of the protein. We tested if branches of the Medicago phylogeny where changes of rhizobial symbionts occurred displayed accelerated rates of amino acid substitutions. Only one branch out of five tested, leading to M. noeana, displays such a pattern. Among other branches, the most likely for having undergone positive selection is not associated with documented shift of rhizobial specificity.ConclusionAdaptive changes in the sequence of the NORK receptor have involved the LRRs, but targeted different sites than in most previous studies of LRR proteins evolution. The fact that positive selection in NORK tends not to be associated to changes in rhizobial specificity indicates that this gene was probably not involved in evolving rhizobial preferences. Other explanations (e.g. coevolutionary arms race) must be tested to explain the adaptive evolution of NORK.

Phylogenetic Reconstruction of Orthology, Paralogy, and Conserved Synteny for Dog and Human

Accurate predictions of orthology and paralogy relationships are necessary to infer human molecular function from experiments in model organisms. Previous genome-scale approaches to predicting these relationships have been limited by their use of protein similarity and their failure to take into account multiple splicing events and gene prediction errors. We have developed PhyOP, a new phylogenetic orthology prediction pipeline based on synonymous rate estimates, which accurately predicts orthology and paralogy relationships for transcripts, genes, exons, or genomic segments between closely related genomes. We were able to identify orthologue relationships to human genes for 93% of all dog genes from Ensembl. Among 1:1 orthologues, the alignments covered a median of 97.4% of protein sequences, and 92% of orthologues shared essentially identical gene structures. PhyOP accurately recapitulated genomic maps of conserved synteny. Benchmarking against predictions from Ensembl and Inparanoid showed that PhyOP is more accurate, especially in its predictions of paralogy. Nearly half (46%) of PhyOP paralogy predictions are unique. Using PhyOP to investigate orthologues and paralogues in the human and dog genomes, we found that the human assembly contains 3-fold more gene duplications than the dog. Species-specific duplicate genes, or “in-paralogues,” are generally shorter and have fewer exons than 1:1 orthologues, which is consistent with selective constraints and mutation biases based on the sizes of duplicated genes. In-paralogues have experienced elevated amino acid and synonymous nucleotide substitution rates. Duplicates possess similar biological functions for either the dog or human lineages. Having accounted for 2,954 likely pseudogenes and gene fragments, and after separating 346 erroneously merged genes, we estimated that the human genome encodes a minimum of 19,700 protein-coding genes, similar to the gene count of nematode worms. PhyOP is a fast and robust approach to orthology prediction that will be applicable to whole genomes from multiple closely related species. PhyOP will be particularly useful in predicting orthology for mammalian genomes that have been incompletely sequenced, and for large families of rapidly duplicating genes.

Longevity and the evolution of the mitochondrial DNA-coded proteins in mammals

The amino acids sequences of the mitochondrial DNA-coded peptides of placental mammals evolved at different rates in different branches of the mammalian phylogenetic tree. Adaptive selection was suggested to account for the faster evolution of some mitochondrial DNA-coded proteins in several branches of the mammalian tree, but the driving force(s) for the accelerated evolution has not been elucidated. Mitochondria generate reactive oxygen species (ROS) that appear to constrain the life span of many species. Therefore, I tested the hypothesis that the evolution of mammalian longevity drives the accelerated evolution of mitochondrial DNA-coded peptides. Using rodents as an outgroup for a clad that included most placental mammals (excluding rodents and hedgehogs) the computed rates of amino acid substitution per site were positively correlated with genus longevity (maximal observed averaged life span) for most of the mitochondrial DNA-coded peptides. The substitution per site of ATP6, the proton conducting subunit of ATPsynthase, CYTB, the core subunit of ubiquinone oxidoreductase that participate in both electron and proton transport, and ND3, a subunit of NADH dehydrogenase, showed the strongest correlations with longevity. Additional confirmation for the hypothesis was obtained by the observation that the genetic distances between placental mammals species that belong to different orders are positively correlated with the sum of longevities of the species pairs. The substitutions per site for the entire amino acid sequence coded by the heavy strand mtDNA were also positively correlated with the average longevities of the placental mammals orders. These results support the hypothesis that the evolution of longevity in mammals drove the accelerated evolution of mtDNA-coded peptide. It is suggested that, in mammals, adaptive selection of mutations that decrease the rate of production of reactive oxygen species, directly or indirectly (e.g. by increasing proton leak), increases longevity.

Clinical and molecular characteristics of patients with non-amyloid light chain deposition disorders, and outcome following treatment with high-dose melphalan and autologous stem cell transplantation

Light chain deposition disease (LCDD) is caused by a clonal plasma cell disorder in which fragments of monoclonal immunoglobulin light chains form non-fibrillary deposits in various tissues resulting in organ dysfunction. Crystal storing histiocytosis (CSH) is another light chain deposition disorder in which monoclonal light chains form intracytoplasmic crystals. Both are uncommon diseases for which there is limited treatment experience. Between 2003 and 2005, five patients with LCDD and one with CSH were treated at Boston University Medical Center with high-dose melphalan and autologous peripheral blood stem cell transplantation (HDM/SCT). Five of the six patients had predominantly renal involvement, and one patient with LCDD had biopsy-proven deposits in the myocardium. Molecular characterization revealed that the pathologic light chains were kappa in four of the six patients, and sequence analysis revealed unusual germline donor genes and high rates of amino-acid substitutions. One light chain sequence encoded a new potential N-linked glycosylation site, and another showed evidence of antigen selection. All patients are alive and five of the six patients are in complete hematologic remission at a median follow-up of 12 months (range 4-29 months) after HDM/SCT. In our experience, HDM/SCT is a feasible and effective treatment approach for these disorders.

The Relationship Between the Rate of Molecular Evolution and the Rate of Genome Rearrangement in Animal Mitochondrial Genomes

Evolution of mitochondrial genes is far from clock-like. The substitution rate varies considerably between species, and there are many species that have a significantly increased rate with respect to their close relatives. There is also considerable variation among species in the rate of gene order rearrangement. Using a set of 55 complete arthropod mitochondrial genomes, we estimate the evolutionary distance from the common ancestor to each species using protein sequences, tRNA sequences, and breakpoint distances (a measure of the degree of genome rearrangement). All these distance measures are correlated. We use relative rate tests to compare pairs of related species in several animal phyla. In the majority of cases, the species with the more highly rearranged genome also has a significantly higher rate of sequence evolution. Species with higher amino acid substitution rates in mitochondria also have more variable amino acid composition in response to mutation pressure. We discuss the possible causes of variation in rates of sequence evolution and gene rearrangement among species and the possible reasons for the observed correlation between the two rates.

Adaptive selection of mitochondrial complex I subunits during primate radiation

Mammalian oxidative phosphorylation (OXPHOS) complexes I, III, IV and V are assembled from both mitochondrial DNA (mtDNA) and nuclear DNA (nDNA) encoded subunits, with complex I encompassing 39 nDNA and seven mtDNA subunits. Yet the sequence variation of the mtDNA genes is more than ten fold greater than that of the nDNA encoded genes of the OXPHOS complexes and the mtDNA proteins have been found to be influenced by positive (adaptive) selection. To maintain a functional complex I, nDNA and mtDNA subunits must interact, implying that certain nDNA complex I genes may also have been influenced by positive selection. To determine if positive selection has influenced nDNA complex I genes, we analyzed the DNA sequences of all of the nDNA and mtDNA encoded complex I subunits from orangutan, gorilla, chimpanzee, human and all available vertebrate sequences. This revealed that three nDNA complex I genes (NDUFC2, NDUFA1, and NDUFA4) had significantly increased amino acid substitution rates by both PAML and Z-test, suggesting that they have been subjected to adaptive selection during primate radiation. Since all three of these subunits reside in the membrane domain of complex I along with the mtDNA subunits, we compared amino acid changes in these three nDNA genes with those of the mtDNA genes across species. Changes in the nDNA NDUFC2 cysteine 39 were found to correlate with those in the mtDNA ND5 cysteine 330. Therefore, adaptive selection has influenced some nDNA complex I genes and nDNA and mtDNA complex I genes may have co-evolved.

Adaptive Evolution of Fertilization Proteins within a Genus: Variation in ZP2 and ZP3 in Deer Mice (Peromyscus)

Rapid evolution of reproductive proteins has been documented in a wide variety of taxa. In internally fertilized species, knowledge about the evolutionary dynamics of these proteins between closely related taxa is primarily limited to accessory gland proteins in the semen of Drosophila. Investigation of additional taxa and functional classes of proteins is necessary in order to determine if there is a general pattern of adaptive evolution of reproductive proteins between recently diverged species. We performed an evolutionary analysis of 2 egg coat proteins, ZP2 and ZP3, in 15 species of deer mice (genus Peromyscus). Both of these proteins are involved in egg-sperm binding, a critical step in maintaining species-specific fertilization. Here, we show that Zp2 and Zp3 gene trees are not consistent with trees based on nonreproductive genes, Mc1r and Lcat, where species formed monophyletic clades. In fact, for both of the reproductive genes, intraspecific amino acid variation was extensive and alleles were sometimes shared across species. We document positive selection acting on ZP2 and ZP3 and identify specific amino acid sites that are likely targets of selection using both maximum likelihood approaches and patterns of parallel amino acid change. In ZP3, positively selected sites are clustered in and around the region implicated in sperm binding in Mus, suggesting changes may impact egg-sperm binding and fertilization potential. Finally, we identify lineages with significantly elevated rates of amino acid substitution using a Bayesian mapping approach. These findings demonstrate that the pattern of adaptive reproductive protein evolution found at higher taxonomic levels can be documented between closely related mammalian species, where reproductive isolation has evolved recently.

The Rate of Adaptive Evolution in Enteric Bacteria

Here we estimate the rate of adaptive substitution in a set of 410 genes that are present in 6 Escherichia coli and 6 Salmonella enterica genomes. We estimate that more than 50% of amino acid substitutions in this set of genes have been fixed by positive selection between the E. coli and S. enterica lineages. We also show that the proportion of adaptive substitutions is uncorrelated with the rate of amino acid substitution or gene function but that it may be correlated with levels of synonymous codon usage bias.

Protein Evolutionary Rates Correlate with Expression Independently of Synonymous Substitutions in Helicobacter pylori

In free-living microorganisms, such as Escherichia coli and Saccharomyces cerevisiae, both synonymous and nonsynonymous substitution frequencies correlate with expression levels. Here, we have tested the hypothesis that the correlation between amino acid substitution rates and expression is a by-product of selection for codon bias and translational efficiency in highly expressed genes. To this end, we have examined the correlation between protein evolutionary rates and expression in the human gastric pathogen Helicobacter pylori, where the absence of selection on synonymous sites enables the two types of substitutions to be uncoupled. The results revealed a statistically significant negative correlation between expression levels and nonsynonymous substitutions in both H. pylori and E. coli. We also found that neighboring genes located on the same, but not on opposite strands, evolve at significantly more similar rates than random gene pairs, as expected by co-expression of genes located in the same operon. However, the two species differ in that synonymous substitutions show a strand-specific pattern in E. coli, whereas the weak similarity in synonymous substitutions for neighbors in H. pylori is independent of gene orientation. These results suggest a direct influence of expression levels on nonsynonymous substitution frequencies independent of codon bias and selective constraints on synonymous sites.