Identification of Single-Nucleotide Polymorphic Loci Associated with Biomass Yield under Water Deficit in Alfalfa (Medicago sativa L.) Using Genome-Wide Sequencing and Association Mapping

Genome-Wide Association and Genomic Prediction of Alfalfa (Medicago sativa L.) Biomass Yield Under Drought Stress

Developing drought-resistant alfalfa (Medicago sativa L.) that maintains high biomass yield is a key breeding goal to enhance productivity in water-limited areas. In this study, 424 alfalfa breeding families were analyzed to identify molecular markers associated with biomass yield under drought stress and to predict high-merit plants. Biomass yield was measured from 18 harvests from 2020 to 2023 in a field trial with deficit irrigation. A total of 131 significant markers were associated with biomass yield, with 80 markers specifically linked to yield under drought stress; among these, 19 markers were associated with multiple harvests. Finally, genomic best linear unbiased prediction (GBLUP) was employed to obtain predictive accuracies (PAs) and genomic estimated breeding values (GEBVs). Removing low-informative SNPs [SNPs with p-values > 0.05 from the additive Genome-Wide Association (GWAS) model] for GBLUP increased PA by 47.3%. The high number of markers associated with yield under drought stress and the highest PA (0.9) represent a significant achievement in improving yield under drought stress in alfalfa.

GWAS for Drought Resilience Traits in Red Clover (Trifolium pratense L.)

Red clover (Trifolium pratense L.) is a well-appreciated grassland crop in temperate climates but suffers from increasingly frequent and severe drought periods. Molecular markers for drought resilience (DR) would benefit breeding initiatives for red clover, as would a better understanding of the genes involved in DR. Two previous studies, as follows, have: (1) identified phenotypic DR traits in a diverse set of red clover accessions; and (2) produced genotypic data using a pooled genotyping-by-sequencing (GBS) approach in the same collection. In the present study, we performed genome-wide association studies (GWAS) for DR using the available phenotypic and genotypic data. Single nucleotide polymorphism (SNP) calling was performed using GBS data and the following two red clover genome assemblies: the recent HEN-17 assembly and the Milvus assembly. SNP positions with significant associations were used to delineate flanking regions in both genome assemblies, while functional annotations were retrieved from Medicago truncatula orthologs. GWAS revealed 19 significant SNPs in the HEN-17-derived SNP set, explaining between 5.3 and 23.2% of the phenotypic variation per SNP-trait combination for DR traits. Among the genes in the SNP-flanking regions, we identified candidate genes related to cell wall structuring, genes encoding sugar-modifying proteins, an ureide permease gene, and other genes linked to stress metabolism pathways. GWAS revealed 29 SNPs in the Milvus-derived SNP set that explained substantially more phenotypic variation for DR traits, between 5.3 and 42.3% per SNP-trait combination. Candidate genes included a DEAD-box ATP-dependent RNA helicase gene, a P-loop nucleoside triphosphate hydrolase gene, a Myb/SANT-like DNA-binding domain protein, and an ubiquitin-protein ligase gene. Most accessions in this study are genetically more closely related to the Milvus genotype than to HEN-17, possibly explaining how the Milvus-derived SNP set yielded more robust associations. The Milvus-derived SNP set pinpointed 10 genomic regions that explained more than 25% of the phenotypic variation for DR traits. A possible next step could be the implementation of these SNP markers in practical breeding programs, which would help to improve DR in red clover. Candidate genes could be further characterized in future research to unravel drought stress resilience in red clover in more detail.

Pearl millet a promising fodder crop for changing climate: a review

The agricultural sector faces colossal challenges amid environmental changes and a burgeoning human population. In this context, crops must adapt to evolving climatic conditions while meeting increasing production demands. The dairy industry is anticipated to hold the highest value in the agriculture sector in future. The rise in the livestock population is expected to result in an increased demand for fodder feed. Consequently, it is crucial to seek alternative options, as crops demand fewer resources and are resilient to climate change. Pearl millet offers an apposite key to these bottlenecks, as it is a promising climate resilience crop with significantly low energy, water and carbon footprints compared to other crops. Numerous studies have explored its potential as a fodder crop, revealing promising performance. Despite its capabilities, pearl millet has often been overlooked. To date, few efforts have been made to document molecular aspects of fodder-related traits. However, several QTLs and candidate genes related to forage quality have been identified in other fodder crops, which can be harnessed to enhance the forage quality of pearl millet. Lately, excellent genomic resources have been developed in pearl millet allowing deployment of cutting-edge genomics-assisted breeding for achieving a higher rate of genetic gains. This review would facilitate a deeper understanding of various aspects of fodder pearl millet in retrospect along with the future challenges and their solution. This knowledge may pave the way for designing efficient breeding strategies in pearl millet thereby supporting sustainable agriculture and livestock production in a changing world.

A genome-wide association study reveals novel loci and candidate genes associated with plant height variation in Medicago sativa

BACKGROUND

Plant height (PH) is an important agronomic trait influenced by a complex genetic network. However, the genetic basis for the variation in PH in Medicago sativa remains largely unknown. In this study, a comprehensive genome-wide association analysis was performed to identify genomic regions associated with PH using a diverse panel of 220 accessions of M. sativa worldwide.

RESULTS

Our study identified eight novel single nucleotide polymorphisms (SNPs) significantly associated with PH evaluated in five environments, explaining 8.59-12.27% of the phenotypic variance. Among these SNPs, the favorable genotype of chr6__31716285 had a low frequency of 16.4%. Msa0882400, located proximal to this SNP, was annotated as phosphate transporter 3;1, and its role in regulating alfalfa PH was supported by transcriptome and candidate gene association analysis. In addition, 21 candidate genes were annotated within the associated regions that are involved in various biological processes related to plant growth and development.

CONCLUSIONS

Our findings provide new molecular markers for marker-assisted selection in M. sativa breeding programs. Furthermore, this study enhances our understanding of the underlying genetic and molecular mechanisms governing PH variations in M. sativa.

An alfalfa MYB-like transcriptional factor MsMYBH positively regulates alfalfa seedling drought resistance and undergoes MsWAV3-mediated degradation

Drought is a major threat to alfalfa (Medicago sativa L.) production. The discovery of important alfalfa genes regulating drought response will facilitate breeding for drought-resistant alfalfa cultivars. Here, we report a genome-wide association study of drought resistance in alfalfa. We identified and functionally characterized an MYB-like transcription factor gene (MsMYBH), which increases the drought resistance in alfalfa. Compared with the wild-types, the biomass and forage quality were enhanced in MsMYBH overexpressed plants. Combined RNA-seq, proteomics and chromatin immunoprecipitation analysis showed that MsMYBH can directly bind to the promoters of MsMCP1, MsMCP2, MsPRX1A and MsCARCAB to improve their expression. The outcomes of such interactions include better water balance, high photosynthetic efficiency and scavenge excess H2O2 in response to drought. Furthermore, an E3 ubiquitin ligase (MsWAV3) was found to induce MsMYBH degradation under long-term drought, via the 26S proteasome pathway. Furthermore, variable-number tandem repeats in MsMYBH promoter were characterized among a collection of germplasms, and the variation is associated with promoter activity. Collectively, our findings shed light on the functions of MsMYBH and provide a pivotal gene that could be leveraged for breeding drought-resistant alfalfa. This discovery also offers new insights into the mechanisms of drought resistance in alfalfa.

Identification of genetic loci associated with five agronomic traits in alfalfa using multi-environment trials

The use of multi-environment trials to test yield-related traits in a diverse alfalfa panel allowed to find multiple molecular markers associated with complex agronomic traits. Yield is one of the most important target traits in alfalfa breeding; however, yield is a complex trait affected by genetic and environmental factors. In this study, we used multi-environment trials to test yield-related traits in a diverse panel composed of 200 alfalfa accessions and varieties. Phenotypic data of maturity stage measured as mean stage by count (MSC), dry matter content, plant height (PH), biomass yield (Yi), and fall dormancy (FD) were collected in three locations in Idaho, Oregon, and Washington from 2018 to 2020. Single-trial and stagewise analyses were used to obtain estimated trait means of entries by environment. The plants were genotyped using a genotyping by sequencing approach and obtained a genotypic matrix with 97,345 single nucleotide polymorphisms. Genome-wide association studies identified a total of 84 markers associated with the traits analyzed. Of those, 29 markers were in noncoding regions and 55 markers were in coding regions. Ten significant SNPs at the same locus were associated with FD and they were linked to a gene annotated as a nuclear fusion defective 4-like (NFD4). Additional SNPs associated with MSC, PH, and Yi were annotated as transcription factors such as Cysteine3Histidine (C3H), Hap3/NF-YB family, and serine/threonine-protein phosphatase 7 proteins, respectively. Our results provide insight into the genetic factors that influence alfalfa maturity, yield, and dormancy, which is helpful to speed up the genetic gain toward alfalfa yield improvement.

Alfalfa genomic selection for different stress-prone growing regions

Alfalfa (Medicago sativa L.) selection for stress-prone regions has high priority for sustainable crop-livestock systems. This study assessed the genomic selection (GS) ability to predict alfalfa breeding values for drought-prone agricultural sites of Algeria, Morocco, and Argentina; managed-stress (MS) environments of Italy featuring moderate or intense drought; and one Tunisian site irrigated with moderately saline water. Additional aims were to investigate genotype × environment interaction (GEI) patterns and the effect on GS predictions of three single-nucleotide polymorphism (SNP) calling procedures, 12 statistical models that exclude or incorporate GEI, and allele dosage information. Our study included 127 genotypes from a Mediterranean reference population originated from three geographically contrasting populations, genotyped via genotyping-by-sequencing and phenotyped based on multi-year biomass dry matter yield of their dense-planted half-sib progenies. The GEI was very large, as shown by 27-fold greater additive genetic variance × environment interaction relative to the additive genetic variance and low genetic correlation for progeny yield responses across environments. The predictive ability of GS (using at least 37,969 SNP markers) exceeded 0.20 for moderate MS (representing Italian stress-prone sites) and the sites of Algeria and Argentina while being quite low for the Tunisian site and intense MS. Predictions of GS were complicated by rapid linkage disequilibrium decay. The weighted GBLUP model, GEI incorporation into GS models, and SNP calling based on a mock reference genome exhibited a predictive ability advantage for some environments. Our results support the specific breeding for each target region and suggest a positive role for GS in most regions when considering the challenges associated with phenotypic selection.

Identification of QTL and candidate genes associated with biomass yield and Feed Quality in response to water deficit in alfalfa (Medicago sativa L.) using linkage mapping and RNA-Seq

Biomass yield and Feed Quality are the most important traits in alfalfa (Medicago sativa L.), which directly affect its economic value. Drought stress is one of the main limiting factors affecting alfalfa production worldwide. However, the genetic and especially the molecular mechanisms for drought tolerance in alfalfa are poorly understood. In this study, linkage mapping was performed in an F1 population by combining 12 phenotypic data (biomass yield, plant height, and 10 Feed Quality-related traits). A total of 48 significant QTLs were identified on the high-density genetic linkage maps that were constructed in our previous study. Among them, nine main QTLs, which explained more than 10% phenotypic variance, were detected for biomass yield (one), plant height (one), CP (two), ASH (one), P (two), K(one), and Mg (one). A total of 31 candidate genes were identified in the nine main QTL intervals based on the RNA-seq analysis under the drought condition. Blast-P was further performed to screen candidate genes controlling drought tolerance, and 22 functional protein candidates were finally identified. The results of the present study will be useful for improving drought tolerance of alfalfa varieties by marker-assisted selection (MAS), and provide promising candidates for further gene cloning and mechanism study.

Toxicity analysis of TNT to alfalfa's mineral nutrition and secondary metabolism

Alfalfa has the ability to degrade TNT. TNT exposure caused root disruption of mineral nutrient metabolism. The exposure of TNT imbalanced basal cell energy metabolism. The mechanism of 2,4,6-trinitrotoluene (TNT) toxicity effects was analyzed in alfalfa (Medicago sativa L.) seedlings by examining the mineral nutrition and secondary metabolism of the plant roots. Exposure to 25-100 mg·L^-1 TNT in a hydroponic solution for 72 h resulted in a TNT absorption rate of 26.8-63.0%. The contents of S, K, and B in root mineral nutrition metabolism increased significantly by 1.70-5.46 times, 1.38-4.01 times, and 1.40-4.03 times, respectively, after TNT exposure. Non-targeted metabolomics analysis of the roots identified 189 significantly upregulated metabolites and 420 significantly downregulated metabolites. The altered metabolites were primarily lipids and lipid-like molecules, and the most significant enrichment pathways were alanine, aspartate, and glutamate metabolism and glycerophospholipid metabolism. TNT itself was transformed in the root system into several intermediate products, including 4-hydroxylamino-2,6-dinitrotoluene, 4-amino-2,6-dinitrotoluene, 2-hydroxylamino-4,6-dinitrotoluene, 2,4',6,6'-tetranitro-2',4-azoxytoluene, 4,4',6,6'-tetranitro-2,2'-azoxytoluene, and 2,4-dinitrotoluene. Overall, TNT exposure disturbed the mineral metabolism balance, and significantly interfered with basic plant metabolism.

Validation of DNA marker-assisted selection for forage biomass productivity under deficit irrigation in alfalfa

Drought and limited irrigation resources threaten agricultural sustainability in many regions of the world. Application of genomic-based breeding strategies may benefit crop variety development for these environments. Here, we provide a first report on the effect of deploying DNA marker-assisted selection (MAS) for the drought resilience quantitative trait in alfalfa (Medicago sativa L.). The goals of this study were to validate the effect of several quantitative trait loci (QTL) associated with alfalfa forage and crown-root (CR) biomass during drought and to determine their potential to improve forage yield of elite germplasm under water-limited conditions. Marker assisted selection was employed to introgress favorable or unfavorable DNA marker alleles affiliated with 10 biomass QTL into three elite backgrounds. Thirty-two populations were developed and evaluated for forage productivity over 3 yr under continuous deficit irrigation management in New Mexico, USA. Significant yield differences (ranging from -13 to 26%) were detected among some MAS-derived populations in all three elite backgrounds. Application of QTL MAS generally resulted in expected phenotypic responses within an elite genetic background that was similar to that in which the QTL were originally identified. However, relative performance of the populations varied substantially across the three genetic backgrounds. These outcomes indicate that QTL MAS can significantly affect forage productivity of elite alfalfa germplasm in drought-stressed environments. However, if biomass QTL are detected in donor germplasm that is genetically dissimilar to targeted elite populations, characterization of donor alleles may be warranted within elite backgrounds of interest to confirm their phenotypic effects prior to implementing MAS-based breeding.

Genome assembly of alfalfa cultivar zhongmu-4 and identification of SNPs associated with agronomic traits

Alfalfa (Medicago sativa L.) is the most important legume forage crop worldwide with high nutritional value and yield. For a long time, the breeding of alfalfa was hampered by lacking reliable information on the autotetraploid genome and molecular markers linked to important agronomic traits. We herein reported the de novo assembly of the allele-aware chromosome-level genome of Zhongmu-4, a cultivar widely cultivated in China, and a comprehensive database of genomic variations based on resequencing of 220 germplasms. Approximate 2.74 Gb contigs (N50 of 2.06 Mb), accounting for 88.39% of the estimated genome, were assembled, and 2.56 Gb contigs were anchored to 32 pseudo-chromosomes. A total of 34,922 allelic genes were identified from the allele-aware genome. We observed the expansion of gene families, especially those related to the nitrogen metabolism, and the increase of repetitive elements including transposable elements, which probably resulted in the increase of Zhongmu-4 genome compared with Medicago truncatula. Population structure analysis revealed that the accessions from Asia and South America had relatively lower genetic diversity than those from Europe, suggesting that geography may influence alfalfa genetic divergence during local adaption. Genome-wide association studies identified 101 single nucleotide polymorphisms (SNPs) associated with 27 agronomic traits. Two candidate genes were predicted to be correlated with fall dormancy and salt response. We believe that the allele-aware chromosome-level genome sequence of Zhongmu-4 combined with the resequencing data of the diverse alfalfa germplasms will facilitate genetic research and genomics-assisted breeding in variety improvement of alfalfa.

Identification of genetic loci associated with forage quality in response to water deficit in autotetraploid alfalfa (Medicago sativa L.)

BACKGROUND

Alfalfa has been cultivated in many regions around the world as an important forage crop due to its nutritive value to livestock and ability to adapt to various environments. However, the genetic basis by which plasticity of quality-relevant traits influence alfalfa adaption to different water conditions remain largely unknown.

RESULTS

In the present study, 198 accessions of alfalfa of the core collection for drought tolerance were evaluated for 26 forage quality traits in a field trial under an imposed deficit irrigation gradient. Regression analysis between quality traits and water stress revealed that values of fiber-related traits were negatively correlated with values of energy-related traits as water deficit increased. More than one hundred significant markers associated with forage quality under different water treatments were identified using genome-wide association studies with genotyping by sequencing. Among them, 131 markers associated with multiple traits in all the water deficit treatments. Most of the associated markers were dependent to the levels of water deficit, suggesting genetic controls for forage quality traits were dependent to the stress treatment. Twenty-four loci associated with forage quality were annotated to functional genes that may play roles in cell development or in response to water stress.

CONCLUSIONS

This study addressed the genetic base of phenotypic variation of forage quality traits under water deficit. The SNP markers identified in this study will be useful in marker-assisted selection for the genetic improvement of alfalfa with enhanced drought tolerance while maintaining forage quality.

High-density linkage map construction and mapping QTL for yield and yield components in autotetraploid alfalfa using RAD-seq

BACKGROUND

Alfalfa (Medicago sativa L.) is an important forage crop grown worldwide. Alfalfa is called the "queen of forage crops" due to its high forage yield and nutritional characteristics. The aim of this study was to undertake quantitative trait loci (QTL) mapping of yield and yield-related traits in an F₁ population of two alfalfa varieties that differ in their yield and yield-related traits.

RESULTS

We constructed a high-density linkage map using single nucleotide polymorphism (SNP) markers generated by restriction-site associated DNA sequencing (RAD-seq). The linkage map contains 4346 SNP and 119 simple sequence repeat (SSR) markers, with 32 linkage groups for each parent. The average marker distances were 3.00 and 1.32 cM, with coverages of 3455 cM and 4381 cM for paternal and maternal linkage maps, respectively. Using these maps and phenotypic data, we identified a total of 21 QTL for yield and yield components, including five for yield, five for plant height, five for branch number, and six for shoot diameter. Among them, six QTL were co-located for more than one trait. Five QTL explained more than 10% of the phenotypic variation.

CONCLUSIONS

We used RAD-seq to construct a linkage map for alfalfa that greatly enhanced marker density compared to previous studies. This high-density linkage map of alfalfa is a useful reference for mapping yield-related traits. Identified yield-related loci could be used to validate their usefulness in developing markers for maker-assisted selection in breeding populations to improve yield potential in alfalfa.

Genetic Loci Associated with Salt Tolerance in Advanced Breeding Populations of Tetraploid Alfalfa Using Genome-Wide Association Studies

Many agricultural lands in the western United States consist of soil with high concentrations of salt, which is detrimental to alfalfa ( L.) growth and production, especially in the region where water resource is limited. Developing alfalfa varieties with salt tolerance is imperative for sustainable production under increasing soil salinity. In the present study, we used advanced alfalfa breeding populations and evaluated five traits related to salt tolerance including biomass dry weight (DW) and fresh weight (FW), plant height (PH), leaf relative water content (RWC), and stomatal conductance (SC) under control and salt stress. Stress susceptibility index (SSI) of each trait and single-nucleotide polymorphism (SNP) markers generated by genotyping-by-sequencing (GBS) were used for genome-wide association studies (GWAS) to identify loci associated with salt tolerance. A total of 53 significant SNPs associated with salt tolerance were identified and they were located at 49 loci through eight chromosomes. A Basic Local Alignment Search Tool (BLAST) search of the regions surrounding the SNPs revealed 21 putative candidate genes associated with salt tolerance. The genetic architecture for traits related to salt tolerance characterized in this report could help in understanding the genetic mechanism by which salt stress affects plant growth and production in alfalfa. The markers and candidate genes identified in the present study would be useful for marker-assisted selection (MAS) in breeding salt-tolerant alfalfa after validation of the markers.

Harnessing the Potential of Forage Legumes, Alfalfa, Soybean, and Cowpea for Sustainable Agriculture and Global Food Security

Substantial improvements in access to food and increased purchasing power are driving many people toward consuming nutrition-rich foods causing an unprecedented demand for protein food worldwide, which is expected to rise further. Forage legumes form an important source of feed for livestock and have potential to provide a sustainable solution for food and protein security. Currently, alfalfa is a commercially grown source of forage and feed in many countries. However, soybean and cowpea also have the potential to provide quality forage and fodder for animal use. The cultivation of forage legumes is under threat from changing climatic conditions, indicating the need for breeding cultivars that can sustain and acclimatize to the negative effects of climate change. Recent progress in genetic and genomic tools have facilitated the identification of quantitative trait loci and genes/alleles that can aid in developing forage cultivars through genomics-assisted breeding. Furthermore, transgenic technology can be utilized to manipulate the genetic makeup of plants to improve forage digestibility for better animal performance. In this article, we assess the genetic potential of three important legume crops, alfalfa, soybean, and cowpea in supplying quality fodder and feed for livestock. In addition, we examine the impact of climate change on forage quality and discuss efforts made in enhancing the adaptation of the plant to the abiotic stress conditions. Subsequently, we suggest the application of integrative approaches to achieve adequate forage production amid the unpredictable climatic conditions.

Ectopic Expression of GsSRK in Medicago sativa Reveals Its Involvement in Plant Architecture and Salt Stress Responses

Receptor-like kinases (RLK) play fundamental roles in plant growth and stress responses. Compared with other RLKs, little information is provided concerning the S-locus LecRLK subfamily, which is characterized by an extracellular G-type lectin domain and an S-locus-glycop domain. Until now, the function of the G-type lectin domain is still unknown. In a previous research, we identified a Glycine soja S-locus LecRLK gene GsSRK, which conferred increased salt stress tolerance in transgenic Arabidopsis. In this study, to investigate the role of the G-type lectin domain and to breed transgenic alfalfa with superior salt stress tolerance, we transformed the full-length GsSRK (GsSRK-f) and a truncated version of GsSRK (GsSRK-t) deleting the G-type lectin domain into alfalfa. Our results showed that overexpression of GsSRK-t, but not GsSRK-f, resulted in changes of plant architecture, as evidenced by more branches but shorter shoots of GsSRK-t transgenic alfalfa, indicating a potential role of the extracellular G-type lectin domain in regulating plant architecture. Furthermore, we also found that transgenic alfalfa overexpressing either GsSRK-f or GsSRK-t showed increased salt stress tolerance, and GsSRK-t transgenic alfalfa displayed better growth (more branches and higher fresh weight) than GsSRK-f lines under salt stress. In addition, our results suggested that both GsSRK-f and GsSRK-t were involved in ion homeostasis, ROS scavenging, and osmotic regulation. Under salt stress, the Na+ content in the transgenic lines was significantly lower, while the K+ content was slightly higher than that in WT. Moreover, the transgenic lines displayed reduced ion leakage and MDA content, but increased SOD activity and proline content than WT. Notably, no obvious difference in these physiological indices was observed between GsSRK-f and GsSRK-t transgenic lines, implying that deletion of the GsSRK G-type lectin domain does not affect its physiological function in salt stress responses. In conclusion, results in this research reveal the dual role of GsSRK in regulating both plant architecture and salt stress responses.