Whole-gene positive selection, elevated synonymous substitution rates, duplication, and indel evolution of the chloroplast clpP1 gene

Background

Synonymous DNA substitution rates in the plant chloroplast genome are generally relatively slow and lineage dependent. Non-synonymous rates are usually even slower due to purifying selection acting on the genes. Positive selection is expected to speed up non-synonymous substitution rates, whereas synonymous rates are expected to be unaffected. Until recently, positive selection has seldom been observed in chloroplast genes, and large-scale structural rearrangements leading to gene duplications are hitherto supposed to be rare.

Methodology/Principle Findings

We found high substitution rates in the exons of the plastid clpP1 gene in Oenothera (the Evening Primrose family) and three separate lineages in the tribe Sileneae (Caryophyllaceae, the Carnation family). Introns have been lost in some of the lineages, but where present, the intron sequences have substitution rates similar to those found in other introns of their genomes. The elevated substitution rates of clpP1 are associated with statistically significant whole-gene positive selection in three branches of the phylogeny. In two of the lineages we found multiple copies of the gene. Neighboring genes present in the duplicated fragments do not show signs of elevated substitution rates or positive selection. Although non-synonymous substitutions account for most of the increase in substitution rates, synonymous rates are also markedly elevated in some lineages. Whereas plant clpP1 genes experiencing negative (purifying) selection are characterized by having very conserved lengths, genes under positive selection often have large insertions of more or less repetitive amino acid sequence motifs.

Conclusions/Significance

We found positive selection of the clpP1 gene in various plant lineages to correlated with repeated duplication of the clpP1 gene and surrounding regions, repetitive amino acid sequences, and increase in synonymous substitution rates. The present study sheds light on the controversial issue of whether negative or positive selection is to be expected after gene duplications by providing evidence for the latter alternative. The observed increase in synonymous substitution rates in some of the lineages indicates that the detection of positive selection may be obscured under such circumstances. Future studies are required to explore the functional significance of the large inserted repeated amino acid motifs, as well as the possibility that synonymous substitution rates may be affected by positive selection.     

Molecular evolution of nuclear genes in Cupressacea, a group of conifer trees

We surveyed the molecular evolutionary characteristics of 11 nuclear genes from 10 conifer trees belonging to the Taxodioideae, the Cupressoideae, and the Sequoioideae. Comparisons of substitution rates among the lineages indicated that the synonymous substitution rates of the Cupressoideae lineage were higher than those of the Taxodioideae. This result parallels the pattern previously found in plastid genes. Likelihood-ratio tests showed that the nonsynonymous-synonymous rate ratio did not change significantly among lineages. In addition, after adjustments for lineage effects, the dispersion indices of synonymous and nonsynonymous substitutions were considerably reduced, and the latter was close to 1. These results indicated that the acceleration of evolutionary rates in the Cupressoideae lineage occurred in both the nuclear and plastid genomes, and that generally, this lineage effect affected synonymous and nonsynonymous substitutions similarly. We also investigated the relationship of synonymous substitution rates with the nonsynonymous substitution rate, base composition, and codon bias in each lineage. Synonymous substitution rates were positively correlated with nonsynonymous substitution rates and GC content at third codon positions, but synonymous substitution rates were not correlated with codon bias. Finally, we tested the possibility of positive selection at the protein level, using maximum likelihood models, assuming heterogeneous nonsynonymous-synonymous rate ratios among codon (amino acid) sites. Although we did not detect strong evidence of positively selected codon sites, the analysis suggested that significant variation in nonsynonymous-synonymous rate ratio exists among the sites. The most likely sites for action of positive selection were found in the ferredoxin gene, which is an important component of the apparatus for photosynthesis.     

Partitioning the variation in mammalian substitution rates

We have used analysis of variance to partition the variation in synonymous and amino acid substitution rates between three effects (gene, lineage, and a gene-by-lineage interaction) in mammalian nuclear and mitochondrial genes. We find that gene effects are stronger for amino acid substitution rates than for synonymous substitution rates and that lineage effects are stronger for synonymous substitution rates than for amino acid substitution rates. Gene-by-lineage interactions, equivalent to overdispersion corrected for lineage effects, are found in amino acid substitutions but not in synonymous substitutions. The variance in the ratio of amino acid and synonymous substitution rates is dominated by gene effects, but there is also a significant gene-by-lineage interaction.     

Molecular evolution of the COX7A gene family in primates.

COX VIIa is one of 10 nuclear-encoded subunits of the COX holoenzyme, and one of three that have isoforms with tissue-specific differences in expression. Analysis of nucleotide substitution rates revealed an accelerated rate of nonsynonymous substitutions relative to that of synonymous substitutions for the heart isoform gene (COX7AH) in six primate lineages. Rate accelerations have been noted for four other COX-related genes in this time period, suggesting that the COX holoenzyme has experienced an episode of adaptive evolution. A third member of the gene family, COX7AR, has recently been described. Although its function is currently unknown, low nonsynonymous substitution/synonymous substitution (N/S) ratios in mammalian evolution suggest that COX7AR is of functional importance. When the COX7A isoforms were divided into domains, examination of nucleotide substitution rates suggested that mitochondrial targeting residues experienced an accelerated nonsynonymous substitution rate in the period following gene duplication. In contrast, paralogous comparisons of the targeting residues of each isoform show they have been relatively conserved in mammalian evolution. This pattern is consistent with the evolution of tissue-specific function.     

Implications of the Plastid Genome Sequence of Typha (Typhaceae, Poales) for Understanding Genome Evolution in Poaceae

Plastid genomes of the grasses (Poaceae) are unusual in their organization and rates of sequence evolution. There has been a recent surge in the availability of grass plastid genome sequences, but a comprehensive comparative analysis of genome evolution has not been performed that includes any related families in the Poales. We report on the plastid genome of Typha latifolia, the first non-grass Poales sequenced to date, and we present comparisons of genome organization and sequence evolution within Poales. Our results confirm that grass plastid genomes exhibit acceleration in both genomic rearrangements and nucleotide substitutions. Poaceae have multiple structural rearrangements, including three inversions, three genes losses (accD, ycf1, ycf2), intron losses in two genes (clpP, rpoC1), and expansion of the inverted repeat (IR) into both large and small single-copy regions. These rearrangements are restricted to the Poaceae, and IR expansion into the small single-copy region correlates with the phylogeny of the family. Comparisons of 73 protein-coding genes for 47 angiosperms including nine Poaceae genera confirm that the branch leading to Poaceae has significantly accelerated rates of change relative to other monocots and angiosperms. Furthermore, rates of sequence evolution within grasses are lower, indicating a deceleration during diversification of the family. Overall there is a strong correlation between accelerated rates of genomic rearrangements and nucleotide substitutions in Poaceae, a phenomenon that has been noted recently throughout angiosperms. The cause of the correlation is unknown, but faulty DNA repair has been suggested in other systems including bacterial and animal mitochondrial genomes.     

Correlated rates of synonymous site evolution across plant genomes

Synonymous substitution rates have been shown to vary among evolutionary lineages of both nuclear and organellar genes across a broad range of taxonomic groups. In animals, rate heterogeneity does not appear to be correlated across nuclear and mitochondrial genes. In this paper, we contrast substitution rates in two plant groups and show that grasses evolve more rapidly than palms at synonymous sites in a mitochondrial, a nuclear, and a plastid gene. Furthermore, we show that the relative rates of synonymous substitution between grasses and palms are similar at the three loci. The correlation in synonymous substitution rates across genes is particularly striking because the three genes evolve at very different absolute rates. In contrast, relative rates of nonsynonymous substitution are not conserved among the three genes.      

Further Examples of Evolution by Gene Duplication Revealed through DNA Sequence Comparisons

To test the theory that evolution by gene duplication occurs as a result of positive Darwinian selection that accompanies the acceleration of mutant substitutions, DNA sequences of recent duplication were analyzed by estimating the numbers of synonymous and nonsynonymous substitutions. For the troponin C family, at the period of differentiation of the fast and slow isoforms, amino acid substitutions were shown to have been accelerated relative to synonymous substitutions. Comparison of the first exon of α-actin genes revealed that amino acid substitutions were accelerated when the smooth muscle, skeletal and cardiac isoforms differentiated. Analysis of members of the heat shock protein 70 gene family of mammals indicates that heat shock responsive genes including duplicated copies are evolving rapidly, contrary to the cognitive genes which have been evolutionarily conservative. For the α(1)-antitrypsin reactive center, the acceleration of amino acid substitution has been found for gene pairs of recent duplication.     

Substitution rates in a new Silene latifolia sex-linked gene, SlssX/Y

Dioecious white campion Silene latifolia has sex chromosomal sex determination, with homogametic (XX) females and heterogametic (XY) males. This species has become popular in studies of sex chromosome evolution. However, the lack of genes isolated from the X and Y chromosomes of this species is a major obstacle for such studies. Here, I report the isolation of a new sex-linked gene, Slss, with strong homology to spermidine synthase genes of other species. The new gene has homologous intact copies on the X and Y chromosomes (SlssX and SlssY, respectively). Synonymous divergence between the SlssX and SlssY genes is 4.7%, and nonsynonymous divergence is 1.4%. Isolation of a homologous gene from nondioecious S. vulgaris provided a root to the gene tree and allowed the estimation of the silent and replacement substitution rates along the SlssX and SlssY lineages. Interestingly, the Y-linked gene has higher synonymous and nonsynonymous substitution rates. The elevated synonymous rate in the SlssY gene, compared with SlssX, confirms our previous suggestion that the S. latifolia Y chromosome has a higher mutation rate, compared with the X chromosome. When differences in silent substitution rate are taken into account, the Y-linked gene still demonstrates significantly faster accumulation of nonsynonymous substitutions, which is consistent with the theoretical prediction of relaxed purifying selection in Y-linked genes, leading to the accumulation of nonsynonymous substitutions and genetic degeneration of the Y-linked genes.     

Relative rates of synonymous substitutions in the mitochondrial, chloroplast and nuclear genomes of seed plants

Previous studies have estimated that, in angiosperms, the synonymous substitution rate of chloroplast genes is three times higher than that of mitochondrial genes and that of nuclear genes is twelve times higher than that of mitochondrial genes. Here we used 12 genes in 27 seed plant species to investigate whether these relative rates of substitutions are common to diverse seed plant groups. We find that the overall relative rate of synonymous substitutions of mitochondrial, chloroplast and nuclear genes of all seed plants is 1:3:10, that these ratios are 1:2:4 in gymnosperms but 1:3:16 in angiosperms and that they go up to 1:3:20 in basal angiosperms. Our results show that the mitochondrial, chloroplast and nuclear genomes of seed plant groups have different synonymous substitutions rates, that these rates are different in different seed plant groups and that gymnosperms have smaller ratios than angiosperms.     

Evaluation of whether accelerated protein evolution in chordates has occurred before, after, or simultaneously with gene duplication

Gene duplication and loss are predicted to be at least of the order of the substitution rate and are key contributors to the development of novel gene function and overall genome evolution. Although it has been established that proteins evolve more rapidly after gene duplication, we were interested in testing to what extent this reflects causation or association. Therefore, we investigated the rate of evolution prior to gene duplication in chordates. Two patterns emerged; firstly, branches, which are both preceded by a duplication and followed by a duplication, display an elevated rate of amino acid replacement. This is reflected in the ratio of nonsynonymous to synonymous substitution (mean nonsynonymous to synonymous nucleotide substitution rate ratio [Ka:Ks]) of 0.44 compared with branches preceded by and followed by a speciation (mean Ka:Ks of 0.23). The observed patterns suggest that there can be simultaneous alteration in the selection pressures on both gene duplication and amino acid replacement, which may be consistent with co-occurring increases in positive selection, or alternatively with concurrent relaxation of purifying selection. The pattern is largely, but perhaps not completely, explained by the existence of certain families that have elevated rates of both gene duplication and amino acid replacement. Secondly, we observed accelerated amino acid replacement prior to duplication (mean Ka:Ks for postspeciation preduplication branches was 0.27). In some cases, this could reflect adaptive changes in protein function precipitating a gene duplication event. In conclusion, the circumstances surrounding the birth of new proteins may frequently involve a simultaneous change in selection pressures on both gene-copy number and amino acid replacement. More precise modeling of the relative importance of preduplication, postduplication, and simultaneous amino acid replacement will require larger and denser genomic data sets from multiple species, allowing simultaneous estimation of lineage-specific fluctuations in mutation rates and adaptive constraints.     

Selective processes and adaptive evolution of the cytochrome b gene in salamanders of the genus Salamandrella

Sequence analysis of the cytochrome b gene fragment in the salamanders of the genus Salamandrella, Siberian salamander and Schrenk salamander was performed with the purpose to elucidate the effect of natural selection on the evolution of mitochondrial DNA (mtDNA) in these species. It was demonstrated that despite of notable influence of negative selection (expressed as very low dN/dS values), speciation and intraspecific divergence in salamanders was accompanied by the appearance of radical amino acid substitutions, caused by the influence of positive (directional) selection. To examine the evolutionary pattern of synonymous mtDNA sites, distribution of conservative and non-conservative substitutions was analyzed. The rates of conservative and non-conservative substitutions were nearly equal, pointing to neutrality of mutation process at synonymous mtDNA sites of salamanders. Analysis of conservative and non-conservative synonymous substitution distributions in different parts of phylogenetic trees showed that the differences between the synonymous groups compared were statistically significant only in one phylogenetic group of Siberian salamander (haplogroup C) (P = 0.02). In the group of single substitutions, located at terminal phylogenetic branches of Siberian salamanders from group C, increased rate of conservative substitutions was observed. Based on these findings, it was suggested that selective processes could have an influence on the formation of the synonymous substitution profile in the Siberian salamander mtDNA fragment examined.     

Explaining the Imperfection of the Molecular Clock of Hominid Mitochondria

The molecular clock of mitochondrial DNA has been extensively used to date various genetic events. However, its substitution rate among humans appears to be higher than rates inferred from human-chimpanzee comparisons, limiting the potential of interspecies clock calibrations for intraspecific dating. It is not well understood how and why the substitution rate accelerates. We have analyzed a phylogenetic tree of 3057 publicly available human mitochondrial DNA coding region sequences for changes in the ratios of mutations belonging to different functional classes. The proportion of non-synonymous and RNA genes substitutions has reduced over hundreds of thousands of years. The highest mutation ratios corresponding to fast acceleration in the apparent substitution rate of the coding sequence have occurred after the end of the Last Ice Age. We recalibrate the molecular clock of human mtDNA as 7990 years per synonymous mutation over the mitochondrial genome. However, the distribution of substitutions at synonymous sites in human data significantly departs from a model assuming a single rate parameter and implies at least 3 different subclasses of sites. Neutral model with 3 synonymous substitution rates can explain most, if not all, of the apparent molecular clock difference between the intra- and interspecies levels. Our findings imply the sluggishness of purifying selection in removing the slightly deleterious mutations from the human as well as the Neandertal and chimpanzee populations. However, for humans, the weakness of purifying selection has been further exacerbated by the population expansions associated with the out-of Africa migration and the end of the Last Ice Age.     

Hotspots of Biased Nucleotide Substitutions in Human Genes

Genes that have experienced accelerated evolutionary rates on the human lineage during recent evolution are candidates for involvement in human-specific adaptations. To determine the forces that cause increased evolutionary rates in certain genes, we analyzed alignments of 10,238 human genes to their orthologues in chimpanzee and macaque. Using a likelihood ratio test, we identified protein-coding sequences with an accelerated rate of base substitutions along the human lineage. Exons evolving at a fast rate in humans have a significant tendency to contain clusters of AT-to-GC (weak-to-strong) biased substitutions. This pattern is also observed in noncoding sequence flanking rapidly evolving exons. Accelerated exons occur in regions with elevated male recombination rates and exhibit an excess of nonsynonymous substitutions relative to the genomic average. We next analyzed genes with significantly elevated ratios of nonsynonymous to synonymous rates of base substitution (d_N/d_S) along the human lineage, and those with an excess of amino acid replacement substitutions relative to human polymorphism. These genes also show evidence of clusters of weak-to-strong biased substitutions. These findings indicate that a recombination-associated process, such as biased gene conversion (BGC), is driving fixation of GC alleles in the human genome. This process can lead to accelerated evolution in coding sequences and excess amino acid replacement substitutions, thereby generating significant results for tests of positive selection.     

Detecting positive selection within genomes: the problem of biased gene conversion

The identification of loci influenced by positive selection is a major goal of evolutionary genetics. A popular approach is to perform scans of alignments on a genome-wide scale in order to find regions evolving at accelerated rates on a particular branch of a phylogenetic tree. However, positive selection is not the only process that can lead to accelerated evolution. Notably, GC-biased gene conversion (gBGC) is a recombination-associated process that results in the biased fixation of G and C nucleotides. This process can potentially generate bursts of nucleotide substitutions within hotspots of meiotic recombination. Here, we analyse the results of a scan for positive selection on genes on branches across the primate phylogeny. We show that genes identified as targets of positive selection have a significant tendency to exhibit the genomic signature of gBGC. Using a maximum-likelihood framework, we estimate that more than 20 per cent of cases of significantly elevated non-synonymous to synonymous substitution rates ratio (d_N/d_S), particularly in shorter branches, could be due to gBGC. We demonstrate that in some cases, gBGC can lead to very high d_N/d_S (more than 2). Our results indicate that gBGC significantly affects the evolution of coding sequences in primates, often leading to patterns of evolution that can be mistaken for positive selection.     

Similarity of Synonymous Substitution Rates Across Mammalian Genomes

Given that a gene has a high (or low) synonymous substitution rate in one mammalian species, will it also have a high (or low) synonymous substitution rate in another mammalian species? Such similarities in the rate of synonymous substitution can reveal both selective pressures and neutral processes acting on mammalian gene sequences; however, the existence of such an effect has been a matter of disagreement. We resolve whether such synonymous substitution rate similarities exist using 7462 ortholog triplets aligned across rat, mouse, and human, a dataset two orders of magnitude larger than previous studies. We find that a gene’s synonymous substitution rate in the rat-mouse branch of the phylogeny is correlated with its rate in the branch connecting human and the rat-mouse ancestor. We confirm this for several different measures of synonymous substitution rate, including corrections for base composition and CpG dinucleotides, and we verify the results in the larger mouse-human-rat-dog phylogeny. This similarity of rates is most apparent for genes in which synonymous sites are well conserved across species, suggesting that a significant component of the effect is due to purifying selection. We observe rate correlations at a resolution as fine as a few hundred kilobases, and the genes with the most similar synonymous substitution rates are enriched for regulatory functions. Genes with above-average substitution rates also exhibit significant, though somewhat weaker, rate correlations, suggesting that some neutral processes may have persisted in the phylogeny as well.     

Codon bias differentiates between the duplicated amylase loci following gene duplication in Drosophila

We examined the pattern of synonymous substitutions in the duplicated Amylase (Amy) genes (called the Amy1- and Amy3-type genes, respectively) in the Drosophila montium species subgroup. The GC content at the third synonymous codon sites of the Amy1-type genes was higher than that of the Amy3-type genes, while the GC content in the 5'-flanking region was the same in both genes. This suggests that the difference in the GC content at third synonymous sites between the duplicated genes is not due to the temporal or regional changes in mutation bias. We inferred the direction of synonymous substitutions along branches of a phylogeny. In most lineages, there were more synonymous substitutions from G/C (G or C) to A/T (A or T) than from A/T to G/C. However, in one lineage leading to the Amy1-type genes, which is immediately after gene duplication but before speciation of the montium species, synonymous substitutions from A/T to G/C were predominant. According to a simple model of synonymous DNA evolution in which major codons are selectively advantageous within each codon family, we estimated the selection intensity for specific lineages in a phylogeny on the basis of inferred patterns of synonymous substitutions. Our result suggested that the difference in GC content at synonymous sites between the two Amy-type genes was due to the change of selection intensity immediately after gene duplication but before speciation of the montium species.     

Nucleic acid composition,codon usage,and the rate of synonymous substitution in protein-coding genes

Summary Based on the rates of synonymous substitution in 42 protein-codin gene pairs from rat and human, a correlation is shown to exist between the frequency of the nucleotides in all positions of the codon and the synonymous substitution rate. The correlation coefficients were positive for A and T and negative for C and G. This means that AT-rich genes accumulate more synonymous substitutions than GC-rich genes. Biased patterns of mutation could not account for this phenomenon. Thus, the variation in synonymous substitution rates and the resulting unequal codon usage must be the consequence of selection against A and T in synonymous positions. Most of the varition in rates of synonymous substitution can be explained by the nucleotide composition in synonymous positions. Codon-anticodon interactions, dinucleotide frequencies, and contextual factors influence neither the rates of synonymous substitution nor codon usage. Interestingly, the nucleotide in the second position of codons (always a nonsynonymous position) was found to affect the rate of synonymous substitution. This finding links the rate of nonsynonymous substitution with the synonymous rate. Consequently, highly conservative proteins are expected to be encoded by genes that evolve slowly in terms of synonymous substitutions, and are consequently highly biased in their codon usage.     

Evolution of the Sry genes

Existing DNA sequence data on the Sry gene, the mammalian sex- determining locus in the Y chromosome, were analyzed for primates, rodents, and bovids. In all three taxonomic groups, the terminal sequences evolved faster than the HMG (high mobility group) boxes, and this applies both to synonymous (Ks) and nonsynonymous (Ka) nucleotide substitutions. Similar intragenic correlation between synonymous and nonsynonymous substitution rates was not found either in other mammalian genes that contain a conservative box (Sox, Msx) or in the MADS-box genes of plants. The rate of nonsynonymous substitutions exceeds significantly that of synonymous substitutions in the terminal Sry sequences of apes. We did not find good support for the hypothesis that the high evolutionary rate of Sry would be associated with a promiscuous mating system.