Comparison of whole chloroplast genome sequences to choose noncoding regions for phylogenetic studies in angiosperms: the tortoise and the hare III

Chloroplast simple sequence repeats (cpSSRs): technical resources and recommendations for expanding cpSSR discovery and applications to a wide array of plant species

Chloroplast microsatellites, or simple sequence repeats (cpSSRs), are typically mononucleotide tandem repeats. When located in the noncoding regions of the chloroplast genome (cpDNA), they commonly show intraspecific variation in repeat number. Despite the growing number of studies applying cpSSRs, studies of economically important plants and their relatives remain over-represented. Thus, the potential of cpSSRs to offer unique insights into ecological and evolutionary processes in wild plant species has yet to be fully realized. This review provides an overview of the technical resources available to aid cpSSR discovery including a list of cpSSR primer sets available and cpDNA sequencing resources. Our updated analysis of 99 whole chloroplast genomes downloaded from GenBank confirms that potentially variable cpSSRs are abundant in the noncoding cpDNA of plants. Overall variation in the frequency of cpSSRs was extreme, ranging from one to 700 per genome (median = 93), while in 81 vascular plants, between 35 and 160 cpSSRs were detected per genome (median = 86). We offer five recommendations to aid wider development and application of cpSSRs: (i) When genus-specific cpSSR primers are available, cross-species amplification can often be fruitful. (ii) While potentially useful, universal cpSSR primers at best provide access to only a small number of variable markers. (iii) De novo sequencing of noncoding cpDNA is the most effective and efficient way to develop cpSSR markers in wild species. (iv) DNA sequencing of cpSSR alleles is essential, given the complex nature of the genetic variation associated with hypervariable cpDNA regions. (v) The reliability of cpSSR length based genetic assays need to be validated in all studies.

A new set of universal de novo sequencing primers for extensive coverage of noncoding chloroplast DNA: new opportunities for phylogenetic studies and cpSSR discovery

We present a new set of universal de novo sequencing primers targeting noncoding chloroplast DNA. The set of 107 polymerase chain reaction (PCR) primers span approximately 86% of the noncoding nucleotides in the large single copy region of Nicotiana tabacum, Oryza sativa and the orchid Phalaenopsis aphrodite. PCR tests confirmed the primers are effective in a wide range of monocots and dicots. More than 19.5 kb of cpDNA sequence was obtained across representative orchid genera with up to 82 chloroplast simple sequence repeats (cpSSRs) detected per genus. This primer set will facilitate both phylogenetic studies and rapid discovery of cpSSRs for plants, such as orchids, where there are limited genomic resources.

Multiplex sequencing of plant chloroplast genomes using Solexa sequencing-by-synthesis technology

Organellar DNA sequences are widely used in evolutionary and population genetic studies, however, the conservative nature of chloroplast gene and genome evolution often limits phylogenetic resolution and statistical power. To gain maximal access to the historical record contained within chloroplast genomes, we have adapted multiplex sequencing-by-synthesis (MSBS) to simultaneously sequence multiple genomes using the Illumina Genome Analyzer. We PCR-amplified approximately 120 kb plastomes from eight species (seven Pinus, one Picea) in 35 reactions. Pooled products were ligated to modified adapters that included 3 bp indexing tags and samples were multiplexed at four genomes per lane. Tagged microreads were assembled by de novo and reference-guided assembly methods, using previously published Pinus plastomes as surrogate references. Assemblies for these eight genomes are estimated at 88-94% complete, with an average sequence depth of 55x to 186x. Mononucleotide repeats interrupt contig assembly with increasing repeat length, and we estimate that the limit for their assembly is 16 bp. Comparisons to 37 kb of Sanger sequence show a validated error rate of 0.056%, and conspicuous errors are evident from the assembly process. This efficient sequencing approach yields high-quality draft genomes and should have immediate applicability to genomes with comparable complexity.

The D4 set: primers that target highly variable intron loops in plant chloroplast genomes

Chloroplast group II introns offer high-quality, rapidly evolving single-copy loci for comparative sequence analysis. These introns feature diagnostic secondary structures with loops that are among the least evolutionarily constrained sequence in plastomes. We exploited these structures to develop universal primers that amplify and sequence the large Domain IV (D4) loop in several angiosperm introns. With a single sequence read, we recover 300-600 nucleotides of highly variable sequence across angiosperms, with rates of change that are equal to or higher than many of the best known intergenic spacers in plant chloroplast genomes.

Evolution and phylogenetic utility of the PHOT gene duplicates in the Verbena complex (Verbenaceae): dramatic intron size variation and footprint of ancestral recombination

A well-resolved species level phylogeny is critically important in studying organismal evolution (e.g., hybridization, polyploidization, adaptive speciation). Lack of appropriate molecular markers that give sufficient resolution to gene trees is one of the major impediments to inferring species level phylogenies. In addition, sampling multiple independent loci is essential to overcome the lineage sorting problem. The availability of nuclear loci has often been a limiting factor in plant species-level phylogenetic studies. Here the two PHOT loci were developed as new sources of nuclear gene trees. The PHOT1 and PHOT2 gene trees of the Verbena complex (Verbenaceae) are well resolved and have good clade support. These gene trees are consistent with each other and previously generated chloroplast and nuclear waxy gene trees in most of the phylogenetic backbone as well as some terminal relationships, but are incongruent in some other relationships. Locus-specific primers were optimized for amplifying and sequencing these two loci in all Lamiales. Comparing intron size in the context of the gene trees shows dramatic variation within the Verbena complex, particularly at the PHOT1 locus. These variations are largely caused by invasions of short transposable elements and frequent long deletions and insertions of unknown causes. In addition, inspection of DNA sequences and phylogenetic analyses unmask a clear footprint of ancestral recombination in one species.

Intraspecific variation in Viola suavis in Europe: parallel evolution of white-flowered morphotypes

BACKGROUND AND AIMS

Viola species are commonly grown for their ornamental flowers, but their evolutionary history and taxonomy are often complicated and have been poorly explored so far. This is a study of the polymorphic, typically blue-flowered species Viola suavis, concentrating on the white-flowered populations of uncertain taxonomic assignment that occur in Spain and central and south-eastern Europe. The aim was to resolve their origin and taxonomic status and to study the intraspecific structure and (post)glacial history of this species.

METHODS

Viola suavis and five close relatives were sampled from multiple locations and subjected to molecular (AFLP, sequencing of nrDNA ITS) and morphometric analyses. Data on ploidy level and pollen fertility were also obtained, to address an assumed hybrid origin of the white-flowered populations.

KEY RESULTS

In V. suavis a strong intraspecific genetic split into two groups was observed, indicating that there has been a long-term isolation and survival in distinct glacial refugia. The white-flowered populations could be placed within the variation range of this species, and it is clear that they evolved independently in two distant areas. Their parallel evolution is supported by both morphological and genetic differentiation. The strongly reduced genetic variation and absence of unique AFLP fragments suggest their derived status and origin from the typical, blue-flowered populations.

CONCLUSIONS

These results suggest that intraspecific variation in V. suavis has been largely shaped by population isolations during the last glaciation and subsequent recolonizations, although cultivation and vegetative spread by humans have affected the present picture as well.

Phylogeny of Nolana (Solanaceae) of the Atacama and Peruvian deserts inferred from sequences of four plastid markers and the nuclear LEAFY second intron

The phylogeny of Nolana (Solanaceae), a genus primarily distributed in the coastal Atacama and Peruvian deserts with a few species in the Andes and one species endemic to the Galápagos Islands, was reconstructed using sequences of four plastid regions (ndhF, psbA-trnH, rps16-trnK and trnC-psbM) and the nuclear LEAFY second intron. The monophyly of Nolana was strongly supported by all molecular data. The LEAFY data suggested that the Chilean species, including Nolana sessiliflora, the N. acuminata group and at least some members of the Alona group, are basally diverged, supporting the Chilean origin of the genus. Three well-supported clades in the LEAFY tree were corroborated by the SINE (short interspersed elements) or SINE-like insertions. Taxa from Peru are grouped roughly into two clades. Nolana galapagensis from the Galápagos Island is most likely to have derived from a Peruvian ancestor. The monophyly of the morphologically well-diagnosed Nolana acuminata group (N. acuminata, N. baccata, N. paradoxa, N. parviflora, N. pterocarpa, N. rupicola and N. elegans) was supported by both plastid and LEAFY data. Incongruence between the plastid and the LEAFY data was detected concerning primarily the positions of N. sessiliflora, N. galapagensis, taxa of the Alona group and the two Peruvian clades. Such incongruence may be due to reticulate evolution or in some cases lineage sorting of plastid DNA. Incongruence between our previous GBSSI trees and the plastid-LEAFY trees was also detected concerning two well-supported major clades in the GBSSI tree. Duplication of the GBSSI gene may have contributed to this incongruence.

Multiple multilocus DNA barcodes from the plastid genome discriminate plant species equally well

A universal barcode system for land plants would be a valuable resource, with potential utility in fields as diverse as ecology, floristics, law enforcement and industry. However, the application of plant barcoding has been constrained by a lack of consensus regarding the most variable and technically practical DNA region(s). We compared eight candidate plant barcoding regions from the plastome and one from the mitochondrial genome for how well they discriminated the monophyly of 92 species in 32 diverse genera of land plants (N = 251 samples). The plastid markers comprise portions of five coding (rpoB, rpoC1, rbcL, matK and 23S rDNA) and three non-coding (trnH-psbA, atpF-atpH, and psbK-psbI) loci. Our survey included several taxonomically complex groups, and in all cases we examined multiple populations and species. The regions differed in their ability to discriminate species, and in ease of retrieval, in terms of amplification and sequencing success. Single locus resolution ranged from 7% (23S rDNA) to 59% (trnH-psbA) of species with well-supported monophyly. Sequence recovery rates were related primarily to amplification success (85-100% for plastid loci), with matK requiring the greatest effort to achieve reasonable recovery (88% using 10 primer pairs). Several loci (matK, psbK-psbI, trnH-psbA) were problematic for generating fully bidirectional sequences. Setting aside technical issues related to amplification and sequencing, combining the more variable plastid markers provided clear benefits for resolving species, although with diminishing returns, as all combinations assessed using four to seven regions had only marginally different success rates (69-71%; values that were approached by several two- and three-region combinations). This performance plateau may indicate fundamental upper limits on the precision of species discrimination that is possible with DNA barcoding systems that include moderate numbers of plastid markers. Resolution to the contentious debate on plant barcoding should therefore involve increased attention to practical issues related to the ease of sequence recovery, global alignability, and marker redundancy in multilocus plant DNA barcoding systems.

Phylogenetic utility of 141 low-copy nuclear regions in taxa at different taxonomic levels in two distantly related families of rosids

Angiosperm systematics has progressed to the point where it is now expected that multiple, independent markers be used in phylogenetic studies. Universal primers for amplifying informative regions of the chloroplast genome are readily available, but in the faster-evolving nuclear genome it is challenging to discover priming sites that are conserved across distantly related taxa. With goals including the identification of informative markers in rosids, and perhaps other angiosperms, we screened 141 nuclear primer combinations for phylogenetic utility in two distinct groups of rosids at different taxonomic levels-Psiguria (Cucurbitaceae) and Geraniaceae. We discovered three phylogenetically informative regions in Psiguria and two in Geraniaceae, but none that were useful in both groups. Extending beyond rosids, we combined our findings with those of another recent effort testing these primer pairs in Asteraceae, Brassicaceae, and Orchidaceae. From this comparison, we identified 32 primer combinations that amplified regions in representative species of at least two of the five distantly related angiosperm families, giving some prior indication about phylogenetic usefulness of these markers in other flowering plants. This reduced set of primer pairs for amplifying low-copy nuclear markers along with a recommended experimental strategy provide a framework for identifying phylogenetically informative regions in angiosperms.

Phylogeny of marattioid ferns (Marattiaceae): inferring a root in the absence of a closely related outgroup

Closely related outgroups are optimal for rooting phylogenetic trees; however, such ideal outgroups are not always available. A phylogeny of the marattioid ferns (Marattiaceae), an ancient lineage with no close relatives, was reconstructed using nucleotide sequences of multiple chloroplast regions (rps4 + rps4-trnS spacer, trnS-trnG spacer + trnG intron, rbcL, atpB), from 88 collections, selected to cover the broadest possible range of morphologies and geographic distributions within the extant taxa. Because marattioid ferns are phylogenetically isolated from other lineages, and internal branches are relatively short, rooting was problematic. Root placement was strongly affected by long-branch attraction under maximum parsimony and by model choice under maximum likelihood. A multifaceted approach to rooting was employed to isolate the sources of bias and produce a consensus root position. In a statistical comparison of all possible root positions with three different outgroups, most root positions were not significantly less optimal than the maximum likelihood root position, including the consensus root position. This phylogeny has several important taxonomic implications for marattioid ferns: Marattia in the broad sense is paraphyletic; the Hawaiian endemic Marattia douglasii is most closely related to tropical American taxa; and Angiopteris is monophyletic only if Archangiopteris and Macroglossum are included.

Reticulate or tree-like chloroplast DNA evolution in Sileneae (Caryophyllaceae)?

Despite sampling of up to 25kb of chloroplast DNA sequence from 24 species in Sileneae a number of nodes in the phylogeny remain poorly supported and it is not expected that additional sequence sampling will converge to a reliable phylogenetic hypothesis in these parts of the tree. The main reason for this is probably a combination of rapid radiation and substitution rate heterogeneity. Poor resolution among closely related species are often explained by low levels of variation in chloroplast data, but the problem with our data appear to be high levels of homoplasy. Tree-like cpDNA evolution cannot be rejected, but apparent incongruent patterns between different regions are evaluated with the possibility of ancient interspecific chloroplast recombination as explanatory model. However, several major phylogenetic relationships, previously not recognized, are confidently resolved, e.g. the grouping of the two SW Anatolian taxa S. cryptoneura and S. sordida strongly disagrees with previous studies on nuclear DNA sequence data, and indicate a possible case of homoploid hybrid origin. The closely related S. atocioides and S. aegyptiaca form a sister group to Lychnis and the rest of Silene, thus suggesting that Silene may be paraphyletic, despite recent revisions based on molecular data.

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.