Assessment of genetic diversity in cultivars of white clover (Trifolium repens L.) detected by SSR polymorphisms

Outlier analyses and genome-wide association study identify glgC and ERD6-like 4 as candidate genes for foliar water-soluble carbohydrate accumulation in Trifolium repens

Increasing water-soluble carbohydrate (WSC) content in white clover is important for improving nutritional quality and reducing environmental impacts from pastoral agriculture. Elucidation of genes responsible for foliar WSC variation would enhance genetic improvement by enabling molecular breeding approaches. The aim of the present study was to identify single nucleotide polymorphisms (SNPs) associated with variation in foliar WSC in white clover. A set of 935 white clover individuals, randomly sampled from five breeding pools selectively bred for divergent (low or high) WSC content, were assessed with 14,743 genotyping-by-sequencing SNPs, using three outlier detection methods: PCAdapt, BayeScan and KGD-F_ST. These analyses identified 33 SNPs as discriminating between high and low WSC populations and putatively under selection. One SNP was located in the intron of ERD6-like 4, a gene coding for a sugar transporter located on the vacuole membrane. A genome-wide association study using a subset of 605 white clover individuals and 5,757 SNPs, identified a further 12 SNPs, one of which was associated with a starch biosynthesis gene, glucose-1-phosphate adenylyltransferase, glgC. Our results provide insight into genomic regions underlying WSC accumulation in white clover, identify candidate genomic regions for further functional validation studies, and reveal valuable information for marker-assisted or genomic selection in white clover.

EST-SSR Primer Development and Genetic Structure Analysis of Psathyrostachys juncea Nevski

Psathyrostachys juncea is a perennial forage grass which plays an important role in soil and water conservation and ecological maintenance in cold and dry areas of temperate regions. In P. juncea, a variety of biotic and abiotic stress related genes have been used in crop improvement, indicating its agronomic, economic, forage, and breeding value. To date, there have been few studies on the genetic structure of P. juncea. Here, the genetic diversity and population structure of P. juncea were analyzed by EST-SSR molecular markers to evaluate the genetic differentiation related to tillering traits in P. juncea germplasm resources. The results showed that 400 simple sequence repeat (SSR) loci were detected in 2,020 differentially expressed tillering related genes. A total of 344 scored bands were amplified using 103 primer pairs, out of which 308 (89.53%) were polymorphic. The Nei's gene diversity of 480 individuals was between 0.092 and 0.449, and the genetic similarity coefficient was between 0.5008 and 0.9111, with an average of 0.6618. Analysis of molecular variance analysis showed that 93% of the variance was due to differences within the population, and the remaining 7% was due to differences among populations. Psathyrostachys juncea materials were clustered into five groups based on population genetic structure, principal coordinate analysis and unweighted pair-group method with arithmetic means (UPGMA) analysis. The results were similar between clustering methods, but a few individual plants were distributed differently by the three models. The clustering results, gene diversity and genetic similarity coefficients showed that the overall genetic relationship of P. juncea individuals was relatively close. A Mantel test, UPGMA and structural analysis also showed a significant correlation between genetic relationship and geographical distribution. These results provide references for future breeding programs, genetic improvement and core germplasm collection of P. juncea.

Genetic diversity and population structure analysis in a large collection of Vicia amoena in China with newly developed SSR markers

Vicia amoena is a high-nutritional quality forage similar to alfalfa. However, studies on the genetic background of V. amoena are scarce. In the present study, the genetic variation of 24 V. amoena populations was assessed with newly developed simple sequence repeat (SSR) markers. A total of 8799 SSRs were identified in the V. amoena genomic-enriched sequences, and the most abundant repeat number was four. A total of 569 sampled individuals were assayed to evaluate the genetic diversity of the V. amoena populations based on 21 polymorphic SSR primers. The polymorphism information content (PIC) ranged from 0.896 to 0.968, with an average of 0.931, which indicated that the markers were highly informative. Based on analysis of molecular variance, 88% of the variance occurred within populations, and the remaining 12% of the variance occurred among populations. The high degree of gene flow (Nm= 4.958) also showed slight differentiation among the V. amoena populations. The V. amoena populations were mainly clustered by steppe and mountain habitats based on principal coordinate analysis (PCoA) and STRUCTURE analysis. This indicated that the elevation and special habitat of geographical origins may be important factors affecting the clustered pattern of V. amoena populations. Neighbour-joining (NJ) analysis did not separate the populations well by geographical origin, which indicated that the genetic structure of V. amoena was complex and needs further study. Overall, our results showed that the newly developed SSR markers could benefit the V. amoena research community by providing genetic background information to help establish a foundation for breeding improvement and germplasm resource conservation.

Genetic trade-offs underlie divergent life history strategies for local adaptation in white clover

Local adaptation is common in plants, yet characterization of its underlying genetic basis is rare in herbaceous perennials. Moreover, while many plant species exhibit intraspecific chemical defence polymorphisms, their importance for local adaptation remains poorly understood. We examined the genetic architecture of local adaptation in a perennial, obligately-outcrossing herbaceous legume, white clover (Trifolium repens). This widespread species displays a well-studied chemical defence polymorphism for cyanogenesis (HCN release following tissue damage) and has evolved climate-associated cyanogenesis clines throughout its range. Two biparental F₂  mapping populations, derived from three parents collected in environments spanning the U.S. latitudinal species range (Duluth, MN, St. Louis, MO and Gainesville, FL), were grown in triplicate for two years in reciprocal common garden experiments in the parental environments (6,012 total plants). Vegetative growth and reproductive fitness traits displayed trade-offs across reciprocal environments, indicating local adaptation. Genetic mapping of fitness traits revealed a genetic architecture characterized by allelic trade-offs between environments, with 100% and 80% of fitness QTL in the two mapping populations showing significant QTL×E interactions, consistent with antagonistic pleiotropy. Across the genome there were three hotspots of QTL colocalization. Unexpectedly, we found little evidence that the cyanogenesis polymorphism contributes to local adaptation. Instead, divergent life history strategies in reciprocal environments were major fitness determinants: selection favoured early investment in flowering at the cost of multiyear survival in the southernmost site versus delayed flowering and multiyear persistence in the northern environments. Our findings demonstrate that multilocus genetic trade-offs contribute to contrasting life history characteristics that allow for local adaptation in this outcrossing herbaceous perennial.

Genetic diversity and population structure analysis in a large collection of white clover (Trifolium repens L.) germplasm worldwide

White clover is an important temperate legume forage with high nutrition. In the present study, 448 worldwide accessions were evaluated for the genetic variation and polymorphisms using 22 simple sequence repeat (SSR) markers. All the markers were highly informative, a total of 341 scored bands were amplified, out of which 337 (98.83%) were polymorphic. The PIC values ranged from 0.89 to 0.97 with an average of 0.95. For the AMOVA analysis, 98% of the variance was due to differences within the population and the remaining 2% was due to differences among populations. The white clover accessions were divided into different groups or subgroups based on PCoA, UPGMA, and STRUCTURE analyses. The existence of genetic differentiation between the originally natural and introduced areas according to the PCoA analysis of the global white clover accessions. There was a weak correlation between genetic relationships and geographic distribution according to UPGMA and STRUCTURE analyses. The results of the present study will provide the foundation for future breeding programs, genetic improvement, core germplasm collection establishment for white clover.

Understanding the Complexity of Cold Tolerance in White Clover using Temperature Gradient Locations and a GWAS Approach

White clover ( L.) is the most important grazing perennial forage legume in temperate climates. However, its limited capacity to survive and restore growth after low temperatures during winter constrains the productivity and wide adoption of the crop. Despite the importance of cold tolerance for white clover cultivar development, the genetic basis of this trait remains largely unknown. Hence, in this study, we performed the first genome-wide association study (GWAS) analyses in white clover to identify quantitative trait loci (QTL) for cold-tolerance-related traits. Seeds from 192 divergent genotypes from six populations in the Patagonia region of South America were collected and seed-derived plants were further clonally propagated. Clonal trials were established in three locations representing temperature gradient associated with elevation. Given the allotetraploid nature of the white clover genome, distinct genetic models (diploid and tetraploid) were tested. Only the tetraploid parameterization was able to detect the 53 loci associated with cold-tolerance traits. Out of the 53 single nucleotide polymorphism (SNP) trait associations, 17 controlled more than one trait or were stable across multiple sites. This work represents the first report of QTL for cold-tolerance-related traits, providing insights into its genetic basis and candidate genomic regions for further functional validation studies.

Continent-Wide Climatic Variation Drives Local Adaptation in North American White Clover

Climate-associated clines in adaptive polymorphisms are commonly cited as evidence of local adaptation within species. However, the contribution of the clinally varying trait to overall fitness is often unknown. To address this question, we examined survival, vegetative growth, and reproductive output in a central US common garden experiment using 161 genotypes of white clover (Trifolium repens L.) originating from 15 locations across North America. White clover is polymorphic for cyanogenesis (hydrogen cyanide release upon tissue damage), a chemical defense against generalist herbivores, and climate-associated cyanogenesis clines have repeatedly evolved across the species range. Over a 12-month experiment, we observed striking correlations between the population of origin and plant performance in the common garden, with climatic distance from the common garden site predicting fitness more accurately than geographic distance. Assessments of herbivore leaf damage over the 2015 growing season indicated marginally lower herbivory on cyanogenic plants; however, this effect did not result in increased fitness in the common garden location. Linear mixed modeling suggested that while cyanogenesis variation had little predictive value for vegetative growth, it is as important as climatic variation for predicting reproductive output in the central United States. Together, our findings suggest that knowledge of climate similarity, as well as knowledge of locally favored adaptive traits, will help to inform transplantation strategies for restoration ecology and other conservation efforts in the face of climate change.

Exploring genetic variability within lentil (Lens culinaris Medik.) and across related legumes using a newly developed set of microsatellite markers

Lentil (Lens culinaris Medik.) is an economically important grain legume, yet the genetic and genomic resources remain largely uncharacterized and unexploited in this crop. Microsatellites have become markers of choice for crop improvement applications. Hence, simple sequence repeat (SSR) markers were developed for lentil through the construction of genomic library enriched for GA/CT motifs. As a result 122 functional SSR primer pairs were developed from 151 microsatellite loci and validated in L. culinaris cv. Precoz. Thirty three SSR markers were utilized for the analysis of genetic relationships between cultivated and wild species of Lens and related legumes. A total of 123 alleles were amplified at 33 loci ranging from 2-5 alleles with an average of 3.73 alleles per locus. Polymorphic information content (PIC) for all the loci ranged from 0.13 to 0.99 with an average of 0.66 per locus. Varied levels of cross genera transferability were obtained ranging from 69.70 % across Pisum sativum to 12.12 % across Vigna radiata. The UPGMA based dendrogram was able to establish the uniqueness of each genotype and grouped them into two major clusters clearly resolving the genetic relationships within lentil and related species. The new set of SSR markers reported here were efficient and highly polymorphic and would add to the existing repertoire of lentil SSR markers to be utilized in molecular breeding. Moreover, the improved knowledge about intra- and inter-specific genetic relationships would facilitate germplasm utilization for lentil improvement.

Aridity shapes cyanogenesis cline evolution in white clover (Trifolium repens L.)

Adaptive differentiation between populations is often proposed to be the product of multiple interacting selective pressures, although empirical support for this is scarce. In white clover, populations show adaptive differentiation in frequencies of cyanogenesis, the ability to produce hydrogen cyanide after tissue damage. This polymorphism arises through independently segregating polymorphisms for the presence/absence of two required cyanogenic components, cyanogenic glucosides and their hydrolysing enzyme. White clover populations worldwide have evolved a series of recurrent, climate-associated clines, with higher frequencies of cyanogenic plants in warmer locations. These clines have traditionally been hypothesized to reflect a fitness trade-off between chemical defence in herbivore-rich areas (warmer climates) and energetic costs of producing cyanogenic components in areas of low herbivore pressure (cooler climates). Recent observational studies suggest that cyanogenic components may also be beneficial in water-stressed environments. We investigated fitness trade-offs associated with temperature-induced water stress in the cyanogenesis system using manipulative experiments in growth chambers and population surveys across a longitudinal precipitation gradient in the central United States. We find that plants producing cyanogenic glucosides have higher relative fitness in treatments simulating a moderate, persistent drought stress. In water-neutral treatments, there are energetic costs to producing cyanogenic components, but only in treatments with nutrient stress. These fitness trade-offs are consistent with cyanogenesis frequencies in natural populations, where we find clinal variation in the proportion of plants producing cyanogenic glucosides along the precipitation gradient. These results suggest that multiple selective pressures interact to maintain this adaptive polymorphism and that modelling adaptation will require knowledge of environment-specific fitness effects.

Recurrent gene deletions and the evolution of adaptive cyanogenesis polymorphisms in white clover (Trifolium repens L.)

Understanding the molecular evolution of genes that underlie intraspecific polymorphisms can provide insights into the process of adaptive evolution. For adaptive polymorphisms characterized by gene presence/absence (P/A) variation, underlying loci commonly show signatures of long-term balancing selection, with gene-presence and gene-absence alleles maintained as two divergent lineages. We examined the molecular evolution of two unlinked P/A polymorphisms that underlie a well-documented adaptive polymorphism for cyanogenesis (hydrogen cyanide release with tissue damage) in white clover. Both cyanogenic and acyanogenic plants occur in this species, and the ecological forces that maintain this chemical defence polymorphism have been studied for several decades. Using a sample of 65 plants, we investigated the molecular evolution of sequences flanking the two underlying cyanogenesis genes: Ac/ac (controlling the presence/absence of cyanogenic glucosides) and Li/li (controlling the presence/absence of their hydrolysing enzyme, linamarase). A combination of genome walking, PCR assays, DNA sequence analysis and Southern blotting was used to test whether these adaptive P/A polymorphisms show evidence of long-term balancing selection, or whether gene-absence alleles have evolved repeatedly through independent deletion events. For both loci, we detect no signatures of balancing selection in the closest flanking genomic sequences. Instead, we find evidence for variation in the size of the deletions characterizing gene-absence alleles. These observations strongly suggest that both of these polymorphisms have been evolving through recurrent gene deletions over time. We discuss the genetic mechanisms that could account for this surprising pattern and the implications of these findings for mechanisms of rapid adaptive evolution in white clover.

Molecular characterisation and interpretation of genetic diversity within globally distributed germplasm collections of tall fescue (Festuca arundinacea Schreb.) and meadow fescue (F. pratensis Huds.)

Allohexaploid tall fescue (Festuca arundinacea Schreb. syn. Lolium arundinaceum [Schreb.] Darbysh.) is an agriculturally important grass cultivated for pasture and turf world-wide. Genetic improvement of tall fescue could benefit from the use of non-domesticated germplasm to diversify breeding populations through the incorporation of novel and superior allele content. However, such potential germplasm must first be characterised, as three major morphotypes (Continental, Mediterranean and rhizomatous) with varying degrees of hybrid interfertility are commonly described within this species. As hexaploid tall fescue is also a member of a polyploid species complex that contains tetraploid, octoploid and decaploid taxa, it is also possible that germplasm collections may have inadvertently sampled some of these sub-species. In this study, 1,040 accessions from the publicly available United States Department of Agriculture tall fescue and meadow fescue germplasm collections were investigated. Sequence of the chloroplast genome-located matK gene and the nuclear ribosomal DNA internal transcribed spacer (rDNA ITS) permitted attribution of accessions to the three previously known morphotypes and also revealed the presence of tall fescue sub-species of varying ploidy levels, as well as other closely related species. The majority of accessions were, however, identified as Continental hexaploid tall fescue. Analysis using 34 simple sequence repeat markers was able to further investigate the level of genetic diversity within each hexaploid tall fescue morphotype group. At least two genetically distinct sub-groups of Continental hexaploid tall fescue were identified which are probably associated with palaeogeographic range expansion of this morphotype. This work has comprehensively characterised a large and complex germplasm collection and has identified genetically diverse accessions which may potentially contribute valuable alleles at agronomic loci for tall fescue cultivar improvement programs.

Molecular breeding of transgenic white clover (Trifolium repens L.) with field resistance to Alfalfa mosaic virus through the expression of its coat protein gene

Viral diseases, such as Alfalfa mosaic virus (AMV), cause significant reductions in the productivity and vegetative persistence of white clover plants in the field. Transgenic white clover plants ectopically expressing the viral coat protein gene encoded by the sub-genomic RNA4 of AMV were generated. Lines carrying a single copy of the transgene were analysed at the molecular, biochemical and phenotypic level under glasshouse and field conditions. Field resistance to AMV infection, as well as mitotic and meiotic stability of the transgene, were confirmed by phenotypic evaluation of the transgenic plants at two sites within Australia. The T(0) and T(1) generations of transgenic plants showed immunity to infection by AMV under glasshouse and field conditions, while the T(4) generation in an agronomically elite 'Grasslands Sustain' genetic background, showed a very high level of resistance to AMV in the field. An extensive biochemical study of the T(4) generation of transgenic plants, aiming to evaluate the level and composition of natural toxicants and key nutritional parameters, showed that the composition of the transgenic plants was within the range of variation seen in non-transgenic populations.

Genetic diversity of Aquilegia (Ranunculaceae) species and cultivars assessed by AFLPs

Species of the genus Aquilegia are exceptionally diverse in their floral morphology and color, commonly known as columbine. They are widely planted ornamentals and are highly attractive for hummingbirds. However, little is known about their genetic diversity. We examined the genetic diversity of the species and cultivars using amplified fragment length polymorphism (AFLP) markers. Sixteen EcoRI/MseI AFLP primer combinations produced 327 informative polymorphic bands, with a mean of 20.4 bands scored per primer. Jaccard's coefficient of similarity varied from 0.61 to 0.93, indicative of high levels of genetic variation. Cluster analysis using the unweighted pair group method with arithmetic mean algorithm placed the 64 accessions into two main clusters, each divided into two sub-clusters. The AFLP variability was significantly associated with the geographic origins, as the Asian species and the North American species grouped into two distinct clusters. The genetic diversity found among Aquilegia demonstrated the potential value of Chinese germplasm for cultivar improvement and for widening the genetic basis of breeding programs and breeding material selection. We concluded that AFLPs are informative and can provide significant insights for genetic diversity research in columbine species.

High levels of nucleotide diversity and fast decline of linkage disequilibrium in rye (Secale cereale L.) genes involved in frost response

BACKGROUND

Rye (Secale cereale L.) is the most frost tolerant cereal species. As an outcrossing species, rye exhibits high levels of intraspecific diversity, which makes it well-suited for allele mining in genes involved in the frost responsive network. For investigating genetic diversity and the extent of linkage disequilibrium (LD) we analyzed eleven candidate genes and 37 microsatellite markers in 201 lines from five Eastern and Middle European rye populations.

RESULTS

A total of 147 single nucleotide polymorphisms (SNPs) and nine insertion-deletion polymorphisms were found within 7,639 bp of DNA sequence from eleven candidate genes, resulting in an average SNP frequency of 1 SNP/52 bp. Nucleotide and haplotype diversity of candidate genes were high with average values π = 5.6 × 10(-3) and Hd = 0.59, respectively. According to an analysis of molecular variance (AMOVA), most of the genetic variation was found between individuals within populations. Haplotype frequencies varied markedly between the candidate genes. ScCbf14, ScVrn1, and ScDhn1 were dominated by a single haplotype, while the other 8 genes (ScCbf2, ScCbf6, ScCbf9b, ScCbf11, ScCbf12, ScCbf15, ScIce2, and ScDhn3) had a more balanced haplotype frequency distribution. Intra-genic LD decayed rapidly, within approximately 520 bp on average. Genome-wide LD based on microsatellites was low.

CONCLUSIONS

The Middle European population did not differ substantially from the four Eastern European populations in terms of haplotype frequencies or in the level of nucleotide diversity. The low LD in rye compared to self-pollinating species promises a high resolution in genome-wide association mapping. SNPs discovered in the promoters or coding regions, which attribute to non-synonymous substitutions, are suitable candidates for association mapping.

Comparison of genome structure between white clover and Medicago truncatula supports homoeologous group nomenclature based on conserved synteny

Computational analysis has been used to align the genetic map of white clover (Trifolium repens L.) with the draft genome sequence of the model legume species Medicago truncatula Gaertn. In silico comparison based on white clover expressed sequence tags that contain simple sequence repeat loci revealed substantial macrosynteny between the genomes of these two species, which are closely related within the Trifolieae tribe of the Fabaceae family. Six of the eight homoeologous chromosome groups (HGs) of allotetraploid white clover show predominant relationships with single M. truncatula (Mt) chromosomes, while the two remaining groups may have participated in an evolutionary reciprocal translocation event. On this basis, a new chromosome nomenclature system for allotetraploid white clover is proposed such that HG A = 3, HG B = 8, HG C = 7, HG D = 4, HG E = 1, HG F = 2, HG G = 5, and HG H = 6. A rationalized linkage map ordering system has also been demonstrated. Improved knowledge of the relationships between agricultural and model forage legume genomes will facilitate prediction of gene location for key agronomic traits for pasture production.

Identification of homologous, homoeologous and paralogous sequence variants in an outbreeding allopolyploid species based on comparison with progenitor taxa

The combination of homologous, homoeologous and paralogous classes of sequence variation presents major challenges for SNP discovery in outbreeding allopolyploid species. Previous in vitro gene-associated SNP discovery studies in the allotetraploid forage legume white clover (Trifolium repens L.) were vulnerable to such effects, leading to prohibitive levels of attrition during SNP validation. Identification of T. occidentale and T. pallescens as the putative diploid progenitors of white clover has permitted discrimination of the different sequence variant categories. Amplicons from selected abiotic stress tolerance-related genes were obtained using mapping family parents and individuals from each diploid species. Following cloning, progenitor comparison allowed tentative assignment of individual haplotypes to one or other sub-genome, as well as to gene copies within sub-genomes. A high degree of coincidence and identity between SNPs and HSVs was observed. Close similarity was observed between the genome of T. occidentale and one white clover sub-genome, but the affinity between T. pallescens and the other sub-genome was weaker, suggesting that a currently uncharacterised taxon may be the true second progenitor. Selected validated SNPs were attributed to individual sub-genomes by assignment to and naming of homoeologous linkage groups, providing the basis for improved genetic trait-dissection studies. The approach described in this study is broadly applicable to a range of allopolyploid taxa of equivocal ancestry.

Validation of in silico-predicted genic SNPs in white clover (Trifolium repens L.), an outbreeding allopolyploid species

White clover (Trifolium repens L.) is an obligate outbreeding allotetraploid forage legume. Gene-associated SNPs provide the optimum genetic system for improvement of such crop species. An EST resource obtained from multiple cDNA libraries constructed from numerous genotypes of a single cultivar has been used for in silico SNP discovery and validation. A total of 58 from 236 selected sequence clusters (24.5%) were fully validated as containing polymorphic SNPs by genotypic analysis across the parents and progeny of several two-way pseudo-testcross mapping families. The clusters include genes belonging to a broad range of predicted functional categories. Polymorphic SNP-containing ESTs have also been used for comparative genomic analysis by comparison with whole genome data from model legume species, as well as Arabidopsis thaliana. A total of 29 (50%) of the 58 clusters detected putative ortholoci with known chromosomal locations in Medicago truncatula, which is closely related to white clover within the Trifolieae tribe of the Fabaceae. This analysis provides access to translational data from model species. The efficiency of in silico SNP discovery in white clover is limited by paralogous and homoeologous gene duplication effects, which are resolved unambiguously by the transmission test. This approach will also be applicable to other agronomically important cross-pollinating allopolyploid plant species.