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.

Saline-alkaline stress is a universal abiotic stress that adversely affects plant growth and productivity. Saline-alkaline conditions results in plant abnormal transcriptome expression finally manifesting as defective phenotypes. Considerable research has revealed the active role of microRNA in various stress conditions. This study was aimed to identify novel miRNAs and the miRNA expression patterns in the leguminous model plant R108 (Medicago truncatula). The miRNA contained in the total RNA extracted from Medicago truncatula seedlings (72 h) that had been treated with solutions mimicking saline and alkaline soils was subjected to miRNA deep sequencing. The Illumina HiSeq sequencing platform was used to analyze nine small RNA libraries of three treatment groups: distilled water, 20 mM NaCl + Na₂SO₄ and 5 mM Na₂CO₃ + NaHCO₃. Sequencing revealed that 876 miRNAs including 664 known miRNAs and 212 potential novel miRNAs were present in all the libraries. The miR159 family, miR156 family, miR2086-3p, miR396, miR166, miR319, miR167, miR5213-5p, miR1510 and miR2643 were among the most expressed miRNAs in all libraries. The results of miRNAs expression under treatments were validated by reverse-transcription quantitative PCR (RT-qPCR). Target gene prediction through computational analysis and pathway annotation analysis revealed that the primary pathways affected by stress were related to plant development, including metabolic processes, single-organism processes and response to the stimulus. Our results provide valuable information towards elucidating the molecular mechanisms of salt/alkali tolerance in Medicago truncatula and provide insight into the putative role of miRNAs in plant stress resistance.

Related to ABSCISIC ACID INSENSITIVE3 (ABI3)/VIVIPAROUS1(VP1)(RAV) transcription factors, which encode a B3 domain and an APETALA2(AP2) domain, belong to the APETALA2/ethylene-responsive element binding factor(AP2/ERF) or B3 superfamily and play an important role in regulating plant growth and development and responding to abiotic stress. Although there have been many functional studies on RAV, the functional differences between RAVs are not clear. Therefore, in this study, the functional differences of RAVs of Medicago truncatula were analyzed. Based on sequence data from the plant transcription factor database and the M. truncatula genome database, we cloned three RAV genes from M. truncatula, named MtRAV1, MtRAV2, and MtRAV3. The cis-acting elements of these genes promoters were predicted, and the expression patterns of MtRAVs under exogenous conditions (4°C, NaCl, Polyethylene Glycol, Abscisic acid) were analyzed. MtRAVs transgenic Arabidopsis thaliana were obtained and subjected to adversity treatment. Subcellular localization results indicated that MtRAVs were located in the nucleus. A much lower expression level was observed for MtRAV3 than the levels of MtRAV1 and MtRAV2 in M. truncatula for growth in normal conditions, but under 4°C or PEG and NaCl treatment, the expression level of MtRAV3 was significantly increased. Only the MtRAV3 overexpression transgenic plants showed strong cold resistance, but the overexpressed MtRAV1 and MtRAV2 transgenic plants showed no difference from wild type plants. MtRAV transgenic plants exhibited similar response to exogenous mannitol, NaCl, and ABA, and the expression of some adverse-related marker genes were up-regulated, such as COLD REGULATED 414 THYLAKOID MEMBRANE 1 (COR414-TM1), Arabidopsis thaliana drought-induced 21 (AtDI21), and Arabidopsis thaliana phosphatidylinositol-specific phospholipase C (ATPLC). MtRAVs transgenic Arabidopsis thaliana exhibited increasing of branch number. These results indicated that there was some function redundancy during MtRAVs proteins of M. truncatula, and MtRAV3 has increased function compared to the other two genes. The results of this study should provide the foundation for future application of MtRAVs in legumes.