Universal primers for the amplification of chloroplast microsatellites in grasses (Poaceae)

Molecular discrimination of tall fescue morphotypes in association with Festuca relatives

Tall fescue (Festuca arundinacea Schreb.) is an important cool-season perennial grass species used as forage and turf, and in conservation plantings. There are three morphotypes in hexaploid tall fescue: Continental, Mediterranean and Rhizomatous. This study was conducted to develop morphotype-specific molecular markers to distinguish Continental and Mediterranean tall fescues, and establish their relationships with other species of the Festuca genus for genomic inference. Chloroplast sequence variation and simple sequence repeat (SSR) polymorphism were explored in 12 genotypes of three tall fescue morphotypes and four Festuca species. Hypervariable chloroplast regions were retrieved by using 33 specifically designed primers followed by sequencing the PCR products. SSR polymorphism was studied using 144 tall fescue SSR primers. Four chloroplast (NFTCHL17, NFTCHL43, NFTCHL45 and NFTCHL48) and three SSR (nffa090, nffa204 and nffa338) markers were identified which can distinctly differentiate Continental and Mediterranean morphotypes. A primer pair, NFTCHL45, amplified a 47 bp deletion between the two morphotypes is being routinely used in the Noble Research Institute's core facility for morphotype discrimination. Both chloroplast sequence variation and SSR diversity showed a close association between Rhizomatous and Continental morphotypes, while the Mediterranean morphotype was in a distant clade. F. pratensis and F. arundinacea var. glaucescens, the P and G1G2 genome donors, respectively, were grouped with the Continental clade, and F. mairei (M1M2 genome) grouped with the Mediterranean clade in chloroplast sequence variation, while both F. pratensis and F. mairei formed independent clade in SSR analysis. Age estimation based on chloroplast sequence variation indicated that the Continental and Mediterranean clades might have been colonized independently during 0.65 ± 0.06 and 0.96 ± 0.1 million years ago (Mya) respectively. The findings of the study will enhance tall fescue breeding for persistence and productivity.

Integral Phylogenomic Approach over Ilex L. Species from Southern South America

The use of molecular markers with inadequate variation levels has resulted in poorly resolved phylogenetic relationships within Ilex. Focusing on southern South American and Asian species, we aimed at contributing informative plastid markers. Also, we intended to gain insights into the nature of morphological and physiological characters used to identify species. We obtained the chloroplast genomes of I.paraguariensis and I. dumosa, and combined these with all the congeneric plastomes currently available to accomplish interspecific comparisons and multilocus analyses. We selected seven introns and nine IGSs as variable non-coding markers that were used in phylogenomic analyses. Eight extra IGSs were proposed as candidate markers. Southern South American species formed one lineage, except for I. paraguariensis, I. dumosa and I. argentina, which occupied intermediate positions among sampled taxa; Euroasiatic species formed two lineages. Some concordant relationships were retrieved from nuclear sequence data. We also conducted integral analyses, involving a supernetwork of molecular data, and a simultaneous analysis of quantitative and qualitative morphological and phytochemical characters, together with molecular data. The total evidence tree was used to study the evolution of non-molecular data, evidencing fifteen non-ambiguous synapomorphic character states and consolidating the relationships among southern South American species. More South American representatives should be incorporated to elucidate their origin.

Comparative assessment of genetic diversity in cytoplasmic and nuclear genome of upland cotton

The importance of the cytoplasmic genome for many economically important traits is well documented in several crop species, including cotton. There is no report on application of cotton chloroplast specific SSR markers as a diagnostic tool to study genetic diversity among improved Upland cotton lines. The complete plastome sequence information in GenBank provided us an opportunity to report on 17 chloroplast specific SSR markers using a cost-effective data mining strategy. Here we report the comparative analysis of genetic diversity among a set of 42 improved Upland cotton lines using SSR markers specific to chloroplast and nuclear genome, respectively. Our results revealed that low to moderate level of genetic diversity existed in both nuclear and cytoplasm genome among this set of cotton lines. However, the specific estimation suggested that genetic diversity is lower in cytoplasmic genome compared to the nuclear genome among this set of Upland cotton lines. In summary, this research is important from several perspectives. We detected a set of cytoplasm genome specific SSR primer pairs by using a cost-effective data mining strategy. We reported for the first time the genetic diversity in the cytoplasmic genome within a set of improved Upland cotton accessions. Results revealed that the genetic diversity in cytoplasmic genome is narrow, compared to the nuclear genome within this set of Upland cotton accessions. Our results suggested that most of these polymorphic chloroplast SSRs would be a valuable complementary tool in addition to the nuclear SSR in the study of evolution, gene flow and genetic diversity in Upland cotton.

Development of chloroplast genomic resources for Cynara

In this study, new chloroplast (cp) resources were developed for the genus Cynara, using whole cp genomes from 20 genotypes, by means of high-throughput sequencing technologies. Our target species included seven globe artichokes, two cultivated cardoons, eight wild artichokes, and three other wild Cynara species (C. baetica, C. cornigera and C. syriaca). One complete cp genome was isolated using short reads from a whole-genome sequencing project, while the others were obtained by means of long-range PCR, for which primer pairs are provided here. A de novo assembly strategy combined with a reference-based assembly allowed us to reconstruct each cp genome. Comparative analyses among the newly sequenced genotypes and two additional Cynara cp genomes ('Brindisino' artichoke and C. humilis) retrieved from public databases revealed 126 parsimony informative characters and 258 singletons in Cynara, for a total of 384 variable characters. Thirty-nine SSR loci and 34 other INDEL events were detected. After data analysis, 37 primer pairs for SSR amplification were designed, and these molecular markers were subsequently validated in our Cynara genotypes. Phylogenetic analysis based on all cp variable characters provided the best resolution when compared to what was observed using only parsimony informative characters, or only short 'variable' cp regions. The evaluation of the molecular resources obtained from this study led us to support the 'super-barcode' theory and consider the total cp sequence of Cynara as a reliable and valuable molecular marker for exploring species diversity and examining variation below the species level.

Analysis of variation in chloroplast DNA sequences

This chapter introduces and reviews methods for analyzing variation in chloroplast DNA, mainly by polymerase chain reaction (PCR) and subsequent revelation of polymorphisms. Sources for chloroplast primers are discussed, as well as methods such as Sanger sequencing, PCR followed by restriction fragment length polymorphism (RFLP), gel electrophoresis, fragment analysis on automated DNA sequencers, denaturing high-performance liquid chromatography (dHPLC), and next-generation sequencing (NGS). A special section deals with peculiarities of chloroplast DNA variation, such as tandem repeats and mini- and microsatellites.

Development of chloroplast simple sequence repeats (cpSSRs) for the intraspecific study of Gracilaria tenuistipitata (Gracilariales, Rhodophyta) from different populations

Background

Gracilaria tenuistipitata is an agarophyte with substantial economic potential because of its high growth rate and tolerance to a wide range of environment factors. This red seaweed is intensively cultured in China for the production of agar and fodder for abalone. Microsatellite markers were developed from the chloroplast genome of G. tenuistipitata var. liui to differentiate G. tenuistipitata obtained from six different localities: four from Peninsular Malaysia, one from Thailand and one from Vietnam. Eighty G. tenuistipitata specimens were analyzed using eight simple sequence repeat (SSR) primer-pairs that we developed for polymerase chain reaction (PCR) amplification.

Findings

Five mononucleotide primer-pairs and one trinucleotide primer-pair exhibited monomorphic alleles, whereas the other two primer-pairs separated the G. tenuistipitata specimens into two main clades. G. tenuistipitata from Thailand and Vietnam were grouped into one clade, and the populations from Batu Laut, Middle Banks and Kuah (Malaysia) were grouped into another clade. The combined dataset of these two primer-pairs separated G. tenuistipitata obtained from Kelantan, Malaysia from that obtained from other localities.

Conclusions

Based on the variations in repeated nucleotides of microsatellite markers, our results suggested that the populations of G. tenuistipitata were distributed into two main geographical regions: (i) populations in the west coast of Peninsular Malaysia and (ii) populations facing the South China Sea. The correct identification of G. tenuistipitata strains with traits of high economic potential will be advantageous for the mass cultivation of seaweeds.

Development of chloroplast microsatellite markers for phylogenetic analysis in Brassicaceae

By employing in silico tools, we devised new chloroplast microsatellite primers for inferring phylogenetic relationships within Brassicaceae. Microsatellite repeats were scanned in 12 chloroplast genomes of Brassicaceae, regions flanking these repeats were aligned and 19 universal primers were designed. Fifteen of these primer pairs are predicted to yield polymorphic amplicons, that are more or less evenly distributed throughout the chloroplast genomes. Finally, using PCR, we have validated three primer pairs on a limited 'test set' of plants, different from those used in computational analysis.

Contrasting patterns of pollen and seed flow influence the spatial genetic structure of sweet vernal grass (Anthoxanthum odoratum) populations

The spatial genetic structure of plant populations is determined by a combination of gene flow, genetic drift, and natural selection. Gene flow in most plants can result from either seed or pollen dispersal, but detailed investigations of pollen and seed flow among populations that have diverged following local adaptation are lacking. In this study, we compared pollen and seed flow among 10 populations of sweet vernal grass (Anthoxanthum odoratum) on the Park Grass Experiment. Overall, estimates of genetic differentiation that were based on chloroplast DNA (cpDNA) and, which therefore resulted primarily from seed flow, were lower (average F(ST) = 0.058) than previously published estimates that were based on nuclear DNA (average F(ST) = 0.095). Unlike nuclear DNA, cpDNA showed no pattern of isolation by adaptation; cpDNA differentiation was, however, inversely correlated with the number of additions (nutrients and lime) that each plot had received. We suggest that natural selection is restricting pollen flow among plots, whereas nutrient additions are increasing seed flow and genetic diversity by facilitating the successful germination and growth of immigrant seeds. This study highlights the importance of considering all potential gene flow mechanisms when investigating determinants of spatial genetic structure, and cautions against the widespread assumption that pollen flow is more important than seed flow for population connectivity in wind-pollinated species.

Evidence for the equal resilience of Triodia spp. (Poaceae), from different functional groups, to frequent fire dating back to the late Pleistocene

Species with different regenerative responses to fire are hypothesised to coexist by utilising the different temporal and spatial niche opportunities created by the stochasticity of the fire regime. This is strongly supported by observations of instability of species' presence and abundance at the local scale while these are stable at the community scale. However, observations of species coexistence in fire-prone communities are limited to several decades only. To improve the robustness of this hypothesis, coalescent analysis, using chloroplast microsatellites, was undertaken on three sympatric species of Triodia from different functional groups in the fire-prone Kimberley region of Western Australia. The results inferred that T. bitextura, an obligate resprouter, Triodia sp., an obligate seeder, and T. epactia, a facultative resprouter, had mean T(mrca) values of 65k, 40k and 111k generations, respectively. Using a mutation rate of 3.2 × 10(-5) and a generation time of 5 years gave T(mrca) values of 436k, 203k and 556 k years, respectively. These results provide evidence for the coexistence of these species to the same fire regime dating back to the late Pleistocene. It also demonstrates the long-term resilience of an obligate seeder, Triodia sp., in a frequently burnt environment at the community scale.

A set of 100 chloroplast DNA primer pairs to study population genetics and phylogeny in monocotyledons

Chloroplast DNA sequences are of great interest for population genetics and phylogenetic studies. However, only a small set of markers are commonly used. Most of them have been designed for amplification in a large range of Angiosperms and are located in the Large Single Copy (LSC). Here we developed a new set of 100 primer pairs optimized for amplification in Monocotyledons. Primer pairs amplify coding (exon) and non-coding regions (intron and intergenic spacer). They span the different chloroplast regions: 72 are located in the LSC, 13 in the Small Single Copy (SSC) and 15 in the Inverted Repeat region (IR). Amplification and sequencing were tested in 13 species of Monocotyledons: Dioscorea abyssinica, D. praehensilis, D. rotundata, D. dumetorum, D. bulbifera, Trichopus sempervirens (Dioscoreaceae), Phoenix canariensis, P. dactylifera, Astrocaryum scopatum, A. murumuru, Ceroxylon echinulatum (Arecaceae), Digitaria excilis and Pennisetum glaucum (Poaceae). The diversity found in Dioscorea, Digitaria and Pennisetum mainly corresponded to Single Nucleotide Polymorphism (SNP) while the diversity found in Arecaceae also comprises Variable Number Tandem Repeat (VNTR). We observed that the most variable loci (rps15-ycf1, rpl32-ccsA, ndhF-rpl32, ndhG-ndhI and ccsA) are located in the SSC. Through the analysis of the genetic structure of a wild-cultivated species complex in Dioscorea, we demonstrated that this new set of primers is of great interest for population genetics and we anticipate that it will also be useful for phylogeny and bar-coding studies.

Evolution and polyploid origins in North American Arctic Puccinellia (Poaceae) based on nuclear ribosomal spacer and chloroplast DNA sequences

The proportion of polyploid plant species increases at higher latitudes, and it has been suggested that original postglacial Arctic immigrants of some large groups, including grasses, were polyploid. We analyzed noncoding nuclear and chloroplast DNA of all North American diploid Puccinellia (Poaceae) and a subset of arctic polyploids to hypothesize evolutionary relationships among diploids and to evaluate the parentage of polyploids. Diploids formed three lineages: one uniting arctic species P. arctica and P. banksiensis; a second comprising arctic species P. tenella, P. alaskana, P. vahliana, and P. wrightii; and a third uniting the two temperate species P. lemmonii and P. parishii. The arctic species P. angustata (hexaploid) and P. andersonii (primarily octoploid) apparently derive from the P. arctica-P. banksiensis lineage based on ITS and chloroplast sequences, and share an ancestor with arctic triploid/tetraploid P. phryganodes based on nrDNA sequences. Sequence comparisons also suggest tetraploid P. bruggemannii evolved from two arctic lineages: P. vahliana-P. wrightii and P. arctica-P. banksiensis. These patterns and the predominance of arctic rather than temperate diploid species support the idea that diploid Puccinellia recolonized the Arctic from northern glacial refugia like Beringia, and also formed stabilized polyploid hybrids during these refugial events or subsequently during postglacial colonization.

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.

Advances in molecular marker techniques and their applications in plant sciences

Detection and analysis of genetic variation can help us to understand the molecular basis of various biological phenomena in plants. Since the entire plant kingdom cannot be covered under sequencing projects, molecular markers and their correlation to phenotypes provide us with requisite landmarks for elucidation of genetic variation. Genetic or DNA based marker techniques such as RFLP (restriction fragment length polymorphism), RAPD (random amplified polymorphic DNA), SSR (simple sequence repeats) and AFLP (amplified fragment length polymorphism) are routinely being used in ecological, evolutionary, taxonomical, phylogenic and genetic studies of plant sciences. These techniques are well established and their advantages as well as limitations have been realized. In recent years, a new class of advanced techniques has emerged, primarily derived from combination of earlier basic techniques. Advanced marker techniques tend to amalgamate advantageous features of several basic techniques. The newer methods also incorporate modifications in the methodology of basic techniques to increase the sensitivity and resolution to detect genetic discontinuity and distinctiveness. The advanced marker techniques also utilize newer class of DNA elements such as retrotransposons, mitochondrial and chloroplast based microsatellites, thereby revealing genetic variation through increased genome coverage. Techniques such as RAPD and AFLP are also being applied to cDNA-based templates to study patterns of gene expression and uncover the genetic basis of biological responses. The review details account of techniques used in identification of markers and their applicability in plant sciences.

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

Although the chloroplast genome contains many noncoding regions, relatively few have been exploited for interspecific phylogenetic and intraspecific phylogeographic studies. In our recent evaluation of the phylogenetic utility of 21 noncoding chloroplast regions, we found the most widely used noncoding regions are among the least variable, but the more variable regions have rarely been employed. That study led us to conclude that there may be unexplored regions of the chloroplast genome that have even higher relative levels of variability. To explore the potential variability of previously unexplored regions, we compared three pairs of single-copy chloroplast genome sequences in three disparate angiosperm lineages: Atropa vs. Nicotiana (asterids); Lotus vs. Medicago (rosids); and Saccharum vs. Oryza (monocots). These three separate sequence alignments highlighted 13 mutational hotspots that may be more variable than the best regions of our former study. These 13 regions were then selected for a more detailed analysis. Here we show that nine of these newly explored regions (rpl32-trnL((UAG)), trnQ((UUG))-5'rps16, 3'trnV((UAC))-ndhC, ndhF-rpl32, psbD-trnT((GGU)), psbJ-petA, 3'rps16-5'trnK((UUU)), atpI-atpH, and petL-psbE) offer levels of variation better than the best regions identified in our earlier study and are therefore likely to be the best choices for molecular studies at low taxonomic levels.

Length polymorphism scanning is an efficient approach for revealing chloroplast DNA variation

Phylogeographic and population genetic screens of chloroplast DNA (cpDNA) provide insights into seed-based gene flow in angiosperms, yet studies are frequently hampered by the low mutation rate of this genome. Detection methods for intraspecific variation can be either direct (DNA sequencing) or indirect (PCR–RFLP), although no single method incorporates the best features of both approaches. We show that screening universal chloroplast amp li cons for length polymorphism provides an accurate and efficient method for identifying cpDNA variation. By sequencing 4500 bp of cpDNA from 17 accessions of Purshia tridentata (bitterbrush), we detected 9 haplotypes, 8 of which were identifiable by unique multilocus length combinations resolvable by automated fragment analysis. In silico estimates of PCR–RFLP for these loci show that 5 haplotypes would be resolved by agarose electrophoresis. A survey of 4 intraspecific data sets from diverse angiosperms revealed that length variation in cpDNA amplicons is nearly ubiquitous, and 61 of 67 haplotypes identified by direct sequencing could be identified by screening length variation. Combined with automated fluorescent detection, length polymorphism screening of universal cpDNA regions offers a simple screen for intraspecific variation that can be used across angiosperms with minimal optimization, providing detection limits that rival direct sequencing at a fraction of the cost.Key words: cpDNA, intraspecific polymorphism, population genetics, phylogeography, indels.

Reinforcement of reproductive isolation between adjacent populations in the Park Grass Experiment

It has been debated, ever since Charles Darwin and Alfred Russell Wallace disagreed about the matter, whether natural selection plays a role in reinforcing reproductive isolation during the earliest stages of speciation. Recent theory suggests that it can do so, but until now the empirical evidence has conspicuously lacked a case in which reinforcement has actually been observed to split a population. We show that this has occurred at least once in populations of the grass Anthoxanthum odoratum growing in the Park Grass Experiment where flowering time has shifted at the boundaries between plots. As a consequence, gene flow via pollen has been severely limited and adjacent populations that had a common origin at the start of the experiment in 1856 have now diverged at neutral marker loci.

Chloroplast microsatellites: measures of genetic diversity and the effect of homoplasy

Chloroplast microsatellites have been widely used in population genetic studies of conifers in recent years. However, their haplotype configurations suggest that they could have high levels of homoplasy, thus limiting the power of these molecular markers. A coalescent-based computer simulation was used to explore the influence of homoplasy on measures of genetic diversity based on chloroplast microsatellites. The conditions of the simulation were defined to fit isolated populations originating from the colonization of one single haplotype into an area left available after a glacial retreat. Simulated data were compared with empirical data available from the literature for a species of Pinus that has expanded north after the Last Glacial Maximum. In the evaluation of genetic diversity, homoplasy was found to have little influence on Nei's unbiased haplotype diversity (H(E)) while Goldstein's genetic distance estimates (D2sh) were much more affected. The effect of the number of chloroplast microsatellite loci for evaluation of genetic diversity is also discussed.