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.

Overexpression of MsDREB1C Modulates Growth and Improves Forage Quality in Tetraploid Alfalfa (Medicago sativa L.).

DREB has been reported to be involved in plant growth and response to environmental factors. However, the function of DREB in growth and development has not been elucidated in alfalfa (Medicago sativa L.), a perennial tetraploid forage cultivated worldwide. In this study, an ortholog of MtDREB1C was characterized from alfalfa and named MsDREB1C accordingly. MsDREB1C was significantly induced by abiotic stress. The transcription factor MsDREB1C resided in the nucleus and had self-transactivation activity. The MsDREB1C overexpression (OE) alfalfa displayed growth retardation under both long-day and short-day conditions, which was supported by decreased MsGA20ox and upregulated MsGA2ox in the OE lines. Consistently, a decrease in active gibberellin (GA) was detected, suggesting a negative effect of MsDREB1C on GA accumulation in alfalfa. Interestingly, the forage quality of the OE lines was better than that of WT lines, with higher crude protein and lower lignin content, which was supported by an increase in the leaf-stem ratio (LSR) and repression of several lignin-synthesis genes (MsNST, MsPAL1, MsC4H, and Ms4CL). Therefore, this study revealed the effects of MsDREB1C overexpression on growth and forage quality via modifying GA accumulation and lignin synthesis, respectively. Our findings provide a valuable candidate for improving the critical agronomic traits of alfalfa, such as overwintering and feeding value of the forage.

Evolutionary genomics of climatic adaptation and resilience to climate change in alfalfa.

Given the escalating impact of climate change on agriculture and food security, gaining insights into the evolutionary dynamics of climatic adaptation and uncovering climate-adapted variation can empower the breeding of climate-resilient crops to face future climate change. Alfalfa (Medicago sativa subsp. sativa), the queen of forages, shows remarkable adaptability across diverse global environments, making it an excellent model for investigating species responses to climate change. In this study, we performed population genomic analyses using genome resequencing data from 702 accessions of 24 Medicago species to unravel alfalfa's climatic adaptation and genetic susceptibility to future climate change. We found that interspecific genetic exchange has contributed to the gene pool of alfalfa, particularly enriching defense and stress-response genes. Intersubspecific introgression between M. sativa subsp. falcata (subsp. falcata) and alfalfa not only aids alfalfa's climatic adaptation but also introduces genetic burden. A total of 1671 genes were associated with climatic adaptation, and 5.7% of them were introgressions from subsp. falcata. By integrating climate-associated variants and climate data, we identified populations that are vulnerable to future climate change, particularly in higher latitudes of the Northern Hemisphere. These findings serve as a clarion call for targeted conservation initiatives and breeding efforts. We also identified pre-adaptive populations that demonstrate heightened resilience to climate fluctuations, illuminating a pathway for future breeding strategies. Collectively, this study enhances our understanding about the local adaptation mechanisms of alfalfa and facilitates the breeding of climate-resilient alfalfa cultivars, contributing to effective agricultural strategies for facing future climate change.

Functional Characterization of the MsFKF1 Gene Reveals Its Dual Role in Regulating the Flowering Time and Plant Height in Medicago sativa L.

Alfalfa (M. sativa), a perennial legume forage, is known for its high yield and good quality. As a long-day plant, it is sensitive to changes in the day length, which affects the flowering time and plant growth, and limits alfalfa yield. Photoperiod-mediated delayed flowering in alfalfa helps to extend the vegetative growth period and increase the yield. We isolated a blue-light phytohormone gene from the alfalfa genome that is an ortholog of soybean FKF1 and named it MsFKF1. Gene expression analyses showed that MsFKF1 responds to blue light and the circadian clock in alfalfa. We found that MsFKF1 regulates the flowering time through the plant circadian clock pathway by inhibiting the transcription of E1 and COL, thus suppressing FLOWERING LOCUS T a1 (FTa1) transcription. In addition, transgenic lines exhibited higher plant height and accumulated more biomass in comparison to wild-type plants. However, the increased fiber (NDF and ADF) and lignin content also led to a reduction in the digestibility of the forage. The key genes related to GA biosynthesis, GA20OX1, increased in the transgenic lines, while GA2OX1 decreased for the inactive GA transformation. These findings offer novel insights on the function of MsFKF1 in the regulation of the flowering time and plant height in cultivated M. sativa. These insights into MsFKF1's roles in alfalfa offer potential strategies for molecular breeding aimed at optimizing flowering time and biomass yield.

The GARP family transcription factor MtHHO3 negatively regulates salt tolerance in Medicago truncatula.

High salinity is one of the detrimental environmental factors restricting plant growth and crop production throughout the world. This study demonstrated that the GARP family transcription factor MtHHO3 is involved in response to salt stress and abscisic acid (ABA) signaling in Medicago truncatula. The transcription of MtHHO3 was repressed by salt, osmotic stress, and ABA treatment. The seed germination assay showed that, overexpression of MtHHO3 in Arabidopsis thaliana caused hypersensitivity to salt and osmotic stress, but increased resistance to ABA inhibition. Overexpression of MtHHO3 in M. truncatula resulted in decreased tolerance of salinity, while loss-of-function mutants mthho3-1 and mthho3-2 were more resistant to salt stress compared with wild-type plants. qRT-PCR analyses showed that MtHHO3 downregulated the expression of genes in stress and ABA responsive pathways. We further demonstrated that MtHHO3 repressed the transcription of the pathogenesis-related gene MtPR2 by binding to its promoter. Overall, these results indicate that MtHHO3 negatively regulates salt stress response in plants and deepen our understanding of the role of the GARP subfamily transcription factors in modulating salt stress and ABA signaling.

A genome-wide study of the lipoxygenase gene families in Medicago truncatula and Medicago sativa reveals that MtLOX24 participates in the methyl jasmonate response.

BACKGROUND: Lipoxygenase (LOX) is a multifunctional enzyme that is primarily related to plant organ growth and development, biotic and abiotic stress responses, and production of flavor-associated metabolites. In higher plants, the LOX family encompasses several isozymes with varying expression patterns between tissues and developmental stages. These affect processes including seed germination, seed storage, seedling growth, fruit ripening, and leaf senescence. LOX family genes have multiple functions in response to hormones such as methyl jasmonate (MeJA) and salicylic acid. RESULTS: In this study, we identified 30 and 95 LOX homologs in Medicago truncatula and Medicago sativa, respectively. These genes were characterized with analyses of their basic physical and chemical properties, structures, chromosomal distributions, and phylogenetic relationships to understand structural variations and their physical locations. Phylogenetic analysis was conducted for members of the three LOX subfamilies (9-LOX, type I 13-LOX, and type II 13-LOX) in Arabidopsis thaliana, Glycine max, M. truncatula, and M. sativa. Analysis of predicted promoter elements revealed several relevant cis-acting elements in MtLOX and MsLOX genes, including abscisic acid (ABA) response elements (ABREs), MeJA response elements (CGTCA-motifs), and antioxidant response elements (AREs). Cis-element data combined with transcriptomic data demonstrated that LOX gene family members in these species were most likely related to abiotic stress responses, hormone responses, and plant development. Gene expression patterns were confirmed via quantitative reverse transcription PCR. Several MtLOX genes (namely MtLOX15, MtLOX16, MtLOX20, and MtLOX24) belonging to the type I 13-LOX subfamily and other LOX genes (MtLOX7, MtLOX11, MsLOX23, MsLOX87, MsLOX90, and MsLOX94) showed significantly different expression levels in the flower tissue, suggesting roles in reproductive growth. Type I 13-LOXs (MtLOX16, MtLOX20, MtLOX21, MtLOX24, MsLOX57, MsLOX84, MsLOX85, and MsLOX94) and type II 13-LOXs (MtLOX5, MtLOX6, MtLOX9, MtLOX10, MsLOX18, MsLOX23, and MsLOX30) were MeJA-inducible and were predicted to function in the jasmonic acid signaling pathway. Furthermore, exogenous MtLOX24 expression in Arabidopsis verified that MtLOX24 was involved in MeJA responses, which may be related to insect-induced abiotic stress. CONCLUSIONS: We identified six and four LOX genes specifically expressed in the flowers of M. truncatula and M. sativa, respectively. Eight and seven LOX genes were induced by MeJA in M. truncatula and M. sativa, and the LOX genes identified were mainly distributed in the type I and type II 13-LOX subfamilies. MtLOX24 was up-regulated at 8 h after MeJA induction, and exogenous expression in Arabidopsis demonstrated that MtLOX24 promoted resistance to MeJA-induced stress. This study provides valuable new information regarding the evolutionary history and functions of LOX genes in the genus Medicago.

FtsH proteases confer protection against salt and oxidative stress in Medicago sativa L.

Plant filamentation temperature-sensitive H (FtsH) proteins are ATP-dependent zinc proteases that play an important role in regulating abiotic stress adaptions. Here we explore their potential role in abiotic stress tolerance in alfalfa, an important legume crop. Genomic analysis revealed seventeen MsFtsH genes in five clusters, which generally featured conserved domains and gene structures. Furthermore, the expression of MsFtsHs was found to be tightly associated with abiotic stresses, including osmotic, salt and oxidative stress. In addition, numerous stress responsive cis-elements, including those related to ABA, auxin, and salicylic acid, were identified in their promoter regions. Moreover, MsFtsH8 overexpression was shown to confer tolerance to salt and oxidative stress which was associated with reduced levels of reactive oxygen species, and enhanced expression and activity of antioxidant enzymes. Our results highlight MsFtsHs as key factors in abiotic stress tolerance, and show their potential usefulness for breeding alfalfa and other crops with improved yield and stress tolerance.

Genome-wide identification and characterization of filamentation temperature-sensitive H (FtsH) genes and expression analysis in response to multiple stresses in Medicago truncatula.

BACKGROUND: Filamentation temperature-sensitive H (FtsH) is an AAA+ ATP-dependent protease that plays a vital role in plant environmental adaption and tolerance. However, little is known about the function of the FtsH gene family in the most important legume model plant, Medicago truncatula. METHODS AND RESULTS: To identify and investigate the potential stress adaptation roles of FtsH gene family in M. truncatula, we conducted a series of genome-wide characterization and expression analyses. Totally, twenty MtFtsH genes were identified, which were unevenly distributed across eight chromosomes and classified into six evolution groups based on their phylogenetic relationships, with each group containing similar structures and motifs. Furthermore, MtFtsH genes exhibited a high degree of collinearity and homology with leguminous plants such as alfalfa and soybean. Multiple cis-elements in the upstream region of MtFtsH genes were also identified that responded to light, abiotic stress, and phytohormones. Public RNA-seq data indicated that MtFtsH genes were induced under both salt and drought stresses, and our transcript expression analysis showed that MtFtsH genes of MtFtsH1, MtFtsH2, MtFtsH4, MtFtsH9, and MtFtsH10 were up-regulated after ABA, H2O2, PEG, and NaCl treatments. These results suggest that MtFtsH genes may play a critical role in drought and high salt stress responses and the adaption processes of plants. CONCLUSIONS: This study provides a systematic analysis of FtsH gene family in M. truncatula, serving as a valuable molecular theoretical basis for future functional investigations. Our findings also extend the pool of potential candidate genes for the genetic improvement of abiotic stress tolerance in legume crops.

Characterization of the Heat Shock Transcription Factor Family in Medicago sativa L. and Its Potential Roles in Response to Abiotic Stresses.

Heat shock transcription factors (HSFs) are important regulatory factors in plant stress responses to various biotic and abiotic stresses and play important roles in growth and development. The HSF gene family has been systematically identified and analyzed in many plants but it is not in the tetraploid alfalfa genome. We detected 104 HSF genes (MsHSFs) in the tetraploid alfalfa genome ("Xinjiangdaye" reference genome) and classified them into three subgroups: 68 in HSFA, 35 in HSFB and 1 in HSFC subgroups. Basic bioinformatics analysis, including genome location, protein sequence length, protein molecular weight and conserved motif identification, was conducted. Gene expression analysis revealed tissue-specific expression for 13 MsHSFs and tissue-wide expression for 28 MsHSFs. Based on transcriptomic data analysis, 21, 11 and 27 MsHSFs responded to drought stress, cold stress and salt stress, respectively, with seven responding to all three. According to RT-PCR, MsHSF27/33 expression gradually increased with cold, salt and drought stress condition duration; MsHSF6 expression increased over time under salt and drought stress conditions but decreased under cold stress. Our results provide key information for further functional analysis of MsHSFs and for genetic improvement of stress resistance in alfalfa.

Transcriptome and GWAS Analyses Reveal Candidate Gene for Root Traits of Alfalfa during Germination under Salt Stress.

Alfalfa growth and production in China are negatively impacted by high salt concentrations in soils, especially in regions with limited water supplies. Few reliable genetic markers are currently available for salt tolerance selection. As a result, molecular breeding strategies targeting alfalfa are hindered. Therefore, with the continuous increase in soil salinity in agricultural lands, it is indispensable that a salt-tolerant variety of alfalfa is produced. We collected 220 alfalfa varieties around the world for resequencing and performed genome-wide association studies (GWASs). Alfalfa seeds were germinated in saline water with different concentrations of NaCl, and the phenotypic differences in several key root traits were recorded. In the phenotypic analysis, the breeding status and geographical origin strongly affected the salt tolerance of alfalfa. Forty-nine markers were significantly associated with salt tolerance, and 103 candidate genes were identified based on linkage disequilibrium. A total of 2712 differentially expressed genes were upregulated and 3570 were downregulated based on transcriptomic analyses. Some candidate genes that affected root development in the seed germination stage were identified through the combination of GWASs and transcriptome analyses. These genes could be used for molecular breeding strategies to increase alfalfa's salt tolerance and for further research on salt tolerance in general.

Overexpression of the elongation factor MtEF1A1 promotes salt stress tolerance in Arabidopsis thaliana and Medicago truncatula.

BACKGROUND: Elongation factor 1 A (EF1A), an essential regulator for protein synthesis, has been reported to participate in abiotic stress responses and environmental adaption in plants. However, the role of EF1A in abiotic stress response was barely studied in Medicago truncatula. Here, we identified elongation factor (EF) genes of M. truncatula and studied the salt stress response function of MtEF1A1 (MTR_6g021805). RESULTS: A total of 34 EF genes were identified in the M. truncatula genome. Protein domains and motifs of EFs were highly conserved in plants. MtEF1A1 has the highest expression levels in root nodules and roots, followed by the leaves and stems. Transgenic Arabidopsis thaliana overexpressing MtEF1A1 was more resistant to salt stress treatment, with higher germination rate, longer roots, and more lateral roots than wild type plant. In addition, lower levels of H2O2 and malondialdehyde (MDA) were also detected in transgenic Arabidopsis. Similarly, MtEF1A1 overexpressing M. truncatula was more resistant to salt stress and had lower levels of reactive oxygen species (ROS) in leaves. Furthermore, the expression levels of abiotic stress-responsive genes (MtRD22A and MtCOR15A) and calcium-binding genes (MtCaM and MtCBL4) were upregulated in MtEF1A1 overexpressing lines of M. truncatula. CONCLUSION: These results suggested that MtEF1A1 play a positive role in salt stress regulation. MtEF1A1 may realize its function by binding to calmodulin (CaM) or by participating in Ca2+-dependent signaling pathway. This study revealed that MtEF1A1 is an important regulator for salt stress response in M. truncatula, and provided potential strategy for salt-tolerant plant breeding.

Genome-Wide Identification and Phylogenetic and Expression Analyses of the PLATZ Gene Family in Medicago sativa L.

The PLATZ family is a novel class of plant-specific zinc finger transcription factors with important roles in plant growth and development and abiotic stress responses. PLATZ members have been identified in many plants, including Oryza sativa, Zea mays, Triticum aestivum, Fagopyrum tataricum, and Arabidopsis thaliana; however, due to the complexity of the alfalfa reference genome, the members of the PLATZ gene family in alfalfa (Medicago sativa L.) have not been systematically identified and analyzed. In this study, 55 Medicago sativa PLATZ genes (MsPLATZs) were identified in the alfalfa "Xinjiangdaye" reference genome. Basic bioinformatic analysis was performed, including the characterization of sequence lengths, protein molecular weights, genomic positions, and conserved motifs. Expression analysis reveals that 7 MsPLATZs are tissue-specifically expressed, and 10 MsPLATZs are expressed in all examined tissues. The transcriptomic expression of these genes is obvious, indicating that these MsPLATZs have different functions in the growth and development of alfalfa. Based on transcriptome data analysis and real-time quantitative PCR (RT-qPCR), we identified 22, 22, and 21 MsPLATZ genes that responded to salt, cold, and drought stress, respectively, with 20 MsPLATZs responding to all three stresses. This study lays a foundation for further exploring the functions of MsPLATZs, and provides ideas for the improvement of alfalfa varieties and germplasm innovation.

Application of machine learning to explore the genomic prediction accuracy of fall dormancy in autotetraploid alfalfa.

Fall dormancy (FD) is an essential trait to overcome winter damage and for alfalfa (Medicago sativa) cultivar selection. The plant regrowth height after autumn clipping is an indirect way to evaluate FD. Transcriptomics, proteomics, and quantitative trait locus mapping have revealed crucial genes correlated with FD; however, these genes cannot predict alfalfa FD very well. Here, we conducted genomic prediction of FD using whole-genome SNP markers based on machine learning-related methods, including support vector machine (SVM) regression, and regularization-related methods, such as Lasso and ridge regression. The results showed that using SVM regression with linear kernel and the top 3000 genome-wide association study (GWAS)-associated markers achieved the highest prediction accuracy for FD of 64.1%. For plant regrowth height, the prediction accuracy was 59.0% using the 3000 GWAS-associated markers and the SVM linear model. This was better than the results using whole-genome markers (25.0%). Therefore, the method we explored for alfalfa FD prediction outperformed the other models, such as Lasso and ElasticNet. The study suggests the feasibility of using machine learning to predict FD with GWAS-associated markers, and the GWAS-associated markers combined with machine learning would benefit FD-related traits as well. Application of the methodology may provide potential targets for FD selection, which would accelerate genetic research and molecular breeding of alfalfa with optimized FD.

Genome-wide identification, characterization, and expression analysis of UDP-glycosyltransferase genes associated with secondary metabolism in alfalfa (Medicago sativa L.).

Uridine diphosphate glycosyltransferases (UGTs) are enzymes that catalyze glycosylation modifications and play an essential role in regulating plant metabolism. Alfalfa (Medicago sativa L.) is the most important legume in the world due to its high yields and protein content; however, the UGT genes in alfalfa have not yet been studied. Identifying UGT genes with metabolic roles in alfalfa is essential for identifying and modifying genetic traits that are relevant to yield and quality. In this study, 90 of the 239 UGT genes identified from the alfalfa "Zhongmu No. 1" genome database were found to be related to secondary metabolism, and a series of gene family characterization analyses were conducted on each. The results demonstrated that all 90 UGT genes were unevenly distributed on eight chromosomes with few introns and that tandem duplications were the crucial driving force expanding the UGT family in alfalfa. Notably, the 90 UGT genes can be clustered into ten evolutionary groups which contain specific PSPG motifs, and genes in these ten groups have specific tissue expressions. This suggests that the UGT genes in each group could have similar glycosylation roles corresponding to analogous secondary metabolites in alfalfa. Additionally, multiple cis-acting elements found in MsUGT promoter regions, such as phytohormone and flavonoids, indicate that 90 UGT members could be induced by these features, which are also related to secondary metabolism. Therefore, our study identified 90 UGT members inten evolutionary groups that are likely related to glycosylation modifications with secondary metabolites in alfalfa. These findings help uncover pivotal regulatory mechanisms associated with secondary metabolism in plant yield and quality and contribute to genetic modification and breeding in alfalfa and other plant species.

Genome-Wide Identification and Expression Analysis of the NAC Gene Family in Alfalfa Revealed Its Potential Roles in Response to Multiple Abiotic Stresses.

NAC (NAM, ATAF1/2, and CUC2) transcription factors compose one of the largest families of plant-specific transcription factors; they are widely involved in plant growth and development and have especially important roles in improving stress resistance in plants. However, NAC gene family members in alfalfa (Medicago sativa L.) have not been systematically identified and analyzed genome-wide due to the complexity of the alfalfa reference genome. In this study, a total of 421 M. sativa NAC genes (MsNACs) were identified from the alfalfa "Xinjiangdaye" reference genome. Basic bioinformatics analysis, including characterization of sequence length, protein molecular weight and genome position and conserved motif analysis, was conducted. Expression analysis showed that 47 MsNACs had tissue-specific expression, and 64 MsNACs were expressed in all tissues. The transcriptomic profiles of the genes were very different, indicating that these MsNACs have various functions in alfalfa growth and development. We identified 25, 42 and 47 MsNACs that respond to cold, drought and salt stress based on transcriptome data analysis and real-time quantitative PCR (RT−qPCR). Furthermore, 22 MsNACs were found to respond to both salt and drought stress, and 15 MsNACs were found to respond to cold, salt and drought stress. The results of this study could provide valuable information for further functional analysis of MsNACs and for the improvement of stress resistance in alfalfa.

A Genome-Wide Association Study Coupled With a Transcriptomic Analysis Reveals the Genetic Loci and Candidate Genes Governing the Flowering Time in Alfalfa (Medicago sativa L.).

The transition to flowering at the right time is very important for adapting to local conditions and maximizing alfalfa yield. However, the understanding of the genetic basis of the alfalfa flowering time remains limited. There are few reliable genes or markers for selection, which hinders progress in genetic research and molecular breeding of this trait in alfalfa. We sequenced 220 alfalfa cultivars and conducted a genome-wide association study (GWAS) involving 875,023 single-nucleotide polymorphisms (SNPs). The phenotypic analysis showed that the breeding status and geographical origin strongly influenced the alfalfa flowering time. Our GWAS revealed 63 loci significantly related to the flowering time. Ninety-five candidate genes were detected at these SNP loci within 40 kb (20 kb up- and downstream). Thirty-six percent of the candidate genes are involved in development and pollen tube growth, indicating that these genes are key genetic mechanisms of alfalfa growth and development. The transcriptomic analysis showed that 1,924, 2,405, and 3,779 differentially expressed genes (DEGs) were upregulated across the three growth stages, while 1,651, 2,613, and 4,730 DEGs were downregulated across the stages. Combining the results of our GWAS and transcriptome analysis, in total, 38 candidate genes (7 differentially expressed during the bud stage, 13 differentially expressed during the initial flowering stage, and 18 differentially expressed during the full flowering stage) were identified. Two SNPs located in the upstream region of the Msa0888690 gene (which is involved in isop renoids) were significantly related to flowering. The two significant SNPs within the upstream region of Msa0888690 existed as four different haplotypes in this panel. The genes identified in this study represent a series of candidate targets for further research investigating the alfalfa flowering time and could be used for alfalfa molecular breeding.

Genome-wide identification, phylogeny and expression analysis of the SPL gene family and its important role in salt stress in Medicago sativa L.

BACKGROUND: SQUAMOSA promoter-binding protein-like (SPL) transcription factors are widely present in plants and are involved in signal transduction, the stress response and development. The SPL gene family has been characterized in several model species, such as A. thaliana and G. max. However, there is no in-depth analysis of the SPL gene family in forage, especially alfalfa (Medicago sativa L.), one of the most important forage crops worldwide. RESULT: In total, 76 putative MsSPL genes were identified in the alfalfa genome with an uneven distribution. Based on their identity and gene structure, these MsSPLs were divided into eight phylogenetic groups. Seventy-three MsSPL gene pairs arose from segmental duplication events, and the MsSPLs on the four subgenomes of individual chromosomes displayed high collinearity with the corresponding M. truncatula genome. The prediction of the cis-elements in the promoter regions of the MsSPLs detected two copies of ABA (abscisic acid)-responsive elements (ABREs) on average, implying their potential involvement in alfalfa adaptation to adverse environments. The transcriptome sequencing of MsSPLs in roots and leaves revealed that 54 MsSPLs were expressed in both tissues. Upon salt treatment, three MsSPLs (MsSPL17, MsSPL23 and MsSPL36) were significantly regulated, and the transcription level of MsSPL36 in leaves was repressed to 46.6% of the control level. CONCLUSION: In this study, based on sequence homology, we identified 76 SPL genes in the alfalfa. The SPLs with high identity shared similar gene structures and motifs. In total, 71.1% (54 of 76) of the MsSPLs were expressed in both roots and leaves, and the majority (74.1%) preferred underground tissues to aerial tissues. MsSPL36 in leaves was significantly repressed under salt stress. These findings provide comprehensive information regarding the SPB-box gene family for improve alfalfa tolerance to high salinity.

RAD-Seq-Based High-Density Linkage Maps Construction and Quantitative Trait Loci Mapping of Flowering Time Trait in Alfalfa (Medicago sativa L.).

Alfalfa (Medicago sativa L.) is a perennial forage crop known as the "Queen of Forages." To dissect the genetic mechanism of flowering time (FT) in alfalfa, high-density linkage maps were constructed for both parents of an F1 mapping population derived from a cross between Cangzhou (P1) and ZhongmuNO.1 (P2), consisting of 150 progenies. The FT showed a transgressive segregation pattern in the mapping population. A total of 13,773 single-nucleotide polymorphism markers was obtained by using restriction-site associated DNA sequencing and distributed on 64 linkage groups, with a total length of 3,780.49 and 4,113.45 cM and an average marker interval of 0.58 and 0.59 cM for P1 and P2 parent, respectively. Quantitative trait loci (QTL) analyses were performed using the least square means of each year as well as the best linear unbiased prediction values across 4 years. Sixteen QTLs for FT were detected for P1 and 22 QTLs for P2, accounting for 1.40-16.04% of FT variation. RNA-Seq analysis at three flowering stages identified 5,039, 7,058, and 7,996 genes that were differentially expressed between two parents, respectively. Based on QTL mapping, DEGs analysis, and functional annotation, seven candidate genes associated with flowering time were finally detected. This study discovered QTLs and candidate genes for alfalfa FT, making it a useful resource for breeding studies on this essential crop.

Identification and characterization of stress responsive homeodomain leucine zipper transcription factors in Medicago truncatula.

BACKGROUND: Homeodomain leucine zipper (HD-ZIP) transcription factors play roles in regulating plant development and responses to abiotic stresses; however, how HD-ZIP genes in Medicago truncatula are involved in abiotic stress response remains elusive. METHODS AND RESULTS: The HD-ZIP I genes in Medicago truncatula were identified and characterized, and their expression patterns in different tissues and under different abiotic stresses were analyzed. A total of 15 Medicago truncatula HD-ZIP I genes were identified and a phylogenetic analysis of HD-ZIP I proteins in Arabidopsis thaliana and Medicago truncatula was conducted. Fifteen HD-ZIP I genes showed diverse tissue preferences. Among them, expressions of MtHB22 and MtHB51 were specially detected in vegetative buds. In addition, they responded to various abiotic stresses, including salinity and osmotic stress and abscisic acid (ABA). For instance, MtHB7 and MtHB12 expression levels were found to be positively associated with salt, osmotic stress and ABA in both shoots and roots, while MtHB13 and MtHB23 were negatively associated with these stresses in Medicago truncatula. CONCLUSION: The HD-ZIP I genes in Medicago truncatula are evolutionarily conserved, but also exhibit gene duplication and gene loss events. Differential expression analysis of Medicago truncatula HD-ZIP I genes indicated their crucial roles in abiotic stress responses. Our genome-wide analysis of the HD-ZIP I transcription factor family in Medicago truncatula provided a valuable reference for further research.

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.