Highly effective sequencing whole chloroplast genomes of angiosperms by nine novel universal primer pairs

Chloroplast genomes supply indispensable information that helps improve the phylogenetic resolution and even as organelle‐scale barcodes. Next‐generation sequencing technologies have helped promote sequencing of complete chloroplast genomes, but compared with the number of angiosperms, relatively few chloroplast genomes have been sequenced. There are two major reasons for the paucity of completely sequenced chloroplast genomes: (i) massive amounts of fresh leaves are needed for chloroplast sequencing and (ii) there are considerable gaps in the sequenced chloroplast genomes of many plants because of the difficulty of isolating high‐quality chloroplast DNA, preventing complete chloroplast genomes from being assembled. To overcome these obstacles, all known angiosperm chloroplast genomes available to date were analysed, and then we designed nine universal primer pairs corresponding to the highly conserved regions. Using these primers, angiosperm whole chloroplast genomes can be amplified using long‐range PCR and sequenced using next‐generation sequencing methods. The primers showed high universality, which was tested using 24 species representing major clades of angiosperms. To validate the functionality of the primers, eight species representing major groups of angiosperms, that is, early‐diverging angiosperms, magnoliids, monocots, Saxifragales, fabids, malvids and asterids, were sequenced and assembled their complete chloroplast genomes. In our trials, only 100 mg of fresh leaves was used. The results show that the universal primer set provided an easy, effective and feasible approach for sequencing whole chloroplast genomes in angiosperms. The designed universal primer pairs provide a possibility to accelerate genome‐scale data acquisition and will therefore magnify the phylogenetic resolution and species identification in angiosperms.     

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.     

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.     

Design and testing of a plant-specific PCR primer for ecological and evolutionary studies

A polymerase chain reaction (PCR) primer, 28KJ (5-GGCGGTAAATTCCGTCC-3), was developed to specifically amplify plant DNA. This primer is located approximately 250 bases downstream of the 5′ end of the 28S ribosomal RNA gene, and it was used in combination with the universal primer 28C. The specificity of this primer combination was tested against 31 angiosperms, 9 conifers, 1 alga and 30 fungi (21 basidiomycetes and 9 ascomycetes). Both herbarium specimens and fresh samples were tested. The 28KJ/28C primer combination successfully amplified all angiosperm and conifer DNAs, but no fungal or algal DNAs. Plant DNA was amplified from plant/fungal symbioses (ectomycorrhizae of conifers and ericoid mycorrhizae of Ericaceous plants), and the plants involved in these symbioses were identified by comparing DNA sequences or restriction enzyme digest patterns of the mycorrhizal DNAs to those of known plant samples. These methods allow rapid and accurate identification of plant associates in complex plant/ fungal systems when the identity of the roots is unclear.     

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.     

16S/18S/ITS扩增子高通量测序引物的生物信息学评估和改进

【背景】在过去的十几年里,基于核糖体RNA基因的扩增子测序技术被广泛用于各种生态系统中微生物群落的多样性检测。扩增子测序的使用极大地促进了土壤、水体、空气等环境中微生物生态的相关研究。【目的】随着高通量测序技术的不断发展和参考数据库的不断更新,针对不同的环境样本的引物选择和改进仍然需要更深入的校验。【方法】本文收集了目前在微生物群落研究中被广泛采用的标记基因扩增通用引物,包括16S rRNA基因扩增常用的8对通用引物和2对古菌引物、9对真菌转录间隔区(internal transcribed spacer,ITS)基因扩增引物,以及18S rRNA基因扩增的4对真核微生物通用引物和1对真菌特异性引物。这些引物中包括了地球微生物组计划(Earth Microbiome Project,EMP)推荐的2对16S rRNA基因扩增引物、1对ITS1基因扩增引物和1对18S rRNA基因扩增引物。采用最近更新的标准数据库对这些引物进行了覆盖度和特异性评价。【结果】EMP推荐的引物依然具有较高的覆盖度,而其他引物在覆盖度及对特定环境或类群的特异性上也各有特点。此外,最近有研究对这些通用引物进行了一些改进,而我们也发现,一个碱基的变化都可能会导致评价结果或扩增产物发生明显变化,简并碱基的引入既可以覆盖更多的物种,但同时也会在一定程度上降低关注物种的特异性。【结论】研究结果为微生态研究中标记基因的引物选择提供了一个广泛的指导,但在关注具体科学问题时,引物的选择仍需数据指导与实验尝试。     

The development and evaluation of consensus chloroplast primer pairs that possess highly variable sequence regions in a diverse array of plant taxa

Although universal or consensus chloroplast primers are available, they are limited by their number and genomic distribution. Therefore, a set of consensus chloroplast primer pairs for simple sequence repeats (ccSSRs) analysis was constructed from tobacco (Nicotiana tabacum L.) chloroplast sequences. These were then tested for their general utility in the genetic analysis of a diverse array of plant taxa. In order to increase the number of ccSSRs beyond that previously reported, the target sequences for SSR motifs was set at A or T (n 7) mononucleotide repeats. Each SSR sequence motif, along with ±200-bp flanking sequences from the first of each mononucleotide base repeat, was screened for homologies with chloroplast DNA sequences of other plant species in GenBank databases using BLAST search procedures. Twenty three putative marker loci that possessed conserved flanking sequence spans were selected for consensus primer pair construction using commercially available computer algorithms. All primer pairs produced amplicons after PCR employing genomic DNA from members of the Cucurbitaceae (six species) and Solanaceae (four species). Sixteen, 22 and 19 of the initial 23 primer pairs were successively amplified by PCR using template DNA from species of the Apiaceae (two species), Brassicaceae (one species) and Fabaceae (two species), respectively. Twenty of 23 primer pairs were also functional in three monocot species of the Liliaceae [onion (Allium cepa L.) and garlic (Allium sativum L.)], and the Poaceae [oat (Avena sativa L.)]. Sequence analysis of selected ccSSR fragments suggests that ccSSR length and sequence variation could be useful as a tool for investigating the genetic relationships within a genus or closely related taxa (i.e., tribal level). In order to provide for a marker system having significant coverage of the cucumber chloroplast genome, ccSSR primers were strategically "recombined" and named recombined consensus chloroplast primers (RCCP) for PCR analysis. Successful amplification after extended-length PCR of 16 RCCP primer pairs from cucumber (Cucumis sativus L.) DNA suggested that the amplicons detected are representative of the cucumber chloroplast genome. These RCCP pairs, therefore, could be useful as an initial molecular tool for investigation of traits related to a chloroplast gene(s) in cucumber, and other closely related species.Communicated by C. Möllers     

Assessment of fungal diversity in deep-sea sediments by multiple primer approach

Increasing evidence of the fungal diversity in deep-sea sediments has come from amplification of environmental DNA with fungal specific or eukaryote primer sets. In order to assess the fungal diversity in deep-sea sediments of the Central Indian Basin (CIB) at ~5,000 m depth, we amplified sediment DNA with four different primer sets. These were fungal-specific primer pair ITS1F/ITS4 (internal transcribed spacers), universal 18S rDNA primers NS1/NS2, Euk18S-42F/Euk18S-1492R and Euk18S-555F/Euk18S-1269R. One environmental library was constructed with each of the primer pairs, and 48 clones were sequenced per library. These sequences resulted in 8 fungal Operational Taxonomic Units (OTUs) with ITS and 19 OTUs with 18S rDNA primer sets respectively by taking into account the 2% sequence divergence cut-off for species delineation. These OTUs belonged to 20 distinct fungal genera of the phyla Ascomycota and Basidiomycota. Seven sequences were found to be divergent by 79–97% from the known sequences of the existing database and may be novel. A majority of the sequences clustered with known sequences of the existing taxa. The phylogenetic affiliation of a few fungal sequences with known environmental sequences from marine and hypersaline habitat suggests their autochthonous nature or adaptation to marine habitat. The amplification of sequences belonging to Exobasidiomycetes and Cystobasidiomycetes from deep-sea is being reported for the first time in this study. Amplification of fungal sequences with eukaryotic as well as fungal specific primers indicates that among eukaryotes, fungi appear to be a dominant group in the sampling site of the CIB.     

Barcoding old Korean lepidopteran specimens using newly designed specific primer pairs

Currently, DNA barcodes are often required to be analyzed using old museum specimens when they are the only available specimens for rare or endangered species, or even type series. In this study, using eight universal primers and newly designed 315 species-specific primers, we tried to recover full-length barcode sequences from 45 dried specimens of 36 butterfly species collected between 1959 and 1980 in Korea. The eight universal primers failed entirely in the PCR amplification and sequencing of all the specimens. On the other hand, 284 primer pairs consisting of the 315 primers, targeting fragments of 71–417 bp, amplified various lengths of barcode sequences from all specimens. The fragments were successfully combined to generate the barcode sequences ranging from 444 bp to 658 bp. Notably, of the 284 primer pairs, 26 primer pairs designed for Limenitis camilla, Argynnis niobe, and Brenthis daphne successfully amplified the barcode sequences of congeneric species, Limenitis doerriesi, Argynnis nerippe, and Brenthis ino, suggesting that the species-specific primers can be available for analyzing barcode sequences of closely related species. Our study reveals that the newly designed species-specific primers will be effective in acquiring COI sequences from old butterfly specimens.     

Barcoding the kingdom Plantae: new PCR primers for ITS regions of plants with improved universality and specificity

The internal transcribed spacer (ITS) of nuclear ribosomal DNA is one of the most commonly used DNA markers in plant phylogenetic and DNA barcoding analyses, and it has been recommended as a core plant DNA barcode. Despite this popularity, the universality and specificity of PCR primers for the ITS region are not satisfactory, resulting in amplification and sequencing difficulties. By thoroughly surveying and analysing the 18S, 5.8S and 26S sequences of Plantae and Fungi from GenBank, we designed new universal and plant‐specific PCR primers for amplifying the whole ITS region and a part of it (ITS1 or ITS2) of plants. In silico analyses of the new and the existing ITS primers based on these highly representative data sets indicated that (i) the newly designed universal primers are suitable for over 95% of plants in most groups; and (ii) the plant‐specific primers are suitable for over 85% of plants in most groups without amplification of fungi. A total of 335 samples from 219 angiosperm families, 11 gymnosperm families, 24 fern and lycophyte families, 16 moss families and 17 fungus families were used to test the performances of these primers. In vitro PCR produced similar results to those from the in silico analyses. Our new primer pairs gave PCR improvements up to 30% compared with common‐used ones. The new universal ITS primers will find wide application in both plant and fungal biology, and the new plant‐specific ITS primers will, by eliminating PCR amplification of nonplant templates, significantly improve the quality of ITS sequence information collections in plant molecular systematics and DNA barcoding.     

双翅目昆虫线粒体基因组结构特点及其测序通用引物的设计和应用

研究双翅目昆虫线粒体基因组的结构特点, 并设计其测序的通用引物, 为今后双翅目昆虫线粒体基因组的研究提供参考和依据。利用比较基因组学和生物信息学方法, 分析了已经完全测序的26个双翅目昆虫线粒体基因组的结构特点、 碱基组成和保守区, 并据此设计了双翅目昆虫基因组测序的通用引物。结果表明： 双翅目昆虫线粒体基因组长14 503~19 517 bp, 其结构保守, 含有37个编码基因, 包括13个蛋白质编码基因, 22个tRNA编码基因和2个rRNA编码基因, 此外还包含一段长度差异很大的非编码区(AT富含区)。基因组内基因排列次序稳定, 除个别基因外, 其余都与黑腹果蝇Drosophila melanogaster基因排列次序一致。基因组的碱基组成不均衡, AT含量在72.59%~85.15%之间, 碱基使用存在偏向性, 偏好使用AC碱基。全基因组的核苷酸和氨基酸序列保守, 共鉴定了11个保守区。在保守区内共设计了26对双翅目线粒体基因组测序通用引物, 扩增的目标片段都在1 200 bp以内。将该套通用引物用于葱蝇Delia antiqua线粒体全基因组测序, 结果证明其高效、 合用。     

Transferability and utility of white oat (Avena sativa) microsatellite markers for genetic studies in black oat (Avena strigosa)

Preservation and use of wild oat species germplasm are essential for further improvement of cultivated oats. We analyzed the transferability and utility of cultivated (white) oat Avena sativa (AACCDD genome) microsatellite markers for genetic studies of black oat A. strigosa (A(s)A(s) genome) genotypes. The DNA of each black oat genotype was extracted from young leaves and amplified by PCR using 24 microsatellite primers developed from white oat. The PCR products were separated on 3% agarose gel. Eighteen microsatellite primer pairs amplified consistent products and 15 of these were polymorphic in A. strigosa, demonstrating a high degree of transferability. Microsatellite primer pairs AM3, AM4, AM21, AM23, AM30, and AM35 consistently amplified alleles only in A. sativa, which indicates that they are putative loci for either the C or D genomes of Avena. Using the data generated by the 15 polymorphic primer pairs, it was possible to separate 40 genotypes of the 44 that we studied. The four genotypes that could not be separated are probably replicates. We conclude that A. sativa microsatellites have a high transferability index and are a valuable resource for genetic studies and characterization of A. strigosa genotypes.     

Development of primer sets designed for use with the PCR to amplify conserved genes from filamentous ascomycetes.

We constructed nine sets of oligonucleotide primers on the basis of the results of DNA hybridization of cloned genes from Neurospora crassa and Aspergillus nidulans to the genomes of select filamentous ascomycetes and deuteromycetes (with filamentous ascomycete affiliations). Nine sets of primers were designed to amplify segments of DNA that span one or more introns in conserved genes. PCR DNA amplification with the nine primer sets with genomic DNA from ascomycetes, deuteromycetes, basidiomycetes, and plants revealed that five of the primer sets amplified a product only from DNA of the filamentous ascomycetes and deuteromycetes. The five primer sets were constructed from the N. crassa genes for histone 3, histone 4, beta-tubulin, and the plasma membrane ATPase. With these five primer sets, polymorphisms were observed in both the size of and restriction enzyme sites in the amplified products from the filamentous ascomycetes. The primer sets described here may provide useful tools for phylogenetic studies and genome analyses in filamentous ascomycetes and deuteromycetes (with ascomycete affiliations), as well as for the rapid differentiation of fungal species by PCR.     

New specific primers for amplification of the Internal Transcribed Spacer region in Clitellata (Annelida)

Nuclear molecular evidence, for example, the rapidly evolving Internal Transcribed Spacer region (ITS), integrated with maternally inherited (mitochondrial) COI barcodes, has provided new insights into the diversity of clitellate annelids. PCR amplification and sequencing of ITS, however, are often hampered by poor specificity of primers used. Therefore, new clitellate‐specific primers for amplifying the whole ITS region (ITS: 29F/1084R) and a part of it (ITS2: 606F/1082R) were developed on the basis of a collection of previously published ITS sequences with flanking rDNA coding regions. The specificity of these and other ITS primers used for clitellates were then tested in silico by evaluating their mismatches with all assembled and annotated sequences (STD, version r127) from EMBL, and the new primers were also tested in vitro for a taxonomically broad sample of clitellate species (71 specimens representing 11 families). The in silico analyses showed that the newly designed primers have a better performance than the universal ones when amplifying clitellate ITS sequences. In vitro PCR and sequencing using the new primers were successful, in particular, for the 606F/1082R pair, which worked well for 65 of the 71 specimens. Thus, using this pair for amplifying the ITS2 will facilitate further molecular systematic investigation of various clitellates. The other pair (29F/1084R), will be a useful complement to existing ITS primers, when amplifying ITS as a whole.     

Optimization of primer screening for evaluation of genetic relationship in rose cultivars

Optimization of primer screening for evaluation of genetic relationship in 34 cultivars of rose through random amplified polymorphic DNA (RAPD) markers was investigated. Four series of decamer primers were used for screening and optimization of RAPD analysis between which A and N series performed good amplification of fragments as compared with other series. The primers OPN-07 and OPN-15 produced maximum number of DNA fragments in Rosa hybrida cv. Anuraag. Some primer either did not produce amplification or produced very poor amplification. Further, ten selected primers were used for genetic analysis of 34 rose cultivars. The primer OPN-15 amplified 21 fragments in all cultivars tested. A total of 162 distinct DNA fragments (bands) ranging from 100 to 3400 base pairs were amplified by using 10 selected random primers. The cluster analysis indicated that these rose cultivars formed nine clusters.     

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.     

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.     

藻类植物叶绿体基因组研究进展

藻类植物的cpDNA结构复杂,普遍缺失反向重复序列IR,且存在IR的藻类植物种类的cpDNA也有IR变短退化迹象.藻类植物的cpDNA包含的基因一般比高等植物要多,编码能力更强.藻类植物cpDNA全序列的测定方法主要是Fosmid文库构建,配合使用Long-PCR技术.该文对国内外有关藻类植物叶绿体基因组结构、叶绿体编码基因、叶绿体基因组在藻类系统发育中的应用以及藻类植物叶绿体基因组的提取和序列测定方法等进行综述,为藻类植物的系统发育和叶绿体起源以及功能基因组学的研究提供理论依据.     

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.