Determining the Relative Rates of Change for Prokaryotic and Eukaryotic Proteins with Anciently Duplicated Paralogs

The relative rates of change for eight sets of ubiquitous proteins were determined by a test in which anciently duplicated paralogs are used to root the universal tree and distances are calculated between each taxonomic group and the last common ancestor. The sets included ATPase subunits, elongation factors, signal recognition particle and its receptor, three sets of tRNA synthetases, transcarbamoylases, and an internal duplication in carbamoyl phosphate synthase. In each case phylogenetic trees were constructed and the distances determined for all pairs. Taken over the period of time since their last common ancestor, average evolutionary rates are remarkably similar for Bacteria and Eukarya, but Archaea exhibit a significantly slower average rate. Received: 30 December 1999 / Accepted: 5 April 2000     

The Root of the Universal Tree of Life Inferred from Anciently Duplicated Genes Encoding Components of the Protein-Targeting Machinery

The key protein of the signal recognition particle (termed SRP54 for Eucarya and Ffh for Bacteria) and the protein (termed SRα for Eucarya and Ftsy for bacteria) involved in the recognition and binding of the ribosome SRP nascent polypeptide complex are the products of an ancient gene duplication that appears to predate the divergence of all extant taxa. The paralogy of the genes encoding the two proteins (both of which are GTP triphosphatases) is argued by obvious sequence similarities between the N-terminal half of SRP54(Ffh) and the C-terminal half of SRα(Ftsy). This enables a universal phylogeny based on either protein to be rooted using the second protein as an outgroup. Phylogenetic trees inferred by various methods from an alignment (220 amino acid positions) of the shared SRP54(Ffh) and SRα(Ftsy) regions generate two reciprocally rooted universal trees corresponding to the two genes. The root of both trees is firmly positioned between Bacteria and Archaea/Eucarya, thus providing strong support for the notion (Iwabe et al. 1989; Gogarten et al. 1989) that the first bifurcation in the tree of life separated the lineage leading to Bacteria from a common ancestor to Archaea and Eucarya. None of the gene trees inferred from the two paralogues support a paraphyletic Archaea with the crenarchaeota as a sister group to Eucarya. Received: 19 March 1998 / Accepted: 5 June 1998     

Evolution Rates of Genes on Leading and Lagging DNA Strands

One of the main causes of bacterial chromosome asymmetry is replication-associated mutational pressure. Different rates of nucleotide substitution accumulation on leading and lagging strands implicate qualitative and quantitative differences in the accumulation of mutations in protein coding sequences lying on different DNA strands. We show that the divergence rate of orthologs situated on leading strands is lower than the divergence rate of those situated on lagging strands. The ratio of the mutation accumulation rate for sequences lying on lagging strands to that of sequences lying on leading strands is rather stable and time-independent. The divergence rate of sequences which changed their positions, with respect to the direction of replication fork movement, is not stable—sequences which have recently changed their positions are the most prone to mutation accumulation. This effect may influence estimations of evolutionary distances between species and the topology of phylogenetic trees. Received: 24 July 2000 / Accepted: 16 January 2001     

High Evolutionary Rates in Nuclear Genes of Squamates

We compared nonsynonymous substitution rates (Ka) of nuclear coding genes between four major groups of living sauropsids (reptiles): birds, squamates, crocodiles, and turtles. Since only 9 orthologous genes are known in all the four taxonomic groups, we searched for orthologous genes known in chicken and at least one of any representative of poikilotherm sauropsids. Thus, we analyzed three additional data sets: 28 genes identified in chicken and various squamates, 24 genes identified in chicken and crocodilians, and 20 genes identified in chicken and turtles. To compare nonsynonymous substitution rates between all lineages of sauropsids, we used the relative-rate test with human genes as the outgroup. We show that 22/28 nuclear coding genes of squamates, especially snakes (15/16), have an higher evolutionary rate than those in chicken (in mean, 30–40% faster). However, no such difference is detected between crocodiles, turtles and chicken. Higher substitution rate in squamates nuclear coding genes than in chicken, and probably than in other sauropsids, could explain some of the difficulties in resolving the molecular phylogeny of reptiles. Received: 5 July 2000 / Accepted: 13 February 2001     

Ascidian and Amphioxus Adh Genes Correlate Functional and Molecular Features of the ADH Family Expansion During Vertebrate Evolution

The alcohol dehydrogenase (ADH) family has evolved into at least eight ADH classes during vertebrate evolution. We have characterized three prevertebrate forms of the parent enzyme of this family, including one from an urochordate (Ciona intestinalis) and two from cephalochordates (Branchiostoma floridae and Branchiostoma lanceolatum). An evolutionary analysis of the family was performed gathering data from protein and gene structures, exon–intron distribution, and functional features through chordate lines. Our data strongly support that the ADH family expansion occurred 500 million years ago, after the cephalochordate/vertebrate split, probably in the gnathostome subphylum line of the vertebrates. Evolutionary rates differ between the ancestral, ADH3 (glutathione-dependent formaldehyde dehydrogenase), and the emerging forms, including the classical alcohol dehydrogenase, ADH1, which has an evolutionary rate 3.6-fold that of the ADH3 form. Phylogenetic analysis and chromosomal mapping of the vertebrate Adh gene cluster suggest that family expansion took place by tandem duplications, probably concurrent with the extensive isoform burst observed before the fish/tetrapode split, rather than through the large-scale genome duplications also postulated in early vertebrate evolution. The absence of multifunctionality in lower chordate ADHs and the structures compared argue in favor of the acquisition of new functions in vertebrate ADH classes. Finally, comparison between B. floridae and B. lanceolatum Adhs provides the first estimate for a cephalochordate speciation, 190 million years ago, probably concomitant with the beginning of the drifting of major land masses from the Pangea. Received: 10 April 2001 / Accepted: 23 May 2001     

Rates of Conservative and Radical Nonsynonymous Nucleotide Substitutions in Mammalian Nuclear Genes

To understand the process and mechanism of protein evolution, it is important to know what types of amino acid substitutions are more likely to be under selection and what types are mostly neutral. An amino acid substitution can be classified as either conservative or radical, depending on whether it involves a change in a certain physicochemical property of the amino acid. Assuming Kimura's two-parameter model of nucleotide substitution, I present a method for computing the numbers of conservative and radical nonsynonymous (amino acid altering) nucleotide substitutions per site and estimate these rates for 47 nuclear genes from mammals. The results are as follows. (1) The average radical/conservative rate ratio is 0.81 for charge changes, 0.85 for polarity changes, and 0.49 when both polarity and volume changes are considered. (2) The radical/conservative rate ratio is positively correlated with the nonsynonymous/synonymous rate ratio for charge changes or when both polarity and volume changes are considered. (3) Both the conservative/synonymous rate ratio and the radical/synonymous rate ratio are lower in the rodent lineage than in the primate or artiodactyl lineage, suggesting more intense purifying selection in the rodent lineage, for both conservative and radical nonsynonymous substitutions. (4) Neglecting transition/transversion bias would cause an underestimation of both radical and conservative rates and the ratio thereof. (5) Transversions induce more dramatic genetic alternations than transitions in that transversions produce more amino acid altering changes and among which, more radical changes. Received: 6 April 1999 / Accepted: 16 August 1999     

Synonymous Nucleotide Divergence and Saturation: Effects of Site-Specific Variations in Codon Bias and Mutation Rates

The synonymous divergence between Escherichia coli and Salmonella typhimurium is explained in a model where there is a large variation between mutation rates at different nucleotide sites in the genome. The model is based on the experimental observation that spontaneous mutation rates can vary over several orders of magnitude at different sites in a gene. Such site-specific variation must be taken into account when studying synonymous divergence and will result in an apparent saturation below the level expected from an assumption of uniform rates. Recently, it has been suggested that codon preference in enterobacteria has a very large site-specific variation and that the synonymous divergence between different species, e.g., E. coli and Salmonella, is saturated. In the present communication it is shown that when site-specific variation in mutation rates is introduced, there is no need to invoke assumptions of saturation and a large variability in codon preference. The same rate variation will also bring average mutation rates as estimated from synonymous sequence divergence into numerical agreement with experimental values. Received: 10 July 1998 / Accepted: 20 August 1998     

Evolution of Duplicated <Emphasis Type="BoldItalic">reggie</Emphasis> Genes in Zebrafish and Goldfish

Invertebrates, tetrapod vertebrates, and fish might be expected to differ in their number of gene copies, possibly due the occurrence of genome duplication events during animal evolution. Reggie (flotillin) genes code for membrane-associated proteins involved in growth signaling in developing and regenerating axons. Until now, there appeared to be only two reggie genes in fruitflies, mammals, and fish. The aim of this research was to search for additional copies of reggie genes in fishes, since a genome duplication might have increased the gene copy number in this group. We report the presence of up to four distinct reggie genes (two reggie-1 and two reggie-2 genes) in the genomes of zebrafish and goldfish. Phylogenetic analyses show that the zebrafish and goldfish sequence pairs are orthologous, and that the additional copies could have arisen through a genome duplication in a common ancestor of bony fish. The presence of novel reggie mRNAs in fish embryos indicates that the newly discovered gene copies are transcribed and possibly expressed in the developing and regenerating nervous system. The intron/exon boundaries of the new fish genes characterized here correspond with those of human genes, both in location and phase. An evolutionary scenario for the evolution of reggie intron-exon structure, where loss of introns appears to be a distinctive trait in invertebrate reggie genes, is presented. Received: 24 January 2001 / Accepted: 27 July 2001     

Evidence for Gene Conversion Between Tandemly Duplicated Cytoplasmic Actin Genes of Helicoverpa armigera (Lepidoptera:Noctuidae)

Tandemly duplicated actin genes have been isolated from a Helicoverpa armigera genomic library. Sequence comparisons with actin genes from other species suggest they encode cytoplasmic actins, being most closely related to the Bombyx mori A3 actin gene. The duplicated H. armigera actin genes, termed A3a and A3b, share 98.3% nucleotide sequence identity over their entire putative coding region. Analysis of the distribution of nucleotide differences shows the first 763 bp are identical between the two coding regions, with the 18 nucleotide changes occurring in the remaining 366 bp. This observation suggests a gene conversion event has taken place between the duplicated H. armigera A3a and A3b actin genes. Translation of the open-reading frames indicates the products of these genes are identical, apart from a single amino acid difference at codon 273. Polymerase chain reaction and northern blot analysis have shown both H. armigera A3a and A3b genes are expressed during pupal development and in the brain of newly eclosed adults. A region 5′ of the H. armigera A3a actin gene start codon has been identified which contains regulatory sequences commonly found in the promoter region of actin genes, including TATA, CAAT, and CArG motifs. Received: 10 January 1996 / Accepted: 12 March 1996     

Evolution of Rh Blood Group Genes Have Experienced Gene Conversions and Positive Selection

There are two tightly linked loci (D and CE) for the human Rh blood group. Their gene products are membrane proteins having 12 transmembrane domains and form a complex with Rh50 glycoprotein on erythrocytes. We constructed phylogenetic networks of human and nonhuman primate Rh genes, and the network patterns suggested the occurrences of gene conversions. We therefore used a modified site-by-site reconstruction method by using two assumed gene trees and detected 9 or 11 converted regions. After eliminating the effect of gene conversions, we estimated numbers of nonsynonymous and synonymous substitutions for each branch of both trees. Whichever gene tree we selected the branch connecting hominoids and Old World monkeys showed significantly higher nonsynonymous than synonymous substitutions, an indication of positive selection. Many other branches also showed higher nonsynonymous than synonymous substitutions; this suggests that the Rh genes have experienced some kind of positive selection. Received: 16 March 1999 / Accepted: 17 June 1999     

Molecular Evolution of Calmodulin and Calmodulin-Like Genes in the Cephalochordate Branchiostoma

Calmodulin is a calcium-binding EF-hand protein that is an activator of many enzymes as well as ion pumps and channels. Due to its multiple targets and its central role in the cell, understanding the evolutionary history of calmodulin genes should provide insights into the origin of genetic complexity in eukaryotes. We have previously isolated and characterized a calmodulin gene from the early-diverging chordate Branchiostoma lanceolatum (CaM1). In this paper, we report the existence of a second calmodulin gene (CaM2) as well as two CaM-like genomic fragments (CaML-2, CaML-3) in B. lanceolatum and a CaM2 and three CaM-like genes (CaML-1, CaML-2, CaML-3) in B. floridae. The CaM-like genes were isolated using low-stringency PCR. Surprisingly, the nucleotide sequences of the B. lanceolatum CaM1 and CaM2 cDNAs differ by 19.3%. Moreover, the CaM2 protein differs at two positions from the amino acid sequence of CaM1; the latter is identical to calmodulins in Drosophila melanogaster, the mollusc Aplysia californica, and the tunicate Halocynthia roretzi. The two B. lanceolatum CaM-like genes are more closely related to the CaM2 than to the CaM1 gene. This relationship is supported by the phylogenetic analyses and the identical exon/intron organization of these three genes, a relationship unique among animal CaM sequences. These data demonstrate the existence of a CaM multigene family in the cephalochordate Branchiostoma, which may have evolved independently from the multigene family in vertebrates. Received: 2 November 1999 / Accepted: 25 April 2000     

Clustering of Tissue-Specific Genes Underlies Much of the Similarity in Rates of Protein Evolution of Linked Genes

Are genes nonrandomly distributed around the genome and might this explain why it was found that, in the mouse genome, proteins of linked genes evolve at similar rates? Anecdotal evidence suggests that the similarity of expression of linked genes might, in part, explain the similarity in their rates of evolution. Immune system genes, for example, are known to evolve at a high rate and sometimes cluster in the genome. Here we develop methods for statistical tests of similarity of expression of linked genes and report that there is a significant tendency for genes of similar expression breadth to be linked. Significantly, when we exclude tissue specific genes from our sample, the similarity in rates of protein evolution of linked genes is greatly diminished, if not abolished. This diminution is not a sampling artifact. In contrast, while half of the immune genes in our sample reside in 1 of 10 immune clusters in the mouse genome, this clustering appears not to affect the extent of local similarity in rates of evolution. The distribution of placentally expressed genes, in contrast, does have an effect.     

Positive Darwinian Selection Promotes Heterogeneity Among Members of the Antifreeze Protein Multigene Family

A variety of organisms have independently evolved proteins exhibiting antifreeze activity that allows survival at subfreezing temperatures. The antifreeze proteins (AFPs) bind ice nuclei and depress the freezing point by a noncolligative absorption–inhibition mechanism. Many organisms have a heterogeneous suite of AFPs with variation in primary sequence between paralogous loci. Here, we demonstrate that the diversification of the AFP paralogues is promoted by positive Darwinian selection in two independently evolved AFPs from fish and beetle. First, we demonstrate an elevated rate of nonsynonymous substitutions compared to synonymous substitutions in the mature protein coding region. Second, we perform phylogeny-based tests of selection to demonstrate a subset of codons is subjected to positive selection. When mapped onto the three-dimensional structure of the fish antifreeze type III antifreeze structure, these codons correspond to amino acid positions that surround but do not interrupt the putative ice-binding surface. The selective agent may be related to efficient binding to diverse ice surfaces or some other aspect of AFP function. Received: 27 February 2001 / Accepted: 12 September 2001     

Detecting Changes in the Functional Constraints of Paralogous Genes

We describe a new procedure to determine whether regional alterations in the evolutionary constraints imposed on paralogous proteins have occurred. We used as models the A and B (alternatively called α and β) subunits of V/F/A-ATPases, originated by a gene duplication more than 3 billion years ago. Changes associated to three major splits (eubacteria versus Archaea-eukaryotes; Archaea versus eukaryotes; and among free-living bacteria and symbiotic mitochondria) were studied. Only in the first case, when we compared eubacterial or mitochondrial F-ATPases versus eukaryotic vacuolar V-ATPases or archaeal A-ATPases, constraint changes were observed. Modifications in the degree of regional constraining were not detected for the other two types of comparisons (V-ATPases versus A-ATPases and within F-ATPases, respectively). When the rates of evolution of the two subunits were compared, it was found that F-ATPases regulatory subunits evolved faster than catalytic subunits, but the opposite was true for A- and V-ATPases. Our results suggest that, even for universal and essential proteins, selective constraints may be occasionally altered. On the other hand, in some cases no changes were detected after periods of more than 2.2 billion years. Received: 24 February 2000 / Accepted: 12 August 2000     

Divergent Human Y-Chromosome Microsatellite Evolution Rates

In this work, we analyze several characteristics influencing the low variability of the microsatellite DYS19 in the major founder Amerindian Y chromosome lineage containing the point mutation DYS199-T. Variation of DYS19 was compared with that of five other Y-linked tetranucleotide repeat loci (DYS389A, DYS389B, DYS390, DYS391, and DYS393) in the DYS199-T lineage. All the other microsatellites showed significantly higher levels of variability than DYS19 as measured by gene diversity and repeat number variance. Moreover, we had previously shown that DYS19 had high diversity in Brazilians and in several other populations worldwide. Thus, the slow DYS19 evolution in the DYS199-T lineage seems to be both locus and allele specific. To understand the slow DYS19 evolutionary rate, the microsatellite loci were compared according to their mapping on the Y chromosome and also on the basis of structural aspects such as the base composition of the repeat motif and flanking regions and the degree of perfection and size (repeat number) of the variable blocks. The only observed difference that might be related to the low DYS19 variability is its small average number of repeats, a value expected to be closer to the founder DYS19 allele in the DYS199-T lineage. These data were also compared to other derived Y lineages. The Tat-C lineage displayed a lower DYS19 variability correlated to a small average repeat number, while in the DYS234-G lineage, a high DYS19 variability was found associated to a larger average repeat number. This approach reveals that evolution of Y microsatellites in lineages defined by slowly evolving markers, such as point mutations, can be greatly influenced by the size (number of repeats of the variable block) of the founder allele in each microsatellite locus. Thus lineage-dating methods using microsatellite variation should be practiced with great care. Received: 7 November 1998 / Accepted: 9 April 1999     

Fast Evolution of Interleukin-2 in Mammals and Positive Selection in Ruminants

Interleukin-2 (IL-2) is a cytokine involved in induction and regulation of the immune response in mammals. There have been numerous reports about the search for IL-2 in species other than mammals, and recently an IL-2-like gene has been isolated in chicken. Using PCR, we searched for IL-2 gene sequences in a wide variety of mammals, including marsupials and monotremes, as well as in birds. Although we can readily amplify IL-2 gene fragments in placental mammals, no amplification was obtained in other species. This is best explained by very high substitution rates. This suggest that strategies to isolate IL-2 homologous genes outside mammals should involve functional assays, as for the chicken gene, and not hybridization-based techniques. Nonsynonymous substitution rates are especially high in ruminants, due to positive selection acting on regions important in term of structure-function. We suggest that, although globally similar, the immune response of various mammals is not identical, mainly at the level of cytokine-mediated regulations. Received: 27 July 1999 / Accepted: 15 April 2000     

Synonymous and Nonsynonymous Substitution Rates in Diatoms: A Comparison Between Chloroplast and Nuclear Genes

Rates of synonymous and nonsynonymous nucleotide substitutions and codon usage bias (ENC) were estimated for a number of nuclear and chloroplast genes in a sample of centric and pennate diatoms. The results suggest that DNA evolution has taken place, on an average, at a slower rate in the chloroplast genes than in the nuclear genes: a rate variation pattern similar to that observed in land plants. Synonymous substitution rates in the chloroplast genes show a negative association with the degree of codon usage bias, suggesting that genes with a higher degree of codon usage bias have evolved at a slower rate. While this relationship has been shown in both prokaryotes and multicellular eukaryotes, it has not been demonstrated before in diatoms. Received: 3 June 1998 / Accepted: 11 August 1998     

The Role of Interelement Selection in Saccharomyces cerevisiae Ty Element Evolution

Retrotransposons are mobile genetic elements that are ubiquitous components of eukaryotic genomes. The evolutionary success of retrotransposons is explained by their ability to replicate faster than the host genomes in which they reside. Elements with higher rates of genomic replication possess a selective advantage over less active elements. Retrotransposon populations, therefore, are shaped largely by selective forces acting at the genomic level between elements. To evaluate rigorously the effects of selective forces acting on retrotransposons, detailed information on the patterns of molecular variation within and between retrotransposon families is needed. The sequencing of the Saccharomyces cerevisiae genome, which includes the entire genomic complement of yeast retrotransposons, provides an unprecedented opportunity to access and analyze such data. In this study, we analyzed in detail the patterns of nucleotide variation within the open reading frames of two parental (Ty1 and Ty2) and one hybrid (Ty1/2) family of yeast retrotransposons. The pattern and distribution of nucleotide changes on the phylogenetic reconstructions of the three families of Ty elements reveal evidence of negative selection on both internal and external branches of the Ty phylogenies. These results indicate that most, if not all, Ty elements examined represent active or recently active retrotransposon lineages. We discuss the relevance of these findings with respect to the coevolutionary dynamic operating between genomic element populations and the host organisms in which they reside. Received: 5 November 1998 / Accepted: 17 March 1999     

Rates of Molecular Evolution in RNA Viruses: A Quantitative Phylogenetic Analysis

The study of rates of nucleotide substitution in RNA viruses is central to our understanding of their evolution. Herein we report a comprehensive analysis of substitution rates in 50 RNA viruses using a recently developed maximum likelihood phylogenetic method. This analysis revealed a significant relationship between genetic divergence and isolation time for an extensive array of RNA viruses, although more rate variation was usually present among lineages than would be expected under the constraints of a molecular clock. Despite the lack of a molecular clock, the range of statistically significant variation in overall substitution rates was surprisingly narrow for those viruses where a significant relationship between genetic divergence and time was found, as was the case when synonymous sites were considered alone, where the molecular clock was rejected less frequently. An analysis of the ecological and genetic factors that might explain this rate variation revealed some evidence of significantly lower substitution rates in vector-borne viruses, as well as a weak correlation between rate and genome length. Finally, a simulation study revealed that our maximum likelihood estimates of substitution rates are valid, even if the molecular clock is rejected, provided that sufficiently large data sets are analyzed. Received: 23 February 2001 / Accepted: 3 July 2001