Evolution of land plant genes encoding L-Ala-D/L-Glu epimerases (AEEs) via horizontal gene transfer and positive selection
- Zefeng Yang†1,
- Yifan Wang†1,
- Yong Zhou1,
- Qingsong Gao1,
- Enying Zhang1,
- Lei Zhu1,
- Yunyun Hu1 and
- Chenwu Xu1Email author
© Yang et al.; licensee BioMed Central Ltd. 2013
Received: 23 January 2013
Accepted: 21 February 2013
Published: 1 March 2013
The L-Ala-D/L-Glu epimerases (AEEs), a subgroup of the enolase superfamily, catalyze the epimerization of L-Ala-D/L-Glu and other dipeptides in bacteria and contribute to the metabolism of the murein peptide of peptidoglycan. Although lacking in peptidoglycan, land plants possess AEE genes that show high similarity to those in bacteria.
Similarity searches revealed that the AEE gene is ubiquitous in land plants, from bryophytas to angiosperms. However, other eukaryotes, including green and red algae, do not contain genes encoding proteins with an L-Ala-D/L-Glu_epimerase domain. Homologs of land plant AEE genes were found to only be present in prokaryotes, especially in bacteria. Phylogenetic analysis revealed that the land plant AEE genes formed a monophyletic group with some bacterial homologs. In addition, land plant AEE proteins showed the highest similarity with these bacterial homologs and shared motifs only conserved in land plant and these bacterial AEEs. Integrated information on the taxonomic distribution, phylogenetic relationships and sequence similarity of the AEE proteins revealed that the land plant AEE genes were acquired from bacteria through an ancient horizontal gene transfer (HGT) event. Further evidence revealed that land plant AEE genes had undergone positive selection and formed the main characteristics of exon/intron structures through gaining some introns during the initially evolutionary period in the ancestor of land plants.
The results of this study clearly demonstrated that the ancestor of land plants acquired an AEE gene from bacteria via an ancient HGT event. Other findings illustrated that adaptive evolution through positive selection has contributed to the functional adaptation and fixation of this gene in land plants.
Horizontal gene transfer (HGT), also known as lateral gene transfer (LGT), refers to the transfer of genetic material between organisms that are reproductively isolated [1, 2]. HGT plays important roles in accelerating the evolution of the acceptor lineages because it can considerably expand the gene pool beyond species barriers. It is believed that HGT is one of the major forces driving the evolution of prokaryotes, leading to the acquisition or modification of certain adaptive traits, such as antibiotic resistance, virulence, and photosynthesis . For example, a genome-wide analysis revealed that 755 out of 4,288 genes have been transferred to the Escherichia coli genome, and at least 234 HGT events contributed to the origin of these transferred genes .
The frequency of HGT into eukaryotic genomes is possibly lower than in prokaryotes, but HGT has also been an important force in the evolution of eukaryotes . Genome-wide identification revealed that 7.6% of the secreted proteome of Phytophthora ramorum has been acquired from fungi via HGT, suggesting that oomycetes became successful plant parasites through multiple acquisitions of genes from fungi . Single-celled organisms have been found to be the dominant agent for genetic transfer, and many microbial eukaryotes and plant mitochondria provide rich in examples of HGT . Gene recruitment has been thought to be difficult in multi-cellular eukaryotes. However, recent investigations in plants, fungi and animals have refuted this conjecture [8, 9]. The classical example of plant HGT is that the ct-DNA sequences in some tobacco nuclear genomes were probably horizontally acquired from Agrobacterium rhizogenes during ancient infections . In addition, evidence has also showed that land plants can recruit genes from species with distinct relationships, such as fungi, bacteria, and other plant species [11–15]. Although land plant genes acquired through HGT are quite rare, they might play critical roles in adaptation to environments. For example, some anciently employed genes have been found to be involved in many plant-specific activities, including xylem formation, plant defense, nitrogen recycling and the biosynthesis of starch, polyamines, hormones and glutathione .
L-Ala-D/L-Glu epimerase (AEE) belongs to the enolase superfamily and catalyzes the epimerization of L-Ala-D/L-Glu and other dipeptides . Studies examining AEEs in E. coli and Bacillus subtilis indicated a probable role in the metabolism of the murein peptide of peptidoglycan, of which L-Ala-D-Glu is a component [16, 17]. However, it has been shown that the AEE family contains members from plants and archaea that lack peptidoglycan, suggesting that the proteins of this family might have other functions . A recent investigation in Thermotoga martima, a species of thermotoga bacteria, resulted in the assignment of epimerase activity for L-Ala-D/L-Phe, L-Ala-D/L-Tyr, and L-Ala-D/L-His to one member of the AEE family . It has also been noted that the genomes of Oryza and Arabidopsis possess genes encoding AEEs . However, the AEE member found in Arabidopsis has long been annotated as a protein of the cytochrome P450 superfamily and was looked at as a pseudogene until it was found to be expressed . Because the Arabidopsis AEE shows similarity with bacterial TfdD, an enzyme in the degradation pathway for chlorinated aromatics, it is assumed to have the potential to degrade aromatic compounds when the bacterial TfdC gene is introduced to the plant .
The ubiquity of the AEE gene in land plants suggests that its functions could include a wide range of selectivity, although its actual function remains unclear. The origin of this gene in land plants remains unknown at present. The increased availability of AEE sequences in public databases allows us to explore the functional diversity from a phylogenetic perspective within the AEE family in land plants. Here, we examined the evolutionary relationship of the land plant AEE genes and their homologs in cellular organisms. Our bioinformatic analyses revealed that land plant AEE genes originated from an ancient HGT event, and the putative donor was bacteria. Further evidence showed that positive selection followed by purifying selection has contributed to the evolution of this gene in land plants.
AEEgenes in land plants
List of AEE genes in 13 representative land plant genomes
6261905 - 6264309
4532099 - 4536153
32914724 - 32919062
7188312 - 7191927
17393161 - 17395652
3148656 - 3153072
148949 - 150726
42573 - 44511
6866003 - 6869858
8094542 - 8098228
2074753 - 2077903
19939526 - 19944315
13336541 - 13339982
1705605 - 1709251
1839471 - 1840852
1851300 - 1852645
98135 - 100510
The origin of land plant AEEgenes
There is no doubt that land plants originated from green algae and that most of the genes in the genomes of land plants were vertically inherited from their common ancestor [22, 23]. We searched the nr and EST databases of NCBI and available eukaryotic genome databases for homologs of land plant AEE genes. To our surprise, the results indicated that there was no homolog in any other eukaryote, including in the genomes of green and red algae. Blast results also revealed that homologs of land plant AEE proteins only existed in prokaryotes, mainly in bacteria. The taxonomic distribution of the AEE genes suggested that their emergence in land plants might have been a result of horizontal gene transfer (HGT) from a prokaryote, and the universality of their distribution in bacteria also suggested that this gene first emerged in bacteria.
HGT has been demonstrated to be one of the major forces driving the evolution of prokaryotes. Recently, accumulating data have indicated that this process has also occurred during the evolution of eukaryotic genomes . However, in eukaryotes, another contributor to the accumulation of nuclear genes is intracellular gene transfer (IGT), which is gene transfer from the genomes of mitochondria or plastids to the nucleus . It has been demonstrated that all mitochondrial genomes originated from an ancient endosymbiotic uptake of an alphaproteobacterium . In addition, the chloroplast genome has been confirmed to have originated from the genome of an ancestor of extant cyanobacteria . Thus, genes showing cyanobacterial and plastid-containing eukaryotic homologs as top hits were mostly considered plastid derived, while those with alphaproteobacterial and other eukaryotic homologs as top hits were considered to likely have been mitochondrion derived . In our analysis, additional searches were performed to exclude the possibility of an IGT origin of land plant AEE genes. First, we searched the nr database, NCBI dbEST database and available eukaryotic genomic databases and found that there was no eukaryotic gene encoding an AEE other than those in land plants. Second, database searches revealed that no AEE protein was encoded by a mitochondrial or chloroplast gene. Third, although there were AEE genes found in both alphaproteobacteria and cyanobacteria, none of the AEE genes in these bacteria fell within the land plant branch in the phylogeny. Therefore, the scenario that land plant AEE genes originated through IGT requires too many independent gene loss events to seem likely. In addition, the hypothesis that this gene was present in the common ancestor of eukaryotes and only retained in land plants also requires numerous independent gene losses. Thus, under the assumption that the chance of the same gene being repeatedly transferred among different organismal groups is relatively low, the most parsimonious explanation is that the origin of land plant AEE genes was the result of an ancient HGT event from bacteria.
Selective constraints on land plant AEEgenes
Likelihood ratio tests of positive selection were applied using ML methods and the codon substitution models of Yang and his colleagues [26–28]. First, we compared models M0 and M3 to evaluate whether there were variations in the d N /d S ratio among codon positions in the AEE genes of land plants (Additional file 2). Overall, the maximum likelihood estimate of the d N /d S value for land plant AEE genes under model M0 was 0.1401, suggesting that relaxed purifying selection was the predominant force for the evolution of the AEE genes in land plants. Interestingly, the log-likelihood differences between models M3 and M0 were statistically significant (LRT=644.6478, p <0.01), illustrating that the overall level of selective constraint has fluctuated. Second, the LRTs employed to compare the fit of the data to model M2a vs. M1a and M8 vs. M7 were used to address whether positive selection promoted the divergence of this family in land plants. To our surprise, neither of these comparisons provided evidence of positive selection. This result revealed that the main constraint on the evolution of AEE genes in land plants was relaxed purifying selection following fixation after an HGT event. To compare the driving forces with the AEE genes in bacteria, 7 bacterial AEE genes that fell within the same branch as those in land plants were also selected to test selective constraints. The results revealed that purifying selection was the predominant force for the evolution of the AEE genes in bacteria, and no positive selection signature was found during their evolution in bacteria.
Parameters of the branch-site models used for the detection of positive selection
p0 = 0.6577, p1 = 0.1957, p2a = 0.1130, p2b = 0.0336
Background: ω0 = 0.0838, ω1 = 1.0000, ω2a = 0.0838, ω2b = 1.0000
Foreground: ω0 = 0.0838, ω1 = 1.0000, ω2a = 1.0000, ω2b = 1.0000
p0 = 0.6537, p1 = 0.1911, p2a = 0.1201, p2b = 0.0351
Background: ω0 = 0.0845, ω1 = 1.0000, ω2a = 0.0845, ω2b = 1.0000
Foreground: ω0 = 0.0845, ω1 = 1.0000, ω2a = 999.0000, ω2b = 999.0000
The origin of land plants has played fundamental roles in the formation of modern terrestrial ecosystems . Numerous lines of evidence have revealed that land plants evolved from water-based green algae. In the transmission from water to land, the ancestor of land plants is expected to have evolved genes with new functions to colonize the land. Among the mechanisms underlying the formation of these genes, HGT is one essential way to acquire new genetic material . Recent investigations have revealed that HGT-derived genes play important roles in plant colonization of land, as some land plant genes that function in plant-specific activities, including plant defense, stress tolerance and the biosynthesis of plant polyamines and hormones, have been demonstrated to have been acquired through HGT . In the present work, using integrated information on the taxonomic distribution, phylogenetic relationships and sequence similarity of the proteins possessing L-Ala-D/L-Glu_epimerase domains, we concluded that an ancient HGT event from bacteria contributed to the origin of AEE genes in land plants. The function of these genes in bacteria was demonstrated to be metabolism of the murein peptide of peptidoglycan [16, 17]. Although no peptidoglycan has been found in land plants, the ubiquity of AEE genes in land plants and the evidence of their expression indicate that they are functional and may play important roles in the growth and development of plants. These genes are expected to exhibit other functions or have evolved new functions in land plants.
In this study, we also noted that all of the AEE genes include introns in their coding regions and that the positions and phases of these introns are quite conserved, illustrating that most of the introns were present in the ancestor of land plants. The vast majority of prokaryotic genes contain no introns, and the only introns that have been shown to be present in prokaryotic genes are self-splicing type II introns, which are functionally quite distinct from the spliceosome-dependent nuclear introns in eukaryotic genes . Because no sequences showed highly similarity with these introns, it is currently unclear where they originated. However, it is conceivable that the introns in land plant AEE genes arose through insertions shortly after the HGT event and before the separation of land plant lineages.
It has been demonstrated that the phenotypic diversity of a gene family is controlled by selection as a function of evolutionary fitness . A rigorous and clear signal of selection pressure in molecular evolution is a significantly higher nonsynonymous (d N ; resulting in amino acid replacement) than synonymous (d S ; silent) substitution rate. The ratio of the two rates, d N /d S , or ω, measures the quantity and direction of selective pressure on a protein, where ω≈1, ω<1, and ω>1, indicate neutral evolution, purifying selection, and positive selection, respectively . Purifying selection is important for the evolution of a gene family because it can help the genes that belong to a family maintain their optimal function. However, positive selection is an important source of evolutionary innovation and is a major force underlying the adaptation of species to a new environment . In our analysis, we found that the dominant driving force for AEE genes was purifying selection in both land plants and bacteria, which would contribute to functional stabilization. However, when we employed the bacterial genes as background, positive selection was found to contribute greatly to the evolution of land plant AEE genes.
In general, positive selection is thought to act on only a few amino acid sites and for a short evolutionary period . Land plant AEE genes originated from bacteria through HGT as well as both the genomic and living conditions differ tremendously between bacteria and land plants. The functional adaptation of the AEE genes to the genomic and living environment of the ancestor of land plants was aided by positive selection. A successful HGT event leading to gene fixation results from providing a benefit to the host. Through positive selection, the AEE gene underwent complete functional innovation during a short evolutionary period in the ancestor of land plants. It is thought that if a transferred protein is not functional, neutral mutation will occur in the gene encoding it. The fate of the transferred gene will therefore be that it will be lost during evolution because of the accumulation of mutations. In addition to facilitating the adaptation of an organism to a particular niche, HGT can also provide a mechanism for genomic innovation and plasticity. After acquiring materials for innovation and adaptation to a new environment, positive selection acting on the transferred gene will modify its sequences to generate new functions. Thus, positive selection will reduce the chances of transferred gene losses caused by the accumulation of mutations.
The gene encoding L-Ala-D/L-Glu epimerase (AEE) was found to be present in all of the available sequenced genomes of land plants, whereas homologs of this gene were not found in any other eukaryotic genome, including those of green and red algae. In this study, we performed extensive analyses of the taxonomic distribution and phylogeny of the AEE protein, which catalyzes the epimerization of L-Ala-D/L-Glu and other dipeptides and plays an important role in the metabolism of the murein peptide of peptidoglycan in bacteria. Our results revealed that the ancestor of land plants acquired the AEE gene from bacteria through an ancient HGT event. We also noted that rapid evolution and drastic sequence variation occurred during the initially short evolutionary period of the AEE gene in land plants following HGT. In addition to generating additional introns in the coding region of the gene, adaptive evolution via positive selection helped the AEE to undergo functional innovation and fixation in the genome of the land plant ancestor.
Sequence data sources
To identify the land plant genes encoding AEEs, BLASTP searches were performed in the Phytozome database  using the amino acid sequence of the B. subtilis YkfB gene  as a query. The CD-search tool in the Conserved Domain Database  of NCBI was used to predict the L-Ala-DL-Glu_epimerase domain (cd03319) for the obtained BLAST hits. The proteins that contained this conserved domain were defined as land plant AEEs. The new AEE sequences detected in land plants were used reiteratively to search the respective sequence database. EST searches for land plant AEE genes were performed using the BLASTN tool against the EST database of NCBI.
To identify the homologs of land plant AEE genes, BLAST searches against the non-redundant (nr) protein sequence database, NCBI EST database and available eukaryotic genome databases (Additional file 3) were performed using the land plant AEE protein sequences as queries. The obtained hits were further analyzed via an NCBI conserved domain search to confirm the presence of the L-Ala-D/L-Glu_epimerase domain in their protein structure. Protein sequences were sampled for further combined phylogenetic analysis from representative groups within each domain of life (bacteria, archaea, and eukaryotes) based on the BLASTP results.
Multiple sequence alignment and phylogenetic tree reconstruction
All of the selected representative protein sequences were aligned using Clustal X . The gaps and ambiguously aligned sites were removed manually. Phylogenetic analyses were performed using a maximum likelihood (ML) approach with PhyML version 3.0  and a neighbor-joining (NJ) method using MEGA . The ML phylogenetic analyses were conducted with the following parameters: JTT model, estimated proportion of invariable sites, 4 rate categories, estimated gamma distribution parameter, and optimized starting BIONJ tree. The JTT model was also employed for the construction of NJ trees. A total of 100 non-parametric bootstrap samplings were carried out to estimate the support level for each internal branch for both the ML and NJ trees. The branch lengths and topologies of all phylogenies were calculated with PhyML. Phylogenetic trees were visualized using the explorer program in MEGA.
Detection of positive selection
A phylogenetically based maximum likelihood method was used to estimate the selective pressure acting on coding regions. The values of the d N /d S ratio (or ω) for the land plant and selected bacterial AEE genes were calculated using the program codeml from PAML v4.4 . The PAL2NAL program  was utilized for conversion of the protein sequence alignment into the corresponding codon-based nucleotide alignment, which, in turn, was input into the codeml program in PAML. Using the codeml program, we detected a variation in ω between sites by employing likelihood ratio tests (LRTs) of M0 vs. M3, M1a vs. M2a, and M7 vs. M8. The LRT for the M0 vs. M3 comparison was used to test the heterogeneity in ω between the codon sites, while the other two LRTs were used to detect the role of positive selection. For one LRT, twice the difference of the log likelihood of the two models was compared with chi-square (χ 2 ) statistics, with degrees of freedom (DFs) equal to the difference in the number of parameters. In our analyses, the DFs were 3 for the M0/M3 test and 2 for the M1a/M2a and M7/M8 tests [28, 42].
An improved branch-site model  was also used to detect the role of positive selection acting on the land plant AEE gene following HGT. For this analysis, we compared the null hypothesis (ω fixed to 1) with the alternative hypothesis (free ω) to test whether positive selection acted on the evolution of land plant AEE genes. A phylogenetic tree was generated using the land plant and bacterial AEE genes with the program PHYML. Here, only the bacterial genes falling within the same branch as the land plant genes were used. The land plant branch was used as the foreground, while the branch containing the genes from bacteria, the putative donors of the land plant AEE gene, was used as the background. The Bayes empirical Bayes procedure  in codeml was used to calculate the posterior probability that each site was subject to positive selection in the foreground branch.
This work was supported by grants from the National Basic Research Program of China (2011CB100100), the Priority Academic Program Development of Jiangsu Higher Education Institutions, the National Natural Science Foundation of China (31200943, 31171187), the Natural Science Foundation of Jiangsu province (BK2012261), the Vital Project of Natural Science of Universities in Jiangsu Province (09KJA210002), the Innovative Research Team of Universities in Jiangsu Province.
- Andersson JO: Lateral gene transfer in eukaryotes. Cell Mol Life Sci. 2005, 62 (11): 1182-1197. 10.1007/s00018-005-4539-z.PubMedView ArticleGoogle Scholar
- Ros VI, Hurst GD: Lateral gene transfer between prokaryotes and multicellular eukaryotes: ongoing and significant?. BMC Biol. 2009, 7: 20. 10.1186/1741-7007-7-20.PubMedPubMed CentralView ArticleGoogle Scholar
- Richardson AO, Palmer JD: Horizontal gene transfer in plants. J Exp Bot. 2007, 58 (1): 1-9.PubMedView ArticleGoogle Scholar
- Lawrence JG, Ochman H: Molecular archaeology of the Escherichia coli genome. Proc Natl Acad Sci USA. 1998, 95 (16): 9413-9417. 10.1073/pnas.95.16.9413.PubMedPubMed CentralView ArticleGoogle Scholar
- Huang J, Gogarten JP: Concerted gene recruitment in early plant evolution. Genome Biol. 2008, 9 (7): R109. 10.1186/gb-2008-9-7-r109.PubMedPubMed CentralView ArticleGoogle Scholar
- Richards TA, Soanes DM, Jones MD, Vasieva O, Leonard G, Paszkiewicz K, Foster PG, Hall N, Talbot NJ: Horizontal gene transfer facilitated the evolution of plant parasitic mechanisms in the oomycetes. Proc Natl Acad Sci USA. 2011, 108 (37): 15258-15263. 10.1073/pnas.1105100108.PubMedPubMed CentralView ArticleGoogle Scholar
- Keeling PJ: Functional and ecological impacts of horizontal gene transfer in eukaryotes. Curr Opin Genet Dev. 2009, 19 (6): 613-619. 10.1016/j.gde.2009.10.001.PubMedView ArticleGoogle Scholar
- Dunning Hotopp JC, Clark ME, Oliveira DC, Foster JM, Fischer P, Munoz Torres MC, Giebel JD, Kumar N, Ishmael N, Wang S: Widespread lateral gene transfer from intracellular bacteria to multicellular eukaryotes. Science. 2007, 317 (5845): 1753-1756. 10.1126/science.1142490.PubMedView ArticleGoogle Scholar
- Ni T, Yue J, Sun G, Zou Y, Wen J, Huang J: Ancient gene transfer from algae to animals: mechanisms and evolutionary significance. BMC Evol Biol. 2012, 12: 83. 10.1186/1471-2148-12-83.PubMedPubMed CentralView ArticleGoogle Scholar
- Suzuki K, Yamashita I, Tanaka N: Tobacco plants were transformed by Agrobacterium rhizogenes infection during their evolution. Plant J. 2002, 32 (5): 775-787. 10.1046/j.1365-313X.2002.01468.x.PubMedView ArticleGoogle Scholar
- Richards TA, Soanes DM, Foster PG, Leonard G, Thornton CR, Talbot NJ: Phylogenomic analysis demonstrates a pattern of rare and ancient horizontal gene transfer between plants and fungi. Plant Cell. 2009, 21 (7): 1897-1911. 10.1105/tpc.109.065805.PubMedPubMed CentralView ArticleGoogle Scholar
- Emiliani G, Fondi M, Fani R, Gribaldo S: A horizontal gene transfer at the origin of phenylpropanoid metabolism: a key adaptation of plants to land. Biol Direct. 2009, 4: 7. 10.1186/1745-6150-4-7.PubMedPubMed CentralView ArticleGoogle Scholar
- Yoshida S, Maruyama S, Nozaki H, Shirasu K: Horizontal gene transfer by the parasitic plant Striga hermonthica. Science. 2010, 328 (5982): 1128. 10.1126/science.1187145.PubMedView ArticleGoogle Scholar
- Yue J, Hu X, Sun H, Yang Y, Huang J: Widespread impact of horizontal gene transfer on plant colonization of land. Nat Commun. 2012, 3: 1152.PubMedPubMed CentralView ArticleGoogle Scholar
- Christin PA, Edwards EJ, Besnard G, Boxall SF, Gregory R, Kellogg EA, Hartwell J, Osborne CP: Adaptive evolution of C(4) photosynthesis through recurrent lateral gene transfer. Curr Biol. 2012, 22 (5): 445-449. 10.1016/j.cub.2012.01.054.PubMedView ArticleGoogle Scholar
- Gulick AM, Schmidt DM, Gerlt JA, Rayment I: Evolution of enzymatic activities in the enolase superfamily: crystal structures of the L-Ala-D/L-Glu epimerases from Escherichia coli and Bacillus subtilis. Biochemistry. 2001, 40 (51): 15716-15724. 10.1021/bi011641p.PubMedView ArticleGoogle Scholar
- Klenchin VA, Schmidt DM, Gerlt JA, Rayment I: Evolution of enzymatic activities in the enolase superfamily: structure of a substrate-liganded complex of the L-Ala-D/L-Glu epimerase from Bacillus subtilis. Biochemistry. 2004, 43 (32): 10370-10378. 10.1021/bi049197o.PubMedView ArticleGoogle Scholar
- Kalyanaraman C, Imker HJ, Fedorov AA, Fedorov EV, Glasner ME, Babbitt PC, Almo SC, Gerlt JA, Jacobson MP: Discovery of a dipeptide epimerase enzymatic function guided by homology modeling and virtual screening. Structure. 2008, 16 (11): 1668-1677. 10.1016/j.str.2008.08.015.PubMedPubMed CentralView ArticleGoogle Scholar
- Saghatelian A, Cravatt BF: Assignment of protein function in the postgenomic era. Nat Chem Biol. 2005, 1 (3): 130-142. 10.1038/nchembio0805-130.PubMedView ArticleGoogle Scholar
- Schuler MA, Duan H, Bilgin M, Ali S: Arabidopsis cytochrome P450s through the looking glass: a window on plant biochemistry. Phytochemistry Rev. 2006, 5 (2–3): 205-237.View ArticleGoogle Scholar
- Liao Y, Zhou X, Yu J, Cao Y, Li X, Kuai B: The key role of chlorocatechol 1,2-dioxygenase in phytoremoval and degradation of catechol by transgenic Arabidopsis. Plant Physiol. 2006, 142 (2): 620-628. 10.1104/pp.106.085936.PubMedPubMed CentralView ArticleGoogle Scholar
- Wodniok S, Brinkmann H, Glockner G, Heidel AJ, Philippe H, Melkonian M, Becker B: Origin of land plants: do conjugating green algae hold the key?. BMC Evol Biol. 2011, 11: 104. 10.1186/1471-2148-11-104.PubMedPubMed CentralView ArticleGoogle Scholar
- Lewis LA, McCourt RM: Green algae and the origin of land plants. Am J Bot. 2004, 91 (10): 1535-1556. 10.3732/ajb.91.10.1535.PubMedView ArticleGoogle Scholar
- de Duve C: The origin of eukaryotes: a reappraisal. Nat Rev Genet. 2007, 8 (5): 395-403. 10.1038/nrg2071.PubMedView ArticleGoogle Scholar
- Sato N: Comparative analysis of the genomes of cyanobacteria and plants. Genome Inform. 2002, 13: 173-182.PubMedGoogle Scholar
- Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24 (8): 1586-1591. 10.1093/molbev/msm088.PubMedView ArticleGoogle Scholar
- Yang Z, Wong WS, Nielsen R: Bayes empirical bayes inference of amino acid sites under positive selection. Mol Biol Evol. 2005, 22 (4): 1107-1118. 10.1093/molbev/msi097.PubMedView ArticleGoogle Scholar
- Nielsen R, Yang Z: Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene. Genetics. 1998, 148 (3): 929-936.PubMedPubMed CentralGoogle Scholar
- Zhang J, Nielsen R, Yang Z: Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level. Mol Biol Evol. 2005, 22 (12): 2472-2479. 10.1093/molbev/msi237.PubMedView ArticleGoogle Scholar
- Long M, Betran E, Thornton K, Wang W: The origin of new genes: glimpses from the young and old. Nat Rev Genet. 2003, 4 (11): 865-875.PubMedView ArticleGoogle Scholar
- Koonin EV: The origin of introns and their role in eukaryogenesis: a compromise solution to the introns-early versus introns-late debate?. Biol Direct. 2006, 1: 22. 10.1186/1745-6150-1-22.PubMedPubMed CentralView ArticleGoogle Scholar
- Klassen JL: Pathway evolution by horizontal transfer and positive selection is accommodated by relaxed negative selection upon upstream pathway genes in purple bacterial carotenoid biosynthesis. J Bacteriol. 2009, 191 (24): 7500-7508. 10.1128/JB.01060-09.PubMedPubMed CentralView ArticleGoogle Scholar
- Yang Z, Gu S, Wang X, Li W, Tang Z, Xu C: Molecular evolution of the CPP-like gene family in plants: insights from comparative genomics of Arabidopsis and rice. J Mol Evol. 2008, 67 (3): 266-277. 10.1007/s00239-008-9143-z.PubMedView ArticleGoogle Scholar
- Kosiol C, Vinar T, da Fonseca RR, Hubisz MJ, Bustamante CD, Nielsen R, Siepel A: Patterns of positive selection in six mammalian genomes. PLoS Genet. 2008, 4 (8): e1000144. 10.1371/journal.pgen.1000144.PubMedPubMed CentralView ArticleGoogle Scholar
- Shen YY, Liang L, Zhu ZH, Zhou WP, Irwin DM, Zhang YP: Adaptive evolution of energy metabolism genes and the origin of flight in bats. Proc Natl Acad Sci USA. 2010, 107 (19): 8666-8671. 10.1073/pnas.0912613107.PubMedPubMed CentralView ArticleGoogle Scholar
- Goodstein DM, Shu S, Howson R, Neupane R, Hayes RD, Fazo J, Mitros T, Dirks W, Hellsten U, Putnam N: Phytozome: a comparative platform for green plant genomics. Nucleic Acids Res. 2012, 40 (Database issue): D1178-D1186.PubMedPubMed CentralView ArticleGoogle Scholar
- Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, Fong JH, Geer LY, Geer RC, Gonzales NR: CDD: a conserved domain database for the functional annotation of proteins. Nucleic Acids Res. 2011, 39 (Database issue): D225-D229.PubMedPubMed CentralView ArticleGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R: Clustal W and clustal X version 2.0. Bioinformatics. 2007, 23 (21): 2947-2948. 10.1093/bioinformatics/btm404.PubMedView ArticleGoogle Scholar
- Guindon S, Delsuc F, Dufayard JF, Gascuel O: Estimating maximum likelihood phylogenies with PhyML. Methods Mol Biol. 2009, 537: 113-137. 10.1007/978-1-59745-251-9_6.PubMedView ArticleGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28 (10): 2731-2739. 10.1093/molbev/msr121.PubMedPubMed CentralView ArticleGoogle Scholar
- Suyama M, Torrents D, Bork P: PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Res. 2006, 34 (Web Server issue): W609-W612.PubMedPubMed CentralView ArticleGoogle Scholar
- Wong WS, Yang Z, Goldman N, Nielsen R: Accuracy and power of statistical methods for detecting adaptive evolution in protein coding sequences and for identifying positively selected sites. Genetics. 2004, 168 (2): 1041-1051. 10.1534/genetics.104.031153.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.