Fifteen novel universal primer pairs for sequencing whole chloroplast genomes and a primer pair for nuclear ribosomal DNAs

Chloroplast genome information helps improve the phylogenetic resolution and can act as organelle-scale barcodes in recently radiated plant groups. Previously we reported that nine universal primer pairs could amplify angiosperm whole chloroplast genomes by long-range polymerase chain reaction and using next-generation sequencing. Although these primers show high universality and efficiency for sequencing whole chloroplast genomes in angiosperms, they did not fully resolve the following two issues surrounding sequencing angiosperm chloroplast genomes: (i) approximately 30% of angiosperms cannot be amplified successfully; and (ii) only fresh leaves can be applied. In this study, we designed another set of 15 universal primer pairs for amplifying angiosperm whole chloroplast genomes to complement the original nine primer pairs. Furthermore, we designed a primer pair for nuclear ribosomal DNAs (nrDNAs). To validate the functionality of the primers, we tested 44 species with silica gel-dried leaves and 15 species with fresh leaves that have been shown to not be amplified with the original nine primer pairs. The result showed that, in 65.9% and 88.6% of the 44 species with silica gel-dried leaves, the whole chloroplast genome and nrDNAs could be amplified, respectively. In addition, all 15 fresh leaf samples could have the whole chloroplast genome successfully amplified. The nrDNAs comprise partial sequences of 18S and 26S, along with the complete sequence of 5.8S and the internal transcribed spacers ITS1 and ITS2. The mean size of nrDNA was 5800 bp. This study shows that the 15 universal primer set is an indispensable tool for amplifying whole chloroplast genomes in angiosperms, and these are an important supplement to the nine reported primer pairs.     

Rapid sequencing of the bamboo mitochondrial genome using Illumina technology and parallel episodic evolution of organelle genomes in grasses

Background

Compared to their counterparts in animals, the mitochondrial (mt) genomes of angiosperms exhibit a number of unique features. However, unravelling their evolution is hindered by the few completed genomes, of which are essentially Sanger sequenced. While next-generation sequencing technologies have revolutionized chloroplast genome sequencing, they are just beginning to be applied to angiosperm mt genomes. Chloroplast genomes of grasses (Poaceae) have undergone episodic evolution and the evolutionary rate was suggested to be correlated between chloroplast and mt genomes in Poaceae. It is interesting to investigate whether correlated rate change also occurred in grass mt genomes as expected under lineage effects. A time-calibrated phylogenetic tree is needed to examine rate change.

Methodology/Principal Findings

We determined a largely completed mt genome from a bamboo, Ferrocalamus rimosivaginus (Poaceae), through Illumina sequencing of total DNA. With combination of de novo and reference-guided assembly, 39.5-fold coverage Illumina reads were finally assembled into scaffolds totalling 432,839 bp. The assembled genome contains nearly the same genes as the completed mt genomes in Poaceae. For examining evolutionary rate in grass mt genomes, we reconstructed a phylogenetic tree including 22 taxa based on 31 mt genes. The topology of the well-resolved tree was almost identical to that inferred from chloroplast genome with only minor difference. The inconsistency possibly derived from long branch attraction in mtDNA tree. By calculating absolute substitution rates, we found significant rate change (∼4-fold) in mt genome before and after the diversification of Poaceae both in synonymous and nonsynonymous terms. Furthermore, the rate change was correlated with that of chloroplast genomes in grasses.

Conclusions/Significance

Our result demonstrates that it is a rapid and efficient approach to obtain angiosperm mt genome sequences using Illumina sequencing technology. The parallel episodic evolution of mt and chloroplast genomes in grasses is consistent with lineage effects.     

被子植物叶绿体基因组的结构变异研究进展

叶绿体是绿色植物特有的细胞器,其基因组信息被广泛应用于植物系统发育和比较基因组学研究。目前,越来越多的物种有了叶绿体全基因组序列,人们对叶绿体基因组的结构及其变异规律有了更深入的了解。该文对近年来国内外有关被子植物叶绿体基因组插入/缺失、短片段倒位与重复、基因组结构重排以及基因丢失等结构变异式样的研究进展进行综述,并分析了叶绿体基因组结构研究中仍存在的问题以及该领域未来的发展趋势。     

Primers for 52 polymorphic regions in the Quercus rubra chloroplast, 47 of which amplify across 11 tracheophyte clades

Postglacial migration studies in Quercus rubra L. (northern red oak) are hampered by low levels of population differentiation in the widely used universal chloroplast primers. We sequenced the large single copy (LSC) regions of the Q. rubra and Quercus ellipsoidalis chloroplasts to enable us to query additional regions for future studies on migration and speciation. Using 454 sequencing of long-range PCR amplicons and Sanger sequencing for gap closure, we report 65 coding sequences from Q. rubra and 59 from Q. ellipsoidalis. Comparison of our de novo assembly of the LSC region sequence for Q. rubra to Q. rubra chloroplast sequence (NCBI Reference Sequence: NC_020152.1) from a different tree revealed 106 polymorphisms, all within intergenic regions, that can serve as tools for postglacial migration studies and taxonomic studies within the Lobatae. Sequence alignment for the 59 complete coding regions in common for theQ. rubrachloroplast reference sequence, our Q. rubra sequence and our Q. ellipsoidalis sequence revealed no sequence polymorphisms and no indels. We also report the 52 primer pairs we used for gap closure, including 53 new primer pairs not previously reported. We tested these 52 primer pairs against 11 species representing the Tracheophyta and detected 47 that produced amplicons in all 11 species. The new universal primers we have identified provide additional tools for resolving the taxonomic relationships among the congeneric taxa of forest trees in the temperate and subtropical forests of the Northern Hemisphere.     

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.     

Mitochondrial genes from 18 angiosperms fill sampling gaps for phylogenomic inferences of the early diversification of flowering plants

The early diversification of angiosperms is thought to have been a rapid process, which may complicate phylogenetic analyses of early angiosperm relationships. Plastid and nuclear phylogenomic studies have raised several conflicting hypotheses regarding overall angiosperm phylogeny, but mitochondrial genomes have been largely ignored as a relevant source of information. Here we sequenced mitochondrial genomes from 18 angiosperms to fill taxon-sampling gaps in Austrobaileyales, magnoliids, Chloranthales, Ceratophyllales, and major lineages of eudicots and monocots. We assembled a data matrix of 38 mitochondrial genes from 107 taxa to assess how well mitochondrial genomic data address current uncertainties in angiosperm relationships. Although we recovered conflicting phylogenies based on different data sets and analytical methods, we also observed congruence regarding deep relationships of several major angiosperm lineages: Chloranthales were always inferred to be the sister group of Ceratophyllales, Austrobaileyales to mesangiosperms, and the unplaced Dilleniales was consistently resolved as the sister to superasterids. Substitutional saturation, GC compositional heterogeneity, and codon-usage bias are possible reasons for the noise/conflict that may impact phylogenetic reconstruction; and angiosperm mitochondrial genes may not be substantially affected by these factors. The third codon positions of the mitochondrial genes appear to contain more parsimony-informative sites than the first and second codon positions, and therefore produced better resolved phylogenetic relationships with generally strong support. The relationships among these major lineages remain incompletely resolved, perhaps as a result of the rapidity of early radiations. Nevertheless, data from mitochondrial genomes provide additional evidence and alternative hypotheses for exploring the early evolution and diversification of the angiosperms.     

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.     

New universal matK primers for DNA barcoding angiosperms

The chloroplast maturase K gene (matK) is one of the most variable coding genes of angiosperms and has been suggested to be a "barcode" for land plants. However, matK exhibits low amplification and sequencing rates due to low universality of currently available primers and mononucleotide repeats. To resolve these technical problems, we evaluated the entire matK region to find a region of 600-800 bp that is highly variable, represents the best of all matK regions with priming sites conservative enough to design universal primers, and avoids the mononucleotide repeats. After careful evaluation, a region in the middle was chosen and a pair of primers named natK472F and matK1248R was designed to amplify and sequence the matK fragment of approximately 776 bp. This region encompasses the most variable sites, represents the entire matK region best, and also exhibits high amplification rates and quality of sequences. The universality of this primer pair was tested using 58 species from 47 families of angiosperm plants. The primers showed a strong amplification (93.1%) and sequencing (92.6%)successes in the species tested. We propose that the new primers will solve, in part, the problems encountered when using matK and promote the adoption of matK as a DNA barcode for angiosperms.     

Molecular evolution and phylogenetic utility of the petD group II intron: a case study in basal angiosperms

Sequences of spacers and group I introns in plant chloroplast genomes have recently been shown to be very effective in phylogenetic reconstruction at higher taxonomic levels and not only for inferring relationships among species. Group II introns, being more frequent in those genomes than group I introns, may be further promising markers. Because group II introns are structurally constrained, we assumed that sequences of a group II intron should be alignable across seed plants. We designed universal amplification primers for the petD intron and sequenced this intron in a representative selection of 47 angiosperms and three gymnosperms. Our sampling of taxa is the most representative of major seed plant lineages to date for group II introns. Through differential analysis of structural partitions, we studied patterns of molecular evolution and their contribution to phylogenetic signal. Nonpairing stretches (loops, bulges, and interhelical nucleotides) were considerably more variable in both substitutions and indels than in helical elements. Differences among the domains are basically a function of their structural composition. After the exclusion of four mutational hotspots accounting for less than 18% of sequence length, which are located in loops of domains I and IV, all sequences could be aligned unambiguously across seed plants. Microstructural changes predominantly occurred in loop regions and are mostly simple sequence repeats. An indel matrix comprising 241 characters revealed microstructural changes to be of lower homoplasy than are substitutions. In showing Amborella first branching and providing support for a magnoliid clade through a synapomorphic indel, the petD data set proved effective in testing between alternative hypotheses on the basal nodes of the angiosperm tree. Within angiosperms, group II introns offer phylogenetic signal that is intermediate in information content between that of spacers and group I introns on the one hand and coding sequences on the other.     

Characterizing homologues of crop domestication genes in poorly described wild relatives by high‐throughput sequencing of whole genomes

Wild crop relatives represent a source of novel alleles for crop genetic improvement. Screening biodiversity for useful or diverse gene homologues has often been based upon the amplification of targeted genes using available sequence information to design primers that amplify the target gene region across species. The crucial requirement of this approach is the presence of sequences with sufficient conservation across species to allow for the design of universal primers. This approach is often not successful with diverse organisms or highly variable genes. Massively parallel sequencing (MPS) can quickly produce large amounts of sequence data and provides a viable option for characterizing homologues of known genes in poorly described genomes. MPS of genomic DNA was used to obtain species‐specific sequence information for 18 rice genes related to domestication characteristics in a wild relative of rice, Microlaena stipoides. Species‐specific primers were available for 16 genes compared with 12 genes using the universal primer method. The use of species‐specific primers had the potential to cover 92% of the sequence of these genes, while traditional universal primers could only be designed to cover 80%. A total of 24 species‐specific primer pairs were used to amplify gene homologues, and 11 primer pairs were successful in capturing six gene homologues. The 23 million, 36‐base pair (bp) paired end reads, equated to an average of 2X genome coverage, facilitated the successful amplification and sequencing of six target gene homologues, illustrating an important approach to the discovery of useful genes in wild crop relatives.     

Microfluidic PCR-based target enrichment: A case study in two rapid radiations of Commiphora (Burseraceae) from Madagascar

Developing effective and cost-efficient multilocus nuclear datasets for angiosperm species is a continuing challenge to the systematics community. Here we describe the development and validation of a novel set of 91 nuclear markers for PCR-based target enrichment. Using microfluidic PCR and Illumina MiSeq, we generated nuclear, subgenomic libraries for 96 species simultaneously and sequenced them for a total cost of ca. $6000 USD. Approximately half of these costs include reusable reagents (primers, barcodes, and custom sequencing primers) and taxon sampling could be increased by an order of magnitude to maximize sequencing depth efficiency. The principle benefit of microfluidic PCR over alternative target enrichment strategies is that it bypasses costly library preparation. After sequencing, we evaluated the ability of the loci to resolve species level relationships within two recently radiated lineages of endemic Madagascan Commiphora Jacq. (Burseraceae) species. Our results demonstrate that (i) effective nuclear markers can be designed for non-model angiosperm taxa from these publicly available datasets; (ii) that microfluidic PCR amplification followed by high throughput sequencing can produce highly complete taxon by locus sequence data matrices with minimal resource investment; and (iii) that these numerous nuclear phylogenomic markers can improve our understanding of phylogenetic relationships withinCommiphora. We provide a synopsis of ongoing activities to enhance this microfluidic PCR-based target enrichment strategy through broader primer assays, multiplexing, and increased efficiency of sequencing depth.     

Interspecific polymorphism at non-coding regions of chloroplast,mitochondrial DNA and rRNA IGS region in Elymus species

Several published universal primers for amplification of non-coding regions of chloroplast, mitochondrial and ribosomal (rRNA) IGS region were tested whether they can amplify respective regions in Elymus species. PCR-RFLP analysis of the chloroplast, mitochondral DNA, and rRNA IGS region of the genus Elymus was used to determine if the method could be employed to detect inter-specific variation in this genus. Published universal primers for amplification of trnK [tRNA-Lys (UUU) exon 1]-trnK [tRNA-Lys (UUU) exon2], and mitochondrial nad1 exon B-nadl exon C intron successfully amplified the respective regions in Elymus species. However, the primers for amplification of chloroplast trnD-trnT intron and rRNA IGS failed to amplify the respective region in Elymus species. New primer pairs were designed and successfully amplified the cpDNA trnD-trnT intron and rRNA IGS region in Elymus species. The amplification products were digested with seven restriction enzymes. The results showed that the investigated regions of chloroplast and mitochondrial genomes are variable in most of the tested taxa and contain multiple variable regions. These regions should serve as useful molecular markers in phylogenetic studies of closely related species, at least at the interspecific level in Elymus. It is likely that further studies, including larger sample sizes, more regions of these genomes and/or more powerful methods for the detection of cpDNA and mt DNA variation will reveal additional variation for this genus. Highly inter- and intra-specific polymorphisms for rRNA IGS region were detected, suggesting the IGS will be a useful molecular marker for population studies of Elymus species.     

Comparative Chloroplast Genomes of Camellia Species

Background

Camellia , comprising more than 200 species, is a valuable economic commodity due to its enormously popular commercial products: tea leaves, flowers, and high-quality edible oils. It is the largest and most important genus in the family Theaceae. However, phylogenetic resolution of the species has proven to be difficult. Consequently, the interspecies relationships of the genus Camellia are still hotly debated. Phylogenomics is an attractive avenue that can be used to reconstruct the tree of life, especially at low taxonomic levels.

Methodology/Principal Findings

Seven complete chloroplast (cp) genomes were sequenced from six species representing different subdivisions of the genus Camellia using Illumina sequencing technology. Four junctions between the single-copy segments and the inverted repeats were confirmed and genome assemblies were validated by PCR-based product sequencing using 123 pairs of primers covering preliminary cp genome assemblies. The length of the Camellia cp genome was found to be about 157kb, which contained 123 unique genes and 23 were duplicated in the IR regions. We determined that the complete Camellia cp genome was relatively well conserved, but contained enough genetic differences to provide useful phylogenetic information. Phylogenetic relationships were analyzed using seven complete cp genomes of six Camellia species. We also identified rapidly evolving regions of the cp genome that have the potential to be used for further species identification and phylogenetic resolution.

Conclusions/Significance

In this study, we wanted to determine if analyzing completely sequenced cp genomes could help settle these controversies of interspecies relationships in Camellia . The results demonstrate that cp genome data are beneficial in resolving species definition because they indicate that organelle-based “barcodes”, can be established for a species and then used to unmask interspecies phylogenetic relationships. It reveals that phylogenomics based on cp genomes is an effective approach for achieving phylogenetic resolution between Camellia species.     

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.     

Development and testing of novel chloroplast microsatellite markers for Lolium perenne and other grasses (Poaceae) from de novo sequencing and in silico sequences

Twelve primers to amplify microsatellite markers from the chloroplast genome of Lolium perenne were designed and optimized using de novo sequencing and in silico sequences. With one exception, each locus was polymorphic with a range from two to nine alleles in L. perenne. The newly developed primer pairs cross‐amplified in different species of Lolium and in 50 other grass species representing nine grass subfamilies.     

Complete chloroplast genome sequence of Elodea canadensis and comparative analyses with other monocot plastid genomes

Elodea canadensis is an aquatic angiosperm native to North America. It has attracted great attention due to its invasive nature when transported to new areas in its non-native range. We have determined the complete nucleotide sequence of the chloroplast (cp) genome of Elodea. Taxonomically Elodea is a basal monocot, and only few monocot cp genomes representing early lineages of monocots have been sequenced so far. The genome is a circular double-stranded DNA molecule 156,700bp in length, and has a typical structure with large (LSC 86,194bp) and small (SSC 17,810bp) single-copy regions separated by a pair of inverted repeats (IRs 26,348bp each). The Elodea cp genome contains 113 unique genes and 16 duplicated genes in the IR regions. A comparative analysis showed that the gene order and organization of the Elodea cp genome is almost identical to that of Amborella trichopoda, a basal angiosperm. The structure of IRs in Elodea is unique among monocot species with the whole cp genome sequenced. In Elodea and another monocot Lemna minor the borders between IRs and LSC are located upstream of rps19 gene and downstream of trnH-GUG gene, while in most monocots, IR has extended to include both trnH and rps19 genes. A phylogenetic analysis conducted using Bayesian method, based on the DNA sequences of 81 chloroplast genes from 17 monocot taxa provided support for the placement of Elodea together with Lemna as a basal monocot and the next diverging lineage of monocots after Acorales. In comparison with other monocots, the Elodea cp genome has gone through only few rearrangements or gene losses. IR of Elodea has a unique structure among the monocot species studied so far as its structure is similar to that of a basal angiosperm Amborella. This result together with phylogenetic analyses supports the placement of Elodea as a basal monocot to the next diverging lineage of monocots after Acorales. So far, only few cp genomes representing early lineages of monocots have been sequenced and, therefore, this study provides valuable information about the course of evolution in divergence of monocot lineages.     

Using phylogenomics to resolve mega-families: An example from Compositae

Next-generation sequencing and phylogenomics hold great promise for elucidating complex relationships among large plant families. Here, we performed targeted capture of low copy sequences followed by next-generation sequencing on the Illumina platform in the large and diverse angiosperm family Compositae (Asteraceae). The family is monophyletic, based on morphology and molecular data, yet many areas of the phylogeny have unresolved polytomies and interpreting phylogenetic patterns has been historically difficult. In order to outline a method and provide a framework and for future phylogenetic studies in the Compositae, we sequenced 23 taxa from across the family in which the relationships were well established as well as a member of the sister family Calyceraceae. We generated nuclear data from 795 loci and assembled chloroplast genomes from off-target capture reads enabling the comparison of nuclear and chloroplast genomes for phylogenetic analyses. We also analyzed multi-copy nuclear genes in our data set using a clustering method during orthology detection, and we applied a network approach to these clusters—analyzing all related locus copies. Using these data, we produced hypotheses of phylogenetic relationships employing both a conservative (restricted to only loci with one copy per targeted locus) and a multigene approach (including all copies per targeted locus). The methods and bioinformatics workflow presented here provide a solid foundation for future work aimed at understanding gene family evolution in the Compositae as well as providing a model for phylogenomic analyses in other plant mega-families.     

Complete plastid genome sequences of Drimys, Liriodendron, and Piper : implications for the phylogenetic relationships of magnoliids

Background  

The magnoliids with four orders, 19 families, and 8,500 species represent one of the largest clades of early diverging angiosperms. Although several recent angiosperm phylogenetic analyses supported the monophyly of magnoliids and suggested relationships among the orders, the limited number of genes examined resulted in only weak support, and these issues remain controversial. Furthermore, considerable incongruence resulted in phylogenetic reconstructions supporting three different sets of relationships among magnoliids and the two large angiosperm clades, monocots and eudicots. We sequenced the plastid genomes of three magnoliids, Drimys (Canellales), Liriodendron (Magnoliales), and Piper (Piperales), and used these data in combination with 32 other angiosperm plastid genomes to assess phylogenetic relationships among magnoliids and to examine patterns of variation of GC content.