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.

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.

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.

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 Analysis of the LATERAL ORGAN BOUNDARIES Domain (LBD) Members in Alfalfa and the Involvement of MsLBD48 in Nitrogen Assimilation.

The LATERAL ORGAN BOUNDARIES DOMAIN (LBD) proteins, a transcription factor family specific to the land plants, have been implicated in multiple biological processes including organ development, pathogen response and the uptake of inorganic nitrogen. The study focused on LBDs in legume forage Alfalfa. The genome-wide analysis revealed that in Alfalfa 178 loci across 31 allelic chromosomes encoded 48 unique LBDs (MsLBDs), and the genome of its diploid progenitor M. sativa spp. Caerulea encoded 46 LBDs. Synteny analysis indicated that the expansion of AlfalfaLBDs was attributed to the whole genome duplication event. The MsLBDs were divided into two major phylogenetic classes, and the LOB domain of the Class I members was highly conserved relative to that of the Class II. The transcriptomic data demonstrated that 87.5% of MsLBDs were expressed in at least one of the six test tissues, and Class II members were preferentially expressed in nodules. Moreover, the expression of Class II LBDs in roots was upregulated by the treatment of inorganic nitrogen such as KNO3 and NH4Cl (0.3 mM). The overexpression of MsLBD48, a Class II member, in Arabidopsis resulted in growth retardance with significantly declined biomass compared with the non-transgenic plants, and the transcription level of the genes involved in nitrogen uptake or assimilation, including NRT1.1, NRT2.1, NIA1 and NIA2 was repressed. Therefore, the LBDs in Alfalfa are highly conserved with their orthologs in embryophytes. Our observations that ectopic expression of MsLBD48 inhibited Arabidopsis growth by repressing nitrogen adaption suggest the negative role of the transcription factor in plant uptake of inorganic nitrogen. The findings imply the potential application of MsLBD48 in Alfalfa yield improvement via gene editing.

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.

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.

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.

Combining QTL mapping and RNA-Seq Unravels candidate genes for Alfalfa (Medicago sativa L.) leaf development.

BACKGROUND: Leaf size affects crop canopy morphology and photosynthetic efficiency, which can influence forage yield and quality. It is of great significance to mine the key genes controlling leaf development for breeding new alfalfa varieties. In this study, we mapped leaf length (LL), leaf width (LW), and leaf area (LA) in an F1 mapping population derived from a cultivar named ZhongmuNo.1 with larger leaf area and a landrace named Cangzhou with smaller leaf area. RESULTS: This study showed that the larger LW was more conducive to increasing LA. A total of 24 significant quantitative trait loci (QTL) associated with leaf size were identified on both the paternal and maternal linkage maps. Among them, nine QTL explained about 11.50-22.45% phenotypic variation. RNA-seq analysis identified 2,443 leaf-specific genes and 3,770 differentially expressed genes. Combining QTL mapping, RNA-seq alalysis, and qRT-PCR, we identified seven candidate genes associated with leaf development in five major QTL regions. CONCLUSION: Our study will provide a theoretical basis for marker-assisted breeding and lay a foundation for further revealing molecular mechanism of leaf development in alfalfa.

Metabolomic changes in crown of alfalfa (Medicago sativa L.) during de-acclimation.

Alfalfa is a high-quality forage legume species that is widely cultivated at high latitudes worldwide. However, a decrease in cold tolerance in early spring seriously affects regrowth and persistence of alfalfa. There has been limited research on the metabolomic changes that occur during de-acclimation. In this study, a liquid chromatography-mass spectrometry system was used to compare the metabolites in two alfalfa cultivars during a simulated overwintering treatment. In four pairwise comparisons, 367 differential metabolites were identified, of which 31 were annotated according to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Many of these metabolites were peptides, carbohydrates, and lipids. At the subclass level, 17 major pathways were revealed to be significantly enriched (P < 0.05). The main differential metabolites included amino acids, peptides and analogs, carbohydrates, and glycerol phosphocholines. A metabolomic analysis showed that the up-regulation of unsaturated fatty acids and amino acids as well as the enhancement of the related metabolic pathways might be an effective strategy for increasing alfalfa cold tolerance. Furthermore, glycerophospholipid metabolism affects alfalfa cold tolerance in early spring. Study results provide new insights about the changes in alfalfa metabolites that occur during de-acclimation, with potential implications for the selection and breeding of cold-tolerant cultivars.

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.

MtPT5 phosphate transporter is involved in leaf growth and phosphate accumulation of Medicago truncatula.

Phosphorus (P) is an indispensable mineral nutrient for plant growth and agricultural production. Plants acquire and redistribute inorganic phosphate (Pi) via Pi transporters (PHT1s/PTs). However, apart from MtPT4, functions of the M. truncatula (Medicago truncatula) PHT1s remain unclear. In this study, we evaluated the function of the PHT1 family transporter MtPT5 in M. truncatula. MtPT5 was closely related to AtPHT1; 1 in Arabidopsis (Arabidopsis thaliana) and GmPT7 in soybean (Glycine max). MtPT5 was highly expressed in leaves in addition to roots and nodules. Ectopic expression of MtPT5 complemented the Pi-uptake deficiency of Arabidopsis pht1;1Δ4Δ double mutant, demonstrating the Pi-transport activity of MtPT5 in plants. When overexpressing MtPT5 in M. truncatula, the transgenic plants showed larger leaves, accompanying with higher biomass and Pi enrichment compared with wild type. All these data demonstrate that MtPT5 is important for leaf growth and Pi accumulation of M. truncatula and provides a target for molecular breeding to improve forage productivity.