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.

Salt and drought stresses are two primary abiotic stresses that inhibit growth and reduce the activity of photosynthetic apparatus in plants. Abscisic acid (ABA) plays a key role in abiotic stress regulation in plants. Some aldo–keto reductases (AKRs) can enhance various abiotic stresses resistance by scavenging cytotoxic aldehydes in some plants. However, there are few comprehensive reports of plant AKR genes and their expression patterns in response to abiotic stresses. In this study, we identified 30 putative AKR genes from Medicago truncatula. The gene characteristics, coding protein motifs, and expression patterns of these MtAKRs were analyzed to explore and identify candidate genes in regulation of salt, drought, and ABA stresses. The phylogenetic analysis result indicated that the 52 AKRs in Medicago truncatula and Arabidopsis thaliana can be divided into three groups and six subgroups. Fifteen AKR genes in M. truncatula were randomly selected from each group or subgroup, to investigate their response to salt (200 mM of NaCl), drought (50 g·L⁻¹ of PEG 6000), and ABA (100 µM) stresses in both leaves and roots. The results suggest that MtAKR1, MtAKR5, MtAKR11, MtAKR14, MtAKR20, and MtAKR29 may play important roles in response to these stresses.

Bioactive gibberellic acids (GAs) are diterpenoid plant hormones that are biosynthesized through complex pathways and control various aspects of growth and development. Although GA biosynthesis has been intensively studied, the downstream metabolic pathways regulated by GAs have remained largely unexplored. We investigated Tnt1 retrotransposon insertion mutant lines of Medicago truncatula with a dwarf phenotype by forward and reverse genetics screening and phylogenetic, molecular, biochemical, proteomic and metabolomic analyses. Three Tnt1 retrotransposon insertion mutant lines of the gibberellin 3-beta-dioxygenase 1 gene (GA3ox1) with a dwarf phenotype were identified, in which the synthesis of GAs (GA₃ and GA₄) was inhibited. Phenotypic analysis revealed that plant height, root and petiole length of ga3ox1 mutants were shorter than those of the wild type (Medicago truncatula ecotype R108). Leaf size was also much smaller in ga3ox1 mutants than that in wild-type R108, which is probably due to cell-size diminution instead of a decrease in cell number. Proteomic and metabolomic analyses of ga3ox1/R108 leaves revealed that in the ga3ox1 mutant, flavonoid isoflavonoid biosynthesis was significantly up-regulated, while nitrogen metabolism was down-regulated. Additionally, we further demonstrated that flavonoid and isoflavonoid biosynthesis was induced by prohexadione calcium, an inhibitor of GA3ox enzyme, and inhibited by exogenous GA₃. In contrast, nitrogen metabolism was promoted by exogenous GA₃ but inhibited by prohexadione calcium. The results of this study further demonstrated that GAs play critical roles in positively regulating nitrogen metabolism and transport and negatively regulating flavonoid biosynthesis through GA-mediated signaling pathways in leaves.

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.

Appropriate root system architecture (RSA) can improve alfalfa yield, yet its genetic basis remains largely unexplored. Here, we evaluated six RSA traits in 171 alfalfa genotypes under controlled greenhouse conditions. We also analyzed five yield-related traits in normal and drought stress environments, finding a significant correlation (0.50) between root dry weight (RDW) and alfalfa dry weight under normal conditions (N_DW). A genome-wide association study (GWAS) was performed using 1,303,374 single nucleotide polymorphisms (SNPs) to explore the relationships between RSA traits. Sixty significant SNPs (-log₁₀(P) ≥ 5) were identified, with genes within 50 kb upstream and downstream ranges primarily enriched in GO terms related to root development, hormone synthesis and signaling, as well as morphological development. Further analysis identified 19 high-confidence candidate genes, including AUXIN RESPONSE FACTORs (ARFs), LATERAL ORGAN BOUNDARIES-DOMAIN (LBD), and WUSCHEL-RELATED HOMEOBOX (WOX). We verified that the forage dry weight under both normal and drought conditions exhibited significant differences among materials with different numbers of favourable haplotypes. Alfalfa containing more favourable haplotypes exhibited higher forage yields, whereas favourable haplotypes were not subjected to human selection during alfalfa breeding. Genomic prediction (GP) utilized SNPs from GWAS and machine learning for each RSA trait, achieving prediction accuracies ranging from 0.70 for secondary root position (SRP) to 0.80 for root length (RL), indicating robust predictive capability across the assessed traits. These findings provide new insights into the genetic underpinnings of root development in alfalfa, potentially informing future breeding strategies aimed at improving yield.