Search

Background : Drought is currently a global environmental problem, which inhibits plant growth and development and seriously restricts crop yields. Many plants exposed to drought stress can generate stress memory, which provides some advantages for resisting recurrent drought. DNA methylation is a mechanism involved in stress memory formation, and many plants can alter methylation levels to form stress memories; however, it remains unclear whether Medicago ruthenica exhibits drought stress memory, as the epigenetic molecular mechanisms underlying this process have not been described in this species. Methods : We conducted methylome and transcriptome sequencing to identify gene methylation and expression changes in plants with a history of two drought stress exposures. Results : Methylation analysis showed that drought stress resulted in an approximately 4.41% decrease in M. ruthenica genome methylation levels. The highest methylation levels were in CG dinucleotide contexts, followed by CHG contexts, with CHH contexts having the lowest levels. Analysis of associations between methylation and transcript levels showed that most DNA methylation was negatively correlated with gene expression except methylation within CHH motifs in gene promoter regions. Genes were divided into four categories according to the relationship between methylation and gene expression; the up-regulation of hypo-methylated gene expression accounted for the vast majority (692 genes) and included genes encoding factors key for abscisic acid (ABA) and proline synthesis. The hypo-methylation of the promoter and body regions of these two gene groups induced increased gene transcription levels. Conclusions : In conclusion, DNA methylation may contribute to drought stress memory formation and maintenance in M. ruthenica by increasing the transcription levels of genes key for ABA and proline biosynthesis.

Alfalfa species Medicago sativa L. (MS) and Medicago falcata L. (MF), globally prominent perennial leguminous forages, hold substantial economic value. However, our comprehension of the molecular mechanisms governing their resistance to cold stress remains limited. To address this knowledge gap, we scrutinized and compared MS and MF cold-stress responses at the molecular level following 24 h and 120 h low-temperature exposure (4 °C). Our study revealed that MF had superior physiological resilience to cold stress compared with MS, and its morphology was healthier under cold stress, and its malondialdehyde content and superoxide dismutase activity increased, first, and then decreased, while the soluble sugar content continued to accumulate. Transcriptome analysis showed that after 120 h of exposure, there were different gene-expression patterns between MS and MF, including 1274 and 2983 genes that were continuously up-regulated, respectively, and a total of 923 genes were included, including star cold-resistant genes such as ICE1 and SIP1 . Gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses revealed numerous inter-species differences in sustained cold-stress responses. Notably, MS-exclusive genes included a single transcription factor (TF) gene and several genes associated with a single DNA repair-related pathway, whereas MF-exclusive genes comprised nine TF genes and genes associated with 14 pathways. Both species exhibited high-level expression of genes encoding TFs belonging to AP2-EREBP, ARR-B, and bHLH TF families, indicating their potential roles in sustaining cold resistance in alfalfa-related species. These findings provide insights into the molecular mechanisms governing cold-stress responses in MS and MF, which could inform breeding programs aimed at enhancing cold-stress resistance in alfalfa cultivars.

Phosphate starvation response (PHR) transcription factors play essential roles in regulating phosphate uptake in plants through binding to the P1BS cis-element in the promoter of phosphate starvation response genes. Recently, PHRs were also shown to positively regulate arbuscular mycorrhizal colonization in rice and lotus by controlling the expression of many symbiotic genes. However, their role in arbuscule development has remained unclear. In Medicago, we previously showed that arbuscule degradation is controlled by two SPX proteins that are highly expressed in arbuscule-containing cells. Since SPX proteins bind to PHRs and repress their activity in a phosphate-dependent manner, we investigated whether arbuscule maintenance is also regulated by PHR. Here, we show that PHR2 is a major regulator of the phosphate starvation response in Medicago. Knockout of phr2 showed reduced phosphate starvation response, symbiotic gene expression, and fungal colonization levels. However, the arbuscules that formed showed less degradation, suggesting a negative role for PHR2 in arbuscule maintenance. This was supported by the observation that overexpression of PHR2 led to enhanced degradation of arbuscules. Although many arbuscule-induced genes contain P1BS elements in their promoters, we found that the P1BS cis-elements in the promoter of the symbiotic phosphate transporter PT4 are not required for arbuscule-containing cell expression. Since both PHR2 and SPX1/3 negatively affect arbuscule maintenance, our results indicate that they control arbuscule maintenance partly via different mechanisms. While PHR2 potentiates symbiotic gene expression and colonization, its activity in arbuscule-containing cells needs to be tightly controlled to maintain a successful symbiosis in Medicago., (© 2024 The Author(s). New Phytologist © 2024 New Phytologist Foundation.)

A perennial leguminous forage, Medicago ruthenica has outstanding tolerance to abiotic stresses. The genome of Medicago ruthenica is large and has a complex genetic background, making it challenging to accurately determine genetic information. However, the chloroplast genome is widely used for researching issues related to evolution, genetic diversity, and other studies. To better understand its chloroplast characteristics and adaptive evolution, chloroplast genomes of 61 Medicago ruthenica were assembled (including 16 cultivated Medicago ruthenica germplasm and 45 wild Medicago ruthenica germplasm). These were used to construct the pan-chloroplast genome of Medicago ruthenica , and the chloroplast genomes of cultivated and wild Medicago ruthenica were compared and analyzed. Phylogenetic and haplotype analyses revealed two main clades of 61 Medicago ruthenica germplasm chloroplast genomes, distributed in eastern and western regions. Meanwhile, based on chloroplast variation information, 61 Medicago ruthenica germplasm can be divided into three genetic groups. Unlike the phylogenetic tree constructed from the chloroplast genome, a new intermediate group has been identified, mainly consisting of samples from the eastern region of Inner Mongolia, Shanxi Province, and Hebei Province. Transcriptomic analysis showed that 29 genes were upregulated and three genes were downregulated. The analysis of these genes mainly focuses on enhancing plant resilience and adapting adversity by stabilizing the photosystem structure and promoting protein synthesis. Additionally, in the analysis of adaptive evolution, the accD , clpP and ycf1 genes showed higher average Ka/Ks ratios and exhibited significant nucleotide diversity, indicating that these genes are strongly positively selected. The editing efficiency of the ycf1 and clpP genes significantly increases under abiotic stress, which may positively contribute to plant adaptation to the environment. In conclusion, the construction and comparative analysis of the complete chloroplast genomes of 61 Medicago ruthenica germplasm from different regions not only revealed new insights into the genetic variation and phylogenetic relationships of Medicago ruthenica germplasm, but also highlighted the importance of chloroplast transcriptome analysis in elucidating the model of chloroplast responses to abiotic stress. These provide valuable information for further research on the adaptive evolution of Medicago ruthenica .

Alfalfa is widely recognized as an important forage crop. To understand the morphological characteristics and genetic basis of seed morphology in alfalfa, we screened 318 Medicago spp., including 244 Medicago sativa subsp. sativa (alfalfa) and 23 other Medicago spp., for seed area size, length, width, length-to-width ratio, perimeter, circularity, the distance between the intersection of length & width (IS) and center of gravity (CG), and seed darkness & red-green-blue (RGB) intensities. The results revealed phenotypic diversity and correlations among the tested accessions. Based on the phenotypic data of M. sativa subsp. sativa, a genome-wide association study (GWAS) was conducted using single nucleotide polymorphisms (SNPs) called against the Medicago truncatula genome. Genes in proximity to associated markers were detected, including CPR1, MON1, a PPR protein, and Wun1(threshold of 1E-04). Machine learning models were utilized to validate GWAS, and identify additional marker-trait associations for potentially complex traits. Marker S7_33375673, upstream of Wun1, was the most important predictor variable for red color intensity and highly important for brightness. Fifty-two markers were identified in coding regions. Along with strong correlations observed between seed morphology traits, these genes will facilitate the process of understanding the genetic basis of seed morphology in Medicago spp., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Aluminum (Al) is the major limiting factor affecting plant productivity in acidic soils. Al 3+ ions exhibit increased solubility at a pH below 5, leading to plant root tip toxicity. Alternatively, plants can perceive very low concentrations of Al 3+ , and Al triggers downstream signaling even at pH 5.7 without causing Al toxicity. The ALUMINUM-ACTIVATED-MALATE-TRANSPORTER (ALMT) family members act as anion channels, with some regulating the secretion of malate from root apices to chelate Al, which is a crucial mechanism for plant Al resistance. To date, the role of the ALMT gene family within the legume Medicago species has not been fully characterized. In this study, we investigated the ALMT gene family in M. sativa and M. truncatula and identified 68 MsALMTs and 18 MtALMTs, respectively. Phylogenetic analysis classified these genes into five clades, and synteny analysis uncovered genuine paralogs and orthologs. The real-time quantitative reverse transcription PCR (qRT-PCR) analysis revealed that MtALMT8, MtALMT9, and MtALMT15 in clade 2-2b are expressed in both roots and root nodules, and MtALMT8 and MtALMT9 are significantly upregulated by Al in root tips. We also observed that MtALMT8 and MtALMT9 can partially restore the Al sensitivity of Atalmt1 in Arabidopsis. Moreover, transcriptome analysis examined the expression patterns of these genes in M. sativa in response to Al at both pH 5.7 and pH 4.6, as well as to protons, and found that Al and protons can independently induce some Al-resistance genes. Overall, our findings indicate that MtALMT8 and MtALMT9 may play a role in Al resistance, and highlight the resemblance between the ALMT genes in Medicago species and those in Arabidopsis., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier GmbH. All rights reserved.)

BBX protein is a class of zinc finger transcription factors that have B-box domains at the N-terminus, and some of these proteins contain a CCT domain at the C-terminus. It plays an important role in plant growth, development, and metabolism. However, the expression pattern of BBX genes in alfalfa under hormonal and salt stresses is still unclear. In this study, we identified a total of 125 BBX gene family members by the available Medicago reference genome in diploid alfalfa ( Medicago sativa spp. Caerulea ), a model plant ( M. truncatula ), and tetraploid alfalfa ( M . sativa ), and divided these members into five subfamilies. We found that the conserved motifs of BBXs of the same subfamily reveal similarities. We analyzed the collinearity relationship and duplication mode of these BBX genes and found that the expression pattern of BBX genes is specific in different tissues. Analysis of the available transcriptome data suggests that some members of the BBX gene family are involved in multiple abiotic stress responses, and the highly expressed genes are often clustered together. Furthermore, we identified different expression patterns of some BBX genes under salt, ethylene, salt and ethylene, salicylic acid, and salt and salicylic acid treatments, verified by qRT-PCR, and analyzed the subcellular localization of MsBBX2 , MsBBX17 , and MsBBX32 using transient expression in tobacco. The results showed that BBX genes were localized in the nucleus. This study systematically analyzed the BBX gene family in Medicago plants, which provides a basis for the study of BBX gene family tolerance to abiotic stresses.

Acyl-CoA-binding proteins (ACBPs) are mainly involved in acyl-CoA ester binding and trafficking in eukaryotic cells, and they function in lipid metabolism, membrane biosynthesis, cellular signaling, stress response, disease resistance, and other biological activities in plants. However, the roles of ACBP family members in Medicago remain unclear. In this study, a total of eight ACBP genes were identified in the genome of Medicago truncatula and Medicago sativa, and they were clustered into four sub-families (Class I-IV). Many cis-acting elements related to abiotic response were identified in the promoter region of these ACBP genes, in particular light-responsive elements. These ACBP genes exhibited distinct expression pattern in various tissues, and the expression level of MtACBP1/MsACBP1 and MtACBP2/MsACBP2 gene pairs were significantly increased under NaCl treatment. Subcellular localization analysis showed that MtACBP1/MsACBP1 and MtACBP2/MsACBP2 were localized in the endoplasmic reticulum of tobacco epidermal cells. Arabidopsis seedlings over-expressing MtACBP2/MsACBP2 displayed increased root length than the wild type under short light, Cu 2+ , ABA, PEG, and NaCl treatments. Over-expression of MtACBP2/MsACBP2 also significantly enhanced Arabidopsis tolerance under NaCl and PEG treatments in mature plants. Collectively, our study identified salt and drought responsive ACBP genes in Medicago and verified their functions in increasing resistance against salt and drought stresses., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024. Published by Elsevier B.V.)

Calcium-dependent protein kinases (CDPKs) are an important class of calcium-sensitive response proteins that play an important regulatory role in response to abiotic stresses. To date, little is known about the CDPK genes in white clover. White clover is a high-quality forage grass with high protein content, but it is susceptible to cold stress. Therefore, we performed a genome-wide analysis of the CDPK gene family in white clover and identified 50 members of the CDPK genes. Phylogenetic analysis using CDPKs from the model plant Arabidopsis divided the TrCDPK genes into four groups based on their sequence similarities. Motif analysis showed that TrCDPKs within the same group had similar motif compositions. Gene duplication analysis revealed the evolution and expansion of TrCDPK genes in white clover. Meanwhile, a genetic regulatory network (GRN) containing TrCDPK genes was reconstructed, and gene ontology (GO) annotation analysis of these functional genes showed that they contribute to signal transduction, cellular response to stimuli, and biological regulation, all of which are important processes in response to abiotic stresses. To determine the function of TrCDPK genes, we analyzed the RNA-seq dataset and found that most TrCDPK genes were highly up-regulated under cold stress, particularly in the early stages of cold stress. These results were validated by qRT-PCR experiments, implying that TrCDPK genes are involved in various gene regulatory pathways in response to cold stress. Our study may help to further investigate the function of TrCDPK genes and their role in response to cold stress, which is important for understanding the molecular mechanisms of cold tolerance in white clover and improving its cold tolerance.

Plant mitochondria are crucial for species evolution, phylogenetics, classification, and identification as maternal genetic material. However, the presence of numerous repetitive sequences, complex structures, and a low number of genes in the mitochondrial genome has hindered its complete assembly and related research endeavors. In this study, we assembled two mitochondrial genomes of alfalfa varieties of Zhongmu No.1 (299,123 bp) and Zhongmu No.4 (306,983 bp), based on a combination of PacBio, Illumina, and Hi-C sequences. The comparison of genome assemblies revealed that the same number of mitochondrial genes, including thirty-three protein-coding genes, sixteen tRNA genes, and three rRNA genes existed in the two varieties. Additionally, large fragments of repetitive sequences were found underlying frequent mitochondrial recombination events. We observed extensive transfer of mitochondrial fragments into the nuclear genome of Zhongmu No.4. Analysis of the cox1 and rrn18s genes in 35 Medicago accessions revealed the presence of population-level deletions and substitutions in the rrn18s gene. We propose that mitochondrial structural reorganizations may contribute to alfalfa evolution.

Lateral roots are crucial for plant growth and development, making them an important target for research aiming to improve crop yields and food security. However, their endogenous ontogeny and, as it were, stochastic appearance challenge their study. Lateral Root Inducible Systems (LRIS) can be used to overcome these challenges by inducing lateral roots massively and synchronously. The combination of LRISs with transcriptomic approaches significantly advanced our insights in the molecular control of lateral root formation, in particular for Arabidopsis. Despite this success, LRISs have been underutilized for other plant species or for lateral root developmental stages later than the initiation. In this study, we developed and/or adapted LRISs in rice, Medicago, and Arabidopsis to perform RNA-sequencing during time courses that cover different developmental stages of lateral root formation and primordium development. As such, our study provides three extensive datasets of gene expression profiles during lateral root development in three different plant species. The three LRISs are highly effective but timing and spatial distribution of lateral root induction vary among the species. Detailed characterization of the stages in time and space in the respective species enabled an interspecies co-expression analysis to identify conserved players involved in lateral root development, as illustrated for the AUX/IAA and LBD gene families. Overall, our results provide a valuable resource to identify potentially conserved regulatory mechanisms in lateral root development, and as such will contribute to a better understanding of the complex regulatory network underlying lateral root development., (© 2023 Society for Experimental Biology and John Wiley & Sons Ltd.)

Caucasian clover ( Trifolium ambiguum M. Bieb.) is an excellent perennial plant in the legume family Fabaceae, with a well-developed rhizome and strong clonal growth. Auxin is one of the most important phytohormones in plants and plays an important role in plant growth and development. Auxin response factor (ARF) can regulate the expression of auxin-responsive genes, thus participating in multiple pathways of auxin transduction signaling in a synergistic manner. No genomic database has been established for Caucasian clover. In this study, 71 TaARF genes were identified through a transcriptomic database of Caucasian clover rhizome development. Phylogenetic analysis grouped the TaARFs into six (1-6) clades. Thirty TaARFs contained a complete ARF structure, including three relatively conserved regions. Physical and chemical property analysis revealed that TaARFs are unstable and hydrophilic proteins. We also analyzed the expression pattern of TaARFs in different tissues (taproot, horizontal rhizome, swelling of taproot, rhizome bud and rhizome bud tip). Quantitative real-time RT-PCR revealed that all TaARFs were responsive to phytohormones (indole-3-acetic acid, gibberellic acid, abscisic acid and methyl jasmonate) in roots, stems and leaves. These results helped elucidate the role of ARFs in responses to different hormone treatments in Caucasian clover.

While shaping of plant microbiome composition through 'host filtering' is well documented in legume-rhizobium symbioses, it is less clear to what extent different varieties and genotypes of the same plant species differentially influence symbiont community diversity and composition. Here, we compared how clover host varieties and genotypes affect the structure of Rhizobium populations in root nodules under conventional field and controlled greenhouse conditions. We first grew four Trifolium repens (white clover) F 2 crosses and one variety in a conventional field trial and compared differences in root nodule Rhizobium leguminosarum symbiovar trifolii (Rlt) genotype diversity using high-throughput amplicon sequencing of chromosomal housekeeping (rpoB and recA) genes and auxiliary plasmid-borne symbiosis genes (nodA and nodD). We found that Rlt nodule diversities significantly differed between clover crosses, potentially due to host filtering. However, variance in Rlt diversity largely overlapped between crosses and was also explained by the spatial distribution of plants in the field, indicative of the role of local environmental conditions for nodule diversity. To test the effect of host filtering, we conducted a controlled greenhouse trial with a diverse Rlt inoculum and several host genotypes. We found that different clover varieties and genotypes of the same variety selected for significantly different Rlt nodule communities and that the strength of host filtering (deviation from the initial Rhizobium inoculant composition) was positively correlated with the efficiency of symbiosis (rate of plant greenness colouration). Together, our results suggest that selection by host genotype and local growth conditions jointly influence white clover Rlt nodule diversity and community composition., (© 2023 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

White clover (Trifolium repens L.; Fabaceae) is an important forage and cover crop in agricultural pastures around the world and is increasingly used in evolutionary ecology and genetics to understand the genetic basis of adaptation. Historically, improvements in white clover breeding practices and assessments of genetic variation in nature have been hampered by a lack of high-quality genomic resources for this species, owing in part to its high heterozygosity and allotetraploid hybrid origin. Here, we use PacBio HiFi and chromosome conformation capture (Omni-C) technologies to generate a chromosome-level, haplotype-resolved genome assembly for white clover totaling 998 Mbp (scaffold N50 = 59.3 Mbp) and 1 Gbp (scaffold N50 = 58.6 Mbp) for haplotypes 1 and 2, respectively, with each haplotype arranged into 16 chromosomes (8 per subgenome). We additionally provide a functionally annotated haploid mapping assembly (968 Mbp, scaffold N50 = 59.9 Mbp), which drastically improves on the existing reference assembly in both contiguity and assembly accuracy. We annotated 78,174 protein-coding genes, resulting in protein BUSCO completeness scores of 99.6% and 99.3% against the embryophyta_odb10 and fabales_odb10 lineage datasets, respectively., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Aims: This study evaluated the red clover (Trifolium pratense) root-associated microbiota to clarify the presence of pathogenic and beneficial microorganisms in 89 Swedish field sites., Methods and Results: 16S rRNA and ITS amplicon sequencing analysis were performed on DNA extracted from the red clover root samples collected to determine the composition of the prokaryotic and eukaryotic root-associated microbe communities. Alpha and beta diversities were calculated and relative abundance of various microbial taxa and their co-occurrence were analyzed. Rhizobium was the most prevalent bacterial genus, followed by Sphingomonas, Mucilaginibacter, Flavobacterium, and the unclassified Chloroflexi group KD4-96. The Leptodontidium, Cladosporium, Clonostachys, and Tetracladium fungal genera known for endophytic, saprotrophic, and mycoparasitic lifestyles were also frequently observed in all samples. Sixty-two potential pathogenic fungi were identified with a bias toward grass pathogens and a higher abundance in samples from conventional farms., Conclusions: We showed that the microbial community was mainly shaped by geographic location and management procedures. Co-occurrence networks revealed that the Rhizobiumleguminosarum bv. trifolii was negatively associated with all fungal pathogenic taxa recognized in this study., (© The Author(s) 2023. Published by Oxford University Press on behalf of Applied Microbiology International.)

Medicago polymorpha L. (bur clover), an invasive plant species of the genus Medicago, has been traditionally used in China as an edible vegetable crop because of its high nutritive value. However, few molecular markers for M. polymorpha have been identified. Using the recently published high-quality reference genome of M. polymorpha , we performed a specific-locus amplified fragment sequencing (SLAF-seq) analysis of 10 M. polymorpha accessions to identify molecular markers and explore genetic diversity. A total of 52,237 high-quality single nucleotide polymorphisms (SNPs) were developed. These SNPs were mostly distributed on pseudochromosome 3, least distributed on pseudochromosome 7, and relatively evenly distributed on five other pseudochromosomes of M. polymorpha . Phenotypic analysis showed that there was a great difference in phenotypic traits among different M. polymorpha accessions. Moreover, clustering all M. polymorpha accessions based on their phenotypic traits revealed three groups. Both phylogenetic analysis and principal component analysis (PCA) of all M. polymorpha accessions based on SNP markers consistently indicated that all M. polymorpha accessions could be divided into three distinct groups (I, II, and III). Subsequent genetic diversity analysis for the 10 M. polymorpha accessions validated the effectiveness of the M. polymorpha germplasm molecular markers in China. Additionally, SSR mining analysis was also performed to identify polymorphic SSR motifs, which could provide valuable candidate markers for the further breeding of M. polymorpha . Since M. polymorpha genetics have not been actively studied, the molecular markers generated from our research will be useful for further research on M. polymorpha resource utilization and marker-assisted breeding., Competing Interests: The authors declare that they have no competing interests., (© 2023 Ren et al.)

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., (© 2022. The Author(s).)

Key Message: High variability for and candidate loci associated with resistance to southern anthracnose and clover rot in a worldwide collection of red clover provide a first basis for genomics-assisted breeding. Red clover (Trifolium pratense L.) is an important forage legume of temperate regions, particularly valued for its high yield potential and its high forage quality. Despite substantial breeding progress during the last decades, continuous improvement of cultivars is crucial to ensure yield stability in view of newly emerging diseases or changing climatic conditions. The high amount of genetic diversity present in red clover ecotypes, landraces, and cultivars provides an invaluable, but often unexploited resource for the improvement of key traits such as yield, quality, and resistance to biotic and abiotic stresses. A collection of 397 red clover accessions was genotyped using a pooled genotyping-by-sequencing approach with 200 plants per accession. Resistance to the two most pertinent diseases in red clover production, southern anthracnose caused by Colletotrichum trifolii, and clover rot caused by Sclerotinia trifoliorum, was assessed using spray inoculation. The mean survival rate for southern anthracnose was 22.9% and the mean resistance index for clover rot was 34.0%. Genome-wide association analysis revealed several loci significantly associated with resistance to southern anthracnose and clover rot. Most of these loci are in coding regions. One quantitative trait locus (QTL) on chromosome 1 explained 16.8% of the variation in resistance to southern anthracnose. For clover rot resistance we found eight QTL, explaining together 80.2% of the total phenotypic variation. The SNPs associated with these QTL provide a promising resource for marker-assisted selection in existing breeding programs, facilitating the development of novel cultivars with increased resistance against two devastating fungal diseases of red clover., (© 2022. The Author(s).)

Pod shattering is a seed dispersal strategy and an important agronomical trait in domesticated crops. The relationship between pod shattering and pod morphology in the genus Medicago is well known; however, the detailed mechanism underlying pod dehiscence in Medicago ruthenica, a perennial legume used for forage production, is unknown. Here, the pod ventral sutures of shatter-resistant and shatter-susceptible M. ruthenica genotypes were examined at 8, 12, 16, and 20 d after flowering. The mechanism of pod shattering was analyzed through microscopic observations, polygalacturonase (PG) and cellulase (CE) activity analyses, and RNA-sequencing (RNA-Seq), and the results were verified via reverse transcriptase-quantitative polymerase chain reaction. Pod shattering at the ventral suture in M. ruthenica occurs via a combination of two mechanisms: degradation of the middle lamella at the abscission layers (ALs) and detachment of lignified cells on either side of the ALs triggered by physical forces. Increased PG and CE activities in the pod ventral suture are essential for AL cell-autolysis in the shatter-susceptible genotype. RNA-Seq revealed that 11 genes encoding PG and CE were highly expressed in the ventral sutures of the shatter-susceptible genotype. The expression levels of auxin biosynthesis-related genes decreased in the AL cells and they were negatively associated with pod dehiscence. These results enhance our understanding of the pod shattering mechanism not only in M. ruthenica but also in other leguminous plants., (Copyright © 2022 Elsevier Masson SAS. All rights reserved.)

Medicago falcata is one of the leguminous forage crops, which grows well in arid and semiarid region. To fully investigate the mechanism of drought resistance response in M. falcata, we challenged the M. falcata plants with 30% PEG-6000, and performed physiological and transcriptome analyses. It was found that, the activities of antioxidant enzymes (eg. SOD, POD, and CAT) and soluble sugar content were all increased in the PEG-treated group, as compared to the control group. Transcriptome results showed that a total of 706 genes were differentially expressed in the PEG-treated plants in comparison with the control. Gene enrichment analyses on differentially expressed genes revealed that a number of genes in various pathway were significantly enriched, including the phenylpropanoid biosynthesis (ko00940) and glycolysis/gluconeogenesis (ko00010), indicating the involvement of these key pathways in drought response. Furthermore, the expression levels of seven differentially expressed genes were verified to be involved in drought response in M. falcata by qPCR. Taken together, these results will provide valuable information related to drought response in M. falcata and lay a foundation for molecular studies and genetic breeding of legume crops in future research., Competing Interests: The authors have declared that no competing interests exist.

Plant nuclear genomes harbor sequence elements derived from the organelles (mitochondrion and plastid) through intracellular gene transfer (IGT). Nuclear genomes also show a dramatic range of repeat content, suggesting that any sequence can be readily amplified. These two aspects of plant nuclear genomes are well recognized but have rarely been linked. Through investigation of 31 Medicago taxa we detected exceptionally high post-IGT amplification of mitochondrial (mt) DNA sequences containing rps10 in the nuclear genome of Medicago polymorpha and closely related species. The amplified sequences were characterized as tandem arrays of five distinct repeat motifs (2157, 1064, 987, 971, and 587 bp) that have diverged from the mt genome (mitogenome) in the M. polymorpha nuclear genome. The mt rps10-like arrays were identified in seven loci (six intergenic and one telomeric) of the nuclear chromosome assemblies and were the most abundant tandem repeat family, representing 1.6-3.0% of total genomic DNA, a value approximately three-fold greater than the entire mitogenome in M. polymorpha. Compared to a typical mt gene, the mt rps10-like sequence coverage level was 691.5-7198-fold higher in M. polymorpha and closely related species. In addition to the post-IGT amplification, our analysis identified the canonical telomeric repeat and the species-specific satellite arrays that are likely attributable to an ancestral chromosomal fusion in M. polymorpha. A possible relationship between chromosomal instability and the mt rps10-like tandem repeat family in the M. polymorpha clade is discussed., (© 2022 Society for Experimental Biology and John Wiley & Sons Ltd.)

Key Message: Transformation of MruGSTU39 in M. ruthenica and alfalfa enhanced growth and survival of transgenic plants by up-regulating GST and glutathione peroxidase activity to detoxify ROS under drought stress. Glutathione S-transferases (GSTs) are ubiquitous supergene family which play crucial roles in detoxification of reactive oxygen species (ROS). Despite studies on GSTs, few studies have focused on them in perennial, wild plant species with high tolerance to environmental stress. Here, we identified 66 MruGST genes from the genome of Medicago ruthenica, a perennial legume species native to temperate grasslands with high tolerance to environmental stress. These genes were divided into eight classes based on their conserved domains, phylogenetic tree and gene structure, with the tau class being the most numerous. Duplication analysis revealed that GST family in M. ruthenica was expanded by segmental and tandem duplication. Several drought-responsive MruGSTs were identified by transcriptomic analyses. Of them, expression of MruGSTU39 was up-regulated much more in a tolerant accession by drought stress. Transformation of MruGSTU39 in M. ruthenica and alfalfa (Medicago sativa) enhanced growth and survival of transgenic seedlings than their wild-type counterparts under drought. We demonstrated that MruGSTU39 can detoxify ROS to reduce its damage to membrane by up-regulating activities of GST and glutathione peroxidase. Our findings provide full-scale knowledge on GST family in the wild legume M. ruthenica with high tolerance to drought, and highlight improvement tolerance of legume forages to drought using genomic information of M. ruthenica., (© 2021. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

The species belonging to the genus Medicago are considered a very important genetic resource at global level both for planet's food security and for sustainable rangelands management. The checklist of the Italian flora (2021) includes a total number of 40 Medicago species for Italy, and 27 for Campania region, with a number of doubtful records or related to species no more found in the wild. In this study, 10 Medicago species native to Campania region, and one archaeophyte ( M. sativa) , identified by means of morphological diagnostic characters, were analyzed in a blind test to assay the efficacy of nine microsatellite markers (five cp-SSRs and four n-SSRs). A total number of 33 individuals from 6 locations were sampled and genotyped. All markers were polymorphic, 40 alleles were obtained with n-SSRs ranging from 8-12 alleles per locus with an average of 10 alleles per marker, PIC values ranged from 0.672 to 0.847, and the most polymorphic SSR was MTIC 564. The cp-SSRs markers were highly polymorphic too; PIC values ranged from 0.644 to 0.891 with an average of 0.776, the most polymorphic cp-SSR was CCMP10. 56 alleles were obtained with cp-SSRs ranging from 7 to 17 alleles per locus with an average of 11. AMOVA analysis with n-SSR markers highlighted a great level of genetic differentiation among the 11 species, with a statistically significant fixation index ( F ST ). UPGMA clustering and Bayesian-based population structure analysis assigned these 11 species to two main clusters, but the distribution of species within clusters was not the same for the two analyses. In conclusion, our results demonstrated that the combination of the used SSRs well distinguished the 11 Medicago species. Moreover, our results demonstrated that the use of a limited number of SSRs might be considered for further genetic studies on other Medicago species.

Background: Although plastomes are highly conserved with respect to gene content and order in most photosynthetic angiosperms, extensive genomic rearrangements have been reported in Fabaceae, particularly within the inverted repeat lacking clade (IRLC) of Papilionoideae. Two hypotheses, i.e., the absence of the IR and the increased repeat content, have been proposed to affect the stability of plastomes. However, this is still unclear for the IRLC species. Here, we aimed to investigate the relationships between repeat content and the degree of genomic rearrangements in plastomes of Medicago and its relatives Trigonella and Melilotus, which are nested firmly within the IRLC., Results: We detected abundant repetitive elements and extensive genomic rearrangements in the 75 newly assembled plastomes of 20 species, including gene loss, intron loss and gain, pseudogenization, tRNA duplication, inversion, and a second independent IR gain (IR ~ 15 kb in Melilotus dentata) in addition to the previous first reported cases in Medicago minima. We also conducted comparative genomic analysis to evaluate plastome evolution. Our results indicated that the overall repeat content is positively correlated with the degree of genomic rearrangements. Some of the genomic rearrangements were found to be directly linked with repetitive sequences. Tandem repeated sequences have been detected in the three genes with accelerated substitution rates (i.e., accD, clpP, and ycf1) and their length variation could be explained by the insertions of tandem repeats. The repeat contents of the three localized hypermutation regions around these three genes with accelerated substitution rates are also significantly higher than that of the remaining plastome sequences., Conclusions: Our results suggest that IR reemergence in the IRLC species does not ensure their plastome stability. Instead, repeat-mediated illegitimate recombination is the major mechanism leading to genome instability, a pattern in agreement with recent findings in other angiosperm lineages. The plastome data generated herein provide valuable genomic resources for further investigating the plastome evolution in legumes., (© 2021. The Author(s).)

Human induced environmental change may require rapid adaptation of plant populations and crops, but the genomic basis of environmental adaptation remain poorly understood. We analysed polymorphic loci from the perennial crop Medicago sativa (alfalfa or lucerne) and the annual legume model species M. truncatula to search for a common set of candidate genes that might contribute to adaptation to abiotic stress in both annual and perennial Medicago species. We identified a set of candidate genes of adaptation associated with environmental gradients along the distribution of the two Medicago species. Candidate genes for each species were detected in homologous genomic linkage blocks using genome-environment (GEA) and genome-phenotype association analyses. Hundreds of GEA candidate genes were species-specific, of these, 13.4% (M. sativa) and 24% (M. truncatula) were also significantly associated with phenotypic traits. A set of 168 GEA candidates were shared by both species, which was 25.4% more than expected by chance. When combined, they explained a high proportion of variance for certain phenotypic traits associated with adaptation. Genes with highly conserved functions dominated among the shared candidates and were enriched in gene ontology terms that have shown to play a central role in drought avoidance and tolerance mechanisms by means of cellular shape modifications and other functions associated with cell homeostasis. Our results point to the existence of a molecular basis of adaptation to abiotic stress in Medicago determined by highly conserved genes and gene functions. We discuss these results in light of the recently proposed omnigenic model of complex traits., (© 2021 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Caucasian clover is the only perennial herb of the genus Leguminous clover with underground rhizomes. However, we know very little about its development process and mechanism. Transcriptome studies were conducted on the roots of Caucasian clover without a rhizome (NR) at the young seedling stage and the fully developed rhizome, including the root neck (R1), main root (R2), horizontal root (R3), and rhizome bud (R4), of the tissues in the mature phase. Compared with the rhizome in the mature phase, NR had 893 upregulated differentially expressed genes (DEGs), most of which were enriched in 'phenylpropanoid biosynthesis', 'phenylalanine metabolism', 'DNA replication' and 'biosynthesis of amino acids'. A higher number of transcription factors (AP2/ERF, C2H2 and FAR1) were found in NR. There were highly expressed genes for R4, such as auxin response factor SAUR, galacturonosyltransferase (GAUT), and sucrose synthase (SUS). Phenylpropanoids are very important for the entire process of rhizome development. We drew a cluster heat map of genes related to the phenylpropanoid biosynthesis pathway, in which the largest number of genes belonged to COMT, and most of them were upregulated in R4., Competing Interests: The authors have declared that no competing interests exist.

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., (© 2021 John Wiley & Sons Ltd.)

Mechanisms that allow plants to survive and reproduce after herbivory are considered to play a key role in plant evolution. In this study, we evaluated how tolerance varies in species with different historic exposure to herbivores considering ontogeny. We exposed the range-restricted species Medicago citrina and its closely related and widespread species M. arborea to one and two herbivory simulations (80 % aerial biomass loss). Physiological and growth parameters related to tolerance capacity were assessed to evaluate constitutive values (without herbivory) and induced tolerance after damage. Constitutive traits were not always related to greater tolerance, and each species compensated for herbivory through different traits. Herbivory damage only led to mortality in M. citrina; adults exhibited root biomass loss and increased oxidative stress after damage, but also compensated aerial biomass. Despite seedlings showed a lower death percentage than adults after herbivory in M. citrina, they showed less capacity to recover control values than adults. Moderate tolerance to M. arborea herbivory and low tolerance to M. citrina is found. Thus, although the constitutive characteristics are maintained in the lineage, the tolerance of plants decreases in M. citrina. That represents how plants respond to the lack of pressure from herbivores in their habitat., (Copyright © 2021 Elsevier B.V. All rights reserved.)

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.

Auxin induced in root culture (AIR12) is a single gene in Arabidopsis and codes for a mono-heme cytochrome b, but it is unknown whether plant AIR12 is involved in abiotic stress responses. MfAIR12 was identified from Medicago falcata that is legume germplasm with great cold tolerance. Transcript levels of MfAIR12 and its homolog MtAIR12 from Medicago truncatula was induced under low temperature. Overexpression of MfAIR12 led to the accumulation of H 2 O 2 in apoplast and enhanced cold tolerance, which was blocked by H 2 O 2 scavengers, indicating that the increased cold tolerance was dependent upon the accumulated H 2 O 2 . In addition, declined cold tolerance was observed in Arabidopsis mutant air12, which could be restored by expressing MfAIR12. Compared to the wild type, higher levels of ascorbic acid and ascorbate redox state, as well as transcripts of the C repeat/dehydration responsive element-binding factor (CBF) transcription factors and their downstream cold-responsive genes, were observed in MfAIR12 transgenic lines, but lower levels of those in air12 mutant. It is suggested AIR12 confers cold tolerance as a result of the altered H 2 O 2 in the apoplast that is signaling in the regulation of CBF cold response pathway and ascorbate homeostasis., (© 2021 John Wiley & Sons Ltd.)

Although most invasive species engage in mutualism, we know little about how mutualism evolves as partners colonize novel environments. Selection on cooperation and standing genetic variation for mutualism traits may differ between a mutualism's invaded and native ranges, which could alter cooperation and coevolutionary dynamics. To test for such differences, we compare mutualism traits between invaded- and native-range host-symbiont genotype combinations of the weedy legume, Medicago polymorpha, and its nitrogen-fixing rhizobium symbiont, Ensifer medicae, which have coinvaded North America. We find that mutualism benefits for plants are indistinguishable between invaded- and native-range symbioses. However, rhizobia gain greater fitness from invaded-range mutualisms than from native-range mutualisms, and this enhancement of symbiont fecundity could increase the mutualism's spread by increasing symbiont availability during plant colonization. Furthermore, mutualism traits in invaded-range symbioses show lower genetic variance and a simpler partitioning of genetic variance between host and symbiont sources, compared to native-range symbioses. This suggests that biological invasion has reduced mutualists' potential to respond to coevolutionary selection. Additionally, rhizobia bearing a locus (hrrP) that can enhance symbiotic fitness have more exploitative phenotypes in invaded-range than in native-range symbioses. These findings highlight the impacts of biological invasion on the evolution of mutualistic interactions., (© 2021 The Authors. Evolution © 2021 The Society for the Study of Evolution.)

Co-expressed genes tend to have regulatory relationships and participate in similar biological processes. Construction of gene correlation networks from microarray or RNA-seq expression data has been widely applied to study transcriptional regulatory mechanisms and metabolic pathways under specific conditions. Furthermore, since transcription factors (TFs) are critical regulators of gene expression, it is worth investigating TFs on the promoters of co-expressed genes. Although co-expressed genes and their related metabolic pathways can be easily identified from previous resources, such as EXPath and EXPath Tool, this information is not simultaneously available to identify their regulatory TFs. EXPath 2.0 is an updated database for the investigation of regulatory mechanisms in various plant metabolic pathways with 1,881 microarray and 978 RNA-seq samples. There are six significant improvements in EXPath 2.0: (i) the number of species has been extended from three to six to include Arabidopsis, rice, maize, Medicago, soybean and tomato; (ii) gene expression at various developmental stages have been added; (iii) construction of correlation networks according to a group of genes is available; (iv) hierarchical figures of the enriched Gene Ontology (GO) terms are accessible; (v) promoter analysis of genes in a metabolic pathway or correlation network is provided; and (vi) user's gene expression data can be uploaded and analyzed. Thus, EXPath 2.0 is an updated platform for investigating gene expression profiles and metabolic pathways under specific conditions. It facilitates users to access the regulatory mechanisms of plant biological processes. The new version is available at http://EXPath.itps.ncku.edu.tw., (© The Author(s) 2020. Published by Oxford University Press on behalf of Japanese Society of Plant Physiologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Fenugreek as a self-pollinated plant is ideal for genome-wide association mapping where traits can be marked by their association with natural mutations. However, fenugreek is poorly investigated at the genomic level due to the lack of information regarding its genome. To fill this gap, we genotyped a collection of 112 genotypes with 153,881 SNPs using double digest restriction site-associated DNA sequencing. We used 38,142 polymorphic SNPs to prove the suitability of the population for association mapping. One significant SNP was associated with both seed length and seed width, and another SNP was associated with seed color. Due to the lack of a comprehensive genetic map, it is neither possible to align the newly developed markers to chromosomes nor to predict the underlying genes. Therefore, systematic targeting of those markers to homologous genomes of other legumes can overcome those problems. A BLAST search using the genomic fenugreek sequence flanking the identified SNPs showed high homology with several members of the Trifolieae tribe indicating the potential of translational approaches to improving our understanding of the fenugreek genome. Using such a comprehensively-genotyped fenugreek population is the first step towards identifying genes underlying complex traits and to underpin fenugreek marker-assisted breeding programs.

Key Message: MfLEA3 is involved in protection of catalase activity and confers multiple abiotic stress tolerance. Late embryogenesis abundant (LEA) proteins are involved in plant growth, development and abiotic stress tolerance. A member of group 3 LEA proteins from Medicago sativa subsp. falcata (L.) Arcang, MfLEA3, was investigated in the study. MfLEA3 transcript was induced in response to cold, dehydration, and abscisic acid (ABA), while the cold-induced transcript of MfLEA3 was blocked by pretreatment with inhibitor of ABA synthesis. Constitutive expression of MfLEA3 led to enhanced tolerance to cold, drought, and high-light stress in transgenic tobacco plants. Compared to accumulated reactive oxygen species (ROS) in the wild-type in response to treatments with low temperature, drought, and high light, ROS were not accumulated in transgenic plants. Superoxide dismutase, catalase (CAT), and ascorbate-peroxidase activities were increased in all plants after treatments with the above stresses, while higher CAT activity was maintained in transgenic plants compared with wild-type. However, transcript level of CAT-encoding genes including CAT1, CAT2, and CAT3 showed no significant difference between transgenic plants and wild-type, indicating that the higher CAT activity was not associated with its gene expression. ABA sensitivity and transcripts of several ABA and stress-responsive genes showed no difference between transgenic plant and wild-type, indicating that ABA signaling was not affected by constitutive expression of MfLEA3. The results suggest that MfLEA3 may be involved in the protection of CAT activity and confers multiple abiotic stress tolerance.

Root system architecture (RSA) influences the effectiveness of resources acquisition from soils but the genetic networks that control RSA remain largely unclear. We used rhizoboxes, X-ray computed tomography, grafting, auxin transport measurements and hormone quantification to demonstrate that Arabidopsis and Medicago CEP (C-TERMINALLY ENCODED PEPTIDE)-CEP RECEPTOR signalling controls RSA, the gravitropic set-point angle (GSA) of lateral roots (LRs), auxin levels and auxin transport. We showed that soil-grown Arabidopsis and Medicago CEP receptor mutants have a narrower RSA, which results from a steeper LR GSA. Grafting showed that CEPR1 in the shoot controls GSA. CEP receptor mutants exhibited an increase in rootward auxin transport and elevated shoot auxin levels. Consistently, the application of auxin to wild-type shoots induced a steeper GSA and auxin transport inhibitors counteracted the CEP receptor mutant's steep GSA phenotype. Concordantly, CEP peptides increased GSA and inhibited rootward auxin transport in wild-type but not in CEP receptor mutants. The results indicated that CEP-CEP receptor-dependent signalling outputs in Arabidopsis and Medicago control overall RSA, LR GSA, shoot auxin levels and rootward auxin transport. We propose that manipulating CEP signalling strength or CEP receptor downstream targets may provide means to alter RSA., (© 2020 The Authors. New Phytologist © 2020 New Phytologist Trust.)

Lotus japonicus is a herbaceous perennial legume that has been used extensively as a genetically tractable model system for deciphering the molecular genetics of symbiotic nitrogen fixation. Our aim is to improve the L. japonicus reference genome sequence, which has so far been based on Sanger and Illumina sequencing reads from the L. japonicus accession MG-20 and contained a large fraction of unanchored contigs. Here, we use long PacBio reads from L. japonicus Gifu combined with Hi-C data and new high-density genetic maps to generate a high-quality chromosome-scale reference genome assembly for L. japonicus. The assembly comprises 554 megabases of which 549 were assigned to six pseudomolecules that appear complete with telomeric repeats at their extremes and large centromeric regions with low gene density. The new L. japonicus Gifu reference genome and associated expression data represent valuable resources for legume functional and comparative genomics. Here, we provide a first example by showing that the symbiotic islands recently described in Medicago truncatula do not appear to be conserved in L. japonicus., (© The Author(s) 2020. Published by Oxford University Press on behalf of Kazusa DNA Research Institute.)

The No apical meristem-Arabidopsis transcription activation factor-Cup-shaped cotyledon (NAC) proteins play vital roles in plant development processes and responses to abiotic stress. In this study, 146 unigenes were identified as NAC genes from wild Medicago falcata L. by RNA sequencing. Among these were 30 full-length NACs, which, except for MfNAC63, MfNAC64 and MfNAC91, contained a complete DNA-binding domain and a variable transcriptional activation region. Sequence analyses of MfNACs along with their Arabidopsis thaliana (L.) Heynh. counterparts allowed these proteins to be phylogenetically classified into nine groups. MfNAC35, MfNAC88, MfNAC79, MfNAC26 and MfNAC95 were found to be stress-responsive genes. The eight MfNAC genes that were chosen for further analysis had different expression abilities in the leaves, stems and roots of M. falcata. Additionally, their expression levels were regulated by salinity, drought and cold stress, and ABA. This study will be useful for understanding the roles of MfNACs in wild M. falcata and could provide important information for the selection of candidate genes associated with stress tolerance.

A field study was conducted on smallholder farmer fields between 2012 to 2014 to evaluate the performance of cv. Agaitti Berseem-2002, against local landraces exchanged between farmers (LBF1) or available from local markets (LBM1). The effects of genotype and harvesting regimen on forage production, quality and seed production were evaluated. Significant differences (P < 0.05) among genotypes and cutting treatments were recorded for forage and seed yields, and forage quality across all research sites in both years. Maximum cumulative fresh forage (89.7 t/ha) and dry matter (DM; 13.4 t/ha) yields were obtained with Agaitti Berseem-2002 when harvesting occurred five times over the season. However, maximum seed yield (1048 kg/ha) with higher 1000-seed weight (3.63 g) were obtained if forage was only harvested three times and the crop then left for seed set. Agaitti Berseem-2002 also produced forage with the higher crude protein content (27%), DM digestibility (69%), digestible organic matter (DM basis; 65%) and metabolizable energy content (10%) compared to the local landraces (LBF1 and LBM1). Therefore, the harvesting regimen for greatest economic return which produced optimum fresh and DM forage yields of highest nutritive values and maximum seed yield, were comprised of taking three forage cuts (at 65, 110 and 150 days after sowing) prior to seed harvest.

Herbicide resistance is a recurrent evolutionary event that has been reported across many species and for all major herbicide modes of action. The synthetic auxinic herbicide 2,4-dichlorophenoxyacetic acid (2,4-D) has been widely used since the 1940s, however the genetic variation underlying naturally evolving resistance remains largely unknown. In this study, we used populations of the forage legume crop red clover (Trifolium pratense L.) that were recurrently selected for 2,4-D resistance to detect genome-wide signatures of adaptation. Four susceptible and six derived resistant populations were sequenced using a less costly approach by combining targeted sequencing (Capture-Seq) with pooled individuals (Pool-Seq). Genomic signatures of selection were identified using: (i) pairwise allele frequency differences; (ii) genome scan for overly differentiated loci; and (iii) genome-wide association. Fifty significant SNPs were consistently detected, most located in a single chromosome, which can be useful for marker assisted selection. Additionally, we searched for candidate genes at these genomic regions to gain insights into potential molecular mechanisms underlying 2,4-D resistance. Among the predicted functions of candidate genes, we found some related to the auxin metabolism, response to oxidative stress, and detoxification, which are also promising for further functional validation studies.

Background: Low temperature is one of the main environmental factors that limits crop growth, development, and production. Medicago falcata is an important leguminous herb that is widely distributed worldwide. M. falcata is related to alfalfa but is more tolerant to low temperature than alfalfa. Understanding the low temperature tolerance mechanism of M. falcata is important for the genetic improvement of alfalfa., Results: In this study, we explored the transcriptomic changes in the roots of low-temperature-treated M. falcata plants by combining SMRT sequencing and NGS technologies. A total of 115,153 nonredundant sequences were obtained, and 8849 AS events, 73,149 SSRs, and 4189 lncRNAs were predicted. A total of 111,587 genes from SMRT sequencing were annotated, and 11,369 DEGs involved in plant hormone signal transduction, protein processing in endoplasmic reticulum, carbon metabolism, glycolysis/gluconeogenesis, starch and sucrose metabolism, and endocytosis pathways were identified. We characterized 1538 TF genes into 45 TF gene families, and the most abundant TF family was the WRKY family, followed by the ERF, MYB, bHLH and NAC families. A total of 134 genes, including 101 whose expression was upregulated and 33 whose expression was downregulated, were differentially coexpressed at all five temperature points. PB40804, PB75011, PB110405 and PB108808 were found to play crucial roles in the tolerance of M. falcata to low temperature. WGCNA revealed that the MEbrown module was significantly correlated with low-temperature stress in M. falcata. Electrolyte leakage was correlated with most genetic modules and verified that electrolyte leakage can be used as a direct stress marker in physiological assays to indicate cell membrane damage from low-temperature stress. The consistency between the qRT-PCR results and RNA-seq analyses confirmed the validity of the RNA-seq data and the analysis of the regulatory mechanism of low-temperature stress on the basis of the transcriptome., Conclusions: The full-length transcripts generated in this study provide a full characterization of the transcriptome of M. falcata and may be useful for mining new low-temperature stress-related genes specific to M. falcata. These new findings could facilitate the understanding of the low-temperature-tolerance mechanism of M. falcata.

The explosive growth of genomic data provides an opportunity to make increased use of sequence variations for phenotype prediction. We have developed a prediction machine for quantitative phenotypes (WhoGEM) that overcomes some of the bottlenecks limiting the current methods. We demonstrated its performance by predicting quantitative disease resistance and quantitative functional traits in the wild model plant species, Medicago truncatula, using geographical locations as covariates for admixture analysis. The method's prediction reliability equals or outperforms all existing algorithms for quantitative phenotype prediction. WhoGEM analysis produces evidence that variation in genome admixture proportions explains most of the phenotypic variation for quantitative phenotypes.

Background: An eukaryotic translation elongation factor-2 (eEF-2) plays an important role in protein synthesis, however, investigation on its role in abiotic stress responses is limited. A cold responsive eEF2 named as MfEF2 was isolated from yellow-flowered alfalfa [Medicago sativa subsp. falcata (L.) Arcang, thereafter M. falcata], a forage legume with great cold tolerance, and transgenic tobacco (Nicotiana tabacum L.) plants overexpressing MfEF2 were analyzed in cold tolerance and proteomic profiling was conducted under low temperature in this study., Results: MfEF2 transcript was induced and peaked at 24 h and remained at the high level during cold treatment up to 96 h. Overexpression of MfEF2 in trasngenic tobacco plants resulted in enhanced cold tolerance. Compared to the wild type, transgenic plants showed higher survival rate after freezing treatment, higher levels of net photosynthetic rate (A), maximum photochemical efciency of photosystem (PS) II (F v /F m ) and nonphotochemical quenching (NPQ) and lower levels of ion leakage and reactive oxygen species (ROS) production after chilling treatment. iTRAQ-based quantitative proteomic analysis identified 336 differentially expressed proteins (DEPs) from leaves of one transgenic line versus the wild type after chilling treatment for 48 h. GO and KEGG enrichment were conducted for analysis of the major biological process, cellular component, molecular function, and pathways of the DEPs involving in. It is interesting that many down-regulated DEPs were grouped into "photosynthesis" and "photosynthesis-antenna", such as subunits of PSI and PSII as well as light harvesting chlorophyll protein complex (LHC), while many up-regulated DEPs were grouped into "spliceosome"., Conclusions: The results suggest that MfEF2 confers cold tolerance through regulating hundreds of proteins synthesis under low temperature conditions. The elevated cold tolerance in MfEF2 transgenic plants was associated with downregulation of the subunits of PSI and PSII as well as LHC, which leads to reduced capacity for capturing sunlight and ROS production for protection of plants, and upregulation of proteins involving in splicesome, which promotes alternative splicing of pre-mRNA under low temperature.

Recent studies on E3 of endoplasmic reticulum (ER)-associated degradation (ERAD) in plants have revealed homologs in yeast and animals. However, it remains unknown whether the plant ERAD system contains a plant-specific E3 ligase. Here, we report that MfSTMIR, which encodes an ER-membrane-localized RING E3 ligase that is highly conserved in leguminous plants, plays essential roles in the response of ER and salt stress in Medicago. MfSTMIR expression was induced by salt and tunicamycin (Tm). mtstmir loss-of-function mutants displayed impaired induction of the ER stress-responsive genes BiP1/2 and BiP3 under Tm treatment and sensitivity to salt stress. MfSTMIR promoted the degradation of a known ERAD substrate, CPY*. MfSTMIR interacted with the ERAD-associated ubiquitin-conjugating enzyme MtUBC32 and Sec61-translocon subunit MtSec61γ. MfSTMIR did not affect MtSec61γ protein stability. Our results suggest that the plant-specific E3 ligase MfSTMIR participates in the ERAD pathway by interacting with MtUBC32 and MtSec61γ to relieve ER stress during salt stress., (© 2019 The Authors. The Plant Journal published by John Wiley & Sons Ltd and Society for Experimental Biology.)

The plant genome comprises a coevolving, integrated genetic system housed in three subcellular compartments: the nucleus, mitochondrion, and the plastid. The typical land plant plastid genome (plastome) comprises the sum of repeating units of 130-160 kb in length. The plastome inverted repeat (IR) divides each plastome monomer into large and small single copy regions, an architecture highly conserved across land plants. There have been varying degrees of expansion or contraction of the IR, and in a few distinct lineages, including the IR-lacking clade of papilionoid legumes, one copy of the IR has been lost. Completion of plastome sequencing and assembly for 19 Medicago species and Trigonella foenum-graceum and comparative analysis with other IR-lacking clade taxa revealed modest divergence with regard to structural organization overall. However, one clade contained unique variation suggesting an ancestor had experienced repeat-mediated changes in plastome structure. In Medicago minima, a novel IR of ∼9 kb was confirmed and the role of repeat-mediated, recombination-dependent replication in IR reemergence is discussed., (© The Author(s) 2019. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The amino acid/auxin permease (AAAP) gene family plays an important role in the long-distance amino acid transport pathway and takes part in various stages of plant growth and development. However, little is known about the AAAP gene family in Medicago truncatula. Here, we identified 86 putative MtAAAP family members using genome sequence information. Based on phylogenetic analysis, these MtAAAP genes were categorized into eight distinct subfamilies. The MtAAAP genes were mapped on 8 chromosomes and duplication events appeared widely, with 19 and 21 pairs of MtAAAP genes showing segment and tandem duplication events, respectively. Ratio of Ka/Ks indicated that duplicated genes underwent purifying selection. Analysis of RNA-seq data showed that MtAAAP genes exhibited specific expression patterns among different tissues and abiotic stress, indicating that MtAAAP members were involved in plant developmental regulation and stress responses. Expression patterns of 16 MtAAAP genes under abiotic stress were verified by qRT-PCR. The present study provides a foundation for the functional analysis of MtAAAPs in developmental regulation and stress responses.

Legume-rhizobial symbiosis plays an important role in agriculture and ecological restoration. However, knowledge of the molecular mechanisms, especially the microstructure and global transcriptional profiling, of the symbiosis process under heavy metal contamination is limited. In this study, a heavy metal-tolerant legume, Medicago lupulina, was treated with different concentrations of copper (Cu). The results showed that the early infection process was inhibited and the nodule ultrastructure was changed under 200 mg kg -1 Cu stress. Most infection threads (ITs) were prevented from entering the nodule cells, and few rhizobia were released into the host cells, in which thickening of the plant cell wall and IT wall was observed, demonstrating that rhizobial invasion was inhibited under Cu stress. RNA-seq analysis indicated that a strong shift in gene expression occurred (3257 differentially expressed genes, DEGs). The most pronounced effect was the upregulation of a set of 71 of 73 DEGs for nodule-specific cysteine-rich peptides, which have been shown to control the terminal differentiation of rhizobia in the nodules and to have antimicrobial activity. Various genes for metal transport, chelation binding and antioxidant defence were regulated. In particular, the DEGs for Cu trafficking and detoxification were induced during nodule formation. The DEGs for ethylene (ET) biosynthesis and signalling were also differentially expressed during nodulation, suggesting that the inhibition of nodulation by Cu occurred partially through ET signalling. Furthermore, the genes related to the cell wall were mostly upregulated and most likely involved in cell wall thickening. These findings provide an integrated understanding of the effects of Cu on legume nodule symbiosis at the molecular and phenotypic levels., (Copyright © 2018. Published by Elsevier B.V.)

Genetic diversity between two ecotypes of Egyptian clover varieties, namely Fahl (mono-cut) and Helaly (multi-cut) have been assessed based on forage yield and yield components as well as molecular marker systems. The two parental genotypes were crossed to produce seeds of F 1 and F 2 progenies. Analyses of variance indicated significant differences between four populations (P 1 (Fahl), P 2 (Helaly), F 1 and F 2 ) for fresh forage yield, number of florets/inflorescence, number of seeds/inflorescence and 1000 seed weight. The mean of F 1 hybrid indicated over-dominance of the higher performance. The phenotypic and genotypic coefficients of variation were high for fresh forage yield, intermediate for 1000-seed weight and low for number of florets/inflorescence and number of seeds/inflorescence. Four molecular marker systems with 80 primers, 30 RAPD, 10 ISSR, 10 SRAP and 30 SSR were used for studying the genetic diversity between the two parents, out of which 64 primers (26 RAPD, 7 ISSR, 7 SRAP and 24 SSR) were polymorphic between the parents. The four molecular marker systems generated unique DNA bands for each parent. Twenty-one primers which produced higher unique bands in both parents were surveyed on bulked DNA from the extremes of four agro-morphological traits within and between the two ecotypes in F 2 generations. Twenty-one primers produced bands distinguish between the bulked extremes for at least one trait within each ecotype or between the two ecotypes. All polymorphic primers were subjected to QTL analysis, out of them 23 only were mapped on three linkage groups with four agro-morphological traits and showed 24 putative QTLs.

To establish and spread in a new location, an invasive species must be able to carry out its life cycle in novel environmental conditions. A key trait underlying fitness is the shift from vegetative to reproductive growth through floral development. In this study, we used a common garden experiment and genotyping-by-sequencing to test whether the latitudinal flowering cline of the North American invasive plant Medicago polymorpha was translocated from its European native range through multiple introductions, or whether the cline rapidly established due to evolution following a genetic bottleneck. Analysis of flowering time in 736 common garden plants showed a latitudinal flowering time cline in both the native and invaded ranges where genotypes from lower latitudes flowered earlier. Genotyping-by-sequencing of 9,658 SNPs in 446 individuals revealed two major subpopulations of M. polymorpha in the native range, only one of which is present in the invaded range. Additionally, native range populations have higher genetic diversity than invaded range populations, suggesting that a genetic bottleneck occurred during invasion. All invaded range individuals are closely related to plants collected from native range populations in Portugal and southern Spain, and population assignment tests assigned invaded range individuals to this same narrow source region. Taken together, our results suggest that latitudinal clinal variation in flowering time has rapidly evolved across the invaded range despite a genetic bottleneck following introduction., (© 2018 John Wiley & Sons Ltd.)

Pea (Pisum sativum) is one of relatively few genetically amenable plant species with compound leaves. Pea leaves have a variety of specialized organs: leaflets, tendrils, pulvini and stipules, which enable the identification of mutations that transform or affect distinct parts of the leaf. Characterization of these mutations offers insights into the development and evolution of novel leaf traits. The previously characterized morphological gene Cochleata, conferring stipule identity, was known to interact with Stipules reduced (St), which conditions stipule size in pea, but the St gene remained unknown. Here we analysed Fast Neutron irradiated pea mutants by restriction site associated DNA sequencing. We identified St as a gene encoding a C2H2 zinc finger transcription factor that is regulated by Cochleata. St regulates both cell division and cell expansion in the stipule. Our approach shows how systematic genome-wide screens can be used successfully for the analysis of traits in species for which whole genome sequences are not available., (© 2018 The Authors New Phytologist © 2018 New Phytologist Trust.)

Medicago flowering, like that of Arabidopsis, is promoted by vernalization and long days, but alternative mechanisms are predicted because Medicago lacks the key regulators CO and FLC. Three Medicago SOC1-like genes, including MtSOC1a, were previously implicated in flowering control, but no legume soc1 mutants with altered flowering were reported. Here, reverse transciption-quantitative PCR (RT-qPCR) indicated that the timing and magnitude of MtSOC1a expression was regulated by the flowering promoter FTa1, while in situ hybridization indicated that MtSOC1a expression increased in the shoot apical meristem during the floral transition. A Mtsoc1a mutant showed delayed flowering and short primary stems. Overexpression of MtSOC1a partially rescued the flowering of Mtsoc1a, but caused a dramatic increase in primary stem height, well before the transition to flowering. Internode cell length correlated with stem height, indicating that MtSOC1a promotes cell elongation in the primary stem. However, application of gibberellin (GA3) caused stem elongation in both the wild type and Mtsoc1a, indicating that the mutant was not defective in gibberellin responsiveness. These results indicate that MtSOC1a may function as a floral integrator gene and promotes primary stem elongation. Overall, this study suggests that apart from some conservation with the Arabidopsis flowering network, MtSOC1a has a novel role in regulating aspects of shoot architecture.