BACKGROUND
chrysanthemyl diphosphate synthase  catalyzes the condensation of two molecules of dimethylallyl diphosphate  to form chrysanthemyl diphosphate and is a key enzyme in biosynthetic pathways involving the production of pyrethrins and irregular monoterpenes  <cit> .
irregular monoterpenes are much less common than other isoprenoids and are confined to plants of the tribe anthemideae in the family asteraceae  <cit> .
these secondary metabolites play an important role in their defense against herbivorous insects  <cit> .
many enzymes involved in secondary plant metabolism are encoded by gene families that originated through gene duplications  <cit> .
evidence shows that the cds genes resulted from duplication of the farnesyl diphosphate synthase  genes that belong to a small family  <cit> , with copy numbers ranging from one in grape and two in arabidopsis thaliana, to five in rice  <cit> .
in artemisia tridentata , three fds genes have been identified: fds <dig>  fds <dig> and cds   <cit> .
despite the high sequence similarity of the three genes, divergent functions have been found among them  <cit> .
fds <dig> and fds <dig> catalyze the sequential head-to-tail condensation of two molecules of isopentenyl diphosphate  with dmapp to produce farnesyl diphosphate , and are involved in the biosynthesis of regular sesquiterpenes  <cit> ; whereas cds catalyzes two molecules of dmapp to form irregular monoterpenes  <cit> .
moreover, fds and its products are found in organisms ranging from prokaryotes to eukaryotes  <cit> , while the products of cds are only present in anthemideae  <cit> .thus, fds seems to be an ancient gene of which cds is a derived or modified orthologous copy, and irregular monoterpenes might be products of a pathway arising from those of other isoprenoids.
gene duplication is prevalent in plant genomes and the duplicated genes face different evolutionary fates, including pseudogenization , retention of the original function, subfunctionalization, and neofunctionalization under the functional view  <cit> .
the duplication origin and distinct function of cds raise a number of interesting questions.
first, when did the duplication event leading to the origin of cds take place?
because the products of cds occur exclusively in anthemideae  <cit> , it would be expected that the cds gene occurred at the same time as the origin of this tribe.
to date, cds sequences have only been cloned from two species of anthemideae   <cit>  and no information is available as to whether cds occurs in any other members of the anthemideae and the related relatives.
second, given the fact that cds uses a new nucleophilic substrate to generate new products, did cds accumulate mutations in sites homologous to the substrate-binding and catalysis sites of fds?
fds contains five conserved regions, two of which are ddxxd motifs  <cit> .
at the three-dimensional structure level, conserved amino-acids in the five regions and the c-terminal of fds are located in a large central cavity surrounded by  <dig> α-helices, which have been identified as substrate-binding and catalysis sites  <cit> .
previous studies demonstrated that cds has a t → g substitution in region iv and a d → n substitution in the first aspartate in region v  <cit> .
finally, what was the evolutionary force leading to the functional divergence of cds?
it is controversial whether the functional divergence of duplicate genes arises from the relaxation of selective constraints or positive selection  <cit> .
generally, neutral evolution with relaxed selective constraints is treated as the null hypothesis  <cit>  and positive selection is invoked if the null hypothesis is rejected.
positive selection has been detected in many genes after their duplication  <cit> .
it has been proposed that positive selection promotes the functional divergence of gene family members encoding enzymes involved in secondary metabolism  <cit> .
because cds carries out the production of irregular monoterpenes that are important secondary metabolites for defense against herbivorous insects in anthemideae, we investigated whether the functional divergence of cds was driven by positive selection.
in this study, we reconstructed the phylogeny of the fds gene family based on cdna and est sequences from the main asteraceae lineages.
we detected two rounds of gene duplication during the evolution of the asteraceae fds gene family, and inferred the possible time of origin of the cds gene.
homology modeling and molecular evolutionary analyses showed that two mutations in cds might be responsible for the fact that cds does not prefer ipp as the nucleophilic substrate like fds, and demonstrated that positive selection has played a role in the functional divergence of cds in anthemideae.
methods
amplification of fds homologs from anthemideae and its relatives
previous studies have well resolved the major clades of asteraceae and their relationships, with the anthemideae and astereae tribes being most closely related  <cit>  .
in this study, four species representing four subtribes of anthemideae  were sampled.
we also sampled one representative species from each of four tribes: aster ageratoides , helianthus annuus , taraxacum mongolicum , and gerbera anandria .
we first amplified fds homologs from genomic dna and then obtained the full-length cdna using gene-specific primers based on the partial dna fragments.
genomic dna was extracted using a plant genomic dna kit .
two cycles of polymerase chain reaction  were conducted to amplify sequences corresponding to conserved regions ii through v of fds and cds.
in the first round of pcr, the primers cdsii  and cdsvb  were used to amplify cds, while the primers fdsii  and fdsvb  were used to amplify fds.
in the second round of pcr, the primers cdsii and cdsva  were used to amplify cds, and fdsii and fdsva  were used to amplify fds.
each amplification reaction  contained 1× buffer,  <dig>  μm of each primer,  <dig> μm of each dntp, and  <dig>  u of taq polymerase, to which 1– <dig>  μl of each genomic dna template was added.
the thermocycling program comprised an initial  <dig> min at 95°c, followed by  <dig> cycles of  <dig> min at 95°c,  <dig> min at 52– <dig> °c depending on the dna template, and  <dig> min at 72°c, with a final extension step of  <dig> min at 72°c.
the amplification products were gel-purified and cloned into pgem-t vectors .
twenty positive clones were screened using restriction enzyme fragment analysis.
all distinct clones with the correct insertion were sequenced and the products were run on an abi <dig> automatic sequencer.
total rna from a mixture of leaves and shoots was extracted using a trizol kit , and 5′ rapid amplification of cdna ends  was performed using the 5′ race system .
after performing first-strand cdna synthesis with gene-specific primers, the original mrna was removed with rnase h and rnase t <dig>  and a polyc tail was added to the 5′-end.
then, two rounds of pcr amplification were performed with nested primers .
for 3′ race, first-strand cdna was produced using superscripttm ii reverse transcriptase .
then, one or two rounds of pcr amplification were carried out with nested primers .
the amplification products  were purified, cloned into pgem-t vectors  and sequenced.
the sequences were deposited in genbank .
est database search
a survey of genbank revealed that more than  <dig> million expressed sequence tag  sequences were available for species from five tribes of asteraceae  <cit> .
most of these sequences were contributed by the compositae genome project.
in the present study, libraries with more than  <dig>  ests were used.
for genera such as helianthus, although est libraries for several species are available, we used only one species to represent each genus.
hence, est libraries for nine asteraceae species representing five tribes were downloaded from genbank.
these species were helianthus exilis , cichorium intybus, lactuca saligna and taraxacum officinale , centaurea solstitialis, carthamus tinctorius and cynara scolymus , gerbera hybrida , and barnadesia spinosa .
a blast database for each est library was constructed using the formatdb program implemented in stand-alone blast  <dig> .
<dig> software  <cit> .
blast searches for sequences similar to fds in each database were conducted using the tblastn program with the fds <dig> protein sequence of a. tridentata <cit>  as the query.
overlapping ests were assembled manually.
detailed information for contigs and singletons  is listed in additional file  <dig>
phylogenetic analysis
phylogenetic analyses were based on the coding sequences.
the fds2-like sequence of lactuca saligna was excluded because of a large stretch of missing bases  at the n-terminal.
sequences were aligned based on the translated amino-acid sequences using clustalw in dambe  <cit> .
five additional fds sequences from species of asterids were retrieved from ncbi nonredundant sequence databases, with one sequence from a gentianaceae species as the outgroup .
maximum likelihood  analysis implemented in phyml version  <dig>   <cit>  and bayesian inference  implemented in mrbayes version  <dig> .
<dig>  <cit>  were used to construct phylogenetic trees.
the best-fit evolutionary model, gtr + i + g, was selected with the akaike information criterion using modeltest  <dig>   <cit> .
for the ml analysis, the starting tree was obtained with bionj, and parameter values were estimated from the data.
branch support was estimated from  <dig> bootstrap replicates .
in the bi analysis, two independent markov chain monte carlo runs were run simultaneously starting from a random tree for  <dig> million generations, sampled every  <dig> generations.
the first 10% of samples were discarded as burn-in, and the remaining trees were used to construct the 50% majority-rule consensus tree.
selection test
the codeml program in the paml 4b package  <cit>  was used to analyze possible positive selection acting on the fds gene family.
to reduce the impact of missing sites, our analyses were limited to fds genes that contained the full-length coding region.
first, branch models allowing the ω ratio  to vary among lineages  <cit>  were used to determine whether the selective pressure differed among different lineages.
the one ratio model  assumes the same ω for all branches and all sites.
the free ratio model  assumes an independent ω parameter for each branch in the tree.
in the phylogenetic analyses, three major clades, fds <dig>  fds <dig>  and cds, were resolved.
we assigned ω <dig>  ω <dig>  and ωc to the lineages ancestral to the fds <dig>  fds <dig>  and cds clades, respectively.
the two ratio models  assumed one ω ratio for branches of interest and the other ratio, ω <dig>  for all other branches; e.g., m2c assumed ωc for the branch ancestral to cds and ω <dig> for all other branches .
the three ratio models  assumed two branches of interest with different ω ratios and all other branches had a ratio of ω <dig>  a more complex four ratio model  assumed four independent ω ratios: one ratio each for the ancestral branches of fds <dig> , fds <dig> , and cds , and one for all other branches .
these models were compared using likelihood ratio tests  of the log likelihood  to check which model fit the data significantly better.
because the branch models average the ω ratio over all sites and were unable to detect a positive signal in many cases, site-specific models allowing ω to vary among sites  <cit>  were subsequently used to determine whether particular amino-acid residues within fds gene families have been subject to positive selection.
in addition to the one ratio model , five site models   <cit>  were used.
the nearly neutral model  assumes two classes of sites: conserved sites under strict constraint  and others under neutral selection .
the positive selection model  is an extension of m <dig> and assumes a third class of positively-selected sites .
the discrete model  uses a general discrete distribution with three site classes.
the beta model  assumes a beta distribution for the ω ratios over sites, while the beta&ω model  adds another site class to m <dig>  allowing the ω values to exceed  <dig>  three lrts of nested models were applied: m <dig> versus m <dig>  m <dig> versus m <dig>  and m <dig> versus m <dig>
as the branch model showed that the ω value for the branch ancestral to cds was significantly different from that for the other branches, we further used branch-site model a to test for sites that were potentially under positive selection on the branch ancestral to the cds subfamily  <cit> .
the model assumes four classes of sites.
the first two have ω <dig>  and ω <dig>  along all lineages in the phylogeny, whereas the third and fourth have ω <dig> along the ancestral cds branch, but ω <dig> and ω <dig> along other background branches.
the branch-site model a was compared with the nearly neutral model .
structure modeling
a homology model was constructed for cds  based on the crystal structure of human fds in complex with zoledronate and isopentenyl diphosphate .
the first  <dig> residues in the n-terminal of cds were cut off because this portion is removed in the mature cds protein  <cit> .
the align 2d structure alignment program  was used to align the sequences, and the modeler module of insightii was used to automatically generate models  <cit> .
to select the best model, all optimized models were evaluated using the profile-3d program  <cit> .
molecular graphics were created with pymol  <cit> .
RESULTS
characterization of the fds gene family
in total, we cloned  <dig> full-length cdna sequences for the fds gene family from  <dig> species of asteraceae.
for cds in leucanthemum vulgare and fds <dig> in aster ageratoides, we were unable to isolate a full-length cdna despite great attempts using different cdna templates, primers, pcr programs and annealing temperatures.
alignment of orf sequences revealed a 2-base insertion in helianthus annuus cds, which was confirmed by repeated pcr amplification, cloning, and sequencing.
this frameshift mutation leads to a premature stop codon, indicative of a nonfunctional pseudogene.
the length of the cfds sequences ranged from  <dig> to  <dig> bp, with the orf ranging from  <dig> to  <dig> bp, encoding proteins of  <dig> to  <dig> amino-acids.
with the exception of the a. ageratoides ccds, which included a large indel in the n-terminal, the ccds sequences varied in length from  <dig> to  <dig> bp, with the orf having a length of  <dig> to  <dig> bp, encoding proteins of  <dig> to  <dig> amino-acids.
compared to fds, cds exhibited an ~50-amino acid extension at its n-terminal.
the extension sequences of cds from different species were highly variable, being rich in serine and threonine residues  and showing a lack of acidic amino-acids.
they were identified as potential chloroplast transit peptides by targetp version  <dig>   <cit> .
these peptides shared little similarity with any other database sequence entries based on blast searches.
sixteen fds contigs and singletons were obtained from nine est libraries, of which thirteen had a length of more than  <dig> bp.
however, there were only three sequences with a full-length coding sequence: the fds <dig> genes from cichorium intybus, taraxacum officinale and cynara scolymus.
phylogeny of the fds gene family
two phylogenetic methods  generated almost the same tree topology and thus only ml tree is shown, with the internal node supports from two methods .
the ml tree clearly showed that all fds sequences from asteraceae form a well-supported clade.
fds sequences from the species of other families fell out of the asteraceae clade.
in this clade, there were three clusters of genes, one consisting of fds <dig> homologs, one fds <dig> homologs and another cds homologs, indicative of two duplication events happened in the evolution of fds gene family of asteraceae.
cds clade was sister to that of fds <dig>  suggesting that the second duplication event led to the origin of cds.
all the anthemideae species  had fds <dig>  fds <dig> and cds genes, whereas only two types of fds homologs were obtained from most species of other tribes.
in aster ageratoides of the tribe astereae and helianthus species of heliantheae, both fds <dig> and cds homologs were found.
in species of cichorieae and cardueae, one type of fds copy formed a cluster with the fds <dig> clade and another with the fds <dig> clade.
interestingly, two fds homologs were cloned from gerbera anandria of mutisieae, the basal tribe of asteraceae.
one of them fell in fds <dig> clade, while another clustered with the fds2+cds ancestral clade.
these findings indicated that the second duplication event giving rise to fds <dig> and cds occurred after the divergence of g. anandria and before the origin of the tribe anthemideae.
barnadesia spinosa, a species from the most basal tribe barnadesieae , had one fds homolog clustered with fds <dig>  suggesting that the first duplication event occurred in the common ancestor of asteraceae.
structure modeling
previous studies demonstrated that the following sixteen fds residues are involved in substrate-binding and catalysis: 56g, 57k, 60r, 96q, 103d, 107d, 112r, 113r, 174d, 200k, 201t, 239f, 240q, 243d, 257k, and 351r    <cit> .
among these sites of fds, four were mutated in homologous sites of cds: t <dig> in fds → g <dig>  in cds; f <dig> in fds → y <dig> in cds; d <dig> in fds → n <dig> in cds; and r <dig> in fds → g <dig> in cds .
of these mutated sites, two  involved in ipp substrate-binding in fds are shown in figure 2b.
the structural feature of fds is the arrangement of  <dig> core helices around a large central cavity, and the highly-conserved amino-acids are all located in this cavity  <cit> .
these and other conserved fds features are preserved in cds, as shown in figure 2c.
selection test
we compared the log likelihood values from different branch models to explore whether the ω ratios varied among different lineages and, particularly, whether the ratios for each ancestral branch to the fds <dig>  fds <dig>  and cds subfamilies differed from those for other branches in the phylogeny.
the results are shown in table  <dig>  the free ratio model  fit the data significantly better than the one ratio model  , suggesting that the ω ratios varied among lineages .
however, the ω < <dig> under mf indicated purifying selection in the gene family.
the two ratio models m2a and m2b, which assigned one ω to the ancestral lineages of fds <dig> or fds <dig> and the other ratio ω <dig> to all other branches, produced log likelihood values very similar to the one ratio model and were not significantly better than the one ratio model .
in contrast, the two ratio model m2c, with ωc =  <dig>  for the lineage ancestral to cds and ω <dig> =  <dig>  for all other branches, was significantly better than the one ratio model , indicating a significant increase in the ω ratio for the ancestral branch for cds.
finally, the three ratio models, m3a, m3b, and m3c, and the four-ratio model m4a were rejected in favor of m2c .
* significant at p < <dig>  level, *** significant at p < <dig>  level.
p, the number of free parameters for the ω ratio.
the lrts for m <dig> vs. m <dig>  suggested that the positive selection model  was not significantly better than the nearly neutral model .
although models m <dig> and m <dig> fit the data significantly better than the null models m <dig> and m <dig> , they did not identify sites with an ω value significantly greater than  <dig>
given that the branch ancestral to cds exhibited an increased ω  and cds was endowed with a new function in the biosynthesis of terpenoids, branch-site model a was further used to test for evidence of positive selection on the branch ancestral to cds .
lrts showed that this model was significantly better than the nearly neutral model m <dig>  the parameter estimates under branch-site model a suggested that  <dig> % of codons along the cds branch had been under positive selection, with ω =  <dig> .
bayes empirical bayes  analyses showed that at the p >50% level, branch-site model a identified  <dig> sites as being potentially subjected to positive selection on the cds branch.
at the p >80% level, the following  <dig> sites were identified: 102m, 103v, 113q, 166e, 191q, 198h, 211i, 218t, 220c, 225q, 235l, 236n, 241q, 273i, 279m, 281y, 285n, 290t, 293d, 295d, 300t, 306e, 307c, 334i, 353k, 355a, 356y, 386c and 394h  .
at the p >95% level, three sites, 103v, 218t and 290t, were identified.
discussion
gene duplication is considered to be a major mechanism in the generation of evolutionary novelty and adaptation  <cit> .
in plants, gene duplication followed by functional divergence is particularly important for the diversification of biochemical metabolites  <cit> .
isoprenoids are a large and diverse class of metabolites  <cit>  derived from five-carbon isoprene units, which can be classified into regular and irregular forms depending on the bond between isoprene units or monoterpenes, sesquiterpenes and diterpenes according to the number of isoprene units  <cit> .
fdss are involved in the biosynthesis of sesquiterpenes, and are encoded by a small gene family.
it seems that this gene family has experienced lineage-specific gene expansions multiple times.
for example, the two copies of arabidopsis thaliana formed a species-specific clade .
fds copies from oryza sative and sorghum bicolor formed a poaceae-specific clade .
based on phylogenetic analyses , we clearly showed that two rounds of gene duplication occurred in the evolutionary history of the asteraceae fds gene family.
the first round of duplication appears to have occurred in the common ancestor of the asteraceae, since the genes from asteraceae formed a monophyletic group separated from the clusters for species from other families , and even from the species of two families closely related to asteraceae .
the fds gene duplications in asteraceae might contribute to the diversity of their sesquiterpenes because of the role of fdss in the biosynthesis of sesquiterpenes.
this is consistent with the large number of sesquiterpenes that have been extracted from asteraceae  <cit> .
the second duplication, which generated the lineage of cds, occurred after the divergence of the mutisieae from the other tribes of asteraceae and before the divergence of the tribe anthemideae.
the evidence includes 1) one fds copy in g. anandria  clustered with fds <dig>  while the other clustered with the ancestor of fds <dig> and cds; 2) all of the sampled anthemideae species had three fds copies ; and 3) cds was also cloned from aster ageratoides  and helianthus annuus  that are close relatives of the tribe anthemideae  <cit> .
after the origin of cds, it developed a new function, involved in the biosynthetic pathway of irregular monoterpenes.
the cds gene is common in the tribe anthemideae, which is consistent with the fact that its products are typically found in anthemideae species.
our results suggested that the duplication and divergence of fds genes has played a major role in determining the novelty of irregular monoterpenes in anthemideae.
after gene duplication, cds accumulated amino-acid changes toward the change of a substrate.
cdss have four substitutions in the substrate-binding and catalysis sites of fds: t <dig> → g <dig> , f <dig> → y <dig>  d <dig> → n <dig>  and r <dig> → g <dig>  substitutions can be divided into either radical or conservative, based on the biochemical properties of the amino-acids  <cit> .
for example, substitutions associated with a change of polarity group are defined as radical and those with the polarity group unaltered as conservative  <cit> .
in the present study, except for t <dig> → s <dig>  all the substitutions are radical, which is consistent with the fact that the evolution of new function requires alterations in the biochemical properties of the amino-acid sequence  <cit> .
f <dig> and r <dig> in fds are involved in binding of the nucleophilic substrate ipp  <cit> .
the radical replacements of these two sites in cds are in good agreement with the finding that cds does not prefer ipp as a nucleophilic reagent.
f <dig> in fds binds ipp through hydrophobic interactions  <cit> .
the corresponding residue is y <dig> in cds, which, owing to the polarity of the hydroxyl, may not interact with ipp through hydrophobic interactions.
r <dig> in fds interacts with the pyrophosphate moiety of ipp through water-mediated hydrogen bonds  <cit>  .
the radical replacement of r <dig> → g <dig> involves changes in the charge  and the molecular volume of the amino-acids , which could affect the ipp binding.
hence, these two substitutions might explain why cds does not prefer ipp as a substrate.
substitutions can change the function of duplicated genes, and may be due to either a relaxation of purifying selection or to the action of positive selection  <cit> .
branch-site model a provided evidence of positive selection acting at  <dig>  sites along the branch ancestral to cds .
interestingly, y <dig>  and n <dig>  noted above were found to be under positive selection by the branch-site model, suggesting the important role of positive selection in the functional evolution of cds.
the biochemical context of substitutions that were under positive selection is consistent with a scenario involving the adaptive evolution of cds.
these sites  were scattered throughout the primary sequences , whereas in the three-dimensional structures , they clustered in the large central cavity.
among these sites, two  were located in conserved region i, and nine  were located in conserved region v. they are all conserved in the fds gene family and important for the precise function of the protein  <cit> .
the mutations at these sites suggested that their importance for the enzymatic activity of fds was altered in cds.
a few sites that were detected to have experienced positive selection with high probability may be responsible for the novel function of cds.
further studies using site-directed mutagenesis are needed to determine whether these positively-selected sites, especially those with high posterior probability , confer an ability on cds to discriminate different substrate types.
it has been proposed that positive selection promotes the functional divergence of gene family members encoding enzymes involved in secondary metabolism because secondary products are considered to be a response to challenges imposed by the environment  <cit> .
for example, the methylthioalkylmalate synthase gene  controls an early step in the biosynthesis of glucosinolates, which play an important role in arabidopsis thaliana and other crucifers’ defense against herbivorous insects  <cit> .
benderoth et al.
<cit>  found that positive selection had driven the evolution of mam <dig> that originated from a lineage-specific duplication of mama in a. thaliana.
another example is the sabath gene family of methyltransferases, which encodes enzymes catalyzing the formation of a variety of secondary metabolites in plants such as those that contribute to floral scent and plant defense.
branch-site analysis suggested that positive selection for a single amino-acid change promoted the substrate discrimination of salicylic acid methyltransferase  <cit> .
here, we provide an additional example in which positive selection has promoted the functional divergence of duplicated genes in a secondary metabolic pathway.
the adaptive evolution of the cds gene at the molecule level is consistent with the adaptive roles of the products of cds , and plays an important role in plant survival.
many models have been proposed to explain the evolutionary fates of duplicated genes, including neofunctionalization, duplication-degeneration-complementation  and escape from adaptive conflict  models  <cit> .
compared with cds gene in which adaptive selection has been detected, fds <dig> gene, as a sister duplicated copy of cds, did not show any signature of positive selection.
so, it seems that the evolution of fds <dig> and cds are not consistent with the predictions of the eac model, where both duplicated copies would evolve under positive selection  <cit> .
particularly, cds has a ~ 50-amino acid extension at its n-terminal, which was identified as a plastidial transit peptide, in agreement with the category ii-c model   <cit> .
however, the peptide of cds shares little similarity with any other sequence database entries by blast searches.
further work including functional analysis and the exploration on the origin of the peptide of cds would provide insights into the evolutionary fate of the fds gene family in asteraceae.
CONCLUSIONS
based on phylogenetic analyses of fds sequences, we demonstrated that two duplication events occurred in the evolution of the asteraceae fds gene family.
the first occurred in the common ancestor of the asteraceae and the second after the divergence of the mutisioideae from the other tribes, but before the birth of the anthemideae tribe.
we found that cds accumulated four mutations in sites homologous to the substrate-binding and catalysis sites of fds: t <dig> → g <dig> in conserved region iv, d <dig> → n <dig> in the first aspartate in conserved region v, f <dig> → y <dig> in region v, and r <dig> → g <dig> in the c-terminal.
of the four replaced sites of fds, f <dig> and r <dig> are involved in the binding of the nucleophilic substrate isopentenyl diphosphate.
likelihood analyses of a branch-site model provided evidence of positive selection acting on  <dig> sites  and  <dig> sites  on the branch ancestral to cds.
positive selection associated with gene duplication has played a major role in the evolution of cds.
abbreviations
bi: bayesian inference; cds: chrysanthemyl diphosphate synthase; dmapp: dimethylallyl diphosphate; eac: escape from adaptive conflict; est: expressed sequence tags; fds: farnesyl diphosphate synthase; ipp: isopentenyl diphosphate; mam: methylthioalkylmalate synthase; ml: maximum likelihood; race: rapid amplification of cdna ends.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
gyr, ypg and pll designed the experiments.
pll and jnw conducted the experiments.
gyr, ypg, sg and pll analyzed and interpreted the data.
gyr, ypg, sg and pll wrote the manuscript.
all authors read and approved the final manuscript.
supplementary material
additional file 1
primers for race.
click here for file
additional file 2
sequences included in this study.
click here for file
additional file 3
maximum likelihood tree of fds gene family from different plants.
numbers next to branches are bootstrap percentages from maximum likelihood analysis, and posterior probabilities from bayesian analysis.
click here for file
acknowledgments
we thank prof. jing-chu luo and dr. hong qu of peking university for their helps in the structural modeling analyses.
we thank dr. iain c. bruce at school of medicine, zhejiang university for linguistic editing.
the anonymous referees are also acknowledged for their valuable comments on the manuscript.
this work is supported by the national natural science foundation of china .

 chrysanthemyl diphosphate synthase catalyzes the condensation of two molecules of dimethylallyl diphosphate to form chrysanthemyl diphosphate and is a key enzyme in biosynthetic pathways involving the production of pyrethrins and irregular monoterpenes. irregular monoterpenes are much less common than other isoprenoids and are confined to plants of the tribe anthemideae in the family asteraceae. these secondary metabolites play an important role in their defense against herbivorous insects. many enzymes involved in secondary plant metabolism are encoded by gene families that originated through gene duplications. evidence shows that the cds genes resulted from duplication of the farnesyl diphosphate synthase genes that belong to a small family, with copy numbers ranging from one in grape and two in arabidopsis thaliana, to five in rice. in artemisia tridentata, three fds genes have been identified: fds fds and cds. despite the high sequence similarity of the three genes, divergent functions have been found among them. fds and fds catalyze the sequential head-to-tail condensation of two molecules of isopentenyl diphosphate with dmapp to produce farnesyl diphosphate, and are involved in the biosynthesis of regular sesquiterpenes; whereas cds catalyzes two molecules of dmapp to form irregular monoterpenes. moreover, fds and its products are found in organisms ranging from prokaryotes to eukaryotes, while the products of cds are only present in anthemideae.thus, fds seems to be an ancient gene of which cds is a derived or modified orthologous copy, and irregular monoterpenes might be products of a pathway arising from those of other isoprenoids. gene duplication is prevalent in plant genomes and the duplicated genes face different evolutionary fates, including pseudogenization, retention of the original function, subfunctionalization, and neofunctionalization under the functional view. the duplication origin and distinct function of cds raise a number of interesting questions. first, when did the duplication event leading to the origin of cds take place? because the products of cds occur exclusively in anthemideae, it would be expected that the cds gene occurred at the same time as the origin of this tribe. to date, cds sequences have only been cloned from two species of anthemideae and no information is available as to whether cds occurs in any other members of the anthemideae and the related relatives. second, given the fact that cds uses a new nucleophilic substrate to generate new products, did cds accumulate mutations in sites homologous to the substrate-binding and catalysis sites of fds? fds contains five conserved regions, two of which are ddxxd motifs. at the three-dimensional structure level, conserved amino-acids in the five regions and the c-terminal of fds are located in a large central cavity surrounded by α-helices, which have been identified as substrate-binding and catalysis sites. previous studies demonstrated that cds has a t → g substitution in region iv and a d → n substitution in the first aspartate in region v. finally, what was the evolutionary force leading to the functional divergence of cds? it is controversial whether the functional divergence of duplicate genes arises from the relaxation of selective constraints or positive selection. generally, neutral evolution with relaxed selective constraints is treated as the null hypothesis and positive selection is invoked if the null hypothesis is rejected. positive selection has been detected in many genes after their duplication. it has been proposed that positive selection promotes the functional divergence of gene family members encoding enzymes involved in secondary metabolism. because cds carries out the production of irregular monoterpenes that are important secondary metabolites for defense against herbivorous insects in anthemideae, we investigated whether the functional divergence of cds was driven by positive selection. in this study, we reconstructed the phylogeny of the fds gene family based on cdna and est sequences from the main asteraceae lineages. we detected two rounds of gene duplication during the evolution of the asteraceae fds gene family, and inferred the possible time of origin of the cds gene. homology modeling and molecular evolutionary analyses showed that two mutations in cds might be responsible for the fact that cds does not prefer ipp as the nucleophilic substrate like fds, and demonstrated that positive selection has played a role in the functional divergence of cds in anthemideae. methods amplification of fds homologs from anthemideae and its relatives previous studies have well resolved the major clades of asteraceae and their relationships, with the anthemideae and astereae tribes being most closely related. in this study, four species representing four subtribes of anthemideae were sampled. we also sampled one representative species from each of four tribes: aster ageratoides, helianthus annuus, taraxacum mongolicum, and gerbera anandria. we first amplified fds homologs from genomic dna and then obtained the full-length cdna using gene-specific primers based on the partial dna fragments. genomic dna was extracted using a plant genomic dna kit. two cycles of polymerase chain reaction were conducted to amplify sequences corresponding to conserved regions ii through v of fds and cds. in the first round of pcr, the primers cdsii and cdsvb were used to amplify cds, while the primers fdsii and fdsvb were used to amplify fds. in the second round of pcr, the primers cdsii and cdsva were used to amplify cds, and fdsii and fdsva were used to amplify fds. each amplification reaction contained 1× buffer, μm of each primer, μm of each dntp, and u of taq polymerase, to which 1– μl of each genomic dna template was added. the thermocycling program comprised an initial min at 95°c, followed by cycles of min at 95°c, min at 52– °c depending on the dna template, and min at 72°c, with a final extension step of min at 72°c. the amplification products were gel-purified and cloned into pgem-t vectors. twenty positive clones were screened using restriction enzyme fragment analysis. all distinct clones with the correct insertion were sequenced and the products were run on an abi automatic sequencer. total rna from a mixture of leaves and shoots was extracted using a trizol kit, and 5′ rapid amplification of cdna ends was performed using the 5′ race system. after performing first-strand cdna synthesis with gene-specific primers, the original mrna was removed with rnase h and rnase t and a polyc tail was added to the 5′-end. then, two rounds of pcr amplification were performed with nested primers. for 3′ race, first-strand cdna was produced using superscripttm ii reverse transcriptase. then, one or two rounds of pcr amplification were carried out with nested primers. the amplification products were purified, cloned into pgem-t vectors and sequenced. the sequences were deposited in genbank. est database search a survey of genbank revealed that more than million expressed sequence tag sequences were available for species from five tribes of asteraceae. most of these sequences were contributed by the compositae genome project. in the present study, libraries with more than ests were used. for genera such as helianthus, although est libraries for several species are available, we used only one species to represent each genus. hence, est libraries for nine asteraceae species representing five tribes were downloaded from genbank. these species were helianthus exilis, cichorium intybus, lactuca saligna and taraxacum officinale, centaurea solstitialis, carthamus tinctorius and cynara scolymus, gerbera hybrida, and barnadesia spinosa. a blast database for each est library was constructed using the formatdb program implemented in stand-alone blast. software. blast searches for sequences similar to fds in each database were conducted using the tblastn program with the fds protein sequence of a. tridentata as the query. overlapping ests were assembled manually. detailed information for contigs and singletons is listed in additional file phylogenetic analysis phylogenetic analyses were based on the coding sequences. the fds2-like sequence of lactuca saligna was excluded because of a large stretch of missing bases at the n-terminal. sequences were aligned based on the translated amino-acid sequences using clustalw in dambe. five additional fds sequences from species of asterids were retrieved from ncbi nonredundant sequence databases, with one sequence from a gentianaceae species as the outgroup. maximum likelihood analysis implemented in phyml version and bayesian inference implemented in mrbayes version. were used to construct phylogenetic trees. the best-fit evolutionary model, gtr + i + g, was selected with the akaike information criterion using modeltest. for the ml analysis, the starting tree was obtained with bionj, and parameter values were estimated from the data. branch support was estimated from bootstrap replicates. in the bi analysis, two independent markov chain monte carlo runs were run simultaneously starting from a random tree for million generations, sampled every generations. the first 10% of samples were discarded as burn-in, and the remaining trees were used to construct the 50% majority-rule consensus tree. selection test the codeml program in the paml 4b package was used to analyze possible positive selection acting on the fds gene family. to reduce the impact of missing sites, our analyses were limited to fds genes that contained the full-length coding region. first, branch models allowing the ω ratio to vary among lineages were used to determine whether the selective pressure differed among different lineages. the one ratio model assumes the same ω for all branches and all sites. the free ratio model assumes an independent ω parameter for each branch in the tree. in the phylogenetic analyses, three major clades, fds fds and cds, were resolved. we assigned ω ω and ωc to the lineages ancestral to the fds fds and cds clades, respectively. the two ratio models assumed one ω ratio for branches of interest and the other ratio, ω for all other branches; e.g., m2c assumed ωc for the branch ancestral to cds and ω for all other branches. the three ratio models assumed two branches of interest with different ω ratios and all other branches had a ratio of ω a more complex four ratio model assumed four independent ω ratios: one ratio each for the ancestral branches of fds, fds, and cds, and one for all other branches. these models were compared using likelihood ratio tests of the log likelihood to check which model fit the data significantly better. because the branch models average the ω ratio over all sites and were unable to detect a positive signal in many cases, site-specific models allowing ω to vary among sites were subsequently used to determine whether particular amino-acid residues within fds gene families have been subject to positive selection. in addition to the one ratio model, five site models were used. the nearly neutral model assumes two classes of sites: conserved sites under strict constraint and others under neutral selection. the positive selection model is an extension of m and assumes a third class of positively-selected sites. the discrete model uses a general discrete distribution with three site classes. the beta model assumes a beta distribution for the ω ratios over sites, while the beta&ω model adds another site class to m allowing the ω values to exceed three lrts of nested models were applied: m versus m m versus m and m versus m as the branch model showed that the ω value for the branch ancestral to cds was significantly different from that for the other branches, we further used branch-site model a to test for sites that were potentially under positive selection on the branch ancestral to the cds subfamily. the model assumes four classes of sites. the first two have ω and ω along all lineages in the phylogeny, whereas the third and fourth have ω along the ancestral cds branch, but ω and ω along other background branches. the branch-site model a was compared with the nearly neutral model. structure modeling a homology model was constructed for cds based on the crystal structure of human fds in complex with zoledronate and isopentenyl diphosphate. the first residues in the n-terminal of cds were cut off because this portion is removed in the mature cds protein. the align 2d structure alignment program was used to align the sequences, and the modeler module of insightii was used to automatically generate models. to select the best model, all optimized models were evaluated using the profile-3d program. molecular graphics were created with pymol.  characterization of the fds gene family in total, we cloned full-length cdna sequences for the fds gene family from species of asteraceae. for cds in leucanthemum vulgare and fds in aster ageratoides, we were unable to isolate a full-length cdna despite great attempts using different cdna templates, primers, pcr programs and annealing temperatures. alignment of orf sequences revealed a 2-base insertion in helianthus annuus cds, which was confirmed by repeated pcr amplification, cloning, and sequencing. this frameshift mutation leads to a premature stop codon, indicative of a nonfunctional pseudogene. the length of the cfds sequences ranged from to bp, with the orf ranging from to bp, encoding proteins of to amino-acids. with the exception of the a. ageratoides ccds, which included a large indel in the n-terminal, the ccds sequences varied in length from to bp, with the orf having a length of to bp, encoding proteins of to amino-acids. compared to fds, cds exhibited an ~50-amino acid extension at its n-terminal. the extension sequences of cds from different species were highly variable, being rich in serine and threonine residues and showing a lack of acidic amino-acids. they were identified as potential chloroplast transit peptides by targetp version. these peptides shared little similarity with any other database sequence entries based on blast searches. sixteen fds contigs and singletons were obtained from nine est libraries, of which thirteen had a length of more than bp. however, there were only three sequences with a full-length coding sequence: the fds genes from cichorium intybus, taraxacum officinale and cynara scolymus. phylogeny of the fds gene family two phylogenetic methods generated almost the same tree topology and thus only ml tree is shown, with the internal node supports from two methods. the ml tree clearly showed that all fds sequences from asteraceae form a well-supported clade. fds sequences from the species of other families fell out of the asteraceae clade. in this clade, there were three clusters of genes, one consisting of fds homologs, one fds homologs and another cds homologs, indicative of two duplication events happened in the evolution of fds gene family of asteraceae. cds clade was sister to that of fds suggesting that the second duplication event led to the origin of cds. all the anthemideae species had fds fds and cds genes, whereas only two types of fds homologs were obtained from most species of other tribes. in aster ageratoides of the tribe astereae and helianthus species of heliantheae, both fds and cds homologs were found. in species of cichorieae and cardueae, one type of fds copy formed a cluster with the fds clade and another with the fds clade. interestingly, two fds homologs were cloned from gerbera anandria of mutisieae, the basal tribe of asteraceae. one of them fell in fds clade, while another clustered with the fds2+cds ancestral clade. these findings indicated that the second duplication event giving rise to fds and cds occurred after the divergence of g. anandria and before the origin of the tribe anthemideae. barnadesia spinosa, a species from the most basal tribe barnadesieae, had one fds homolog clustered with fds suggesting that the first duplication event occurred in the common ancestor of asteraceae. structure modeling previous studies demonstrated that the following sixteen fds residues are involved in substrate-binding and catalysis: 56g, 57k, 60r, 96q, 103d, 107d, 112r, 113r, 174d, 200k, 201t, 239f, 240q, 243d, 257k, and 351r. among these sites of fds, four were mutated in homologous sites of cds: t in fds → g in cds; f in fds → y in cds; d in fds → n in cds; and r in fds → g in cds. of these mutated sites, two involved in ipp substrate-binding in fds are shown in figure 2b. the structural feature of fds is the arrangement of core helices around a large central cavity, and the highly-conserved amino-acids are all located in this cavity. these and other conserved fds features are preserved in cds, as shown in figure 2c. selection test we compared the log likelihood values from different branch models to explore whether the ω ratios varied among different lineages and, particularly, whether the ratios for each ancestral branch to the fds fds and cds subfamilies differed from those for other branches in the phylogeny. the results are shown in table the free ratio model fit the data significantly better than the one ratio model, suggesting that the ω ratios varied among lineages. however, the ω < under mf indicated purifying selection in the gene family. the two ratio models m2a and m2b, which assigned one ω to the ancestral lineages of fds or fds and the other ratio ω to all other branches, produced log likelihood values very similar to the one ratio model and were not significantly better than the one ratio model. in contrast, the two ratio model m2c, with ωc = for the lineage ancestral to cds and ω = for all other branches, was significantly better than the one ratio model, indicating a significant increase in the ω ratio for the ancestral branch for cds. finally, the three ratio models, m3a, m3b, and m3c, and the four-ratio model m4a were rejected in favor of m2c. * significant at p < level, *** significant at p < level. p, the number of free parameters for the ω ratio. the lrts for m vs. m suggested that the positive selection model was not significantly better than the nearly neutral model. although models m and m fit the data significantly better than the null models m and m, they did not identify sites with an ω value significantly greater than given that the branch ancestral to cds exhibited an increased ω and cds was endowed with a new function in the biosynthesis of terpenoids, branch-site model a was further used to test for evidence of positive selection on the branch ancestral to cds. lrts showed that this model was significantly better than the nearly neutral model m the parameter estimates under branch-site model a suggested that % of codons along the cds branch had been under positive selection, with ω =. bayes empirical bayes analyses showed that at the p >50% level, branch-site model a identified sites as being potentially subjected to positive selection on the cds branch. at the p >80% level, the following sites were identified: 102m, 103v, 113q, 166e, 191q, 198h, 211i, 218t, 220c, 225q, 235l, 236n, 241q, 273i, 279m, 281y, 285n, 290t, 293d, 295d, 300t, 306e, 307c, 334i, 353k, 355a, 356y, 386c and 394h. at the p >95% level, three sites, 103v, 218t and 290t, were identified. discussion gene duplication is considered to be a major mechanism in the generation of evolutionary novelty and adaptation. in plants, gene duplication followed by functional divergence is particularly important for the diversification of biochemical metabolites. isoprenoids are a large and diverse class of metabolites derived from five-carbon isoprene units, which can be classified into regular and irregular forms depending on the bond between isoprene units or monoterpenes, sesquiterpenes and diterpenes according to the number of isoprene units. fdss are involved in the biosynthesis of sesquiterpenes, and are encoded by a small gene family. it seems that this gene family has experienced lineage-specific gene expansions multiple times. for example, the two copies of arabidopsis thaliana formed a species-specific clade. fds copies from oryza sative and sorghum bicolor formed a poaceae-specific clade. based on phylogenetic analyses, we clearly showed that two rounds of gene duplication occurred in the evolutionary history of the asteraceae fds gene family. the first round of duplication appears to have occurred in the common ancestor of the asteraceae, since the genes from asteraceae formed a monophyletic group separated from the clusters for species from other families, and even from the species of two families closely related to asteraceae. the fds gene duplications in asteraceae might contribute to the diversity of their sesquiterpenes because of the role of fdss in the biosynthesis of sesquiterpenes. this is consistent with the large number of sesquiterpenes that have been extracted from asteraceae. the second duplication, which generated the lineage of cds, occurred after the divergence of the mutisieae from the other tribes of asteraceae and before the divergence of the tribe anthemideae. the evidence includes 1) one fds copy in g. anandria clustered with fds while the other clustered with the ancestor of fds and cds; 2) all of the sampled anthemideae species had three fds copies; and 3) cds was also cloned from aster ageratoides and helianthus annuus that are close relatives of the tribe anthemideae. after the origin of cds, it developed a new function, involved in the biosynthetic pathway of irregular monoterpenes. the cds gene is common in the tribe anthemideae, which is consistent with the fact that its products are typically found in anthemideae species. our results suggested that the duplication and divergence of fds genes has played a major role in determining the novelty of irregular monoterpenes in anthemideae. after gene duplication, cds accumulated amino-acid changes toward the change of a substrate. cdss have four substitutions in the substrate-binding and catalysis sites of fds: t → g, f → y d → n and r → g substitutions can be divided into either radical or conservative, based on the biochemical properties of the amino-acids. for example, substitutions associated with a change of polarity group are defined as radical and those with the polarity group unaltered as conservative. in the present study, except for t → s all the substitutions are radical, which is consistent with the fact that the evolution of new function requires alterations in the biochemical properties of the amino-acid sequence. f and r in fds are involved in binding of the nucleophilic substrate ipp. the radical replacements of these two sites in cds are in good agreement with the finding that cds does not prefer ipp as a nucleophilic reagent. f in fds binds ipp through hydrophobic interactions. the corresponding residue is y in cds, which, owing to the polarity of the hydroxyl, may not interact with ipp through hydrophobic interactions. r in fds interacts with the pyrophosphate moiety of ipp through water-mediated hydrogen bonds. the radical replacement of r → g involves changes in the charge and the molecular volume of the amino-acids, which could affect the ipp binding. hence, these two substitutions might explain why cds does not prefer ipp as a substrate. substitutions can change the function of duplicated genes, and may be due to either a relaxation of purifying selection or to the action of positive selection. branch-site model a provided evidence of positive selection acting at sites along the branch ancestral to cds. interestingly, y and n noted above were found to be under positive selection by the branch-site model, suggesting the important role of positive selection in the functional evolution of cds. the biochemical context of substitutions that were under positive selection is consistent with a scenario involving the adaptive evolution of cds. these sites were scattered throughout the primary sequences, whereas in the three-dimensional structures, they clustered in the large central cavity. among these sites, two were located in conserved region i, and nine were located in conserved region v. they are all conserved in the fds gene family and important for the precise function of the protein. the mutations at these sites suggested that their importance for the enzymatic activity of fds was altered in cds. a few sites that were detected to have experienced positive selection with high probability may be responsible for the novel function of cds. further studies using site-directed mutagenesis are needed to determine whether these positively-selected sites, especially those with high posterior probability, confer an ability on cds to discriminate different substrate types. it has been proposed that positive selection promotes the functional divergence of gene family members encoding enzymes involved in secondary metabolism because secondary products are considered to be a response to challenges imposed by the environment. for example, the methylthioalkylmalate synthase gene controls an early step in the biosynthesis of glucosinolates, which play an important role in arabidopsis thaliana and other crucifers’ defense against herbivorous insects. benderoth et al. found that positive selection had driven the evolution of mam that originated from a lineage-specific duplication of mama in a. thaliana. another example is the sabath gene family of methyltransferases, which encodes enzymes catalyzing the formation of a variety of secondary metabolites in plants such as those that contribute to floral scent and plant defense. branch-site analysis suggested that positive selection for a single amino-acid change promoted the substrate discrimination of salicylic acid methyltransferase. here, we provide an additional example in which positive selection has promoted the functional divergence of duplicated genes in a secondary metabolic pathway. the adaptive evolution of the cds gene at the molecule level is consistent with the adaptive roles of the products of cds, and plays an important role in plant survival. many models have been proposed to explain the evolutionary fates of duplicated genes, including neofunctionalization, duplication-degeneration-complementation and escape from adaptive conflict models. compared with cds gene in which adaptive selection has been detected, fds gene, as a sister duplicated copy of cds, did not show any signature of positive selection. so, it seems that the evolution of fds and cds are not consistent with the predictions of the eac model, where both duplicated copies would evolve under positive selection. particularly, cds has a ~ 50-amino acid extension at its n-terminal, which was identified as a plastidial transit peptide, in agreement with the category ii-c model. however, the peptide of cds shares little similarity with any other sequence database entries by blast searches. further work including functional analysis and the exploration on the origin of the peptide of cds would provide insights into the evolutionary fate of the fds gene family in asteraceae.  based on phylogenetic analyses of fds sequences, we demonstrated that two duplication events occurred in the evolution of the asteraceae fds gene family. the first occurred in the common ancestor of the asteraceae and the second after the divergence of the mutisioideae from the other tribes, but before the birth of the anthemideae tribe. we found that cds accumulated four mutations in sites homologous to the substrate-binding and catalysis sites of fds: t → g in conserved region iv, d → n in the first aspartate in conserved region v, f → y in region v, and r → g in the c-terminal. of the four replaced sites of fds, f and r are involved in the binding of the nucleophilic substrate isopentenyl diphosphate. likelihood analyses of a branch-site model provided evidence of positive selection acting on sites and sites on the branch ancestral to cds. positive selection associated with gene duplication has played a major role in the evolution of cds. abbreviations bi: bayesian inference; cds: chrysanthemyl diphosphate synthase; dmapp: dimethylallyl diphosphate; eac: escape from adaptive conflict; est: expressed sequence tags; fds: farnesyl diphosphate synthase; ipp: isopentenyl diphosphate; mam: methylthioalkylmalate synthase; ml: maximum likelihood; race: rapid amplification of cdna ends. competing interests the authors declare that they have no competing interests. authors’ contributions gyr, ypg and pll designed the experiments. pll and jnw conducted the experiments. gyr, ypg, sg and pll analyzed and interpreted the data. gyr, ypg, sg and pll wrote the manuscript. all authors read and approved the final manuscript. supplementary material additional file 1 primers for race. click here for file additional file 2 sequences included in this study. click here for file additional file 3 maximum likelihood tree of fds gene family from different plants. numbers next to branches are bootstrap percentages from maximum likelihood analysis, and posterior probabilities from bayesian analysis. click here for file acknowledgments we thank prof. jing-chu luo and dr. hong qu of peking university for their helps in the structural modeling analyses. we thank dr. iain c. bruce at school of medicine, zhejiang university for linguistic editing. the anonymous referees are also acknowledged for their valuable comments on the manuscript. this work is supported by the national natural science foundation of china.
these compounds are confined to plants of the tribe anthemideae of the asteraceae, and play an important role in defending the plants against herbivorous insects. chrysanthemyl diphosphate synthase is a key enzyme in biosynthetic pathways producing pyrethrins and irregular monoterpenes. positive selection associated with gene duplication played a major role in the evolution of cds. two duplication events were detected in the evolutionary history of the fds gene family in the asteraceae, and the second duplication led to the origin of cds. maximum likelihood analyses showed that some sites in cds were under positive selection and were scattered throughout primary sequences, whereas in the three-dimensional structure model they clustered in the large central cavity. however, the duplication time toward the origin of cds and the evolutionary force behind the functional divergence of the cds gene are still unknown. after its origin, cds accumulated four mutations in sites homologous to the substrate-binding and catalysis sites of fds. cds occurred after the divergence of the tribe mutisieae from other tribes of asteraceae but before the birth of the anthemideae tribe. of these, two sites were involved in the binding of the nucleophilic substrate isopentenyl diphosphate in fds. it has been proposed that the cds genes arose from duplication of the farnesyl diphosphate synthase gene and have different function from fdss.
 chrysanthemyl diphosphate synthase is a key enzyme in biosynthetic pathways producing pyrethrins and irregular monoterpenes. these compounds are confined to plants of the tribe anthemideae of the asteraceae, and play an important role in defending the plants against herbivorous insects. it has been proposed that the cds genes arose from duplication of the farnesyl diphosphate synthase gene and have different function from fdss. however, the duplication time toward the origin of cds and the evolutionary force behind the functional divergence of the cds gene are still unknown.  two duplication events were detected in the evolutionary history of the fds gene family in the asteraceae, and the second duplication led to the origin of cds. cds occurred after the divergence of the tribe mutisieae from other tribes of asteraceae but before the birth of the anthemideae tribe. after its origin, cds accumulated four mutations in sites homologous to the substrate-binding and catalysis sites of fds. of these, two sites were involved in the binding of the nucleophilic substrate isopentenyl diphosphate in fds. maximum likelihood analyses showed that some sites in cds were under positive selection and were scattered throughout primary sequences, whereas in the three-dimensional structure model they clustered in the large central cavity.  positive selection associated with gene duplication played a major role in the evolution of cds. 