Skip to main content

The complete plastid genome sequence of Welwitschia mirabilis: an unusually compact plastome with accelerated divergence rates

Abstract

Background

Welwitschia mirabilis is the only extant member of the family Welwitschiaceae, one of three lineages of gnetophytes, an enigmatic group of gymnosperms variously allied with flowering plants or conifers. Limited sequence data and rapid divergence rates have precluded consensus on the evolutionary placement of gnetophytes based on molecular characters. Here we report on the first complete gnetophyte chloroplast genome sequence, from Welwitschia mirabilis, as well as analyses on divergence rates of protein-coding genes, comparisons of gene content and order, and phylogenetic implications.

Results

The chloroplast genome of Welwitschia mirabilis [GenBank: EU342371] is comprised of 119,726 base pairs and exhibits large and small single copy regions and two copies of the large inverted repeat (IR). Only 101 unique gene species are encoded. The Welwitschia plastome is the most compact photosynthetic land plant plastome sequenced to date; 66% of the sequence codes for product. The genome also exhibits a slightly expanded IR, a minimum of 9 inversions that modify gene order, and 19 genes that are lost or present as pseudogenes. Phylogenetic analyses, including one representative of each extant seed plant lineage and based on 57 concatenated protein-coding sequences, place Welwitschia at the base of all seed plants (distance, maximum parsimony) or as the sister to Pinus (the only conifer representative) in a monophyletic gymnosperm clade (maximum likelihood, bayesian). Relative rate tests on these gene sequences show the Welwitschia sequences to be evolving at faster rates than other seed plants. For these genes individually, a comparison of average pairwise distances indicates that relative divergence in Welwitschia ranges from amounts about equal to other seed plants to amounts almost three times greater than the average for non-gnetophyte seed plants.

Conclusion

Although the basic organization of the Welwitschia plastome is typical, its compactness, gene content and high nucleotide divergence rates are atypical. The current lack of additional conifer plastome sequences precludes any discrimination between the gnetifer and gnepine hypotheses of seed plant relationships. However, both phylogenetic analyses and shared genome features identified here are consistent with either of the hypotheses that link gnetophytes with conifers, but are inconsistent with the anthophyte hypothesis.

Background

Welwitschia mirabilis Hook f. (Welwitschiaceae) is a morphologically unique gymnosperm ("without parallel among all living vascular plants" [1]) of the Namib Desert of southwestern Africa (Namibia and Angola) and is the only extant member of its genus and family. The species is dioecious, and each adult plant consists of a giant taproot, a very short woody stem, and two permanent 'strap-shaped' leaves [1, 2]. Welwitschiaceae is one of three families in the Gnetophyta (each family would be placed in its own order and class when the group is recognized at the level of phylum), the other two being Ephedraceae, comprised of the genus Ephedra, and Gnetaceae, comprised of Gnetum. Welwitschiaceae diverged from other gnetophytes prior to the early Cretaceous [3].

Gnetophytes are of intense interest not only due to their peculiar morphology but also because controversies regarding seed plant phylogeny revolve around their placement (reviewed in [4, 5]). In morphologically based cladistic analyses [e.g., [6, 7]], gnetophytes are the extant sister to angiosperms (the "anthophyte" hypothesis). However the anthophyte hypothesis has rarely been recovered in analyses based on molecular data [4, 5, 8]. In recent work, multiple outcomes are often supported, even within the same paper, depending on which genes, which sites and which methods are employed [e.g., [5, 9, 10]]. Gnetophytes have been placed sister to all other seed plants [e.g., [11]], sister to conifers (the "gnetifer" hypothesis) [e.g., [12]], or nesting within conifers as the sister to Pinaceae (the "gne-pine" hypothesis) [e.g., [13–15]].

Placement of the gnetophytes based on molecular data has been problematic, apparently due to accelerated rates of evolution in the lineage, which can lead to Long Branch Attraction or LBA [16, 17]. Increased taxon sampling to break up long branches is commonly regarded as an effective approach for overcoming LBA, as is adding more sequence data, selecting slower markers, selecting slower positions and representing lineages with slowly evolving exemplars [18]. Furthermore, using genomic characters, e.g., inversions and gene or intron losses, which are less vulnerable to LBA, can also be helpful [19, 20]. In plants, the chloroplast genome is the primary target when attempting to generate large amounts of sequence data and genomic characters for phylogenetics, but unfortunately gymnosperms are very poorly represented among currently available chloroplast genome sequences. This paucity of gymnosperm chloroplast genome sequences (only three gymnosperm plastomes are available, and two are species from the same genus: Pinus thunbergii [21], P. koraiensis, and Cycas taitungensis [22]) currently limits the ability to construct genome level data matrices for seed plant phylogenetics.

The typical seed plant plastid genome [23] contains two copies of a large inverted repeat (IR) separated by large (LSC) and small (SSC) single copy regions. The genome is usually comprised of 150,000 to 160,000 base pairs (bp), includes approximately 120 different genes and is highly conserved in both gene order and content [24–26]. This general form (excepting some minor variation in IR boundaries and some differences in gene content between gymnosperms and angiosperms) is found in plastomes of Cycas [22], early-diverging angiosperms (such as Amborella [27] and Nuphar [28]), magnollids [29, 30] and various eudicots [31].

In contrast, some seed plant lineages contain plastid genomes that vary from this typical form [25]. Some plant groups contain genomes that lack one copy of the large IR [32, 33] or have greatly expanded IR regions [34]. Gene loss can occur, and the most extreme examples can be found in the plastid genomes of parasitic angiosperms [35]. Rearrangements have affected the gene order of plastid genomes, in some lineages slightly [e.g., [36–38]] but in a few cases greatly [34, 39–41]. And, although the rate of evolutionary change of plastid genes is largely conservative [42, 43], it is elevated in some lineages [44, 45]. Here, in describing the completely sequenced Welwitschia mirabilis plastome, we add to these examples of atypical plastid genomes and discuss phylogenetic implications.

Results

General Characteristics of the Genome

The 119,726 bp Welwitschia mirabilis plastid genome [GenBank: EU342371], similar to other chloroplast genomes, is A+T rich overall and in all compartments except for the RNA genes [see Additional file 1 and Additional file 2 for details]. The plastome consists of a large single copy region (LSC) of 68,556 base pairs (bp) and a small single copy region (SSC) of 11,156 bp, separated by two copies of the large inverted repeat (IRa and IRb) of 20,007 bp each (Fig 1). This quadripartite structure is typical among most land plant and some algal chloroplast genomes [46, 47]. The small size of the genome is unexpected; Welwitschia possesses the smallest plastid genome of any published non-parasitic land plant that still contains the large IR. The genome is similar in size to two other publicly available chloroplast genomes from gymnosperms, Pinus koraiensis (116,866 bp) and Pinus thunbergii (119,707). However, both Pinus plastomes have lost all but 475 or 495 bp (respectively) of their IR [48]. The Welwitschia plastome is less than 75% the size of the plastid genome reported for Cycas taitungensis [22].

Figure 1
figure 1

Gene map of the Welwitschia mirabilis plastome. Genes shown on the inside of the circle are transcribed counterclockwise and those on the outside clockwise. Gene boxes are color coded by functional group as shown in the key in the figure. The genome shows a structure typical to most chloroplast genomes: two copies of the inverted repeat region separating the large and small single copy regions. There are, however, multiple gene-order changes and gene losses relative to the ancestral genome organization, as well as slight expansions of the IR at the LSC and SSC boundaries.

The Welwitschia plastome is also unusually compact (i.e., a higher proportion of the genome is comprised of coding sequence and less of non-coding sequence). Coding regions constitute 66% of the Welwitschia plastome, making this the most compact of any non-parasitic, land plant chloroplast genome sequenced so far. Pinus thunbergii is the most similar, with 61% coding (Table 1). However, the NCBI value for Pinus thunbergii (and perhaps some of the other genomes) includes non-conserved ORFs and pseudogenes as coding sequence, which we excluded from coding in our calculations for Welwitschia. The NCBI coding percentage values are automatically generated based on the submitted annotation. Therefore, different methods of annotation (e.g., whether or not all ORFs are annotated) lead to ambiguities. For example, we recalculated compactness for the Pinus thunbergii plastid genome, using the same approach we applied to the Welwitschia genome, and found P. thunbergii to be 57.3% coding (rather than the 61% reported at NCBI). Thus our conclusions are conservative, and in no instance would the values reported in Table 1 be increased if calculated using the method applied to Welwitschia. Therefore, the compactness of our genome may be even more remarkable when taking into account that the genomes in the NCBI database may be less compact, by our methods of calculation, than reported. Even so, overall, the average coding percentage of the available land plant genomes, using NCBI values and excluding Welwitschia, is only 50.0% (Table 1), compared to 66% in Welwitschia.

Table 1 Compactness measures for vascular plant plastid genomes.

The LSC boundary of the IR in Welwitschia is located in the 3' end of psbA. The LSC-end of the IR includes the genes ycf2-trnH-trnI-3'psbA, with ycf2, trnH-GUG, and psbA on one strand and trnI-CAU on the other. Two events are required to explain the expanded IR seen in Welwitschia, from the ancestral seed plant condition (i.e., ycf2 and trnH duplicated in the IR ancestrally) proposed by Wu et al [22]. First, trnI-CAU was duplicated into the IR via an expansion of IRb at the JLB boundary. Second, IRa expanded at the JLA boundary to include a portion of psbA. Interestingly, the remnant IR found in Pinus thunbergii, which includes trnI and 3'psbA, matches the LSC end of the IR in Welwitschia (almost exactly: in the Welwitschia plastome only 77 bp more of psbA is duplicated). The duplication of trnI in this context and followed by a partial duplication of psbA is found in no other completely sequenced chloroplast genome, and other methods (targeted PCR and sequencing) show this "motif" to exist only in plastomes of gnetophytes and conifers [22, 49]. Presumably the loss of the IR in pines (and perhaps all conifers) is a further modification of a Welwitschia-like IR.

The gene order in the Welwitschia chloroplast genome is rearranged compared to more "typical" seed plant plastomes such as Cycas, Amborella, or Nicotiana. There are 14 locations (i.e., breakpoints) where gene adjacencies differ in Welwitschia compared to Cycas taitungensis, excluding differences due to gene losses and IR boundary shifts. A minimum of nine inversions (i.e., reversals on the chromosomal scale that change the order and orientation of one or more genes) would be required to convert the Cycas (ancestral) gene order to that of Welwitschia (derived), in addition to gene losses and IR extent changes. No gene order changes occur in the IR, but 4 inversions (7 breakpoints) are proposed for the LSC and 5 inversions (7 breakpoints) for the SSC (Fig 2). Although both Pinus and Welwitschia plastomes have undergone inversions, they share none in common. Over the entire plastome, there is only one clear instance where both genomes are disrupted at the same general location (trnT-GGU and trnE-UUC, immediately adjacent in unrearranged land plant plastomes, have been moved apart in both genomes). However, a shared inversion would require two shared points of disruption. There is a second possible shared disruption in the region of ndhF. However, even if there is truly a shared point of disruption (hard to define due to gene losses), these two endpoints (trnT-trnE and ndhF) could not have been used together in a single inversion; if these two breakpoints flanked one inversion in Welwitschia (with one endpoint in the LSC and the other in the SSC) the IR would have become a direct rather than inverted repeat. Since the repeat is still inverted we can be confident that such an inversion did not take place.

Figure 2
figure 2

Rearrangement scenario showing one possible explanation of differences observed in the Welwitschia SSC region. Gene order and orientation are represented in the diagram, but genes and IGS are not shown to scale. Map 1 shows the SSC and flanking ends of the IR as the region appears in Cycas or Ginkgo, maps 2–6 are hypothetical, and map 7 illustrates the region as it appears in the Welwitschia genome. Event a is an inversion that reverses the orientation of trnN. One endpoint of this inversion may have disrupted ndhF, leading to its loss. Event b is an inversion that reverses chlN. One endpoint may have disrupted chlL. Event c is the copy correction of the second IR copy to reflect the gene order change (event a) in the other copy. Events d, e, and f are inversions modifying the order and orientation of blocks of genes within the SSC. Again, in step e, inversion breakpoints may have disrupted genes. Events g, h, and i are additional gene losses (chlN is still detectable as a pseudogene) not directly related to inversion breakpoints. Event j is an extension of the inverted repeat into the SSC to copy rpl32 into the IR. The positions of the inversion endpoints are defined by the gene adjacencies in Welwitschia as compared to the ancestral condition, however exactly how those endpoints are combined into inversion events is speculative. Thus the endpoints indicated for events b-f may have been combined in different ways and the events may have occurred in an alternative order than that represented in this model.

Gene Content

The Welwitschia genome contains 101 distinct, presumably functional, genes: 31 unique tRNA genes, four rRNA gene species, and 66 different protein-coding genes (including 5 widely conserved ORFs or ycf genes). Four rRNA genes, eight tRNA genes and four protein-coding genes are fully or in part duplicated in the IR. Of the 18 genes usually found to contain introns in land plant plastid genomes, only 12 are still present in the Welwitschia plastome and two, petD and clpP, lack introns. In addition to genes duplicated in the IR, Welwitschia also has a duplication of the trnR-CCG gene, as well as a partial duplication of trnF-GAA. Of the genes normally present in land plant chloroplast genomes, four (Fig 3) are pseudogenes (detectable but truncated and containing various frameshift mutations resulting in numerous premature stop codons) and 15 are completely absent (at least we were unable to detect any remnants of them in the genome).

Figure 3
figure 3

Harr plots comparing sequence from the Welwitschia plastome with homologous regions from Podocarpus and Ginkgo . The extent of genes found in each region of the Podocarpus and Ginkgo plastomes are shown at the top of each graph as colored boxes. The Welwitschia (WEMI) sequence is represented as the Y-axis. Diagonal lines indicate regions of similarity between the two sequences being compared. Most portions of these genes found in Podocarpus and Ginkgo lack equivalent sequence in Welwitschia and the sections exhibiting similarity are disjointed.

The entire complement of ndh genes (11 subunits encoding NADH dehydrogenase) is non-functional, with 10 being absent and one, ndhB, being a pseudogene. The ndhB remnant is a small section of the 5' exon (Fig 3), only 121 bp long, whereas the fully functional Nicotiana tabacum ndhB gene is 1,530 bp long. The loss of the ndh genes also has been reported in Pinus [21], where more of the genes (seven) remain detectable as pseudogenes. In Welwitschia elevated rates and selection for compactness could accelerate the loss of detectable gene remnants following the initial loss of function. Presumably the remnant ndhB gene is still recognizable because it resides in the IR where rates are reduced [50]. Within completely sequenced land plant plastomes, the ndh genes have also been lost from the chloroplast genome of the orchid Phalaenopsis [45] and the non-photosynthetic parasite Epifagus [51]. The genes are found in some but not all green algal plastid genomes, and are lacking from plastomes of the red algal lineage and Cyanophora [52].

The genes for the three chl subunits, encoding the enzyme protochlorophyllide reductase, are also missing or non-functional, with one gene (chlL) being completely absent and the other two (chlN and chlB) being pseudogenes. The chlN and chlB pseudogenes are truncated, to 697 bp and 604 bp respectively, from a full length (as represented by Pinus thunbergii) of 1,530 bp and 1,401 bp. Based on earlier work [53], we expected chlL to be missing or highly divergent in Welwitschia. However results of that previous study, using Southern hybridization, indicated that chlL is present in Ephedra and Gnetum, conflicting with the results of Wu et al [22] for Gnetum. An examination of the Welwitschia plastome suggests that the loss of the chl genes may have been initiated by an inversion. One breakpoint may have been located within the chlL gene itself (Fig 2), so that the inversion split the gene apart, disrupting it and causing it to become non-functional. Of course, it is also possible that the gene was inactivated by simple base substitution. Once the chlL (or other) subunit was inactivated, no functional enzyme could be produced, eliminating any selection to retain the other two subunits as intact genes. It seems suggestive that the chlL gene is completely missing whereas chlN and chlB remain as pseudogenes, consistent with the hypothesis that chlL was lost first. These three genes are also missing from all angiosperm plastomes [31] as well as from the Psilotum chloroplast genome [53]. The product of these chl genes encodes an enzyme that allows chlorophyll to "green in the dark". An alternative, nuclear-encoded protein that requires light can also affect chlorophyll maturation. Thus the loss of these genes presumably would not be lethal even if the plastid genes were not successfully transferred to the nuclear genome.

Beyond the loss of the ndh and chl gene families, we were also unable to detect the genes accD, psaM, rpl23, rps16 or trnP-GGG in the Welwitschia plastid genome. The gene accD has been lost independently in numerous prokaryote and eukaryote lineages [41]. The gene rpl23 has been reported missing from the plastome of the angiosperms Spinacia [54, 55] and Trachelium [31]. The gene rps16 has experienced numerous independent losses in land plants [31, 56, 57]. Overall, the gene content of Welwitschia appears to be very close to what Wu et al [22] report for Gnetum, except that rps15, which is reported as absent in Gnetum, is present in Welwitschia.

Divergence for Protein-Coding Genes

In most cases Welwitschia plastid genes are more divergent than the genes of other seed plants. We attempted to measure relative divergence with a ratio derived from a simple comparison of pairwise distances, an extension of the approach used by Hajibabaei, Xia, and Drouin [58]. For each of 57 genes (Table 2), we calculated the average pairwise distance among non-gnetophyte seed plants and the average pairwise distance between Welwitschia and each representative non-gnetophyte seed plant. We then determined the ratio of these two averages to determine a "Relative Divergence Factor" for each gene. If the Relative Divergence Factor is less than one, then that gene in the Welwitschia plastome shows less divergence than average, if it is one, the gene is equally divergent, and, if greater than one, the Welwitschia gene exhibits above average divergence.

Table 2 Comparison of Relative Divergence Factors calculated using different reference sets of non-gnetophyte seed plants.

We selected from among 10 taxa to calculate the non-gnetophyte average: Ginkgo, Cycas, Pinus, Podocarpus, Amborella, Nuphar, Nymphaea, Calycanthus, Ranunculus, and Acorus. We used A) all 10 taxa, as well as B) Ginkgo and Cycas plus Ranunculus and Pinus, C) Ginkgo and Cycas plus Amborella and Pinus, D) Ginkgo and Cycas plus Amborella and Podocarpus, and E) Ginkgo and Cycas plus Amborella, Pinus and Podocarpus. Phylogenetic structure within the non-gnetophyte taxa could confound this approach. For example, in Set A, which contains multiple angiosperms, Nuphar and Nymphaea, especially, are rather closely related. Having such "small" distances included in the non-gnetophyte average could bias the result. However, we did not obtain very different results when using single angiosperm exemplars. In terms of the higher-level structure in the non-gnetophytes, internal branches are probably so short (Fig 4) relative to the terminals that this structure would have minimal impact. Also, if the true phylogenetic position of Welwitschia is anywhere within the other seed plants, then the bias should be against the Welwitschia to non-gnetophyte distances. In summary, this seems a reasonable enough way to investigate which genes are more or less divergent.

Figure 4
figure 4

Phylogenetic hypotheses obtained using different optimization criteria based on 57 protein-coding plastid genes. Values associated with branches indicate the level of bootstrap support or, for the Bayesian analysis, posterior probability. Branch lengths are shown proportional to inferred amount of change in all trees.

Across all the comparisons and all genes, the Relative Divergence Factor ranged from 0.87 to 3.45 (Table 2). [See Additional files 3, 4, 5, 6, 7 for details of the calculations for each reference set.] For about 20% of the comparisons, the difference between the non-gnetophyte average and the Welwitschia to non-gnetophtye average was not significant, either because the averages were basically equal or, more often, because the variation around the average was great enough to make them statistically indistinguishable. About 25% of the genes had an average Relative Divergence Factor above 2. When the Relative Divergence Factor was determined as an average of the five different calculations, only one gene (psaJ) had a divergence rate about equal to the non-gnetophytes (factor of 0.97), whereas 38 genes had Relative Divergence Factors above 1.5 and 14 above 2.

The average Relative Divergence Factor over all genes for each reference set ranged from 1.50 to 1.69 (all significant, p < 0.0001). The average Relative Divergence Factor from all genes, over all calculations, is 1.68. For all the genes together, we also applied Tajima's Relative Rate Test (Table 3). The Welwitschia genes have very significantly (p < 0.00001) elevated rates based on any of the comparisons conducted, whether the three taxon statement was compatible with commonly produced trees or not. In these tests, the proportion of unique sites identified in Welwitschia relative to the comparison taxa suggests that the Relative Divergence Factor method might be conservative.

Table 3 Tajima's Relative Rate Tests between Welwitschia and various other seed plants.

Phylogenetic Inference

Using the same 57 genes (47,858 aligned nt) as for our Relative Divergence Factor and Relative Rate comparisons, we inferred phylogenetic trees using various methods. Here, we also included genes from the pteridophytes Adiantum and Psilotum. The trees shown (Fig 4) include only Amborella as the representative angiosperm. However, analyses using other angiosperms (from those listed as "set A", above) as single exemplars or including all six angiosperm taxa yielded identical topologies (data not shown; however we will note that when the six angiosperms were included together, Amborella was strongly supported as the earliest diverging angiosperm by all methods). Neighbor-joining, Minimum Evolution (ME) and Maximum Parsimony (MP) trees yielded topologies with Welwitschia diverging first among the seed plants. In contrast, Maximum Likelihood (ML) and Bayesian (BI) analyses resulted in topologies that included a monophyletic gymnosperm clade and a Welwitschia-Pinus sister group. Bootstrap (BS) values and posterior probabilities (PP) indicate that each method, except ML, very strongly supports its particular solution. However, the SH test indicates that, for these data, the three topologies shown in Fig 4 do not have significantly different likelihoods. In contrast, the best ML topology (shown in the figure) is significantly better than the best anthophyte topology (p < 0.001).

It is interesting to note that although taxon sampling in our phylogenetic analyses is very limited, the trees obtained are mostly robust (e.g., high BS, PP values) and consistent with expectation. In other studies (e.g., [58] and [5] along with studies reviewed therein), based on total evidence from fewer genes but more taxa, distance methods and MP commonly place gnetophytes at the base of all seed plants, whereas ML and BI support a topology including gnetifers or gnepines. Here we see exactly the same pattern. By all methods, Welwitschia is inferred to represent a long branch, consistent with the relative rate tests and the relative divergence factor analyses. Since ME and MP are more likely to be confounded by rate heterogeneity, we would guess that the placement of Welwitschia in these trees is an artifact. However, when MP analyses were conducted on nucleotide data with the third position excluded or on deduced amino acid data, both of which should reduce the effects of rate heterogeneity, the same MP topology is recovered. As mentioned previously, according to SH test results, these data (for these taxa) cannot discriminate among the three topologies of Fig. 4. Although the addition of more sequence data allows us (with a sampling of only five seed plant taxa) to do no better than earlier work based on better taxon sampling, likewise our results are no worse. This suggests that a significant increase in number of genes (of a level available from whole plastome sequences) will allow a rather limited taxon sampling to provide important insights into seed plant phylogeny. Of course that sampling will need to be less limited that what is available to us and must include at least one conifer II representative.

Discussion

The morphology of Welwitschia makes it a very unusual organism. The plastid genome also has unusual features. Although the overall organization of the genome is typical of most land plants (i.e., having two single copy regions separated by two inverted repeat regions), the extent of the IR, gene content, gene order, rate of nucleotide divergence, and compactness, all are atypical. The phylogenetic significance of any of these atypical features can, at this point, be interpreted in only a limited manner. Until more gnetophyte plastome sequences are available we are uncertain whether features seen in Welwitschia are unique to that species or are characteristic of gnetophytes. Likewise, traits of the Pinus plastid genomes provide little detail about the nature of conifer genomes in general. Derived characters shared between Pinus and Welwitschia are compatible with either the gnetifer or the gne-pine hypotheses; however, they are inconsistent with the anthophyte hypothesis.

There are some differences in gene content that are shared between Welwitschia and angiosperms. Angiosperms and Welwitschia (but not Ephedra) lack the three chl genes. Angiosperms and Welwitschia (but not Gnetum) share the loss of trnP-GGG. These loss events could be shared due to common ancestry only if the gnetophytes were not monophyletic, in conflict with almost all recent work. Plastome gene loss has been a common pattern over the history of the plastid as an endosymbiont [52, 59]. The same genes are often lost independently in unrelated lineages [52]. Even within angiosperms, where gene content is largely stable in photosynthetic representatives, some genes have been lost in multiple instances (reviewed recently in [31]).

The most distinctive structural feature common to both Welwitschia and Pinus is the shared IR extent. The Pinus remnant IR (a 495 bp sequence containing trnI and a 3' portion of psbA [48]) matches the LSC end of the IR in Welwitschia and Gnetum [22] suggesting that the Pinus IR represents a reduction from a Welwitschia/Gnetum-like ancestor. Welwitschia and Pinus plastomes also both lack functional copies of all eleven ndh genes, the rps16 gene, and both introns of clpP. However, due to the reasons mentioned above, gene losses must be used cautiously as phylogenetic markers. Phylogenetic analyses, using ML and BI, based on 57 plastid genes also link Welwitschia and Pinus (as analyses based on 56 genes linked Gnetum and Pinus [22]). Of course, a representative of the second (non-Pinaceae) lineage of conifers (i.e., "conifer II" or Cupressophyta [60]) is necessary to distinguish between the gnetifer and gne-pine hypotheses.

Based on the Welwitschia plastome gene organization, we can speculate (due to placement of breakpoints) that inversions in the LSC or the SSC (Fig 2) may have initially destroyed some ndh genes leading to the loss of all the subunits, as was also a possible explanation in the case of the chlL gene loss. Inversion endpoints are located in the areas of the ndhCKJ cluster in the LSC, as well as ndhF and the ndhHAIGE cluster in the SSC. If any one of the ndh genes was disabled by an inversion (or by any other type of mutation) the remaining subunits would gradually decay to extinction. A gene-disrupting inversion in ndhF also may have initiated ndh gene loss in Pinaceae. Although an initial disruption of ndhF is a viable hypothesis to explain ndh loss in both Pinus and Welwitschia, the gene disrupting inversions could not be held in common between the gnetophyte and Pinaceae plastomes because, as discussed earlier, none of the inversions can be shared. Of course the initial gene disruptions could have been point mutations, not inversions; in which case, the losses of the ndh genes could possibly represent a synapomorphy supporting the gnepine clade (as ndh genes are reported missing from other gnetophytes [22, 61] and yet are commonly amplified from members of Cupressophyta [61]). Unfortunately gene losses lack complex characteristics to aid in determining homology.

We detected elevated levels of sequence divergence in most Welwitschia genes analyzed. Earlier work, based on limited numbers of genes [15, 62–64], has consistently found gnetophytes to have higher rates of sequence evolution in genes from each of the three compartments (plastid, mitochondrial and nuclear genomes) and in both ribosomal and protein-coding genes. Hajibabaei, Xia, and Drouin [58] showed, for nine genes (four plastid, three nuclear, and two mitochondrial), that the average pairwise distances between gnetophytes and non-gnetophyes were significantly higher than average pairwise distances among the non-gnetophytes. Here we expanded their approach to 57 plastid genes and found that almost all these genes exhibit above average divergence. We also found that, using the Relative Rate Test on the 57 genes, rates in Welwitschia were significantly higher in a variety of comparisons. Rates in Gnetum also were analyzed [22] using relative rate tests and, although transition rates in Gnetum were not always significantly different, transversion rates were. Thus these expanded studies (ours and Wu et al [22]) indicate that plastome-wide rate elevation has probably taken place in the gnetophtyes. Rate elevation is also suggested by the phylogenetic analyses, both in the estimated branch lengths and, presumably, in the conflict in outcomes among the different optimization methods.

The Welwitschia chloroplast genome was found to be unusually compact. The compactness of the genome can be interpreted in numerous ways, but might suggest that the small compact Welwitschia genome is the result of selective pressure to more rapidly replicate the genome by reducing intergenic space and by losing "non-essential" nucleotides. It is interesting to note that increased rate of replication is one factor hypothesized to lead to increased mutation rate [65]. The small genome size can be attributed to the high incidence of gene loss in addition to compactness. This suggests that a small, compact genome was more important to the success of Welwitschia (or, more likely, its ancestors) than the function of any genes lost.

Conclusion

Here we describe the first completely sequenced plastid genome of a gnetophyte. The Welwitschia plastome provides insight into the rates of sequence evolution of this highly divergent group of plants, as well as illustrating the possible gnetophyte pattern of gene loss and rearrangement. The most distinctive, potentially phylogenetically informative, feature of the Welwitschia chloroplast genome, the IR extent, supports a relationship of gnetophytes with conifers while being inconsistent with the anthophyte hypothesis. Phylogenetic hypotheses supported in analyses of 57 plastid genes, but minimal taxonomic representation, also lack any support for an anthophyte clade. The availability of the Welwitschia mirabilis plastome will provide important information for use in further phylogenetic studies resolving major questions about the evolution of seed plants. As plastome sequences for gymnosperms accumulate, genome level phylogenetic analyses should contribute to the resolution of controversies of seed plant phylogeny; currently however the number of gymnosperm plastomes is very limited. Additionally, further gnetophyte plastomes need to be sequenced in order to determine whether atypical characteristics seen in Welwitschia are shared by all gnetophytes, or are unusual genomic features of a very unique plant.

Methods

Chloroplast Extraction and DNA Sequencing

Chloroplasts were extracted from leaf tissue of Welwitschia mirabilis using the sucrose gradient method as described in Jansen et al [66]. Genomic RCA product (prepared using the Operon Repli-G Kit) was used to prepare the template for shotgun sequencing. Shotgun sequencing was performed at DOE Joint Genome Institute [67]. Assembly and its assessment were conducted using Phrap as implemented in Consed 15.1 [68]. Additional targeted sequencing was conducted on PCR products to attain quality in the finished sequence of a level of Q50 or higher for every nucleotide [69]. Each of the four IR-SC boundaries was also confirmed with independent PCR and sequencing reactions.

Annotation

The Welwitschia genome was annotated with the aid of DOGMA (Dual Organellar GenoMe Annotator) [70, 71]. Each region in which genes were not detected using DOGMA was investigated using TBLASTX [72] and ORFfinder [73]. In all cases gene boundaries (start and stop codons as well as intron/exon boundaries) were determined through comparison with other plastome annotations rather than via experimental evidence.

Coding Percentages

The coding percentage of the Welwitschia chloroplast genome sequence was calculated as the proportion of nucleotides that would be represented in mature gene products relative to the total number of nucleotides in the genome. Introns and pseudogenes were not included in the coding percentage calculation. The percent coding of all other genomes reported in Table 1 were obtained from the NCBI database.

Pseudogenes

Pseudogenes were detected and their extent determined using a variety of approaches. Any region that was found using DOGMA to contain fragmentary sequences of genes was aligned (using ClustalW) to functional versions of genes from various previously sequenced plants. In addition, Harr or dot plots were constructed using Pipmaker [74], under default stringency settings, to compare potentially homologous regions of the Welwitschia genome with gene sequences from the plastomes of Podocarpus macrophyllus and Ginkgo biloba (Raubeson et al, unpublished). The Harr plots were visually analyzed in order to determine presence of significant sequence similarity. Pseudogenes were annotated as the most extensive region of detectable similarity.

Gene Order and Loss

Gene order rearrangements were investigated manually and using GRIMM [75, 76]. When using GRIMM (which requires gene content to be identical), genomes were simplified such that IR expansions or contractions and gene losses were not taken into account. GRIMM reports the minimum number of inversions required to "convert" one genome into another and returns one model of inversions. The specific inversion scenario suggested by GRIMM is only one of several possible models of equal length that could explain the observed pattern. As GRIMM is unable to take into account IR expansion and contraction, visual analyses were also performed in order to take into account these events. These visual analyses were also used to hypothesize the loss of the ndh and chl genes.

Divergence

For 57 shared protein-coding genes (these are the 61 plastid genes first used for phylogenetic inference by Goremykin et al [27], minus four (clpP, petL, rpl2, and rpl32) that were problematic to align), pairwise distance analyses were performed initially between 10 taxa (Ginkgo, Cycas, Pinus, Podocarpus, Amborella, Nuphar, Nymphaea, Calycanthus, Ranunculus, and Acorus) and Welwitschia using MEGA3.1 [77]. The 10 taxa were chosen to represent all of the major clades of extant non-gnetophtye seed plants and constituted reference Set A. Distances were calculated using both the LogDet (data not shown) and the Kimura two-parameter models, which gave very similar results. Default MEGA3.1 parameters were used in all calculations. Standard error for each distance was calculated via bootstrapping with 100 replicates. The Relative Divergence Factor was calculated by dividing the average distance calculated from Welwitschia to each non-gnetophytes by the average distance from comparisons among the 10 non-gnetophytes. We repeated these calculations using Kimura two-parameter distance and different subsets of the ten taxa as reference sets. For each of the five sets of calculations, we determined whether or not the difference between the two average distances was significant using the two-tailed t-test. In addition, we preformed Tajima's Relative Rate test, using MEGA 4.0 [78], on the 57 concatenated genes.

Phylogenetic Inference

We added Welwitschia plastome genes to the alignment/data matrix of Leebens-Mack et al [79]. As mentioned above, four of the "standard" 61 genes were excluded from our analyses. The pteridophytes Psilotum and Adiantum were used to root the gymnosperm tree, although all analyses were run as and trees generated as "unrooted" topologies, i.e., we did not force Psilotum and Adiantum to any particular position on any tree. Bootstrap (BS) trees (100 heuristic replicates) were generated under Minimum Evolution (LogDet distances), Maximum Parsimony (MP), and Maximum Likelihood (ML), using a GTR + I + Γ model, optimizations using PAUP4.0b10 [80]. Individual MP analyses were conducted using branch and bound searches, always yielding a single most parsimonious tree identical in topology to the BS tree. Multiple individual ML heuristic searches were also conducted using estimation of all model parameters and obtaining starting trees by 10 replicates of random stepwise addition. Mr. Bayes 3.1 [81] was used, with default settings except as noted, to conduct a Bayesian analysis again using GTR + I + Γ (the model indicated by ModelTest [82] as the best fit to the data). Two analyses were run in parallel for 100,000 generations. After 6,000 generations, the standard deviation between the two hot chains was 0.000000 and remained at that level for the remainder of the run. Visual analysis using AWTY [83] also indicated that the chains had converged very early in the run. The first ten percent of the trees were discarded as "burn-in". PAUP4b10 was used to construct the consensus tree for the 901 trees retained. The Shimodaira-Hasegawa, or SH, Test [84] was implemented in PAUP4b10 using 1000 replicates and RELL approximation.

Abbreviations

Anthophyte:

hypothesized lineage composed of gnetophytes and angiosperms (and extinct taxa), BI: Bayesian Inference, bp: base pairs, BS: Booststrap, Conifer II: conifers are composed of two major clades, Pinaceae and the other families

the non-Pinaceae clade has been designated 'conifer II' or more recently Cupressophyta [60]:

Gnepine: hypothesized group composed of gnetophytes and Pinaceae, Gnetifer: hypothesized clade composed of gnetophytes as sister to Pinaceae + ConiferII, GRIMM: Genome rearrangement algorithms (genome rearrangements in mouse and man) [75], IR: Inverted Repeat (in plastid genomes usually a large, 10–25 kb, region that includes the rRNA genes), LBA: Long Branch Attraction (a phenomenon that confounds phylogenetic reconstruction as rapidly evolving lineages may share more characteristics due to chance rather than descent), LSC: Large Single Copy (a region of the plastid genome), ME: Minimum Evolution, ML: Maximum Likelihood, MP: Maximum Parsimony, ORF: Open Reading Frame, PP: Posterior Probability, RCA: Rolling Circle Amplification, SSC: Small Single Copy (a region of the plastid genome), WEMI: We lwitschia mi rabilis

References

  1. Gifford EM, Foster AS: Morphology and evolution of vascular plants. 1989, New York: W.H. Freeman and Company, 3

    Google Scholar 

  2. Chamberlain CJ: Gymnosperms: Structure and Function. 1935, Chicago: University of Chicago Press

    Google Scholar 

  3. Dilcher DL, Bernardes-De-Oliveira ME, Pons D, Lott TA: Welwitschiaceae from the Lower Cretaceous of Northeastern Brazil. Am J Bot. 2005, 92: 1294-1310. 10.3732/ajb.92.8.1294.

    Article  PubMed  Google Scholar 

  4. Magallon S, Sanderson MJ: Relationships among seed plants inferred from highly conserved genes: Sorting conflicting phylogenetic signals among ancient lineages. Am J Bot. 2002, 89: 1991-2006. 10.3732/ajb.89.12.1991.

    Article  CAS  PubMed  Google Scholar 

  5. Burleigh JG, Mathews S: Phylogenetic signal in nucleotide data from seed plants: Implications for resolving the seed plant tree of life. Am J Bot. 2004, 91: 1599-1613. 10.3732/ajb.91.10.1599.

    Article  CAS  PubMed  Google Scholar 

  6. Doyle JA, Donoghue MJ: Seed plant phylogeny and the origin of angiosperms: an experimental cladistic approach. Bot Rev. 1986, 52: 321-431. 10.1007/BF02861082.

    Article  Google Scholar 

  7. Rothwell GW, Serbet R: Lignophyte phylogeny and the evolution of spermatophytes: a numerical cladistic analysis. Syst Bot. 1994, 19: 443-482. 10.2307/2419767.

    Article  Google Scholar 

  8. Donoghue MJ, Doyle JA: Seed plant phylogeny: Demise of the anthophyte hypothesis?. Curr Biol. 2000, 10: R106-R109. 10.1016/S0960-9822(00)00304-3.

    Article  CAS  PubMed  Google Scholar 

  9. Rydin C, Kallersjo M, Friist EM: Seed plant relationships and the systematic position of Gnetales based on nuclear and chloroplast DNA: Conflicting data, rooting problems, and the monophyly of conifers. Int J Plant Sci. 2002, 163: 197-214. 10.1086/338321.

    Article  CAS  Google Scholar 

  10. Soltis DE, Soltis PS, Zanis MJ: Phylogeny of seed plants based on evidence from eight genes. Am J Bot. 2002, 89: 1670-1681. 10.3732/ajb.89.10.1670.

    Article  CAS  PubMed  Google Scholar 

  11. Rai HS, O'Brien HE, Reeves PA, Olmstead RG, Graham SW: Inference of higher-order relationships in the cycads from a large chloroplast data set. Mol Phylogenet Evol. 2003, 29: 350-359. 10.1016/S1055-7903(03)00131-3.

    Article  CAS  PubMed  Google Scholar 

  12. Chaw SM, Zharkikh A, Sung HM, Lau TC, Li WH: Molecular phylogeny of extant gymnosperms and seed plant evolution: Analysis of nuclear 18S rRNA sequences. Mol Biol Evol. 1997, 14: 56-68.

    Article  CAS  PubMed  Google Scholar 

  13. Qiu YL, Lee JH, Bernasconi-Quadroni F, Soltis DE, Soltis PS, Zanis M, Zimmer EA, Chen ZD, Savolainen V, Chase MW: The earliest angiosperms: evidence from mitochondrial, plastid and nuclear genomes. Nature. 1999, 402: 404-407. 10.1038/46536.

    Article  CAS  PubMed  Google Scholar 

  14. Chaw S, Parkinson CL, Cheng Y, Vincent TM, Palmer JD: Seed plant phylogeny inferred from all three plant genomes: Monophyly of extant gymnosperms and origin of Gnetales from conifers. Proc Natl Acad Sci USA. 2000, 97: 4086-4091. 10.1073/pnas.97.8.4086.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Bowe LM, Coat G, dePamphilis CW: Phylogeny of seed plants based on all three genomic compartments: Extant gymnosperms are monophyletic and Gnetales' closest relatives are conifers. Proc Natl Acad Sci USA. 2000, 97: 4092-4097. 10.1073/pnas.97.8.4092.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Felsenstein J: Cases in which parsimony or compatibility methods will be positively misleading. Syst Zool. 1978, 27: 401-410. 10.2307/2412923.

    Article  Google Scholar 

  17. Philippe H, Zhou Y, Brinkmann H, Rodrigue N, Delsuc F: Heterotachy and long-branch attraction in phylogenetics. BMC Evol Biol. 2005, 5: 50-10.1186/1471-2148-5-50.

    Article  PubMed Central  PubMed  Google Scholar 

  18. Moreira D, Phillippe H: Molecular phylogeny: pitfalls and progress. Internatl Microbiol. 2000, 3: 9-16.

    CAS  Google Scholar 

  19. Rokas A, Holland PWH: Rare genomic changes as a tool for phylogenetics. TREE. 2000, 15: 454-459.

    PubMed  Google Scholar 

  20. Boore JL: The use of genome-level characters for phylogenetic reconstruction. TREE. 2006, 21: 439-446.

    PubMed  Google Scholar 

  21. 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. 10.1073/pnas.91.21.9794.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Wu CS, Wang YN, Liu SM, Chaw SM: Chloroplast Genome (cpDNA) of Cycas taitungensis and 56 cp protein-coding genes of Gnetum parvifolium: insights into cpDNA evolution and phylogeny of extant seed plants. Mol Biol Evol. 2007, 24: 1366-1379. 10.1093/molbev/msm059.

    Article  CAS  PubMed  Google Scholar 

  23. Raubeson LA, Jansen RK: Chloroplast genomes of plants. Diversity and Evolution of Plants – Genotypic and phenotypic variation in higher plants. Edited by: Henry R. 2005, London: CABI Publishing, 45-68.

    Chapter  Google Scholar 

  24. Shinozaki K, Ohme M, Tanaka M, Wakasugi T, Hayashida N, Matsubayashi T, Zaita N, Chunwongse J, Obokata J, Yamaguchishinozaki 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 Central  CAS  PubMed  Google Scholar 

  25. 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.

    Chapter  Google Scholar 

  26. Palmer JD: Plastid chromosomes: structure and evolution. The molecular biology of plastids. Cell culture and somatic cell genetics of plants. 7A. Edited by: Hermann RG. 1991, Vienna:Springer-Verlag, 5-53.

    Google Scholar 

  27. 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.

    Article  CAS  PubMed  Google Scholar 

  28. Raubeson LA, Peery R, Chumley TW, Dziubek C, Fourcade HM, Boore JL, Jansen RK: Comparative chloroplast genomics: analyses including new sequences from the angiosperms Nuphar advena and Ranunculus macaranthus. BMC Genomics. 2007, 8: 174-10.1186/1471-2164-8-174.

    Article  PubMed Central  PubMed  Google Scholar 

  29. Cai ZQ, Penaflor C, Kuehl JV, Leebens-Mack J, Carlson JE, dePamphilis CW, Boore JL, Jansen RK: Complete plastid genome sequences of Drimys, Liriodendron, and Piper: implications for the phylogenetic relationships of magnoliids. BMC Evol Biol. 2006, 6: 77-10.1186/1471-2148-6-77.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Hansen DR, Dastidar SG, Cai Z, Penaflor C, Kuehl JV, Boore JL, Jansen RK: Phylogenetic and evolutionary implications of complete chloroplast genome sequences of four early diverging angiosperms: Buxus (Buxaceae), Chloranthus (Chloranthaceae), Dioscorea (Dioscoreaceae), and Illicium (Schisandraceae). Mol Phylog Evol. 2007, 45: 547-563. 10.1016/j.ympev.2007.06.004.

    Article  CAS  Google Scholar 

  31. Jansen RK, Cai Z, Raubeson LA, Daniell H, dePamphilis CW, Leebens-Mack J, Muller KF, Guisinger-Bellian M, Haberle RC, Hansen AK, Chumley TW, Lee S-B, Peery R, McNeal JR, Kuehl JV, Boore JL: Analysis of 81 genes from 64 chloroplast genomes resolves relationships in angiosperms and identifies genome-scale evolutionary patterns. Proc Natl Acad Sci USA. 2007, 104: 19369-19374. 10.1073/pnas.0709121104.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Strauss SH, Palmer JD, Howe GT, Doerksen AH: Chloroplast genomes of two conifers lack a large inverted repeat and are extensively rearranged. Proc Natl Acad Sci USA. 1988, 85: 3898-3902. 10.1073/pnas.85.11.3898.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Palmer JD, Osorio B, Aldrich J, Thompson WF: Chloroplast DNA evolution among legumes – loss of a large inverted repeat occurred prior to other sequence rearrangements. Curr Genet. 1987, 11: 275-286. 10.1007/BF00355401.

    Article  CAS  Google Scholar 

  34. Chumley TW, Palmer JD, Mower JP, Fourcade HM, Calie PJ, Boore JL, Jansen RK: The complete chloroplast genome sequence of Pelargonium x hortorum: Organization and evolution of the largest and most highly rearranged chloroplast genome of land plants. Mol Biol Evol. 2006, 23: 2175-2190. 10.1093/molbev/msl089.

    Article  CAS  PubMed  Google Scholar 

  35. Wolfe KH, Morden CW, Palmer JD: Function and evolution of a minimal plastid genome from a nonphotosynthetic parasitic plant. Proc Natl Acad Sci USA. 1992, 89: 10648-10652. 10.1073/pnas.89.22.10648.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Raubeson LA, Jansen RK: Chloroplast DNA evidence on the ancient evolutionary split in vascular land plants. Science. 1992, 255: 1697-1699. 10.1126/science.255.5052.1697.

    Article  CAS  PubMed  Google Scholar 

  37. Jansen RK, Palmer JD: A chloroplast DNA inversion marks an ancient evolutionary split in the sunflower family (Asteraceae). Proc Natl Acad Sci USA. 1987, 84: 5818-5822. 10.1073/pnas.84.16.5818.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Palmer JD, Osorio B, Thompson WF: Evolutionary significance of inversions in legume chloroplast DNAs. Curr Genet. 1988, 14: 65-74. 10.1007/BF00405856.

    Article  CAS  Google Scholar 

  39. 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.

    Article  CAS  PubMed  Google Scholar 

  40. Cosner ME, Raubeson LA, Jansen RK: Chloroplast DNA rearrangements in Campanulaceae: phylogenetic utility of highly rearranged genomes. BMC Evol Biol. 2004, 4: 27-10.1186/1471-2148-4-27.

    Article  PubMed Central  PubMed  Google Scholar 

  41. Lee HL, Jansen RK, Chumley TW, Kim KJ: Gene relocations within chloroplast genomes of Jasminum and Menodora (Oleaceae) are due to multiple, overlapping inversions. Mol Biol Evol. 2007, 24: 1161-1180. 10.1093/molbev/msm036.

    Article  CAS  PubMed  Google Scholar 

  42. Clegg MT, Gaut BS, Learn GH, Morton BR: Rates and patterns of chloroplast DNA evolution. Proc Natl Acad Sci USA. 1994, 91: 6795-6801. 10.1073/pnas.91.15.6795.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. Gaut BS, Muse SV, Clegg MT: Relative rates of nucleotide substitution in the chloroplast genome. Mol Phylogenet Evol. 1993, 2: 89-96. 10.1006/mpev.1993.1009.

    Article  CAS  PubMed  Google Scholar 

  44. Gaut BS, Morton BR, McCaig BC, Clegg MT: Substitution rate comparisons between grasses and palms: Synonymous rate differences at the nuclear gene Adh parallel rate differences at the plastid gene rbcL. Proc Natl Acad Sci USA. 1996, 93: 10274-10279. 10.1073/pnas.93.19.10274.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Chang CC, Lin HC, Lin IP, Chow TY, Chen HH, Chen WH, Cheng CH, Lin CY, Liu SM, Chang CC, Chaw SM: 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.

    Article  CAS  PubMed  Google Scholar 

  46. Raubeson LA, Jansen RK: Chloroplast genomes of plants. Plant Diversity and Evolution: Genotypic and Phenotypic Variation in Higher Plants. Edited by: Henry RJ. 2005, CAB International, 45-68.

    Chapter  Google Scholar 

  47. Turmel M, Otis C, Lemieux C: The complete chloroplast DNA sequence of the green alga Nephroselmis olivacea: Insights into the architecture of ancestral chloroplast genomes. Proc Natl Acad Sci USA. 1999, 96: 10248-10253. 10.1073/pnas.96.18.10248.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  48. Tsudzuki J, Nakashima K, Tsudzuki T, Hiratsuka J, Shibata M, Wakasugi T, Sugiura M: Chloroplast DNA of black pine retains a residual inverted repeat lacking rRNA genes: nucleotide sequences of trnQ, trnK, psbA, trnI and trnH and the absence of rps16. Mol Gen Genet. 1992, 232: 206-214.

    CAS  PubMed  Google Scholar 

  49. Raubeson LA: Chloroplast DNA structural similarities shared by conifers and Gnetales: coincidence or common ancestry? [abstract]. Am J Bot. 1998, 85: S153-10.2307/2446303.

    Article  Google Scholar 

  50. Perry AS, Wolfe KH: Nucleotide substitution rates in legume chloroplast DNA depend on the presence of the inverted repeat. J Mol Evol. 2002, 55: 501-508. 10.1007/s00239-002-2333-y.

    Article  CAS  PubMed  Google Scholar 

  51. Wolfe KH, Morden CW, Palmer JD: Function and evolution of a minimal plastid genome from a nonphotosynthetic parasitic plant. Proc Natl Acad Sci USA. 1992, 89: 10648-10652. 10.1073/pnas.89.22.10648.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  52. Martin W, Rujan T, Richly E, Hansen A, Cornelsen S, Lins T, Leister D, Stoebe B, Hasegawa M, Penny D: Evolutionary analysis of Arabidopsis, cyanobacterial, and chloroplast genomes reveals plastid phylogeny and thousands of cyanobacterial genes in the nucleus. Proc Natl Acad Sci USA. 2002, 99: 12246-12251. 10.1073/pnas.182432999.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  53. Burke DH, Raubeson LA, Alberti M, Hearst J, Jordon ET, Kirsch S, Conant DS, Valinski AEC, Stein DB: The chlL (frxC) Gene – Phylogenetic distribution in vascular plants and DNA sequence from Polystichum acrostichoides (Pteridophyta) and Synechococcus pp 7002 (Cyanobacteria). Plant Syst Evol. 1993, 187: 89-102. 10.1007/BF00994092.

    Article  CAS  Google Scholar 

  54. Thomas F, Massenet O, Dorne AM, Briat JF, Mache R: Expression of the rpl23, rpl2, and rps19 gene in spinach chloroplasts. Nucl Acids Res. 1988, 25: 2461-2472. 10.1093/nar/16.6.2461.

    Article  Google Scholar 

  55. 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.

    Article  CAS  Google Scholar 

  56. 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. 10.2307/2419496.

    Article  Google Scholar 

  57. Ohyama K: Chloroplast and mitochondrial genomes from a liverwort, Marchantia polymorpha – gene organization and molecular evolution. Bioscience Biotechnology and Biochemistry. 1996, 60: 16-24.

    Article  CAS  Google Scholar 

  58. Hajibabaei M, Xia J, Drouin G: Seed plant phylogeny: gnetophytes are derived conifers and a sister group to Pinaceae. Mol Phylogenet Evol. 2006, 40: 208-217. 10.1016/j.ympev.2006.03.006.

    Article  CAS  PubMed  Google Scholar 

  59. Martin W, Stoebe B, Goremykin V, Hansmann S, Hasegawa M, Kowallik KV: Gene transfer to the nucleus and the evolution of chloroplasts. Nature. 1998, 393: 162-165. 10.1038/30234.

    Article  CAS  PubMed  Google Scholar 

  60. Cantino PD, Doyle JA, Graham SW, Judd WS, Olmstead RG, Soltis DE, Soltis PS, Donoghue MJ: Towards a phylogenetic nomenclature of Tracheophyta. Taxon. 2007, 56: E1-E44.

    Article  Google Scholar 

  61. Graham SW, Olmstead RG: Utility of 17 chloroplast genes for inferring the phylogeny of the basal angiosperms. Am J Bot. 2000, 87: 1712-1730. 10.2307/2656749.

    Article  CAS  PubMed  Google Scholar 

  62. Chaw SM, Parkinson CL, Cheng YC, Vincent TM, Palmer JD: Seed plant phylogeny inferred from all three plant genomes: monophyly of extant gymnosperms and origin of Gnetales from conifers. Proc Natl Acad Sci USA. 2000, 97: 4086-4091. 10.1073/pnas.97.8.4086.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  63. Sanderson MJ, Wojciechowski MF, Hu JM, Khan TS, Brady SG: Error, bias, and long-branch attraction in data for two chloroplast photosystem genes in seed plants. Mol Biol Evol. 2000, 17: 782-797.

    Article  CAS  PubMed  Google Scholar 

  64. Sanderson MJ: Estimating absolute rates of molecular evolution and divergence times: a penalized likelihood approach. Mol Biol Evol. 2002, 19: 101-109.

    Article  CAS  PubMed  Google Scholar 

  65. Mindell DP, Thacker CE: Rates of Molecular Evolution: Phylogenetic Issues and Applications. Ann Rev Ecol Syst. 1996, 27: 279-303. 10.1146/annurev.ecolsys.27.1.279.

    Article  Google Scholar 

  66. 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 whole chloroplast genome sequences. Methods Enzymol. 2005, 395: 348-384.

    Article  CAS  PubMed  Google Scholar 

  67. JGI – Protocols in Production Sequencing. [http://www.jgi.doe.gov/sequencing/protocols/prots_production.html]

  68. Gordon D, Abajian C, Green P: Consed: A graphical tool for sequence finishing. Genome Res. 1998, 8: 195-202.

    Article  CAS  PubMed  Google Scholar 

  69. Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8: 186-194.

    Article  CAS  PubMed  Google Scholar 

  70. Wyman SK, Jansen RK, Boore JL: Automatic annotation of organellar genomes with DOGMA. Bioinformatics. 2004, 20: 3252-3255. 10.1093/bioinformatics/bth352.

    Article  CAS  PubMed  Google Scholar 

  71. DOGMA (Dual Organellar Genome Annotator). [http://dogma.ccbb.utexas.edu/]

  72. BLAST (Basic Local Alignment Tool). [http://www.ncbi.nlm.nih.gov/BLAST/]

  73. ORF Finder. [http://www.ncbi.nlm.nih.gov/gorf/gorf.html]

  74. Schwartz S, Zhang Z, Frazer KA, Smit A, Riemer C, Bouck J, Gibbs R, Hardison R, Miller W: PipMaker – A Web server for aligning two genomic DNA sequences. Genome Res. 2000, 10: 577-586. 10.1101/gr.10.4.577.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  75. Tesler G: GRIMM: genome rearrangements web server. Bioinformatics. 2002, 18: 492-493. 10.1093/bioinformatics/18.3.492.

    Article  CAS  PubMed  Google Scholar 

  76. GRIMM – Genome rearrangement algorithms. [http://nbcr.sdsc.edu/GRIMM/mgr.cgi]

  77. Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform. 2004, 5: 150-163. 10.1093/bib/5.2.150.

    Article  CAS  PubMed  Google Scholar 

  78. Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.

    Article  CAS  PubMed  Google Scholar 

  79. 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.

    Article  CAS  PubMed  Google Scholar 

  80. Swofford DL: PAUP*: Phylogenetic analysis using parsimony (* and other methods), ver. 4.0. 2003, Sunderland MA: Sinauer Associates

    Google Scholar 

  81. Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogeny. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.

    Article  CAS  PubMed  Google Scholar 

  82. Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.

    Article  CAS  PubMed  Google Scholar 

  83. Wilgenbusch JC, Warren DL, Swofford DL: AWTY: A system for graphical exploration of MCMC convergence in Bayesian phylogenetic inference. [http://ceb.csit.fsu.edu/awty]

  84. Shimodaira H, Hasegawa M: Multiple comparisons of log-likelihoods with applications to phylogenetic inference. Mol Biol Evol. 1999, 16: 1114-1116.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We thank Claude dePamphilis and Anthony Omeis for providing plant material from the Penn State Biology Greenhouse, Bob Jansen's lab for preparing the sequencing template, Josh Leigh for initial work on an earlier draft sequence, Rhiannon Peery for help with NCBI submission and use of Mulan, Kai Müller and Claude dePamphilis for help with the alignment, and Bob Jansen, Mary Guisinger and Tim Chumley for critically reading an earlier version of the manuscript. The comments of four anonymous reviewers also helped us improve the manuscript. Funding for this work was provided by the CWU Science Honors Research Program, the M.J. Murdock Charitable Trust, the CWU Faculty Development and Research Committee, and the National Science Foundation (DEB0075700 to LAR, DEB0120709 to LAR and JLB, and EF0629607 to LAR). This work was performed partly under the auspices of the US Department of Energy's Office of Science, Biological and Environmental Research Program, and by the University of California, Lawrence Berkeley National Laboratory under Contract No. DE-AC02-05CH11231.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Linda A Raubeson.

Additional information

Authors' contributions

SRM annotated the genome, performed most analyses, and wrote the initial draft as his undergraduate thesis; JVK and JLB generated the draft genome sequence and initial assembly; JLB also assisted with manuscript preparation; LAR conceived of the project, helped with finishing, generated the final assembly, performed phylogenetic work and relative rate tests, assisted with other analyses, and revised the thesis for submission.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

McCoy, S.R., Kuehl, J.V., Boore, J.L. et al. The complete plastid genome sequence of Welwitschia mirabilis: an unusually compact plastome with accelerated divergence rates. BMC Evol Biol 8, 130 (2008). https://doi.org/10.1186/1471-2148-8-130

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2148-8-130

Keywords