A Trifolium repens flavodoxin-like quinone reductase 1 (TrFQR1) improves plant adaptability to high temperature associated with oxidative homeostasis and lipids remodeling.

Maintenance of stable mitochondrial respiratory chains could enhance adaptability to high temperature, but the potential mechanism was not elucidated clearly in plants. In this study, we identified and isolated a TrFQR1 gene encoding the flavodoxin-like quinone reductase 1 (TrFQR1) located in mitochondria of leguminous white clover (Trifolium repens). Phylogenetic analysis indicated that amino acid sequences of FQR1 in various plant species showed a high degree of similarities. Ectopic expression of TrFQR1 protected yeast (Saccharomyces cerevisiae) from heat damage and toxic levels of benzoquinone, phenanthraquinone and hydroquinone. Transgenic Arabidopsis thaliana and white clover overexpressing TrFQR1 exhibited significantly lower oxidative damage and better photosynthetic capacity and growth than wild-type in response to high-temperature stress, whereas AtFQR1-RNAi A. thaliana showed more severe oxidative damage and growth retardation under heat stress. TrFQR1-transgenic white clover also maintained better respiratory electron transport chain than wild-type plants, as manifested by significantly higher mitochondrial complex II and III activities, alternative oxidase activity, NAD(P)H content, and coenzyme Q10 content in response to heat stress. In addition, overexpression of TrFQR1 enhanced the accumulation of lipids including phosphatidylglycerol, monogalactosyl diacylglycerol, sulfoquinovosyl diacylglycerol and cardiolipin as important compositions of bilayers involved in dynamic membrane assembly in mitochondria or chloroplasts positively associated with heat tolerance. TrFQR1-transgenic white clover also exhibited higher lipids saturation level and phosphatidylcholine:phosphatidylethanolamine ratio, which could be beneficial to membrane stability and integrity during a prolonged period of heat stress. The current study proves that TrFQR1 is essential for heat tolerance associated with mitochondrial respiratory chain, cellular reactive oxygen species homeostasis, and lipids remodeling in plants. TrFQR1 could be selected as a key candidate marker gene to screen heat-tolerant genotypes or develop heat-tolerant crops via molecular-based breeding.

Genome-wide identification and characterization of NAC transcription factor family members in Trifolium pratense and expression analysis under lead stress.

BACKGROUND: The NAC TF family is widely involved in plant responses to various types of stress. Red clover (Trifolium pratense) is a high-quality legume, and the study of NAC genes in red clover has not been comprehensive. The aim of this study was to analyze the NAC gene family of red clover at the whole-genome level and explore its potential role in the Pb stress response. RESULTS: In this study, 72 TpNAC genes were identified from red clover; collinearity analysis showed that there were 5 pairs of large fragment replicators of TpNAC genes, and red clover was found to be closely related to Medicago truncatula. Interestingly, the TpNAC genes have more homologs in Arabidopsis thaliana than in soybean (Glycine max). There are many elements in the TpNAC genes promoters that respond to stress. Gene expression analysis showed that all the TpNAC genes responded to Pb stress. qRT-PCR showed that the expression levels of TpNAC29 and TpNAC42 were significantly decreased after Pb stress. Protein interaction network analysis showed that 21 TpNACs and 23 other genes participated in the interaction. In addition, the TpNAC proteins had three possible 3D structures, and the secondary structure of these proteins were mainly of other types. These results indicated that most TpNAC members were involved in the regulation of Pb stress in red clover. CONCLUSION: These results suggest that most TpNAC members are involved in the regulation of Pb stress in red clover. TpNAC members play an important role in the response of red clover to Pb stress.

Genome-wide identification of Aux/IAA gene family in white clover (Trifolium repens L.) and functional verification of TrIAA18 under different abiotic stress.

BACKGROUND: White clover (Trifolium repens L.) is an excellent leguminous cool-season forage with a high protein content and strong nitrogen-fixing ability. Despite these advantages, its growth and development are markedly sensitive to environmental factors. Indole-3-acetic acid (IAA) is the major growth hormone in plants, regulating plant growth, development, and response to adversity. Nevertheless, the specific regulatory functions of Aux/IAA genes in response to abiotic stresses in white clover remain largely unexplored. RESULTS: In this study, we identified 47 Aux/IAA genes in the white clover genome, which were categorized into five groups based on phylogenetic analysis. The TrIAAs promoter region co-existed with different cis-regulatory elements involved in developmental and hormonal regulation, and stress responses, which may be closely related to their diverse regulatory roles. Collinearity analysis showed that the amplification of the TrIAA gene family was mainly carried out by segmental duplication. White clover Aux/IAA genes showed different expression patterns in different tissues and under different stress treatments. In addition, we performed a yeast two-hybrid analysis to investigate the interaction between white clover Aux/IAA and ARF proteins. Heterologous expression indicated that TrIAA18 could enhance stress tolerance in both yeast and transgenic Arabidopsis thaliana. CONCLUSION: These findings provide new scientific insights into the molecular mechanisms of growth hormone signaling in white clover and its functional characteristics in response to environmental stress.

Starmerella kisarazuensis f.a., sp. nov., a novel yeast isolated from Trifolium pratense flowers.

Three yeast isolates, NBRC 115909T, NBRC 115910 and NBRC 116270, were isolated from Trifolium pratense (red clover) flowers collected from Kisarazu, Chiba, Japan. Analysis of the sequences of the D1/D2 domains of the large subunit (LSU) rRNA gene and the internal transcribed spacer (ITS) regions revealed that these isolates represent a single novel species within the genus Starmerella. Also, no ascospore formation was observed. The yeast isolates were closely related to Starmerella vitae UWOPS 00-107.2T and Starmerella bombi NRRL Y-17081T. They differed from S. vitae, the most closely related species with a validly published name, by ten nucleotide substitutions with two gaps in the D1/D2 domains and 20 nucleotide substitutions in the ITS region. Moreover, the three isolates exhibited distinct phenotypic characteristics from the closely related species. Therefore, we suggest that these three isolates represent a novel species, designated as Starmerella kisarazuensis f.a., sp. nov. The holotype is NBRC 115909T (isotype: CBS 18485T).

High-quality chromosome-level de novo assembly of the Trifolium repens.

BACKGROUND: White clover (Trifolium repens L.), an excellent perennial legume forage, is an allotetraploid native to southeastern Europe and southern Asia. It has high nutritional, ecological, genetic breeding, and medicinal values and exhibits excellent resistance to cold, drought, trample, and weed infestation. Thus, white clover is widely planted in Europe, America, and China; however, the lack of reference genome limits its breeding and cultivation. This study generated a white clover de novo genome assembly at the chromosomal level and annotated its components. RESULTS: The PacBio third-generation Hi-Fi assembly and sequencing methods generated a 1096 Mb genome size of T. repens, with contigs of N50 = 14 Mb, scaffolds of N50 = 65 Mb, and BUSCO value of 98.5%. The newly assembled genome has better continuity and integrity than the previously reported white clover reference genome; thus provides important resources for the molecular breeding and evolution of white clover and other forage. Additionally, we annotated 90,128 high-confidence gene models from the genome. White clover was closely related to Trifolium pratense and Trifolium medium but distantly related to Glycine max, Vigna radiata, Medicago truncatula, and Cicer arietinum. The expansion, contraction, and GO functional enrichment analysis of the gene families showed that T. repens gene families were associated with biological processes, molecular function, cellular components, and environmental resistance, which explained its excellent agronomic traits. CONCLUSIONS: This study reports a high-quality de novo assembly of white clover genome obtained at the chromosomal level using PacBio Hi-Fi sequencing, a third-generation sequencing. The generated high-quality genome assembly of white clover provides a key basis for accelerating the research and molecular breeding of this important forage crop. The genome is also valuable for future studies on legume forage biology, evolution, and genome-wide mapping of quantitative trait loci associated with the relevant agronomic traits.

Rhizobium nodule diversity and composition are influenced by clover host selection and local growth conditions.

While shaping of plant microbiome composition through 'host filtering' is well documented in legume-rhizobium symbioses, it is less clear to what extent different varieties and genotypes of the same plant species differentially influence symbiont community diversity and composition. Here, we compared how clover host varieties and genotypes affect the structure of Rhizobium populations in root nodules under conventional field and controlled greenhouse conditions. We first grew four Trifolium repens (white clover) F2 crosses and one variety in a conventional field trial and compared differences in root nodule Rhizobium leguminosarum symbiovar trifolii (Rlt) genotype diversity using high-throughput amplicon sequencing of chromosomal housekeeping (rpoB and recA) genes and auxiliary plasmid-borne symbiosis genes (nodA and nodD). We found that Rlt nodule diversities significantly differed between clover crosses, potentially due to host filtering. However, variance in Rlt diversity largely overlapped between crosses and was also explained by the spatial distribution of plants in the field, indicative of the role of local environmental conditions for nodule diversity. To test the effect of host filtering, we conducted a controlled greenhouse trial with a diverse Rlt inoculum and several host genotypes. We found that different clover varieties and genotypes of the same variety selected for significantly different Rlt nodule communities and that the strength of host filtering (deviation from the initial Rhizobium inoculant composition) was positively correlated with the efficiency of symbiosis (rate of plant greenness colouration). Together, our results suggest that selection by host genotype and local growth conditions jointly influence white clover Rlt nodule diversity and community composition.

Variable expression of cyanide detoxification and tolerance genes in cyanogenic and acyanogenic white clover (Trifolium repens).

PREMISE: ß-Cyanoalanine synthase (ß-CAS) and alternative oxidase (AOX) play important roles in the ability of plants to detoxify and tolerate hydrogen cyanide (HCN). These functions are critical for all plants because HCN is produced at low levels during basic metabolic processes, and especially for cyanogenic species, which release high levels of HCN following tissue damage. However, expression of ß-CAS and Aox genes has not been examined in cyanogenic species, nor compared between cyanogenic and acyanogenic genotypes within a species. METHODS: We used a natural polymorphism for cyanogenesis in white clover to examine ß-CAS and Aox gene expression in relation to cyanogenesis-associated HCN exposure. We identified all ß-CAS and Aox gene copies present in the genome, including members of the Aox1, Aox2a, and Aox2d subfamilies previously reported in legumes. Expression levels were compared between cyanogenic and acyanogenic genotypes and between damaged and undamaged leaf tissue. RESULTS: ß-CAS and Aox2a expression was differentially elevated in cyanogenic genotypes, and tissue damage was not required to induce this increased expression. Aox2d, in contrast, appeared to be upregulated as a generalized wounding response. CONCLUSIONS: These findings suggest a heightened constitutive role for HCN detoxification (via elevated ß-CAS expression) and HCN-toxicity mitigation (via elevated Aox2a expression) in plants that are capable of cyanogenesis. As such, freezing-induced cyanide autotoxicity is unlikely to be the primary selective factor in the evolution of climate-associated cyanogenesis clines.

Biochemical changes after cold acclimation in Nordic red clover (Trifolium pratense L.) accessions with contrasting levels of freezing tolerance.

The ability to tolerate low freezing temperatures is an important component of winter survival and persistence of red clover. Cold acclimation (CA) allows plants to acquire higher levels of freezing tolerance. However, the biochemical responses to cold and the importance of such changes for the plant to acquire adequate freezing tolerance have not been investigated in red clover of Nordic origin, which has a distinct genetic background. To shed light on this, we selected five freezing tolerant (FT) and five freezing susceptible (FS) accessions and studied the effect of CA on the contents of carbohydrates, amino acids, and phenolic compounds in the crowns. Among those compounds which increased during CA, FT accessions had higher contents of raffinose, pinitol, arginine, serine, alanine, valine, phenylalanine, and one phenolic compound (a pinocembrin hexoside derivative) than FS accessions, suggesting a role for these compounds in the freezing tolerance in the selected accessions. These findings, together with a description of the phenolic profile of red clover crowns, significantly add to the current knowledge of the biochemical changes during CA and their role in freezing tolerance in Nordic red clover.

Red clover root-associated microbiota is shaped by geographic location and choice of farming system.

AIMS: This study evaluated the red clover (Trifolium pratense) root-associated microbiota to clarify the presence of pathogenic and beneficial microorganisms in 89 Swedish field sites. METHODS AND RESULTS: 16S rRNA and ITS amplicon sequencing analysis were performed on DNA extracted from the red clover root samples collected to determine the composition of the prokaryotic and eukaryotic root-associated microbe communities. Alpha and beta diversities were calculated and relative abundance of various microbial taxa and their co-occurrence were analyzed. Rhizobium was the most prevalent bacterial genus, followed by Sphingomonas, Mucilaginibacter, Flavobacterium, and the unclassified Chloroflexi group KD4-96. The Leptodontidium, Cladosporium, Clonostachys, and Tetracladium fungal genera known for endophytic, saprotrophic, and mycoparasitic lifestyles were also frequently observed in all samples. Sixty-two potential pathogenic fungi were identified with a bias toward grass pathogens and a higher abundance in samples from conventional farms. CONCLUSIONS: We showed that the microbial community was mainly shaped by geographic location and management procedures. Co-occurrence networks revealed that the Rhizobiumleguminosarum bv. trifolii was negatively associated with all fungal pathogenic taxa recognized in this study.

Optimization of Protocols for the Induction of Callus and Plant Regeneration in White Clover (Trifolium repens L.).

White clover is a widely grown temperate legume forage with high nutritional value. Research on the functional genomics of white clover requires a stable and efficient transformation system. In this study, we successfully induced calluses from the cotyledons and leaves of 10 different white clover varieties. The results showed that the callus formation rate in the cotyledons did not vary significantly among the varieties, but the highest callus formation rate was observed in 'Koala' leaves. Subsequently, different concentrations of antioxidants and hormones were tested on the browning rate and differentiation ability of the calluses, respectively. The results showed that the browning rate was the lowest on MS supplemented with 20 mg L-1 AgNO3 and 25 mg L-1 VC, respectively, and the differentiation rate was highest on MS supplemented with 1 mg L-1 6-BA, 1 mg L-1 KT and 0.5 mg L-1 NAA. In addition, the transformation system for Agrobacterium tumefaciens-mediated transformation of 4-day-old leaves was optimized to some extent and obtained a positive callus rate of 8.9% using green fluorescent protein (GFP) as a marker gene. According to our data, by following this optimized protocol, the transformation efficiency could reach 2.38%. The results of this study will provide the foundation for regenerating multiple transgenic white clover from a single genetic background.

Identification and Evolutionary Analysis of the Auxin Response Factor (ARF) Family Based on Transcriptome Data from Caucasian Clover and Analysis of Expression Responses to Hormones.

Caucasian clover (Trifolium ambiguum M. Bieb.) is an excellent perennial plant in the legume family Fabaceae, with a well-developed rhizome and strong clonal growth. Auxin is one of the most important phytohormones in plants and plays an important role in plant growth and development. Auxin response factor (ARF) can regulate the expression of auxin-responsive genes, thus participating in multiple pathways of auxin transduction signaling in a synergistic manner. No genomic database has been established for Caucasian clover. In this study, 71 TaARF genes were identified through a transcriptomic database of Caucasian clover rhizome development. Phylogenetic analysis grouped the TaARFs into six (1-6) clades. Thirty TaARFs contained a complete ARF structure, including three relatively conserved regions. Physical and chemical property analysis revealed that TaARFs are unstable and hydrophilic proteins. We also analyzed the expression pattern of TaARFs in different tissues (taproot, horizontal rhizome, swelling of taproot, rhizome bud and rhizome bud tip). Quantitative real-time RT-PCR revealed that all TaARFs were responsive to phytohormones (indole-3-acetic acid, gibberellic acid, abscisic acid and methyl jasmonate) in roots, stems and leaves. These results helped elucidate the role of ARFs in responses to different hormone treatments in Caucasian clover.

A de novo assembled high-quality chromosome-scale Trifolium pratense genome and fine-scale phylogenetic analysis.

BACKGROUND: Red clover (Trifolium pratense L.) is a diploid perennial temperate legume with 14 chromosomes (2n = 14) native to Europe and West Asia, with high nutritional and economic value. It is a very important forage grass and is widely grown in marine climates, such as the United States and Sweden. Genetic research and molecular breeding are limited by the lack of high-quality reference genomes. In this study, we used Illumina, PacBio HiFi, and Hi-C to obtain a high-quality chromosome-scale red clover genome and used genome annotation results to analyze evolutionary relationships among related species. RESULTS: The red clover genome obtained by PacBio HiFi assembly sequencing was 423 M. The assembly quality was the highest among legume genome assemblies published to date. The contig N50 was 13 Mb, scaffold N50 was 55 Mb, and BUSCO completeness was 97.9%, accounting for 92.8% of the predicted genome. Genome annotation revealed 44,588 gene models with high confidence and 52.81% repetitive elements in red clover genome. Based on a comparison of genome annotation results, red clover was closely related to Trifolium medium and distantly related to Glycine max, Vigna radiata, Medicago truncatula, and Cicer arietinum among legumes. Analyses of gene family expansions and contractions and forward gene selection revealed gene families and genes related to environmental stress resistance and energy metabolism. CONCLUSIONS: We report a high-quality de novo genome assembly for the red clover at the chromosome level, with a substantial improvement in assembly quality over those of previously published red clover genomes. These annotated gene models can provide an important resource for molecular genetic breeding and legume evolution studies. Furthermore, we analyzed the evolutionary relationships among red clover and closely related species, providing a basis for evolutionary studies of clover leaf and legumes, genomics analyses of forage grass, the improvement of agronomic traits.

Asynapsis and unreduced gamete formation in a Trifolium interspecific hybrid.

BACKGROUND: Unreduced gametes, a driving force in the widespread polyploidization and speciation of flowering plants, occur relatively frequently in interspecific or intergeneric hybrids. Studies of the mechanisms leading to 2n gamete formation, mainly in the wheat tribe Triticeae have shown that unreductional meiosis is often associated with chromosome asynapsis during the first meiotic division. The present study explored the mechanisms of meiotic nonreduction leading to functional unreduced gametes in an interspecific Trifolium (clover) hybrid with three sub-genomes from T. ambiguum and one sub-genome from T. occidentale. RESULTS: Unreductional meiosis leading to 2n gametes occurred when there was a high frequency of asynapsis during the first meiotic division. In this hybrid, approximately 39% of chromosomes were unpaired at metaphase I. Within the same cell at anaphase I, sister chromatids of univalents underwent precocious separation and formed laggard chromatids whereas paired chromosomes segregated without separation of sister chromatids as in normal meiosis. This asynchrony was frequently accompanied by incomplete or no movement of chromosomes toward the poles and restitution leading to unreduced chromosome constitutions. Reductional meiosis was restored in progeny where asynapsis frequencies were low. Two progeny plants with approximately 5 and 7% of unpaired chromosomes at metaphase I showed full restoration of reductional meiosis. CONCLUSIONS: The study revealed that formation of 2n gametes occurred when asynapsis (univalent) frequency at meiosis I was high, and that normal gamete production was restored in the next generation when asynapsis frequencies were low. Asynapsis-dependent 2n gamete formation, previously supported by evidence largely from wheat and its relatives and grasshopper, is also applicable to hybrids from the dicotyledonous plant genus Trifolium. The present results align well with those from these widely divergent organisms and strongly suggest common molecular mechanisms involved in unreduced gamete formation.

Identification of loci controlling timing of stem elongation in red clover using genotyping by sequencing of pooled phenotypic extremes.

MAIN CONCLUSION: Through selective genotyping of pooled phenotypic extremes, we identified a number of loci and candidate genes putatively controlling timing of stem elongation in red clover. We have identified candidate genes controlling the timing of stem elongation prior to flowering in red clover (Trifolium pratense L.). This trait is of ecological and agronomic significance, as it affects fitness, competitivity, climate adaptation, forage and seed yield, and forage quality. We genotyped replicate pools of phenotypically extreme individuals (early and late-elongating) within cultivar Lea using genotyping-by-sequencing in pools (pool-GBS). After calling and filtering SNPs and GBS locus haplotype polymorphisms, we estimated allele frequencies and searched for markers with significantly different allele frequencies in the two phenotypic groups using BayeScan, an FST-based test utilizing replicate pools, and a test based on error variance of replicate pools. Of the three methods, BayeScan was the least stringent, and the error variance-based test the most stringent. Fifteen significant markers were identified in common by all three tests. The candidate genes flanking the markers include genes with potential roles in the vernalization, autonomous, and photoperiod regulation of floral transition, hormonal regulation of stem elongation, and cell growth. These results provide a first insight into the potential genes and mechanisms controlling transition to stem elongation in a perennial legume, which lays a foundation for further functional studies of the genetic determinants regulating this important trait.

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 F2  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 trade-offs and unexpected life history traits shape local adaptation in Trifolium repens.

Local adaptation has evolved in numerous taxa across the tree of life in response to divergent selection acting on populations that inhabit different environments (Briscoe Runquist et al., 2019; Hargreaves et al., 2019; Hereford, 2009). The genetic basis of local adaptation has intrigued researchers for decades. In their foundational reciprocal transplant study, Clausen and Hiesey (1960) evaluated the genetics of adaptation primarily using hybrid lines derived from crosses of low elevation and alpine Potentilla glandulosa (Rosaceae) ecotypes. Their work revealed that transgressive segregation can lead to a wider range of trait values than expressed by the parents, complex traits are often genetically correlated and evolve in tandem, and local adaptation is typically polygenic, that is, controlled by many loci of small effect (Clausen & Hiesey, 1960). Within the past 15 years, a burgeoning literature has investigated whether local adaptation evolves through genetic trade-offs, such that local alleles at a QTL (quantitative trait locus) or candidate gene have a fitness advantage in their home environment, but suffer a fitness cost when transplanted into a contrasting habitat type (Mitchell-Olds et al., 2007). Alternatively, local adaptation could arise through conditional neutrality, in which an allele native to one habitat type has elevated fitness in its home site relative to foreign alleles, but is not at a fitness disadvantage elsewhere. Conditional neutrality could maintain local adaptation if gene flow is spatially restricted (Hall et al., 2010). In a From the Cover article in this issue of Molecular Ecology, Wright et al. (2022) examined the genetic basis of local adaptation in white clover (Trifolium repens), discovering strong signatures of local adaptation, and revealing that both genetic trade-offs and conditional neutrality contribute to local adaptation. The most surprising and intriguing result to emerge from this study was that variation in a key antiherbivore defence did not appear to influence contemporary patterns of local adaptation. Rather, divergence in life history strategies was crucial, with early reproduction favoured in the southern garden and delayed reproduction and longer lifespans emerging in the north. These findings highlight the challenges of identifying the multivariate targets of divergent selection in locally-adapted systems, and reveal that not all traits that vary across populations contribute to adaptive differentiation. As studies continue to investigate local adaptation, experiments that manipulate environmental conditions and quantify the magnitude and direction of selection on traits will shed light on the processes that drive local adaptation.

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.

Major effect loci for plant size before onset of nitrogen fixation allow accurate prediction of yield in white clover.

KEY MESSAGE: Accurate genomic prediction of yield within and across generations was achieved by estimating the genetic merit of individual white clover genotypes based on extensive genetic replication using cloned material. White clover is an agriculturally important forage legume grown throughout temperate regions as a mixed clover-grass crop. It is typically cultivated with low nitrogen input, making yield dependent on nitrogen fixation by rhizobia in root nodules. Here, we investigate the effects of clover and rhizobium genetic variation by monitoring plant growth and quantifying dry matter yield of 704 combinations of 145 clover genotypes and 170 rhizobium inocula. We find no significant effect of rhizobium variation. In contrast, we can predict yield based on a few white clover markers strongly associated with plant size prior to nitrogen fixation, and the prediction accuracy for polycross offspring yield is remarkably high. Several of the markers are located near a homolog of Arabidopsis thaliana GIGANTUS 1, which regulates growth rate and biomass accumulation. Our work provides fundamental insight into the genetics of white clover yield and identifies specific candidate genes as breeding targets.

Phenotypic variation and quantitative trait loci for resistance to southern anthracnose and clover rot in red clover.

KEY MESSAGE: High variability for and candidate loci associated with resistance to southern anthracnose and clover rot in a worldwide collection of red clover provide a first basis for genomics-assisted breeding. Red clover (Trifolium pratense L.) is an important forage legume of temperate regions, particularly valued for its high yield potential and its high forage quality. Despite substantial breeding progress during the last decades, continuous improvement of cultivars is crucial to ensure yield stability in view of newly emerging diseases or changing climatic conditions. The high amount of genetic diversity present in red clover ecotypes, landraces, and cultivars provides an invaluable, but often unexploited resource for the improvement of key traits such as yield, quality, and resistance to biotic and abiotic stresses. A collection of 397 red clover accessions was genotyped using a pooled genotyping-by-sequencing approach with 200 plants per accession. Resistance to the two most pertinent diseases in red clover production, southern anthracnose caused by Colletotrichum trifolii, and clover rot caused by Sclerotinia trifoliorum, was assessed using spray inoculation. The mean survival rate for southern anthracnose was 22.9% and the mean resistance index for clover rot was 34.0%. Genome-wide association analysis revealed several loci significantly associated with resistance to southern anthracnose and clover rot. Most of these loci are in coding regions. One quantitative trait locus (QTL) on chromosome 1 explained 16.8% of the variation in resistance to southern anthracnose. For clover rot resistance we found eight QTL, explaining together 80.2% of the total phenotypic variation. The SNPs associated with these QTL provide a promising resource for marker-assisted selection in existing breeding programs, facilitating the development of novel cultivars with increased resistance against two devastating fungal diseases of red clover.

Transcriptome Analysis of Air Space-Type Variegation Formation in Trifolium pratense.

Air space-type variegation is the most diverse among the species of known variegated leaf plants and is caused by conspicuous intercellular spaces between the epidermal and palisade cells and among the palisade cells at non-green areas. Trifolium pratense, a species in Fabaceae with V-shaped air space-type variegation, was selected to explore the application potential of variegated leaf plants and accumulate basic data on the molecular regulatory mechanism and evolutionary history of leaf variegation. We performed comparative transcriptome analysis on young and adult leaflets of variegated and green plants and identified 43 candidate genes related to air space-type variegation formation. Most of the genes were related to cell-wall structure modification (CESA, CSL, EXP, FLA, PG, PGIP, PLL, PME, RGP, SKS, and XTH family genes), followed by photosynthesis (LHCB subfamily, RBCS, GOX, and AGT family genes), redox (2OG and GSH family genes), and nitrogen metabolism (NodGS family genes). Other genes were related to photooxidation, protein interaction, and protease degradation systems. The downregulated expression of light-responsive LHCB subfamily genes and the upregulated expression of the genes involved in cell-wall structure modification were important conditions for air space-type variegation formation in T. pratense. The upregulated expression of the ubiquitin-protein ligase enzyme (E3)-related genes in the protease degradation systems were conducive to air space-type variegation formation. Because these family genes are necessary for plant growth and development, the mechanism of the leaf variegation formation in T. pratense might be a widely existing regulation in air space-type variegation in nature.