Large number of replacement polymorphisms in rapidly evolving genes of Drosophila. Implications for genome-wide surveys of DNA polymorphism

We present a survey of nucleotide polymorphism of three novel, rapidly evolving genes in populations of Drosophila melanogaster and D. simulans. Levels of silent polymorphism are comparable to other loci, but the number of replacement polymorphisms is higher than that in most other genes surveyed in D. melanogaster and D. simulans. Tests of neutrality fail to reject neutral evolution with one exception. This concerns a gene located in a region of high recombination rate in D. simulans and in a region of low recombination rate in D. melanogaster, due to an inversion. In the latter case it shows a very low number of polymorphisms, presumably due to selective sweeps in the region. Patterns of nucleotide polymorphism suggest that most substitutions are neutral or nearly neutral and that weak (positive and purifying) selection plays a significant role in the evolution of these genes. At all three loci, purifying selection of slightly deleterious replacement mutations appears to be more efficient in D. simulans than in D. melanogaster, presumably due to different effective population sizes. Our analysis suggests that current knowledge about genome-wide patterns of nucleotide polymorphism is far from complete with respect to the types and range of nucleotide substitutions and that further analysis of differences between local populations will be required to understand the forces more completely. We note that rapidly diverging and nearly neutrally evolving genes cannot be expected only in the genome of Drosophila, but are likely to occur in large numbers also in other organisms and that their function and evolution are little understood so far.     

Evaluation of Sequence Variation and Selection in the Bindin Locus of the Red Sea Urchin, Strongylocentrotus franciscanus

Recent evidence suggests that gamete recognition proteins may be subjected to directed evolutionary pressure that enhances sequence variability. We evaluated whether diversity enhancing selection is operating on a marine invertebrate fertilization protein by examining the intraspecific DNA sequence variation of a 273-base pair region located at the 5′ end of the sperm bindin locus in 134 adult red sea urchins (Strongylocentrotus franciscanus). Bindin is a sperm recognition protein that mediates species-specific gamete interactions in sea urchins. The region of the bindin locus examined was found to be polymorphic with 14 alleles. Mean pairwise comparison of the 14 alleles indicates moderate sequence diversity (p-distance = 1.06). No evidence of diversity enhancing selection was found. It was not possible to reject the null hypothesis that the sequence variation observed in S. franciscanus bindin is a result of neutral evolution. Statistical evaluation of expected proportions of replacement and silent nucleotide substitutions, observed versus expected proportions of radical replacement substitutions, and conformance to the McDonald and Kreitman test of neutral evolution all indicate that random mutation followed by genetic drift created the polymorphisms observed in bindin. Observed frequencies were also highly similar to results expected for a neutrally evolving locus, suggesting that the polymorphism observed in the 5′ region of S. franciscanus bindin is a result of neutral evolution. Received: 19 June 1998 / Accepted: 2 August 2000     

Sequence variation in two lignin biosynthesis genes,cinnamoyl CoA reductase (CCR) and cinnamyl alcohol dehydrogenase 2 (CAD2)

Cinnamoyl CoA reductase (CCR) and cinnamyl alcohol dehydrogenase 2 (CAD2) are genes which may influence variation in lignin content and composition within plants. Sequence variation within these genes may be responsible for changes in enzyme activity and/or specificity, which could cause variation in lignin content or composition. This study examines sequence variation within these two genes in Eucalyptus globulus, an important species used in pulp and paper-making. Twenty-one single nucleotide polymorphisms (SNPs) were identified in the exons of CCR, of which nine were neutral mutations and 12 were missense mutations. Six of the missense mutations affected highly conserved amino acids within the protein sequence of CCR. Eight SNPs were identified in the CAD2 exons, six of which were neutral mutations and two which were missense mutations. One of the missense mutations affected a highly conserved amino acid within the protein sequence. In addition, 32 SNPs were identified in the CCR introns along with four insertion/deletions and two polyA length variation regions. Polymorphism affecting highly conserved amino acids may alter enzyme function and this molecular variation may be linked to variation in lignin profiles. Selecting positive alleles which produce favourable lignin profiles would be advantageous in tree breeding programs.     

Genetic diversity and evolution of reduced sulfur storage during domestication of maize

The domestication of maize has spanned a period of over 9000 years, during which time its wild relative teosinte underwent natural and artificial selection. We hypothesize that environmental conditions could have played a major role in this process. One factor of environmental variation is soil composition, which includes sulfur availability. Sulfur is reduced during photosynthesis and is used to synthesize cysteine and methionine, which drive the accumulation of δ10 (Zm00001d045937), δ18 (Zm00001d037436), β15 (Zm00001d035760), γ16 (Zm00001d005793), γ27 (Zm00001d020592), and γ50 (Zm00001d020591) zeins, representing the zein2 fraction (z2) of storage proteins in maize seeds. In this study, polymorphisms and haplotypes were detected based on six z2 genes in 60 maize and teosintes lines. Haplotypes were unevenly distributed, and abundant genetic diversity was found in teosintes. Polymorphism was highest in z2δ18, whereas for z2β15 single nucleotide polymorphism (SNP) density and insertion/deletion (indel) abundance were the lowest, indicating differential roles in seed evolution. Indels showed a clustered distribution, and most of these derived from teosintes. The indels not only led to tandem repeat polymorphisms, but also to frameshift mutations, which could also be used as null variants. In addition, neutral evolutionary tests, phylogenetic analyses, and population structures indicated that z2δ10 and z2γ50 had undergone natural selection. Indeed, a natural selection imprint could also be found with z2γ27 and z2γ16, whereas z2δ18 and z2β15 tended to be under neutral evolution. These results suggested that genetic diversity and evolution of a subset of sulfur‐rich zeins could be under environmental adaptation during maize domestication.     

Nucleotide Sequence Variation Does Not Relate to Differences in Kinetic Properties of Neutral Trehalase from the Insect Pathogenic Fungus <Emphasis Type="Italic">Metarhizium anisopliae</Emphasis>

Genetic variability in a putative virulence factor, the neutral trehalase (Ntl) gene, was examined in strains of the insect pathogenic fungi Metarhizium anisopliae and Metarhizium flavoviride by restriction fragment length polymorphism (RFLP). The Ntl gene was sequenced from four of these strains that showed dissimilar RFLP patterns. Enzyme kinetic experiments were also performed on the partially purified neutral trehalase in order to assess whether nucleotide changes in these strains related to differences in enzyme catalytic function (i.e., K _m , V _max, and K _cat). Finally, the Metarhizium strains were assessed in bioassays against waxworm larvae in order to relate nucleotide variation with Ntl enzyme kinetics and insect virulence. The greatest RFLP variation was observed with Rsa1. M. flavoviride was found to be most dissimilar in RFLP patterns when compared with the M. anisopliae strains. RFLP patterns for Ntl were diagnostic markers for previously studied genetic groups of M. anisopliae. Comparisons of Ntl sequences showed that the introns were found to be more variable (6.2%) than the exons (3.1%). Comparisons of the translated nucleotide codons showed high levels (91%) of synonymous sequence variation between strains. Another fraction of the remaining mutations was neutral, resulting in amino acid substitutions with similar functions. The neutral trehalase was partially purified by preparative isoelectric focus, revealing a single band of enzyme activity as assessed by analytical isoelectric focusing (pI ca. 5). Kinetic properties of the neutral trehalases revealed no differences between the M. anisopliae strains, while the M. flavovoride had a lower K _cat/K _m . However, there was lower virulence in one strain that showed Ntl enzyme kinetic properties that were similar to the other strains, suggesting that factors other than neutral trehalase may be responsible for delimiting virulence in this insect pathogenic fungi. Although there is nucleotide variation in genes involved in pathogenicity, this variation is mostly neutral in nature, and there is strong stabilizing selection to maintain enzyme function.     

Local Patterns of Nucleotide Polymorphism Are Highly Variable in the Selfing Species Arabidopsis thaliana

Neighboring genes predictably share similar evolutionary histories to an extent delineated by recombination. This correlation should extend across multiple linked genes in a selfing species such as Arabidopsis thaliana due to its low effective recombination rate. To test this prediction, we performed a molecular population genetics analysis of nucleotide polymorphism and divergence in chromosomal regions surrounding four low-diversity loci. Three of these loci, At1g67140, At3g03700, and TERMINAL FLOWER1 (TFL1), have been previously implicated as targets of selection and we would predict stronger correlations in polymorphism between neighboring loci due to genetic hitchhiking around these loci. The remaining locus, At1g04300, was identified in a study of linkage disequilibrium surrounding the CRYPTOCHROME2 (CRY2) locus. Although we found broad valleys of reduced nucleotide variation around two of our focal genes, At1g67140 and At3g03700, all chromosomal regions exhibited extreme variation in the patterns of polymorphism and evolution between neighboring loci. Although three of our four regions contained potential targets of selection, application of the composite-likelihood-ratio test of selection in conjunction with a goodness-of-fit test supports the selection hypothesis only for the region containing At3g03700. The degree of discordance in evolutionary histories between linked loci within each region generally correlated with estimates of recombination and linkage disequilibrium for that region, with the exception of the region containing At1g04300. We discuss the implications of these data for future population genetics analyses and genomics studies in A. thaliana. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The mean and variance of the number of segregating sites since the last hitchhiking event

 Tight linkage may cause a reduction of nucleotide diversity in a chromosomal region if an advantageous mutation appears in that region which is driven to fixation by directional selection. This process is usually called genetic hitchhiking. If selection is strong, the entire process takes place during a time period of length 2s ln (2N) that is very short relative to 2N generations [s is the selection coefficient of the advantageous mutation and N the effective diploid population size]. On the time scale of 2N generations, which is characteristic for neutral evolution, we may therefore call this process a hitchhiking event. Using coalescent methods, we analyzed a model in which a hitchhiking event occurred in a chromosomal region of zero-recombination in the past at time x. Such a hitchhiking “catastrophe” wipes out completely genetic variation that existed in a population before that time. Standing variation observed at present must therefore be due to mutations that have arisen since time point x. Assuming that all newly arising mutations are neutral, we derived expressions for the expectation, variance and also for the higher moments of the number of nucleotide sites segregating in a sample of n genes as a function of x. The result for the first moment is then used to estimate the time back to the last hitchhiking event based on DNA polymorphism data from Drosophila. Assuming that directional selection is the sole determinant of the level of genetic variation in the gene regions surveyed, we obtained estimates of x that were typically in the order of 0.1N generations. Received 14 May 1996; received in revised form 26 August 1996     

Rapid Concerted Evolution via Gene Conversion at the <Emphasis Type="BoldItalic">Drosophila hsp70</Emphasis> Genes

We analyzed nucleotide variation in the hsp70 genes of Drosophila melanogaster (five genes) and D. simulans (four genes) to characterize the homogenizing and diversifying roles of gene conversion in their evolution. Gene conversion within and between the 87A7 and 87C1 gene clusters homogenize the hsp70 coding regions; in both D. melanogaster and D. simulans, same-cluster paralogues are virtually identical, and large intercluster conversion tracts diminish 87A7/87C1 divergence. Same-cluster paralogues share many polymorphisms, consistent with frequent intracluster conversion. Shared polymorphism is highly biased toward silent variation; homogenizing conversion interacts with purifying selection. In contrast to the coding regions, some hsp70 flanking regions show conversion-mediated diversification. Strong reductions of nucleotide variability and linkage disequilibria among conversion-mediated sites in hsp70Ab and hsp70Bb alleles sampled from a single natural population are consistent with a selective sweep. Comparison of the D. melanogaster and D. simulans hsp70 genes reveals whole-family fixed differences, consistent with rapid propagation of novel mutations among duplicate genes. These results suggest that the homogenizing and diversifying roles of conversion interact to drive dynamic concerted evolution of the hsp70 genes. Received: 25 June 2001 / Accepted: 10 October 2001     

Evidence for Pervasive Adaptive Protein Evolution in Wild Mice

The relative contributions of neutral and adaptive substitutions to molecular evolution has been one of the most controversial issues in evolutionary biology for more than 40 years. The analysis of within-species nucleotide polymorphism and between-species divergence data supports a widespread role for adaptive protein evolution in certain taxa. For example, estimates of the proportion of adaptive amino acid substitutions (α) are 50% or more in enteric bacteria and Drosophila. In contrast, recent estimates of α for hominids have been at most 13%. Here, we estimate α for protein sequences of murid rodents based on nucleotide polymorphism data from multiple genes in a population of the house mouse subspecies Mus musculus castaneus, which inhabits the ancestral range of the Mus species complex and nucleotide divergence between M. m. castaneus and M. famulus or the rat. We estimate that 57% of amino acid substitutions in murids have been driven by positive selection. Hominids, therefore, are exceptional in having low apparent levels of adaptive protein evolution. The high frequency of adaptive amino acid substitutions in wild mice is consistent with their large effective population size, leading to effective natural selection at the molecular level. Effective natural selection also manifests itself as a paucity of effectively neutral nonsynonymous mutations in M. m. castaneus compared to humans.     

A Comparison of Intraspecific Patterns of DNA Sequence Variation in Mitochondrial DNA, α-Enolase,and MHC Class II B Loci in Auklets (Charadriiformes: Alcidae)

Patterns of DNA sequence variation can be used to learn about mechanisms of organismal evolution, but only if mechanisms of sequence evolution are well understood. Although theories of molecular evolution are well developed, few empirical studies have addressed patterns and mechanisms of sequence evolution in nuclear genes within species. In the present study, we compared DNA sequences among three loci with different evolutionary constraints to determine the influences of effective population size, balancing selection, and linkage on intraspecific patterns of sequence variation. Specifically, we assessed the degree and nature of polymorphism in a 307-base pair (bp) fragment of the mitochondrial cytochrome b gene, intron VIII of the gene for -enolase (a presumably neutral nuclear gene), and an ~600-bp fragment of an MHC class II B gene, including 155 bp of the hypervariable peptide binding region (a nuclear locus thought to be under balancing selection) for least and crested auklets (Aethia pusilla and A. cristatella; Charadriiformes: Alcidae). Transspecies polymorphism was found in both -enolase and the MHC but not cytochrome b and, given estimates of effective population size, probably represents retained ancestral variation. Biases in nucleotide composition suggested that mutational bias, tRNA availability, and the secondary structure of mRNA and/or DNA may influence base usage. Several lines of evidence indicated that balancing selection may be acting on the MHC II B exon 2. However, no evidence of balancing selection was observed in the intron and exon sequences immediately downstream of MHC II B exon 2. Current address (Hollie E. Walsh): Department of Zoology, University of Washington, Box 351800, Seattle, WA 98195-1800, USA     

Genetic Constraints on Protein Evolution

ABSTRACT

Evolution requires the generation and optimization of new traits (“adaptation”) and involves the selection of mutations that improve cellular function. These mutations were assumed to arise by selection of neutral mutations present at all times in the population. Here we review recent evidence that indicates that deleterious mutations are more frequent in the population than previously recognized and that these mutations play a significant role in protein evolution through continuous positive selection. Positively selected mutations include adaptive mutations, i.e. mutations that directly affect enzymatic function, and compensatory mutations, which suppress the pleiotropic effects of adaptive mutations. Compensatory mutations are by far the most frequent of the two and would allow potentially adaptive but deleterious mutations to persist long enough in the population to be positively selected during episodes of adaptation. Compensatory mutations are, by definition, context-dependent and thus constrain the paths available for evolution. This provides a mechanistic basis for the examples of highly constrained evolutionary landscapes and parallel evolution reported in natural and experimental populations. The present review article describes these recent advances in the field of protein evolution and discusses their implications for understanding the genetic basis of disease and for protein engineering in vitro.     

The Underlying Structure of Adaptation under Strong Selection in 12 Experimental Yeast Populations

The aims of this study were to determine (i) whether adaptation under strong selection occurred through mutations in a narrow target of one or a few nucleotide sites or a broad target of numerous sites and (ii) whether the programs of adaptation previously observed from three experimental populations were unique or shared among populations that underwent parallel evolution. We used archived population samples from a previous study, representing 500 generations of experimental evolution in 12 populations under strong selection, 6 populations in a high-salt environment and 6 populations in a low-glucose environment. Each set of six populations included four with sexual reproduction and two with exclusively asexual reproduction. Populations were sampled as resequenced genomes of 115 individuals and as bulk samples from which frequencies of mutant alleles were estimated. In a high-salt environment, a broad target of 11 mutations within the proton exporter, PMA1, was observed among the six populations, in addition to expansions of the ENA gene cluster. This pattern was shared among populations that underwent parallel evolution. In a low-glucose environment, two programs of adaptation were observed. The originally observed pattern of mutation in MDS3/MKT1 in population M8 was a narrow target of a single nucleotide, unique to this population. Among the other five populations, the three mutations were shared in a broad target, sensing/signaling genes RAS1 and RAS2. RAS1/RAS2 mutations were not observed in the high-salt populations; PMA1 mutations were observed only in a high-salt environment.     

Is the fixation of observable mutations distributed randomly among the three nucleotide positions of the codon?

Summary The distribution among the three nucleotide positions of the codons of 642 mutations fixed during the descent of 49 sequences of cytochromec was examined. This was compared to the distribution expected if the number of ways of getting a selectively acceptable amino acid alternative from a single nucleotide replacement at each coding position were random,i.e. proportional to the total number of ways of changing the encoded amino acid by a single nucleotide replacement at each coding position. It was found that the observed distribution was significantly different from random, there being 40% more mutations in the first coding position than in the second whereas one would have expected 10% more in the second than in the first. The probability of the result occurring by chance is < 10^–6.The same test was made on the distribution of 347 mutations fixed in the descent of 19 sequences of alpha hemoglobin and 286 mutations fixed in the descent of 16 beta and 4 delta hemoglobins. The result for the alpha hemoglobins was a similar non-randomness but the probability of its occurring by chance rose to 0.005. The result for the beta-delta hemoglobins was in the same direction but was not significant (p = 0.3). The degree of non-randomness among the three genes in the distribution of fixations over the three nucleotide positions of their codons appears to be correlated (negatively) with their rates of evolution, the plasticity required of the molecule to adapt to new environments, and the recency of exploitation of opportunities for change in functional specificity provided by such processes as gene duplication.     

The genetic basis of adaptation: lessons from concealing coloration in pocket mice

Recent studies on the genetics of adaptive coat-color variation in pocket mice (Chaetodipus intermedius) are reviewed in the context of several on-going debates about the genetics of adaptation. Association mapping with candidate genes was used to identify mutations responsible for melanism in four different populations of C. intermedius. Here, I review four main results (i) a single gene, the melanocortin-1-receptor (Mc1r), appears to be responsible for most of the phenotypic variation in color in one population, the Pinacate site; (ii) four or fewer nucleotide changes at Mc1r appear to be responsible for the difference in receptor function; (iii) studies of migration-selection balance suggest that the selection coefficient associated with the dark Mc1r allele at the Pinacate site is large; and (iv) different (unknown) genes underlie the evolution of melanism on three other lava flows. These findings are discussed in light of the evolution of convergent phenotypes, the average size of phenotypic effects underlying adaptation, the evolution of dominance, and the distinction between adaptations caused by changes in gene dosage versus gene structure.     

球孢白僵菌种内比较线粒体基因组学分析

【目的】明确球孢白僵菌种内线粒体基因组的分化程度。【方法】从GenBank下载已知的球孢白僵菌6个菌株线粒体基因组序列,详细分析基因组的组成结构,比较外显子区、内含子区和基因间区的碱基变异情况,分析菌株间的系统发育关系。【结果】球孢白僵菌不同菌株的线粒体基因组大小为28.8–32.3 kb,都有14个常见的核心蛋白编码基因、2个rRNA基因和25个tRNA基因,具有很强的共线性关系。但是,不同菌株含有的线粒体内含子数目存在差异(2–5个/菌株),在cox1、cox2和nad1基因中表现出内含子插入/缺失多态性,这是导致线粒体基因组大小变化的主要因素。对外显子、内含子和基因间区的碱基变异情况进行分析,发现内含子和基因间区相对变异较大,而外显子区相对变异较小。系统发育分析发现,这些球孢白僵菌菌株以很高的支持度聚在一起,具有相同内含子分布规律的菌株也具有较近的聚类关系。【结论】本研究首次报道球孢白僵菌因内含子数目不同、插入缺失突变和单核苷酸变异等在线粒体基因组上表现出较大程度的遗传分化,为认识真菌种内线粒体基因组分化提供了新的证据。     

Molecular population genetics of the male and female mitochondrial DNA molecules of the California sea mussel, Mytilus californianus

The presence of two gender-associated mitochondrial genomes in marine mussels provides a unique opportunity to investigate the dynamics of mtDNA evolution without complications inherent in interspecific comparisons. Here, we assess the relative importance of selection, mutation, and differential constraint in shaping the patterns of polymorphism within and divergence between the male (M) and female (F) mitochondrial genomes of the California sea mussel, Mytilus californianus. Partial sequences were obtained from homologous regions of four genes (nad2, cox1, atp6, and nad5) totaling 2307 bp in length. The M and F mtDNA molecules of M. californianus exhibited extensive levels of nucleotide polymorphism and were more highly diverged than observed in other mytilids (overall Tamura-Nei distances >40%). Consistent with previous studies, the M molecule had significantly higher levels of silent and replacement polymorphism relative to F. Both genomes possessed large numbers of singleton and low-frequency mutations that gave rise to significantly negative Tajima's D values. Mutation-rate scalars estimated for silent and replacement mutations were elevated in the M genome but were not sufficient to account for its higher level of polymorphism. McDonald-Kreitman tests were highly significant at all loci due to excess numbers of fixed replacement mutations between molecules. Strong purifying selection was evident in both genomes in keeping the majority of replacement mutations at low population frequencies but appeared to be slightly relaxed in M. Our results suggest that a reduction in selective constraint acting on the M genome remains the best explanation for its greater levels of polymorphism and faster rate of evolution.     

Molecular characterization of a novel,nuclear-encoded,NAD+-dependent glyceraldehyde-3-phosphate dehydrogenase in plastids of the gymnosperm Pinus sylvestris L.

Angiosperms and algae possess two distinct glyceraldehyde-3-phosphate dehydrogenase (GAPDH) enzymes, an NAD⁺-dependent tetramer involved in cytosolic glycolysis and an NADP⁺-dependent enzyme of the Calvin cycle in chloroplasts. We have found that the gymnosperm Pinus sylvestris possesses, in addition to these, a nuclear-encoded, plastid-specific, NAD⁺-dependent GAPDH, designated GapCp, which has not previously been described from any plant. Several independent full-size cDNAs for this enzyme were isolated which encode a functional transit peptide and mature subunit very similar to that of cytosolic GAPDH of angiosperms and algae. A molecular phylogeny reveals that chloroplast GapCp and cytosolic GapC arose through gene duplication early in chlorophyte evolution. The GapCp gene is expressed as highly as that for GapC in light-grown pine seedlings. These findings suggest that aspects of compartmentalized sugar phosphate metabolism may differ in angiosperms and gymnosperms and furthermore underscore the contributions of endosymbiotic gene transfer and gene duplication to the nuclear complement of genes for enzymes of plant primary metabolism.     

The comparative amino acid sequences, substrate specificities and gene or cDNA nucleotide sequences of some prokaryote and eukaryote amidinotransferases: implications for evolution

The amino acid sequences of the amidinotransferases and the nucleotide sequences of their genes or cDNA from four Streptomyces species (seven genes) and from the kidneys of rat, pig, human and human pancreas were compared. The overall amino acid and nucleotide sequences of the prokaryotes and eukaryotes were very similar and further, three regions were identified that were highly identical. Evidence is presented that there is virtually zero chance that the overall and high identity regions of the amino acid sequence similarities and the overall nucleotide sequence similarities between Streptomyces and mammals represent random match. Both rat and lamprey amidinotransferases were able to use inosamine phosphate, the amidine group acceptor of Streptomyces. We have concluded that the structure and function of the amidinotransferases and their genes has been highly conserved through evolution from prokaryotes to eukaryotes. The evolution has occurred with: (1) a high degree of retention of nucleotide and amino acid sequences; (2) a high degree of retention of the primitive Streptomyces guanine+cytosine (G+C) third codon position composition in certain high identity regions of the eukaryote cDNA; (3) a decrease in the specificities for the amidine group acceptors; and (4) most of the mutations silent in the regions suggested to code for active sites in the enzymes.     

Molecular evolution of paclitaxel biosynthetic genes <Emphasis Type="Italic">TS</Emphasis> and <Emphasis Type="Italic">DBAT</Emphasis> of <Emphasis Type="Italic">Taxus</Emphasis> species

Evolutionary patterns of sequence divergence were analyzed in genes from the conifer genus Taxus (yew), encoding paclitaxel biosynthetic enzymes taxadiene synthase (TS) and 10-deacetylbaccatin III-10β-O-acetyltransferase (DBAT). N-terminal fragments of TS, full-length DBAT and internal transcribed spacer (ITS) were amplified from 15 closely related Taxus species and sequenced. Premature stop codons were not found in TS and DBAT sequences. Codon usage bias was not found, suggesting that synonymous mutations are selectively neutral. TS and DBAT gene trees are not consistent with the ITS tree, where species formed monophyletic clades. In fact, for both genes, alleles were sometimes shared across species and parallel amino acid substitutions were identified. While both TS and DBAT are, overall, under purifying selection, we identified a number of amino acids of TS under positive selection based on inference using maximum likelihood models. Positively selected amino acids in the N-terminal region of TS suggest that this region might be more important for enzyme function than previously thought. Moreover, we identify lineages with significantly elevated rates of amino acid substitution using a genetic algorithm. These findings demonstrate that the pattern of adaptive paclitaxel biosynthetic enzyme evolution can be documented between closely related Taxus species, where species-specific taxane metabolism has evolved recently.