Phylogenetic analyses of Vitis (Vitaceae) based on complete chloroplast genome sequences: effects of taxon sampling and phylogenetic methods on resolving relationships among rosids
© Jansen et al; licensee BioMed Central Ltd. 2006
Received: 09 December 2005
Accepted: 09 April 2006
Published: 09 April 2006
The Vitaceae (grape) is an economically important family of angiosperms whose phylogenetic placement is currently unresolved. Recent phylogenetic analyses based on one to several genes have suggested several alternative placements of this family, including sister to Caryophyllales, asterids, Saxifragales, Dilleniaceae or to rest of rosids, though support for these different results has been weak. There has been a recent interest in using complete chloroplast genome sequences for resolving phylogenetic relationships among angiosperms. These studies have clarified relationships among several major lineages but they have also emphasized the importance of taxon sampling and the effects of different phylogenetic methods for obtaining accurate phylogenies. We sequenced the complete chloroplast genome of Vitis vinifera and used these data to assess relationships among 27 angiosperms, including nine taxa of rosids.
The Vitis vinifera chloroplast genome is 160,928 bp in length, including a pair of inverted repeats of 26,358 bp that are separated by small and large single copy regions of 19,065 bp and 89,147 bp, respectively. The gene content and order of Vitis is identical to many other unrearranged angiosperm chloroplast genomes, including tobacco. Phylogenetic analyses using maximum parsimony and maximum likelihood were performed on DNA sequences of 61 protein-coding genes for two datasets with 28 or 29 taxa, including eight or nine taxa from four of the seven currently recognized major clades of rosids. Parsimony and likelihood phylogenies of both data sets provide strong support for the placement of Vitaceae as sister to the remaining rosids. However, the position of the Myrtales and support for the monophyly of the eurosid I clade differs between the two data sets and the two methods of analysis. In parsimony analyses, the inclusion of Gossypium is necessary to obtain trees that support the monophyly of the eurosid I clade. However, maximum likelihood analyses place Cucumis as sister to the Myrtales and therefore do not support the monophyly of the eurosid I clade.
Phylogenies based on DNA sequences from complete chloroplast genome sequences provide strong support for the position of the Vitaceae as the earliest diverging lineage of rosids. Our phylogenetic analyses support recent assertions that inadequate taxon sampling and incorrect model specification for concatenated multi-gene data sets can mislead phylogenetic inferences when using whole chloroplast genomes for phylogeny reconstruction.
The estimation of phylogenetic relationships among angiosperms has received considerable attention during the past decade with the rapid increase in availability of DNA sequence data from a wide diversity of markers and taxa [reviewed in ]. Most previous molecular phylogenetic studies of flowering plants have relied on one to several genes from the chloroplast, mitochondrial, and/or nuclear genomes, though most of these analyses were based on chloroplast markers. These efforts have resolved the relationships among many of the major lineages of angiosperms but a number of outstanding issues remain . Completely sequenced chloroplast genomes provide a rich source of data that can be used to address phylogenetic questions at deep nodes in the angiosperm tree [2–6]. The use of DNA sequences from all of the shared chloroplast genes provides many more characters for phylogeny reconstruction compared to previous studies that have relied on only one or a few genes to address the same questions. However, the whole genome approach can result in misleading estimates of relationships because of limited taxon sampling [5, 7–10] and the use of incorrect models of sequence evolution in concatenated datasets [4, 11]. Thus, there is a growing interest in expanding the taxon sampling of complete chloroplast genome sequences and developing new evolutionary models for phylogenetic analysis of chloroplast sequences  to overcome these concerns.
The rosids represent the largest of the eight major clades of core eudicots and include nearly one third of all flowering plants. Single and multi-gene phylogenies of rosids have identified seven major clades, however, relationships among these clades remain unresolved [13–16]. One of these unresolved clades includes the Vitaceae, which includes grape, an important crop plant. The phylogenetic position of Vitaceae has been controversial for many years. Some previous classifications place the family within the Rhamnales in the subclass Rosidae . More recent molecular phylogenies based on one to four genes provided weak support for the placement of Vitaceae sister to the Caryophylales , asterids , Saxifragales , Dilleniaceae , or to the rosids [14–16]. Thus, the phylogenetic relationship of the grape family to core eudicots remains unresolved.
In this article, we report on the complete sequence of the chloroplast genome of grape (Vitis vinifera, Vitaceae). In addition to describing the organization of the chloroplast genome, we present results of phylogenetic analyses of DNA sequences for 61 genes from grape and 26 other angiosperm chloroplast genomes, including eight other members of the rosid clade. The phylogenetic analyses provide insights into the relationship of Vitaceae to other rosids and illustrate the importance of taxon sampling and analytical method on addressing phylogenetic questions using whole genome sequences. The complete chloroplast genome sequence of Vitis also provides valuable data for using chloroplast genetic engineering for this economically important crop plant .
Size, gene content, order and organization of the grape chloroplast genome
Location of repeats in the grape chloroplast genome. Repeats 1 to 15 are direct, and 16 to 36 are inverted. Table includes repeats at least 30 bp in size, with a sequence identity ≥ 90%. IGS = Intergenic spacer. See Figure 1 for location of repeats on the gene map. Repeats indicated by * (palindrome) and ** (tandem) are only shown once on the circular map in Figure 1.
ycf3 intron, IGS
IGS (4 bp) – trnS-GCU, IGS (4 bp) – trnS-UGA
ycf3 intron, IGS, ndhA intron
psaB exon – psaA exon
IGS, ndhA intron
IGS (3 bp) – trnS-GCU, trnS-GGA
IGS (2 bp) – trnS-UGA, trnS-GGA
ycf3 intron, IGS
ndhA intron, IGS
RNA variable sites in grape chloroplast transcripts
Differences observed by comparison of grape chloroplast genome sequences with EST sequences obtained by BLAST search in Genbank.
Gene size (bp)
Number of variable sites
Amino acid change
Taxa included in phylogenetic analyses with GenBank accession numbers and references.
GenBank Accession Numbers
Gymnosperms – Outgroups
Wakasugi et al. 1994 
Leebens-Mack et al 2005 
Goremykin et al. 2003 
Leebens-Mack et al 2005 
Goremykin et al. 2004 
Leebens-Mack et al 2005 
Hiratsuka et al. 1989 
Asano et al. 2004 
Ikeo and Ogihara, unpublished
Leebens-Mack et al 2005 
Leebens-Mack et al 2005 
Maier et al. 1995 
Goremykin et al. 2003 
Sato et al. 1999 
Schmitz-Linneweber et al. 2002 
Plader et al. unpublished
Steane 2005 
Saski et al. 2005 
Lee et al. 
Kato et al. 2000 
Lin et al., unpublished
Shinozaki et al. 1986 
Hupfer et al. 2000 
Kim and Lee 2004 
Leebens-Mack et al 2005 
Daniell et al. 
Daniell et al. 
Schmitz-Linneweber et al. 2001 
MP analysis of the second dataset of 29 taxa including Gossypium resulted in a single most parsimonious tree with a length of 51,056, a consistency index of 0.46 (excluding uninformative characters) and a retention index of 0.61 (Fig. 4A). Bootstrap analyses indicated that 24 of the 26 nodes were supported by values ≥ 95%, and all but four of these nodes had a bootstrap value of 100%. ML analysis resulted in a tree with a – lnL = 296670.545 (Fig. 4B). ML bootstrap analyses indicated that 22 of the 26 nodes were supported by values ≥ 95% and all but two of these nodes had a bootstrap value of 100%. Both MP and ML analyses provided strong support for Vitis as the earliest diverging lineage of rosids, monophyly of Myrtales, and sister relationship of Brassicales and Malvales. The ML and MP trees had three important topological differences. The first two differences concerned the position of Calycanthus and the basal angiosperms, which were identical to those described above for the analyses that excluded Gossypium. The other difference concerned relationships among rosids. The MP tree (Fig. 4A) showed strong support (100% bootstrap) for monophyly of the eurosid I clade because of the sister relationship between the Fabales and Cucurbitales. In contrast, the ML tree indicated that eurosids I are paraphyletic because Cucurbitales were sister to Mrytales rather than Fabales (Fig. 4B); bootstrap support for this relationship was also strong (92%).
Grapes are an important crop plant grown for wine, juice, raisins, and as fresh fruit. In 2004, the world's grape harvest area in 89 grape-producing countries was 7.5 million hectares, and in the United States grapes were grown in 380,000 hectares . The total production of grapes in the US in 2004 was 5,418,160 metric tons and this generated $2.5 billion . There is considerable interest in using chloroplast genetic engineering as an environmentally friendly approach for engineering disease resistance to powdery and downy mildew, two fungal diseases that have a negative impact on the grape industry. Chloroplast genetic engineering offers a number of unique advantages, including a high-level of transgene expression , multi-gene engineering in a single transformation event [23–26], transgene containment via maternal inheritance [27–29] or cytoplasmic male sterility , and lack of gene silencing, position effect, pleiotropic effects, and undesirable foreign DNA [20, 31–35]. Thus far, transgenes have been stably integrated and expressed via the chloroplast genome to confer several useful agronomic traits, including insect resistance [36, 37, 23], herbicide resistance [27, 38], disease resistance , drought tolerance , salt tolerance , and phytoremediation . The complete grape chloroplast genome sequence reported in this paper provides valuable characterization of spacer regions for potential integration of transgenes at optimal sites via homologous recombination, as well as endogenous regulatory sequences for optimal expression of transgenes.
Genome organization and evolution
The organization of the Vitis genome with two copies of an IR separating the SSC and LSC regions is identical to most sequenced angiosperm chloroplast genomes [reviewed in ]. The size of the genome at 160,928 bp is also within the known size range for angiosperms, which generally vary from 150,519 (Lotus ) to 162,686 bp (Amborella ) among photosynthetic genomes from dicots that have both copies of the IR. Size of the Vitis IR at 26,358 bp is also well within the size range of other sequenced dicot genomes, which range from 23,302 (Calycanthus ) to 27,807 bp (Oenothera ). Gene content and order of the Vitis chloroplast genome is virtually identical to tobacco and many other unrearranged angiosperm chloroplast genomes. Several previously sequenced rosid chloroplast genomes have lost the rpl22 gene, including legumes [45–49]. The distribution of this loss on the chloroplast phylogeny (arrows in Fig. 4A) indicates that there have been at least two independent losses of rpl22 in rosids. Multiple, independent gene losses in angiosperms have been demonstrated for other genes including infA , rps16  and accD [51, 52]. Thus, it is evident that gene losses are not always reliable indicators of phylogenetic relationships.
It is increasingly evident that chloroplast genomes contain repeated sequences other than the IR . Several studies have identified a higher incidence of dispersed repeats in genomes that have experienced extensive rearrangements [53, 54]. However, dispersed repeats are also being detected in unrearranged genomes. In most cases, these repeats are more common in intergenic spacers and introns, which is also true for the Vitis genome. Repeats have been located in a number of other rosids  in the same regions as those identified in the Vitis genome. One of these, a 32 bp repeat in the trnS gene, is in the same location in Gossypium hirsutum , indicating that this repeat may be shared among rosids. Although repeats have been implicated in playing a role in chloroplast genome rearrangements , their effect if any in unrearranged chloroplast genomes is unknown.
Based on previous studies of Atropa  and tobacco , posttranscriptional RNA editing events, as well as deamination-facilitating attacks on nucleotides' exocyclic amino groups, yield primarily C-to-U alterations. Analyses of the Vitis chloroplast genome and the corresponding ESTs indicate that the five C-to-U changes likely represent mRNA edits. However, the remaining six differences could be either sequencing errors in the genomic DNA or EST sequences or due to the use of different cultivars and/or plants/tissues used for sequencing. Our methods eliminate the latter explanation since we only compared DNA and EST sequences from leaves of the chardonnay variety of Vitis vinifera. In view of the high depth of coverage (8X) of our genomic DNA sequences, we believe that the non C-to-U changes represent EST sequencing errors.
Evolutionary loss of RNA editing sites has been observed in earlier studies and could be attributed to a decrease in the effect of RNA-editing enzymes . Additionally, conversions other than C-to-U in Vitis and other plants suggest that chloroplast genomes may be accumulating a considerable number of nucleotide substitutions, and some genes might accumulate more changes than others, such as the petL and ndh genes that have a high frequency of RNA editing . Therefore, despite high levels sequence conservation in chloroplast genomes, variations do occur posttranscriptionally, promoting translational efficiency due to transcript-protein complex binding and/or changes in the chloroplasts microenvironment (e.g., like redox potential or light intensity [61, 62]).
Phylogenetic analyses of 28 (Fig. 3) or 29 (Fig. 4) angiosperms based on 61 protein-coding genes identified many of the major lineages recognized in previous phylogenetic hypotheses of flowering plants [reviewed in ]. Two groups, Amborella and Nymphaelaes (represented by Nuphar and Nymphaea) are basal, with Amborella forming the first diverging lineage in MP analyses and Amborella/Nymphaelaes together forming the most basal clade in ML trees. These results are congruent with recent 61-gene analyses by Leebens-Mack et al.  and support their contention that limited taxon sampling in earlier whole chloroplast genome phylogenies led some previous workers to suggest that Amborella may not be among the earliest diverging angiosperm lineages [2, 3]. Monophyly of the monocots is strongly supported, and they are sister to the remaining angiosperms. Calycanthus, the sole representative of the magnolids, is weakly supported as sister to eudicots in the MP analyses (Figs. 3A and 4A) but the genus is weakly supported as sister to a clade that includes both monocots and eudicots in ML trees (Figs. 3B and 4B). Monophyly of eudicots is strongly supported (100% bootstrap values), in agreement with phylogenies based on both pollen [63, 64] and other molecular data [13, 14, 18, 19, 65–67]. Within eudicots, Ranunculales diverge first and are sister to a strongly supported eudicot clade that includes two moderately to well-supported groups comprising the rosids and asterids. The early divergence of Ranunculales among eudicots is in agreement with many recent molecular phylogenies [see chapter 5 in ]. Although previous studies have clearly indicated that Carylophyllales belong in the core eudicot clade , resolution of the relationships of Caryophyllales to other major clades of eudicots remains uncertain. This order has been considered to be closely allied to rosids, asterids, or simply as an unresolved major eudicot clade sister to the Dilleniaceae . Although taxon sampling is limited in our 61 gene phylogeny, there is moderate to strong support for a sister relationship between the Caryophyllales and asterids (Figs. 3 and 4).
The rosid clade is very diverse, including nearly 140 families representing approximately 39% of the species of angiosperms. The most recent phylogenies of this group [summarized in chapter 8 in ] indicate that there are seven major clades whose relationships still remain unresolved. Eight (Fig. 3) or nine (Fig. 4, includes Gossypium) representatives of four of these major clades are included in our phylogenetic analyses, including members of eurosids I, eurosids II, Myrtales, and Vitaceae. Phylogenetic analyses of both datasets using MP and ML clearly indicate that the Vitaceae is sister to the remaining rosids, and therefore represents an early diverging member of the rosid clade. Previous molecular phylogenetic comparisons that included Vitaceae could not resolve its relationship. Phylogenetic analyses of rbcL sequences alone placed the Vitaceae as sister to either the Caryophyllales or asterid clade with weak support . Phylogenies based on atpB provided only weak support for a sister relationship of Vitaceae to Saxifragales . Several phylogenies based on two to four genes suggested that the Vitaceae are sister to rest of the rosids, with relatively weak support (50–75%; [14–16]). However, phylogenies based on the chloroplast gene matK did not place Vitaceae sister to rosids but instead positioned the family as sister to Dilleniaceae with weak support . In short, the phylogenetic position of Vitaceae is equivocal, though our results strongly support the earlier findings that Vitaceae represent an early diverging clade within rosids (Figs. 3 and 4).
The two datasets we examined differed by only one taxon but the results of MP analyses differed dramatically regarding the placement of three of the four rosid clades examined (compare Figs. 3A and 4A). The 28-taxon dataset (excluding Gossypium) showed relationships that are incongruent with recent molecular phylogenies of rosids  by placing the eurosids II (represented by only Brassicales) sister to the Fabales in eurosids I. This made eurosids I paraphyletic because the other representative of this clade is Cucumis (Cucurbitales), which is sister to the Brassicales in molecular phylogenies of rosids . The addition of Gossypium in the 29-taxon dataset generates MP trees (Fig. 4A) that are congruent with previous angiosperm phylogenies. The Brassicales and Malvales are sister and there is strong support for the monophyly of eurosid II. The addition of Gossypium also makes the eurosid I clade strongly monophyletic in the MP tree by placing the Cucurbitales sister to the Fabales, both of which are members of the nitrogen-fixing clade [see chapter 8 in ]. In contrast to the MP trees, relationships among the major rosid clades do not differ in the ML trees when Gossypium is added. In both the 28 and 29-taxon data sets the ML trees do not support the monophyly of eurosids I since Cucurbitales (eurosid I) are sister to the Myrtales and Brassicales (eurosid II) are sister to the Fabales (eurosid I). Therefore, the ML analyses are incongruent with currently accepted relationships among rosids , though the strongest support for the monophyly of eurosid I clade is only 77% (jackknife support) in a three-gene analysis . Thus, our results suggest that additional phylogenetic studies are needed to assess the monophyly of eurosids I and their relationship to other rosids.
There has been considerable debate regarding the utility of whole genome sequences for phylogeny reconstruction [5, 7–10]. Some have argued that the use of more genes from whole genomes has great potential for providing much more data for resolving phylogenetic relationships [2, 68], whereas others have suggested that problems with limited taxon sampling available for whole genomes [5, 7, 8, 10] and model misspecification [4, 11] overshadows any potential advantages. One example that highlighted each of these concerns centered around the controversy regarding identification of basal angiosperms. Leebens-Mack et al.  demonstrated that inadequate taxon sampling clearly played a role in misleading some previous studies, and Goremykin et al.  demonstrated that ML analyses of whole chloroplast genome data sets can be sensitive to model specification. It is well known that ML methods fail when model parameters are misspecified [69–71]. The phylogenetic analyses in this study provide yet another example of these phenomena. Addition of the Gossypium genome to our parsimony analyses generated trees that are congruent with current understanding of relationships among the major rosid clades. However, the ML analyses are incongruent with the MP trees regarding the monophyly and relationships of the rosid clades and support for the alternative relationships was very strong in each case (compare Figs. 4A and 4B). It is possible, if not likely, that the use of a single "average" model (GTR + I + Γ) in the ML analyses is inappropriate for a data set of 61 concatenated genes [see  for a discussion of this issue]. Future phylogenetic analyses of complete chloroplast genome sequences should consider using methods in which different models can be applied to different partitions of the data (e.g., genes, codon positions, functional groups) . Development of more appropriate models of evolution of chloroplast sequences  may also improve the accuracy of phylogenies based on these genomes. Thus, we need more extensive sampling of whole chloroplast genomes from the major lineages of flowering plants and more rigorous phylogenetic analyses before the full potential of this approach can be realized. Ongoing projects by several labs [see  for a list of some of these] should greatly enhance our taxon sampling so that we can generate reliable phylogenies based on whole chloroplast genomes.
The chloroplast genome of Vitis has a very similar size and organization to other previously sequenced, unrearranged angiosperm chloroplast genomes. These sequences will provide a valuable resource for developing transgenes for this important crop plant using the more environmentally friendly chloroplast genetic engineering technology . Phylogenetic analyses of a dataset of sequences of 61 shared protein-coding genes of Vitis and 26 other angiosperm genomes demonstrated the importance of taxon sampling and methods of phylogenetic analysis for phylogenomic studies. Furthermore, trees generated by both parsimony and likelihood methods provided support for the resolution of relationships among eudicots. This included support for the position of the Ranunculales as the earliest diverging lineage of eudicots, the placement of the Caryophyllales as the sister clade to the asterids, and the position of the Vitaceae sister to all other rosids. However, resolution of relationships among the remaining rosid clades based on complete chloroplast genome sequences remains unresolved due to limited taxon sampling and differences in trees generated by MP and ML analyses.
The bacterial artificial chromosome (BAC) library of grape was constructed by ligating size-fractionated partial Hind III digests of total cellular, high molecular weight DNA with the pINDIGOBAC vector. The average insert size of the grape library is 144 kb. BAC-related resources for this public library can be obtained online from the Clemson University Genomics Institute BAC/EST Resource Center .
BAC clones containing the chloroplast genome inserts were isolated by screening the library with a soybean chloroplast probe. The first 96 positive clones from screening were pulled from the library, arrayed in a 96-well microtitre plate, copied, and archived. Selected clones were then subjected to Hind III fingerprinting and Not I digests. End-sequences were determined and localized on the chloroplast genome of Arabidopsis thaliana to deduce the relative positions of the clones, then clones that covered the entire chloroplast genome of grape were chosen for sequencing.
DNA sequencing and genome assembly
The nucleotide sequences of the BAC clones were determined by the bridging shotgun method. The purified BAC DNA was subjected to hydroshearing, end repair, and then size-fractionated by agarose gel electrophoresis. Fractions of approximately 3.0–5.0 kb were eluted and ligated into the vector pBLUESCRIPT IIKS+. The libraries were plated and arrayed into 40 96-well microtitre plates for the sequencing reactions.
Sequencing was performed using the Dye-terminator cycle sequencing kit (Perkin Elmer Applied Biosystems, USA). Sequence data from the forward and reverse priming sites of the shotgun clones were accumulated. Sequence data equivalent to eight times the size of the genome was assembled using Phred/Phrap programs .
The Vitis genome was annotated using DOGMA (Dual Organellar GenoMe Annotator) , after uploading a FASTA-formatted file of the complete plastid genome to the program's server. BLASTX and BLASTN searches against a custom database of previously published plastid genomes identified Vitis' putative protein-coding genes, and tRNAs or rRNAs. For genes with low sequence identity, manual annotation was performed, after identifying the position of the start and stop codons, as well as the translated amino acid sequence, using the plastid/bacterial genetic code.
Examination of repeat structure
REPuter  was used in order to locate and count the direct (forward) and inverted (palindromic) repeats within the grape chloroplast genome. For repeat identification, the following constraints were set to REPuter: (i) minimum repeat size of 30 bp, and (ii) 90% or greater sequence identity, based on Hamming distance of 3 . Manual verification of the identified repeats was performed in EditSeq, while performing intragenomic blast search of the identified repeat sequence.
Variation between coding sequences and cDNAs
Each of the gene sequences from the grape chloroplast genome was used to perform a BLAST search of expressed sequence tags (ESTs) from Genbank. In order to incorporate specificity into our analyses, the matching ESTs had to meet all of the following criteria: (1) belong to a Vitis vinifera cv., (2) belong to the chardonnay variety, and (3) come from leaf tissue. Due to length variations between the screened ESTs and the related gene, the retrieved EST with the highest bit score was selected for further analyses. The retrieved Vitis vinifera EST was aligned with the corresponding annotated gene using ClustalX , followed by screening for nucleotide and amino acid changes using Megalign and the plastid/bacterial genetic code. Because of variations in the length between an EST and the related gene, the length of the analyzed sequence was recorded.
The 61 genes included in the analyses of Goremykin et al. [2, 3] and Leebens-Mack et al.  were extracted from our new chloroplast genome sequences of Vitis using the organellar genome annotation program DOGMA . The same set of 61 genes was extracted from chloroplast genome sequences of six other recently sequenced angiosperm chloroplast genomes, including tomato, potato, soybean, cotton, cucumber, and Eucalyptus (see Table 3 for complete list of genomes examined). In general, alignment of the DNA sequences was straightforward and simply involved adding the 61 genes for the new angiosperms to the aligned data matrix from Leebens-Mack et al. . In some cases, small in-frame insertions or deletions were required for correct alignment. For two genes, ccsA and matK, the DNA sequences were more divergent, requiring alignment using ClustalX  followed by manual adjustments. The complete nucleotide alignment is available online at .
Phylogenetic analyses using maximum parsimony (MP) and maximum likelihood (ML) were performed using PAUP* version 4.10  on two data sets, one including 28 taxa and a second including 29 taxa by the addition of Gossypium. Phylogenetic analyses excluded gap regions. All MP searches included 100 random addition replicates and TBR branch swapping with the Multrees option. Modeltest 3.7  was used to determine the most appropriate model of DNA sequence evolution for the combined 61-gene dataset. Hierarchical likelihood ratio tests and the Akaikle information criterion were used to assess which of the 56 models best fit the data, which was determined to be GTR + I + Γ by both criteria. For ML analyses we performed an initial parsimony search with 100 random addition sequence replicates and TBR branch swapping, which resulted in a single tree. Model parameters were optimized onto the parsimony tree. We fixed these parameters and performed a ML analysis with three random addition sequence replicates and TBR branch swapping. The resulting ML tree was used to re-optimize model parameters, which then were fixed for another ML search with three random addition sequence replicates and TBR branch swapping. This successive approximation procedure was repeated until the same tree topology and model parameters were recovered in multiple, consecutive iterations. This tree was accepted as the final ML tree (Figs. 3B, 4B). Successive approximation has been shown to perform as well as full-optimization analyses for a number of empirical and simulated datasets . Non-parametric bootstrap analyses  were performed for MP analyses with 1000 replicates with TBR branch swapping, 1 random addition replicate, and the Multrees option and for ML analyses with 100 replicates with NNI branch swapping, 1 random addition replicate, and the Multrees option.
small single copy
large single copy, bp, base pair
hypothetical chloroplast reading frame
maximum likelihood, EST, expressed sequence tags
Investigations reported in this article were supported in part by grants from USDA 3611-21000-017-00D and NIH R01 GM 63879 to Henry Daniell and from NSF DEB 0120709 to Robert K. Jansen. We thank R. Haberle for comments on an earlier version of the manuscript.
- Soltis DE, Soltis PS, Endress PK, Chase MW: Phylogeny and evolution of Angiosperms. 2005, Sunderland Massachusetts: Sinauer Associates IncGoogle Scholar
- Goremykin VV, Hirsch-Ernst KI, Wolfl S, Hellwig FH: The chloroplast genome of Nymphaea alba: whole-genome analyses and the problem of identifying the most basal angiosperm. Mol Biol Evol. 2004, 21: 1445-1454. 10.1093/molbev/msh147.View ArticlePubMedGoogle Scholar
- Goremykin VV, Hirsch-Ernst KI, Wolfl S, Hellwig FH: Analysis of the Amborella trichopoda chloroplast genome sequence suggests that Amborella is not a basal angiosperm. Mol Biol Evol. 2003, 20: 1499-1505. 10.1093/molbev/msg159.View ArticlePubMedGoogle Scholar
- Goremykin VV, Holland B, Hirsch-Ernst KI, Hellwig FH: Analysis of Acorus calamus chloroplast genome and its phylogenetic implications. Mol Biol Evol. 2005, 22: 1813-1822. 10.1093/molbev/msi173.View ArticlePubMedGoogle Scholar
- Leebens-Mack J, Raubeson LA, Cui L, Kuehl J, Fourcade M, Chumley T, Boore JL, Jansen RK, dePamphilis CW: Identifying the basal angiosperms in chloroplast genome phylogenies: sampling one's way out of the Felsenstein zone. Mol Biol Evol. 2005, 22: 1948-1963. 10.1093/molbev/msi191.View ArticlePubMedGoogle Scholar
- Chang C-C, Lin H-C, Lin I-P, Chow T-Y, Chen H-H, Chen W-H, Cheng C-H, Lin C-Y, Liu S-M, Chang C-C, Chaw S-M: The chloroplast genome of Phalaenopsis aphrodite (Orchidaceae): comparative analysis of evolutionary rate with that of grasses and its phylogenetic implications. Mol Biol Evol. 2006, 23: 279-291. 10.1093/molbev/msj029.View ArticlePubMedGoogle Scholar
- Soltis DE, Soltis PS: Amborella not a "basal angiosperm"? Not so fast. Amer J Bot. 2004, 91: 997-1001.View ArticleGoogle Scholar
- Stefanovic S, Rice DW, Palmer JD: Long branch attraction, taxon sampling, and the earliest angiosperms: Amborella or monocots?. BMC Evol Biol. 2004, 4: 35-10.1186/1471-2148-4-35.PubMed CentralView ArticlePubMedGoogle Scholar
- Martin W, Deusch O, Stawski N, Grunheit N, Goremykin V: Chloroplast genome phylogenetics: why we need independent approaches to plant molecular evolution. Trends Plant Sci. 2005, 10: 203-209. 10.1016/j.tplants.2005.03.007.View ArticlePubMedGoogle Scholar
- Soltis DE, Albert VA, Savolainen V, Hilu K, Qiu Y-Q, Chase MW, Farris JS, Stefanović S, Rice DW, Palmer JD, Soltis PS: Genome-scale data, angiosperm relationships, and 'ending incongruence': a cautionary tale in phylogenetics. Trends Plant Sci. 2004, 9: 477-483. 10.1016/j.tplants.2004.08.008.View ArticlePubMedGoogle Scholar
- Lockhart PJ, Penny D: The place of Amborella within the radiation of angiosperms. Trends Plant Sci. 2005, 10: 201-202. 10.1016/j.tplants.2005.03.006.View ArticlePubMedGoogle Scholar
- Ané C, Burleigh JG, McMahon MM, Sanderson MJ: Covarion structure in plastid genome evolution: a new statistical test. Mol Biol Evol. 2005, 22: 914-924. 10.1093/molbev/msi076.View ArticlePubMedGoogle Scholar
- Savolainen V, Fay MF, Albach DC, Backlund A, van der Bank M, Cameron KM, Johnson SA, Lledó MD, Pintaud J-C, Powell M, Sheahan MC, Soltis DE, Soltis PS, Weston P, Whitton WM, Wurdack KJ, Chase MW: Phylogeny of the eudicots: a nearly complete familial analysis based on rbcL gene sequences. Kew Bulletin. 2000, 55: 257-309.View ArticleGoogle Scholar
- Savolainen V, Chase MW, Morton CM, Soltis DE, Bayer C, Fay MF, De Bruijn A, Sullivan S, Qiu Y-L: Phylogenetics of flowering plants based upon a combined analysis of plastid atpB and rbcL gene sequences. Syst Biol. 2000, 49: 306-362. 10.1080/10635159950173861.View ArticlePubMedGoogle Scholar
- Soltis DE, Soltis PS, Chase MW, Mort ME, Albach DC, Zanis M, Savolainen V, Hahn WJ, Hoot SB, Fay MF, Axtell M, Swensen SM, Prince LM, Kress WJ, Nixon KC, Farris JS: Angiosperm phylogeny inferred from 18S rDNA, rbcL, and atpB sequences. Bot J Linn Soc. 2000, 133: 381-461. 10.1006/bojl.2000.0380.View ArticleGoogle Scholar
- Soltis DE, Senters AE, Zanis MJ, Kim S, Thompson JD, Soltis PS, Ronse De Craene LP, Endress PK, Farris KS: Gunnerales are sister to other core eudicots: implications for the evolution of pentamery. Amer J Bot. 2003, 90: 461-470.View ArticleGoogle Scholar
- Cronquist A: An integrated system of classification of flowering plants. 1981, Boston Massachusetts: Columbia University PressGoogle Scholar
- Chase M, Soltis D, Olmstead R, Morgan D, Les D, Mishler B, Duvall M, Price R, Hills H, Qui Y-L, Kron K, Rettig J, Conti E, Palmer J, Manhart J, Sytsma K, Michaels H, Kress J, Karol K, Clark D, Hedren M, Gaut B, Jansen R, Kim K-J, Wimpee C, Smith J, Furnier G, Straus S, Xiang Q-Y, Plunkett G, Soltis P, Swensen S, Williams S, Gadek P, Quinn C, Equiarte L, Golenberg E, Learn G, Graham S, Barrett S, Dayanandan S, Albert V: Phylogenetics of seed plants: an analysis of nucleotide sequences from the plastid gene rbc L. Ann Missouri Bot Gard. 1993, 80: 528-580. 10.2307/2399846.View ArticleGoogle Scholar
- Hilu KW, Borsch T, Muller K, Soltis DE, Soltis PS, Savolainen V, Chase M, Powell M, Alice L, Evans R, Sauquet H, Neinhuis C, Slotta T, Rohwer J, Chatrou L: Inference of angiosperm phylogeny based on matK sequence information. Amer J Bot. 2003, 90: 1758-1776.View ArticleGoogle Scholar
- Daniell H, Kumar S, Duformantel N: Breakthrough in chloroplast genetic engineering of agronomically important crops. Trends Biotechnol. 2005, 23: 238-245. 10.1016/j.tibtech.2005.03.008.PubMed CentralView ArticlePubMedGoogle Scholar
- The Food and Agriculture Organization of the United Nations: FAO Statitstical Databases – Agricultural Production. [http://faostat.fao.org/]
- Foreign Agricultural Services, US Department of Agriculture: World Horticultural Trade and U.S. Export Opportunities. 2005, [http://www.fas.usda.gov/htp/horticulture/Grapes/]Google Scholar
- DeCosa B, Moar W, Lee SB, Miller M, Daniell H: Overexpression of the Bt Cry2Aa2 operon in chloroplasts leads to formation of insecticidal crystals. Nat Biotechnol. 2001, 9: 71-74.Google Scholar
- Ruiz ON, Hussein H, Terry N, Daniell H: Phytoremediation of organomercurial compounds via chloroplast genetic engineering. Plt Phys. 2003, 32: 1344-1352. 10.1104/pp.103.020958.View ArticleGoogle Scholar
- Lossl A, Eibl C, Harloff HJ, Jung C, Koop HU: Polyester synthesis in transplastomic tobacco (Nicotiana tabacum L.): significant contents of polyhydroxybutyrate are associated with growth reduction. Plant Cell Rep. 2003, 21: 891-899.PubMedGoogle Scholar
- Quesada-Vargas T, Ruiz ON, Daniell H: Characterization of heterologous multigene operons in transgenic chloroplasts: transcription, processing, translation. Plt Physiol. 2005, 138: 1746-1762. 10.1104/pp.105.063040.View ArticleGoogle Scholar
- Daniell H, Datta R, Varma S, Gray S, Lee SB: Containment of herbicide resistance through genetic engineering of the chloroplast genome. Nat Biotechnol. 1998, 16: 345-348. 10.1038/nbt0498-345.View ArticlePubMedGoogle Scholar
- Daniell H: Molecular strategies for gene containment in transgenic crops. Nat Biotechnol. 2002, 20: 581-586. 10.1038/nbt0602-581.PubMed CentralView ArticlePubMedGoogle Scholar
- Hagemann R: The Sexual Inheritance of Plant Organelles. Molecular Biology and Biotechnology of Plant Organelles. Edited by: Daniell H, Chase C. 2004, The Netherlands: Springer Publishers, 93-113.View ArticleGoogle Scholar
- Ruiz ON, Daniell H: Engineering cytoplasmic male sterility via the chloroplast genome. Plt Physiol. 2005, 138: 1232-1246. 10.1104/pp.104.057729.View ArticleGoogle Scholar
- Lee SB, Kwon HB, Kwon SJ, Park SC, Jeong MJ, Han SE, Daniell H: Accumulation of trehalose within transgenic chloroplasts confers drought tolerance. Mol Breed. 2003, 11: 1-13. 10.1023/A:1022100404542.View ArticleGoogle Scholar
- Dhingra A, Portis AR, Daniell H: Enhanced translation of a chloroplast expressed rbcS gene restores SSU levels and photosynthesis in nuclear antisense rbcS plants. Proc Natl Acad Sci USA. 2004, 101: 6315-6320. 10.1073/pnas.0400981101.PubMed CentralView ArticlePubMedGoogle Scholar
- Daniell H, Lee SB, Panchal T, Wiebe PO: Expression of cholera toxin B subunit gene and assembly as functional oligomers in transgenic tobacco chloroplasts. J Mol Biol. 2001, 311: 1001-1009. 10.1006/jmbi.2001.4921.PubMed CentralView ArticlePubMedGoogle Scholar
- Leelavathi S, Gupta N, Maiti S, Ghosh A, Reddy VS: Overproduction of an alkali- and thermo-stable xylanase in tobacco chloroplasts and efficiency recovery of the enzyme. Mol Breed. 2003, 11: 59-67. 10.1023/A:1022168321380.View ArticleGoogle Scholar
- Grevich JJ, Daniell H: Chloroplast genetic engineering: recent advances and future perspectives. Crit Rev Plt Sci. 2005, 24: 83-108. 10.1080/07352680590935387.View ArticleGoogle Scholar
- McBride KE, Svab Z, Schaaf DJ, Hogan PS, Stalker DM, Maliga P: Amplification of a chimeric Bacillus gene in chloroplasts leads to an extraordinary level of an insecticidal protein in tobacco. Bio/Technology. 1995, 13: 362-365. 10.1038/nbt0495-362.View ArticlePubMedGoogle Scholar
- Kota M, Daniel H, Varma S, Garczynski SF, Gould F, William MJ: Overexpression of the Bacillus thuringiensis (Bt) Cry2Aa2 protein in chloroplasts confers resistance to plants against susceptible and Bt-resistant insects. Proc Natl Acad Sci USA. 1999, 96: 1840-1845. 10.1073/pnas.96.5.1840.PubMed CentralView ArticlePubMedGoogle Scholar
- Iamtham S, Day A: Removal of antibiotic resistance genes from transgenic tobacco plastids. Nat Biotechnol. 2000, 18: 1172-1176. 10.1038/81161.View ArticlePubMedGoogle Scholar
- DeGray G, Rajasekaran K, Smith F, Sanford J, Daniell H: Expression of an antimicrobial peptide via the chloroplast genome to control phytopathogenic bacteria and fungi. Plt Phys. 2001, 127: 852-862. 10.1104/pp.127.3.852.View ArticleGoogle Scholar
- Kumar S, Dhingra A, Daniell H: Plastid expressed betaine aldehyde dehydrogenase gene in carrot cultured cells, roots and leaves confers enhanced salt tolerance. Plant Physiol. 2004, 136: 2843-2854. 10.1104/pp.104.045187.PubMed CentralView ArticlePubMedGoogle Scholar
- Raubeson LA, Jansen RK: Chloroplast genomes of plants. Diversity and Evolution of Plants-Genotypic and Phenotypic Variation in Higher Plants. Edited by: Henry H. 2005, Wallingford: CABI Publishing, 45-68.View ArticleGoogle Scholar
- Kato T, Kaneko T, Sato S, Nakamura Y, Tabata S: Complete structure of the chloroplast genome of a legume, Lotus japonicus. DNA Res. 2000, 7: 323-330. 10.1093/dnares/7.6.323.View ArticlePubMedGoogle Scholar
- Goremykin VV, Hirsch-Ernst KI, Wolfl S, Hellwig FH: The chloroplast genome of the "basal" angiosperm Calycanthus fertilis – structural and phylogenetic analyses. Plt Syst Evol. 2003, 242: 119-135. 10.1007/s00606-003-0056-4.View ArticleGoogle Scholar
- Hupfer H, Swaitek M, Hornung S, Herrmann RG, Maier RM, Chiu WL, Sears B: Complete nucleotide sequence of the Oenothera elata plastid chromosome, representing plastome 1 of the five distinguishable Euoenthera plastomes. Mol Gen Genet. 2000, 263: 581-585.PubMedGoogle Scholar
- Spielmann A, Roux E, von Allmen J, Stutz E: The soybean chloroplast genome: completed sequence of the rps19 gene, including flanking parts containing exon 2 of rpl2 (upstream), but lacking rpl22 (downstream). Nucl Acids Res. 1988, 16: 1199-PubMed CentralView ArticlePubMedGoogle Scholar
- Milligan BG, Hampton JN, Palmer JD: Dispersed repeats and structural reorganization in subclover chloroplast DNA. Mol Biol Evol. 1989, 6: 355-368.PubMedGoogle Scholar
- Gantt JS, Baldauf SL, Calie PJ, Weeden NF, Palmer JD: Transfer of rpl22 to the nucleus greatly preceded its loss from the chloroplast and involved the gain of an intron. EMBO J. 1991, 10: 3073-3078.PubMed CentralPubMedGoogle Scholar
- Doyle JJ, Doyle JL, Palmer JD: Multiple independent losses of two genes and one intron from legume chloroplast genomes. Syst Bot. 1995, 20: 272-294.View ArticleGoogle Scholar
- Saski C, Lee S, Daniell H, Wood T, Tomkins J, Kim H-G, Jansen RK: Complete chloroplast genome sequence of Glycine max and comparative analyses with other legume genomes. Plt Mol Biol. 2005, 59: 309-322. 10.1007/s11103-005-8882-0.View ArticleGoogle Scholar
- Millen RS, Olmstead RG, Adams KL, Palmer JD, Lao NT, Heggie L, Kavanagh TA, Hibberd JM, Gray JC, Morden CW, Calie PJ, Jermlin LS, Wolfe KH: Many parallel losses of infA from chloroplast DNA during angiosperm evolution with multiple independent transfers to the nucleus. The Plant Cell. 2001, 13: 645-658. 10.1105/tpc.13.3.645.PubMed CentralView ArticlePubMedGoogle Scholar
- Cosner ME, Jansen RK, Palmer JD, Downie SR: The highly rearranged chloroplast genome of Trachelium caeruleum (Campanulaceae): multiple inversions, inverted repeat expansion and contraction, transposition, insertions/deletions, and several repeat families. Curr Genet. 1997, 31: 419-429. 10.1007/s002940050225.View ArticlePubMedGoogle Scholar
- Downie SR, Palmer JD: Use of chloroplast DNA rearrangements in reconstructing plant phylogeny. Molecular Systematics of Plants. Edited by: Soltis PS, Soltis DE, Doyle JJ. 1992, New York: Chapman and Hall, 14-35.View ArticleGoogle Scholar
- Maul JE, Lilly JW, Cui L, dePamphilis CW, Miller W, Harris EH, Stern DB: The Chlamydomonas reinhardtii plastid chromosome: islands of genes in a sea of repeats. The Plt Cell. 2002, 14: 1-22. 10.1105/tpc.140110.View ArticleGoogle Scholar
- Pombert J-F, Otis C, Lemieux C, Turmel M: The chloroplast genome sequence of the green alga Pseudendoclonium akinetum (Ulvophyceae) reveals unusual structural features and new insights into the branching order of chlorophyte lineages. Mol Biol Evol. 2005, 22: 1903-1918. 10.1093/molbev/msi182.View ArticlePubMedGoogle Scholar
- Lee S-B, Kaittanis C, Jansen RK, Hostetler JB, Tallon LJ, Town CD, Daniell H: The complete chloroplast genome sequence of Gossypium hirsutum: organization and phylogenetic relationships to other angiosperms. BMC Genomics. 2006, 7: 61-10.1186/1471-2164-7-61.PubMed CentralView ArticlePubMedGoogle Scholar
- Palmer JD: Plastid chromosomes: structure and evolution. The Molecular Biology of Plastids. Edited by: Bogorad L, Vasil K. 1991, San Diego: Academic Press, 5-53.View ArticleGoogle Scholar
- Schmitz-Linneweber C, Regel R, Du TG, Hupfer H, Herrmann RG, Maier RM: The plastid chromosome of Atropa belladonna and its comparison with that of Nicotiana tabacum: the role of RNA editing in generating divergence in the process of plant speciation. Mol Biol Evol. 2002, 19: 1602-1612.View ArticlePubMedGoogle Scholar
- Hirose T, Kusumegi T, Tsudzuki T, Sugiura M: RNA editing sites in tobacco chloroplast transcripts: editing as a possible regulator of chloroplast RNA polymerase activity. Mol Gen Genet. 1999, 262: 462-467. 10.1007/s004380051106.View ArticlePubMedGoogle Scholar
- Wolf PG, Rowe CA, Hasebe M: High levels of RNA editing in a vascular plant chloroplast genome: analysis of transcripts from the fern Adiantum capillus-veneris. Gene. 2004, 339: 89-97. 10.1016/j.gene.2004.06.018.View ArticlePubMedGoogle Scholar
- Fiebig A, Stegemann S, Bock R: Rapid evolution of RNA editing sites in a small non-essential plastid gene. Nucl Acids Res. 2004, 32: 3615-3622. 10.1093/nar/gkh695.PubMed CentralView ArticlePubMedGoogle Scholar
- Monde RA, Schusterc G, Stern DB: Processing and degradation of chloroplast mRNA. Biochimie. 2000, 82: 573-582. 10.1016/S0300-9084(00)00606-4.View ArticlePubMedGoogle Scholar
- Rochaix JD: Posttranscriptional control of chloroplast gene expression. From RNA to photosynthetic complex. Plt Phys. 2001, 125: 142-144. 10.1104/pp.125.1.142.View ArticleGoogle Scholar
- Donoghue MJ, Doyle JA: Phylogenetic studies of seed plants and angiosperms based on morphological characters. The Hierarchy of Life: Molecules and Morphology in Phylogenetic Studies. Edited by: Bremer K, Jornvall H. 1989, Amersterdam: Elsevier Science Publishers, 181-193.Google Scholar
- Doyle JA, Hotton CL: Diversification of early angiosperm pollen in a cladistic context. Pollen and Spores: Patterns of Diversification. Edited by: Blackmore S, Barnes SH. 1991, Oxford: Clarendon, 169-195.Google Scholar
- Soltis DE, Soltis PS, Nickrent DL, Johnson LA, Hahn WJ, Hoot SB, Sweere JA, Kuzoff RK, Kron KA, Chase M, Swensen SM, Zimmer E, Shaw SM, Gillespie LJ, Kress WJ, Sytsma K: Angiosperm phylogeny inferred from 18S ribosomal DNA sequences. Ann Missouri Bot Gard. 1997, 84: 1-49. 10.2307/2399952.View ArticleGoogle Scholar
- Hoot SB, Magallón S, Crane PR: Phylogeny of basal eudicots based on three molecular data sets: atpB, rbcL, and 18S nuclear ribosomal DNA sequences. Ann Missouri Bot Gard. 1999, 86: 1-32. 10.2307/2666215.View ArticleGoogle Scholar
- Soltis PS, Soltis DE, Chase MW: Angiosperm phylogeny inferred from multiple genes as a tool for comparative biology. Nature. 1999, 402: 402-40. 10.1038/46528.View ArticlePubMedGoogle Scholar
- Rokas A, Williams BL, King N, Carrol SB: Genome-scale approaches to resolving incongruence in molecular phylogenies. Nature. 2003, 425: 798-804. 10.1038/nature02053.View ArticlePubMedGoogle Scholar
- Swofford DL, Olsen GJ, Waddell PJ, Hillis DM: Phylogenetic inference. Molecular Systematics. Edited by: Hillis DM, Moritz C, Mable BK. 1996, Massachusetts: Sinauer Associates Inc, 407-514.Google Scholar
- Bruno WJ, Halpern AL: Topological bias and inconsistency of maximum likelihood using wrong models. Mol Biol Evol. 1999, 16: 564-566.View ArticlePubMedGoogle Scholar
- Swofford DL, Waddell PJ, Huelsenbeck JP, Foster PG, Lewis PO, Rogers JS: Bias in phylogenetic estimation and its relevance to the choice between parsimony and likelihood methods. Syst Biol. 2001, 50: 525-539. 10.1080/106351501750435086.View ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Ronquist F, Nielson R, Bolback JP: Bayesian inference of phylogeny and its impact on evolutionary biology. Science. 2001, 294: 2310-2314. 10.1126/science.1065889.View ArticlePubMedGoogle Scholar
- Jansen RK, Raubeson LA, Boore JL, dePamphilis CW, Chumley TW, Haberle RC, Wyman SK, Alverson AJ, Peery R, Herman SJ, Fourcade HM, Kuehl JV, McNeal JR, Leebens-Mack J, Cui L: Methods for obtaining and analyzing chloroplast genome sequences. Meth Enzym. 2005, 395: 348-384. 10.1016/S0076-6879(05)95020-9.View ArticlePubMedGoogle Scholar
- Clemson University BAC/EST Resource Center: [http://www.genome.clemson.edu]
- Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8: 186-194.View ArticlePubMedGoogle Scholar
- Wyman SK, Boore JL, Jansen RK: Automatic annotation of organellar genomes with DOGMA. Bioinformatics. 2004, 20: 3252-3255. 10.1093/bioinformatics/bth352. [http://bugmaster.jgi-psf.org/dogma]View ArticlePubMedGoogle Scholar
- Kurtz S, Choudhuri JV, Ohlebusch E, Schleiermacher C, Stoye J, Giegerich R: REPuter: the manifold applications of repeat analysis on a genomic scale. Nucl Acids Res. 2001, 29: 4633-4642. 10.1093/nar/29.22.4633.PubMed CentralView ArticlePubMedGoogle Scholar
- Higgins DG, Thompson JD, Gibson TJ: Using CLUSTAL for multiple sequence alignments. Meth Enzy. 1996, 266: 383-402.View ArticleGoogle Scholar
- Cui L, Veeraraghavan N, Richer A, Wall K, Jansen RK, Leebens-Mack J, Makalowska I, de Pamphillis CW: ChloroplastDB: the chloroplast genome database. Nucl Acids Res. 2006, 34: D692-D696. 10.1093/nar/gkj055. [http://chloroplast.cbio.psu.edu/]PubMed CentralView ArticlePubMedGoogle Scholar
- Swofford DL: PAUP*: Phylogenetic analysis using parsimony (*and other methods), ver. 4.0. 2003, Sunderland MA: Sinauer AssociatesGoogle Scholar
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.View ArticlePubMedGoogle Scholar
- Sullivan J, Abdo Z, Joyce P, Swofford DL: Evaluating the performance of a successive-approximations approach to parameter optimization in maximum-likelihood phylogeny estimation. Mol Biol Evol. 2005, 22: 1386-1392. 10.1093/molbev/msi129.View ArticlePubMedGoogle Scholar
- Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985, 39: 783-791.View ArticleGoogle Scholar
- Wakasugi T, Tsudzuki J, Ito S, Nakashima K, Tsudzuki T, Sugiura M: Loss of all ndh genes as determined by sequencing the entire chloroplast genome of the black pine Pinus thunbergii. Proc Natl Acad Sci USA. 1994, 91: 9794-9798.PubMed CentralView ArticlePubMedGoogle Scholar
- Hiratsuka J, Shimada H, Whittier R, Ishibashi T, Sakamoto M, Mori M, Kondo C, Honji Y, Sun CR, Meng BY, Li YQ, Kanno A, Nishizawa Y, Hirai A, Shinozaki K, Sugiura M: The complete sequence of the rice (Oryza sativa) chloroplast genome: intermolecular recombination between distinct tRNA genes accounts for a major plastid DNA inversion during the evolution of the cereals. Mol Gen Genet. 1989, 217: 185-194.View ArticlePubMedGoogle Scholar
- Asano T, Tsudzuki T, Takahashi S, Shimada H, Kadowaki K: Complete nucleotide sequence of the sugarcane (Saccharum officinarum) chloroplast genome: a comparative analysis of four monocot chloroplast genomes. DNA Res. 2004, 11: 93-99. 10.1093/dnares/11.2.93.View ArticlePubMedGoogle Scholar
- Maier RM, Neckermann K, Igloi GL, Kossel H: Complete sequence of the maize chloroplast genome: gene content, hotspots of divergence and fine tuning of genetic information by transcript editing. J Mol Biol. 1995, 251: 614-628. 10.1006/jmbi.1995.0460.View ArticlePubMedGoogle Scholar
- Sato S, Nakamura Y, Kaneko T, Asamizu E, Tabata S: Complete structure of the chloroplast genome of Arabidopsis thaliana. DNA Res. 1999, 6: 283-290. 10.1093/dnares/6.5.283.View ArticlePubMedGoogle Scholar
- Steane DA: Complete nucleotide sequence of the chloroplast genome from the Tasmanian blue gum, Eucalyptus globulus (Myrtaceae). DNA Res. 2005, 12: 215-220.View ArticlePubMedGoogle Scholar
- Shinozaki K, Ohme M, Tanaka M, Wakasugi T, Hayashida N, Matsubayashi T, Zaita N, Chunwongse J, Obokata J, Yamaguchi-Shinozaki K, Ohto C, Torazawa K, Meng BY, Sugita M, Deno H, Kamogashira T, Yamada K, Kusuda J, Takaiwa F, Kato A, Tohdoh N, Shimada H, Sugiura M: The complete nucleotide sequence of the tobacco chloroplast genome: its gene organization and expression. EMBO J. 1986, 5: 2043-2049.PubMed CentralPubMedGoogle Scholar
- Kim K-J, Lee H-L: Complete chloroplast genome sequence from Korean Ginseng (Panax schiseng Nees) and comparative analysis of sequence evolution among 17 vascular plants. DNA Res. 2004, 11: 247-261. 10.1093/dnares/11.4.247.View ArticlePubMedGoogle Scholar
- Daniell H, Lee S-B, Grevich J, Saksi C, Quesada-Vargas T, Guda C, Tomkins J, Jansen RK: Complete chloroplast genome sequences of Solanum bulbocastanum, Solanum lycopersicum and comparative analyses with other Solanaceae genomes. Theor Appl. 2006,Google Scholar
- Schmitz-Linneweber C, Maier RM, Alcaraz JP, Cottet A, Herrmann RG, Mache R: The plastid chromosome of spinach (Spinacia oleracea): complete nucleotide sequence and gene organization. Plt Mol Biol. 2001, 45: 307-315. 10.1023/A:1006478403810.View ArticleGoogle Scholar
- APG II 2002: An update of the Angiosperm Phylogeny Group classification for the orders and families of flowering plants APG II. Bot J Linn Soc. 2003, 141: 399-436. 10.1046/j.1095-8339.2003.t01-1-00158.x.View 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.