The Chromosome-Level Genome Sequence of the Autotetraploid Alfalfa and Resequencing of Core Germplasms Provide Genomic Resources for Alfalfa Research

Advances in basic biology of alfalfa (Medicago sativa L.): a comprehensive overview

Alfalfa (Medicago sativa L.), a perennial legume forage, has been broadly cultivated owing to a variety of favorable characteristics, including comprehensive ecological adaptability, superior nutritive value and palatability, and nitrogen fixation capacity. The productivity traits of alfalfa, specifically its biomass yield and forage quality, are significantly influenced by a series of determinants, including internal developmental factors and external environmental cues. However, the regulatory mechanisms underlying the fundamental biological problems of alfalfa remain elusive. Here, we conducted a comprehensive review focusing on the genomics of alfalfa, advancements in gene-editing technologies, and the identification of genes that control pivotal agronomic characteristics, including biomass formation, nutritional quality, flowering time, and resistance to various stresses. Moreover, a molecular design roadmap for the 'ideal alfalfa' has been proposed and the potential of pangenomes, self-incompatibility mechanisms, de novo domestication, and intelligent breeding strategies to enhance alfalfa's yield, quality, and resilience were further discussed. This review will provide comprehensive information on the basic biology of alfalfa and offer new insights for the cultivation of ideal alfalfa.

Complete mitochondrial genome of Medicago sativa ssp. falcata (Papilionoideae, Fabaceae): characterization and phylogenetic analysis

MAIN CONCLUSION

The first mitogenome of Medicago sativa ssp. falcata complete assembly and a comparative analysis of the four Medicago species base on mitogenome reveal taxonomic insights. Medicago sativa ssp. falcata is primarily distributed in the northern part of the geographical range where alfalfa grows and is a subspecies of the Medicago sativa complex (also called Mediacago falacta). However, compared to M. sativa, M. falcata has better performance in cold resistance and drought tolerance, making it a high-quality gene source for the breeding improvement of Medicago species. We sequenced and assembled the mitochondrial genome of M. falcata with a length of 307,026 bp and successfully annotated 50 genes, of which nad2 exhibited high nucleotide polymorphism in four Medicago species. A total of 197 RNA-editing sites were predicted across 24 protein-coding genes, with alterations at these editing sites resulting in a substantial number of leucine-coding sites, which is consistent with the results of codon usage bias. In addition, we conducted a horizontal comparison of four types of Medicago, including Medicago truncatula, and found that repetitive sequences in their mitogenomes exhibited consistent distribution characteristics. Phylogenetic trees generated through two methods indicated the independent genetic status of M. falcata within the Medicago genus and its partial kinship relationships within the Fabaceae family. The analysis of non-synonymous and synonymous substitution rates of shared protein-coding genes in different plants, along with gene transfer results, suggests that the mitogenome of M. falcata evolved smoothly without showing phases of intense change. This study provides useful information for further understanding the genetic background of M. falcata, with the expectation of contributing to the genomic mining and utilization of germplasm resources in the Medicago genus.

Genome-wide identification of the Medicago sativa L. MYB family and its transcriptional dynamics during pollen development

BACKGROUND

The myeloblastosis (MYB) gene family plays crucial roles in the development of anthers and the establishment of pollen morphology during plant growth. However, little is known about the role of MYB transcription factors in pollen development in alfalfa (Medicago sativa L.).

RESULTS

In this study, we identified 161 MsMYBs in the alfalfa genome, including 34 1R-MYBs, 123 R2R3-MYBs, 3 3R-MYBs, and 1 4R-MYBs (categorized by the number of repeats). These were classified into six subfamilies based on the phylogenetic analysis, conserved structural domains, and gene structures. All MsMYBs were predicted to be hydrophilic and localized in the cell nucleus. The promoter regions contained three classes of cis-regulatory elements related to pollen development, as well as a variable set of functionally diverse elements, including hormone responsiveness, growth and development, and stress responsiveness elements. A transcriptome and qRT-PCR analysis revealed 12 MsMYBs with anther-specific expression and exhibited distinct expression patterns. Some MsMYBs showed a close phylogenetic relationship with Arabidopsis MYBs related to pollen development, such as MsMYB49 and MsMYB100, were found to be localized in the nucleus upon subcellular localization analysis. This genetic proximity suggests a potential role for these MsMYBs in the developmental processes of pollen.

CONCLUSIONS

This study provides a comprehensive understanding of MsMYBs in alfalfa and elucidates their potential roles and expression patterns in pollen development.

Pan-genomic analysis highlights genes associated with agronomic traits and enhances genomics-assisted breeding in alfalfa

Alfalfa (Medicago sativa L.), a globally important forage crop, is valued for its high nutritional quality and nitrogen-fixing capacity. Here, we present a high-quality pan-genome constructed from 24 diverse alfalfa accessions, encompassing a wide range of genetic backgrounds. This comprehensive analysis identified 433,765 structural variations and characterized 54,002 pan-gene families, highlighting the pivotal role of genomic diversity in alfalfa domestication and adaptation. Key structural variations associated with salt tolerance and quality traits were discovered, with functional analysis implicating genes such as MsMAP65 and MsGA3ox1. Notably, overexpression of MsGA3ox1 led to a reduced stem-leaf ratio and enhanced forage quality. The integration of genomic selection and marker-assisted breeding strategies improved genomic estimated breeding values across multiple traits, offering valuable genomic resources for advancing alfalfa breeding. These findings provide insights into the genetic basis of important agronomic traits and establish a solid foundation for future crop improvement.

Systematic characterization of the calmodulin-like (CML) gene family in alfalfa and functional analysis of MsCML70 under salt stress

Calmodulin-like proteins (CMLs), which are widely involved in various abiotic stress responses, are important calcium ion sensors in plants. However, systematic identification and functional analysis of these proteins have not been performed in alfalfa. Here, a total of 211 MsCMLs were identified in the alfalfa genome. Conserved domain analysis revealed that most MsCMLs contained three EF-hand domains. A total of 17 tandem duplication events and 292 segmental duplication events were identified, indicating that segmental duplications were the major factor in the expansion of MsCMLs. There were 28, 36 and 18 MsCMLs that responded to drought, salt and cold stress, respectively, in alfalfa. In addition, MsCML70 overexpression significantly increased salt tolerance in Arabidopsis. MsCML70 participates in the plant salt stress response through various biological pathways, including transcriptional regulation, protein modification, plant hormone metabolism and secondary metabolism. Moreover, MsCML70 significantly increased the expression of HKT1 (high-affinity K+transporter 1), DREB19 (dehydration responsive element binding protein 19), PRX32 (peroxidase 32), JAL10 (jacalin-associated lectins 10), HB17 (homeobox 17), and NPF2.3 (nitrate transporter 2.3) under salt stress to promote tolerance to salt stress in Arabidopsis. The results of this study help elucidate the function of alfalfa CML genes and provide a new gene resource for the breeding of stress-resistant alfalfa.

Pectins play a central role in enhancing Al tolerance of alfalfa via looseing fibre-microfiber arrangement of cell wall in root tips

Cell wall greatly affects Al tolerance of plants, but the precise mechanisms by which the cell wall modulating Al tolerance remains largely unknown. In the present study, Al tolerant alfalfa varieties (WL525 and WL903) accumulated less Al in root tips, cell wall and pectins, averagely decreased by 23.8 %, 41.7 %, and 80.3 %, respectively, than Al-sensitive varieties (WL440 and GN3) under Al treatment. Pectins were more susceptible to Al stress than hemicellulose in cell wall, in which water-soluble pectin (WSP) was most susceptible among the three pectin forms (WSP, CSP and ASP). The decreased pectin and WSP contents exhibited a positive relationship with PG activity and five MsPGs expressions, leading to a smaller WSP polymer, a looser fibre-microfiber arrangement, and a larger porosity of cell wall in the root tips of WL525 than WL440 under Al treatment. Consequently, the Al-induced inhibitor rates of root growth were lower, averagely decreased by 1.9 times, in Al tolerant varieties than Al sensitive varieties. Similarly, the pectin hydrolysis and root length increased in overexpressed MsPG4, but decreased in downregulated MsPG4 transgenic alfalfa plants under Al stress. These results revealed that pectins, especially WSP, play a central role in sustaining cell wall extensibility and increasing Al tolerance of alfalfa.

Anther-specific expression of MsMYB35 transcription factor in alfalfa (Medicago sativa L.) and its crucial role in pollen development

Alfalfa (Medicago sativa L.) is a high-quality forage crop and an essential resource for livestock. Understanding the molecular mechanisms underlying male sterility in alfalfa is pivotal for the development of superior forage varieties. Despite the critical role of anther development in plant reproduction, its molecular regulation-particularly the involvement of transcription factors in M. sativa-remains insufficiently explored. This study bridges this gap by isolating and characterizing an R2R3-MYB transcription factor, MsMYB35, and unveiling its regulatory role in anther development. Quantitative RT-PCR (qRT-PCR) revealed that MsMYB35 is predominantly expressed during early anther development and is homologous to AtMYB35. MsMYB35 was found to localize in both the cytoplasm and nucleus. DNA affinity purification sequencing (DAP-seq) identified 3647 target genes of MsMYB35, with enrichment analysis uncovering three recognition motifs. Integrated DAP-seq and RNA-seq analyses revealed that MsMYB35 directly regulates two key anther development-related genes. Functional analyses showed that overexpression of MsMYB35 promotes anther development, while silencing MsMYB35 leads to defective anther sacs and wrinkled pollen grains. Proper MsMYB35 expression ensures the formation of viable and fertile pollen grains, solidifying its role as a critical regulator of anther development. These findings provide a novel perspective on the molecular mechanisms regulating anther development in M. sativa and offer valuable insights for improving molecular breeding and hybrid seed production strategies. By advancing the fundamental understanding of transcriptional regulation in anther development, this study sets the stage for innovative approaches to alfalfa crop improvement.

The Involvement of Amino Acid Metabolism in the Mechanisms of Salt Tolerance Adaptation in Medicago sativa and Medicago truncatula

Amino acid metabolism constitutes a major metabolic pathway in plants, playing an important role in the modulation of plant responses to stress. In this study, we investigated the amino acid metabolism responses of M. sativa (Medicago sativa L.) and M. truncatula (Medicago truncatula L.) plants under salt stress using transcriptomic and proteomic approaches to elucidate their salt stress tolerance mechanisms in relation to the regulation of amino acid homeostasis. Transcriptome and proteome sequencing followed by Kyoto Gene and Genome Encyclopedia enrichment analysis revealed 34 differentially expressed genes and 45 differentially expressed proteins involved in valine, leucine, and isoleucine degradation, tyrosine metabolism, and glutathione metabolism. Significant differences were observed in the expression of glutathione S-transferase (GST) within the glutathione metabolic pathway between M. sativa and M. truncatula. The induction of valine, leucine, and isoleucine metabolism, aldehyde dehydrogenases (ALDHs), and alanine-glyoxylate aminotransferases (AGXTs), involved in intracellular reactive oxygen species scavenging, also significantly differed under salt stress. Significant differences were identified in the expression of tyrosine decarboxylases (TDCs) involved in tyrosine metabolism, which are responsible for tyramine biosynthesis and can enhance plant tolerance to salt stress. This study delved into the effects of amino acid metabolism on the salt tolerance mechanisms of M. sativa and M. truncatula, which is crucial in guiding the future breeding of salt-tolerant alfalfa varieties.

Genome-wide identification of the HIPPs gene family and functional validation of MsHIPP12 in enhancing cadmium tolerance in Medicago sativa

Heavy metal-associated isoprenylated plant proteins (HIPPs) are crucial for metal ion homeostasis and stress responses in plants exposed to heavy metals. They bind heavy metal ions via their HMA domains, sequestering them to prevent cellular toxicity. The C-terminal isoprenylation enhances interactions with membrane proteins, aiding in ion transport and compartmentalization. In Medicago sativa (alfalfa), we identified 23 MsHIPP genes containing conserved HMA domains and C-terminal isoprenylation motifs using bioinformatics tools. Phylogenetic analysis classified these genes into five clades, indicating functional diversity and evolutionary divergence. Promoter analysis revealed cis-regulatory elements associated with responses to light, drought, cold, abscisic acid (ABA), salicylic acid (SA), and auxin, suggesting roles in environmental adaptation. Expression profiling under drought, cold, salt, cadmium (Cd), ABA, and indole-3-acetic acid (IAA) treatments demonstrated the involvement of MsHIPPs in abiotic stress responses. Notably, overexpression of MsHIPP12 in Arabidopsis thaliana enhanced Cd tolerance by increasing antioxidant enzyme activities (APX, CAT, SOD), reducing malondialdehyde (MDA) levels, and attenuating chlorophyll degradation. These findings identify MsHIPP12 as a promising candidate gene for enhancing stress tolerance in M. sativa. This work provides valuable insights for molecular breeding strategies and phytoremediation approaches to address cadmium-contaminated soils.

Optimization of high-throughput marker systems for genomic prediction in alfalfa family bulks

Alfalfa (Medicago sativa L.) is a perennial forage legume esteemed for its exceptional quality and dry matter yield (DMY); however, alfalfa has historically exhibited low genetic gain for DMY. Advances in genotyping platforms paved the way for a cost-effective application of genomic prediction in alfalfa family bulks. In this context, the optimization of marker density holds potential to reallocate resources within genomic prediction pipelines. This study aimed to (i) test two genotyping platforms for population structure discrimination and predictive ability (PA) of genomic prediction models (G-BLUP) for DMY, and (ii) explore optimal levels of marker density to predict DMY in family bulks. For this, 160 nondormant alfalfa families were phenotyped for DMY across 11 harvests and genotyped via targeted sequencing using Capture-seq with 17K probes and the DArTag 3K panel. Both platforms discriminated similarly against the population structure and resulted in comparable PA for DMY. For genotyping optimization, different levels of marker density were randomly extracted from each platform. In both cases, a plateau was achieved around 500 markers, yielding similar PA as the full set of markers. For phenotyping optimization, models with 500 markers built with data from five harvests resulted in similar PA compared to the full set of 11 harvests and full set of markers. Altogether, genotyping and phenotyping efforts were optimized in terms of number of markers and harvests. Capture-seq and DArTag yielded similar results and have the flexibility to adjust their panels to meet breeders' needs in terms of marker density.

Selection and validation of reference genes in alfalfa based on transcriptome sequence data

Compared with the traditional gene expression techniques of quantitative analysis, RT-qPCR is most widely used because of its low cost, time-saving, and accuracy. It is essential to introduce suitable internal reference genes as reference corrections in RT-qPCR experiments to reduce the RNA quality and reverse transcription efficiency of different samples. In our study, we chose the candidate internal reference genes of alfalfa from transcriptome sequence datasets (162 RNA-seq sequencing data) through comparative analysis. Finally, 10 candidate reference genes were selected. These candidate reference gene expressions were determined by RT-qPCR under five abiotic stresses of drought, alkali, high temperature, low temperature, and acid treatments. The stability index of these candidate genes was calculated and evaluated correspondently using specific software and different lgorithms, such as GeNorm, Normfinder, Bestkeeper, ΔCt method, and an online analysis tool RefFinder. Then the appropriate candidate genes were screened strictly; and validated by the phyA gene. GAPDH and Actin are taken as traditional reference genes on gene expression of quantitative analysis commonly used in alfalfa, Our results showed GAPDH and Actin aren't the most appropriate reference genes of alfalfa under different abiotic stresses, under alkaline stress, the optimal reference gene is UBL-2a, and the optimal combination of reference genes is MS.65,463 (some candidate reference genes haven't been annotated yet, using gene ID abbreviation number of Medicago sativa L. instead )and UBL-2a; Under drought stress, the optimal reference gene is Ms.33,066, and the optimal combination of reference genes is MS.65,463 and UBL-2a; Under high-temperature stress, the optimal reference gene is Actin, and the optimal combination of reference genes is Rer1, Actin, MS.00617, MS.74,923, UBL-2a, MS.33,066, MS.99,505, and MS.65,463; Under low-temperature stress, the optimal reference gene is Actin, and the optimal combination of reference genes is Rer1, Actin, MS.99,505, MS.073307, and UBL-2a. The optimal reference genes and their combinations need further validation under acid stress. This paper provides scientific evidence for quantitative analysis of the genes of alfalfa.

Integrated Analyses of the Mechanism of Flower Color Formation in Alfalfa (Medicago sativa)

BACKGROUND

Alfalfa (Medicago sativa) is one of the most valuable forages in the world. As an outcrossing species, it needs bright flowers to attract pollinators to deal with self-incompatibility. Although various flower colors have been observed and described in alfalfa a long time ago, the biochemical and molecular mechanism of its color formation is still unclear.

METHODS

By analyzing alfalfa lines with five contrasting flower colors including white (cream-colored), yellow, lavender (purple), dark purple and dark blue, various kinds and levels of anthocyanins, carotenoids and other flavonoids were detected in different colored petals, and their roles in color formation were revealed.

RESULTS

Notably, the content of delphinidin-3,5-O-diglucoside in lines 3, 4 and 5 was 58.88, 100.80 and 94.07 times that of line 1, respectively. Delphinidin-3,5-O-diglucoside was the key factor for purple and blue color formation. Lutein and β-carotene were the main factors for the yellow color formation. By analyzing differentially expressed genes responsible for specific biochemical pathways and compounds, 27 genes were found to be associated with purple and blue color formation, and 14 genes were found to play an important role in yellow color formation.

CONCLUSIONS

The difference in petal color between white, purple and blue petals was mainly caused by the accumulation of delphinidin-3,5-O-diglucoside. The difference in petal color between white and yellow petals was mainly affected by the production of lutein and β-carotene. These findings provide a basis for understanding the biochemical and molecular mechanism of alfalfa flower color formation.

Whole-genome sequencing revealed genetic basis of diterpenoid alkaloid difference in Aconitum vilmorinianum

BACKGROUND

Aconitum is an important medicinal genus widely used in traditional Chinese medicine, which produces types of diterpenoid alkaloids (DA) among different species. We performed whole genome resequencing (WGS) research in Aconitum spp., and wish to find diterpenoid alkaloids related genetic variations.

RESULTS

In this study, we re-sequenced 150 Aconitum vilmorinianum (A. vilmorinianum) including 102 from the cultivation garden and 48 from the wild, as well as nine wild samples of Aconitum weixiense. The intra-population differentiation of A. vilmorinianum was detected by evolutionary tree and population structure inference. We identify 47 DA biosynthesis genes that might be highly associated with the specialization of DA based on whole-genome resequencing. Of 616 significant SNPs and 105 significant InDels among these genes could be developed as polymorphic molecular markers capable of effectively recognizing A. vilmorinianum from A. weixiense. Furthermore, the significant SNPs and InDels were almost homozygous alternates in A. weixiense, whereas they tended to be homozygous references in the A. vilmorinianum.

CONCLUSIONS

Our results discussed the difference in genetic background in A. vilmorinianum compared to A. weixiense and these high-quality DA biosynthesis-associated polymorphic locus provided useful genetic information for discrimination of A. vilmorinianum and could serve as a vehicle to study the mechanism of DA differentiation in Aconitum.

Medicago2035: Genomes, functional genomics, and molecular breeding

Medicago, a genus in the Leguminosae or Fabaceae family, includes the most globally significant forage crops, notably alfalfa (Medicago sativa). Its close diploid relative Medicago truncatula serves as an exemplary model plant for investigating legume growth and development, as well as symbiosis with rhizobia. Over the past decade, advances in Medicago genomics have significantly deepened our understanding of the molecular regulatory mechanisms that underlie various traits. In this review, we comprehensively summarize research progress on Medicago genomics, growth and development (including compound leaf development, shoot branching, flowering time regulation, inflorescence development, floral organ development, and seed dormancy), resistance to abiotic and biotic stresses, and symbiotic nitrogen fixation with rhizobia, as well as molecular breeding. We propose avenues for molecular biology research on Medicago in the coming decade, highlighting those areas that have yet to be investigated or that remain ambiguous.

The mixed auto-/allooctoploid genome of Japanese knotweed (Reynoutria japonica) provides insights into its polyploid origin and invasiveness

Reynoutria japonica Houtt. (Polygonaceae), a traditional Chinese medicine, is one of the top 100 most destructive invasive species worldwide due to its aggressive growth and strong adaptability. Here, we report an 8.04 Gb chromosome-scale assembly of R. japonica with 88 chromosomes across eight homologous sets. Through a combined phylogenetic and genomic analysis, we demonstrate that R. japonica is a mixed auto-/allooctoploid (AAAABBBB). Subgenome A (SubA) exhibited a close phylogenetic relationship with the related species Fallopia multiflora. We also unveiled the origin and evolutionary history of octoploid R. japonica based on resequencing data from Reynoutria species with different ploidy. Comparative genomics analysis revealed the genetic basis of R. japonica's invasivity and adaptability. The auxin response factor (ARF) gene family was significantly expanded in R. japonica, and these genes were highly expressed in rhizomes. We also investigated the collaboration and differentiation of the duplicated genes resulting from auto- and allo-polyploidization at the genomic variation, gene family evolution, and gene expression levels in R. japonica. Transcriptomic analysis of stem internodes and apices at different developmental stages revealed that the octuplication and significant expansion of the SAUR19 and SAUR63 subfamilies due to tandem replication in SubA, and the high expression of these genes in stems, likely contribute to the rapid growth of R. japonica. Our study provides important clues into adaptive evolution and polyploidy dominant traits in invasive plants, and will also provide important guidance for the breeding of polyploid crops.

Genome-wide characterization and expression analysis of the ADF gene family in response to salt and drought stress in alfalfa (Medicago sativa)

The microfilament cytoskeleton, formed by the process of actin polymerization, serves not only to support the morphology of the cell, but also to regulate a number of cellular activities. Actin-depolymerizing factors (ADFs) represent a significant class of actin-binding proteins that regulate the dynamic alterations in the microfilament framework, thereby playing a pivotal role in plant growth and development. Additionally, they are instrumental in modulating stress responses in plants. The ADF gene family has been explored in various plants, but there was a paucity of knowledge regarding the ADF gene family in alfalfa (Medicago sativa), which is one of the most significant leguminous forage crops globally. In this study, a total of nine ADF genes (designated MsADF1 through MsADF9) were identified in the alfalfa genome and mapped to five different chromosomes. A phylogenetic analysis indicated that the MsADF genes could be classified into four distinct groups, with members within the same group exhibiting comparable gene structures and conserved motifs. The analysis of the Ka/Ks ratios indicated that the MsADF genes underwent purity-based selection during its evolutionary expansion. The promoter region of these genes was found to contain multiple cis-acting elements related to hormone responses, defence, and stress, indicating that they may respond to a variety of developmental and environmental stimuli. Gene expression profiles analyzed by RT-qPCR experiments demonstrated that MsADF genes exhibited distinct expression patterns among different organs. Furthermore, the majority of MsADF genes were induced by salt and drought stress by more than two-fold, with MsADF1, 2/3, 6, and 9 being highly induced, suggesting their critical role in resistance to abiotic stress. These results provide comprehensive information on the MsADF gene family in alfalfa and lay a solid foundation for elucidating their biological function.

Morphological diversity variation of seed traits among 587 germplasm resources of Medicago Genus and 32 germplasm resources of Trigonella Genus

Germplasm resources within the Medicago genus are highly regarded for their value as forage crops and their critical roles in nitrogen cycling, ecosystem restoration, and soil structure improvement. Therefore, understanding the diversity of seed morphology in this genus is essential for advancing its development and utilization. This study analyzed seed samples from 587 germplasm accessions representing 77 species within Medicago genus, as well as 32 accessions from 21 species within the closely related genus Trigonella. A statistical analysis was conducted on twelve traits, including seven quantitative traits-straight length (SL), straight width (SW), width-to-length ratio (WL), perimeter (PE), radicle length (RL), hilum length (HL), and 100-seed weight (SY)-and five qualitative traits, including seed coat condition, radicle characteristics, seed size, shape, and color. The results revealed that: (1) there was significant diversity (P < 0.05) in SL, SW, WL, PE, RL, HL, and SY across Medicago species; (2) principal component analysis of the 587 Medicago accessions identified SL, SW, PE, HL, RL, and SY as the primary contributors to morphological diversity; and (3) high-resolution images of seeds from various accessions were captured for future research. This study provides a solid foundation for the establishment of seed banks and the enhancement of germplasm resources through the systematic analysis of these morphological traits.

Comprehensive analysis of epigenetic modifications in alfalfa under cadmium stress

Epigenetics plays an important role in plant growth and development and in environmental adaptation. Alfalfa, an important forage crop, is rich in nutrients. However, little is known about the molecular regulatory mechanisms underlying the response of alfalfa to cadmium (Cd) stress. Here, we performed DNA methylation (5mC), RNA methylation (m6A) and transcriptomic sequencing analyses of alfalfa roots under Cd stress. Whole-genome methylation sequencing and transcriptomic sequencing revealed that Cd stress reduced DNA methylation levels. Moreover, a reduced 5mC methylation level was associated with decreased expression of several DNA methyltransferase genes. Compared with those under normal (CK) conditions, the m6A modification levels under Cd stress were greater and were positively correlated with gene expression in alfalfa roots. We also found a negative correlation between the 5mC level and the m6A level, especially in CG and CHG contexts. In yeast, the overexpression of MsNARMP5 (natural resistance-associated macrophage protein) and MsPCR2 (plant cadmium resistance 2), which are modified by 5mC or m6A, significantly increased Cd stress tolerance. These results provide candidate genes for future studies on the mechanism of Cd stress tolerance in alfalfa roots and valuable information for studying heavy metal stress in alfalfa breeding.

Multi-omic analysis on the molecular mechanisms of rapid growth in 'Deqin' alfalfa after space mutagenesis

BACKGROUND

Space-induced plant mutagenesis, driven by cosmic radiation, offers a promising approach for the selective breeding of new plant varieties. By leveraging the unique environment of outer space, we successfully induced mutagenesis in 'Deqin' alfalfa and obtained a fast-growing mutant. However, the molecular mechanisms underlying its rapid growth remain poorly unexplored.

RESULTS

Comparative analyses of transcriptomics, proteomics, and hormone profiles were conducted in root, stem, and leaf tissues of both mutant and non-mutagenic materials. Targeted plant hormone showed notable increases in the levels of 3-indolebutyric, indole-3-acetic, and 3-indolepropionic acids in the mutant, with percentage increases of 33.55%, 32.49%, and 30.39%, respectively. Zeatin-riboside and dihydrozeatin riboside levels increased by 164.92% and 25.86%, while giberellin (GA) 7, GA3, and GA1 levels increased by 219.52%, 68.74%, and 40.98%. Non-mutagenic materials sprayed with exogenous 3-indolebutyric acid, zeatin-riboside, and GA7 exhibited significant growth acceleration. Transcriptomics identified 49,095 annotated genes, with 2,009, 1,889, and 1,760 upregulated and 2,082, 2,035, and 1,499 downregulated in the leaves, stems, and roots, respectively. Twenty-two genes related to plant hormone biosynthesis showed significant alterations. Screening through weighted correlation network analysis revealed ten candidate genes, four of which were associated with photosynthesis and starch and sucrose metabolism. Integrated analysis of targeted plant hormone metabolomics and transcriptomics indicated that plant hormone signal transduction played a crucial role. Proteomics revealed 479 differentially accumulated proteins, of which 174 were upregulated and 305 were downregulated. Integrated proteomics and transcriptomics showed that photosynthesis, starch and sucrose metabolism, carbon metabolism, and carbon fixation in photosynthetic organisms promoted the rapid growth of the mutants. By integrating multi-omics data, we elucidated the synergistic effects of pathways such as hormone signal transduction and tryptophan metabolism on the rapid growth of the mutants.

CONCLUSION

This study demonstrated the significance of plant hormones in the rapid growth of the mutants and identified key genes and metabolic pathways. Our findings provide valuable information for the genetic improvement of alfalfa varieties and serve as a reference for achieving rapid growth in other plants.

Pre-breeding in alfalfa germplasm develops highly differentiated populations, as revealed by genome-wide microhaplotype markers

Plant genebanks contain large numbers of germplasm accessions that likely harbor useful alleles or genes absent in commercial plant breeding programs. Broadening the genetic base of commercial alfalfa germplasm with these valuable genetic variations can be achieved by screening the extensive genetic diversity in germplasm collections and enabling maximal recombination among selected genotypes. In this study, we assessed the genetic diversity and differentiation of germplasm pools selected in northern U.S. latitudes (USDA Plant Hardiness Zone 7 or below) originating from Eurasian germplasm. The germplasm evaluated included four BASE populations (C0) from different geographical origins (Central Asia, Northeastern Europe, Balkans-Turkey-Black Sea, and Siberia/Mongolia), 20 cycle-one populations (C1) derived from each of the four BASE populations selected across five locations in the U.S. and Canada, and four commercial cultivars. Using a panel of 3,000 Diversity Array Technologies (DArTag) marker loci, we retrieved 2,994 target SNPs and approximately 12,000 microhaplotypes. Microhaplotypes exhibited higher genetic diversity values than target SNPs. Principal component analysis and discriminant analysis of principal components revealed significant population structure among the alfalfa populations based on geographical origin, while the check cultivars formed a central cluster. Inbreeding coefficients (FIS) ranged from - 0.1 to 0.006, with 27 out of 28 populations showing negative FIS values, indicating an excess of heterozygotes. Interpopulation genetic distances were calculated using Rho pairwise distances (FST adapted for autotetraploid species) and analysis of molecular variance (AMOVA) parameters. All BASE populations showed lower Rho values compared to C1 populations and check cultivars. AMOVA revealed that most of the genetic diversity was among individuals within populations, especially in BASE populations (92.7%). This study demonstrates that individual plants in BASE populations possess high genetic diversity, low interpopulation distances, and minimal inbreeding, characteristics that are essential for base-broadening selection. The populations developed in this project serve as valuable sources of novel alleles for North American alfalfa breeding programs, offering breeders access to diverse, regionally adapted pools for improving various alfalfa traits.

Oligo-FISH barcode chromosome identification system provides novel insights into the natural chromosome aberrations propensity in the autotetraploid cultivated alfalfa

Alfalfa is one of the most economically valuable forage crops in the world. However, molecular cytogenetic studies in alfalfa lag far behind other cash crops and have reached a bottleneck. Here, we developed a novel chromosome identification system by designing 21 oligo probes in specific regions of each chromosome, which can be used as a barcode to simultaneously distinguish all chromosomes in a cell. Using this system, we revealed the chromosome karyotype features and evolutionary differences among 10 cultivated alfalfa varieties. Interestingly, we also found two chromosomal variation types, i.e. aneuploidy and large chromosomal segment deletions in the seeds of three alfalfa varieties. Variation frequency analysis showed that only 7/173 seeds in those three alfalfa varieties had chromosome aberrations, which indicated that the inheritance and meiosis of alfalfa had evolved to a relatively stable state. Remarkably, 4/7 variation seeds were chromosome 2 aberrations, suggesting that chromosome 2 appears to be more susceptible to natural chromosomal aberrations than other chromosomes during inheritance. DNA sequence variation analysis showed that the difference of presence and absence variations (PAVs) among homologous copies of chromosome 2 was larger than that of the other seven chromosomes. We suggest that such large PAV divergence among homologous copies may provide the physical basis for natural chromosome 2 aberrations propensity. Our study provides a valuable and efficient tool for alfalfa's molecular cytogenetics and sheds new insights into the propensity for natural chromosome aberrations during autopolyploid inheritance.

Refining polyploid breeding in sweet potato through allele dosage enhancement

Allele dosage plays a key role in the phenotypic variation of polyploids. Here we present a genome-wide variation map of hexaploid sweet potato that captures allele dosage information, constructed from deep sequencing of 294 hexaploid accessions. Genome-wide association studies identified quantitative trait loci with dosage effects on 23 agronomic traits. Our analyses reveal that sweet potato breeding has progressively increased the dosage of favourable alleles to enhance trait performance. Notably, the Mesoamerican gene pool has evolved towards higher dosages of favourable alleles at multiple loci, which have been increasingly introgressed into modern Chinese cultivars. We substantiated the breeding-driven dosage accumulation through transgenic validation of IbEXPA4, an expansin gene influencing tuberous root weight. In addition, we explored causative sequence variations that alter the expression of the Orange gene, which regulates flesh colour. Our findings illuminate the breeding history of sweet potato and establish a foundation for leveraging allele dosages in polyploid breeding practices.

Identification of dehydrin family genes in three Medicago species and insights into their tolerant mechanism to salt stress

KEY MESSAGE

All ten dehydrin genes from three Medicago species are responsive to different kinds of abiotic stress, and CAS31 confers transgenic plants salt tolerance by down-regulating HKT1 expression. Dehydrins are protective proteins playing crucial roles in the tolerance of plants to abiotic stresses. However, a full-scale and systemic analysis of total dehydrin genes in Medicago at the genome level is still lacking. In this study, we identified ten dehydrin genes from three Medicago species (M. truncatula, M. ruthenica, and M. sativa), categorizing the coding proteins into four types. Genome collinearity analysis among the three Medicago species revealed six orthologous gene pairs. Promoter regions of dehydrin genes contained various phytohormone- and stress-related cis-elements, and transcriptome analysis showed up-regulation of all ten dehydrin genes under different stress conditions. Transformation of dehydrin gene CAS31 increased the tolerance of transgenic seedlings compared with wild-type seedlings under salt stress. Our study demonstrated that transgenic seedlings maintained the more chlorophyll, accumulated more proline and less hydrogen peroxide and malondialdehyde than wild-type seedlings under salt stress. Further study revealed that CAS31 reduced Na+ accumulation by down-regulating HKT1 expression under salt stress. These findings enhance our understanding of the dehydrin gene family in three Medicago species and provide insights into their mechanisms of tolerance.

Genome-Wide Identification of the GPAT Family in Medicago sativa L. and Expression Profiling Under Abiotic Stress

Glycerol-3-phosphate acyltransferase (GPAT), as a rate-limiting enzyme engaged in lipid synthesis pathways, exerts an important role in plant growth and development as well as environmental adaptation throughout diverse growth stages. Alfalfa (Medicago sativa L.) is one of the most significant leguminous forages globally; however, its growth process is frequently exposed to environmental stress, giving rise to issues such as impeded growth and decreased yield. At present, the comprehension of the GPAT genes in alfalfa and their reactions to abiotic stresses is conspicuously deficient. This study identified 15 GPATs from the genome of "Zhongmu No. 1" alfalfa, which were phylogenetically categorized into three major groups (Groups I ~ III). Furthermore, Group III is further subdivided into three distinct subgroups. MsGPATs belonging to the same subfamily exhibited similar protein conserved motifs and gene structural characteristics, in which groups with simple conserved motifs had more complex gene structures. A multitude of regulatory cis-elements pertinent to hormones and responses to environmental stress were detected in their promoter regions. In addition, a spatial-temporal expression analysis showed that MsGPATs have significant tissue specificity. Furthermore, the transcriptomic analysis of ABA treatment and the qRT-PCR results under drought, salt, and cold stress demonstrated that the majority of MsGPATs respond to abiotic stress with pronounced timely characteristics. It was also ascertained that these GPAT genes might assume a crucial role in salt and drought stress. This research can further constitute a fundamental basis for the exploration of the alfalfa GPAT family, the screening of key GPATs, and the investigation of their functionalities.

Saliva of Therioaphis trifolii (Drepanosiphidae) Activates the SA Plant Hormone Pathway, Inhibits the JA Plant Hormone Pathway, and Improves Aphid Survival Rate

The spotted alfalfa aphid (Therioaphis trifolii) is a kind of destructive pest of cultivated alfalfa (Medicago sativa) that results in significant financial losses for the livestock sector. To understand how T. trifolii navigates the biochemical defenses of its host, we investigated the effects of susceptible and resistant aphid strains on two alfalfa cultivars. T. trifolii was reared for over 50 generations on two cultivars-WL343, which is susceptible to T. trifolii, and Zhongmu No. 1, which is resistant-resulting in the development of a resistant aphid strain (R-aphid) and a susceptible aphid strain (S-aphid). The results showed that the survival rate of R-aphids was significantly higher than that of S-aphids (p = 0.039) on the resistant cultivar Zhongmu No. 1, while there was no significant difference in survival rates between the two aphid strains on WL343 (p = 0.139). This suggests that S- and R-aphids differ in their ability to modulate plant defense mechanisms, influencing their survival rates. The application of saliva from R-aphids reared on Zhongmu No. 1 (R-saliva) reduced the repellency and toxicity of treated plants, improving aphid survival. Furthermore, R-aphid infestation and R-saliva application activated the salicylic acid (SA) signaling pathway while suppressing the jasmonic acid (JA) pathway, enhancing the host plant's suitability for aphid colonization. We propose that R-aphids may use their saliva to activate the SA pathway, which in turn inhibits JA synthesis, weakening the plant's defenses. These findings provide valuable insights into how T. trifolii interacts with host plant defense systems.

Functions and Mechanisms of Brassinosteroids in Regulating Crop Agronomic Traits

Brassinosteroids (BRs) perform crucial functions controlling plant growth and developmental processes, encompassing many agronomic traits in crops. Studies of BR-related genes involved in agronomic traits have suggested that BRs could serve as a potential target for crop breeding. Given the pleiotropic effect of BRs, a systematic understanding of their functions and molecular mechanisms is conducive for application in crop improvement. Here, we summarize the functions and underlying mechanisms by which BRs regulate the several major crop agronomic traits, including plant architecture, grain size, as well as the specific trait of symbiotic nitrogen fixation in legume crops. For plant architecture, we discuss the roles of BRs in plant height, branching number and leaf erectness, and propose how progress in these fields may contribute to designing crops with optimal agronomic traits and improved grain yield by accurately modifying BR levels and signaling pathways.

Genome-Wide Identification and Expression Profiling Analysis of the Mitochondrial Calcium Uniporter Gene Family Under Abiotic Stresses in Medicago sativa

The mitochondrial calcium uniporters (MCUs) are a family of calcium unidirectional transporters important for cytoplasmic Ca2+ signals. Though the MCU proteins in several plant species have been investigated, genome-wide analysis of MCUs in alfalfa is lacking. Here, via genome-wide analysis, a total of 5, 20, and 6 MCU genes were identified in three different alfalfa cultivars, namely Medicago truncatula Jemalong A17, Medicago sativa XinJiangDaYe, and M. sativa Zhongmu No. 1, respectively. They were further phylogenetically classified into three subfamilies. Most MCU genes have only one intron, and gene duplication events of MCU genes were observed within each alfalfa accession and between different accessions. All alfalfa MCU proteins contained a highly conserved MCU domain and 10 conserved motifs, featuring two transmembrane domains and a DI/VME motif. According to the tissue expression data of M. sativa Zhongmu No. 1, MsMCU6.2 was the most abundant transcript with the highest expression in the leaf, and MsMCU5 and MsMCU1.2 showed higher expression levels in the stem than other tissues. We analyzed the expression profiles of five MCU genes (MsMCU1.1/1.2/5/6.1/6.2) under salt, drought, and cold stresses via qRT-PCR assays. All five MCU genes were induced by drought stress, except MsMCU5, whose expression was up-regulated by salt stress, while cold stress slightly altered MsMCU expression. Nine potential interacting proteins and three miRNAs targeting MtMCUs were predicted. These results provide detailed knowledge of alfalfa MCU genes and suggest their potential functions in abiotic stress response.

The Effects of Plant-Microbe-Environment Interactions on Mineral Weathering Patterns in a Granular Basalt

The importance of biota to soil formation and landscape development is widely recognized. As biotic complexity increases during early succession via colonization by soil microbes followed by vascular plants, effects of biota on mineral weathering and soil formation become more complex. Knowledge of the interactions among groups of organisms and environmental conditions will enable us to better understand landscape evolution. Here, we used experimental columns of unweathered granular basalt to investigate how early successional soil microbes, vascular plants (alfalfa; Medicago sativa), and soil moisture interact to affect both plant performance and mineral weathering. We found that the presence of soil microbes reduced plant growth rates, total biomass, and survival, which suggests that plants and microbes were competing for nutrients in this environment. However, we also found considerable genotype-specific variation in plant-microbial interactions, which underscores the importance of within-species genetic variation on biotic interactions. We also found that the presence of vascular plants reduced variability in pH and electrical conductivity, suggesting that plants may homogenize weathering reactions across the soil column. We also show that there is heterogeneity in the abiotic conditions in which microbes, plants, or their combination have the strongest effect on weathering, and that many of these relationships are sensitive to soil moisture. Our findings highlight the importance of interdependent effects of environmental and biotic factors on weathering during initial landscape formation.

A stepwise guide for pangenome development in crop plants: an alfalfa (Medicago sativa) case study

BACKGROUND

The concept of pangenomics and the importance of structural variants is gaining recognition within the plant genomics community. Due to advancements in sequencing and computational technology, it has become feasible to sequence the entire genome of numerous individuals of a single species at a reasonable cost. Pangenomes have been constructed for many major diploid crops, including rice, maize, soybean, sorghum, pearl millet, peas, sunflower, grapes, and mustards. However, pangenomes for polyploid species are relatively scarce and are available in only few crops including wheat, cotton, rapeseed, and potatoes.

MAIN BODY

In this review, we explore the various methods used in crop pangenome development, discussing the challenges and implications of these techniques based on insights from published pangenome studies. We offer a systematic guide and discuss the tools available for constructing a pangenome and conducting downstream analyses. Alfalfa, a highly heterozygous, cross pollinated and autotetraploid forage crop species, is used as an example to discuss the concerns and challenges offered by polyploid crop species. We conducted a comparative analysis using linear and graph-based methods by constructing an alfalfa graph pangenome using three publicly available genome assemblies. To illustrate the intricacies captured by pangenome graphs for a complex crop genome, we used five different gene sequences and aligned them against the three graph-based pangenomes. The comparison of the three graph pangenome methods reveals notable variations in the genomic variation captured by each pipeline.

CONCLUSION

Pangenome resources are proving invaluable by offering insights into core and dispensable genes, novel gene discovery, and genome-wide patterns of variation. Developing user-friendly online portals for linear pangenome visualization has made these resources accessible to the broader scientific and breeding community. However, challenges remain with graph-based pangenomes including compatibility with other tools, extraction of sequence for regions of interest, and visualization of genetic variation captured in pangenome graphs. These issues necessitate further refinement of tools and pipelines to effectively address the complexities of polyploid, highly heterozygous, and cross-pollinated species.

Genome-Wide Identification of the RALF Gene Family and Expression Pattern Analysis in Zea mays (L.) under Abiotic Stresses

Rapid Alkalization Factor (RALF) is a signaling molecule in plants that plays a crucial role in growth and development, reproductive processes, and responses to both biotic and abiotic stresses. Although RALF peptides have been characterized in Arabidopsis and rice, a comprehensive bioinformatics analysis of the ZmRALF gene family in maize is still lacking. In this study, we identified 20 RALF genes in the maize genome. Sequence alignment revealed significant structural variation among the ZmRALF family genes. Phylogenetic analysis indicates that RALF proteins from Arabidopsis, rice, and maize can be classified into four distinct clades. Duplication events suggest that the expansion of the RALF gene family in maize primarily relies on whole-genome duplication. ZmRALF genes are widely expressed across various tissues; ZmRALF1/15/18/19 are highly expressed in roots, while ZmRALF6/11/14/16 are predominantly expressed in anthers. RNA-seq and RT-qPCR demonstrated that the expression levels of ZmRALF7, ZmRALF9, and ZmRALF13 were significantly up-regulated and down-regulated in response to PEG and NaCl stresses, respectively. Overall, our study provides new insights into the role of the RALF gene family in abiotic stress.

The mutagenic effect of cold plasma on Medicago sativa L

To explore the feasibility of using cold plasma as a mutagenesis breeding technology for forage crops, in this study we used the Medicago sativa L. cultivar, Zhongcao No. 3, as the experimental material. The effects of plasma treatments on Medicago sativa L. were analyzed through the use of plasma and activated water. Treatments with plasma and activated water inhibited plant height but promoted root growth. By creating a closed environment, adding a dielectric barrier plate, and combining these two treatment methods, the greatest impact can be had on the growth of Medicago sativa L. seeds. After treatment, the plant heights were approximately half those of the control group, and the root lengths were approximately 1.6 times those of the control group. Through emission spectroscopy, it was found that active particles such as O, NO2, and N2* were present and could be considered to have produced plasma-activated water through contact with the water surface, thus affecting the survival and growth of the seeds. Whole-genome resequencing (WGRS) was performed on the wild-type and selected mutants after treatment, with an average sequencing depth of 115.93×, an average genome alignment rate of 91.72 %, and an average genome coverage rate of 91.85 %. Various types of mutations were detected and annotated. After filtering, 7,822,324 SNP (single nucleotide polymorphisms) sites, 2,161,917 indel sites, 200,544 SV sites, 238 CNV (copy number variation) sites. The SNPs, indels (insertions/deletions), and SVs (structural variations) were mainly heterozygous, with heterozygosity rates of 87.13 %, 92.16 %, and 83.49 %, respectively. The CNVs were dominated by low copy numbers, accounting for 53.77 %. These results indicate that plasma treatment has significant effects on plant growth and genome of Medicago sativa L.

Genome-Wide Identification, Phylogenetic, and Expression Analysis of Jasmonate ZIM-Domain Gene Family in Medicago Sativa L

JASMONATE ZIM domain (JAZ) proteins, inhibitors of the jasmonic acid (JA) signaling pathway, are identified in different plants, such as rice and Arabidopsis. These proteins are crucial for growth, development, and abiotic stress responses. However, limited information is available regarding the JAZ family in alfalfa. This study identified 11 JAZ genes (MsJAZs) in the "Zhongmu No.1" reference genome of alfalfa. The physical and chemical properties, chromosome localization, phylogenetic relationships, gene structure, cis-acting elements, and collinearity of the 11 MsJAZ genes were subsequently analyzed. Tissue-specific analysis revealed distinct functions of different MsJAZ genes in growth and development. The expression patterns of MsJAZ genes under salt stress conditions were validated using qRT-PCR. All MsJAZ genes responded to salt stress, with varying levels of upregulation over time, highlighting their role in stress responses. Furthermore, heterogeneous expression of MsJAZ1 in Arabidopsis resulted in significantly lower seed germination and survival rates in OE-2 and OE-4 compared to the WT under 150 mM NaCl treatment. This study establishes a foundation for further exploration of the function of the JAZ family and provides significant insights into the genetic improvement of alfalfa.

Chromosome-scale and haplotype-resolved genome assembly of the autotetraploid Misgurnus anguillicaudatus

In nature, diploids and tetraploids are two common types of polyploid evolution. Misgurnus anguillicaudatus (mud loach) is a remarkable fish species that exhibits both diploid and tetraploid forms. However, reconstructing the four haplotypes of its autotetraploid genome remains unresolved. Here, we generated the first haplotype-resolved, chromosome-level genome of autotetraploid M. anguillicaudatus with a size of 4.76 Gb, contig N50 of 6.78 Mb, and scaffold N50 of 44.11 Mb. We identified approximately 2.9 Gb (61.03% of genome) of repetitive sequences and predicted 91,485 protein-coding genes. Moreover, allelic gene expression levels indicated the absence of significant dominant haplotypes within the autotetraploid loach genome. This genome will provide a valuable biological model for unraveling the mechanisms of polyploid formation and evolution, adaptation to environmental changes, and benefit for aquaculture applications and biodiversity conservation.

Identification of a Novel Gene MtbZIP60 as a Negative Regulator of Leaf Senescence through Transcriptome Analysis in Medicago truncatula

Leaves are the primary harvest portion in forage crops such as alfalfa (Medicago sativa). Delaying leaf senescence is an effective strategy to improve forage biomass production and quality. In this study, we employed transcriptome sequencing to analyze the transcriptional changes and identify key senescence-associated genes under age-dependent leaf senescence in Medicago truncatula, a legume forage model plant. Through comparing the obtained expression data at different time points, we obtained 1057 differentially expressed genes, with 108 consistently up-regulated genes across leaf growth and senescence. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses showed that the 108 SAGs mainly related to protein processing, nitrogen metabolism, amino acid metabolism, RNA degradation and plant hormone signal transduction. Among the 108 SAGs, seven transcription factors were identified in which a novel bZIP transcription factor MtbZIP60 was proved to inhibit leaf senescence. MtbZIP60 encodes a nuclear-localized protein and possesses transactivation activity. Further study demonstrated MtbZIP60 could associate with MtWRKY40, both of which exhibited an up-regulated expression pattern during leaf senescence, indicating their crucial roles in the regulation of leaf senescence. Our findings help elucidate the molecular mechanisms of leaf senescence in M. truncatula and provide candidates for the genetic improvement of forage crops, with a focus on regulating leaf senescence.

Genome-wide investigation of the TIFY transcription factors in alfalfa (Medicago sativa L.): identification, analysis, and expression

BACKGROUND

Alfalfa (Medicago sativa L.) is an essential leguminous forage with high nutrition and strong adaptability. The TIFY family is a plant-specific transcription factor identified in many plants. However, few reports have been reported on the phylogenetic analysis and gene expression profiling of TIFY family genes in alfalfa.

RESULT

A total of 84 TIFY genes belonging to 4 categories were identified in alfalfa, including 58 MsJAZs, 18 MsZMLs, 4 MsTIFYs and 4 MsPPDs, respectively. qRT-PCR data from 8 genes in different tissues revealed that most MsTIFY genes were highly expressed in roots. The expression of MsTIFY14 was up-regulated after different times in both thrips-resistant and susceptible alfalfa after thrips feeding, and the expression of the remaining MsTIFYs had a strong correlation with the time of thrips feeding. Different abiotic stresses, including drought, salt, and cold, could induce or inhibit the expression of MsTIFY genes to varying degrees. In addition, the eight genes were all significantly up-regulated by JA and/or SA. Interestingly, MsTIFY77 was induced considerably by all the biotic, abiotic, or plant hormones (JA or SA) except ABA.

CONCLUSION

Our study identified members of the TIFY gene family in alfalfa and analyzed their structures and possible functions. It laid the foundation for further research on the molecular functions of TIFYs in alfalfa.

Genome-wide identification of SHMT family genes in alfalfa (Medicago sativa) and its functional analyses under various abiotic stresses

BACKGROUND

Alfalfa (Medicago sativa L.) is the most widely planted legume forage and one of the most economically valuable crops in the world. Serine hydroxymethyltransferase (SHMT), a pyridoxal phosphate-dependent enzyme, plays crucial roles in plant growth, development, and stress responses. To date, there has been no comprehensive bioinformatics investigation conducted on the SHMT genes in M. sativa.

RESULTS

Here, we systematically analyzed the phylogenetic relationship, expansion pattern, gene structure, cis-acting elements, and expression profile of the MsSHMT family genes. The result showed that a total of 15 SHMT members were identified from the M. sativa genome database. Phylogenetic analysis demonstrated that the MsSHMTs can be divided into 4 subgroups and conserved with other plant homologues. Gene structure analysis found that the exons of MsSHMTs ranges from 3 to 15. Analysis of cis-acting elements found that each of the MsSHMT genes contained different kinds of hormones and stress-related cis-acting elements in their promoter regions. Expression and function analysis revealed that MsSHMTs expressed in all plant tissues. qRT-PCR analysis showed that MsSHMTs induced by ABA, Salt, and drought stresses.

CONCLUSIONS

These results provided definite evidence that MsSHMTs might involve in growth, development and adversity responses in M. sativa, which laid a foundation for future functional studies of MsSHMTs.

Transcriptomic analyses provide molecular insight into the cold stress response of cold-tolerant alfalfa

BACKGROUND

Daye No.3 is a novel cultivar of alfalfa (Medicago sativa L.) that is well suited for cultivation in high-altitude regions such as the Qinghai‒Tibet Plateau owing to its high yield and notable cold resistance. However, the limited availability of transcriptomic information has hindered our investigation into the potential mechanisms of cold tolerance in this cultivar. Consequently, we conducted de novo transcriptome assembly to overcome this limitation. Subsequently, we compared the patterns of gene expression in Daye No. 3 during cold acclimatization and exposure to cold stress at various time points.

RESULTS

A total of 15 alfalfa samples were included in the transcriptome assembly, resulting in 141.97 Gb of clean bases. A total of 441 DEGs were induced by cold acclimation, while 4525, 5016, and 8056 DEGs were identified at 12 h, 24 h, and 36 h after prolonged cold stress at 4 °C, respectively. The consistency between the RT‒qPCR and transcriptome data confirmed the accuracy and reliability of the transcriptomic data. KEGG enrichment analysis revealed that many genes related to photosynthesis were enriched under cold stress. STEM analysis demonstrated that genes involved in nitrogen metabolism and the TCA cycle were consistently upregulated under cold stress, while genes associated with photosynthesis, particularly antenna protein genes, were downregulated. PPI network analysis revealed that ubiquitination-related ribosomal proteins act as hub genes in response to cold stress. Additionally, the plant hormone signaling pathway was activated under cold stress, suggesting its vital role in the cold stress response of alfalfa.

CONCLUSIONS

Ubiquitination-related ribosomal proteins induced by cold acclimation play a crucial role in early cold signal transduction. As hub genes, these ubiquitination-related ribosomal proteins regulate a multitude of downstream genes in response to cold stress. The upregulation of genes related to nitrogen metabolism and the TCA cycle and the activation of the plant hormone signaling pathway contribute to the enhanced cold tolerance of alfalfa.

Chromosome-level reference genome and resequencing of 322 accessions reveal evolution, genomic imprint and key agronomic traits in adzuki bean

Adzuki bean (Vigna angularis) is an important legume crop cultivated in over 30 countries worldwide. We developed a high-quality chromosome-level reference genome of adzuki bean cultivar Jingnong6 by combining PacBio Sequel long-read sequencing with short-read and Hi-C technologies. The assembled genome covers 97.8% of the adzuki bean genome with a contig N50 of approximately 16 Mb and a total of 32 738 protein-coding genes. We also generated a comprehensive genome variation map of adzuki bean by whole-genome resequencing (WGRS) of 322 diverse adzuki beans accessions including both wild and cultivated. Furthermore, we have conducted comparative genomics and a genome-wide association study (GWAS) on key agricultural traits to investigate the evolution and domestication. GWAS identified several candidate genes, including VaCycA3;1, VaHB15, VaANR1 and VaBm, that exhibited significant associations with domestication traits. Furthermore, we conducted functional analyses on the roles of VaANR1 and VaBm in regulating seed coat colour. We provided evidence for the highest genetic diversity of wild adzuki (Vigna angularis var. nipponensis) in China with the presence of the most original wild adzuki bean, and the occurrence of domestication process facilitating transition from wild to cultigen. The present study elucidates the genetic basis of adzuki bean domestication traits and provides crucial genomic resources to support future breeding efforts in adzuki bean.

Effects of Selenium Content on Growth, Antioxidant Activity, and Key Selenium-Enriched Gene Expression in Alfalfa Sprouts

To enhance the selenium (Se) intake of the general public, the present study implemented biofortification techniques in alfalfa sprouts. Alfalfa sprouts possess unique nutritional value and provide an optimal Se-enriched supplemental Se source. The impact of sodium selenite (Na2SeO3) on alfalfa shoot germination, shoot length, and biomass was assessed experimentally, and changes in the antioxidant capacity of sprouts treated with optimal Se concentrations were investigated. In addition, the transcriptome of alfalfa sprouts treated with the optimal Na2SeO3 concentration was sequenced. Gene co-expression networks, constructed through differential gene analysis and weighted gene co-expression network analysis, were used to identify the core genes responsible for Se enrichment in alfalfa sprouts. The findings of the present study offer novel insights into the effects of Se treatment on the nutrient composition of alfalfa sprouts, in addition to introducing novel methods and references that could facilitate production of Se-enriched alfalfa sprouts and associated products.

Genetic and epigenetic signatures for improved breeding of cultivated blueberry

Blueberry belongs to the Vaccinium genus and is a highly popular fruit crop with significant economic importance. It was not until the early twentieth century that they began to be domesticated through extensive interspecific hybridization. Here, we collected 220 Vaccinium accessions from various geographical locations, including 154 from the United States, 14 from China, eight from Australia, and 29 from Europe and other countries, comprising 164 Vaccinium corymbosum, 15 Vaccinium ashei, 10 lowbush blueberries, seven half-high blueberries, and others. We present the whole-genome variation map of 220 accessions and reconstructed the hundred-year molecular history of interspecific hybridization of blueberry. We focused on the two major blueberry subgroups, the northern highbush blueberry (NHB) and southern highbush blueberry (SHB) and identified candidate genes that contribute to their distinct traits in climate adaptability and fruit quality. Our analysis unveiled the role of gene introgression from Vaccinium darrowii and V. ashei into SHB in driving the differentiation between SHB and NHB, potentially facilitating SHB's adaptation to subtropical environments. Assisted by genome-wide association studies, our analysis suggested VcTBL44 as a pivotal gene regulator governing fruit firmness in SHB. Additionally, we conducted whole-genome bisulfite sequencing on nine NHB and 12 SHB cultivars, and characterized regions that are differentially methylated between the two subgroups. In particular, we discovered that the β-alanine metabolic pathway genes were enriched for DNA methylation changes. Our study provides high-quality genetic and epigenetic variation maps for blueberry, which offer valuable insights and resources for future blueberry breeding.

Comparative Phenotypic and Transcriptomic Analyses Provide Novel Insights into the Molecular Mechanism of Seed Germination in Response to Low Temperature Stress in Alfalfa

Low temperature is the most common abiotic factor that usually occurs during the seed germination of alfalfa (Medicago sativa L.). However, the potential regulatory mechanisms involved in alfalfa seed germination under low temperature stress are still ambiguous. Therefore, to determine the relevant key genes and pathways, the phenotypic and transcriptomic analyses of low-temperature sensitive (Instict) and low-temperature tolerant (Sardi10) alfalfa were conducted at 6 and 15 h of seed germination under normal (20 °C) and low (10 °C) temperature conditions. Germination phenotypic results showed that Sardi10 had the strongest germination ability under low temperatures, which was manifested by the higher germination-related indicators. Further transcriptome analysis indicated that differentially expressed genes were mainly enriched in galactose metabolism and carbon metabolism pathways, which were the most commonly enriched in two alfalfa genotypes. Additionally, fatty acid metabolism and glutathione metabolism pathways were preferably enriched in Sardi10 alfalfa. The Weighted Gene Co-Expression Network Analysis (WGCNA) suggested that genes were closely related to galactose metabolism, fatty acid metabolism, and glutathione metabolism in Sardi10 alfalfa at the module with the highest correlation (6 h of germination under low temperature). Finally, qRT-PCR analysis further validated the related genes involved in the above pathways, which might play crucial roles in regulating seed germination of alfalfa under low temperature conditions. These findings provide new insights into the molecular mechanisms of seed germination underlying the low temperature stress in alfalfa.

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.

A comparative plastome approach enhances the assessment of genetic variation in the Melilotus genus

BACKGROUND

Melilotus, a member of the Fabaceae family, is a pivotal forage crop that is extensively cultivated in livestock regions globally due to its notable productivity and ability to withstand abiotic stress. However, the genetic attributes of the chloroplast genome and the evolutionary connections among different Melilotus species remain unresolved.

RESULTS

In this study, we compiled the chloroplast genomes of 18 Melilotus species and performed a comprehensive comparative analysis. Through the examination of protein-coding genes, we successfully established a robust phylogenetic tree for these species. This conclusion is further supported by the phylogeny derived from single-nucleotide polymorphisms (SNPs) across the entire chloroplast genome. Notably, our findings revealed that M. infestus, M. siculus, M. sulcatus, and M. speciosus formed a distinct subgroup within the phylogenetic tree. Additionally, the chloroplast genomes of these four species exhibit two shared inversions. Moreover, inverted repeats were observed to have reemerged in six species within the IRLC. The distribution patterns of single-nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) within protein-coding genes indicated that ycf1 and ycf2 accumulated nonconservative alterations during evolutionary development. Furthermore, an examination of the evolutionary rate of protein-coding genes revealed that rps18, rps7, and rpl16 underwent positive selection specifically in Melilotus.

CONCLUSIONS

We present a comparative analysis of the complete chloroplast genomes of Melilotus species. This study represents the most thorough and detailed exploration of the evolution and variability within the genus Melilotus to date. Our study provides valuable chloroplast genomic information for improving phylogenetic reconstructions and making biogeographic inferences about Melilotus and other Papilionoideae species.

Biological nitrogen fixation maintains carbon/nitrogen balance and photosynthesis at elevated CO2

Understanding crop responses to elevated CO2 is necessary to meet increasing agricultural demands. Crops may not achieve maximum potential yields at high CO2 due to photosynthetic downregulation, often associated with nitrogen limitation. Legumes have been proposed to have an advantage at elevated CO2 due to their ability to exchange carbon for nitrogen. Here, the effects of biological nitrogen fixation (BNF) on the physiological and gene expression responses to elevated CO2 were examined at multiple nitrogen levels by comparing alfalfa mutants incapable of nitrogen fixation to wild-type. Elemental analysis revealed a role for BNF in maintaining shoot carbon/nitrogen (C/N) balance under all nitrogen treatments at elevated CO2, whereas the effect of BNF on biomass was only observed at elevated CO2 and the lowest nitrogen dose. Lower photosynthetic rates at were associated with the imbalance in shoot C/N. Genome-wide transcriptional responses were used to identify carbon and nitrogen metabolism genes underlying the traits. Transcription factors important to C/N signalling were identified from inferred regulatory networks. This work supports the hypothesis that maintenance of C/N homoeostasis at elevated CO2 can be achieved in plants capable of BNF and revealed important regulators in the underlying networks including an alfalfa (Golden2-like) GLK ortholog.

The chromosome-level assembly of the wild diploid alfalfa genome provides insights into the full landscape of genomic variations between cultivated and wild alfalfa

Alfalfa (Medicago sativa L.) is one of the most important forage legumes in the world, including autotetraploid (M. sativa ssp. sativa) and diploid alfalfa (M. sativa ssp. caerulea, progenitor of autotetraploid alfalfa). Here, we reported a high-quality genome of ZW0012 (diploid alfalfa, 769 Mb, contig N50 = 5.5 Mb), which was grouped into the Northern group in population structure analysis, suggesting that our genome assembly filled a major gap among the members of M. sativa complex. During polyploidization, large phenotypic differences occurred between diploids and tetraploids, and the genetic information underlying its massive phenotypic variations remains largely unexplored. Extensive structural variations (SVs) were identified between ZW0012 and XinJiangDaYe (an autotetraploid alfalfa with released genome). We identified 71 ZW0012-specific PAV genes and 1296 XinJiangDaYe-specific PAV genes, mainly involved in defence response, cell growth, and photosynthesis. We have verified the positive roles of MsNCR1 (a XinJiangDaYe-specific PAV gene) in nodulation using an Agrobacterium rhizobia-mediated transgenic method. We also demonstrated that MsSKIP23_1 and MsFBL23_1 (two XinJiangDaYe-specific PAV genes) regulated leaf size by transient overexpression and virus-induced gene silencing analysis. Our study provides a high-quality reference genome of an important diploid alfalfa germplasm and a valuable resource of variation landscape between diploid and autotetraploid, which will facilitate the functional gene discovery and molecular-based breeding for the cultivars in the future.

In silico analysis identified bZIP transcription factors genes responsive to abiotic stress in Alfalfa (Medicago sativa L.)

BACKGROUND

Alfalfa (Medicago sativa L.) is the most cultivated forage legume around the world. Under a variety of growing conditions, forage yield in alfalfa is stymied by biotic and abiotic stresses including heat, salt, drought, and disease. Given the sessile nature of plants, they use strategies including, but not limited to, differential gene expression to respond to environmental cues. Transcription factors control the expression of genes that contribute to or enable tolerance and survival during periods of stress. Basic-leucine zipper (bZIP) transcription factors have been demonstrated to play a critical role in regulating plant growth and development as well as mediate the responses to abiotic stress in several species, including Arabidopsis thaliana, Oryza sativa, Lotus japonicus and Medicago truncatula. However, there is little information about bZIP transcription factors in cultivated alfalfa.

RESULT

In the present study, 237 bZIP genes were identified in alfalfa from publicly available sequencing data. Multiple sequence alignments showed the presence of intact bZIP motifs in the identified sequences. Based on previous phylogenetic analyses in A. thaliana, alfalfa bZIPs were similarly divided and fell into 10 groups. The physico-chemical properties, motif analysis and phylogenetic study of the alfalfa bZIPs revealed high specificity within groups. The differential expression of alfalfa bZIPs in a suite of tissues indicates that bZIP genes are specifically expressed at different developmental stages in alfalfa. Similarly, expression analysis in response to ABA, cold, drought and salt stresses, indicates that a subset of bZIP genes are also differentially expressed and likely play a role in abiotic stress signaling and/or tolerance. RT-qPCR analysis on selected genes further verified these differential expression patterns.

CONCLUSIONS

Taken together, this work provides a framework for the future study of bZIPs in alfalfa and presents candidate bZIPs involved in stress-response signaling.

An inducible CRISPR activation tool for accelerating plant regeneration

The inducible CRISPR activation (CRISPR-a) system offers unparalleled precision and versatility for regulating endogenous genes, making it highly sought after in plant research. In this study, we developed a chemically inducible CRISPR-a tool for plants called ER-Tag by combining the LexA-VP16-ER inducible system with the SunTag CRISPR-a system. We systematically compared different induction strategies and achieved high efficiency in target gene activation. We demonstrated that guide RNAs can be multiplexed and pooled for large-scale screening of effective morphogenic genes and gene pairs involved in plant regeneration. Further experiments showed that induced activation of these morphogenic genes can accelerate regeneration and improve regeneration efficiency in both eudicot and monocot plants, including alfalfa, woodland strawberry, and sheepgrass. Our study expands the CRISPR toolset in plants and provides a powerful new strategy for studying gene function when constitutive expression is not feasible or ideal.

MrERF039 transcription factor plays an active role in the cold response of Medicago ruthenica as a sugar molecular switch

Cold stress severely restricts plant development, causing significant agricultural losses. We found a critical transcription factor network in Medicago ruthenica was involved in plant adaptation to low-temperature. APETALA2/ethylene responsive factor (AP2/ERF) transcription factor MrERF039 was transcriptionally induced by cold stress in M. ruthenica. Overexpression of MrERF039 significantly increased the glucose and maltose content, thereby improving the tolerance of M. ruthenica. MrERF039 could bind to the DRE cis-acting element in the MrCAS15A promoter. Additionally, the methyl group of the 14th amino acid in MrERF039 was required for binding. Transcriptome analysis showed that MrERF039 acted as a sugar molecular switch, regulating numerous sugar transporters and sugar metabolism-related genes. In addition, we found that MrERF039 could directly regulate β-amylase gene, UDP glycosyltransferase gene, and C2H2 zinc finger protein gene expression. In conclusion, these findings suggest that high expression of MrERF039 can significantly improve the cold tolerance of M. ruthenica root tissues during cold acclimation. Our results provide a new theoretical basis and candidate genes for breeding new legume forage varieties with high resistance.

Genome organization and botanical diversity

The rich diversity of angiosperms, both the planet's dominant flora and the cornerstone of agriculture, is integrally intertwined with a distinctive evolutionary history. Here, we explore the interplay between angiosperm genome organization and botanical diversity, empowered by genomic approaches ranging from genetic linkage mapping to analysis of gene regulation. Commonality in the genetic hardware of plants has enabled robust comparative genomics that has provided a broad picture of angiosperm evolution and implicated both general processes and specific elements in contributing to botanical diversity. We argue that the hardware of plant genomes-both in content and in dynamics-has been shaped by selection for rather substantial differences in gene regulation between plants and animals such as maize and human, organisms of comparable genome size and gene number. Their distinctive genome content and dynamics may reflect in part the indeterminate development of plants that puts strikingly different demands on gene regulation than in animals. Repeated polyploidization of plant genomes and multiplication of individual genes together with extensive rearrangement and differential retention provide rich raw material for selection of morphological and/or physiological variations conferring fitness in specific niches, whether natural or artificial. These findings exemplify the burgeoning information available to employ in increasing knowledge of plant biology and in modifying selected plants to better meet human needs.