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

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.

Haplotype-Resolved, Chromosome-Level Assembly of White Clover (Trifolium repens L., Fabaceae)

White clover (Trifolium repens L.; Fabaceae) is an important forage and cover crop in agricultural pastures around the world and is increasingly used in evolutionary ecology and genetics to understand the genetic basis of adaptation. Historically, improvements in white clover breeding practices and assessments of genetic variation in nature have been hampered by a lack of high-quality genomic resources for this species, owing in part to its high heterozygosity and allotetraploid hybrid origin. Here, we use PacBio HiFi and chromosome conformation capture (Omni-C) technologies to generate a chromosome-level, haplotype-resolved genome assembly for white clover totaling 998 Mbp (scaffold N50 = 59.3 Mbp) and 1 Gbp (scaffold N50 = 58.6 Mbp) for haplotypes 1 and 2, respectively, with each haplotype arranged into 16 chromosomes (8 per subgenome). We additionally provide a functionally annotated haploid mapping assembly (968 Mbp, scaffold N50 = 59.9 Mbp), which drastically improves on the existing reference assembly in both contiguity and assembly accuracy. We annotated 78,174 protein-coding genes, resulting in protein BUSCO completeness scores of 99.6% and 99.3% against the embryophyta_odb10 and fabales_odb10 lineage datasets, respectively.

Genome-wide analysis of the WRKY genes and their important roles during cold stress in white clover

Background

White clover (Trifolium repens L) is a high-quality forage grass with a high protein content, but it is vulnerable to cold stress, which can negatively affect its growth and development. WRKY transcription factor is a family of plant transcription factors found mainly in higher plants and plays an important role in plant growth, development, and stress response. Although WRKY transcription factors have been studied extensively in other plants, it has been less studied in white clover.

Methods and Results

In the present research, we have performed a genome-wide analysis of the WRKY gene family of white clover, in total, there were 145 members of WRKY transcription factors identified in white clover. The characterization of the TrWRKY genes was detailed, including conserved motif analysis, phylogenetic analysis, and gene duplication analysis, which have provided a better understanding of the structure and evolution of the TrWRKY genes in white clover. Meanwhile, the genetic regulation network (GRN) containing TrWRKY genes was reconstructed, and Gene Ontology (GO) annotation analysis of these function genes showed they contributed to regulation of transcription process, response to wounding, and phosphorylay signal transduction system, all of which were important processes in response to abiotic stress. To determine the TrWRKY genes function under cold stress, the RNA-seq dataset was analyzed; most of TrWRKY genes were highly upregulated in response to cold stress, particularly in the early stages of cold stress. These results were validated by qRT-PCR experiment, implying they are involved in various gene regulation pathways in response to cold stress.

Conclusion

The results of this study provide insights that will be useful for further functional analyses of TrWRKY genes in response to biotic or abiotic stresses in white clover. These findings are likely to be useful for further research on the functions of TrWRKY genes and their role in response to cold stress, which is important to understand the molecular mechanism of cold tolerance in white clover and improve its cold tolerance.

A genome-wide association study of freezing tolerance in red clover (Trifolium pratense L.) germplasm of European origin

Improvement of persistency is an important breeding goal in red clover (Trifolium pratense L.). In areas with cold winters, lack of persistency is often due to poor winter survival, of which low freezing tolerance (FT) is an important component. We conducted a genome wide association study (GWAS) to identify loci associated with freezing tolerance in a collection of 393 red clover accessions, mostly of European origin, and performed analyses of linkage disequilibrium and inbreeding. Accessions were genotyped as pools of individuals using genotyping-by-sequencing (pool-GBS), generating both single nucleotide polymorphism (SNP) and haplotype allele frequency data at accession level. Linkage disequilibrium was determined as a squared partial correlation between the allele frequencies of pairs of SNPs and found to decay at extremely short distances (< 1 kb). The level of inbreeding, inferred from the diagonal elements of a genomic relationship matrix, varied considerably between different groups of accessions, with the strongest inbreeding found among ecotypes from Iberia and Great Britain, and the least found among landraces. Considerable variation in FT was found, with LT50-values (temperature at which 50% of the plants are killed) ranging from -6.0°C to -11.5°C. SNP and haplotype-based GWAS identified eight and six loci significantly associated with FT (of which only one was shared), explaining 30% and 26% of the phenotypic variation, respectively. Ten of the loci were found within or at a short distance (<0.5 kb) from genes possibly involved in mechanisms affecting FT. These include a caffeoyl shikimate esterase, an inositol transporter, and other genes involved in signaling, transport, lignin synthesis and amino acid or carbohydrate metabolism. This study paves the way for a better understanding of the genetic control of FT and for the development of molecular tools for the improvement of this trait in red clover through genomics assisted breeding.

Genome-wide association study identifies GhSAL1 affects cold tolerance at the seedling emergence stage in upland cotton (Gossypium hirsutum L.)

Genomic analysis of upland cotton revealed that cold tolerance was associated with ecological distribution. GhSAL1 on chromosome D09 negatively regulated cold tolerance of upland cotton. Cotton can undergo low-temperature stress at the seedling emergence stage, which adversely affects growth and yield; however, the regulatory mechanism underlying cold tolerance remains nebulous. Here, we analyze the phenotypic and physiological parameters in 200 accessions from 5 ecological distributions under constant chilling (CC) and diurnal variation of chilling (DVC) stresses at the seedling emergence stage. All accessions were clustered into four groups, of which Group IV, with most germplasms from the northwest inland region (NIR), had better phenotypes than Groups I-III under the two kinds of chilling stresses. A total of 575 significantly associated single-nucleotide polymorphism (SNP) were identified, and 35 stable genetic quantitative trait loci (QTL) were obtained, of which 5 were associated with traits under CC and DVC stress, respectively, while the remaining 25 were co-associated. The accumulation of dry weight (DW) of seedling was associated with the flavonoid biosynthesis process regulated by Gh_A10G0500. The emergence rate (ER), DW, and total length of seedling (TL) under CC stress were associated with the SNPs variation of Gh_D09G0189 (GhSAL1). GhSAL1^HapB was the elite haplotype, which increased ER, DW, and TL by 19.04%, 11.26%, and 7.69%, respectively, compared with that of GhSAL1^HapA. The results of virus-induced gene silencing (VIGS) experiment and determination of metabolic substrate content preliminarily illustrated that GhSAL1 negatively regulated cotton cold tolerance through IP₃-Ca²⁺ signaling pathway. The elite haplotypes and candidate genes identified in this study could be used to improve cold tolerance at the seedling emergence stage in future upland cotton breeding.

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.

Low Temperature Stress Tolerance: An Insight Into the Omics Approaches for Legume Crops

The change in climatic conditions is the major cause for decline in crop production worldwide. Decreasing crop productivity will further lead to increase in global hunger rate. Climate change results in environmental stress which has negative impact on plant-like deficiencies in growth, crop yield, permanent damage, or death if the plant remains in the stress conditions for prolonged period. Cold stress is one of the main abiotic stresses which have already affected the global crop production. Cold stress adversely affects the plants leading to necrosis, chlorosis, and growth retardation. Various physiological, biochemical, and molecular responses under cold stress have revealed that the cold resistance is more complex than perceived which involves multiple pathways. Like other crops, legumes are also affected by cold stress and therefore, an effective technique to mitigate cold-mediated damage is critical for long-term legume production. Earlier, crop improvement for any stress was challenging for scientific community as conventional breeding approaches like inter-specific or inter-generic hybridization had limited success in crop improvement. The availability of genome sequence, transcriptome, and proteome data provides in-depth sight into different complex mechanisms under cold stress. Identification of QTLs, genes, and proteins responsible for cold stress tolerance will help in improving or developing stress-tolerant legume crop. Cold stress can alter gene expression which further leads to increases in stress protecting metabolites to cope up the plant against the temperature fluctuations. Moreover, genetic engineering can help in development of new cold stress-tolerant varieties of legume crop. This paper provides a general insight into the "omics" approaches for cold stress in legume crops.

Genetic control of leaf angle in sorghum and its effect on light interception

Developing sorghum genotypes adapted to different light environments requires understanding of a plant's ability to capture light, determined through leaf angle specifically. This study dissected the genetic basis of leaf angle in 3 year field trials at two sites, using a sorghum diversity panel (729 accessions). A wide range of variation in leaf angle with medium heritability was observed. Leaf angle explained 36% variation in canopy light extinction coefficient, highlighting the extent to which variation in leaf angle influences light interception at the whole-canopy level. This study also found that the sorghum races of Guinea and Durra consistently having the largest and smallest leaf angle, respectively, highlighting the potential role of leaf angle in adaptation to distinct environments. The genome-wide association study detected 33 quantitative trait loci (QTLs) associated with leaf angle. Strong synteny was observed with previously detected leaf angle QTLs in maize (70%) and rice (40%) within 10 cM, among which the overlap was significantly enriched according to χ2 tests, suggesting a highly consistent genetic control in grasses. A priori leaf angle candidate genes identified in maize and rice were found to be enriched within a 1-cM window around the sorghum leaf angle QTLs. Additionally, protein domain analysis identified the WD40 protein domain as being enriched within a 1-cM window around the QTLs. These outcomes show that there is sufficient heritability and natural variation in the angle of upper leaves in sorghum which may be exploited to change light interception and optimize crop canopies for different contexts.

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

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.

Into the weeds: new insights in plant stress

Weeds, plants that thrive in the face of disturbance, have eluded human's attempts at control for >12 000 years, positioning them as a unique group of extreme stress tolerators. The most successful weeds have a suite of traits that enable them to rapidly adapt to environments typified by stress, growing in hostile conditions or subject to massive destruction from agricultural practices. Through their ability to persist and adapt, weeds illuminate principles of evolution and provide insights into weed management and crop improvement. Here we highlight why the time is right to move beyond traditional model systems and leverage weeds to gain a deeper understanding of the mechanisms, adaptations, and genetic and physiological bases for stress tolerance.

Understanding the Effectiveness of Genomic Prediction in Tetraploid Potato

Use of genomic prediction (GP) in tetraploid is becoming more common. Therefore, we think it is the right time for a comparison of GP models for tetraploid potato. GP models were compared that contrasted shrinkage with variable selection, parametric vs. non-parametric models and different ways of accounting for non-additive genetic effects. As a complement to GP, association studies were carried out in an attempt to understand the differences in prediction accuracy. We compared our GP models on a data set consisting of 147 cultivars, representing worldwide diversity, with over 39 k GBS markers and measurements on four tuber traits collected in six trials at three locations during 2 years. GP accuracies ranged from 0.32 for tuber count to 0.77 for dry matter content. For all traits, differences between GP models that utilised shrinkage penalties and those that performed variable selection were negligible. This was surprising for dry matter, as only a few additive markers explained over 50% of phenotypic variation. Accuracy for tuber count increased from 0.35 to 0.41, when dominance was included in the model. This result is supported by Genome Wide Association Study (GWAS) that found additive and dominance effects accounted for 37% of phenotypic variation, while significant additive effects alone accounted for 14%. For tuber weight, the Reproducing Kernel Hilbert Space (RKHS) model gave a larger improvement in prediction accuracy than explicitly modelling epistatic effects. This is an indication that capturing the between locus epistatic effects of tuber weight can be done more effectively using the semi-parametric RKHS model. Our results show good opportunities for GP in 4x potato.

Prospects for Trifolium Improvement Through Germplasm Characterisation and Pre-breeding in New Zealand and Beyond

Trifolium is the most used pastoral legume genus in temperate grassland systems, and a common feature in meadows and open space areas in cities and parks. Breeding of Trifolium spp. for pastoral production has been going on for over a century. However, the breeding targets have changed over the decades in response to different environmental and production pressures. Relatively small gains have been made in Trifolium breeding progress. Trifolium breeding programmes aim to maintain a broad genetic base to maximise variation. New Zealand is a global hub in Trifolium breeding, utilising exotic germplasm imported by the Margot Forde Germplasm Centre. This article describes the history of Trifolium breeding in New Zealand as well as the role and past successes of utilising genebanks in forage breeding. The impact of germplasm characterisation and evaluation in breeding programmes is also discussed. The history and challenges of Trifolium breeding and its effect on genetic gain can be used to inform future pre-breeding decisions in this genus, as well as being a model for other forage legumes.

Genome-wide association study of multiple yield traits in a diversity panel of polyploid sugarcane (Saccharum spp.)

Sugarcane (Saccharum spp.) is an important economic crop, contributing up to 80% of sugar and approximately 60% of biofuel globally. To meet the increased demand for sugar and biofuel supplies, it is critical to breed sugarcane cultivars with robust performance in yield traits. Therefore, dissection of causal DNA sequence variants is of great importance, as it provides genetic resources and fundamental information for crop improvement. In this study, we analyzed nine yield traits in a sugarcane diversity panel consisting of 308 accessions primarily selected from the World Collection of Sugarcane and Related Grasses. By genotyping the diversity panel via target enrichment sequencing, we identified a large number of sequence variants. Genome-wide association studies between the markers and traits were conducted, taking dosages and gene actions into consideration. In total, 217 nonredundant markers and 225 candidate genes were identified to be significantly associated with the yield traits, which can serve as a comprehensive genetic resource database for future gene identification, characterization, and selection for sugarcane improvement. We further investigated runs of homozygosity (ROH) in the sugarcane diversity panel. We characterized 282 ROHs and found that the occurrence of ROHs in the genome were nonrandom and probably under selection. The ROHs were associated with total weight and dry weight, and high ROHs resulted in a decrease in the two traits. This study suggests that genomic inbreeding has led to negative impacts on sugarcane yield.