Next Article in Journal
Experimental Evaluation of a New Giant Reed (Arundo Donax L.) Composite Using Citric Acid as a Natural Binder
Previous Article in Journal
Conventional Industrial Robotics Applied to the Process of Tomato Grafting Using the Splicing Technique
Previous Article in Special Issue
Factors Underlying Seed Yield in Red Clover: Review of Current Knowledge and Perspectives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

DNA-Based Assessment of Genetic Diversity in Grassland Plant Species: Challenges, Approaches, and Applications

Molecular Plant Breeding, Institute of Agricultural Sciences, ETH Zurich, Universitaetstrasse 2, 8092 Zurich, Switzerland
*
Author to whom correspondence should be addressed.
Submission received: 14 November 2019 / Revised: 6 December 2019 / Accepted: 10 December 2019 / Published: 12 December 2019

Abstract

:
Grasslands are wide-spread, multi-species ecosystems that provide many valuable services. Plant genetic diversity (i.e., the diversity within species) is closely linked to ecosystem functioning in grasslands and constitutes an important reservoir of genetic resources that can be used to breed improved cultivars of forage grass and legume species. Assessing genetic diversity in grassland plant species is demanding due to the large number of different species and the level of resolution needed. However, recent methodological advances could help in tackling this challenge at a larger scale. In this review, we outline the methods that can be used to measure genetic diversity in plants, highlighting their strengths and limitations for genetic diversity assessments of grassland plant species, with a special focus on forage plants. Such methods can be categorized into DNA fragment, hybridization array, and high-throughput sequencing (HTS) methods, and they differ in terms of resolution, throughput, and multiplexing potential. Special attention is given to HTS approaches (i.e., plastid genome skimming, whole genome re-sequencing, reduced representation libraries, sequence capture, and amplicon sequencing), because they enable unprecedented large-scale assessments of genetic diversity in non-model organisms with complex genomes, such as forage grasses and legumes. As no single method may be suited for all kinds of purposes, we also provide practical perspectives for genetic diversity analyses in forage breeding and genetic resource conservation efforts.

1. Introduction

Grasslands are among the most widespread land ecosystems of the planet, covering more than 50 million square kilometres, roughly 40% of the land surface [1]. In addition to their fundamental role as sources of roughage for ruminant livestock, grasslands provide important ecosystems services related to water (e.g., water provisioning and water flow regulation), climate (e.g., carbon sequestration and storage), soil quality (e.g., erosion prevention), biodiversity (e.g., habitats for pollinating insects) and cultural heritage [2]. They can be categorized according to their management intensity into natural, semi-natural, and improved grasslands. Natural grasslands are dominated by species-rich, herbaceous vegetation and developed with no or minimal human intervention (i.e., prevention of woody overgrowth caused by natural succession). Semi-natural grasslands depend on minimal human management, such as moderate grazing and occasional mowing [3,4]. Improved grasslands, on the other hand, are usually intensively managed and tend to harbour only high yielding forage crop cultivars [2]. Over the last century, the surface area covered by natural and semi-natural grassland has shrunk to make place for arable land or, in the case of some semi-natural grasslands, due to the abandonment of cultivation [2,5]. Demographic growth and changing dietary preferences are driving an increase in the demand for livestock products, contributing to the rise of intensively managed grassland worldwide [2]. The combination of the decline of natural and semi-natural grasslands and the rise of intensified grasslands are thus threatening plant biodiversity in temperate regions of the world, where natural and semi-natural grasslands are important biodiversity hotspots [5]. Such a loss of biodiversity can have great implications for the stability of grassland ecosystems and for food production chains, which calls for targeted conservation efforts [6].
Plant biodiversity is commonly measured as the diversity between species (i.e., species richness) and the diversity within species (i.e., genetic diversity). Species richness can be routinely monitored in natural and semi-natural grasslands with well-established methods such as botanical surveys. Such monitoring methods enable the study of plant species richness in semi-natural grasslands and its changes over time [6,7,8]. In contrast, monitoring genetic diversity is far less straightforward: technical shortcomings limit measuring genetic diversity at large scales and for all relevant grassland species [9].
The genetic diversity of a species consists of the genetic differences among individual plants and encompasses the allelic diversity of genes. The ability of grassland plant species to colonize diverse geographical and climatic zones has given rise to a large pool of genetic resources, into which plant breeders can tap to improve forage crop cultivars [10,11]. Due to local adaptation, plants from such natural or ecotype populations have been shown to produce higher biomass yields when compared to otherwise highly productive commercial cultivars [12,13]. Thus, such genetic resources are particularly important for the development of forage crop cultivars adapted to specific environmental conditions. In addition to its value for breeding, the genetic diversity of plants can also affect ecosystem functioning in grasslands. High genetic diversity was shown to provide an insurance effect against environmental stresses, such as drought or frost [14,15]. Furthermore, grassland plant genetic diversity may also stabilize and, in some cases, increase population biomass yields [16].
Utilizing the above-mentioned benefits of plant genetic diversity largely depends on the ability to assess it in a standardized way. Assessing genetic diversity could assist plant breeders in ensuring that the diversity of the breeding germplasm is sufficient for further genetic gain and in making informed selection decisions [17]. It could also aid in managing and protecting the valuable genetic resources that are present in natural populations or in landraces (i.e., in situ conservation), as well as in gene banks (i.e., ex situ conservation) [17]. Genetic diversity assessments could further help in guiding management practices that influence plant genetic diversity (e.g., grazing and mowing) [18]. Standardized genetic diversity assessments could also help forage breeders to comply with the distinctness, uniformity, and stability (DUS) requirements of cultivar registration, which are difficult to achieve with the highly diverse populations of outcrossing grassland plant species [19].
Genetic diversity assessments are nowadays dominated by DNA-based methods. They can be categorized into fragment analysis-based methods (i.e., methods in which DNA polymorphisms are detected by comparing the sizes of DNA fragments), hybridization array-based methods (i.e., methods in which DNA polymorphisms are detected by hybridizing the sample DNA to arrayed probes), and sequencing-based methods (i.e., methods in which DNA polymorphisms are detected by sequencing). In contrast to alternative methodologies such as morphological and isozyme-based assessments, with DNA-based methods, potentially thousands of DNA polymorphisms can be assessed simultaneously in multiple plants, which greatly reduces costs and labour. Some DNA-based methods can target both coding and non-coding regions of the genome, which provides a wider view of genetic diversity and could potentially be used to detect genomic loci under selection. Furthermore, recent developments in DNA sequencing technologies enable unprecedented large-scale assessments of genetic diversity in non-model organisms with complex genomes, such as grassland plant species.
Optimally exploiting the technical capabilities of such methods to study the genetic diversity of grassland plant populations will require a careful balance between the scope of the assessments, the biological features of the plant species under study, and the features of the method. In this review, we describe DNA-based methods that can be used to assess plant genetic diversity, highlighting their strengths and limitations, as well as their potential for application in grassland plant species. Whenever available, we describe exemplary studies on grasses and legumes, the two major grassland plant families used as forages. We finish by providing practical perspectives on using large-scale genetic diversity monitoring to inform forage breeding, improve genetic resource conservation efforts, and to deepen the study of grassland ecology.

2. DNA Fragment Size Analysis

Genetic diversity assessments based on DNA fragment analysis rely on the size variation among DNA fragments generated from different individuals or populations. Such fragments can be produced by restriction enzyme digestion (e.g., restriction fragment length polymorphism, RFLP [20]), by random PCR amplification (e.g., random amplified polymorphic DNA, RAPD [21]), or by a combination of restriction enzyme digestion, adapter ligation and PCR amplification (e.g., amplified fragment length polymorphism, AFLP [22]). Another approach to generate polymorphic DNA fragments is using primers that target general sequence patterns of coding regions (e.g., sequence-related amplified polymorphism, SRAP [23]). Polymorphic DNA fragments can also be generated by amplifying repetitive DNA, such as retrotransposons (e.g., inter-primer binding site amplification, iPBS [24]) or microsatellites (e.g., simple sequence repeats, SSR [25]), or by using the SSR sequences themselves as primers to amplify the sequence between two SSR (e.g., inter simple sequence repeats, ISSR [26]).
DNA fragment analysis methods can be divided into dominant and co-dominant marker systems, depending on whether only one (dominant) or all (co-dominant) alleles of a locus can be detected. An important advantage of dominant marker systems, such as AFLP, RAPD, iPBS, SRAP, and ISSR, is that markers can be produced at low cost and with no sequence information of the target species needed. This makes them the systems of choice for many grassland plant species where genomic sequence information is scarce. Although a major drawback of dominant DNA marker systems is that they are not able to detect heterozygotes, they have been successfully used to study the diversity of grassland plant species, such as orchardgrass (Dactylis glomerata L.), perennial ryegrass (Lolium perenne L.), Italian ryegrass (L. multiflorum Lam.), tall oat-grass (Arrhenatherum elatius L.), tall fescue (Festuca arundinacea Schreb., syn. Lolium arundinaceum and Schenodurus arundinaceus), meadow fescue (F pratensis Huds., syn. Lolium pratense and Schenodurus pratensis), Kentucky bluegrass (Poa pratensis L.), bird’s-foot trefoil (Lotus corniculatus L.), red clover (Trifolium pratense L.) and white clover (T. repens L.; Table 1).
Co-dominant marker systems, such as SSR, can detect heterozygotes and enable diversity analyses based on allele frequencies. Developing SSR requires significant input and a priori sequence information. However, once obtained, SSR primer sets can be easily standardized and used by multiple laboratories. SSR are highly polymorphic, and usually, 10 to 30 SSR primer combinations are used to assess the diversity of a grassland plant population (Table 1). Genetic diversity estimations obtained with SSR are generally higher than those of dominant marker systems [42]. SSR-based genetic diversity analyses have been reported for orchardgrass, perennial and Italian ryegrass, tall and meadow fescue, timothy (Phleum pratense L.), and white clover (Table 1).
Most economically relevant grassland plants are outbreeding, i.e., they are unable to self-pollinate. Consequently, grassland plant populations tend to be highly genetically diverse, and multiple plants are needed to obtain representative samples of the original population. Sample size is dictated by the frequency of the alleles to be detected. For outbreeding species, the frequency of a given allele in a sample follows a binomial distribution [68]. Thus, to sample alleles that are present in the original population with a frequency of >10%, approximately 30 plants would be necessary at a probability of p = 0.95 [68]. Mainly due to practical limitations, DNA fragment-based studies usually analysed 10 to 30 plants per population, which gave them enough power to detect alleles with a frequency of 25% to 10%, respectively. In tall fescue [69] and Kentucky bluegrass, as few as 16 or three plants per population, respectively, have been found sufficient to detect genetic differentiation among populations [54,58]. Such results suggest that frequent alleles are driving genetic differentiation among genetically distant populations. More plants may be necessary to assess genetic diversity and differentiation in closely related grassland plant populations and to detect rare alleles [68,70]. A useful guideline could be using samples of at least 60 plants, the recommended sample size to characterize new ryegrass cultivars [71]. The number of analysed fragments may also influence genetic diversity estimations. In principle, the more loci are analysed, the more likely it is that polymorphisms will be detected. However, as the number of polymorphisms reaches saturation, each additional fragment stops adding information. For example, 120 AFLP fragments were found sufficient to characterize genetic diversity in 15 grassland species [72].
An important limitation of DNA fragment-based methods is throughput. Throughput is defined as the amount of data that is produced simultaneously by a method. Typically, 5–10 polymorphic RAPD fragments or 20–30 AFLP can be reliably scored per primer combination used, while SSR primers yield from one to as many fragments (i.e., alleles) as the ploidy number of the species being analysed. Producing consistent fragment patterns across samples and laboratories is a critical issue for RAPD, and to some degree, for the other dominant marker systems [73]. While reproducibility is not a significant issue for SSR, SSR primers are usually species-specific and can normally be applied only on a single or a few related species (e.g., the Lolium-Festuca complex) [74]. Another limitation of fragment-based methods is their limited multiplexing potential, which reflects the maximum number of plants or markers that can be analysed in a single run. DNA fragments are usually detected by gel or capillary electrophoresis on instruments with up to 96 capillaries, in each of which up to five differently labelled samples can be analysed simultaneously. Bulked DNA approaches can further increase the multiplexing level by up to twenty plants per sample [33,45,56,62].
Despite the above-mentioned limitations, DNA fragment analysis methods offer a cost-effective alternative for genetic diversity assessments. Such methods are flexible enough to be applied in basic laboratory settings, for example, in field or breeding stations. They are also useful pilot-study approaches that set the guidelines for higher throughput analyses, which can detect genetic differentiation at a higher resolution.

3. Hybridization Array Methods

DNA hybridization arrays are high throughput methods that rely on the hybridization of labelled target DNA to a matrix of probes, which are bound to a solid surface. After removing non-hybridized target DNA, successful hybridizations are visualized by fluorescence or chemiluminescence-based methods. DNA hybridization arrays have a considerably larger throughput than fragment-based, allowing many thousands of loci to be assessed simultaneously. They have been widely applied in crop species to assess gene expression, for gene mapping or for analysing genomic loci related to quantitative traits (i.e., Quantitative Trait Loci, QTL) [75]. Hybridization array-based methods, such as SNP arrays and the Diversity Array Technology (DArT, [76,77]), have been used to measure genetic variation in grassland plant populations (Table 2). Other hybridization-based, sequencing-independent methods, such as competitive allele-specific PCR (KASP®) [78] or TaqMan® assays, have been mainly used as genotyping tools for trait selection studies and have limited potential for diversity analyses [79].

3.1. SNP Arrays

Single nucleotide polymorphisms (SNPs) are co-dominantly inherited mutations that occur at the single-base level in coding or non-coding regions of the genome [73,80]. SNPs are the most common kind of polymorphism and are identified by comparing the sequences of multiple plants from a population. They can be used to study the genetic diversity that is linked to gene function and phenotypic traits (i.e., functional genetic diversity), which may have direct applications in plant breeding [79].
SNP arrays consist of a set of SNP-containing oligonucleotides that are bound to a surface in a grid. The throughput of SNP arrays can reach hundreds of thousands of SNPs [79]. However, SNP arrays are available for only a few grassland plant species, mainly due to their high development costs. For example, an array with 2185 informative SNPs was developed and used to assess the genetic diversification of 90 ecotype accessions of perennial ryegrass from different European geographic regions [81]. Partitioning of the genetic variation showed the greatest amount of diversity to be attributed to individuals within accessions. In addition, the study found a strong differentiation of allele frequencies according to geographical location [81]. The same array was used to show that amenity and forage cultivars of perennial ryegrass are genetically differentiated from each other and that both represent a limited amount of the genetic diversity of European perennial ryegrass ecotypes [82]. In alfalfa, a 7476-SNP array was used to assess genetic differentiation of 280 genotypes, representing four cultivated subspecies of different ploidy levels, as well as related species. The study found strong differentiation according to alfalfa subspecies, ploidy level, and fall dormancy level [83].
An important limitation of SNP arrays for within-species diversity assessments is that they target a fixed set of variants, which may not be useful for applications where rare alleles are to be detected [75]. Such a bias towards frequent alleles is a kind of ascertainment bias and can be overcome to some extent by increasing the amount of SNPs on an array [84]. A reference genome can also be used to map SNPs and produce haplotypes, thus reducing the effect of ascertainment bias. Such a strategy was used to analyze the genetic relationships of 4506 wheat accessions from all over the world with a 280,226-SNP array, revealing three large gene pools that group most of the wheat diversity: an Eastern European and Mediterranean pool, an Asian pool, and a Western European pool [85]. Nevertheless, high-density SNP arrays depend on expensive oligonucleotide synthesis procedures, which makes them accessible only for major crop plants and substantially limits their use in grassland species where comprehensive genome information is often lacking.

3.2. Diversity Arrays (DArT)

In contrast to SNP arrays, diversity arrays (DArT) do not require a priori sequence information, which greatly reduces their development costs. DArT arrays are produced using “genomic representations” of the populations under analysis. Such genomic representations are generated by cutting the pooled DNA of multiple plants with restriction enzymes and then enriching certain fragments using selective PCR primers. The fragments are then cloned into a plasmid library and subsequently tethered to a surface (e.g., microscope slides) [76].
DArT arrays are particularly useful for organisms in which SNP information is scarce. For example, a DArT array based on 7–40 genotypes each of five species of the Festuca-Lolium complex produced 3884 polymorphic markers and was used to analyse the genetic relationships among a group of 184 samples of two ryegrass species (L. perenne and L. multiflorum) and three fescues (F. pratensis, F. arundinacea., and F. glaucescens Hegetschw.) [86]. The Festuca-Lolium DArT array (DArTFest) was able to differentiate each species and also capture the genetic diversity in each of those species, finding that fescues have lower levels of polymorphism than ryegrasses, which is in accordance to DNA fragment-based studies [32,42]. Low levels of genetic diversity were also found in turf-type F. arundinacea using the same array [87]. The DArTFest array has also been used for genetic mapping and to study the genetic constitution of Festuca-Lolium hybrids [88,89]. In crop plants, diversity arrays have been widely used to study genetic diversity in barley, sorghum and wheat, three cereal species with complex genomes [90,91,92].
Major drawbacks of diversity arrays include locus redundancy and ascertainment bias. Diversity arrays are made with DNA probes of several hundred base pairs in length [84,93]. Such probes are not sequence-characterized and usually show some degree of redundancy (i.e., the same fragment or portions of it are represented more than once) [84]. Ascertainment bias is also present in diversity arrays since they are produced from a reduced sample of the total diversity of a species and still can miss rare alleles, which can impact diversity estimates [84].

4. Assessing Genetic Diversity with Sequencing-Based Methods

Over the past two decades, DNA sequencing technologies have greatly increased in throughput. The output of high-throughput sequencing (HTS) technologies now reaches more than one terabase (1012 bp) per run [94], a large difference when compared to the ~105 bases of a 96-capillary Sanger sequencing run. In turn, the costs to obtain a megabase (106 bp) of a sequence have declined from several thousand US dollars in the early 2000s to less than ten cents since 2011 [95]. Such technological advances have enabled new approaches to assess the genetic diversity of plant populations with sequence-level precision. HTS technologies allow for larger multiplexing levels, which greatly increases the number of plants that are assessed simultaneously, especially when combined with bulked DNA approaches. Furthermore, ascertainment bias is reduced with HTS approaches, as SNPs are directly discovered on the sequences obtained, which enables rare SNPs to be identified.

4.1. Plastid Genome Skimming

In plant cells, DNA is mainly contained in the nucleus. However, organelles such as mitochondria or plastids also contain DNA that codes for proteins involved in processes to maintain essential cell functions or photosynthesis. In contrast to the nuclear genome, plastid genomes are circular molecules of ~150 kb that harbour ~100–120 genes [106]. The genes in plastid genomes are under high selective pressure and do not usually undergo sexual recombination, so they tend to accumulate mutations at a slower rate than the nuclear genome and are therefore highly conserved among plants, with few exceptions [106]. The high conservation of plastid genes favours their use as phylogenetic markers and for species identification (Table 2; [106]). Furthermore, organellar genomes are overrepresented several times in a cell when compared to the nuclear genome. Such overrepresentation of organelle genomes is the basis of plastid genome skimming, the recovery of plastid genomes from shotgun sequencing of plant DNA (Figure 1, green arrows).
The analysis of entire plastid genomes has been boosted by plastid genome skimming approaches based on HTS. The availability of plastid genomes has, in turn, enhanced phylogenetic resolution in many groups of plants, including the Pooideae, a group that includes all major cereal crops and forage grasses [107]. Plastid genomes also enhance species-level taxonomic identifications and can be used to produce the information needed to develop SSR [106,108]. In the model grass species Brachypodium distachyon (L.) Beauv., a total of 298 polymorphic sites and 32 haplotypes were found in 53 accessions, allowing for genetic differentiation among major clades [100]. To our knowledge, the utility of plastid genome skimming for genetic diversity assessments in forages species has not yet been thoroughly assessed.

4.2. Whole Genome Re-Sequencing

As sequencing costs per megabase keep dropping, analyzing genetic variation at the whole-genome level is becoming a reality. Once a high-quality reference genome is available for a species, additional genomes can be resequenced at a shallower depth, mapped to the reference genome and have millions of SNPs called on them [109]. Whole genome re-sequencing (WGR; purple arrows in Figure 1) was applied to assess the genomic diversity of Arabidopsis halleri L., an outcrossing plant species using a bulked DNA approach also known as Pool-Seq [97], in which the DNA from multiple genotypes is equimolarly pooled and sequenced [110,111]. Genome-wide SNPs had a good correlation with microsatellite allele richness, although SNPs are more evenly distributed and are likely to show a lower ascertainment bias. Such an even distribution of SNPs makes it possible that a few thousand randomly chosen SNPs can still reliably estimate genetic diversity and differentiation patterns [97].
An important feature of WGR is that it enables the study of pangenomes, i.e., the set of genes that exist in a species. A typical reference genome is sequenced from a single individual and, thus, contains a very narrow set of the total genomic variation that can be found in a species. Pangenome sequencing aims at complementing the reference genome of a species with genome sequences from additional individuals of the same or closely related species, sequenced at an intermediate coverage [112]. Pangenome analyses greatly increase the discovery of associations between agronomical traits and previously uncharacterized genes [112]. Furthermore, pangenomics also allows the characterization of other kinds of genomic variation, such as genomic duplications, deletions, or translocations [112]. For example, in wheat, the core genome (i.e., the genes shared by all specimens sequenced) of 18 wheat cultivars spanned 81,070 ± 1631 genes, while its pangenome spanned 140,500 ± 102 genes [113]. The variable fraction of the pangenome was enriched in genes related to agronomically relevant traits. In addition, 2.87 million SNPs were identified in the pangenome but not in the Chinese Spring wheat reference genome. In B. distachyon, 45% of the pangenome of 54 lines were not represented in the reference genome [102].
Pangenome approaches would be particularly useful to assess genetic diversity in grassland plant species, which are highly diverse and are distributed throughout contrasting geographical regions. Although the genomes of grassland plant species are large and will require a large amount of sequencing capacity, some reference genomes are already available for such species (e.g., red clover [114,115], perennial ryegrass [116]), and Italian ryegrass [117], which can help in assembling additional genomes with low coverage data. Furthermore, pangenomic information is becoming available for crop species with similarly large genomes, such as wheat, thanks to the joint efforts of international consortia (e.g., the International Wheat Genome Sequencing Consortium, [118]).

4.3. Reduced Representation Libraries

Reducing the complexity of an entire genome to a set of loci is a useful strategy to assess genetic diversity, particularly for species with a large genome size like many forage grasses and legumes. Such complexity reduction may be achieved through reduced representation libraries (RRLs; red arrows in Figure 1). RRLs are produced in a similar way as the genomic representations used for DArT (see above), except that, instead of being cloned into a plasmid library, the resulting fragments are used to produce next-generation sequencing libraries by PCR amplification. Such fragments, also known as restriction associated DNA-tags, are short sequences that are adjacent to a restriction enzyme cutting site. RRLs are at the core of methods, such as RADseq (restriction-site associated DNA sequencing [119]) or GBS (genotyping-by-sequencing [120]), among many others, which can produce tens of thousands of SNPs at very low cost [121].
In crop species, RRLs have mainly been used for marker-trait association analyses and genotyping [79]; however, applications of RRLs to study the genetic diversity of plant populations have been reported as well (Table 2). For example, GBS has been applied to assess the genetic diversity of natural populations of the self-compatible, model grass species B. distachyon [101]. Individual GBS libraries were produced from accessions of B. distachyon from different European countries, resulting in >50,000 high quality SNPs, which grouped the accessions into three genetically differentiated populations. Pooled GBS, an adaptation of GBS based on bulked DNA samples, has been used to assess segregation distortion in barley by pooling 76 genotypes to construct a single GBS library [98]. The 76 samples were also used to generate individual GBS libraries for validation and SNP allele frequencies obtained with the pooled samples was an accurate estimation of the frequencies obtained with individual samples. However, pooling of equal amounts of DNA or leaf material and sufficient sequencing depth are crucial for achieving consistent results [98].
Since no previous sequence information is needed to generate RRLs, they are suitable to analyse non-model species with a high density of SNPs. For example, pooled GBS was applied to monitor temporal changes of genetic diversity in experimental populations of perennial ryegrass by using multiple pools of 40 plants to construct GBS libraries [96]. Over the course of four years, swards sown with cultivar mixtures of perennial ryegrass changed their initial composition, while no considerable changes were detected in single-cultivar swards. Pooled GBS was also used to generate Genome Wide Allele Frequency Fingerprints (GWAFFs) in eight cultivars of perennial ryegrass [122]. The GWAFFs, which were based on the frequencies of >30,000 SNPs, successfully captured the genetic differentiation among the cultivars and always clustered together the replicates of the same cultivar [122]. Another remarkable pooled GBS application enabled the genotyping of 995 F2 families of perennial ryegrass, each one composed of 200–500 plants [123]. After filtering the SNPs with low and high allele frequencies (0.02 and 0.98, respectively), 728,359 high-quality SNPs were obtained and used to assess the heritability of crown rust resistance and heading date, finding that the most accurate estimates of heritability were produced with the SNPs with the highest sequencing depth [123].
An important limitation of RRLs is the high prevalence of missing data [124]. A missing data point arises when a given SNP is not found in all samples, which complicates downstream analyses. High rates of missing data in RRLs lead to underestimation of the true genetic diversity of a population [125,126]. Nevertheless, the reproducibility of SNP calls can be enhanced by increasing sequencing depth and stringent SNP filtering [122]. Missing data can also be imputed by substituting missing values with an approximation [127]. Furthermore, although a certain ascertainment bias is also associated with RRLs, it is lower than with other methods, such as DArT arrays [84].

4.4. Sequence Capture

Sequence capture (blue arrows in Figure 1) is another complexity reduction method in which a set of short, biotinylated RNA or DNA molecules are used as probes (or baits) to hybridize with the sample DNA. After hybridization, the conjugates are pulled down with streptavidin-coated beads, resuspended and sequenced. Although also relying on the hybridization of probes and target DNA, sequence capture differs from hybridization arrays in that sequence capture hybridization occurs in an aqueous solution instead of over a surface and polymorphisms are detected by DNA sequencing [128].
Sequence capture approaches are highly customizable. Probes can be directed to specific parts of the genome (e.g., exons, intergenic regions, centromeres), so that different components of genetic diversity can be analyzed (i.e., neutral and non-neutral genetic diversity or chromosomal rearrangements). They can also be designed with enough sequence degeneracy to capture phylogenetically diverse samples [129]. Furthermore, the use of conserved orthologous genes and ultra-conserved noncoding elements as targets is a way to overcome the lack of sequence information for orphan species [130]. Sequence capture approaches commonly result in lower amounts of missing loci than RRLs [131] and has been successfully applied to study the intraspecific diversity of northern switchgrass (Panicum virginatum L.), Sabal spp. palms and Medicago spp. (Table 2).
Combining sequence capture with RRLs has been reported to reduce the amount of missing data on the latter [132,133]. In one study, RAD libraries were subjected to sequence capture, using baits to enrich a subset of 500 pre-defined RAD-tags in 288 fish samples. With this approach >16× coverage was achieved in all the target loci and samples [133]. A similar outcome was obtained by following a modified approach, achieving 99.8% recovery of 964 target loci at >4× coverage for over 90% of 96 tree samples [132].
A commonly cited disadvantage of sequence capture approaches is that they are relatively expensive [131]. However, library pooling strategies, proper bait design, and in-house bait production can help to reduce costs per sequence capture reaction [134,135]. Furthermore, reducing costs may be possible by using standardized sets of baits, with proved efficiency in a wide taxonomic spectrum.

4.5. Amplicon Sequencing

Once a polymorphic locus is identified within a species, primers can be designed to target its flanking regions to amplify and sequence it in individual samples. Amplicon sequencing is a complexity reduction method that has been applied to assess genetic diversity in various plant species (Table 2; orange arrows in Figure 1). Amplicon sequencing is a low-cost alternative for genetic diversity assessments since only standard PCR reagents are needed. However, the method also needs a discovery step to find the genomic loci that are suitable for PCR. For example, sequence capture was used as a discovery tool for conserved genomic loci in five genotypes of 16 species of forage grasses and legumes [136]. Such sequences were then used to identify loci showing large within-species variation and to design multi-species primers. Primer sequences used for other SNP-probing methods, such as the TaqMan™ assay, can also be used to generate within-species variable amplicons [137]. Under certain conditions, publicly available gene sequences can also serve as templates for primer design. For example, the sequences of the flowering time habit locus VrnH1 from two barley lines were obtained from a public repository and used for primer design [138]. In the same study, a single reference sequence of the grain protein gene HvNAM1 was also used for primer design since that gene is known to be polymorphic in lines with differing grain protein content [138]. Still, a discovery phase may not always be necessary for primer design. ISSR can be amplified from a wide range of species by using microsatellite repeats as primers. ISSR-sequencing or MIG-seq (multiplexed ISSR genotyping by sequencing) can produce ~1000 SNPs per sample, which proved to be enough to assess genetic differentiation in populations of an orchid species [105].
The capacity of PCR-based methods can be further enhanced by multiplexing more than one pair of primers in the same reaction. Such approach is commonly called “multiplex PCR” and, for example, was applied as a genotyping tool for 192 loci on 2068 DNA samples of the fish Oncorhynchus mykiss in a single Illumina HiSeq 1500 run [137]. A multiplex PCR of >3000 primer pairs yielded a set of 2655 reproducible SNPs on two pooled samples of fugu fish (Takifugu rubripes), each with the same 326 individually tagged fishes [139].

5. Outlook for Genetic Diversity Assessments in Grassland Plant Species

A wide spectrum of tools is available to assess genetic diversity in plant populations. Such tools differ in terms of resolution (i.e., the fraction of the genome that is analysed), throughput (i.e., the amount of data that is produced simultaneously by a method), and multiplexing potential (i.e., the number of plants or markers that can be assessed simultaneously). Choosing a suitable method requires a careful balance between such methodological features and the exact purpose of the assessment. Adequate sampling of plant populations is also key to reduce ascertainment bias in genetic diversity estimations.
Genetic diversity analysis can assist forage breeders in many aspects. For example, it can be used to ensure sufficient diversity in the breeding gene pool, since low levels of genetic diversity may decrease the probability of further genetic gain [17]. Genetic diversity analysis can also be applied to select genetically diverse parents to optimize heterosis [140]. Many studies have demonstrated that the genetic differentiation between geographically distant ecotypes can be addressed with DNA fragment analysis. About 15–20 SSR or 2–10 AFLP primer combinations can produce enough polymorphic fragments to capture a meaningful proportion of the genetic diversity of a forage plant population (Table 1). Nevertheless, DNA fragment analysis has a rather low throughput (i.e., usually <20 polymorphic fragments per sample) and a low multiplexing potential (i.e., pools of 20–30 plants). In contrast, low-resolution, high-throughput sequencing approaches (e.g., amplicon sequencing) can be applied to characterize 100–1000 SNPs in highly multiplexed samples (i.e., pools of >100 plants). Such a method could, therefore, be useful as a cost-efficient, routine tool to assess genetic differentiation in the starting germplasm for forage breeding.
Genetic diversity assessments based on amplicon sequencing could be applied, for example, as a monitoring tool for breeding in mixtures. Mixture breeding aims at leveraging biodiversity effects and the interactions of multiple cultivars to increase yields and the resilience of grassland systems [141]. Keeping track of the genetic diversity when assembling and testing mixtures could greatly facilitate evaluating their performance [141]. In addition, such diversity assessments could also assist strategies to breed for several species in mixtures simultaneously to optimize interaction traits essential in multi-species grasslands [11].
Cultivar registration is another area where genetic diversity assessments could support traditional phenotypic assessments. Plant breeders need to characterize new varieties in terms of distinctness, uniformity, and stability (DUS). However, ensuring the distinctness of new cultivars of outcrossing forage cultivars is particularly difficult. Within-cultivar variation is usually high, and alleles can be widely spread among cultivars, which leads to lower genetic differentiation between cultivars. DNA fragment-based genetic diversity analyses often detect low levels of differentiation between gene pools and cultivars from diverse origins [142]. High-resolution, high-throughput approaches (e.g., RRLs), in contrast, can produce tens of thousands of SNPs on multiplexed samples, which could be used to generate GWAFFs to thoroughly characterize a novel forage variety at a low cost [143].
Large-scale genetic diversity assessments can also help in setting up and monitoring genetic resource conservation programs. Valuable genetic resources are at risk of disappearing as grassland intensification rates increase through some of the grassland biodiversity hotspots of the globe [144,145,146]. Programs to protect such genetic resources are in place in every continent, either in the form of protected areas (i.e., in situ conservation) or in large collections of germplasm (i.e., ex situ conservation). Genetic diversity information would provide an insight into the effectiveness of such conservation programs. In particular, HTS-based genetic diversity assessments (e.g., WGR, RRLs, sequence capture, and some amplicon sequencing approaches) could become informative and cost-effective tools to manage in situ and ex situ genetic resource conservation efforts, as they can be applied without much previous information [17]. WGR and pangenomic approaches could be especially relevant for assessing and incorporating the diversity of wild relatives of forage crops [147]. WGR approaches may add more genomic context information to SNPs, making it possible to detect variations in linkage disequilibrium (e.g., selective sweeps). Adding more genomic context information to SNPs could, in principle, also increase the power of genome-wide association studies [148].
In conclusion, there are many promising approaches to tackle the challenge of assessing the genetic diversity of grassland plant species. There is no one-size-fits-all solution for all kinds of assessments: the methods differ in their resolution, throughput and multiplexing potential; cost optimization will ultimately depend on the scope and purpose of the assessment. HTS approaches (e.g., RRLs, sequence capture, amplicon sequencing, and whole genome re-sequencing) hold great potential for routine genetic diversity assessments, as long as they can provide information of multiple plants at a low cost. Further increases in sequencing throughput could reduce analysis costs and enable routine monitoring programs that are currently prohibitively expensive [17]. Such large-scale monitoring programs will ultimately be of great value for forage breeding and genetic resource conservation efforts.

Author Contributions

Conceptualization, M.L.-S., R.K. and B.S.; investigation, M.L.-S.; writing—original draft preparation, M.L.-S.; writing—review and editing, M.L.-S., R.K. and B.S.; project administration, R.K. and B.S.; funding acquisition, R.K. and B.S.

Funding

This work was funded by the Swiss Federal Office of Agriculture (FOAG).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Reynolds, S.G. Chapter 1 Introduction. In Grasslands of the World; Suttie, J.M., Reynolds, S.G., Batello, C., Eds.; Food and Agriculture Organization of the United Nations: Rome, Italy, 2005. [Google Scholar]
  2. Bengtsson, J.; Bullock, J.M.; Egoh, B.; Everson, C.; Everson, T.; O’Connor, T.; O’Farrell, P.J.; Smith, H.G.; Lindborg, R. Grasslands-More Important for Ecosystem Services than You Might Think. Ecosphere 2019, 10, e02582. [Google Scholar] [CrossRef]
  3. Huber, S.; Huber, B.; Stahl, S.; Schmid, C.; Reisch, C. Species Diversity of Remnant Calcareous Grasslands in South Eastern Germany Depends on Litter Cover and Landscape Structure. Acta Oecol. 2017, 83, 48–55. [Google Scholar] [CrossRef]
  4. Tälle, M.; Deák, B.; Poschlod, P.; Valkó, O.; Westerberg, L.; Milberg, P. Grazing vs. Mowing: A Meta-Analysis of Biodiversity Benefits for Grassland Management. Agric. Ecosyst. Environ. 2016, 222, 200–212. [Google Scholar] [CrossRef]
  5. Carbutt, C.; Henwood, W.D.; Gilfedder, L.A. Global Plight of Native Temperate Grasslands: Going, Going, Gone? Biodivers. Conserv. 2017, 26, 2911–2932. [Google Scholar] [CrossRef]
  6. Wesche, K.; Krause, B.; Culmsee, H.; Leuschner, C. Fifty Years of Change in Central European Grassland Vegetation: Large Losses in Species Richness and Animal-Pollinated Plants. Biol. Conserv. 2012, 150, 76–85. [Google Scholar] [CrossRef]
  7. Feßel, C.; Meier, I.C.; Leuschner, C. Relationship between Species Diversity, Biomass and Light Transmittance in Temperate Semi-Natural Grasslands: Is Productivity Enhanced by Complementary Light Capture? J. Veg. Sci. 2016, 27, 144–155. [Google Scholar] [CrossRef]
  8. Mariotte, P.; Buttler, A.; Kohler, F.; Gilgen, A.K.; Spiegelberger, T. How Do Subordinate and Dominant Species in Semi-Natural Mountain Grasslands Relate to Productivity and Land-Use Change? Basic Appl. Ecol. 2013, 14, 217–224. [Google Scholar] [CrossRef] [Green Version]
  9. Mimura, M.; Yahara, T.; Faith, D.P.; Vázquez-Domínguez, E.; Colautti, R.I.; Araki, H.; Javadi, F.; Núñez-Farfán, J.; Mori, A.S.; Zhou, S.; et al. Understanding and Monitoring the Consequences of Human Impacts on Intraspecific Variation. Evol. Appl. 2017, 10, 121–139. [Google Scholar] [CrossRef]
  10. Barrett, B.A.; Faville, M.J.; Nichols, S.N.; Simpson, W.R.; Bryan, G.T.; Conner, A.J. Breaking through the Feed Barrier: Options for Improving Forage Genetics. Anim. Prod. Sci. 2015, 55, 883–892. [Google Scholar] [CrossRef] [Green Version]
  11. Litrico, I.; Violle, C. Diversity in Plant Breeding: A New Conceptual Framework. Trends Plant Sci. 2015, 20, 604–613. [Google Scholar] [CrossRef] [Green Version]
  12. Boller, B.; Peter-Schmid, M.K.I.; Tresch, E.; Tanner, P.; Schubiger, F.X. Ecotypes of Italian Ryegrass from Swiss Permanent Grassland Outperform Current Recommended Cultivars. Euphytica 2009, 170, 53–65. [Google Scholar] [CrossRef]
  13. Schröder, R.; Prasse, R. Do Cultivated Varieties of Native Plants Have the Ability to Outperform Their Wild Relatives? PLoS ONE 2013, 8. [Google Scholar] [CrossRef] [Green Version]
  14. Malyshev, A.V.; Arfin Khan, M.A.S.; Beierkuhnlein, C.; Steinbauer, M.J.; Henry, H.A.L.; Jentsch, A.; Dengler, J.; Willner, E.; Kreyling, J. Plant Responses to Climatic Extremes: Within-Species Variation Equals among-Species Variation. Glob. Chang. Biol. 2016, 22, 449–464. [Google Scholar] [CrossRef] [PubMed]
  15. Prieto, I.; Violle, C.; Barre, P.; Durand, J.-L.; Ghesquiere, M.; Litrico, I. Complementary Effects of Species and Genetic Diversity on Productivity and Stability of Sown Grasslands. Nat. Plants 2015, 1, 15033. [Google Scholar] [CrossRef] [PubMed]
  16. Meilhac, J.; Durand, J.-L.; Beguier, V.; Litrico, I. Increasing the Benefits of Species Diversity in Multispecies Temporary Grasslands by Increasing Within-Species Diversity. Ann. Bot. 2019, 123, 891–900. [Google Scholar] [CrossRef] [PubMed]
  17. Mascher, M.; Schreiber, M.; Scholz, U.; Graner, A.; Reif, J.C.; Stein, N. Genebank Genomics Bridges the Gap between the Conservation of Crop Diversity and Plant Breeding. Nat. Genet. 2019, 51, 1076–1081. [Google Scholar] [CrossRef]
  18. Völler, E.; Bossdorf, O.; Prati, D.; Auge, H. Evolutionary Responses to Land Use in Eight Common Grassland Plants. J. Ecol. 2017, 105, 1290–1297. [Google Scholar] [CrossRef] [Green Version]
  19. Jamali, S.H.; Cockram, J.; Hickey, L.T. Insights into Deployment of DNA Markers in Plant Variety Protection and Registration. Theor. Appl. Genet. 2019, 132, 1911–1929. [Google Scholar] [CrossRef]
  20. Botstein, D.; White, R.L.; Skolnick, M.; Davis, R.W. Construction of a Genetic Linkage Map in Man Using Restriction Fragment Length Polymorphisms. Am. J. Hum. Genet. 1980, 32, 314–331. [Google Scholar]
  21. Williams, J.G.K.; Kubelik, A.R.; Livak, K.J.; Rafalski, J.A.; Tingey, S.V. DNA Polymorphisms Amplified by Arbitrary Primers Are Useful as Genetic Markers. Nucleic Acids Res. 1990, 18, 6531–6535. [Google Scholar] [CrossRef] [Green Version]
  22. Vos, P.; Hogers, R.; Bleeker, M.; Reijans, M.; Lee, T.V.D.; Hornes, M.; Friters, A.; Pot, J.; Paleman, J.; Kuiper, M.; et al. AFLP: A New Technique for DNA Fingerprinting. Nucleic Acids Res. 1995, 23, 4407–4414. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Li, G.; Quiros, C.F. Sequence-Related Amplified Polymorphism (SRAP), a New Marker System Based on a Simple PCR Reaction: Its Application to Mapping and Gene Tagging in Brassica. Theor. Appl. Genet. 2001, 103, 455–461. [Google Scholar] [CrossRef]
  24. Kalendar, R.; Antonius, K.; Smýkal, P.; Schulman, A.H. IPBS: A Universal Method for DNA Fingerprinting and Retrotransposon Isolation. Theor. Appl. Genet. 2010, 121, 1419–1430. [Google Scholar] [CrossRef] [PubMed]
  25. Morgante, M.; Olivieri, A.M. PCR-Amplified Microsatellites as Markers in Plant Genetics. Plant J. 1993, 3, 175–182. [Google Scholar] [CrossRef] [PubMed]
  26. Zietkiewicz, E.; Rafalski, A.; Labuda, D. Genome Fingerprinting by Simple Sequence Repeat (SSR)-Anchored Polymerase Chain Reaction Amplification. Genomics 1994, 20, 176–183. [Google Scholar] [CrossRef]
  27. Last, L.; Lüscher, G.; Widmer, F.; Boller, B.; Kölliker, R. Indicators for Genetic and Phenotypic Diversity of Dactylis glomerata in Swiss Permanent Grassland. Ecol. Indic. 2014, 38, 181–191. [Google Scholar] [CrossRef]
  28. Last, L.; Widmer, F.; Fjellstad, W.; Stoyanova, S.; Kölliker, R. Genetic Diversity of Natural Orchardgrass (Dactylis glomerata L.) Populations in Three Regions in Europe. BMC Genet. 2013, 14, 102. [Google Scholar] [CrossRef] [Green Version]
  29. Xie, W.G.; Lu, X.F.; Zhang, X.Q.; Huang, L.K.; Cheng, L. Genetic Variation and Comparison of Orchardgrass (Dactylis glomerata L.) Cultivars and Wild Accessions as Revealed by SSR Markers. Genet. Mol. Res. 2012, 11, 425–433. [Google Scholar] [CrossRef]
  30. Peng, Y.A.N.; Zhang, X.; Deng, Y.; Ma, X. Evaluation of Genetic Diversity in Wild Orchardgrass (Dactylis glomerata L.) Based on AFLP Markers. Hereditas 2008, 145, 174–181. [Google Scholar] [CrossRef]
  31. Costa, R.; Pereira, G.; Garrido, I.; Tavares-De-Sousa, M.M.; Espinosa, F. Comparison of RAPD, ISSR, and AFLP Molecular Markers to Reveal and Classify Orchardgrass (Dactylis glomerata L.) Germplasm Variations. PLoS ONE 2016, 11, e0152972. [Google Scholar] [CrossRef] [Green Version]
  32. Kölliker, R.; Stadelmann, F.J.; Reidy, B.; Nösberger, J. Genetic Variability of Forage Grass Cultivars: A Comparison of Festuca pratensis Huds., Lolium perenne L., and Dactylis glomerata L. Euphytica 1999, 106, 261–270. [Google Scholar] [CrossRef]
  33. Zhang, C.; Sun, M.; Zhang, X.; Chen, S.; Nie, G.; Peng, Y.; Huang, L.; Ma, X. AFLP-Based Genetic Diversity of Wild Orchardgrass Germplasm Collections from Central Asia and Western China, and the Relation to Environmental Factors. PLoS ONE 2018, 13, e0195273. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Zeng, B.; Zhang, X.; Lan, Y.; Yang, W.-Y. Evaluation of Genetic Diversity and Relationships in Orchardgrass (Dactylis glomerata L.) Germplasm Based on SRAP Markers. Can. J. Plant Sci. 2008, 88, 53–60. [Google Scholar] [CrossRef] [Green Version]
  35. Ghariani, S.; Elazreg, H.; Chtourou-Ghorbel, N.; Chakroun, M.; Trifi-Farah, N. Genetic Diversity Analysis in Tunisian Perennial Ryegrass Germplasm as Estimated by RAPD, ISSR, and Morpho-Agronomical Markers. Genet. Mol. Res. 2015, 14, 18523–18533. [Google Scholar] [CrossRef]
  36. Van Treuren, R.; Bas, N.; Goossens, P.J.; Jansen, J.; Van Soest, L.J.M. Genetic Diversity in Perennial Ryegrass and White Clover among Old Dutch Grasslands as Compared to Cultivars and Nature Reserves. Mol. Ecol. 2005, 14, 39–52. [Google Scholar] [CrossRef]
  37. Vieira, E.A.; Castro, C.M.; de Oliveira, A.C.; de Carvalho, F.I.F.; Zimmer, P.D.; Martins, L.F. Genetic Structure of Annual Ryegrass (Lolium multiflorum) Populations Estimated by RAPD. Sci. Agric. 2004, 61, 407–413. [Google Scholar] [CrossRef] [Green Version]
  38. Statkevičiūtė, G.; Aleliūnas, A.; Kemešytė, V.; Pašakinskienė, I.; Brazauskas, G. AFLP Analysis of Genetic Diversity in an Association Mapping Panel of Lolium perenne L. In Quantitative Traits Breeding for Multifunctional Grasslands and Turf; Sokolović, D., Huyghe, C., Radović, J., Eds.; Springer: Dordrecht, The Netherlands, 2014; pp. 285–289. [Google Scholar] [CrossRef]
  39. Cresswell, A.; Sackville Hamilton, N.R.; Roy, A.K.; Viegas, B.M.F. Use of Amplified Fragment Length Polymorphism Markers to Assess Genetic Diversity of Lolium Species from Portugal. Mol. Ecol. 2001, 10, 229–241. [Google Scholar] [CrossRef]
  40. Jonavičienė, K.; Statkevičiūtė, G.; Kemešytė, V. Genetic and Phenotypic Diversity for Drought Tolerance in Perennial Ryegrass (Lolium perenne L.). Zemdirb. -Agric. 2014, 101, 411–418. [Google Scholar] [CrossRef] [Green Version]
  41. Diekmann, K.; Hodkinson, T.R.; Barth, S. New Chloroplast Microsatellite Markers Suitable for Assessing Genetic Diversity of Lolium perenne and Other Related Grass Species. Ann. Bot. 2012, 110, 1327–1339. [Google Scholar] [CrossRef] [Green Version]
  42. Peter-Schmid, M.K.I.; Boller, B.; Kölliker, R. Habitat and Management Affect Genetic Structure of Festuca pratensis but Not Lolium multiflorum Ecotype Populations. Plant Breed. 2008, 127, 510–517. [Google Scholar] [CrossRef]
  43. Roldán-Ruiz, I.; Calsyn, E.; Gilliland, T.J.; Coll, R.; van Eijk, M.J.T.; De Loose, M. Estimating Genetic Conformity between Related Ryegrass (Lolium) Varieties. 2. AFLP Characterization. Mol. Breed. 2000, 6, 593–602. [Google Scholar] [CrossRef]
  44. Barth, S.; McGrath, S.K.; Arojju, S.K.; Hodkinson, T.R. An Irish Perennial Ryegrass Genetic Resource Collection Clearly Divides into Two Major Gene Pools. Plant Genet. Resour. 2017, 15, 269–278. [Google Scholar] [CrossRef]
  45. Liu, S.; Feuerstein, U.; Luesink, W.; Schulze, S.; Asp, T.; Studer, B.; Becker, H.C.; Dehmer, K.J. DArT, SNP, and SSR Analyses of Genetic Diversity in Lolium perenne L. Using Bulk Sampling. BMC Genet. 2018, 19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Michalski, S.G.; Durka, W.; Jentsch, A.; Kreyling, J.; Pompe, S.; Schweiger, O.; Willner, E.; Beierkuhnlein, C. Evidence for Genetic Differentiation and Divergent Selection in an Autotetraploid Forage Grass (Arrhenatherum elatius). Theor. Appl. Genet. 2010, 120, 1151–1162. [Google Scholar] [CrossRef] [PubMed]
  47. Fjellheim, S.; Rognli, O.A. Molecular Diversity of Local Norwegian Meadow Fescue (Festuca pratensis Huds.) Populations and Nordic Cultivars-Consequences for Management and Utilisation. Theor. Appl. Genet. 2005, 111, 640–650. [Google Scholar] [CrossRef] [PubMed]
  48. Fedorenko, O.M.; Gritskikh, M.V.; Malysheva, I.E.; Nikolaevskaya, T.S. Genetic Diversity of Insular Natural Populations of Festuca pratensis Huds.: RAPD Analysis. Russ. J. Genet. 2009, 45, 1134–1138. [Google Scholar] [CrossRef]
  49. Duncan, D.S.; Krohn, A.L.; Jackson, R.D.; Casler, M.D. Conservation Implications of the Introduction History of Meadow Fescue (Festuca pratensis Huds.) to the Driftless Area of the Upper Mississippi Valley, USA. Plant Ecol. Divers. 2015, 8, 91–99. [Google Scholar] [CrossRef]
  50. Hand, M.L.; Cogan, N.O.I.; Forster, J.W. Molecular Characterisation and Interpretation of Genetic Diversity within Globally Distributed Germplasm Collections of Tall Fescue (Festuca arundinacea Schreb.) and Meadow Fescue (F. pratensis Huds.). Theor. Appl. Genet. 2012, 124, 1127–1137. [Google Scholar] [CrossRef]
  51. Kölliker, R.; Stadelmann, F.J.; Reidy, B.; Nösberger, J. Fertilization and Defoliation Frequency Affect Genetic Diversity of Festuca pratensis Huds. in Permanent Grasslands. Mol. Ecol. 1998, 7, 1557–1567. [Google Scholar] [CrossRef]
  52. Cuyeu, R.; Rosso, B.; Pagano, E.; Soto, G.; Fox, R.; Ayub, N.D. Genetic Diversity in a World Germplasm Collection of Tall Fescue. Genet. Mol. Biol. 2013, 36, 237–242. [Google Scholar] [CrossRef] [Green Version]
  53. Sun, X.; Xie, Y.; Bi, Y.; Liu, J.; Amombo, E.; Hu, T.; Fu, J. Comparative Study of Diversity Based on Heat Tolerant-Related Morpho-Physiological Traits and Molecular Markers in Tall Fescue Accessions. Sci. Rep. 2015, 5, 18213. [Google Scholar] [CrossRef] [PubMed]
  54. Xu, W.W.; Sleper, D.A.; Krause, G.F. Genetic Diversity of Tall Fescue Germplasm Based on RFLPs. Crop Sci. 1994, 34, 246–252. [Google Scholar] [CrossRef]
  55. Tehrani, M.S.; Mardi, M.; Sahebi, J.; Catalán, P.; Díaz-Pérez, A. Genetic Diversity and Structure among Iranian Tall Fescue Populations Based on Genomic-SSR and EST-SSR Marker Analysis. Plant Syst. Evol. 2009, 282, 57–70. [Google Scholar] [CrossRef]
  56. Mian, M.A.R.; Hopkins, A.A.; Zwonitzer, J.C. Determination of Genetic Diversity in Tall Fescue with AFLP Markers. Crop Sci. 2002, 42, 944–950. [Google Scholar] [CrossRef]
  57. Fard, J.R.; Zamani, Z.; Moghaddam, M.R.F.; Kafi, M. Evaluation of Genetic Diversity among Some Genotypes of Kentucky Bluegrass by RAPD Molecular Markers. Hortic. Environ. Biotechnol. 2012, 53, 298–303. [Google Scholar] [CrossRef]
  58. Curley, J.; Jung, G. RAPD-Based Genetic Relationships in Kentucky Bluegrass: Comparison of Cultivars, Interspecific Hybrids, and Plant Introductions. Crop Sci. 2004, 44, 1299–1306. [Google Scholar] [CrossRef]
  59. Tanhuanpää, P.; Manninen, O. High SSR Diversity but Little Differentiation between Accessions of Nordic Timothy (Phleum pratense L.). Hereditas 2012, 149, 114–127. [Google Scholar] [CrossRef] [Green Version]
  60. Abraham, E.M.; Ganopoulos, I.; Giagourta, P.; Osathanunkul, M.; Bosmali, I.; Tsaftaris, A.; Papaioannou, A.; Madesis, P. Genetic Diversity of Lotus Corniculatus in Relation to Habitat Type, Species Composition and Species Diversity. Biochem. Syst. Ecol. 2015, 63, 59–67. [Google Scholar] [CrossRef]
  61. Kempf, K.; Mora-Ortiz, M.; Smith, L.M.J.; Kölliker, R.; Skøt, L. Characterization of Novel SSR Markers in Diverse Sainfoin (Onobrychis viciifolia) Germplasm. BMC Genet. 2016, 17. [Google Scholar] [CrossRef] [Green Version]
  62. Herrmann, D.; Boller, B.; Widmer, F.; Kölliker, R. Optimization of Bulked AFLP Analysis and Its Application for Exploring Diversity of Natural and Cultivated Populations of Red Clover. Genome 2005, 48, 474–486. [Google Scholar] [CrossRef]
  63. Hargreaves, S.; Maxted, N.; Hirano, R.; Abberton, M.; Skøt, L.; Ford-Lloyd, B.V. Islands as Refugia of Trifolium repens Genetic Diversity. Conserv. Genet. 2010, 11, 1317–1326. [Google Scholar] [CrossRef]
  64. Collins, R.P.; Helgadottir, A.; Frankow-Lindberg, B.E.; Skot, L.; Jones, C.; Skot, K.P. Temporal Changes in Population Genetic Diversity and Structure in Red and White Clover Grown in Three Contrasting Environments in Northern Europe. Ann. Bot. 2012, 110, 1341–1350. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Annicchiarico, P.; Carelli, M. Origin of Ladino White Clover as Inferred from Patterns of Molecular and Morphophysiological Diversity. Crop Sci. 2014, 54, 2696. [Google Scholar] [CrossRef]
  66. Grauda, D.; Avotiņš, K.; Fokina, O.; Kolodinska-Brantestam, A.; Rashal, I. Genetic Diversity of White Clover (Trifolium repens L.) from the Urban Area of Rīga Baltā Āboliņa. Proc. Latv. Acad. Sci. Sect. B. Nat. Exactappl. Sci. 2015, 69, 132–134. [Google Scholar] [CrossRef] [Green Version]
  67. George, J.; Dobrowolski, M.P.; van Zijll de Jong, E.; Cogan, N.O.; Smith, K.F.; Forster, J.W. Assessment of Genetic Diversity in Cultivars of White Clover (Trifolium repens L.) Detected by SSR Polymorphisms. Genome 2006, 49, 919–930. [Google Scholar] [CrossRef] [PubMed]
  68. Crossa, J. Methodologies for Estimating the Sample Size Required for Genetic Conservation of Outbreeding Crops. Theor. Appl. Genet. 1989, 77, 153–161. [Google Scholar] [CrossRef] [PubMed]
  69. Darbyshire, S.J. Realignment of Festuca Subgenus Schedonorus with the Genus Lolium (Poaceae). Novon 1993, 3, 239. [Google Scholar] [CrossRef]
  70. Bolaric, S.; Barth, S.; Melchinger, A.E.; Posselt, U.K. Genetic Diversity in European Perennial Ryegrass Cultivars Investigated with RAPD Markers. Plant Breed. 2005, 124, 161–166. [Google Scholar] [CrossRef]
  71. UPOV. Guidelines for the Conduct of Tests for Distinctness, Uniformity and Stability for Ryegrass; International Union for the Protection of New Varieties of Plants: Geneva, Switzerland, 2006. [Google Scholar]
  72. Leipold, M.; Tausch, S.; Hirtreiter, M.; Poschlod, P.; Reisch, C. Sampling for Conservation Genetics: How Many Loci and Individuals Are Needed to Determine the Genetic Diversity of Plant Populations Using AFLP? Conserv. Genet. Resour. 2018, 1–10. [Google Scholar] [CrossRef]
  73. Nybom, H.; Weising, K.; Rotter, B. DNA Fingerprinting in Botany: Past, Present, Future. Investig. Genet. 2014, 5, 1. [Google Scholar] [CrossRef] [Green Version]
  74. Saha, M.C.; Cooper, J.D.; Mian, M.A.R.; Chekhovskiy, K.; May, G.D. Tall Fescue Genomic SSR Markers: Development and Transferability across Multiple Grass Species. Theor. Appl. Genet. 2006, 113, 1449–1458. [Google Scholar] [CrossRef] [PubMed]
  75. Voss-Fels, K.; Snowdon, R.J. Understanding and Utilizing Crop Genome Diversity via High-Resolution Genotyping. Plant Biotechnol. J. 2016, 14, 1086–1094. [Google Scholar] [CrossRef] [PubMed]
  76. Jaccoud, D.; Peng, K.; Feinstein, D.; Kilian, A. Diversity Arrays: A Solid State Technology for Sequence Information Independent Genotyping. Nucleic Acids Res. 2001, 29, e25. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Kilian, A.; Wenzl, P.; Huttner, E.; Carling, J.; Xia, L.; Blois, H.; Caig, V.; Heller-Uszynska, K.; Jaccoud, D.; Hopper, C.; et al. Diversity Arrays Technology: A Generic Genome Profiling Technology on Open Platforms. In Data Production and Analysis in Population Genomics: Methods and Protocols; Pompanon, F., Bonin, A., Eds.; Humana Press: Totowa, NJ, USA, 2012; pp. 67–89. [Google Scholar] [CrossRef]
  78. Smith, S.M.; Maughan, P.J. SNP Genotyping Using KASPar Assays. In Plant Genotyping: Methods and Protocols; Batley, J., Ed.; Springer: New York, NY, USA, 2015; pp. 243–256. [Google Scholar] [CrossRef]
  79. Thomson, M.J. High-Throughput SNP Genotyping to Accelerate Crop Improvement. Plant Breed. Biotechnol. 2014, 2, 195–212. [Google Scholar] [CrossRef]
  80. Gupta, P.K.; Rustgi, S.; Mir, R.R. Array-Based High-Throughput DNA Markers for Crop Improvement. Heredity 2008, 101, 5–18. [Google Scholar] [CrossRef] [Green Version]
  81. Blackmore, T.; Thomas, I.; McMahon, R.; Powell, W.; Hegarty, M. Genetic-Geographic Correlation Revealed across a Broad European Ecotypic Sample of Perennial Ryegrass (Lolium perenne) Using Array-Based SNP Genotyping. Theor. Appl. Genet. 2015, 128, 1917–1932. [Google Scholar] [CrossRef] [Green Version]
  82. Blackmore, T.; Thorogood, D.; Skøt, L.; McMahon, R.; Powell, W.; Hegarty, M. Germplasm Dynamics: The Role of Ecotypic Diversity in Shaping the Patterns of Genetic Variation in Lolium perenne. Sci. Rep. 2016, 6, 22603. [Google Scholar] [CrossRef]
  83. Li, X.; Han, Y.; Wei, Y.; Acharya, A.; Farmer, A.D.; Ho, J.; Monteros, M.J.; Brummer, E.C. Development of an Alfalfa SNP Array and Its Use to Evaluate Patterns of Population Structure and Linkage Disequilibrium. PLoS ONE 2014, 9, e84329. [Google Scholar] [CrossRef] [Green Version]
  84. Heslot, N.; Rutkoski, J.; Poland, J.; Jannink, J.-L.; Sorrells, M.E. Impact of Marker Ascertainment Bias on Genomic Selection Accuracy and Estimates of Genetic Diversity. PLoS ONE 2013, 8, e74612. [Google Scholar] [CrossRef] [Green Version]
  85. Balfourier, F.; Bouchet, S.; Robert, S.; DeOliveira, R.; Rimbert, H.; Kitt, J.; Choulet, F.; Paux, E. Worldwide Phylogeography and History of Wheat Genetic Diversity. Sci. Adv. 2019, 5. [Google Scholar] [CrossRef] [Green Version]
  86. Kopecký, D.; Bartoš, J.; Lukaszewski, A.J.; Baird, J.H.; Černoch, V.; Kölliker, R.; Rognli, O.; Blois, H.; Caig, V.; Lübberstedt, T.; et al. Development and Mapping of DArT Markers within the Festuca-Lolium Complex. BMC Genom. 2009, 10, 473. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Baird, J.H.; Kopecký, D.; Lukaszewski, A.J.; Green, R.L.; Bartoš, J.; Doležel, J. Genetic Diversity of Turf-Type Tall Fescue Using Diversity Arrays Technology. Crop Sci. 2012, 52, 408. [Google Scholar] [CrossRef]
  88. Bartoš, J.; Sandve, S.R.; Kölliker, R.; Kopecký, D.; Christelová, P.; Stočes, Š.; Østrem, L.; Larsen, A.; Kilian, A.; Rognli, O.A.; et al. Genetic Mapping of DArT Markers in the Festuca-Lolium Complex and Their Use in Freezing Tolerance Association Analysis. Theor. Appl. Genet. 2011, 122, 1133–1147. [Google Scholar] [CrossRef] [PubMed]
  89. Kopecký, D.; Bartoš, J.; Christelová, P.; Černoch, V.; Kilian, A.; Doležel, J. Genomic Constitution of Festuca × Lolium Hybrids Revealed by the DArTFest Array. Theor. Appl. Genet. 2011, 122, 355–363. [Google Scholar] [CrossRef]
  90. Akbari, M.; Wenzl, P.; Caig, V.; Carling, J.; Xia, L.; Yang, S.; Uszynski, G.; Mohler, V.; Lehmensiek, A.; Kuchel, H.; et al. Diversity Arrays Technology (DArT) for High-Throughput Profiling of the Hexaploid Wheat Genome. Theor. Appl. Genet. 2006, 113, 1409–1420. [Google Scholar] [CrossRef]
  91. Mace, E.S.; Xia, L.; Jordan, D.R.; Halloran, K.; Parh, D.K.; Huttner, E.; Wenzl, P.; Kilian, A. DArT Markers: Diversity Analyses and Mapping in Sorghum Bicolor. BMC Genom. 2008, 9. [Google Scholar] [CrossRef] [Green Version]
  92. Zhang, L.Y.; Marchand, S.; Tinker, N.A.; Belzile, F. Population Structure and Linkage Disequilibrium in Barley Assessed by DArT Markers. Theor. Appl. Genet. 2009, 119, 43–52. [Google Scholar] [CrossRef]
  93. Wenzl, P.; Carling, J.; Kudrna, D.; Jaccoud, D.; Huttner, E.; Kleinhofs, A.; Kilian, A. Diversity Arrays Technology (DArT) for Whole-Genome Profiling of Barley. Proc. Natl. Acad. Sci. USA 2004, 101, 9915–9920. [Google Scholar] [CrossRef] [Green Version]
  94. Shendure, J.; Balasubramanian, S.; Church, G.M.; Gilbert, W.; Rogers, J.; Schloss, J.A.; Waterston, R.H. DNA Sequencing at 40: Past, Present and Future. Nature 2017, 550, 345–353. [Google Scholar] [CrossRef]
  95. Wetterstrand, K.A. DNA Sequencing Costs: Data from the NHGRI Genome Sequencing Program (GSP). Available online: https://www.genome.gov/about-genomics/fact-sheets/DNA-Sequencing-Costs-Data (accessed on 8 October 2019).
  96. Verwimp, C.; Ruttink, T.; Muylle, H.; Van Glabeke, S.; Cnops, G.; Quataert, P.; Honnay, O.; Roldán-Ruiz, I. Temporal Changes in Genetic Diversity and Forage Yield of Perennial Ryegrass in Monoculture and in Combination with Red Clover in Swards. PLoS ONE 2018, 13, e0206571. [Google Scholar] [CrossRef]
  97. Fischer, M.C.; Rellstab, C.; Leuzinger, M.; Roumet, M.; Gugerli, F.; Shimizu, K.K.; Holderegger, R.; Widmer, A. Estimating Genomic Diversity and Population Differentiation—An Empirical Comparison of Microsatellite and SNP Variation in Arabidopsis halleri. BMC Genom. 2017, 18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  98. Bélanger, S.; Esteves, P.; Clermont, I.; Jean, M.; Belzile, F.; Belzile, F. Genotyping-by-Sequencing on Pooled Samples and Its Use in Measuring Segregation Bias during the Course of Androgenesis in Barley. Plant Genome 2016, 9. [Google Scholar] [CrossRef] [Green Version]
  99. Bélanger, S.; Clermont, I.; Esteves, P.; Belzile, F. Extent and Overlap of Segregation Distortion Regions in 12 Barley Crosses Determined via a Pool-GBS Approach. Theor. Appl. Genet. 2016, 129, 1393–1404. [Google Scholar] [CrossRef] [PubMed]
  100. Sancho, R.; Cantalapiedra, C.P.; López-Alvarez, D.; Gordon, S.P.; Vogel, J.P.; Catalán, P.; Contreras-Moreira, B. Comparative Plastome Genomics and Phylogenomics of Brachypodium: Flowering Time Signatures, Introgression and Recombination in Recently Diverged Ecotypes. New Phytol. 2018, 218, 1631–1644. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  101. Tyler, L.; Lee, S.J.; Young, N.D.; DeIulio, G.A.; Benavente, E.; Reagon, M.; Sysopha, J.; Baldini, R.M.; Troìa, A.; Hazen, S.P.; et al. Population Structure in the Model Grass Brachypodium distachyon Is Highly Correlated with Flowering Differences across Broad Geographic Areas. Plant Genome 2016, 9. [Google Scholar] [CrossRef] [Green Version]
  102. Gordon, S.P.; Contreras-Moreira, B.; Woods, D.P.; Des Marais, D.L.; Burgess, D.; Shu, S.; Stritt, C.; Roulin, A.C.; Schackwitz, W.; Tyler, L.; et al. Extensive Gene Content Variation in the Brachypodium distachyon Pan-Genome Correlates with Population Structure. Nat. Commun. 2017, 8, 2184. [Google Scholar] [CrossRef]
  103. De Sousa, F.; Bertrand, Y.J.K.; Nylinder, S.; Oxelman, B.; Eriksson, J.S.; Pfeil, B.E. Phylogenetic Properties of 50 Nuclear Loci in Medicago (Leguminosae) Generated Using Multiplexed Sequence Capture and Next-Generation Sequencing. PLoS ONE 2014, 9, e109704. [Google Scholar] [CrossRef] [Green Version]
  104. Evans, J.; Crisovan, E.; Barry, K.; Daum, C.; Jenkins, J.; Kunde-Ramamoorthy, G.; Nandety, A.; Ngan, C.Y.; Vaillancourt, B.; Wei, C.L.; et al. Diversity and Population Structure of Northern Switchgrass as Revealed through Exome Capture Sequencing. Plant J. 2015, 84, 800–815. [Google Scholar] [CrossRef]
  105. Suyama, Y.; Matsuki, Y. MIG-Seq: An Effective PCR-Based Method for Genome-Wide Single-Nucleotide Polymorphism Genotyping Using the next-Generation Sequencing Platform. Sci. Rep. 2015, 5, 16963. [Google Scholar] [CrossRef] [Green Version]
  106. Tonti-Filippini, J.; Nevill, P.G.; Dixon, K.; Small, I. What Can We Do with 1000 Plastid Genomes? Plant J. 2017, 90, 808–818. [Google Scholar] [CrossRef] [Green Version]
  107. Orton, L.M.; Burke, S.V.; Duvall, M.R. Plastome Phylogenomics and Characterization of Rare Genomic Changes as Taxonomic Markers in Plastome Groups 1 and 2 Poeae (Pooideae; Poaceae). PeerJ 2019, 7, e6959. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  108. Gryta, H.; Van de Paer, C.; Manzi, S.; Holota, H.; Roy, M.; Besnard, G. Genome Skimming and Plastid Microsatellite Profiling of Alder Trees (Alnus Spp., Betulaceae): Phylogenetic and Phylogeographical Prospects. Tree Genet. Genomes 2017, 13, 118. [Google Scholar] [CrossRef]
  109. Fuentes-Pardo, A.P.; Ruzzante, D.E. Whole-Genome Sequencing Approaches for Conservation Biology: Advantages, Limitations and Practical Recommendations. Mol. Ecol. 2017, 26, 5369–5406. [Google Scholar] [CrossRef] [PubMed]
  110. Fischer, M.C.; Rellstab, C.; Tedder, A.; Zoller, S.; Gugerli, F.; Shimizu, K.K.; Holderegger, R.; Widmer, A. Population Genomic Footprints of Selection and Associations with Climate in Natural Populations of Arabidopsis halleri from the Alps. Mol. Ecol. 2013, 22, 5594–5607. [Google Scholar] [CrossRef] [Green Version]
  111. Rellstab, C.; Zoller, S.; Tedder, A.; Gugerli, F.; Fischer, M.C. Validation of SNP Allele Frequencies Determined by Pooled Next-Generation Sequencing in Natural Populations of a Non-Model Plant Species. PLoS ONE 2013, 8, e80422. [Google Scholar] [CrossRef]
  112. Golicz, A.A.; Batley, J.; Edwards, D. Towards Plant Pangenomics. Plant Biotechnol. J. 2016, 14, 1099–1105. [Google Scholar] [CrossRef]
  113. Montenegro, J.D.; Golicz, A.A.; Bayer, P.E.; Hurgobin, B.; Lee, H.T.; Chan, C.K.K.; Visendi, P.; Lai, K.; Doležel, J.; Batley, J.; et al. The Pangenome of Hexaploid Bread Wheat. Plant J. 2017, 90, 1007–1013. [Google Scholar] [CrossRef] [Green Version]
  114. De Vega, J.J.; Ayling, S.; Hegarty, M.; Kudrna, D.; Goicoechea, J.L.; Ergon, Å.; Rognli, O.A.; Jones, C.; Swain, M.; Geurts, R.; et al. Red Clover (Trifolium pratense L.) Draft Genome Provides a Platform for Trait Improvement. Sci. Rep. 2015, 5, 17394. [Google Scholar] [CrossRef]
  115. Ištvánek, J.; Jaroš, M.; Křenek, A.; Řepková, J. Genome Assembly and Annotation for Red Clover (Trifolium pratense; Fabaceae). Am. J. Bot. 2014, 101, 327–337. [Google Scholar] [CrossRef]
  116. Byrne, S.L.; Nagy, I.; Pfeifer, M.; Armstead, I.; Swain, S.; Studer, B.; Mayer, K.; Campbell, J.D.; Czaban, A.; Hentrup, S.; et al. A Synteny-Based Draft Genome Sequence of the Forage Grass Lolium perenne. Plant J. 2015, 84, 816–826. [Google Scholar] [CrossRef]
  117. Knorst, V.; Yates, S.; Byrne, S.; Asp, T.; Widmer, F.; Studer, B.; Kölliker, R. First Assembly of the Gene-Space of Lolium multiflorum and Comparison to Other Poaceae Genomes. Grassl. Sci. 2019, 65, 125–134. [Google Scholar] [CrossRef] [Green Version]
  118. Alaux, M.; Rogers, J.; Letellier, T.; Flores, R.; Alfama, F.; Pommier, C.; Mohellibi, N.; Durand, S.; Kimmel, E.; Michotey, C.; et al. Linking the International Wheat Genome Sequencing Consortium Bread Wheat Reference Genome Sequence to Wheat Genetic and Phenomic Data. Genome Biol. 2018, 19, 111. [Google Scholar] [CrossRef] [PubMed]
  119. Baird, N.A.; Etter, P.D.; Atwood, T.S.; Currey, M.C.; Shiver, A.L.; Lewis, Z.A.; Selker, E.U.; Cresko, W.A.; Johnson, E.A. Rapid SNP Discovery and Genetic Mapping Using Sequenced RAD Markers. PLoS ONE 2008, 3, e3376. [Google Scholar] [CrossRef] [PubMed]
  120. Elshire, R.J.; Glaubitz, J.C.; Sun, Q.; Poland, J.A.; Kawamoto, K.; Buckler, E.S.; Mitchell, S.E. A Robust, Simple Genotyping-by-Sequencing (GBS) Approach for High Diversity Species. PLoS ONE 2011, 6, e19379. [Google Scholar] [CrossRef] [Green Version]
  121. Andrews, K.R.; Good, J.M.; Miller, M.R.; Luikart, G.; Hohenlohe, P.A. Harnessing the Power of RADseq for Ecological and Evolutionary Genomics. Nat. Rev. Genet. 2016, 17, 81–92. [Google Scholar] [CrossRef] [Green Version]
  122. Byrne, S.; Czaban, A.; Studer, B.; Panitz, F.; Bendixen, C.; Asp, T. Genome Wide Allele Frequency Fingerprints (GWAFFs) of Populations via Genotyping by Sequencing. PLoS ONE 2013, 8, e57438. [Google Scholar] [CrossRef] [Green Version]
  123. Ashraf, B.H.; Byrne, S.; Fé, D.; Czaban, A.; Asp, T.; Pedersen, M.G.; Lenk, I.; Roulund, N.; Didion, T.; Jensen, C.S.; et al. Estimating Genomic Heritabilities at the Level of Family-Pool Samples of Perennial Ryegrass Using Genotyping-by-Sequencing. Theor. Appl. Genet. 2016, 129, 45–52. [Google Scholar] [CrossRef] [Green Version]
  124. Lowry, D.B.; Hoban, S.; Kelley, J.L.; Lotterhos, K.E.; Reed, L.K.; Antolin, M.F.; Storfer, A. Breaking RAD: An Evaluation of the Utility of Restriction Site-Associated DNA Sequencing for Genome Scans of Adaptation. Mol. Ecol. Resour. 2017, 17, 142–152. [Google Scholar] [CrossRef]
  125. Cariou, M.; Duret, L.; Charlat, S. How and How Much Does RAD-Seq Bias Genetic Diversity Estimates? BMC Evol. Biol. 2016, 16, 240. [Google Scholar] [CrossRef] [Green Version]
  126. Arnold, B.; Corbett-Detig, R.B.; Hartl, D.; Bomblies, K. RADseq Underestimates Diversity and Introduces Genealogical Biases Due to Nonrandom Haplotype Sampling. Mol. Ecol. 2013, 22, 3179–3190. [Google Scholar] [CrossRef]
  127. Nazzicari, N.; Biscarini, F.; Cozzi, P.; Brummer, E.C.; Annicchiarico, P. Marker Imputation Efficiency for Genotyping-by-Sequencing Data in Rice (Oryza sativa) and Alfalfa (Medicago sativa). Mol. Breed. 2016, 36, 69. [Google Scholar] [CrossRef]
  128. Gasc, C.; Peyretaillade, E.; Peyret, P. Sequence Capture by Hybridization to Explore Modern and Ancient Genomic Diversity in Model and Nonmodel Organisms. Nucleic Acids Res. 2016, 44, 4504–4518. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  129. Faircloth, B.C. Identifying conserved genomic elements and designing universal bait sets to enrich them. Methods Ecol. Evol. 2017, 8, 1103–1112. [Google Scholar] [CrossRef] [Green Version]
  130. Kadlec, M.; Bellstedt, D.U.; Le Maitre, N.C.; Pirie, M.D. Targeted NGS for Species Level Phylogenomics: “Made to Measure” or “One Size Fits All”? PeerJ 2017, 5, e3569. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  131. Harvey, M.G.; Smith, B.T.; Glenn, T.C.; Faircloth, B.C.; Brumfield, R.T. Sequence Capture versus Restriction Site Associated DNA Sequencing for Shallow Systematics. Syst. Biol. 2016, 65, 910–924. [Google Scholar] [CrossRef] [Green Version]
  132. Hoffberg, S.L.; Kieran, T.J.; Catchen, J.M.; Devault, A.; Faircloth, B.C.; Mauricio, R.; Glenn, T.C. RADcap: Sequence Capture of Dual-Digest RADseq Libraries with Identifiable Duplicates and Reduced Missing Data. Mol. Ecol. Resour. 2016, 16, 1264–1278. [Google Scholar] [CrossRef]
  133. Ali, O.A.; O’Rourke, S.M.; Amish, S.J.; Meek, M.H.; Luikart, G.; Jeffres, C.; Miller, M.R. Rad Capture (Rapture): Flexible and Efficient Sequence-Based Genotyping. Genetics 2016, 202, 389–400. [Google Scholar] [CrossRef] [Green Version]
  134. Zarza, E.; Faircloth, B.C.; Tsai, W.L.E.; Bryson, R.W.; Klicka, J.; McCormack, J.E. Hidden Histories of Gene Flow in Highland Birds Revealed with Genomic Markers. Mol. Ecol. 2016, 25, 5144–5157. [Google Scholar] [CrossRef]
  135. Peñalba, J.V.; Smith, L.L.; Tonione, M.A.; Sass, C.; Hykin, S.M.; Skipwith, P.L.; Mcguire, J.A.; Bowie, R.C.K.; Moritz, C. Sequence Capture Using PCR-Generated Probes: A Cost-Effective Method of Targeted High-Throughput Sequencing for Nonmodel Organisms. Mol. Ecol. Resour. 2014, 14, 1000–1010. [Google Scholar] [CrossRef]
  136. Loera-Sánchez, M.A.; Studer, B.; Kölliker, R. Towards a High-Throughput Method to Measure Genetic Diversity in Grassland. In Grassland Science in Europe, Vol. 24—Improving Sown Grasslands through Breeding and Management; Huguenin-Elie, O., Studer, B., Kölliker, R., Reheul, D., Probo, M., Barre, P., Feuerstein, U., Roldán-Ruiz, I., Mariotte, P., Hopkins, A., Eds.; European Grassland Federation: Zurich, Switzerland, 2019; pp. 12–19. [Google Scholar] [CrossRef]
  137. Campbell, N.R.; Harmon, S.; Narum, S.R. Genotyping-in-Thousands by Sequencing (GT-Seq): A Cost Effective SNP Genotyping Method Based on Custom Amplicon Sequencing. Mol. Ecol. Resour. 2014, 15, 855–867. [Google Scholar] [CrossRef]
  138. Fritsch, L.; Soeur, R.; Hansen, C.; Fischer, R.; Schillberg, S.; Schröper, F. Next-Generation Sequencing of Amplicons Is a Rapid and Reliable Method for the Detection of Polymorphisms Relevant for Barley Breeding. Mol. Breed. 2016, 36, 83. [Google Scholar] [CrossRef]
  139. Sato, M.; Hosoya, S.; Yoshikawa, S.; Ohki, S.; Kobayashi, Y.; Itou, T.; Kikuchi, K. A Highly Flexible and Repeatable Genotyping Method for Aquaculture Studies Based on Target Amplicon Sequencing Using Next-Generation Sequencing Technology. Sci. Rep. 2019, 9, 6904. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  140. Kölliker, R.; Boller, B.; Widmer, F. Marker Assisted Polycross Breeding to Increase Diversity and Yield in Perennial Ryegrass (Lolium perenne L.). Euphytica 2005, 146, 55–65. [Google Scholar] [CrossRef]
  141. Barot, S.; Allard, V.; Cantarel, A.; Enjalbert, J.; Gauffreteau, A.; Goldringer, I.; Lata, J.C.; Le Roux, X.; Niboyet, A.; Porcher, E. Designing Mixtures of Varieties for Multifunctional Agriculture with the Help of Ecology. A Review. Agron. Sustain. Dev. 2017, 37, 13. [Google Scholar] [CrossRef] [Green Version]
  142. Roldán-Ruiz, I.; Dendauw, J.; Van Bockstaele, E.; Depicker, A.; De Loose, M. AFLP Markers Reveal High Polymorphic Rates in Ryegrasses (Lolium spp.). Mol. Breed. 2000, 6, 125–134. [Google Scholar] [CrossRef]
  143. Annicchiarico, P.; Nazzicari, N.; Ananta, A.; Carelli, M.; Wei, Y.; Brummer, E.C. Assessment of Cultivar Distinctness in Alfalfa: A Comparison of Genotyping-by-Sequencing, Simple-Sequence Repeat Marker, and Morphophysiological Observations. Plant Genome 2016, 9. [Google Scholar] [CrossRef] [Green Version]
  144. Dengler, J.; Janišová, M.; Török, P.; Wellstein, C. Biodiversity of Palaearctic Grasslands: A Synthesis. Agric. Ecosyst. Environ. 2014, 182, 1–14. [Google Scholar] [CrossRef] [Green Version]
  145. Bartolome, J.W.; Allen-Diaz, B.H.; Barry, S.; Ford, L.D.; Hammond, M.; Hopkinson, P.; Ratcliff, F.; Spiegal, S.; White, M.D. Grazing for Biodiversity in Californian Mediterranean Grasslands. Rangelands 2014, 36, 36–43. [Google Scholar] [CrossRef] [Green Version]
  146. Graf, R.; Müller, M.; Korner, P.; Jenny, M.; Jenni, L. 20% Loss of Unimproved Farmland in 22 Years in the Engadin, Swiss Alps. Agric. Ecosyst. Environ. 2014, 185, 48–58. [Google Scholar] [CrossRef]
  147. Henry, R.J.; Nevo, E. Exploring Natural Selection to Guide Breeding for Agriculture. Plant Biotechnol. J. 2014, 12, 655–662. [Google Scholar] [CrossRef]
  148. Josephs, E.B.; Stinchcombe, J.R.; Wright, S.I. What Can Genome-Wide Association Studies Tell Us about the Evolutionary Forces Maintaining Genetic Variation for Quantitative Traits? New Phytol. 2017, 214, 21–33. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Schematic representation of high-throughput-sequencing methods for genetic diversity assessments. Complexity reduction methods are highlighted in the grey oval. The plus signs indicate the stages where library preparation is usually performed (e.g., Illumina adapter and indexing incorporation). Reduced representation libraries (red arrows) are produced by digesting the starting DNA material with restriction enzymes. Amplicon sequencing (orange arrows) begins by amplifying the targeted loci from the starting DNA. Priming sites are shown in grey. In the particular case of amplicon sequencing, library preparation can start either during PCR (i.e., using adapter-bound primers) or after PCR (i.e., through blunt-end ligation of the adapters). Sequence capture (blue arrows) starts by hybridizing modified RNA or DNA probes (yellow hexagons) to DNA fragments containing target loci; the hybridized DNA fragments are then pulled using magnetic beads that bind to the RNA or DNA probes. Whole genome re-sequencing (WGR; purple arrows) is performed by sequencing random fragments of the starting DNA. Plastid genome skimming (green arrows) recovers the plastid genome sequences from the data produced by, for example, sequence capture or WGR.
Figure 1. Schematic representation of high-throughput-sequencing methods for genetic diversity assessments. Complexity reduction methods are highlighted in the grey oval. The plus signs indicate the stages where library preparation is usually performed (e.g., Illumina adapter and indexing incorporation). Reduced representation libraries (red arrows) are produced by digesting the starting DNA material with restriction enzymes. Amplicon sequencing (orange arrows) begins by amplifying the targeted loci from the starting DNA. Priming sites are shown in grey. In the particular case of amplicon sequencing, library preparation can start either during PCR (i.e., using adapter-bound primers) or after PCR (i.e., through blunt-end ligation of the adapters). Sequence capture (blue arrows) starts by hybridizing modified RNA or DNA probes (yellow hexagons) to DNA fragments containing target loci; the hybridized DNA fragments are then pulled using magnetic beads that bind to the RNA or DNA probes. Whole genome re-sequencing (WGR; purple arrows) is performed by sequencing random fragments of the starting DNA. Plastid genome skimming (green arrows) recovers the plastid genome sequences from the data produced by, for example, sequence capture or WGR.
Agronomy 09 00881 g001
Table 1. Overview of DNA fragment analysis-based genetic diversity studies on major grassland plant species from temperate grasslands.
Table 1. Overview of DNA fragment analysis-based genetic diversity studies on major grassland plant species from temperate grasslands.
SpeciesType and Number of Populations Analysed Plants Per PopulationMarker System aRef.
Orchardgrass (Dactylis glomerata L.)Ecotypes (20)32SSR (29)[27]
Ecotypes (59)32SSR (29)[28]
Ecotypes and cultivars (6)20SSR (15)[29]
Ecotypes (3)2–20AFLP (9)[30]
Subspecies (3)2–49RAPD (26), ISSR (22), AFLP (4)[31]
Cultivars (3)28RAPD (12)[32]
Ecotypes (14)15AFLP (6)[33]
Ecotypes and cultivars (60)25SRAP (21)[34]
Perennial and Italian ryegrasses (Lolium perenne L. and L. multiflorum Lam.)Cultivars (3)28RAPD (12)[32]
Ecotypes and cultivars (8)10ISSR (10)[35]
Ecotypes and cultivars (24)36AFLP (2)[36]
Ecotypes (4)84–99RAPD (6)[37]
Ecotypes and cultivars (3)38–90AFLP (3)[38]
Species and hybrids (16)8AFLP (3)[39]
Ecotypes and cultivars (104)1AFLP (3)[40]
Ecotypes and cultivars (60)1–15SSR (9)[41]
Ecotypes and cultivars (32)23SSR (34)[42]
Breeding lines (12)52–60AFLP (56)[43]
Ecotypes (40)~24SSR (8)[44]
Gene bank accessions (297)24SSR (48)[45]
Tall oat-grass (Arrhenaterum elatius L.)Gene bank accessions (46)11AFLP (4)[46]
Meadow fescue (Festuca pratensis Huds., syn. Lolium pratense and Schenodurus pratensis)Cultivars (3)28RAPD (12)[32]
Ecotypes and cultivars (32)23SSR (34)[42]
Ecotypes and cultivars (43)14–20AFLP (3)[47]
Ecotypes (3)20RAPD (5)[48]
Ecotypes (71)~33SSR (24)[49]
Gene bank accessions (189)1SSR (34)[50]
Ecotypes and cultivars (9)36RAPD (13)[51]
Tall fescue (Festuca arundinacea Schreb. ., syn. Lolium arundinaceum and Schenodurus arundinaceus)Gene bank accessions (161)30 leavesSSR (15)[52]
Ecotypes and cultivars (100)1SSR (90)[53]
Cultivars (16)20RFLP[54]
Ecotypes (9)2–12SSR (11)[55]
Gene bank accessions (851)1SSR (28)[50]
Gene bank accessions and cultivars (18)16AFLP (6)[56]
Kentucky bluegrass (Poa pratensis L.)Ecotypes, accessions and cultivars (16)1RAPD (17)[57]
Cultivars, ecotypes and breeding lines (123)3RAPD (12)[58]
Timothy (Phleum pratense L.)Gene bank accessions (96)15–20SSR (13)[59]
Bird’s-foot trefoil (Lotus corniculatus L.)Ecotypes (8)20ISSR (7)[60]
Sainfoin (Onobrychis viciifolia Scop.)Cultivars and landraces (32)1SSR (400)[61]
Red and white clover (Trifolium pratense L. and T. repens L.)Ecotypes and cultivars (23)36AFLP (2)[36]
Ecotypes, landraces and cultivars (120)40AFLP (12)[62]
Landraces and cultivars (29)12–32AFLP (2)[63]
Cultivars and breeding lines (19)31–96AFLP (4)[64]
Ecotypes and landraces (7)25SSR (32)[65]
Ecotypes and cultivars (11)23–37iPBS (2)[66]
Cultivars (16)8–22SSR (15)[67]
a—Between parenthesis are shown: the number of simple sequence repeats (SSR) or inter SRR (ISSR) primer pairs, the number of amplified fragment length polymorphism (AFLP) or sequence-related amplified polymorphism (SRAP) primer combinations, or the number of inter-primer binding site amplification (iPBS) or random amplified polymorphic DNA (RAPD) primers analysed in each study.
Table 2. Selected applications of hybridization array and high throughput sequencing (HTS) approaches to assess genetic diversity in plants.
Table 2. Selected applications of hybridization array and high throughput sequencing (HTS) approaches to assess genetic diversity in plants.
SpeciesMethodNo. of PolymorphismsType of StudyMain Findings Related to Genetic DiversityRef.
Fescues and ryegrasses (Festuca spp. and Lolium spp.)SNP array2185 SNPsArray development and genetic diversity analysisThe array had 99% call rate and could detect population structure in European ecotypes of perennial ryegrass.
[81]
DArT7680 probesArray development and genetic diversity analysisThere were 3884 polymorphic probes on 40 genotypes of 5 species of fescues and ryegrasses. Ecotypes from the same country were generally more closely related.
[86]
DArT7680 probesGenetic diversity analysisThere were 190 polymorphic probes on 93 tall fescue accessions. Turf-type accessions had lower genetic diversity than forage-type tall fescue.
[87]
Pooled GBS22,324 SNPsGenetic diversity analysisPooled GBS was validated on pools of 40 perennial ryegrass plants. The allele frequencies of mixed-cultivar swards change over time.
[96]
DArT1384 probesGenetic diversity analysisPools of 30 plants from 297 accessions of perennial ryegrass were characterized with DArT, SNPs and SSR. DArT-based pairwise genetic distances had the highest reproducibility and consistency.
[45]
Arabidopsis halleri L.Whole genome re-sequencing2 million SNPsGenomic diversity analysisSNPs provide an unbiased estimation of genetic diversity of 180 A. halleri plants. A few thousand random SNPs are enough to assess genetic differentiation.
[97]
Barley (Hordeum vulgare L.)DArT2304 probesGenetic diversity analysis942 polymorphic probes were found in 170 barley cultivars used in Canada. In general, two-row and six-row barleys were genetically differentiated.
[92]
DArT1920 probesWhole-genome profilingThe DNA from 33 cultivars and two wild accessions was analysed. The DNA from the two wild accessions and two Southeast Asian landraces were underrepresented on the array, and they were genetically differentiated from the rest of the cultivars, which clustered according to expected relationships.
[93]
Pooled GBS1984 SNPsSegregation bias analysisPooled GBS was validated by comparing the allele frequencies from individual plants and pools of such plants. Pooled GBS could track changes of allele frequencies and segregation bias during androgenesis.
[98]
Pooled GBS674–1744 SNPsSegregation bias analysisPooled GBS allowed an unprecedented analysis of segregation bias on 12 segregating barley populations.
[99]
Brachypodium spp.Plastid genome skimming57 plastomes (not polymorphisms)Comparative genomics and phylogenomicsThe plastomes of Brachypodium spp. showed structural variations. Within Brachypodium distachyon (L.) Beauv., two main lineages with differing flowering time traits were detected.
[100]
GBS50,000 SNPsGenetic diversity analysis84 accessions of B. distachyon were genetically differentiated according to flowering time. Seven B. hybridum samples proved to be a mosaic of B. distachyon and B. stacei.
[101]
Whole genome re-sequencing~4 million SNPsGenetic diversity and pan-genome analysis54 lines of B. distachyon were resequenced at a median coverage of 94×. The resulting pangenome was 58% larger (430 Mb vs. 272 Mb) and contained 40% more genes. Core genes were related to essential cell functions, while non-core genes were related to environmental adaptations.
[102]
Medicago spp.SNP array9277Genetic diversity and population structure analysisPrincipal component analysis of 280 alfalfa (Medicago sativa L.) genotypes revealed grouping according to species and ploidy level. Tetraploids had higher heterozygosity levels.
[83]
Sequence capture50 genesPhylogenetic analysisThe sequences of 50 genes were recovered from samples of four Medicago spp. The alignment of such genes revealed known phylogenetic relationships between the Medicago spp. and an outgroup.
[103]
Northern switchgrass (Panicum virginatum L.)Sequence capture1,590,653 SNPsGenetic diversity and population structure analysis537 individuals from 45 upland and 21 lowland populations were genotyped. Population structure analysis revealed seven groups that match geographical origin and up- or lowland distribution.
[104]
Orchid (Cypripedium macranthos var. rebunense)Amplicon sequencing~1000 SNPsGenetic diversity analysisEight samples from two populations of a Japanese orchid were genotyped with 16 MIG-seq primers. A principal component analysis of 209 SNPs was able to differentiate the two populations. Remarkably, the MIG-seq primers were also functional in a wide spectrum of organisms, including a pine tree, bamboo, mushroom, a copepod, a snail, a sea cucumber, and a lizard.
[105]
Sorghum (Sorghum bicolor [L.] Moench.)DArT12,000 probesGenetic diversity analysis90 accessions of sorghum were genotyped with DArT. Cluster analysis revealed 13 groups of genotypes that match their race and origin.
[91]
Wheat (Triticum aestivum L.)SNP array280,226 SNPsGenetic diversity and population structure analysis4506 wheat landraces and cultivars were genotyped with the SNP array. 11 high-confidence groups were detected based on the genetic structure of the samples.
[85]
DArT5137 probesGenetic diversity analysis and genetic mapping62 wheat cultivars, plus a double haploid population were genotyped with DArT. The DArT markers can capture genetic differentiation among cultivars.
[90]
GBS
DArT
38,412 SNPs
1544 probes
Genetic diversity analysis365 soft winter wheat breeding lines were genotyped with GBS and DArT. The genetic diversity estimates of the two methods differed, with DArT overestimating genetic diversity, likely due to ascertainment bias.
[84]

Share and Cite

MDPI and ACS Style

Loera-Sánchez, M.; Studer, B.; Kölliker, R. DNA-Based Assessment of Genetic Diversity in Grassland Plant Species: Challenges, Approaches, and Applications. Agronomy 2019, 9, 881. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy9120881

AMA Style

Loera-Sánchez M, Studer B, Kölliker R. DNA-Based Assessment of Genetic Diversity in Grassland Plant Species: Challenges, Approaches, and Applications. Agronomy. 2019; 9(12):881. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy9120881

Chicago/Turabian Style

Loera-Sánchez, Miguel, Bruno Studer, and Roland Kölliker. 2019. "DNA-Based Assessment of Genetic Diversity in Grassland Plant Species: Challenges, Approaches, and Applications" Agronomy 9, no. 12: 881. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy9120881

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop