Harnessing cold adaptation for postglacial colonisation: Galactinol synthase expression and raffinose accumulation in a polyploid and its progenitors.

Allotetraploid white clover (Trifolium repens) formed during the last glaciation through hybridisation of two European diploid progenitors from restricted niches: one coastal, the other alpine. Here, we examine which hybridisation-derived molecular events may have underpinned white clover's postglacial niche expansion. We compared the transcriptomic frost responses of white clovers (an inbred line and an alpine-adapted ecotype), extant descendants of its progenitor species and a resynthesised white clover neopolyploid to identify genes that were exclusively frost-induced in the alpine progenitor and its derived subgenomes. From these analyses we identified galactinol synthase, the rate-limiting enzyme in biosynthesis of the cryoprotectant raffinose, and found that the extant descendants of the alpine progenitor as well as the neopolyploid white clover rapidly accumulated significantly more galactinol and raffinose than the coastal progenitor under cold stress. The frost-induced galactinol synthase expression and rapid raffinose accumulation derived from the alpine progenitor likely provided an advantage during early postglacial colonisation for white clover compared to its coastal progenitor.

A novel mutation in IAA16 is associated with dicamba resistance in Chenopodium album.

BACKGROUND: Resistance to dicamba in Chenopodium album was first documented over a decade ago, however, the molecular basis of dicamba resistance in this species has not been elucidated. In this research, the resistance mechanism in a dicamba-resistant C. album phenotype was investigated using a transcriptomics (RNA-sequence) approach. RESULTS: The dose-response assay showed that the resistant (R) phenotype was nearly 25-fold more resistant to dicamba than a susceptible (S) phenotype of C. album. Also, dicamba treatment significantly induced transcription of the known auxin-responsive genes, Gretchen Hagen 3 (GH3), small auxin-up RNAs (SAURs), and 1-aminocyclopropane-1-carboxylate synthase (ACS) genes in the susceptible phenotype. Comparing the transcripts of auxin TIR/AFB receptors and auxin/indole-3-acetic acid (AUX/IAA) proteins identified from C. album transcriptomic analysis revealed that the R phenotype contained a novel mutation at the first codon of the GWPPV degron motif of IAA16, resulting in an amino acid substitution of glycine (G) with aspartic acid (D). Sequencing the IAA16 gene in other R and S individuals further confirmed that all the R individuals contained the mutation. CONCLUSION: In this research, we describe the dicamba resistance mechanism in the only case of dicamba-resistant C. album reported to date. Prior work has shown that the dicamba resistance allele confers significant growth defects to the R phenotype investigated here, suggesting that dicamba-resistant C. album carrying this novel mutation in the IAA16 gene may not persist at high frequencies upon removal of dicamba application. © 2024 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

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.

Quantitative genetic analysis reveals potential to breed for improved white clover growth in symbiosis with nitrogen-fixing Rhizobium bacteria.

White clover (Trifolium repens) is integral to mixed pastures in New Zealand and temperate agriculture globally. It provides quality feed and a sustainable source of plant-available nitrogen (N) via N-fixation through symbiosis with soil-dwelling Rhizobium bacteria. Improvement of N-fixation in white clover is a route to enhancing sustainability of temperate pasture production. Focussing on seedling growth critical for crop establishment and performance, a population of 120 half-sibling white clover families was assessed with either N-supplementation or N-fixation via inoculation with a commercial Rhizobium strain (TA1). Quantitative genetic analysis identified significant (p < 0.05) family additive genetic variance for Shoot and Root Dry Matter (DM) and Symbiotic Potential (SP), and Root to Shoot ratio. Estimated narrow-sense heritabilities for above-ground symbiotic traits were moderate (0.24-0.33), and the strong (r ≥ 0.97) genetic correlation between Shoot and Root DM indicated strong pleiotropy or close linkage. The moderate (r = 0.47) phenotypic correlation between Shoot DM under symbiosis vs. under N-supplementation suggested plant growth with mineral-N was not a strong predictor of symbiotic performance. At 5% among-family selection pressure, predicted genetic gains per selection cycle of 19 and 17% for symbiotic traits Shoot DM and Shoot SP, respectively, highlighted opportunities for improved early seedling establishment and growth under symbiosis. Single and multi-trait selection methods, including a Smith-Hazel index focussing on an ideotype of high Shoot DM and Shoot SP, showed commonality of top-ranked families among traits. This study provides a platform for proof-of-concept crosses to breed for enhanced seedling growth under Rhizobium symbiosis and is informative for other legume crops.

Pinoxaden resistance in Lolium perenne L. is due to both target-site and non-target-site mechanisms.

Application of herbicides inhibiting acetyl CoA carboxylase (ACCase) has been one of the main strategies for selectively controlling grass weed species such as perennial ryegrass (Lolium perenne L.) in wheat and barley crops in New Zealand. In this study, we have confirmed and characterized resistance to pinoxaden, an ACCase-inhibiting herbicide, in a population of L. perenne. Dose-response experiments were conducted to assess the level of pinoxaden resistance, and based on the LD50 values, the studied population was 41.4-times more resistant to pinoxaden than a susceptible population. Application of malathion, an inhibitor of the cytochrome P450s, preceding pinoxaden treatment reduced the level of resistance to 9.7-fold. However, pre-treatment with the glutathione S-transferase (GST) inhibitor 4-chloro7- nitrobenzoxadiazole prior to pinoxaden treatment did not affect pinoxaden resistance. Partial sequencing of the ACCase gene revealed that the resistant population had an isoleucine to valine replacement at position 2041. These results suggest that both cytochrome P450-based and target-site mechanisms are jointly associated with this instance of pinoxaden resistance in L. perenne. The pinoxaden-resistant L. perenne individuals were also resistant to quizalofop-p-ethyl (108.6-fold), but they were susceptible to clethodim, which can, therefore, be used to manage this pinoxaden-resistant L. perenne. This is the first report of a L. perenne population in which a rare target-site mutation works in concert with enhanced cytochrome P-450 activity to confer pinoxaden resistance. Evolution of resistance to ACCase-inhibiting herbicides in this L. perenne population indicates that integrated weed management practices are required to prevent widespread resistance developing in New Zealand cereal crop systems.

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-FST. 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.

Deterministic and stochastic modelling of impacts from genomic selection and phenomics on genetic gain for perennial ryegrass dry matter yield.

Increasing the efficiency of current forage breeding programs through adoption of new technologies, such as genomic selection (GS) and phenomics (Ph), is challenging without proof of concept demonstrating cost effective genetic gain (∆G). This paper uses decision support software DeltaGen (tactical tool) and QU-GENE (strategic tool), to model and assess relative efficiency of five breeding methods. The effect on ∆G and cost ($) of integrating GS and Ph into an among half-sib (HS) family phenotypic selection breeding strategy was investigated. Deterministic and stochastic modelling were conducted using mock data sets of 200 and 1000 perennial ryegrass HS families using year-by-season-by-location dry matter (DM) yield data and in silico generated data, respectively. Results demonstrated short (deterministic)- and long-term (stochastic) impacts of breeding strategy and integration of key technologies, GS and Ph, on ∆G. These technologies offer substantial improvements in the rate of ∆G, and in some cases improved cost-efficiency. Applying 1% within HS family GS, predicted a 6.35 and 8.10% ∆G per cycle for DM yield from the 200 HS and 1000 HS, respectively. The application of GS in both among and within HS selection provided a significant boost to total annual ∆G, even at low GS accuracy rA of 0.12. Despite some reduction in ∆G, using Ph to assess seasonal DM yield clearly demonstrated its impact by reducing cost per percentage ∆G relative to standard DM cuts. Open-source software tools, DeltaGen and QuLinePlus/QU-GENE, offer ways to model the impact of breeding methodology and technology integration under a range of breeding scenarios.

A PCR plus restriction enzyme-based technique for detecting target-enzyme mutations at position Pro-106 in glyphosate-resistant Lolium perenne.

The first step in managing herbicide-resistant weeds is to confirm their resistance status. It is, therefore, crucial to have a rapid, reliable and cost-effective technique to assess samples for herbicide resistance. We designed and evaluated three derived cleaved amplified polymorphic sequence (dCAPS) markers for detecting glyphosate resistance in Lolium perenne. conferred by non-synonymous mutations at codon-106 in the enolpyruvylshikimate-3-phosphate synthase (EPSPS) gene. The dCAPS markers involve amplification of the target region, digestion of the amplified products with restriction enzymes and gel-based visualisation of the digested products. The results showed that all three dCAPS markers could successfully detect mutations at codon-106 in the target enzyme. The dCAPS markers can also inform us of the zygosity state of the resistance allele and was confirmed by sequencing the target region of the EPSPS gene. The markers described here are effective quick tests for the monitoring and evaluation of the target-enzyme mechanism of glyphosate resistance in Lolium perenne.

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.

Protocol: a versatile, inexpensive, high-throughput plant genomic DNA extraction method suitable for genotyping-by-sequencing.

BACKGROUND: The recent development of next-generation sequencing DNA marker technologies, such as genotyping-by-sequencing (GBS), generates thousands of informative single nucleotide polymorphism markers in almost any species, regardless of genomic resources. This enables poorly resourced or "orphan" crops/species access to high-density, high-throughput marker platforms which have revolutionised population genetics studies and plant breeding. DNA quality underpins success of GBS methods as the DNA must be amenable to restriction enzyme digestion and sequencing. A barrier to implementing GBS technologies is access to inexpensive, high-throughput extraction methods that yield sequencing-quality genomic DNA (gDNA) from plants. Several high-throughput DNA extraction methods are available, but typically provide low yield or poor quality gDNA, or are costly (US$6-$9/sample) for consumables. RESULTS: We modified a non-organic solvent protocol to extract microgram quantities (1-13 µg) of sequencing-quality high molecular weight gDNA inexpensively in 96-well plates from either fresh, freeze-dried or silica gel-dried plant tissue. The protocol was effective for several easy and difficult-to-extract forage, crop, horticultural and common model species including Trifolium, Medicago, Lolium, Secale, Festuca, Malus, Oryza, and Arabidopsis. The extracted DNA was of high molecular weight and digested readily with restriction enzymes. Contrasting with other extraction protocols we assessed, Illumina-based sequencing of GBS libraries developed from this gDNA had very uniform high quality base-calls to the end of sequence reads. Furthermore, DNA extracted using this method has been sequenced successfully with the PacBio long-read platform. The protocol is scalable, readily automated without requirement for fume hoods, requires approximately three hours to process 192 samples (384-576 samples/day), and is inexpensive at US$0.62/sample for consumables. CONCLUSIONS: This versatile, scalable and simple protocol yields high molecular weight genomic DNA suitable for restriction enzyme digestion and next-generation sequencing applications including GBS and long-read sequencing platforms such as PacBio. The low cost, high-throughput, and extraction of high quality gDNA from a range of fresh and dried source plant material makes this method suitable for many sequencing and genotyping applications including large-scale sample screening underpinning breeding programmes.

Predictive ability of genomic selection models in a multi-population perennial ryegrass training set using genotyping-by-sequencing.

KEY MESSAGE: Genomic prediction models for multi-year dry matter yield, via genotyping-by-sequencing in a composite training set, demonstrate potential for genetic gain improvement through within-half sibling family selection. Perennial ryegrass (Lolium perenne L.) is a key source of nutrition for ruminant livestock in temperate environments worldwide. Higher seasonal and annual yield of herbage dry matter (DMY) is a principal breeding objective but the historical realised rate of genetic gain for DMY is modest. Genomic selection was investigated as a tool to enhance the rate of genetic gain. Genotyping-by-sequencing (GBS) was undertaken in a multi-population (MP) training set of five populations, phenotyped as half-sibling (HS) families in five environments over 2 years for mean herbage accumulation (HA), a measure of DMY potential. GBS using the ApeKI enzyme yielded 1.02 million single-nucleotide polymorphism (SNP) markers from a training set of n = 517. MP-based genomic prediction models for HA were effective in all five populations, cross-validation-predictive ability (PA) ranging from 0.07 to 0.43, by trait and target population, and 0.40-0.52 for days-to-heading. Best linear unbiased predictor (BLUP)-based prediction methods, including GBLUP with either a standard or a recently developed (KGD) relatedness estimation, were marginally superior or equal to ridge regression and random forest computational approaches. PA was principally an outcome of SNP modelling genetic relationships between training and validation sets, which may limit application for long-term genomic selection, due to PA decay. However, simulation using data from the training experiment indicated a twofold increase in genetic gain for HA, when applying a prediction model with moderate PA in a single selection cycle, by combining among-HS family selection, based on phenotype, with within-HS family selection using genomic prediction.

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.

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

BACKGROUND: White clover (Trifolium repens L.) is a temperate forage legume with an allotetraploid genome (2n=4×=32) estimated at 1093 Mb. Several linkage maps of various sizes, marker sources and completeness are available, however, no integrated map and marker set has explored consistency of linkage analysis among unrelated mapping populations. Such integrative analysis requires tools for homoeologue matching among populations. Development of these tools provides for a consistent framework map of the white clover genome, and facilitates in silico alignment with the model forage legume, Medicago truncatula. RESULTS: This is the first report of integration of independent linkage maps in white clover, and adds to the literature on methyl filtered GeneThresher®-derived microsatellite (simple sequence repeat; SSR) markers for linkage mapping. Gene-targeted SSR markers were discovered in a GeneThresher® (TrGT) methyl-filtered database of 364,539 sequences, which yielded 15,647 SSR arrays. Primers were designed for 4,038 arrays and of these, 465 TrGT-SSR markers were used for parental consensus genetic linkage analysis in an F1 mapping population (MP2). This was merged with an EST-SSR consensus genetic map of an independent population (MP1), using markers to match homoeologues and develop a multi-population integrated map of the white clover genome. This integrated map (IM) includes 1109 loci based on 804 SSRs over 1274 cM, covering 97% of the genome at a moderate density of one locus per 1.2 cM. Eighteen candidate genes and one morphological marker were also placed on the IM. Despite being derived from disparate populations and marker sources, the component maps and the derived IM had consistent representations of the white clover genome for marker order and genetic length. In silico analysis at an E-value threshold of 1e-20 revealed substantial co-linearity with the Medicago truncatula genome, and indicates a translocation between T. repens groups 2 and 6 relative to M. truncatula. CONCLUSIONS: This integrated genetic linkage analysis provides a consistent and comprehensive linkage analysis of the white clover genome, with alignment to a model forage legume. Associated marker locus information, particularly the homoeologue-specific markers, offers a new resource for forage legume research to enable genetic analysis and improvement of this forage and grassland species.

A lineage-specific centromeric satellite sequence in the genus Trifolium.

We report the molecular structure, genomic organization, chromosomal distribution and evolutionary dynamics of TrR350, a satellite DNA isolated from the forage legume white clover (Trifolium repens L.; 2n = 4 x = 32). The basic repeating unit is an A+T rich 350 bp Hin dIII fragment with a complex dimeric structure consisting of an internal direct repeat of 156 bp packed between unrelated flanking sequences. Each 156 bp repeat has a conserved 24 bp motif repeating at two places. Most of the 24 bp short repeating units enclose a pentanucleotide CAAAA motif, presumed to be involved in breakage-reunion mechanism of tandemly repeating arrays. The dimers share high sequence homology among themselves while monomers within dimers show significant sequence divergence. Genomic Southern hybridization and/or fluorescence in situ hybridization (FISH) on 17 Trifolium species/subspecies revealed that it is a lineage-specific repeat confined to several species within the section Lotoidea originating in the Mediterranean region. The uniform length of the basic repeating unit and the centromeric localization in most of the species harbouring it reflects its extensive conservation in the lineage. However, the Hin dIII restriction profile in seven species also indicated independent evolution of this repeat.