Isolation and structural analysis of chloroplast DNA

Characterization of the complete chloroplast genome of Brunfelsia brasiliensis (Spreng.) L.B.Sm. & Downs

Brunfelsia brasiliensis (B. brasiliensis) (Spreng.) L.B.Sm. & Downs is a perennial evergreen shrub that is widely cultivated as an ornamental plant in tropical and subtropical regions. This study presents the complete chloroplast genome of B. brasiliensis, having a length of 169,062 bp. The chloroplast genome contains 137 genes, including 92 protein-coding, 37 tRNA genes, and 8 rRNA genes. Phylogenetic analysis reveals that Petunia hybrida and P. exserta are closely linked to B. brasiliensis in evolutionary terms. This first chloroplast genome assembly is a valuable resource for future genetic and molecular biology studies of the genus Brunfelsia.

Complete Chloroplast Genome of Corethrodendron fruticosum (Papilionoideae: Fabaceae): Comparative and Phylogenetic Analysis

Corethrodendron fruticosum is an endemic forage grasses in China with high ecological value. In this study, the complete chloroplast genome of C. fruticosum was sequenced using Illumina paired-end sequencing. The C. fruticosum chloroplast genome was 123,100 bp and comprised 105 genes, including 74 protein-coding genes, 4 rRNA-coding genes, and 27 tRNA-coding genes. The genome had a GC content of 34.53%, with 50 repetitive sequences and 63 simple repeat repetitive sequences that did not contain reverse repeats. The simple repeats included 45 single-nucleotide repeats, which accounted for the highest proportion and primarily comprised A/T repeats. A comparative analysis of C. fruticosum, C. multijugum, and four Hedysarum species revealed that the six genomes were highly conserved, with differentials primarily located in the conserved non-coding regions. Moreover, the accD and clpP genes in the coding regions exhibited high nucleotide variability. Accordingly, these genes may serve as molecular markers for the classification and phylogenetic analysis of Corethrodendron species. Phylogenetic analysis further revealed that C. fruticosum and C. multijugum appeared in different clades than the four Hedysarum species. The newly sequenced chloroplast genome provides further insights into the phylogenetic position of C. fruticosum, which is useful for the classification and identification of Corethrodendron.

The Chloroplast Phylogenomics and Systematics of Zoysia (Poaceae)

The genus Zoysia Willd. (Chloridoideae) is widely distributed from the temperate regions of Northeast Asia—including China, Japan, and Korea—to the tropical regions of Southeast Asia. Among these, four species—Zoysia japonica Steud., Zoysia sinica Hance, Zoysia tenuifolia Thiele, and Zoysia macrostachya Franch. & Sav.—are naturally distributed in the Korean Peninsula. In this study, we report the complete plastome sequences of these Korean Zoysia species (NCBI acc. nos. MF953592, MF967579~MF967581). The length of Zoysia plastomes ranges from 135,854 to 135,904 bp, and the plastomes have a typical quadripartite structure, which consists of a pair of inverted repeat regions (20,962~20,966 bp) separated by a large (81,348~81,392 bp) and a small (12,582~12,586 bp) single-copy region. In terms of gene order and structure, Zoysia plastomes are similar to the typical plastomes of Poaceae. The plastomes encode 110 genes, of which 76 are protein-coding genes, 30 are tRNA genes, and four are rRNA genes. Fourteen genes contain single introns and one gene has two introns. Three evolutionary hotspot spacer regions—atpB~rbcL, rps16~rps3, and rpl32~trnL-UAG—were recognized among six analyzed Zoysia species. The high divergences in the atpB~rbcL spacer and rpl16~rpl3 region are primarily due to the differences in base substitutions and indels. In contrast, the high divergence between rpl32~trnL-UAG spacers is due to a small inversion with a pair of 22 bp stem and an 11 bp loop. Simple sequence repeats (SSRs) were identified in 59 different locations in Z. japonica, 63 in Z. sinica, 62 in Z. macrostachya, and 63 in Z. tenuifolia plastomes. Phylogenetic analysis showed that the Zoysia (Zoysiinae) forms a monophyletic group, which is sister to Sporobolus (Sporobolinae), with 100% bootstrap support. Within the Zoysia clade, the relationship of (Z. sinica, Z japonica), (Z. tenuifolia, Z. matrella), (Z. macrostachya, Z. macrantha) was suggested.

Genomic Profiling: The Strengths and Limitations of Chloroplast Genome-Based Plant Variety Authentication

Genomic profiling is a suitable tool for variety authentication and has applications in both operational quality and regulatory raw material control. It can be used to differentiate species or varieties and to identify admixtures as well as field contaminants. To establish a molecular profile, reliable and very accurate sequence data are required. As a result of the influence of the pollinator plant, nuclear genome-based authentication is in most cases not suitable for a direct application on the fruit. Sequences must be used that come exclusively from the localized mother plant. Parts of the fruit of maternal origin, e.g., components derived from the blossom, are suitable as a basis for this. Alternatively, DNA from cell organelles that are maternally inherited, such as mitochondria or chloroplasts, can be used. The latter will be discussed in this review in closer detail. Although individual gene segments on the chloroplast genome are already used for species differentiation in barcoding studies on plants, little is known about the usefulness of the entire chloroplast genome for intraspecies differentiation in general and for differentiation between modern varieties in particular. Results from the literature as well as from our own work suggest that chloroplast genome sequences are indeed very well-suited for the differentiation of old varieties. On the other hand, they are less or not suitable for the genetic differentiation of modern cultivars, because they are often too closely related.

Novel microsatellite markers for evaluation of genetic diversity in the tetraploid flame azalea, Rhododendron calendulaceum (Ericaceae)

Microsatellite markers have become a popular and useful tool for investigating evolutionary processes at shallow taxonomic scales such as within a species or between extremely closely related species. Rhododendron sect. Pentanthera is a closely related group of deciduous azaleas that demonstrate both naturally occurring and horticulturally derived hybridization. Two species, flame azalea and Cumberland azalea, represent a particularly recalcitrant evolutionary problem, which will benefit from the development of rapidly evolving molecular markers. Microsatellite markers were specifically developed for Rhododendron calendulaceum, the flame azalea, for use in studies of genetic structure and potential hybridization with its close relative Rhododendron cumberlandense, the Cumberland azalea. Forty-eight primer pairs designed from paired-end Illumina MiSeq data were screened for robust amplification. Sixteen of these pairs were PCR-amplified in the presence of fluorescently labeled primers and genotyped in 66 flame azalea individuals from three geographically dispersed populations. Fifteen primer pairs were both reliable and polymorphic and exhibit ample variability for use in downstream population-level investigations. Cross-amplification in all other members of Rhododendron sect. Pentanthera was highly successful, suggesting broad utility across the entire clade. The novel microsatellite markers presented here functioned well within the target species and amplified with high success in the remaining members of the clade. They represent a significant improvement to the genetic toolkit available for Rhododendron sect. Pentanthera, and particularly for the flame/Cumberland azalea evolutionary problem.

Isolation of Intact Chloroplast for Sequencing Plastid Genomes of Five Festuca Species

Isolation of good quality chloroplast DNA (cpDNA) is a challenge in different plant species, although several methods for isolation are known. Attempts were undertaken to isolate cpDNA from Festuca grass species by using available standard protocols; however, they failed due to difficulties separating intact chloroplasts from the polysaccharides, oleoresin, and contaminated nuclear DNA that are present in the crude homogenate. In this study, we present a quick and inexpensive protocol for isolating intact chloroplasts from seven grass varieties/accessions of five Festuca species using a single layer of 30% Percoll solution. This protocol was successful in isolating high quality cpDNA with the least amount of contamination of other DNA. We performed Illumina MiSeq paired-end sequencing (2 × 300 bp) using 200 ng of cpDNA of each variety/accession. Chloroplast genome mapping showed that 0.28%–11.37% were chloroplast reads, which covered 94%–96% of the reference plastid genomes of the closely related grass species. This improved method delivered high quality cpDNA from seven grass varieties/accessions of five Festuca species and could be useful for other grass species with similar genome complexity.

Fourteen polymorphic microsatellite markers for the widespread Labrador tea (Rhododendron groenlandicum)

PREMISE

Microsatellite markers were developed for Labrador tea (Rhododendron groenlandicum, Ericaceae) to facilitate downstream genetic investigation of this species and the extremely closely related, circumboreal Rhododendron subsect. Ledum.

METHODS AND RESULTS

Forty-eight primer pairs were designed using Illumina data and screened for excellent amplification. Sixteen successful pairs were developed as microsatellite markers using fluorescently labeled amplification to generate chromatogram data. These data were evaluated for intrapopulation and interpopulation variability in three populations from Alaska and Maine, USA, and the Northwest Territories, Canada. Fourteen polymorphic markers genotyped reliably, each with one to eight alleles. Cluster analysis indicates that across the range, populations can be easily discriminated. Cross-amplification in other Rhododendron subsect. Ledum species shows broad application of the developed markers within this small, well-supported clade.

CONCLUSIONS

These microsatellite markers exhibit significant variability and will be useful in population genetics within R. groenlandicum and for investigation of species boundaries across Rhododendron subsect. Ledum.

Microsatellite markers for the biogeographically enigmatic sandmyrtle (Kalmia buxifolia, Phyllodoceae: Ericaceae)

PREMISE

Microsatellite markers were developed for sandmyrtle, Kalmia buxifolia (Ericaceae), to facilitate phylogeographic studies in this taxon and possibly many of its close relatives.

METHODS AND RESULTS

Forty-eight primer pairs designed from paired-end Illumina MiSeq data were screened for robust amplification. Sixteen pairs were amplified again, but with fluorescently labeled primers to facilitate genotyping. Resulting chromatograms were evaluated for variability using three populations from Tennessee, North Carolina, and New Jersey, USA. Eleven primer pairs were reliable and polymorphic (mean 3.92 alleles), one was reliable but monomorphic, and four were not reliable. The markers exhibited lower heterozygosity (mean 0.246) than expected (mean 0.464). Cross-amplification in the remaining nine Kalmia species exhibited a phylogenetic pattern, suggesting broad applicability of the markers across the genus.

CONCLUSIONS

These microsatellite markers will be useful in population genetics and species boundaries studies of K. buxifolia, K. procumbens, and likely all other Kalmia species.

Evolution of six novel ORFs in the plastome of Mankyua chejuense and phylogeny of eusporangiate ferns

In this paper, three plastomes of Mankyua chejuense, Helminthostachys zeylanica, and Botrychium ternatum in Ophioglossaceae were completely sequenced in order to investigate the plastome evolution and phylogeny of eusporangiate ferns. They were similar to each other in terms of length and the gene orders; however, six unknown open reading frames (ORFs) were found between rps4 and trnL-UAA genes in M. chejuense. Similar sequence regions of six ORFs of M. chejuense were found at the plastomes of Ophioglossum californicum and H. zeylanica, as well as the mitochondrial genome (mitogenome) of H. zeylanica, but not in B. ternatum. Interestingly, the translated amino acid sequences of three ORFs were more similar to the proteins of distantly related taxa such as algae and bacteria than they were to proteins in land plants. It is likely that the six ORFs region arose from endosymbiotic gene transfer (EGT) or horizontal gene transfer (HGT), but further study is needed to verify this. Phylogenetic analyses suggested that Mankyua was resolved as the earliest diverging lineage and that Ophioglossum was subsequently diverged in Ophioglossaceae. This result supports why the plastome of M. chejuense have contained the most ancestral six ORFs in the family.

The complete plastome of tropical fruit Garcinia mangostana (Clusiaceae)

The complete plastome sequence of Garcinia mangostana L. (Clusiaceae) is completed in this study (NCBI acc. no. KX822787). This is a first complete plastome sequence from the Clusiaceae. The complete plastome size is 158,179 bp in length and consists of a large single copy of 86,458 bp and a small single copy of 17,703 bp, separated by two inverted repeats of 27,009 bp. The G. mangostana plastome shows four minor structural modifications including infA gene loss, rpl32 gene loss, ycf3 gene intron loss and a 363 bp inversion between trnV-UAC and atpE gene. The plastome contains 111 genes, of which 77 are protein-coding genes, 30 are tRNA genes and four are rRNA genes. The average A-T content of the plastome is 63.9%. A total of 110 simple sequence loci are identified from the genome. Phylogenetic analysis reveals that G. mangostana is a sister group of Erythroxylum novogranatense (Erythroxylaceae) with 78% bootstrap support.

The complete plastome sequences of Mangifera indica L. (Anacardiaceae)

In this study, we determined the complete plastome sequence of Mangifera indica L. (Anacardiaceae) (NCBI acc. no. KX871231). The complete plastome is 157,780 bp in length, and consists of a large single copy of 86,673 bp and a small single copy of 18,349 bp, separated by two inverted repeats of 25,792 bp. The plastome contains 112 genes, of which 78 are protein-coding genes, 30 are tRNA genes, and four are rRNA genes. Sixteen genes contain one intron and two genes have two introns. The average A-T content of the plastome is 62.1%. The M. indica plastome has approximately 15 kb inversion between trnT-UGU and trnT-GGU. We identify a palindromic repeat of 18 bp (ATTCTTTTTTTTTTTTTT/AAAAAAAAAAAAAAGAAT) near the inversion breakpoints of M. indica plastome. Phylogenetic analysis revealed that M. indica is a sister group of Rhus chinensis with 100% bootstrap support. Anacardiaceae clade is a sister group of Boswellia sacra (Burseraceae) with 100% bootstrap support.

SIX INDEPENDENT LOSSES OF THE CHLOROPLAST DNA rpl2 INTRON IN DICOTYLEDONS: MOLECULAR AND PHYLOGENETIC IMPLICATIONS

Previous studies have shown that in several angiosperms and the liverwort Marchantia the chloroplast gene rpl2, encoding ribosomal protein L2, is interrupted by an intron, but that in spinach (Spinacia oleracea, Caryophyllales) this intron has been lost. We have determined the distribution of the rpl2 intron for 390 species representing 116 angiosperm families. Filter hybridizations reveal that the intron is absent from the chloroplast genomes of all examined families of the Caryophyllales, suggesting that the intron was lost in the common ancestor of the order. Sequencing of the rpl2 gene in five genera of the Caryophyllales and in Rumex (Polygonales) not only confirms the filter hybridization results, but also shows that for all taxa lacking the intron, the rpl2 gene has undergone a precise deletion of the intron. In all cases, it is the original rpl2 gene that has sustained loss of its intron. This implies that in chloroplast DNA, integration of exogenous genes (e.g., a reverse transcript of a spliced mRNA) occurs mainly by homologous, replacement recombination, rather than by illegitimate recombination elsewhere in the genome. Filter hybridizations also reveal that the rpl2 intron was lost independently in the common ancestors of at least five other lineages of dicotyledons: Saxifragaceae (s.s.), Convolvulaceae (including Cuscuta), Menyanthaceae, two genera of Geraniaceae, and one genus of Droseraceae. The molecular and phylogenetic implications of these independent intron losses are discussed.

Application of a simplified method of chloroplast enrichment to small amounts of tissue for chloroplast genome sequencing

PREMISE OF THE STUDY

High-throughput sequencing of genomic DNA can recover complete chloroplast genome sequences, but the sequence data are usually dominated by sequences from nuclear/mitochondrial genomes. To overcome this deficiency, a simple enrichment method for chloroplast DNA from small amounts of plant tissue was tested for eight plant species including a gymnosperm and various angiosperms.

METHODS

Chloroplasts were enriched using a high-salt isolation buffer without any step gradient procedures, and enriched chloroplast DNA was sequenced by multiplexed high-throughput sequencing.

RESULTS

Using this simple method, significant enrichment of chloroplast DNA-derived reads was attained, allowing deep sequencing of chloroplast genomes. As an example, the chloroplast genome of the conifer Callitris sulcata was assembled, from which polymorphic microsatellite loci were isolated successfully.

DISCUSSION

This chloroplast enrichment method from small amounts of plant tissue will be particularly useful for studies that use sequencers with relatively small throughput and that cannot use large amounts of tissue (e.g., for endangered species).

The complete plastome sequence of Carissa macrocarpa (Eckl.) A. DC. (Apocynaceae)

In this study, we determined the complete plastome sequence of Carissa macrocarpa (Eckl.) A. DC. (Apocynaceae) (NCBI acc. no. KX364402). The gene order and structure of the C. macrocarpa plastome are similar to those of a typical angiosperm. The complete plastome is 155,297 bp in length, and consists of a large single-copy region of 85,586 bp and a small single-copy region of 18,131 bp, which are separated by two inverted repeats of 25,792 bp. The plastome contains 113 genes, of which 79 are protein-coding genes, 30 are tRNA genes and 4 are rRNA genes. Sixteen genes contained one intron and two genes have two introns. The average A–T content of the plastome is 62.0%. A total of 31 simple sequence repeat loci were identified within the genome. Phylogenetic analysis revealed that C. macrocarpa is a member of the paraphyletic subfamily Rauvolfioideae of Apocynaceae. The sister group relationship of C. macrocarpa to the Apocynoideae–Asclepiadoideae clade is supported by 100% bootstrap values.

The first complete plastome sequence from the family Sapotaceae, Pouteria campechiana (Kunth) Baehni

In this study, we determined the complete plastome sequence of Pouteria campechiana (Kunth) Baehni (Sapotaceae) (NCBI acc. no. KX426215). This is the first time a plastome from the Sapotaceae has been sequenced. The gene order and structure of the P. campechiana plastome are collinear with those of the typical plastome of land plants. The complete plastome size is 157,922 bp in length and consists of a large single-copy region of 87,122 bp and a small single-copy region of 18,559 bp, which are separated by a pair of 26,120 bp-long inverted repeat regions. The overall A-T content of the plastome sequence is 63.2%. The plastome contains 113 genes, of which 79 are protein-coding genes, 30 are tRNA genes, and 4 are rRNA genes. Sixteen genes contain one intron and two genes have two introns. A total of 91 simple sequence repeat loci were identified within the genome. Phylogenetic analysis revealed that P. campechiana is a sister group of the Primulaceae-Ebenaceae clade with 100% bootstrap support.

The complete plastome sequence of Pentactina rupicola Nakai (Rosaceae), a genus endemic to Korea

The complete plastid genome (plastome) of Pentactina rupicola Nakai, the sole member of genus Pentactina, endemic to Korea, was determined in this study. The plastome of P. rupicola is 156,612 bp in length and is composed of a pair of 26,351 bp inverted repeat regions (IRa and IRb) separating large (LSC) and small (SSC) single-copy regions of 84,970 and 18,940 bp, respectively (NCBI acc. no. NC 016921). The plastome encodes 129 genes, of which 112 are unique, including 78 protein-coding genes, 30 tRNA genes, and 4 rRNA genes. Seventeen genes contain one intron and the ycf3 gene has two introns. The second intron of clpP is absent in the P. rupicola plastome. The AT content of P. rupicola is 63% overall, and in the LSC, SSC, and IR regions is 65%, 69%, and 58%, respectively. A total of 63 simple sequence repeats (SSR) are distributed among the noncoding regions of the genome. Phylogenetic analysis of the combined 82-gene data set for 35 plastomes suggests that P. rupicola is sister to the Pyrusmalus clade.

Complete plastome sequence of Psidium guajava L. (Myrtaceae)

In this study, we determined the complete plastome sequence of Psidium guajava L. (Myrtaceae) (NCBI acc. no. KX364403). The gene order and structure of the P. guajava plastome are similar to those of a typical angiosperm. The complete plastome is 158,841 bp in length, and consists of a large single copy of 87,675 bp and a small single copy of 18,464 bp, separated by two inverted repeats of 26,351 bp. The overall AT content of the sequence is 63.0%. The plastome contains 112 genes, of which 78 are protein-coding genes, 30 are tRNA genes, and four are rRNA genes. Sixteen genes contain one intron and two genes have two introns. A total of 100 simple sequence loci were identified from the genome. Phylogenetic analysis revealed that P. guajava is a sister group of Eugenia uniflora with 100% bootstrap support.

Complete plastome sequence of Averrhoa carambola L. (Oxalidaceae)

In this study, we determined the complete plastome sequence of Averrhoa carambola L. (Oxalidaceae) (NCBI acc. no. KX364202). To the best of our knowledge, this is the first reported complete plastome sequence from the order Oxalidales. The gene order and structure of the A. carambola plastome are collinear with the typical plastome of land plants. The complete plastome size is 155,965 bp in length and consists of a large single copy region of 87,217 bp and a small single copy region of 17,496 bp, which are separated by a pair of 25,626-bp-long inverted repeat regions. The overall A-T content of the plastome sequence is 61.2%. The plastome contains 111 genes, of which 77 are protein-coding genes, 30 are tRNA genes, and 4 are rRNA genes. Sixteen genes contain one intron and two genes have two introns. A total of 91 simple sequence loci were identified from the genome. Phylogenetic analysis revealed that A. carambola is a sister group of Euonymus japonicus with 100% bootstrap support.

The complete plastome sequence of Diospyros blancoi A. DC. (Ebenaceae)

The plastome sequences of Diospyros blancoi A. DC. (Ebenaceae) were completed in this study (NCBI acc. no. KX426216). The gene order and structure of the D. blancoi plastome are collinear with the typical plastome of land plants. The complete plastome size is 157,745 bp in length and consists of a large single-copy region of 87,246 bp and a small single-copy region of 18,323 bp, which are separated by a pair of 26,088 bp-long inverted repeat regions. The overall A-T content of the plastome sequence is 62.6%. The plastome contains 113 genes, of which 79 are protein-coding genes, 30 are tRNA genes, and 4 are rRNA genes. Sixteen genes contain one intron and two genes have two introns. A total of 45 simple sequence loci were identified from the genome. Phylogenetic analysis revealed that D. blancoi is a sister group of Primulaceae with 100% bootstrap support.

Correlation between sequence divergence and polymorphism reveals similar evolutionary mechanisms acting across multiple timescales in a rapidly evolving plastid genome

Background

Although the plastid genome is highly conserved across most angiosperms, multiple lineages have increased rates of structural rearrangement and nucleotide substitution. These lineages exhibit an excess of nonsynonymous substitutions (i.e., elevated dN/dS ratios) in similar subsets of plastid genes, suggesting that similar mechanisms may be leading to relaxed and/or positive selection on these genes. However, little is known regarding whether these mechanisms continue to shape sequence diversity at the intraspecific level.

Results

We examined patterns of interspecific divergence and intraspecific polymorphism in the plastid genome of Campanulastrum americanum, and across plastid genes found a significant correlation between dN/dS and pN/pS (i.e., the within-species equivalent of dN/dS). A number of genes including ycf1, ycf2, clpP, and ribosomal protein genes exhibited high dN/dS ratios. McDonald-Kreitman tests detected little evidence for positive selection acting on these genes, likely due to the presence of substantial intraspecific divergence.

Conclusions

These results suggest that mechanisms leading to increased nucleotide substitution rates in the plastid genome are continuing to act at the intraspecific level. Accelerated plastid genome evolution may increase the likelihood of intraspecific cytonuclear genetic incompatibilities, and thereby contribute to the early stages of the speciation process.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0268-y) contains supplementary material, which is available to authorized users.

Chloroplast genome differences between Asian and American Equisetum arvense (Equisetaceae) and the origin of the hypervariable trnY-trnE intergenic spacer

Comparative analyses of complete chloroplast (cp) DNA sequences within a species may provide clues to understand the population dynamics and colonization histories of plant species. Equisetum arvense (Equisetaceae) is a widely distributed fern species in northeastern Asia, Europe, and North America. The complete cp DNA sequences from Asian and American E. arvense individuals were compared in this study. The Asian E. arvense cp genome was 583 bp shorter than that of the American E. arvense. In total, 159 indels were observed between two individuals, most of which were concentrated on the hypervariable trnY-trnE intergenic spacer (IGS) in the large single-copy (LSC) region of the cp genome. This IGS region held a series of 19 bp repeating units. The numbers of the 19 bp repeat unit were responsible for 78% of the total length difference between the two cp genomes. Furthermore, only other closely related species of Equisetum also show the hypervariable nature of the trnY-trnE IGS. By contrast, only a single indel was observed in the gene coding regions: the ycf1 gene showed 24 bp differences between the two continental individuals due to a single tandem-repeat indel. A total of 165 single-nucleotide polymorphisms (SNPs) were recorded between the two cp genomes. Of these, 52 SNPs (31.5%) were distributed in coding regions, 13 SNPs (7.9%) were in introns, and 100 SNPs (60.6%) were in intergenic spacers (IGS). The overall difference between the Asian and American E. arvense cp genomes was 0.12%. Despite the relatively high genetic diversity between Asian and American E. arvense, the two populations are recognized as a single species based on their high morphological similarity. This indicated that the two regional populations have been in morphological stasis.

The dynamic history of plastid genomes in the Campanulaceae sensu lato is unique among angiosperms

Why have some plants lost the organizational stability in plastid genomes (plastomes) that evolved in their algal ancestors? During the endosymbiotic transformation of a cyanobacterium into the eukaryotic plastid, most cyanobacterial genes were transferred to the nucleus or otherwise lost from the plastome, and the resulting plastome architecture in land plants confers organizational stability, as evidenced by the conserved gene order among bryophytes and lycophytes, whereas ferns, gymnosperms, and angiosperms share a single, 30-kb inversion. Although some additional gene losses have occurred, gene additions to angiosperm plastomes were previously unknown. Plastomes in the Campanulaceae sensu lato have incorporated dozens of large ORFs (putative protein-coding genes). These insertions apparently caused many of the 125+ large inversions now known in this small eudicot clade. This phylogenetically restricted phenomenon is not biogeographically localized, which indicates that these ORFs came from the nucleus or (less likely) a cryptic endosymbiont.

Chloroplast genome evolution in early diverged leptosporangiate ferns

In this study, the chloroplast (cp) genome sequences from three early diverged leptosporangiate ferns were completed and analyzed in order to understand the evolution of the genome of the fern lineages. The complete cp genome sequence of Osmunda cinnamomea (Osmundales) was 142,812 base pairs (bp). The cp genome structure was similar to that of eusporangiate ferns. The gene/intron losses that frequently occurred in the cp genome of leptosporangiate ferns were not found in the cp genome of O. cinnamomea. In addition, putative RNA editing sites in the cp genome were rare in O. cinnamomea, even though the sites were frequently predicted to be present in leptosporangiate ferns. The complete cp genome sequence of Diplopterygium glaucum (Gleicheniales) was 151,007 bp and has a 9.7 kb inversion between the trnL-CAA and trnVGCA genes when compared to O. cinnamomea. Several repeated sequences were detected around the inversion break points. The complete cp genome sequence of Lygodium japonicum (Schizaeales) was 157,142 bp and a deletion of the rpoC1 intron was detected. This intron loss was shared by all of the studied species of the genus Lygodium. The GC contents and the effective numbers of codons (ENCs) in ferns varied significantly when compared to seed plants. The ENC values of the early diverged leptosporangiate ferns showed intermediate levels between eusporangiate and core leptosporangiate ferns. However, our phylogenetic tree based on all of the cp gene sequences clearly indicated that the cp genome similarity between O. cinnamomea (Osmundales) and eusporangiate ferns are symplesiomorphies, rather than synapomorphies. Therefore, our data is in agreement with the view that Osmundales is a distinct early diverged lineage in the leptosporangiate ferns.

A recurring syndrome of accelerated plastid genome evolution in the angiosperm tribe Sileneae (Caryophyllaceae)

In flowering plants, plastid genomes are generally conserved, exhibiting slower rates of sequence evolution than the nucleus and little or no change in structural organization. However, accelerated plastid genome evolution has occurred in scattered angiosperm lineages. For example, some species within the genus Silene have experienced a suite of recent changes to their plastid genomes, including inversions, shifts in inverted repeat boundaries, large indels, intron losses, and rapid rates of amino acid sequence evolution in a subset of protein genes, with the most extreme divergence occurring in the protease gene clpP. To investigate the relationship between the rates of sequence and structural evolution, we sequenced complete plastid genomes from three species (Silene conoidea, S. paradoxa, and Lychnis chalcedonica), representing independent lineages within the tribe Sileneae that were previously shown to have accelerated rates of clpP evolution. We found a high degree of parallel evolution. Elevated rates of amino acid substitution have occurred repeatedly in the same subset of plastid genes and have been accompanied by a recurring pattern of structural change, including cases of identical inversions and intron loss. This "syndrome" of changes was not observed in the closely related outgroup Agrostemma githago or in the more slowly evolving Silene species that were sequenced previously. Although no single mechanism has yet been identified to explain the correlated suite of changes in plastid genome sequence and structure that has occurred repeatedly in angiosperm evolution, we discuss a possible mixture of adaptive and non-adaptive forces that may be responsible.

An improved chloroplast DNA extraction procedure for whole plastid genome sequencing

BACKGROUND

Chloroplast genomes supply valuable genetic information for evolutionary and functional studies in plants. The past five years have witnessed a dramatic increase in the number of completely sequenced chloroplast genomes with the application of second-generation sequencing technology in plastid genome sequencing projects. However, cost-effective high-throughput chloroplast DNA (cpDNA) extraction becomes a major bottleneck restricting the application, as conventional methods are difficult to make a balance between the quality and yield of cpDNAs.

METHODOLOGY/PRINCIPAL FINDINGS

We first tested two traditional methods to isolate cpDNA from the three species, Oryza brachyantha, Leersia japonica and Prinsepia utihis. Both of them failed to obtain properly defined cpDNA bands. However, we developed a simple but efficient method based on sucrose gradients and found that the modified protocol worked efficiently to isolate the cpDNA from the same three plant species. We sequenced the isolated DNA samples with Illumina (Solexa) sequencing technology to test cpDNA purity according to aligning sequence reads to the reference chloroplast genomes, showing that the reference genome was properly covered. We show that 40-50% cpDNA purity is achieved with our method.

CONCLUSION

Here we provide an improved method used to isolate cpDNA from angiosperms. The Illumina sequencing results suggest that the isolated cpDNA has reached enough yield and sufficient purity to perform subsequent genome assembly. The cpDNA isolation protocol thus will be widely applicable to the plant chloroplast genome sequencing projects.

Recent acceleration of plastid sequence and structural evolution coincides with extreme mitochondrial divergence in the angiosperm genus Silene

The angiosperm genus Silene exhibits some of the most extreme and rapid divergence ever identified in mitochondrial genome architecture and nucleotide substitution rates. These patterns have been considered mitochondrial specific based on the absence of correlated changes in the small number of available nuclear and plastid gene sequences. To better assess the relationship between mitochondrial and plastid evolution, we sequenced the plastid genomes from four Silene species with fully sequenced mitochondrial genomes. We found that two species with fast-evolving mitochondrial genomes, S. noctiflora and S. conica, also exhibit accelerated rates of sequence and structural evolution in their plastid genomes. The nature of these changes, however, is markedly different from those in the mitochondrial genome. For example, in contrast to the mitochondrial pattern, which appears to be genome wide and mutationally driven, the plastid substitution rate accelerations are restricted to a subset of genes and preferentially affect nonsynonymous sites, indicating that altered selection pressures are acting on specific plastid-encoded functions in these species. Indeed, some plastid genes in S. noctiflora and S. conica show strong evidence of positive selection. In contrast, two species with more slowly evolving mitochondrial genomes, S. latifolia and S. vulgaris, have correspondingly low rates of nucleotide substitution in plastid genes as well as a plastid genome structure that has remained essentially unchanged since the origin of angiosperms. These results raise the possibility that common evolutionary forces could be shaping the extreme but distinct patterns of divergence in both organelle genomes within this genus.

Development and characterization of polymorphic microsatellite markers for Conopholis americana (Orobanchaceae)

PREMISE OF THE STUDY

Conopholis americana is an obligate root parasite with highly reduced morphology. To investigate the population structure, genetic diversity, and mating system of this predominantly eastern North American species, we developed polymorphic microsatellite markers for C. americana.

METHODS AND RESULTS

Using an enrichment cloning protocol, we report the isolation and characterization of 11 microsatellite markers. Product size varied from 198-370 bp. These loci show moderate levels of allelic variation (averaging 4.182 alleles per locus) and very low levels of heterozygosity (average observed heterozygosity = 0.054).

CONCLUSIONS

These microsatellite markers will be useful in obtaining estimates of population-level genetic diversity and in phylogeographic studies of C. americana.

Implications of the plastid genome sequence of typha (typhaceae, poales) for understanding genome evolution in poaceae

Plastid genomes of the grasses (Poaceae) are unusual in their organization and rates of sequence evolution. There has been a recent surge in the availability of grass plastid genome sequences, but a comprehensive comparative analysis of genome evolution has not been performed that includes any related families in the Poales. We report on the plastid genome of Typha latifolia, the first non-grass Poales sequenced to date, and we present comparisons of genome organization and sequence evolution within Poales. Our results confirm that grass plastid genomes exhibit acceleration in both genomic rearrangements and nucleotide substitutions. Poaceae have multiple structural rearrangements, including three inversions, three genes losses (accD, ycf1, ycf2), intron losses in two genes (clpP, rpoC1), and expansion of the inverted repeat (IR) into both large and small single-copy regions. These rearrangements are restricted to the Poaceae, and IR expansion into the small single-copy region correlates with the phylogeny of the family. Comparisons of 73 protein-coding genes for 47 angiosperms including nine Poaceae genera confirm that the branch leading to Poaceae has significantly accelerated rates of change relative to other monocots and angiosperms. Furthermore, rates of sequence evolution within grasses are lower, indicating a deceleration during diversification of the family. Overall there is a strong correlation between accelerated rates of genomic rearrangements and nucleotide substitutions in Poaceae, a phenomenon that has been noted recently throughout angiosperms. The cause of the correlation is unknown, but faulty DNA repair has been suggested in other systems including bacterial and animal mitochondrial genomes.

Comparative analysis of the complete sequence of the plastid genome of Parthenium argentatum and identification of DNA barcodes to differentiate Parthenium species and lines

BACKGROUND

Parthenium argentatum (guayule) is an industrial crop that produces latex, which was recently commercialized as a source of latex rubber safe for people with Type I latex allergy. The complete plastid genome of P. argentatum was sequenced. The sequence provides important information useful for genetic engineering strategies. Comparison to the sequences of plastid genomes from three other members of the Asteraceae, Lactuca sativa, Guitozia abyssinica and Helianthus annuus revealed details of the evolution of the four genomes. Chloroplast-specific DNA barcodes were developed for identification of Parthenium species and lines.

RESULTS

The complete plastid genome of P. argentatum is 152,803 bp. Based on the overall comparison of individual protein coding genes with those in L. sativa, G. abyssinica and H. annuus, we demonstrate that the P. argentatum chloroplast genome sequence is most closely related to that of H. annuus. Similar to chloroplast genomes in G. abyssinica, L. sativa and H. annuus, the plastid genome of P. argentatum has a large 23 kb inversion with a smaller 3.4 kb inversion, within the large inversion. Using the matK and psbA-trnH spacer chloroplast DNA barcodes, three of the four Parthenium species tested, P. tomentosum, P. hysterophorus and P. schottii, can be differentiated from P. argentatum. In addition, we identified lines within P. argentatum.

CONCLUSION

The genome sequence of the P. argentatum chloroplast will enrich the sequence resources of plastid genomes in commercial crops. The availability of the complete plastid genome sequence may facilitate transformation efficiency by using the precise sequence of endogenous flanking sequences and regulatory elements in chloroplast transformation vectors. The DNA barcoding study forms the foundation for genetic identification of commercially significant lines of P. argentatum that are important for producing latex.

ANALYSIS OF A PLASTID MULTIGENE DATA SET AND THE PHYLOGENETIC POSITION OF THE MARINE MACROALGA CAULERPA FILIFORMIS (CHLOROPHYTA)(1)

Molecular phylogenetic relationships within the Chlorophyta have relied heavily on rRNA data. These data have revolutionized our insight in green algal evolution, yet some class relationships have never been well resolved. A commonly used class within the Chlorophyta is the Ulvophyceae, although there is not much support for its monophyly. The relationships among the Ulvophyceae, Trebouxiophyceae, and Chlorophyceae are also contentious. In recent years, chloroplast genome data have shown their utility in resolving relationships between the main green algal clades, but such studies have never included marine macroalgae. We provide partial chloroplast genome data (∼30,000 bp, 23 genes) of the ulvophycean macroalga Caulerpa filiformis (Suhr) K. Herig. We show gene order conservation for some gene combinations and rearrangements in other regions compared to closely related taxa. Our data also revealed a pseudogene (ycf62) in Caulerpa species. Our phylogenetic results, based on analyses of a 23-gene alignment, suggest that neither Ulvophyceae nor Trebouxiophyceae are monophyletic, with Caulerpa being more closely related to the trebouxiophyte Chlorella than to Oltmannsiellopsis and Pseudendoclonium.

An optimized chloroplast DNA extraction protocol for grasses (Poaceae) proves suitable for whole plastid genome sequencing and SNP detection

Background

Obtaining chloroplast genome sequences is important to increase the knowledge about the fundamental biology of plastids, to understand evolutionary and ecological processes in the evolution of plants, to develop biotechnological applications (e.g. plastid engineering) and to improve the efficiency of breeding schemes. Extraction of pure chloroplast DNA is required for efficient sequencing of chloroplast genomes. Unfortunately, most protocols for extracting chloroplast DNA were developed for eudicots and do not produce sufficiently pure yields for a shotgun sequencing approach of whole plastid genomes from the monocot grasses.

Methodology/Principal Findings

We have developed a simple and inexpensive method to obtain chloroplast DNA from grass species by modifying and extending protocols optimized for the use in eudicots. Many protocols for extracting chloroplast DNA require an ultracentrifugation step to efficiently separate chloroplast DNA from nuclear DNA. The developed method uses two more centrifugation steps than previously reported protocols and does not require an ultracentrifuge.

Conclusions/Significance

The described method delivered chloroplast DNA of very high quality from two grass species belonging to highly different taxonomic subfamilies within the grass family (Lolium perenne, Pooideae; Miscanthus×giganteus, Panicoideae). The DNA from Lolium perenne was used for whole chloroplast genome sequencing and detection of SNPs. The sequence is publicly available on EMBL/GenBank.

Extensive rearrangements in the chloroplast genome of Trachelium caeruleum are associated with repeats and tRNA genes

Chloroplast genome organization, gene order, and content are highly conserved among land plants. We sequenced the chloroplast genome of Trachelium caeruleum L. (Campanulaceae), a member of an angiosperm family known for highly rearranged genomes. The total genome size is 162,321 bp, with an inverted repeat (IR) of 27,273 bp, large single-copy (LSC) region of 100,114 bp, and small single-copy (SSC) region of 7,661 bp. The genome encodes 112 different genes, with 17 duplicated in the IR, a tRNA gene (trnI-cau) duplicated once in the LSC region, and a protein-coding gene (psbJ) with two duplicate copies, for a total of 132 putatively intact genes. ndhK may be a pseudogene with internal stop codons, and clpP, ycf1, and ycf2 are so highly diverged that they also may be pseudogenes. ycf15, rpl23, infA, and accD are truncated and likely nonfunctional. The most conspicuous feature of the Trachelium genome is the presence of 18 internally unrearranged blocks of genes inverted or relocated within the genome relative to the ancestral gene order of angiosperm chloroplast genomes. Recombination between repeats or tRNA genes has been suggested as a mechanism of chloroplast genome rearrangements. The Trachelium chloroplast genome shares with Pelargonium and Jasminum both a higher number of repeats and larger repeated sequences in comparison to eight other angiosperm chloroplast genomes, and these are concentrated near rearrangement endpoints. Genes for tRNAs occur at many but not all inversion endpoints, so some combination of repeats and tRNA genes may have mediated these rearrangements.

Phylogenetic and evolutionary implications of complete chloroplast genome sequences of four early-diverging angiosperms: Buxus (Buxaceae), Chloranthus (Chloranthaceae), Dioscorea (Dioscoreaceae), and Illicium (Schisandraceae)

We have determined the complete chloroplast genome sequences of four early-diverging lineages of angiosperms, Buxus (Buxaceae), Chloranthus (Chloranthaceae), Dioscorea (Dioscoreaceae), and Illicium (Schisandraceae), to examine the organization and evolution of plastid genomes and to estimate phylogenetic relationships among angiosperms. For the most part, the organization of these plastid genomes is quite similar to the ancestral angiosperm plastid genome with a few notable exceptions. Dioscorea has lost one protein-coding gene, rps16; this gene loss has also happened independently in four other land plant lineages, liverworts, conifers, Populus, and legumes. There has also been a small expansion of the inverted repeat (IR) in Dioscorea that has duplicated trnH-GUG. This event has also occurred multiple times in angiosperms, including in monocots, and in the two basal angiosperms Nuphar and Drimys. The Illicium chloroplast genome is unusual by having a 10 kb contraction of the IR. The four taxa sequenced represent key groups in resolving phylogenetic relationships among angiosperms. Illicium is one of the basal angiosperms in the Austrobaileyales, Chloranthus (Chloranthales) remains unplaced in angiosperm classifications, and Buxus and Dioscorea are early-diverging eudicots and monocots, respectively. We have used sequences for 61 shared protein-coding genes from these four genomes and combined them with sequences from 35 other genomes to estimate phylogenetic relationships using parsimony, likelihood, and Bayesian methods. There is strong congruence among the trees generated by the three methods, and most nodes have high levels of support. The results indicate that Amborella alone is sister to the remaining angiosperms; the Nymphaeales represent the next-diverging clade followed by Illicium; Chloranthus is sister to the magnoliids and together this group is sister to a large clade that includes eudicots and monocots; and Dioscorea represents an early-diverging lineage of monocots just internal to Acorus.

Complete chloroplast genome sequences of Hordeum vulgare, Sorghum bicolor and Agrostis stolonifera, and comparative analyses with other grass genomes

Comparisons of complete chloroplast genome sequences of Hordeum vulgare, Sorghum bicolor and Agrostis stolonifera to six published grass chloroplast genomes reveal that gene content and order are similar but two microstructural changes have occurred. First, the expansion of the IR at the SSC/IRa boundary that duplicates a portion of the 5' end of ndhH is restricted to the three genera of the subfamily Pooideae (Agrostis, Hordeum and Triticum). Second, a 6 bp deletion in ndhK is shared by Agrostis, Hordeum, Oryza and Triticum, and this event supports the sister relationship between the subfamilies Erhartoideae and Pooideae. Repeat analysis identified 19-37 direct and inverted repeats 30 bp or longer with a sequence identity of at least 90%. Seventeen of the 26 shared repeats are found in all the grass chloroplast genomes examined and are located in the same genes or intergenic spacer (IGS) regions. Examination of simple sequence repeats (SSRs) identified 16-21 potential polymorphic SSRs. Five IGS regions have 100% sequence identity among Zea mays, Saccharum officinarum and Sorghum bicolor, whereas no spacer regions were identical among Oryza sativa, Triticum aestivum, H. vulgare and A. stolonifera despite their close phylogenetic relationship. Alignment of EST sequences and DNA coding sequences identified six C-U conversions in both Sorghum bicolor and H. vulgare but only one in A. stolonifera. Phylogenetic trees based on DNA sequences of 61 protein-coding genes of 38 taxa using both maximum parsimony and likelihood methods provide moderate support for a sister relationship between the subfamilies Erhartoideae and Pooideae.

A comparative analysis of the Lactuca and Helianthus (Asteraceae) plastid genomes: identification of divergent regions and categorization of shared repeats

We have sequenced two complete chloroplast genomes in the Asteraceae, Helianthus annuus (sunflower), and Lactuca sativa (lettuce), which belong to the distantly related subfamilies, Asteroideae and Cichorioideae, respectively. The Helianthus chloroplast genome is 151 104 bp and the Lactuca genome is 152 772 bp long, which is within the usual size range for chloroplast genomes in flowering plants. When compared to tobacco, both genomes have two inversions: a large 22.8-kb inversion and a smaller 3.3-kb inversion nested within it. Pairwise sequence divergence across all genes, introns, and spacers in Helianthus and Lactuca has resulted in the discovery of new, fast-evolving DNA sequences for use in species-level phylogenetics, such as the trnY-rpoB, trnL-rpl32, and ndhC-trnV spacers. Analysis and categorization of shared repeats resulted in seven classes useful for future repeat studies: double tandem repeats, three or more tandem repeats, direct repeats dispersed in the genome, repeats found in reverse complement orientation, hairpin loops, runs of A's or T's in excess of 12 bp, and gene or tRNA similarity. Results from BLAST searches of our genomic sequence against expressed sequence tag (EST) databases for both genomes produced eight likely RNA edited sites (C → U changes). These detailed analyses in Asteraceae contribute to a broader understanding of plastid evolution across flowering plants.

Paternal, maternal, and biparental inheritance of the chloroplast genome in Passiflora (Passifloraceae): implications for phylogenetic studies

Patterns of inheritance of the chloroplast genome in Passiflora were analyzed by examining the progeny from both interspecific and intraspecific crosses. Artificial crosses of field-collected material were performed in greenhouses at The University of Texas at Austin. DNA from fresh leaf material was analyzed by Southern blot techniques to identify the donor of the chloroplast genome. Initially, single progeny were analyzed for 11 crosses; two intraspecific crosses demonstrated maternal inheritance, whereas the nine interspecific crosses had paternal inheritance. Subsequently, the donor of the chloroplast genome was determined for multiple progeny in seven crosses. Passiflora oerstedii × P. retipetala showed strict paternal inheritance in all of 17 progeny. A series of five crosses and backcrosses between P. oerstedii and P. menispermifolia demonstrated strictly paternal inheritance. Finally, when 15 progeny were analyzed for the P. costaricensis × P. costaricensis cross, 12 of the 15 showed maternal inheritance, whereas the remaining three were biparental. Interestingly, all interspecific crosses had primarily paternal inheritance, whereas all intraspecific crosses had primarily maternal inheritance. The implications of heteroplasmy on phylogenetic analyses of chloroplast DNA are discussed.

The phylogenetic history and biogeography of the frankincense and myrrh family (Burseraceae) based on nuclear and chloroplast sequence data

Generalized hypotheses for the vicariant, Gondwanan origin of pantropically distributed eudicotyledon families must be refined to accommodate recently revised dates that indicate major continental rifting events predate the evolution of many tricolpate angiosperm clades. Here, we use molecular phylogenies of an eudicotyledon family previously hypothesized to have a Gondwanan origin, the Burseraceae, to test this and other alternative biogeographical hypotheses in light of recalibrated geological events. Phylogenies based on nuclear and chloroplast data were reconstructed for 13 of the 18 genera (50 spp. total) of Burseraceae using parsimony, maximum likelihood, and Bayesian methods. Ages of all lineages were estimated using penalized likelihood and semiparametric rate smoothing [Bioinformatics 2003 (19) 301], which allows the user to calibrate phylogenies based on non-clock-like DNA sequence data with fossil information. Biogeographical hypotheses were tested by comparing ages of species and more inclusive lineages with their extant and most parsimonious ancestral distributions. Our data support a North American Paleocene origin for the Burseraceae followed by dispersal of ancestral lineages to eastern Laurasia and Southern Hemisphere continents.

Chloroplast haplotype diversity patterns inPackera pauciflora(Asteraceae) are affected by geographical isolation, hybridization, and breeding system

Packera pauciflora (Pursh) Löve & Löve is a wide ranging alpine or subalpine species that is disjunct in three regions of North America: the northern Rocky mountains, Quebec–Labrador, and the Sierra Nevada of California. It is one of two Packera species known to be self-fertile. We assessed intrapopulational chloroplast haplotype diversity (n) and haplotype nucleotide diversity ([Formula: see text]) and found it was lower and population structure (Φst), higher than all but one of the five other Packera species examined to date. In addition, distinct differences exist between populations from California and those from the Rocky mountains and eastern Canada. Our results suggest that geographical isolation, hybridization, and breeding system have all played a role in the formation of haplotype diversity patterns, and that, in spite of the wide disjunctions, the species is monophyletic.Key words: Asteraceae, cpDNA, haplotypes, Packera, phylogeography, self-fertility.

Two chloroplast DNA inversions originated simultaneously during the early evolution of the sunflower family (Asteraceae)

The chloroplast DNA (cpDNA) inversion in the Asteraceae has been cited as a classic example of using genomic rearrangements for defining major lineages of plants. We further characterize cpDNA inversions in the Asteraceae using extensive sequence comparisons among 56 species, including representatives of all major clades of the family and related families. We determine the boundaries of the 22-kb (now known as 22.8 kb) inversion that defines a major split within the Asteraceae, and in the process, we characterize the second and a new, smaller 3.3-kb inversion that occurs at one end of the larger inversion. One end point of the smaller inversion is upstream of the trnE-UUC gene, and the other end point is located between the trnC-GCA and rpoB genes. Although a diverse sampling of Asteraceae experienced substantial length variation and base substitution during the long evolutionary history subsequent to the inversion events, the precise locations of the inversion end points are identified using comparative sequence alignments in the inversion regions. The phylogenetic distribution of two inversions is identical among the members of Asteraceae, suggesting that the inversion events likely occurred simultaneously or within a short time period shortly after the origin of the family. Estimates of divergence times based on ndhF and rbcL sequences suggest that two inversions originated during the late Eocene (38-42 MYA). The divergence time estimates also suggest that the Asteraceae originated in the mid Eocene (42-47 MYA).

Methods for obtaining and analyzing whole chloroplast genome sequences

During the past decade, there has been a rapid increase in our understanding of plastid genome organization and evolution due to the availability of many new completely sequenced genomes. There are 45 complete genomes published and ongoing projects are likely to increase this sampling to nearly 200 genomes during the next 5 years. Several groups of researchers including ours have been developing new techniques for gathering and analyzing entire plastid genome sequences and details of these developments are summarized in this chapter. The most important developments that enhance our ability to generate whole chloroplast genome sequences involve the generation of pure fractions of chloroplast genomes by whole genome amplification using rolling circle amplification, cloning genomes into Fosmid or bacterial artificial chromosome (BAC) vectors, and the development of an organellar annotation program (Dual Organellar GenoMe Annotator [DOGMA]). In addition to providing details of these methods, we provide an overview of methods for analyzing complete plastid genome sequences for repeats and gene content, as well as approaches for using gene order and sequence data for phylogeny reconstruction. This explosive increase in the number of sequenced plastid genomes and improved computational tools will provide many insights into the evolution of these genomes and much new data for assessing relationships at deep nodes in plants and other photosynthetic organisms.

Molecular genetic evidence for interspecific hybridization among endemic Hispaniolan Bursera (Burseraceae)

Historically, genetic introgression among species as well as hybrid origins for species of the diploid tree genus Bursera (Burseraceae) have been proposed based on the supposition that individuals morphologically intermediate between sympatric "parent" species must be derived from hybridization. This study reports the first molecular genetic evidence for both unidirectional and reciprocal interspecific hybridization within Bursera. Phylogenies of hybrids and other species in B. subgenus Bursera are reconstructed based on nuclear and chloroplast sequence data. Compelling evidence supports the hybrid origin of three endemic Hispaniolan species: B. brunea (B. nashii × B. simaruba), B. gracilipes (B. spinescens × B. simaruba), and B. ovata (B. simaruba × B. spinescens). Cloning studies of nuclear markers from B. ovata suggests that this species is an introgressed or later backcross generation hybrid and thus reproduces sexually.

Molecular and cytological examination of Calopogon (Orchidaceae, Epidendroideae): circumscription, phylogeny, polyploidy, and possible hybrid speciation

The orchid genus Calopogon R.Br. (Orchidaceae), native to eastern North America and the northern Caribbean, currently contains five species and up to three varieties. Using nuclear internal transcribed spacer (ITS) ribosomal DNA sequences, amplified fragment length polymorphisms (AFLPs), chloroplast DNA restriction fragments, and chromosome counts, we present a phylogenetic and taxonomic study of the genus. Calopogon multiflorus and C. pallidus are consistently sister species, but the relationships of C. barbatus, C. oklahomensis, and C. tuberosus are not as clear. In the ITS analysis C. oklahomensis is sister to C. barbatus, whereas it is sister to C. tuberosus in the plastid restriction fragment analysis. Furthermore, all species were found to have chromosome numbers of 2n = 38 and 40, with the exception of the putatively hybrid-derived C. oklahomensis with 2n = 114 and 120. The hexaploidy of the latter, plus the discrepancy in its position between the ITS and plastid restriction fragment trees, could suggest that it is of hybrid origin. However, the presence of unique morphological and molecular characters might indicate that it is either an ancient hybrid or not of hybrid derivation at all. Finally, using these molecular methods all taxa appear to generally be discrete groups, with the exception of C. tuberosus vars. latifolius and tuberosus, the former of which is best combined with the latter.