Characterization of the microRNA408-LACCASE5 module as a regulatory axis for photosynthetic efficiency in Medicago ruthenica: implications for forage yield enhancement.

Medicago ruthenica is closely related to Medicago sativa, a commonly cultivated forage. Characterized by its high tolerance to environmental stress, M. ruthenica is a valuable genetic resource. However, low yield limits its large-scale utilization. Leaf morphology, an important agronomic trait, is closely related to forage yield and photosynthetic efficiency. In the presented study, "Correlation of Leaf Morphology and Photosynthetic Performance with Forage Yield in Medicago ruthenica: The Underlying Molecular Mechanisms," comprehensive data analysis revealed a significant positive association between leaf width and leaf area with forage yield in Medicago ruthenica (p < 0.05). The specific cultivar "Mengnong No.1 (MN No.1) had a large leaf area, and its physiological parameters related to photosynthetic characteristics were superior. Anatomical examination revealed that the leaves of MN No.1 had strong palisade tissue and compact cell structure. Subsequent investigations, utilizing small RNA and transcriptome sequencing, discerned critical miRNA-target gene networks that underpin the high photosynthetic efficiency in M. ruthenica. A total of 63 differentially expressed miRNAs (DEMs) were identified, inclusive of several well-characterized miRNAs such as miR408, miR171, and miR398. These miRNAs were predicted to target 55 genes (mRNAs), of which 6 miRNA-target gene pairs, particularly those involving miR408and miR171, exhibited inverse expression patterns. Among the six postulated miRNA-target gene pairs, the targeted cleavage of LACCASE5 (LAC5) by miR408 was conclusively validated through degradome sequencing, with the cleavage site pinpointed between the 9th and 10th nucleotides from the 5'end of miR408 via the 5'-RLM-RACE assay. Therefore, it is posited that the miR408-MrLAC5 module constitutes a central mechanism in fostering high photosynthetic efficiency in M. ruthenica. Moreover, these findings also provide valuable information for further study of the regulatory genes and miRNA functions of forage yield in legume forage.

Leaf transcriptome analysis of Medicago ruthenica revealed its response and adaptive strategy to drought and drought recovery.

BACKGROUND: Drought is one of the main causes of losses in forage crop yield and animal production. Medicago ruthenica (L.) cv. Zhilixing is a high-yielding alfalfa cultivar also known for its high tolerance to drought. We analyzed the transcriptome profile of this cultivar throughout drought stress and recovery and we were able to describe its phased response through the expression profiles of overlapping gene networks and drought-specific genes. RESULTS: The ABA and auxin signal transduction pathways are overlapping pathways in response to drought and drought recovery in forage crops. Medicago ruthenica (L.) cv. Zhilixing adopts different strategies at different degrees of drought stress. On the 9th day of drought, transcriptional regulations related to osmoregulation are enhanced mainly through increased activities of carbohydrate and amino acid metabolism, while photosynthetic activities were reduced to slow down growth. With drought prolonging, on the 12th day of drought, the synthesis of proline and other stored organic substances was suppressed in general. After recovery, Medicago ruthenica synthesizes flavonoids through the flavonoid biosynthesis pathway to remove accumulated ROS and repair the oxidative damage from water stress. In addition, the regulation of circadian rhythm seems to accelerate the drought recovery process. CONCLUSIONS: Medicago ruthenica adapts to drought by regulating the osmoregulatory system and photosynthesis, which appears to involve the ABA and auxin signaling pathways as key regulators. Furthermore, the synthesis of flavonoids and the regulation of the circadian rhythm can accelerate the recovery process. These results enriched our knowledge of molecular responses to drought and drought recovery in Medicago ruthenica and provide useful information for the development of new legume forage grass varieties with improved adaptability to drought stress.

The MicroRNA397a-LACCASE17 module regulates lignin biosynthesis in Medicago ruthenica (L.).

Mechanical strength is essential for the upright growth habit, which is one of the most important characteristics of terrestrial plants. Lignin, a phenylpropanoid-derived polymer mainly present in secondary cell walls plays critical role in providing mechanical support. Here, we report that the prostrate-stem cultivar of the legume forage Medicago ruthenica cultivar 'Mengnong No. 1' shows compromised mechanical strength compared with the erect-stem cultivar 'Zhilixing'. The erect-stem cultivar, 'Zhilixing' has significantly higher lignin content, leading to higher mechanical strength than the prostrate-stem cultivar. The low abundance of miRNA397a in the Zhiixing cultivar causes reduced cleavage of MrLAC17 transcript, which results in enhanced expression level of MrLAC17 compared to that in the prostrate-stem cultivar Mengnong No. 1. Complementation of the Arabidopsis lac4 lac17 double mutants with MrLAC17 restored the lignin content to wild-type levels, confirming that MrLAC17 perform an exchangeable role with Arabidopsis laccases. LAC17-mediated lignin polymerization is therefore increased in the 'Zhilixing', causing the erect stem phenotype. Our data reveal the importance of the miR397a in the lignin biosynthesis and suggest a strategy for molecular breeding targeting plant architecture in legume forage.

Targeted Metabolomics Provide Chemotaxonomic Insights of Medicago ruthenica, with Coupled Transcriptomics Elucidating the Mechanism Underlying Floral Coloration.

Medicago ruthenica, a wild legume forage widely distributed in the Eurasian steppe, demonstrates high genetic and phenotypic variation. M. ruthenica with a purely yellow flower (YFM), differing from the general phenotype of M. ruthenica with a purple flower (PFM), was recently discovered. The similar characteristics of YFM with Medicago falcata have led to conflicting opinions on its taxonomy using traditional morphological methods. The lack of chemotaxonomy information about M. ruthenica species and the unclear flower coloration mechanisms have hampered their study. Here, we investigated M. ruthenica using targeted metabolomics based on the chemotaxonomy method and elaborated the floral coloration mechanisms using transcriptomics. The identified flavonoids were the same types, but there were different contents in YFM and PFM, especially the contents of cyanidin-3-O-glucoside (C3G), an anthocyanin that causes the purple-reddish color of flowers. The over-accumulation of C3G in PFM was 1,770 times more than YFM. Nineteen anthocyanin-related genes were downregulated in YFM compared with their expression in PFM. Thus, YFM could be defined as a variety of M. ruthenica rather than a different species. The loss of purple flower coloration in YFM was attributed to the downregulation of these genes, resulting in reduced C3G accumulation. The taxonomic characteristics and molecular and physiological characteristics of this species will contribute to further research on other species with similar external morphologies.

MrERF, MrbZIP, and MrSURNod of Medicago ruthenica Are Involved in Plant Growth and Abiotic Stress Response.

Abiotic stresses affect plant growth and productivity. The outstanding stress resistance of Medicago ruthenica makes it a desirable gene resource to improve the stress tolerance of other plants. The roles of three differently expressed genes [(DEGs) (MrERF, MrbZIP, and MrSURNod)] from M. ruthenica in stress resistance have not been fully elucidated. Therefore, we constructed their expression vectors, transformed them into tobacco, and subjected transgenic lines to abiotic stresses. Through comprehensive bioinformatics, transcriptomic, morphological, and physiological analyses of transgenic lines, we have revealed the critical role of these three DEGs in plant growth and abiotic stress response. The upregulation of genes enhanced the germination rate, biomass, root length number, etc. Additionally, the accumulation of osmolytes increased the activity of antioxidant enzymes. These genes are also associated with improved seed yield, increased branching, and early flowering, thereby shortening the growth period. Potentially, this is one of the ways for tobacco to cope with stress. Furthermore, the resistance of transgenic tobacco expressing MrERF or MrbZIP was better than that with MrSURNod. MrERF and MrbZIP can improve drought and salt tolerance of plants, whereas MrSURNod is beneficial in improving drought and cold resistance. Moreover, MrERF or MrbZIP can promote root elongation and increase the root number, whereas MrSURNod mainly promotes root elongation. This may be the reason why stress resistance conferred by MrSURNod is weaker than that associated with the other two genes. Overall, MrERF, MrbZIP, and MrSURNod positively modulate plant growth and stress tolerance.

Genome-wide identification and analysis of the NLR gene family in Medicago ruthenica.

Medicago ruthenica, important forage in the legume family, possesses high nutritional value and carries abundant tolerance genes. This study used whole-genome data of M. ruthenica to perform a genome-wide analysis of the nucleotide-binding site-leucine-rich repeat receptor (NLR) gene family, which is the largest family of plant disease resistance genes (R genes). A total of 338 NLR genes were identified in the M. ruthenica genome, including 160 typical genes that contained 80 coiled-coil (CC)-NBS-LRR (CNL) genes, 76 toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) genes, four resistance to powdery mildew 8 (RPW8)-NBS-LRR (RNL) subclass genes, and 178 atypical NLR genes encoding proteins without at least one important domain. Among its eight chromosomes, M. ruthenica chromosomes 3 and 8 contained most of the NLR genes. More than 40% of all NLR genes were located on these two chromosomes, mainly in multigene clusters. The NLR proteins of M. ruthenica had six highly conserved motifs: P-loop, GLPL, RNBS-D, kinase-2, RNBS-C, and MHDV. Phylogenetic analysis revealed that the NLR genes of M. ruthenica formed three deeply separated clades according to the N-terminal domain of the proteins encoded by these genes. Gene duplication and syntenic analysis suggested four gene duplication types in the NLR genes of M. ruthenica, namely, tandem, proximal, dispersed, and segmental duplicates, which involved 189, 49, 59, and 41 genes, respectively. A total of 41 segmental duplication genes formed 23 NLR gene pairs located on syntenic chromosomal blocks mainly between chromosomes 6 and 7. In addition, syntenic analysis between M. truncatula and M. ruthenica revealed 193 gene pairs located on syntenic chromosomal blocks of the two species. The expression analysis of M. ruthenica NLR genes showed that 303 (89.6%) of the NLR genes were expressed in different varieties. Overall, this study described the full NLR profile of the M. ruthenica genome to provide an important resource for mining disease-resistant genes and disease-resistant breeding.

Genome-wide identification of microRNAs associated with osmotic stress and elucidation of the role of miR319 in Medicago ruthenica seedlings.

Drought is a major environmental stress that affects plant growth, development, and productivity. Medicago ruthenica, a leguminous forage, has garnered attention owing to its resistance to abiotic stress. The purpose of the current study was to explore genes conferring drought resistance to M. ruthenica. MicroRNAs (miRNAs) are post-transcriptional regulators of gene expression in plants and are associated with developmental plasticity and abiotic/biotic stress responses. Here, high-throughput small RNA, mRNA, and degradome sequencing analyses were performed to analyze miRNAs and their potential target genes in the leaves of M. ruthenica seedlings under osmotic stress conditions. In total, 591 miRNAs were identified. A comparison of the expression levels showed that 15 miRNAs (14 upregulated and 1 downregulated) were significantly differentially expressed following PEG6000 treatment compared with those in the control (0 h). Most miRNAs are highly conserved between M. ruthenica and Medicago truncatula. Using TargetFinder, 11 target genes were predicted; the expression of these target genes negatively correlated with that of five miRNAs related to osmotic stress response. miR319 downregulated the expression of teosinte branched/cycloidea/proliferating cell factor 4 (TCP4), which encodes plant-specific transcription factors, more significantly in the leaves than in the roots. These results were confirmed using quantitative real-time polymerase chain reaction, northern blotting, RLM 5'RACE, and a Nicotiana benthamiana transient expression system. The miR319-TCP4 module may act as a homeostasis factor in M. ruthenica roots following drought injury, and it is conserved among plant species.

The complete chloroplast genome of Medicago ruthenica cv. 'Taihang' (Fabaceae) from Shanxi, China.

Medicago ruthenica is an important perennial forage with multiple characteristics of resistance. In this study, we sequenced and characterized the complete chloroplast genome of M. ruthenica 'Taihang', which is 124, 254 bp in length. A total of 108 genes were identified, including 74 protein-coding, 30 tRNA, and four rRNA genes. Phylogenetic analysis based on 27 chloroplast genomes showed that M. ruthenica 'Taihang' has a close relationship with M. ruthenica from Qinghai Province, China. The data are useful in better understanding the genetic diversity and stress resistance of Medicago and contribute to the phylogenetic study of Trifolieae.

Utilization of Transcriptome, Small RNA, and Degradome Sequencing to Provide Insights Into Drought Stress and Rewatering Treatment in Medicago ruthenica.

Drought is a major limiting factor in foraging grass yield and quality. Medicago ruthenica (M. ruthenica) is a high-quality forage legume with drought resistance, cold tolerance, and strong adaptability. In this study, we integrated transcriptome, small RNA, and degradome sequencing in identifying drought response genes, microRNAs (miRNAs), and key miRNA-target pairs in M. ruthenica under drought and rewatering treatment conditions. A total of 3,905 genes and 50 miRNAs (45 conserved and 5 novel miRNAs) were significantly differentially expressed in three test conditions (CK: control, DS: plants under drought stress, and RW: plants rewatering after drought stress). The degradome sequencing (AllenScore < 4) analysis revealed that 104 miRNAs (11 novel and 93 conserved miRNAs) were identified with 263 target transcripts, forming 296 miRNA-target pairs in three libraries. There were 38 differentially expressed targets from 16 miRNAs in DS vs. CK, 31 from 11 miRNAs in DS vs. RW, and 6 from 3 miRNAs in RW vs. CK; 21, 18, and 3 miRNA-target gene pairs showed reverse expression patterns in DS vs. CK, DS vs. RW, and RW vs. CK comparison groups, respectively. These findings provide valuable information for further functional characterization of genes and miRNAs in response to abiotic stress, in general, and drought stress in M. ruthenica, and potentially contribute to drought resistance breeding of forage in the future.

Environment-Driven Adaptations of Leaf Cuticular Waxes Are Inheritable for Medicago ruthenica.

Cuticular waxes covering the plant surface play pivotal roles in helping plants adapt to changing environments. However, it is still not clear whether the responses of plant cuticular waxes to their growing environments are inheritable. We collected seeds of Medicago ruthenica (a perennial legume) populations from 30 growing sites in northern China and examined the variations of leaf cuticular waxes in a common garden experiment. Four wax genes, MrFAR3-1, MrFAR3-2, MrCER1, and MrKCS1, involved in biosynthesis of predominant wax classes (primary alcohol and alkane) and wax precursors, were isolated to test the contributions of genetic variations of the coding sequences (CDS) and the promoter sequences and epigenetic modifications. The plasticity responses of the cuticular waxes were further validated by two stress-modeling experiments (drought and enhancing ultraviolet B). Great variations in total wax coverage and abundance of wax classes or wax compounds were observed among M. ruthenica populations in a common garden experiment. Stress-modeling experiments further validated that M. ruthenica would alter leaf wax depositions under changed growing conditions. The transcriptional levels of the wax genes were positively or negatively correlated with amounts of cuticular waxes. However, the analysis of promoter methylation showed that the methylation level of the promoter region was not associated with their expressions. Although both promoter sequences and CDS showed a number of polymorphic sites, the promoters were not naturally selected and insignificant difference could be observed in the numbers and types of acting elements of the four wax genes among populations. In contrast, the CDS of the wax genes were naturally selected, with a number of missense mutations resulting in alterations of the amino acid as well as their isoelectric points and polarities, which could impact on enzyme function/activity. We conclude that long-term adaptation under certain environments would induce genetic mutation of wax biosynthesis genes, resulting in inheritable alterations of cuticular wax depositions.

The genome of a wild Medicago species provides insights into the tolerant mechanisms of legume forage to environmental stress.

BACKGROUND: Medicago ruthenica, a wild and perennial legume forage widely distributed in semi-arid grasslands, is distinguished by its outstanding tolerance to environmental stress. It is a close relative of commonly cultivated forage of alfalfa (Medicago sativa). The high tolerance of M. ruthenica to environmental stress makes this species a valuable genetic resource for understanding and improving traits associated with tolerance to harsh environments. RESULTS: We sequenced and assembled genome of M. ruthenica using an integrated approach, including PacBio, Illumina, 10×Genomics, and Hi-C. The assembled genome was 904.13 Mb with scaffold N50 of 99.39 Mb, and 50,162 protein-coding genes were annotated. Comparative genomics and transcriptomic analyses were used to elucidate mechanisms underlying its tolerance to environmental stress. The expanded FHY3/FAR1 family was identified to be involved in tolerance of M. ruthenica to drought stress. Many genes involved in tolerance to abiotic stress were retained in M. ruthenica compared to other cultivated Medicago species. Hundreds of candidate genes associated with drought tolerance were identified by analyzing variations in single nucleotide polymorphism using accessions of M. ruthenica with varying tolerance to drought. Transcriptomic data demonstrated the involvements of genes related to transcriptional regulation, stress response, and metabolic regulation in tolerance of M. ruthenica. CONCLUSIONS: We present a high-quality genome assembly and identification of drought-related genes in the wild species of M. ruthenica, providing a valuable resource for genomic studies on perennial legume forages.

Genomic analysis of Medicago ruthenica provides insights into its tolerance to abiotic stress and demographic history.

Medicago ruthenica has been recently cultivated as a new forage crop and has been recognized as a source of genes to improve abiotic stress tolerance in cultivated alfalfa because of its remarkable tolerance to drought, salinity-alkalinity, and cold and snowy winters. Here, we reveal a chromosome-scale genome sequence of M. ruthenica based on Illumina, PacBio, and Hi-C data. The assembled genome consists of 903.56 Mb with 50,268 annotated protein-coding genes, which is larger and contains relatively more genes than Medicago truncatula (420 Mb and 44,623 genes) and Medicago sativa spp. caerulea (793 Mb and 47,202 genes). All three species shared the ancestral Papilionoideae whole-genome duplication event before their divergence. The more recent expansion of repetitive elements compared to that in the other two species was determined to have contributed greatly to the larger genome size of M. ruthenica. We further found that multiple gene and transcription factor families (e.g., SOS homologous genes, NAC, C2H2, and CAMTA) have expanded in M. ruthenica, which might have led to its enhanced tolerance to abiotic stress. In addition, M. ruthenica harbors more genes involved in the lignin and cellulose biosynthesis pathways than the other two species. Finally, population genomic analyses revealed two genetic lineages, reflecting the west and east of its geographical distribution, respectively. The two lineages probably diverged during the last glaciation and survived in multiple refugia at the last glacial maximum, followed by recent expansion. Our genomic data provide a genetic basis for further molecular breeding research on M. ruthenica and alfalfa.

Complete chloroplast genome of a high-quality forage in north China, Medicago ruthenica (Fabaceae：Trifolieae).

Medicago ruthenica is a well-known high-quality forage due to its good palatability and strong tolerance to drought, cold and saline-alkali stress. Here, the complete chloroplast genome sequence of M. ruthenica was reported. The chloroplast genome is 126,939 bp in length. This chloroplast genome has no inverted repeat (IR) regions, which is very common in the family Fabaceae. The M. ruthenica chloroplast genome encodes 107 genes, including 73 protein-coding genes, 30 tRNA genes, and 4 rRNA genes. Phylogenetic analysis result strongly suggested that M. ruthenica is a distinct lineage in Medicago, being sister to highly supported clade composed of three species (M. hybrida, M. papillosa and M. sativa).