A set of 35 consensus primer pairs amplifying genes and introns of plant mitochondrial DNA

Being generally uniparentally inherited, the plant mitochondrial genome is a source of original markers potentially useful for studies of phylogeny and population genetics. We designed 24 new pairs of plant mitochondrial DNA primers that allow amplification of either introns, intergenic regions or genes. They have been defined for consensus over angiosperms and were tested along with 11 previously described mitochondrial primer pairs on 28 plant species representing 19 families of higher plants. The total set allows amplification of 40 kb (～11%) of the mtDNA genome of Arabidopsis thaliana. The amplification rate ranged between 76% and 100% depending on the species.     

In silico mining for simple sequence repeat loci in a pineapple expressed sequence tag database and cross-species amplification of EST-SSR markers across Bromeliaceae

A collection of 5,659 expressed sequence tags (ESTs) from pineapple [Ananas comosus (L.) Merr.] was screened for simple sequence repeats (EST-SSRs) with motif lengths between 1 and 6 bp. Lower thresholds of 15, 7 and 5 repeat units were used to define microsatellites of the mono-, di-, and tri- to hexanucleotide repeat type, respectively. Based on these criteria, 696 SSRs were identified among 3,389 EST unigenes, together representing 2,840 kb. This corresponds to an average density of one SSR every 4.1 kb of non-redundant EST sequences. Dinucleotide repeats were most abundant (38.4% of all SSRs) followed by trinucleotide repeats (38.1%). Flanking primer pairs were designed for 537 EST-SSR loci, and 49 of these were screened for their functionality in 12 accessions of A. comosus, 14 accessions of 5 additional Ananas species and 1 species of Pseudananas. Distinct PCR products of the expected size range were obtained with 36 primer pairs. Eighteen loci analyzed in more detail were all polymorphic in pineapple, and primer pairs flanking these loci also generated PCR products from a wide range of genera and species from six subfamilies of the Bromeliaceae. The potential to reveal polymorphism in a heterologous target species was demonstrated in Deuterocohnia brevifolia (subfamily Pitcairnioideae).     

High‐throughput multiplex cpDNA resequencing clarifies the genetic diversity and genetic relationships among Brassica napus,Brassica rapa and Brassica oleracea

Brassica napus (rapeseed) is a recent allotetraploid plant and the second most important oilseed crop worldwide. The origin of B. napus and the genetic relationships with its diploid ancestor species remain largely unresolved. Here, chloroplast DNA (cpDNA) from 488 B. napus accessions of global origin, 139 B. rapa accessions and 49 B. oleracea accessions were populationally resequenced using Illumina Solexa sequencing technologies. The intraspecific cpDNA variants and their allelic frequencies were called genomewide and further validated via EcoTILLING analyses of the rpo region. The cpDNA of the current global B. napus population comprises more than 400 variants (SNPs and short InDels) and maintains one predominant haplotype (Bncp1). Whole‐genome resequencing of the cpDNA of Bncp1 haplotype eliminated its direct inheritance from any accession of the B. rapa or B. oleracea species. The distribution of the polymorphism information content (PIC) values for each variant demonstrated that B. napus has much lower cpDNA diversity than B. rapa; however, a vast majority of the wild and cultivated B. oleracea specimens appeared to share one same distinct cpDNA haplotype, in contrast to its wild C‐genome relatives. This finding suggests that the cpDNA of the three Brassica species is well differentiated. The predominant B. napus cpDNA haplotype may have originated from uninvestigated relatives or from interactions between cpDNA mutations and natural/artificial selection during speciation and evolution. These exhaustive data on variation in cpDNA would provide fundamental data for research on cpDNA and chloroplasts.     

Exploiting BAC-end sequences for the mining,characterization and utility of new short sequences repeat (SSR) markers in Citrus

The aim of this study was to develop a large set of microsatellite markers based on publicly available BAC-end sequences (BESs), and to evaluate their transferability, discriminating capacity of genotypes and mapping ability in Citrus. A set of 1,281 simple sequence repeat (SSR) markers were developed from the 46,339 Citrus clementina BAC-end sequences (BES), of them 20.67% contained SSR longer than 20 bp, corresponding to roughly one perfect SSR per 2.04 kb. The most abundant motifs were di-nucleotide (16.82%) repeats. Among all repeat motifs (TA/AT)n is the most abundant (8.38%), followed by (AG/CT)n (4.51%). Most of the BES-SSR are located in the non-coding region, but 1.3% of BES-SSRs were found to be associated with transposable element (TE). A total of 400 novel SSR primer pairs were synthesized and their transferability and polymorphism tested on a set of 16 Citrus and Citrus relative’s species. Among these 333 (83.25%) were successfully amplified and 260 (65.00%) showed cross-species transferability with Poncirus trifoliata and Fortunella sp. These cross-species transferable markers could be useful for cultivar identification, for genomic study of Citrus, Poncirus and Fortunella sp. Utility of the developed SSR marker was demonstrated by identifying a set of 118 markers each for construction of linkage map of Citrus reticulata and Poncirus trifoliata. Genetic diversity and phylogenetic relationship among 40 Citrus and its related species were conducted with the aid of 25 randomly selected SSR primer pairs and results revealed that citrus genomic SSRs are superior to genic SSR for genetic diversity and germplasm characterization of Citrus spp.     

A chloroplast DNA hypervariable region in eucalypts

Eucalypt chloroplast DNA (cpDNA) provides useful markers for phylogenetic and population research including gene flow and maternity studies. All cpDNA studies in Eucalyptus to date have been based on the RFLP technique, which requires relatively large amounts of clean DNA. The objective of this study was to develop PCR-based cpDNA markers for Eucalyptus. The chloroplast genome of Eucalyptus, like that of most angiosperms, possesses inverted repeats (IR). The two junctions between the IRs and the large single copy (LSC) regions are defined as J_LA andJ_LB. The region surrounding the J_LA junction was sequenced from 26 Eucalyptus DNA samples (21 of E. globulus, plus 5 other species), andtheJ_LB region was sequenced using 5 of these samples. The samples were chosen to represent all major haplotypes identified in previous cpDNA RFLP studies. The J_LA products were 150–210 bp in size, while those fromJ_LB were approximately 500 bp in size. Eighteen mutations were scored in total. Extensive variation was found in the IR within the intergenic spacer between rpl2 and the IR/LSC junctions. Many of these characters were complex indels. One sample of E. globulus possessed a relatively large (38 bp) insertion which caused displacement of the IR/LSC junctions. This is the first report of intraspecific variation in the position of IR/LSC junctions; interspecific variation was also found. The LSC region near J_LB had a low rate of mutation and none were informative. The LSC region near J_LA possessed 2 informative mutations. The variation revealed from these sequences reflects the differentiation previously uncovered in an extensive RLFP analysis on the same samples. These results suggest that the J_LA region will provide very useful cpDNA polymorphisms for future studies in Eucalyptus. Received: 20 March 1999 / Accepted: 16 December 1999     

A multiple-method approach reveals a declining amount of chloroplast DNA during development in Arabidopsis

Background  

A decline in chloroplast DNA (cpDNA) during leaf maturity has been reported previously for eight plant species, including Arabidopsis thaliana. Recent studies, however, concluded that the amount of cpDNA during leaf development in Arabidopsis remained constant.     

Organellar genome analysis of rye (Secale cereale) representing diverse geographic regions.

Rye (Secale cereale) is an important diploid (2n = 14, RR) crop species of the Triticeae and a better understanding of its organellar genome variation can aid in its improvement. Previous genetic analyses of rye focused on the nuclear genome. In the present study, the objective was to investigate the organellar genome diversity and relationships of 96 accessions representing diverse geographic regions using chloroplast (cp) and mitochondrial (mt) DNA PCR-RFLPs. Seven cpDNA and 4 mtDNA coding and noncoding regions were amplified using universal cpDNA and mtDNA primer pairs. Each amplified fragment was digested with 13 different restriction enzymes. mtDNA analysis indicated that the number of polymorphic loci (20) was low and genetic differentiation (GST) was 0.60, excluding the outgroups (hexaploid wheat, Triticum aestivum, 2n = 6x = 42, AABBDD; triticale, xTriticosecale Wittmack, 2n = 6x = 42, AABBRR). cpDNA analysis revealed a low level of polymorphism (40%) among the accessions, and GST was 0.39. Of the 96 genotypes studied, 70 could not be differentiated using cpDNA PCR-RFLPs even though they are from different geographic regions. This is most likely due to germplasm exchange, indicating that genotypes might have a common genetic background. Two cpDNA and 3 mtDNA fragments were significantly correlated to the site of germplasm collection. However, there was no clear trend. These results indicate that the level of organellar polymorphism is low among the cultivated rye genotypes. The cpDNA and mtDNA PCR-RFLP markers used in the present study could be used as molecular markers in rye genetics and breeding programs.     

Comparative analysis of chloroplast DNA in Pyrus species: physical map and gene localization

A physical map of chloroplast DNA (cpDNA) of pear [Pyrus ussuriensis var. hondoensis (Nakai et Kikuchi) Rehder] was constructed using five restriction enzymes, SalI, XhoI, BamHI, SacI and PstI. This information will make it possible to investigate the phylogenetic relationships between Pyrus species. Pear cpDNA was found to be a circular molecule with a total size of about 156 kb in which two inverted repeats of 24.8 kb divide the molecule into small (17 kb) and large (90 kb) single-copy regions. The endonuclease recognition sites in the physical map were determined by single and double digestion of 13 lambda phage clones which covered the entire sequence of the pear cpDNA. Twenty nine genes were localized on the physical map of the pear cpDNA. The structure of pear cpDNA was almost the same in terms of genome size and gene order as that of tobacco cpDNA. RFLP analysis was carried out on cpDNAs from five Pyrus species (Pyrus pyrifolia, Pyrus ussuriensis, Pyrus calleryana, Pyrus elaeagrifolia and Pyrus communis). Two mutations, a recognition-site mutation and a length mutation (deletion), were found only in the cpDNA of P. pyrifolia cultivars. These mutations were localized on the physical map of pear cpDNA. The number of mutations of cpDNA in Pyrus species are small in comparison with those of other angiosperms, suggesting a high degree of genome conservatism in Pyrus species.     

The tortoise and the hare II: relative utility of 21 noncoding chloroplast DNA sequences for phylogenetic analysis

Chloroplast DNA sequences are a primary source of data for plant molecular systematic studies. A few key papers have provided the molecular systematics community with universal primer pairs for noncoding regions that have dominated the field, namely trnL-trnF and trnK/matK. These two regions have provided adequate information to resolve species relationships in some taxa, but often provide little resolution at low taxonomic levels. To obtain better phylogenetic resolution, sequence data from these regions are often coupled with other sequence data. Choosing an appropriate cpDNA region for phylogenetic investigation is difficult because of the scarcity of information about the tempo of evolutionary rates among different noncoding cpDNA regions. The focus of this investigation was to determine whether there is any predictable rate heterogeneity among 21 noncoding cpDNA regions identified as phylogenetically useful at low levels. To test for rate heterogeneity among the different cpDNA regions, we used three species from each of 10 groups representing eight major phylogenetic lineages of phanerogams. The results of this study clearly show that a survey using as few as three representative taxa can be predictive of the amount of phylogenetic information offered by a cpDNA region and that rate heterogeneity exists among noncoding cpDNA regions.     

Sequence-based novel genomic microsatellite markers for robust genotyping purposes in foxtail millet [<Emphasis Type="Italic">Setaria italica</Emphasis> (L.) P. Beauv.]

The unavailability of microsatellite markers and saturated genetic linkage map has restricted the genetic improvement of foxtail millet [Setaria italica (L.) P. Beauv.], despite the fact that in recent times it has been documented as a new model species for biofuel grasses. With the objective to generate a good number of microsatellite markers in foxtail millet cultivar ‘Prasad’, 690 clones were sequenced which generated 112.95 kb high quality sequences obtained from three genomic libraries each enriched with different microsatellite repeat motifs. Microsatellites were identified in 512 (74.2%) of the 690 positive clones and 172 primer pairs (pp) were successfully designed from 249 (48.6%) unique SSR-containing clones. The efficacies of the microsatellite containing genomic sequences were established by superior primer designing ability (69%), PCR amplification efficiency (85.5%) and polymorphic potential (52%) in the parents of F₂ mapping population. Out of 172 pp, functional 147 markers showed high level of cross-species amplification (~74%) in six grass species. Higher polymorphism rate and broad range of genetic diversity (0.30–0.69 averaging 0.58) obtained in constructed phylogenetic tree using 52 microsatellite markers, demonstrated the utility of markers in germplasm characterizations. In silico comparative mapping of 147 foxtail millet microsatellite containing sequences against the mapping data of sorghum (~18%), maize (~16%) and rice (~5%) indicated the presence of orthologous sequences of the foxtail millet in the respective species. The result thus demonstrates the applicability of microsatellite markers in various genotyping applications, determining phylogenetic relationships and comparative mapping in several important grass species.     

Deep genomic analysis of Coelastrella saipanensis (Scenedesmaceae,Chlorophyta): comparative chloroplast genomics of Scenedesmaceae

Many species belonging to the coccoid green algae genus Coelastrella are considered potential candidates for the large-scale production of natural pigments and biofuels. However, little is known about the structural, functional and molecular aspects of the chloroplast genomes (cpDNAs) of this genus. In the present study, the complete sequence of the cpDNA of strain FACHB-2138, which was further identified as Coelastrella saipanensis Hanagata based on morphological and molecular analyses, was elucidated. The 196 140 bp cpDNA sequence that was assembled as a circular map was found to possess the typical quadripartite structure. The two identical copies of 11 897 bp inverted repeat (IR) sequences were separated from one another by single copy regions. The large single copy region (LSC) was 104 949 bp, whereas the small single copy region (SSC) was 67 397 bp. The cpDNA encoded a total of 96 unique genes, which included 67 protein-coding genes, three rRNA genes and 26 tRNA genes. A total of 19 group I introns were annotated in this genome. Comparative analyses with three species from the family Scenedesmaceae showed C. saipanensis had a slightly expanded genome, higher GC content and less skewed distribution of its genes between the two DNA strands than that of the other three species. The cpDNA data deduced from the present study helps to expand our present understanding of plant systematics and phylogenetic reconstruction, and identify the possible biotechnological applications of the species belonging to the studied taxa.     

Comparison among three methods for obtaining chloroplast genome sequences from the conifer Pinus massoniana

The chloroplast genome (CPG) is a powerful tool for phylogenetic studies. Many CPGs have been determined using NGS. However, the large nuclear-genome and difficult CPG-DNA separation in conifers limit their application in related research. In this study, three methods (PCR + Sanger, PCR + HiSeq, cpDNA+HiSeq) for obtaining the CPGs of Pinus massoniana were compared for sequence accuracy, time and cost. PCR + Sanger obtained the most accurate CPGs with advantages in cost (3.08$/kb) and time (2–3 days); PCR + HiSeq generated some DNA fragments with low depth, and the SNPs false-positive-rate (0.44) and sequencing error-rate (0.0265) of this method were higher than those of the cpDNA+HiSeq. Moreover, the cost (~6.17$/kb) and time (4–5 weeks) would significantly increase when HiSeq sequencing were outsourced to sequencing service company. Thus, for the study of intraspecific and interspecies variation in CPGs, CPG sequences can be obtained by comprehensive methods to bridge the method shortcomings. Scuh as sequence accuracy, cost and time.     

Development of Juglans Regia SSR Markers by Data Mining of the EST Database

Walnut (Juglans regia), an economically important woody plant, is widely cultivated in temperate regions for its timber and nutritional fruits. Despite abundant studies in germplasm, systemic molecular evaluations of walnut are sparsely reported mainly due to the limited molecular markers available. Expressed sequence tags (EST) provide a valuable resource for developing simple sequence repeat (SSR) markers. In this study, a total of 5,025 walnut ESTs (covering 16.41 Mb) were retrieved from the National Center for Biotechnology Information database. The SSR motifs were then analyzed by the SSRHunter software. In total, 398 SSRs were obtained with an average frequency of 1/4.08 kb. Dinucleotide (di-) repeat motifs accounted for 69.85% of all SSRs, followed by trinucleotide (tri-) with a frequency of 27.64%, while low frequency (2.51%) of tetranucleotide (tetra-) to hexanucleotide (hexa-) was observed. Meanwhile, GCA and TC motifs were prevalent among di- and tri- loci, respectively. Subsequently, a total of 123 primer pairs were designed from the non-redundant SSR-containing unigenes with the selection threshold of SSR length set to 10 bp or more. To examine the efficiency of candidate markers, seven DNA pools were collected from geographically different accessions. Results demonstrated that 41 SSR primer sets could generate high polymorphic amplification products (33.3%), and these polymorphic loci were mainly located in the 3′-untranslated region. Annotation analysis revealed that only two of these 41 loci were located inside open reading frames of characterized proteins (E ≤ 1E−30).     

Transferability of microsatellite markers from Eucalyptus spp. to Eugenia dysenterica (Myrtaceae family)

Cagaita (Eugenia dysenterica) is a widespread plant found in the Brazilian cerrado. Its fruit is used for popular consumption and for industrial purposes. A battery of 346 primer pairs developed for Eucalyptus spp. was tested on cagaita. Only 10 primer pairs were found to be transferable between the two species. Using a polyacrilamide gel, an average of 10.4 alleles per locus was detected, in a sample of 116 individuals from 10 natural populations of cagaita. Seven polymorphic loci allowed estimation of genetic parameters, including expected average heterozigosity H_E = 0.442, diversity among populations, R_ST = 0.268 and gene flow Nm = 0.680.     

A comparative analysis of complete plastid genomes from Prangos fedtschenkoi and Prangos lipskyi (Apiaceae)

Prangos fedtschenkoi (Regel & Schmalh.) Korovin and P. lipskyi Korovin (Apiaceae) are rare plant species endemic to mountainous regions of Middle Asia. Both are edificators of biotic communities and valuable resource plants. The results of recent phylogenetic analyses place them in Prangos subgen. Koelzella (M. Hiroe) Lyskov & Pimenov and suggest they may possibly represent sister species. To aid in development of molecular markers useful for intraspecific phylogeographic and population‐level genetic studies of these ecologically and economically important plants, we determined their complete plastid genome sequences and compared the results obtained to several previously published plastomes of Apiaceae. The plastomes of P. fedtschenkoi and P. lipskyi are typical of Apiaceae and most other higher plant plastid DNAs in their sizes (153,626 and 154,143 bp, respectively), structural organization, gene arrangement, and gene content (with 113 unique genes). A total of 49 and 48 short sequence repeat (SSR) loci of 10 bp or longer were detected in P. fedtschenkoi and P. lipskyi plastomes, respectively, representing 42–43 mononucleotides and 6 AT dinucleotides. Seven tandem repeats of 30 bp or longer with a sequence identity ≥90% were identified in each plastome. Further comparisons revealed 319 polymorphic sites between the plastomes (IR, 21; LSC, 234; SSC, 64), representing 43.8% transitions (Ts), 56.1% transversions (Tv), and a Ts/Tv ratio of 0.78. Within genic regions, two indel events were observed in rpoA (6 and 51 bp) and ycf1 (3 and 12 bp), and one in ndhF (6 bp). The most variable intergenic spacer region was that of accD/psaI, with 21.1% nucleotide divergence. Each Prangos species possessed one of two separate inversions (either 5 bp in ndhB intron or 9 bp in petB intron), and these were predicted to form hairpin structures with flanking repeat sequences of 18 and 19 bp, respectively. Both species have also incorporated novel DNA in the LSC region adjacent to the LSC/IRa junction, and BLAST searches revealed it had a 100 bp match (86% sequence identity) to noncoding mitochondrial DNA. Prangos‐specific primers were developed for the variable accD/psaI intergenic spacer and preliminary PCR‐surveys suggest that this region will be useful for future phylogeographic and population‐level studies.     

Chloroplast DNA variation and phylogeny of theRanunculaceae

Restriction site mapping of chloroplast DNA from 31 species representing 26 genera of theRanunculaceae was performed using eleven restriction endonucleases. The chloroplast genome varies in length from approximately 152 to 160 kb. Length variants are frequent in theRanunculaceae and range from usually less than 300 bp to rarely 1.5 kb. The inverted repeat is extended into the large single copy (LSC) region by 4–4.5 kb inAnemone, Clematis, Clematopsis, Hepatica, Knowltonia, andPulsatilla. Several inversions are present in the LSC-region of the cpDNA in all these genera and inAdonis. The frequency of restriction site mutations varies within the chloroplast genome in theRanunculaceae between 4 and 32 mutations per kilobase, and is lowest in the inverted repeat and the regions containing the ATPase-genes and the genespsaA, psaB, psbA, rpoB, andrbcL. A total of 547 phylogenetically informative restriction sites was utilized in cladistic analyses of the family using Wagner, Dollo, and weighted parsimony. These three parsimony analyses result in different tree topologies. Four, six, and one equally most parsimonious trees were obtained with Wagner, Dollo, and weighted parsimony, respectively. The amount of support for the monophyletic groups was evaluated using bootstrapping and decay analysis. All three parsimony methods suggest thatHydrastis is the sister group to the remainder of theRanunculaceae, and that theAnemone-Clematis group, which shares several derived cpDNA rearrangements, is monophyletic. Only a few of the traditional groups in theRanunculaceae are supported by cpDNA restriction side data. Only Dollo parsimony provides support for the hypothesis thatThalictroideae andRanunculoideae are monophyletic.     

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.     

Molecular identification of the economically important freshwater mussels (Mollusca–Bivalvia–Unionoida) of Thailand: developing species‐specific markers from AFLPs

Shells of certain freshwater mussel (Unionoida) species are highly demanded and serve as raw material for a range of decorative and pharmaceutical products. In Thailand, most animals for this purpose are currently harvested from wild populations, with unionoid culture still being in its infancy. Whilst reliable species identification is a prerequisite for developing a large‐scale industry, identification by morphological means is hampered by extensive phenotypic plasticity and poor knowledge of species delimitations. To facilitate alternative molecular identification, we developed species‐specific markers for the three Thai unionoids with considerable economic potential (CEP): that is, Chamberlainia hainesiana, Hyriopsis desowitzi and Hyriopsis myersiana. For this purpose, amplified fragment length polymorphism (AFLP) fingerprints using 24 specific primer pairs were generated for eight samples of each CEP species and four samples of the closely related, non‐CEP species Contradens contradens. Cloning and sequencing of 13 CEP species‐specific AFLP bands revealed fragment collision at three occasions. In total, 16 species‐specific primer pairs were designed and tested on 92 Thai specimens spanning seven species and four genera. Thereby, specificity of (1) three primers to C. hainesiana, (2) one primer to H. desowitzi + Hyriopsis bialata, (3) one primer to H. myersiana + H. bialata and (4) four primers to all three Hyriopsis species tested was confirmed. Respective multiplex PCR protocols are provided. The developed primers enable cheap, quick and reliable identification of the Thai CEP species by one to three PCRs and offer a tool for a range of additional applications within mussel culture and ecological and evolutionary research on these important organisms.