An integrated genetic linkage map for white clover (Trifolium repens L.) with alignment to Medicago

An improved reference genome for Trifolium subterraneum L. provides insight into molecular diversity and intra-specific phylogeny

Subterranean clover (Trifolium subterraneum L., Ts) is a geocarpic, self-fertile annual forage legume with a compact diploid genome (n = x = 8, 544 Mb/1C). Its resilience and climate adaptivity have made it an economically important species in Mediterranean and temperate zones. Using the cultivar Daliak, we generated higher resolution sequence data, created a new genome assembly TSUd_3.0, and conducted molecular diversity analysis for copy number variant (CNV) and single-nucleotide polymorphism (SNP) among 36 cultivars. TSUd_3.0 substantively improves prior genome assemblies with new Hi-C and long-read sequence data, covering 531 Mb, containing 41,979 annotated genes and generating a 94.4% BUSCO score. Comparative genomic analysis among select members of the tribe Trifolieae indicated TSUd 3.0 corrects six assembly-error inversion/duplications and confirmed phylogenetic relationships. Its synteny with T. pratense, T. repens, Medicago truncatula and Lotus japonicus genomes were assessed, with the more distantly related T. repens and M. truncatula showing higher levels of co-linearity with Ts than between Ts and its close relative T. pratense. Resequencing of 36 cultivars discovered 7,789,537 SNPs subsequently used for genomic diversity assessment and sequence-based clustering. Heterozygosity estimates ranged from 1% to 21% within the 36 cultivars and may be influenced by admixture. Phylogenetic analysis supported subspecific genetic structure, although it indicates four or five groups, rather than the three recognized subspecies. Furthermore, there were incidences where cultivars characterized as belonging to a particular subspecies clustered with another subspecies when using genomic data. These outcomes suggest that further investigation of Ts sub-specific classification using molecular and morpho-physiological data is needed to clarify these relationships. This upgraded reference genome, complemented with comprehensive sequence diversity analysis of 36 cultivars, provides a platform for future gene functional analysis of key traits, and genome-based breeding strategies for climate adaptation and agronomic performance. Pangenome analysis, more in-depth intra-specific phylogenomic analysis using the Ts core collection, and functional genetic and genomic studies are needed to further augment knowledge of Trifolium genomes.

R2R3-MYB gene family: Genome-wide identification provides insight to improve the content of proanthocyanidins in Trifolium repens

The R2R3-MYB family is one of largest transcription factor families in plants playing significant roles in regulating anthocyanin and proanthocyanidin biosynthesis. Proanthocyanidins are one of major objectives to improve the quality of white clover (Trifolium repens L.), which have a beneficial effect on ruminant to prevent the lethal pasture bloat. A total of 133 TrR2R3-MYB genes were identified and distributed on all 16 chromosomes based on the whole genome information of white clover. Also, by exploring the gene structure, motifs and duplication events of TrR2R3-MYBs, as well as the evolutionary relationship with TrR2R3-MYB genes of other species, 10 TrR2R3-MYB genes with the potential to regulate the anthocyanins and proanthocyanidins biosynthesis were screened. These TrR2R3-MYB genes responded significantly to low temperature in white clover. In addition, they have different expression patterns in leaves, petioles and inflorescences of white clover. Importantly, TrMYB116 and TrMYB118 may positively regulate anthocyanin accumulation and low temperature response in white clover. TrMYB118 may also be associated with anthocyanin pigmentation pattern in Purple leaves. This study provides a basis for verifying the function of TrR2R3-MYB and breeding white clover cultivars with high proanthocyanidins.

Sugar Transporters in Plasmodiophora brassicae: Genome-Wide Identification and Functional Verification

Plasmodiophora brassicae, an obligate intracellular pathogen, can hijack the host’s carbohydrates for survival. When the host plant is infected by P. brassicae, a large amount of soluble sugar accumulates in the roots, especially glucose, which probably facilitates the development of this pathogen. Although a complete glycolytic and tricarboxylic acid cycle (TCA) cycle existed in P. brassicae, very little information about the hexose transport system has been reported. In this study, we screened 17 putative sugar transporters based on information about their typical domains. The structure of these transporters showed a lot of variation compared with that of other organisms, especially the number of transmembrane helices (TMHs). Phylogenetic analysis indicated that these sugar transporters were far from the evolutionary relationship of other organisms and were unique in P. brassicae. The hexose transport activity assay indicated that eight transporters transported glucose or fructose and could restore the growth of yeast strain EBY.VW4000, which was deficient in hexose transport. The expression level of these glucose transporters was significantly upregulated at the late inoculation time when resting spores and galls were developing and a large amount of energy was needed. Our study provides new insights into the mechanism of P. brassicae survival in host cells by hijacking and utilizing the carbohydrates of the host.

Chromosome and Genome Diversity in the Genus Trifolium (Fabaceae)

Trifolium L. is an economically important genus that is characterized by variable karyotypes relating to its ploidy level and basic chromosome numbers. The advent of genomic resources combined with molecular cytogenetics provides an opportunity to develop our understanding of plant genomes in general. Here, we summarize the current state of knowledge on Trifolium genomes and chromosomes and review methodologies using molecular markers that have contributed to Trifolium research. We discuss possible future applications of cytogenetic methods in research on the Trifolium genome and chromosomes.

Dual-species origin of an adaptive chemical defense polymorphism

Allopolyploid speciation and chemical defense diversification are two of the most characteristic features of plant evolution; although the former has likely shaped the latter, this has rarely been documented. Here we document allopolyploidy-mediated chemical defense evolution in the origin of cyanogenesis (HCN release upon tissue damage) in white clover (Trifolium repens). We combined linkage mapping of the loci that control cyanogenesis (Ac, controlling production of cyanogenic glucosides; and Li, controlling production of their hydrolyzing enzyme linamarase) with genome sequence comparisons between white clover, a recently evolved allotetraploid, and its diploid progenitors (Trifolium pallescens, Trifolium occidentale). The Ac locus (a three-gene cluster comprising the cyanogenic glucoside pathway) is derived from T. occidentale; it maps to linkage group 2O (occidentale subgenome) and is orthologous to a highly similar cluster in the T. occidentale reference genome. By contrast, Li maps to linkage group 4P (pallescens subgenome), indicating an origin in the other progenitor species. These results indicate that cyanogenesis evolved in white clover as a product of the interspecific hybridization that created the species. This allopolyploidization-derived chemical defense, together with subsequent selection on intraspecific cyanogenesis variation, appears to have contributed to white clover's ecological success as a globally distributed weed species.

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.

Fingerprint identification of white clover cultivars based on SSR molecular markers

White clover (Trifolium repens L.) is an important perennial legume forage with high productivity and quality. To strengthen the basic research on the genetic characteristics, fingerprint identification and adaptability of white clover germplasm resources, Simple sequence repeat (SSR) molecular markers were applied to 10 white clover cultivars to assess the genetic diversity and related lines of white clover at the molecular level in order to lay a theoretical foundation for the selection of high-quality seeds and cultivars of white clover. A total of 120 different bands were amplified by 29 pairs of SSR primers with good polymorphism, of which 103 (89.5%) were polymorphic. Meanwhile, the PIC of each primer was 0.181-0.588, with an average of 0.329. Analysis of molecular variance revealed that 57% of the genetic variation occurred within cultivars and 43% occurred among cultivars. The results of cluster analysis and the principal coordinate analysis revealed that the parental relationships of the 10 cultivars, with the 'Purple' cultivar very distantly related to the other 9 cultivars and the closest parental relationship between 'Ladino' and 'Sulky'. The fingerprints constructed by three representative primers (gtrs679, gtrs319, and gtrs678) have a strong identification ability. In summary, the SSR markers had good polymorphism and could be used for DNA fingerprint analysis of white clover cultivars.

Breaking Free: The Genomics of Allopolyploidy-Facilitated Niche Expansion in White Clover

The merging of distinct genomes, allopolyploidization, is a widespread phenomenon in plants. It generates adaptive potential through increased genetic diversity, but examples demonstrating its exploitation remain scarce. White clover (Trifolium repens) is a ubiquitous temperate allotetraploid forage crop derived from two European diploid progenitors confined to extreme coastal or alpine habitats. We sequenced and assembled the genomes and transcriptomes of this species complex to gain insight into the genesis of white clover and the consequences of allopolyploidization. Based on these data, we estimate that white clover originated ∼15,000 to 28,000 years ago during the last glaciation when alpine and coastal progenitors were likely colocated in glacial refugia. We found evidence of progenitor diversity carryover through multiple hybridization events and show that the progenitor subgenomes have retained integrity and gene expression activity as they traveled within white clover from their original confined habitats to a global presence. At the transcriptional level, we observed remarkably stable subgenome expression ratios across tissues. Among the few genes that show tissue-specific switching between homeologous gene copies, we found flavonoid biosynthesis genes strongly overrepresented, suggesting an adaptive role of some allopolyploidy-associated transcriptional changes. Our results highlight white clover as an example of allopolyploidy-facilitated niche expansion, where two progenitor genomes, adapted and confined to disparate and highly specialized habitats, expanded to a ubiquitous global presence after glaciation-associated allopolyploidization.

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.

Accounting for Errors in Low Coverage High-Throughput Sequencing Data When Constructing Genetic Maps Using Biparental Outcrossed Populations

Next-generation sequencing is an efficient method that allows for substantially more markers than previous technologies, providing opportunities for building high-density genetic linkage maps, which facilitate the development of nonmodel species' genomic assemblies and the investigation of their genes. However, constructing genetic maps using data generated via high-throughput sequencing technology (e.g., genotyping-by-sequencing) is complicated by the presence of sequencing errors and genotyping errors resulting from missing parental alleles due to low sequencing depth. If unaccounted for, these errors lead to inflated genetic maps. In addition, map construction in many species is performed using full-sibling family populations derived from the outcrossing of two individuals, where unknown parental phase and varying segregation types further complicate construction. We present a new methodology for modeling low coverage sequencing data in the construction of genetic linkage maps using full-sibling populations of diploid species, implemented in a package called GUSMap. Our model is based on the Lander-Green hidden Markov model but extended to account for errors present in sequencing data. We were able to obtain accurate estimates of the recombination fractions and overall map distance using GUSMap, while most existing mapping packages produced inflated genetic maps in the presence of errors. Our results demonstrate the feasibility of using low coverage sequencing data to produce genetic maps without requiring extensive filtering of potentially erroneous genotypes, provided that the associated errors are correctly accounted for in the model.

Anthocyanin leaf markings are regulated by a family of R2R3-MYB genes in the genus Trifolium

Anthocyanin pigments accumulate to form spatially restricted patterns in plants, particularly in flowers, but also occur in vegetative tissues. Spatially restricted anthocyanin leaf markings are poorly characterised in plants, but are common in forage legumes. We hypothesised that the molecular basis for anthocyanin leaf markings in Trifolium spp. is due to the activity of a family of R2R3-MYB genes. R2R3-MYB genes were identified that are associated with the two classic pigmentation loci in T. repens. The R locus patterns 'red leaf', 'red midrib' and 'red fleck' are conditioned by a single MYB gene, RED LEAF. The 'diffuse red leaf' trait is regulated by the RED LEAF DIFFUSE MYB gene. The V locus was identified through mapping two V-linked traits, 'V-broken yellow' (Vby) and 'red leaflet' (Vrl). Two highly similar R2R3-MYB genes, RED V-a and RED V-b, mapped to the V locus and co-segregated with the RED V pigmentation pattern. Functional characterisation of RED LEAF and RED V was performed, confirming their function as anthocyanin regulators and identifying a C-terminal region necessary for transactivation. The mechanisms responsible for generating anthocyanin leaf markings in T. repens provide a valuable system to compare with mechanisms that regulate complex floral pigmentation.

Breaking through the feed barrier: options for improving forage genetics

Pasture based on perennial ryegrass (Lolium perenne L.) and white clover (Trifolium repens L.) is the foundation for production and profit in the Australasian pastoral sectors. The improvement of these species offers direct opportunities to enhance sector performance, provided there is good alignment with industry priorities as quantified by means such as the forage value index. However, the rate of forage genetic improvement must increase to sustain industry competitiveness. New forage technologies and breeding strategies that can complement and enhance traditional approaches are required to achieve this. We highlight current and future research in plant breeding, including genomic and gene technology approaches to improve rate of genetic gain. Genomic diversity is the basis of breeding and improvement. Recent advances in the range and focus of introgression from wild Trifolium species have created additional specific options to improve production and resource-use-efficiency traits. Symbiont genetic resources, especially advances in grass fungal endophytes, make a critical contribution to forage, supporting pastoral productivity, with benefits to both pastures and animals in some dairy regions. Genomic selection, now widely used in animal breeding, offers an opportunity to lift the rate of genetic gain in forages as well. Accuracy and relevance of trait data are paramount, it is essential that genomic breeding approaches be linked with robust field evaluation strategies including advanced phenotyping technologies. This requires excellent data management and integration with decision-support systems to deliver improved effectiveness from forage breeding. Novel traits being developed through genetic modification include increased energy content and potential increased biomass in ryegrass, and expression of condensed tannins in forage legumes. These examples from the wider set of research emphasise forage adaptation, yield and energy content, while covering the spectrum from exotic germplasm and symbionts through to advanced breeding strategies and gene technologies. To ensure that these opportunities are realised on farm, continuity of industry-relevant delivery of forage-improvement research is essential, as is sustained research input from the supporting pasture and plant sciences.

A saturated genetic linkage map of autotetraploid alfalfa (Medicago sativa L.) developed using genotyping-by-sequencing is highly syntenous with the Medicago truncatula genome

A genetic linkage map is a valuable tool for quantitative trait locus mapping, map-based gene cloning, comparative mapping, and whole-genome assembly. Alfalfa, one of the most important forage crops in the world, is autotetraploid, allogamous, and highly heterozygous, characteristics that have impeded the construction of a high-density linkage map using traditional genetic marker systems. Using genotyping-by-sequencing (GBS), we constructed low-cost, reasonably high-density linkage maps for both maternal and paternal parental genomes of an autotetraploid alfalfa F₁ population. The resulting maps contain 3591 single-nucleotide polymorphism markers on 64 linkage groups across both parents, with an average density of one marker per 1.5 and 1.0 cM for the maternal and paternal haplotype maps, respectively. Chromosome assignments were made based on homology of markers to the M. truncatula genome. Four linkage groups representing the four haplotypes of each alfalfa chromosome were assigned to each of the eight Medicago chromosomes in both the maternal and paternal parents. The alfalfa linkage groups were highly syntenous with M. truncatula, and clearly identified the known translocation between Chromosomes 4 and 8. In addition, a small inversion on Chromosome 1 was identified between M. truncatula and M. sativa. GBS enabled us to develop a saturated linkage map for alfalfa that greatly improved genome coverage relative to previous maps and that will facilitate investigation of genome structure. GBS could be used in breeding populations to accelerate molecular breeding in alfalfa.

Fast magnetic isolation of simple sequence repeat markers in microfluidic channels

Simple sequence repeat (SSR) markers are widely used for genome mapping, genetic diversity characterization and medical diagnosis. The fast isolation by AFLP of sequence containing repeats (FIASCO) is a powerful method for SSR marker isolation, but it is laborious, costly, and time consuming and requires multiple rounds of washing. Here, we report a superparamagnetic bead (SPMB)-based FIASCO method in a magnetic field controllable microfluidic chip (MFCM-Chip). This method dramatically reduces the assay time by 4.25-fold and reduces the quantity of magnetic beads and probes by 10-fold through the magnetic capture of (AG)n-containing fragments from Herba Leonuri, followed by washing and eluting on a microchip. The feasibility of this method was further evaluated by PCR and sequencing, and the results showed that the proportion of fragments containing SSRs was 89%, confirming that this platform is a fast and efficient method for SSR marker isolation. This cost-effective platform will make the powerful FIASCO technique more accessible for routine use with a wide variety of materials.