SEARCH · Search PubMed

Results for “synonymous codon”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 847 records · Page 47Linked to original sources

PubMed2004

Antigenic and/or phase variation of PorA protein in non-subtypable Neisseria meningitidis strains isolated in Spain.

The PorA protein is a potential candidate as a vaccine component against meningococcal disease. However, this protein experiences antigenic variation and is subject to phase variations to evade immune selective pressure. In this study, the mechanisms responsible for altered expression of the PorA protein were analysed in 50 non-subtypable strains isolated from patients with meningococcal disease in Spain. The porA gene was amplified from 47 of the 50 strains. The majority of isolates were not recognized by the subtyping panel, as a result of non-synonymous base changes in the variable regions of the porA gene. Two of these strains revealed a premature stop codon before the variable region VR1 of PorA due to a single base-pair substitution at position 109 of the porA coding region. Another two presented a homopolymeric tract of eight adenine residues in the coding region, producing a DNA strand-slippage mechanism and PorA phase variation.

Antigenic Variation↗

PubMed2005

Evolutionary selection across the nuclear hormone receptor superfamily with a focus on the NR1I subfamily (vitamin D, pregnane X, and constitutive androstane receptors).

BACKGROUND: The nuclear hormone receptor (NR) superfamily complement in humans is composed of 48 genes with diverse roles in metabolic homeostasis, development, and detoxification. In general, NRs are strongly conserved between vertebrate species, and few examples of molecular adaptation (positive selection) within this superfamily have been demonstrated. Previous studies utilizing two-species comparisons reveal strong purifying (negative) selection of most NR genes, with two possible exceptions being the ligand-binding domains (LBDs) of the pregnane X receptor (PXR, NR1I2) and the constitutive androstane receptor (CAR, NR1I3), two proteins involved in the regulation of toxic compound metabolism and elimination. The aim of this study was to apply detailed phylogenetic analysis using maximum likelihood methods to the entire complement of genes in the vertebrate NR superfamily. Analyses were carried out both across all vertebrates and limited to mammals and also separately for the two major domains of NRs, the DNA-binding domain (DBD) and LBD, in addition to the full-length sequences. Additional functional data is also reported for activation of PXR and the vitamin D receptor (VDR; NR1I1) to gain further insight into the evolution of the NR1I subfamily. RESULTS: The NR genes appear to be subject to strong purifying selection, particularly in the DBDs. Estimates of the ratio of the non-synonymous to synonymous nucleotide substitution rates (the omega ratio) revealed that only the PXR LBD had a sub-population of codons with an estimated omega ratio greater than 1. CAR was also unusual in showing high relative omega ratios in both the DBD and LBD, a finding that may relate to the recent appearance of the CAR gene (presumably by duplication of a pre-mammalian PXR gene) just prior to the evolution of mammals. Functional analyses of the NR1I subfamily show that human and zebrafish PXRs show similar activation by steroid hormones and early bile salts, properties not shared by sea lamprey, mouse, or human VDRs, or by Xenopus laevis PXRs. CONCLUSION: NR genes generally show strong sequence conservation and little evidence for positive selection. The main exceptions are PXR and CAR, genes that may have adapted to cross-species differences in toxic compound exposure.

Journal Article↗

PubMed1997

Messenger RNA release from ribosomes during 5'-translational blockage by consecutive low-usage arginine but not leucine codons in Escherichia coli.

In '5'-translational blockage', significantly reduced yields of proteins are synthesized in Escherichia coli when consecutive low-usage codons are inserted near translation starts of messages (with reduced or no effect when these same codons are inserted downstream). We tested the hypothesis that ribosomes encountering these low-usage codons near the translation start prematurely release the mRNA. RNA from polysome gradients was fractionated into pools of polysomes and monosomes and a ribosome-free pool. New hybridization probes, called 'molecular beacons', and standard slot blots were used to detect test messages containing either consecutive low-usage AGG (arginine) or synonymous high-usage CGU insertions near the 5' end. The results show an approximately twofold increase in the ratio of free to bound mRNA when the low-usage codons were present in the message compared with when high-usage codons were present. In contrast, there was no difference in the ratio of free to bound mRNA when consecutive low-usage CUA or high-usage CUG (leucine) codons were inserted or when the arginine codons were inserted near the 3' end. These data indicate that at least some mRNA is released from ribosomes during 5'-translational blockage by arginine but not leucine codons, and they support proposals that premature termination of translation can occur in some conditions in vivo in the absence of a stop codon.

Arginine↗

PubMed2000

Gene expression, amino acid conservation, and hydrophobicity are the main factors shaping codon preferences in Mycobacterium tuberculosis and Mycobacterium leprae.

Mycobacterium tuberculosis and Mycobacterium leprae are the ethiological agents of tuberculosis and leprosy, respectively. After performing extensive comparisons between genes from these two GC-rich bacterial species, we were able to construct a set of 275 homologous genes. Since these two bacterial species also have a very low growth rate, translational selection could not be so determinant in their codon preferences as it is in other fast-growing bacteria. Indeed, principal-components analysis of codon usage from this set of homologous genes revealed that the codon choices in M. tuberculosis and M. leprae are correlated not only with compositional constraints and translational selection, but also with the degree of amino acid conservation and the hydrophobicity of the encoded proteins. Finally, significant correlations were found between GC3 and synonymous distances as well as between synonymous and nonsynonymous distances.

Amino Acid Sequence↗

PubMed2005

Comparative analysis of protein coding sequences from human, mouse and the domesticated pig.

BACKGROUND: The availability of abundant sequence data from key model organisms has made large scale studies of molecular evolution an exciting possibility. Here we use full length cDNA alignments comprising more than 700,000 nucleotides from human, mouse, pig and the Japanese pufferfish Fugu rubrices in order to investigate 1) the relationships between three major lineages of mammals: rodents, artiodactyls and primates, and 2) the rate of evolution and the occurrence of positive Darwinian selection using codon based models of sequence evolution. RESULTS: We provide evidence that the evolutionary splits among primates, rodents and artiodactyls happened shortly after each other, with most gene trees favouring a topology with rodents as outgroup to primates and artiodactyls. Using an unrooted topology of the three mammalian species we show that since their diversification, the pig and mouse lineages have on average experienced 1.44 and 2.86 times as many synonymous substitutions as humans, respectively, whereas the rates of non-synonymous substitutions are more similar. The analysis shows the highest average dN/dS ratio in the human lineage, followed by the pig and then the mouse lineages. Using codon based models we detect signals of positive Darwinian selection in approximately 5.3%, 4.9% and 6.0% of the genes on the human, pig and mouse lineages respectively. Approximately 16.8% of all the genes studied here are not currently annotated as functional genes in humans. Our analyses indicate that a large fraction of these genes may have lost their function quite recently or may still be functional genes in some or all of the three mammalian species. CONCLUSIONS: We present a comparative analysis of protein coding genes from three major mammalian lineages. Our study demonstrates the usefulness of codon-based likelihood models in detecting selection and it illustrates the value of sequencing organisms at different phylogenetic distances for comparative studies.

Animals↗

PubMed1993

Population genetics and phylogenetics of DNA sequence variation at multiple loci within the Drosophila melanogaster species complex.

Two regions of the genome, a 1-kbp portion of the zeste locus and a 1.1-kbp portion of the yolk protein 2 locus, were sequenced in six individuals from each of four species: Drosophila melanogaster, D. simulans, D. mauritiana, and D. sechellia. The species and strains were the same as those of a previous study of a 1.9-kbp region of the period locus. No evidence was found for recent balancing or directional selection or for the accumulation of selected differences between species. Yolk protein 2 has a high level of amino acid replacement variation and a low level of synonymous variation, while zeste has the opposite pattern. This contrast is consistent with information on gene function and patterns of codon bias. Polymorphism levels are consistent with a ranking of effective population sizes, from low to high, in the following order: D. sechellia, D. melanogaster, D.mauritiana, and D. simulans. The apparent species relationships are very similar to those suggested by the period locus study. In particular, D. simulans appears to be a large population that is still segregating variation that arose before the separation of D. mauritiana and D. sechellia. It is estimated that the separation of ancestral D. melanogaster from the other species occurred 2.5-3.4 Mya. The separations of D. sechellia and D. mauritiana from ancestral D. simulans appear to have occurred 0.58-0.86 Mya, with D. mauritiana having diverged from ancestral D. simulans 0.1 Myr more recently than D. sechellia.

Animals↗

PubMed2005

Hepatitis C virus F protein sequence reveals a lack of functional constraints and a variable pattern of amino acid substitution.

Hepatitis C virus (HCV) is an important human pathogen that affects 170 million people worldwide. The HCV genome is an RNA molecule that is approximately 9.6 kb in length and encodes a polyprotein that is cleaved proteolytically to generate at least 10 mature viral proteins. Recently, a new HCV protein named F has been described, which is synthesized as a result of a ribosomal frameshift. Little is known about the biological properties of this protein, but the possibility that the F protein may participate in HCV morphology or replication has been raised. In this work, the presence of functional constraints in the F protein was investigated. It was found that the rate of amino acid substitutions along the F protein was significantly higher than the rate of synonymous substitutions, and comparisons involving genes that represented independent phylogenetic lineages yielded very different divergence/conservation patterns. The distribution of stop codons in the F protein across all HCV genotypes was also investigated; genotypes 2 and 3 were found to have more stop codons than genotype 1. The results of this work suggest strongly that the pattern of divergence in the F protein is not affected by functional constraints.

Amino Acid Substitution↗

PubMed2006

Identifying sites under positive selection with uncertain parameter estimates.

Codon-based substitution models are routinely used to measure selective pressures acting on protein-coding genes. To this effect, the nonsynonymous to synonymous rate ratio (dN/dS = omega) is estimated. The proportion of amino-acid sites potentially under positive selection, as indicated by omega > 1, is inferred by fitting a probability distribution where some sites are permitted to have omega > 1. These sites are then inferred by means of an empirical Bayes or by a Bayes empirical Bayes approach that, respectively, ignores or accounts for sampling errors in maximum-likelihood estimates of the distribution used to infer the proportion of sites with omega > 1. Here, we extend a previous full-Bayes approach to include models with high power and low false-positive rates when inferring sites under positive selection. We propose some heuristics to alleviate the computational burden, and show that (i) full Bayes can be superior to empirical Bayes when analyzing a small data set or small simulated data, (ii) full Bayes has only a small advantage over Bayes empirical Bayes with our small test data, and (iii) Bayesian methods appear relatively insensitive to mild misspecifications of the random process generating adaptive evolution in our simulations, but in practice can prove extremely sensitive to model specification. We suggest that the codon model used to detect amino acids under selection should be carefully selected, for instance using Akaike information criterion (AIC).

Amino Acids↗

PubMed2006

Analysis of transitions at two-fold redundant sites in mammalian genomes. Transition redundant approach-to-equilibrium (TREx) distance metrics.

BACKGROUND: The exchange of nucleotides at synonymous sites in a gene encoding a protein is believed to have little impact on the fitness of a host organism. This should be especially true for synonymous transitions, where a pyrimidine nucleotide is replaced by another pyrimidine, or a purine is replaced by another purine. This suggests that transition redundant exchange (TREx) processes at the third position of conserved two-fold codon systems might offer the best approximation for a neutral molecular clock, serving to examine, within coding regions, theories that require neutrality, determine whether transition rate constants differ within genes in a single lineage, and correlate dates of events recorded in genomes with dates in the geological and paleontological records. To date, TREx analysis of the yeast genome has recognized correlated duplications that established a new metabolic strategies in fungi, and supported analyses of functional change in aromatases in pigs. TREx dating has limitations, however. Multiple transitions at synonymous sites may cause equilibration and loss of information. Further, to be useful to correlate events in the genomic record, different genes within a genome must suffer transitions at similar rates. RESULTS: A formalism to analyze divergence at two fold redundant codon systems is presented. This formalism exploits two-state approach-to-equilibrium kinetics from chemistry. This formalism captures, in a single equation, the possibility of multiple substitutions at individual sites, avoiding any need to "correct" for these. The formalism also connects specific rate constants for transitions to specific approximations in an underlying evolutionary model, including assumptions that transition rate constants are invariant at different sites, in different genes, in different lineages, and at different times. Therefore, the formalism supports analyses that evaluate these approximations. Transitions at synonymous sites within two-fold redundant coding systems were examined in the mouse, rat, and human genomes. The key metric (f2), the fraction of those sites that holds the same nucleotide, was measured for putative ortholog pairs. A transition redundant exchange (TREx) distance was calculated from f2 for these pairs. Pyrimidine-pyrimidine transitions at these sites occur approximately 14% faster than purine-purine transitions in various lineages. Transition rate constants were similar in different genes within the same lineages; within a set of orthologs, the f2 distribution is only modest overdispersed. No correlation between disparity and overdispersion is observed. In rodents, evidence was found for greater conservation of TREx sites in genes on the X chromosome, accounting for a small part of the overdispersion, however. CONCLUSION: The TREx metric is useful to analyze the history of transition rate constants within these mammals over the past 100 million years. The TREx metric estimates the extent to which silent nucleotide substitutions accumulate in different genes, on different chromosomes, with different compositions, in different lineages, and at different times.

Animals↗

PubMed1992

High expression vectors for the production of recombinant single-chain urinary plasminogen activator from Escherichia coli.

An expression cassette containing a synonymous gene for human single-chain urokinase-type plasminogen activator (Rscu-PA) 5'-flanked by a trp promoter and the Shine-Dalgarno sequence of the xyl A operon of Bacillus subtilis and terminated by the terminators trp A and Tn10 was constructed and inserted into a pBR322 derivative to yield pBF160. When compared to pUK54 trp 207-1 containing the natural scu-PA gene without the Shine-Dalgarno sequence and terminator, the expression efficiency of pBF160 in Escherichia coli strains was improved by one order of magnitude. Replacement of the trp by the tac promoter (pBF171) did not affect expression. Inserting the Shine-Dalgarno sequence and Tn10 terminator into pUK54 trp 207-1 (pWH1320) slightly increased the expression level, whereas elimination of the Shine-Dalgarno sequence and the terminators from pBF160 with almost complete conservation of the synonymous structural gene (pBF191) significantly reduced the expression. Variation of the distance between the Shine-Dalgarno sequence and the start codon between 8 and 10 bp (pBF163) proved irrelevant. In conclusion, poor expression of mammalian genes in E. coli may result from both improperly designed regulatory elements and structural features of the coding region and therefore de-novo synthesis of the gene may be required to obtain satisfactory expression.

Amino Acid Sequence↗

PubMed1996

Sequence analysis of a polymorphic Mhc class II gene in Pacific salmon.

Polymorphism of the nucleotide sequences encoding 149 amino acids of linked major histocompatibility complex (Mhc) class II B1 and B2 peptides, and of the intervening intron (548-773 base pairs), was examined within and among seven Pacific salmon (Oncorhynchus) species. Levels of nucleotide diversity were higher for the B1 sequence than for B2 or the intron in comparisons both within and between species. For the codons of the peptide binding region of the B1 sequence, the level of nonsynonymous nucleotide substitution (dN) exceeded the level of synonymous substitution (dS) by a factor of ten for within-species comparisons, and by a factor of four for between-species comparisons. The excess of dN indicates that balancing selection maintains diversity at this salmonid Mhc class II locus, as is common for Mhc loci in other vertebrates. Levels of nucleotide diversity for both the exon and intron sequences were greater among than within species, and there were numerous species-specific nucleotides present in both the coding and noncoding regions. Thus, neighbor-joining analysis of both the intron and exon regions provided phylogenies in which the sequences clustered strongly by species. There was little evidence of shared ancestral (trans-species) polymorphism in the exon phylogeny, and the intron phylogeny depicted standard relationships among the Pacific salmon species. The lack of shared allelic B1 lineages in these closely related species may result from severe bottlenecks that occurred during speciation or during the ice ages that glaciated the rim of the north Pacific Ocean approximately every 100000 years in the Pleistocene.

Amino Acid Sequence↗

PubMed2005

A specific haplotype in the 3' end of estrogen-receptor alpha gene is associated with low bone mineral density in premenopausal women and increased risk of postmenopausal osteoporosis.

It is well established that the development of postmenopausal osteoporosis is under genetic influence. We have recently identified a synonymous single nucleotide polymorphism (SNP) in exon 8 of estrogen receptor-alpha (ERalpha) gene in the vicinity of the stop codon (G2014A) that is associated with an increased risk of postmenopausal osteoporosis. In the present study, we attempted to locate SNPs in the 3'-unstranslated region (3'UTR) of the ERalpha gene that are in linkage disequilibrium with the exon 8 SNP and assessed their utilization in the risk assessment of postmenopausal osteoporosis in 352 Thai postmenopausal women. The association with bone mineral density (BMD) in premenopausal women was also investigated in 202 premenopausal women. A C to A SNP 1,748 nucleotides distal to the end of the stop codon (C3768A) was identified. The C3768A SNP was not overrepresented in subjects with osteoporosis. However, the presence of the A-C haplotype allele based on the A2014 and C3768 alleles was significantly related to the risk of osteoporosis independently of age, body weight, the G2014A and C3768A SNPs (odds ratio 2.36, 95% CI 1.42-3.91). Moreover, the presence of the A-C haplotype allele was associated with lower femoral neck BMD in premenopausal women ( P =0.05). We concluded that a specific haplotype in the 3' end of the ERalpha gene is associated with lower BMD in premenopausal women and is associated with a higher risk of osteoporosis in postmenopausal women. It is likely that the haplotype allele exerts its influence on bone as early as during young adulthood to increase the risk of osteoporosis later in life.

3' Untranslated Regions↗

PubMed1998

Extensive polymorphism of the FUT2 gene in an African (Xhosa) population of South Africa.

The human secretor type alpha(1,2)fucosyltransferase gene (FUT2) polymorphism was investigated in Xhosa and Caucasian populations of South Africa by polymerase chain reaction-restriction fragment length polymorphism and DNA sequencing. Six new base substitutions were found in the coding region of FUT2. A single base (C) deletion at nucleotide 778, which led to a frame shift and produced a stop codon at codon 275, was responsible for the enzyme inactivation. Three nonsynonymous base substitutions, A40G (lle14Val), C379T (Arg127Cys), and G481A (Asp161Asn), and two synonymous base substitutions, A375G (Glu125) and C480T (His160), were also identified in functional alleles. As a result, seven new alleles, Se40, Se481, Se40,481, Se357,480, Se357,379,480, Se375, and se357,480,778 were identified. Population studies revealed that an allele containing a nonsense mutation G428A (Trp143stop) (se428) was the common null allele in both Xhosa and Caucasian populations, whereas an allele containing a missense A385T (Ile129Phe) mutation (se357,385), which is the common null allele in Orientals, was found to be absent from both populations. The heterozygosity rates of FUT2 genotypes were as high as 0.75 in the Xhosa population and 0.65 in the Caucasian population. Therefore, the extensive polymorphism and race specificity of the FUT2 gene make it suitable for application as a new tool in genetic studies of modern human evolutionary history.

Alleles↗

PubMed1991

Nucleotide sequence of mouse Tcp-1a cDNA.

We have isolated complete cDNA clones encoding the mouse t-complex polypeptides 1A and 1B (TCP-1A and TCP-1B) from t-haplotype and wild-type (wt) mice, respectively. The complete nucleotide (nt) sequence of the Tcp-1a cDNA was determined. The Tcp-1a cDNA has an open reading frame (ORF) encoding a 60-kDa protein of 556 amino acids (aa). A comparison of nt sequences between the Tcp-1a and Tcp-1b cDNAs revealed that the 1786-bp regions upstream from their polyadenylation signals differed by 17 substitutions and that Tcp-1a had different polyadenylation sites from Tcp-1b. In these ORFs, 15 bp were substituted between the two alleles, occurring in 14 codons and resulting in eleven single-aa substitutions. Among these 15 substitutions, twelve were nonsynonymous (aa change) and three were synonymous (no aa change). The aa substitution in TCP-1 has occurred at least 20 times faster between t-haplotype and wt than between mouse and human or mouse and Drosophila.

Amino Acid Sequence↗

PubMed2006

Influenza A H5N1 hemagglutinin cleavable signal sequence substitutions.

Eleven influenza A H5N1 hemagglutinin N-terminal cleavable signal sequences, coded by single nucleotide substitutions relative to reference A/Viet Nam/1203/2004, were identified by BLASTN search of GenBank and were characterized by molecular modeling. The signal sequences statistically segregated into two classes of states. Members of one class were uncharged and conformationally compact while members of the second class each carried a 2+ electric charge and were conformationally extended. Virtual signal sequences, not found on GenBank and based upon hypothetical transversions in the third codon, had molecular characteristics intermediate to those of the two classes of actual signal sequences. The high incidence of non-synonymous substitutions (63.6%), the high transition/transversion ratio (10/1) and the results of molecular modeling all suggest that the N-terminal cleavable signal sequence is mutationally evolving more rapidly than proteins which must assume specific conformational states in the mature influenza virion.

Evolution, Molecular↗

PubMed2002

Accuracy and power of bayes prediction of amino acid sites under positive selection.

Bayes prediction quantifies uncertainty by assigning posterior probabilities. It was used to identify amino acids in a protein under recurrent diversifying selection indicated by higher nonsynonymous (d(N)) than synonymous (d(S)) substitution rates or by omega = d(N)/d(S) > 1. Parameters were estimated by maximum likelihood under a codon substitution model that assumed several classes of sites with different omega ratios. The Bayes theorem was used to calculate the posterior probabilities of each site falling into these site classes. Here, we evaluate the performance of Bayes prediction of amino acids under positive selection by computer simulation. We measured the accuracy by the proportion of predicted sites that were truly under selection and the power by the proportion of true positively selected sites that were predicted by the method. The accuracy was slightly better for longer sequences, whereas the power was largely unaffected by the increase in sequence length. Both accuracy and power were higher for medium or highly diverged sequences than for similar sequences. We found that accuracy and power were unacceptably low when data contained only a few highly similar sequences. However, sampling a large number of lineages improved the performance substantially. Even for very similar sequences, accuracy and power can be high if over 100 taxa are used in the analysis. We make the following recommendations: (1) prediction of positive selection sites is not feasible for a few closely related sequences; (2) using a large number of lineages is the best way to improve the accuracy and power of the prediction; and (3) multiple models of heterogeneous selective pressures among sites should be applied in real data analysis.

Amino Acid Substitution↗

PubMed1997

Selection on the protein-coding genes of the TBE1 family of transposable elements in the ciliates Oxytricha fallax and O. trifallax.

TBE1s are "cut-and-paste" transposable elements found in high copy number in the germline genomes of the ciliates Oxytricha fallax and O. trifallax. TBE1 "family" sequence (sequence of mixed polymerase chain reaction products generated using primers that match roughly half the TBE1s in host whole-cell DNA) was obtained from both host species. Although family sequence autoradiograms represent thousands of different elements, they are as legible as those representing corresponding sequences of a single TBE1, implying that ideal polymorphisms are rare within the genes examined. Nucleotide polymorphisms among TBE1s (indicated by ambiguities in family sequence) are far more common at third than at first or second positions of codons of genes, implying that selection has conserved the amino acid sequences of these genes in the majority of TBE1s. Portions of the transposase gene and another TBE1 gene have been sequenced from 10 individual TBE1s. None of these portions is interrupted by stop codons or frameshifts, and, for both genes, pairwise comparisons of these sequences show that nonsynonymous differences are significantly less common than synonymous differences, again implicating conservative selection Phylogenetic analysis shows that multiple divergent lineages of TBE1s have evolved under this selection within O. fallax. All these results are unexpected for cut-and-paste transposons in eukaryotic hosts: since transposase encoded by intact elements presumably acts in trans, it can duplicate mutant copies (those that do not encode functional transposase) found in the same genome, and thus no selection is expected to maintain the transposase gene. The selection demonstrated here could act at transposition (if functional TBE1s are preferentially transposed) or at the level of the host (if the host's fitness depends on functional TBE1 genes). TBE1-encoded proteins might be responsible for the precise excision of TBE1s that occurs during development of the host somatic nucleus; selection on hosts for uninterrupted somatic genes would then translate into selection for TBE1 protein-coding competence. We suggest a method for distinguishing between these two classes of explanations by finding and analyzing divergent alleles of ancestral transposable element insertions.

Animals↗

PubMed1998

Evolutionary rates for tuf genes in endosymbionts of aphids.

The gene encoding elongation factor Tu (tuf) in aphid endosymbionts (genus Buchnera) evolves at rates of 1.3 x 10(-10) to 2.5 x 10(-10) nonsynonymous substitutions and 3.9 x 10(-9) to 8.0 x 10(-9) synonymous substitutions per position per year. These rates, which are at present among the most reliable substitution rates for protein-coding genes of bacteria, have been obtained by calibrating the nodes in the phylogenetic tree produced from the Buchnera EF-Tu sequences using divergence times for the corresponding ancestral aphid hosts. We also present data suggesting that the rates of nonsynonymous substitutions are significantly higher in the endosymbiont lineages than in the closely related free-living bacteria Escherichia coli and Salmonella typhimurium. Synonymous substitution rates for Buchnera approximate estimated mutation rates for E. coli and S. typhimurium, as expected if synonymous changes act as neutral mutations in Buchnera. We relate the observed differences in substitution frequencies to the absence of selective codon preferences in Buchnera and to the influence of Muller's ratchet on small asexual populations.

Amino Acid Sequence↗