Genome-wide characterization, transcriptome profiling, and functional analysis of the ALMT gene family in Medicago for aluminum resistance

Aluminum (Al) is the major limiting factor affecting plant productivity in acidic soils. Al3+ ions exhibit increased solubility at a pH below 5, leading to plant root tip toxicity. Alternatively, plants can perceive very low concentrations of Al3+, 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.

Evolutionary genomics of climatic adaptation and resilience to climate change in alfalfa

Given the escalating impact of climate change on agriculture and food security, gaining insights into the evolutionary dynamics of climatic adaptation and uncovering climate-adapted variation empower the breeding of climate-resilience crops to face future climate change. Alfalfa (Medicago sativa subsp. sativa), the queen of forages with remarkable adaptability across diverse global environments, is an excellent model for investigating species' responses to climate change. We conducted population genomic analyses to unravel alfalfa's climatic adaptation and genetic susceptibility to future climate change, utilizing genome resequencing data from 702 accessions of 24 Medicago species. We found that interspecific genetic exchange has fueled the gene pool of alfalfa, particularly enriching defense and stress response genes. Inter-subspecific introgression between Medicago sativa subsp. falcata (subsp. falcata) and alfalfa not only aids alfalfa's climatic adaptation but also introduces genetic burden. A total of 1671 genes were associated with climatic adaptation, and 5.7% of them were introgression from subsp. falcata. Integrating climate-associated variants and climate data, we identified vulnerable populations to future climate change, particularly in higher latitudes of the northern hemisphere, serving as a clarion call for targeted conservation initiatives and breeding efforts. Moreover, we unveil pre-adaptive populations demonstrating heightened resilience to climate fluctuations, illuminating a pathway for future breeding strategies. This study enhances our understanding of alfalfa's local adaptation and facilitates breeding of climate-resilient cultivars, contributing to effective agricultural strategies facing future climate change.

The single-cell transcriptome program of nodule development cellular lineages in Medicago truncatula

Legumes establish a symbiotic relationship with nitrogen-fixing rhizobia by developing nodules. Nodules are modified lateral roots that undergo changes in their cellular development in response to bacteria, but the transcriptional reprogramming that occurs in these root cells remains largely uncharacterized. Here, we describe the cell-type-specific transcriptome response of Medicago truncatula roots to rhizobia during early nodule development in the wild-type genotype Jemalong A17, complemented with a hypernodulating mutant (sunn-4) to expand the cell population responding to infection and subsequent biological inferences. The analysis identifies epidermal root hair and stele sub-cell types associated with a symbiotic response to infection and regulation of nodule proliferation. Trajectory inference shows cortex-derived cell lineages differentiating to form the nodule primordia and, posteriorly, its meristem, while modulating the regulation of phytohormone-related genes. Gene regulatory analysis of the cell transcriptomes identifies new regulators of nodulation, including STYLISH 4, for which the function is validated.

Gene-edited Mtsoc1 triple mutant Medicago plants do not flower

Optimized flowering time is an important trait that ensures successful plant adaptation and crop productivity. SOC1-like genes encode MADS transcription factors, which are known to play important roles in flowering control in many plants. This includes the best-characterized eudicot model Arabidopsis thaliana (Arabidopsis), where SOC1 promotes flowering and functions as a floral integrator gene integrating signals from different flowering-time regulatory pathways. Medicago truncatula (Medicago) is a temperate reference legume with strong genomic and genetic resources used to study flowering pathways in legumes. Interestingly, despite responding to similar floral-inductive cues of extended cold (vernalization) followed by warm long days (VLD), such as in winter annual Arabidopsis, Medicago lacks FLC and CO which are key regulators of flowering in Arabidopsis. Unlike Arabidopsis with one SOC1 gene, multiple gene duplication events have given rise to three MtSOC1 paralogs within the Medicago genus in legumes: one Fabaceae group A SOC1 gene, MtSOC1a, and two tandemly repeated Fabaceae group B SOC1 genes, MtSOC1b and MtSOC1c. Previously, we showed that MtSOC1a has unique functions in floral promotion in Medicago. The Mtsoc1a Tnt1 retroelement insertion single mutant showed moderately delayed flowering in long- and short-day photoperiods, with and without prior vernalization, compared to the wild-type. In contrast, Mtsoc1b Tnt1 single mutants did not have altered flowering time or flower development, indicating that it was redundant in an otherwise wild-type background. Here, we describe the generation of Mtsoc1a Mtsoc1b Mtsoc1c triple mutant lines using CRISPR-Cas9 gene editing. We studied two independent triple mutant lines that segregated plants that did not flower and were bushy under floral inductive VLD. Genotyping indicated that these non-flowering plants were homozygous for the predicted strong mutant alleles of the three MtSOC1 genes. Gene expression analyses using RNA-seq and RT-qPCR indicated that these plants remained vegetative. Overall, the non-flowering triple mutants were dramatically different from the single Mtsoc1a mutant and the Arabidopsis soc1 mutant; implicating multiple MtSOC1 genes in critical overlapping roles in the transition to flowering in Medicago.

A genome-wide study of the lipoxygenase gene families in Medicago truncatula and Medicago sativa reveals that MtLOX24 participates in the methyl jasmonate response

BACKGROUND

Lipoxygenase (LOX) is a multifunctional enzyme that is primarily related to plant organ growth and development, biotic and abiotic stress responses, and production of flavor-associated metabolites. In higher plants, the LOX family encompasses several isozymes with varying expression patterns between tissues and developmental stages. These affect processes including seed germination, seed storage, seedling growth, fruit ripening, and leaf senescence. LOX family genes have multiple functions in response to hormones such as methyl jasmonate (MeJA) and salicylic acid.

RESULTS

In this study, we identified 30 and 95 LOX homologs in Medicago truncatula and Medicago sativa, respectively. These genes were characterized with analyses of their basic physical and chemical properties, structures, chromosomal distributions, and phylogenetic relationships to understand structural variations and their physical locations. Phylogenetic analysis was conducted for members of the three LOX subfamilies (9-LOX, type I 13-LOX, and type II 13-LOX) in Arabidopsis thaliana, Glycine max, M. truncatula, and M. sativa. Analysis of predicted promoter elements revealed several relevant cis-acting elements in MtLOX and MsLOX genes, including abscisic acid (ABA) response elements (ABREs), MeJA response elements (CGTCA-motifs), and antioxidant response elements (AREs). Cis-element data combined with transcriptomic data demonstrated that LOX gene family members in these species were most likely related to abiotic stress responses, hormone responses, and plant development. Gene expression patterns were confirmed via quantitative reverse transcription PCR. Several MtLOX genes (namely MtLOX15, MtLOX16, MtLOX20, and MtLOX24) belonging to the type I 13-LOX subfamily and other LOX genes (MtLOX7, MtLOX11, MsLOX23, MsLOX87, MsLOX90, and MsLOX94) showed significantly different expression levels in the flower tissue, suggesting roles in reproductive growth. Type I 13-LOXs (MtLOX16, MtLOX20, MtLOX21, MtLOX24, MsLOX57, MsLOX84, MsLOX85, and MsLOX94) and type II 13-LOXs (MtLOX5, MtLOX6, MtLOX9, MtLOX10, MsLOX18, MsLOX23, and MsLOX30) were MeJA-inducible and were predicted to function in the jasmonic acid signaling pathway. Furthermore, exogenous MtLOX24 expression in Arabidopsis verified that MtLOX24 was involved in MeJA responses, which may be related to insect-induced abiotic stress.

CONCLUSIONS

We identified six and four LOX genes specifically expressed in the flowers of M. truncatula and M. sativa, respectively. Eight and seven LOX genes were induced by MeJA in M. truncatula and M. sativa, and the LOX genes identified were mainly distributed in the type I and type II 13-LOX subfamilies. MtLOX24 was up-regulated at 8 h after MeJA induction, and exogenous expression in Arabidopsis demonstrated that MtLOX24 promoted resistance to MeJA-induced stress. This study provides valuable new information regarding the evolutionary history and functions of LOX genes in the genus Medicago.

The Extraordinary Alkali Bee, Nomia melanderi (Halictidae), the World's Only Intensively Managed Ground-Nesting Bee

Among the ground-nesting bees are several proven crop pollinators, but only the alkali bee (Nomia melanderi) has been successfully managed. In <80 years, it has become the world's most intensely studied ground-nesting solitary bee. In many ways, the bee seems paradoxical. It nests during the torrid, parched midsummer amid arid valleys and basins of the western United States, yet it wants damp nesting soil. In these basins, extensive monocultures of an irrigated Eurasian crop plant, alfalfa (lucerne), subsidize millions of alkali bees. Elsewhere, its polylectic habits and long foraging range allow it to stray into neighboring crops contaminated with insecticides. Primary wild floral hosts are either non-native or poorly known. Kleptoparasitic bees plague most ground nesters, but not alkali bees, which do, however, host other well-studied parasitoids. Building effective nesting beds requires understanding the hydraulic conductivity of silty nesting soils and its important interplay with specific soil mineral salts. Surprisingly, some isolated populations endure inhospitably cold climates by nesting amid hot springs. Despite the peculiarities and challenges associated with its management, the alkali bee remains the second most valuable managed solitary bee for US agriculture and perhaps the world.

Alleviation of gadolinium stress on Medicago by elevated atmospheric CO2 is mediated by changes in carbohydrates, Anthocyanin, and proline metabolism

Rare earth elements (REE) like Gadolinium (Gd), are increasingly used in industry and agriculture and this is concomitant with the increasingly leaking of Gd into the environment. Under a certain threshold concentration, REE can promote plant growth, however, beyond this concentration, they exert negative effects on plant growth. Moreover, the effect of Gd on plants growth and metabolism under a futuristic climate with increasingly atmospheric CO2 has not yet been studied. To this end, we investigated the effect of soil contamination with Gd (150 mg/kg soil) on the growth, carbohydrates, proline, and anthocyanin metabolism of Medicago plants grown under ambient (aCO2, 410 ppm) or elevated CO2 (eCO2, 720 ppm) concentration. Gd negatively affected the growth and photosynthesis of plants and imposed oxidative stress i.e., increased H2O2 and lipid peroxidation (MDA) level. As defense lines, the level and metabolism of osmoprotectants (soluble sugars and proline) and antioxidants (phenolics, anthocyanins, and tocopherols) were increased under Gd treatment. High CO2 positively affected the growth and metabolism of Medicago plants. Moreover, eCO2 mitigated the negative impacts of Gd on Medicago growth. It further induced the levels of osmoprotectants and antioxidants. In line with increased proline and anthocyanins, their metabolic enzymes (e.g. OAT, P5CS, PAL, and CS) were also increased. This study advances our understanding of how Gd adversely affects Medicago plant growth and metabolism. It also sheds light on the biochemical mechanisms underlying the Gd stress-reducing impact of eCO2.

Identification and evolution analysis of YUCCA genes of Medicago sativa and Medicago truncatula and their expression profiles under abiotic stress

The YUCCAs (YUC) are functionally identified flavin-containing monooxidases (FMOs) in plants that act as an important rate-limiting enzyme functioning in the auxin synthesis IPA (indole-3-pyruvic acid) pathway. In this study, 12 MsYUCs and 15 MtYUCs containing characteristic conserved motifs were identified in M. sativa (Medicago sativa L.) and M. truncatula (Medicago truncatula Gaertn.), respectively. Phylogenetic analysis revealed that YUC proteins underwent an evolutionary divergence. Both tandem and segmental duplication events were presented in MsYUC and MtYUC genes. Comparative syntenic maps of M. sativa with M. truncatula, Arabidopsis (Arabidopsis thaliana), or rice (Oryza sativa L.) were constructed to illustrate the evolution relationship of the YUC gene family. A large number of cis-acting elements related to stress response and hormone regulation were revealed in the promoter sequences of MsYUCs. Expression analysis showed that MsYUCs had a tissue-specific, genotype-differential expression and a differential abiotic stress response pattern based on transcriptome data analysis of M. sativa online. In addition, RT-qPCR confirmed that salt stress significantly induced the expression of MsYUC1/MsYUC10 but significantly inhibited MsYUC2/MsYUC3 expression and the expression of MsYUC10/MsYUC11/MsYUC12 was significantly induced by cold treatment. These results could provide valuable information for functional analysis of YUC genes via gene engineering of the auxin synthetic IPA pathway in Medicago.

A simplified protocol for Agrobacterium-mediated transformation of cell suspension cultures of the model species Medicago truncatula A17

This manuscript describes a unique protocol for the rapid transformation of Medicago truncatula A17 cell suspension cultures mediated by Agrobacterium tumefaciens. Medicago cells were collected on day 7 of the growth curve, which corresponded to the beginning of the exponential phase. They were then co-cultured with Agrobacterium for 3 days before being spread onto a petri dish with appropriate antibiotic selection. The Receptor Binding Domain of the Spike protein of SARS-CoV-2 was used as a model to develop this protocol. The presence of the transgene was assessed using PCR, and the integrity of the product was evaluated by SDS-PAGE and Western-blotting.

Diversification of Sinorhizobium populations associated with Medicago polymorpha and Medicago lupulina in purple soil of China

The double selection of environment adaptation and host specificity forced the diversification of rhizobia in nature. In the tropical region of China, Medicago polymorpha and Medicago lupulina are widely distributed, particularly in purple soil. However, the local distribution and diversity of rhizobia associated with these legumes has not been systematically investigated. To this end, root nodules of M. polymorpha and M. lupulina grown in purple soil at seven locations in Yunnan Province of China were collected for rhizobial isolation. The obtained rhizobia were characterized by RFLP of 16S-23S rRNA intergenic spacer, BOXAIR fingerprinting, and phylogeny of housekeeping and symbiosis genes. As result, a total of 91 rhizobial strains were classified into species Sinorhizobium medicae and S. meliloti, while three nodC gene types were identified among them. S. medicae containing nodC of type I was dominant in farmlands associated with M. polymorpha; while S. meliloti harboring nodC of type III was dominant in wild land nodulated by M. lupulina. For both rhizobial species, greater genetic diversity was detected in the populations isolated from their preferred host plant. A high level of genetic differentiation was observed between the two Sinorhizobium species, and gene flow was evident within the populations of the same species derived from different soil types, indicating that rhizobial evolution is likely associated with the soil features. To examine the effects of environmental features on rhizobial distribution, soil physicochemical traits and rhizobial genotypes were applied for constrained analysis of principle coordinates, which demonstrated that soil features like pH, nitrogen and sodium were the principle factors governing the rhizobial geographical distribution. Altogether, both S. medicae and S. meliloti strains could naturally nodulate with M. polymorpha and M. lupulina, but the rhizobium-legume symbiosis compatibility determined by both the host species and soil factors was also highlighted.

Highly Efficient Gene Knockout in Medicago truncatula Genotype R108 Using CRISPR-Cas9 System and an Optimized Agrobacterium Transformation Method

Medicago truncatula is the model plant species for studying symbioses with nitrogen-fixing rhizobia and arbuscular mycorrhizae, where edited mutants are invaluable for elucidating the contributions of known genes in these processes. Streptococcus pyogenes Cas9 (SpCas9)-based genome editing is a facile means of achieving loss of function, including where multiple gene knockouts are desired in a single generation. We describe how the user can customize our vector to target single or multiple genes, then how the vector is used to make M. truncatula transgenic plants containing target site mutations. Finally, obtaining transgene-free homozygous mutants is covered.

Paired Medicago receptors mediate broad-spectrum resistance to nodulation by Sinorhizobium meliloti carrying a species-specific gene

Plants have evolved the ability to distinguish between symbiotic and pathogenic microbial signals. However, potentially cooperative plant-microbe interactions often abort due to incompatible signaling. The Nodulation Specificity 1 (NS1) locus in the legume Medicago truncatula blocks tissue invasion and root nodule induction by many strains of the nitrogen-fixing symbiont Sinorhizobium meliloti. Controlling this strain-specific nodulation blockade are two genes at the NS1 locus, designated NS1a and NS1b, which encode malectin-like leucine-rich repeat receptor kinases. Expression of NS1a and NS1b is induced upon inoculation by both compatible and incompatible Sinorhizobium strains and is dependent on host perception of bacterial nodulation (Nod) factors. Both presence/absence and sequence polymorphisms of the paired receptors contribute to the evolution and functional diversification of the NS1 locus. A bacterial gene, designated rns1, is required for activation of NS1-mediated nodulation restriction. rns1 encodes a type I-secreted protein and is present in approximately 50% of the nearly 250 sequenced S. meliloti strains but not found in over 60 sequenced strains from the closely related species Sinorhizobium medicae. S. meliloti strains lacking functional rns1 are able to evade NS1-mediated nodulation blockade.

Herbivory modifies plant symbiont number and impact on host plant performance in the field

Species interactions are a unifying theme in ecology and evolution. Both fields are currently moving beyond their historical focus on isolated pairwise relationships to understand how ecological communities affect focal interactions. Additional species can modify both the number of interactions and the fitness consequences of each interaction (i.e., selection). Although only selection affects the evolution of the focal interaction, the two are often conflated, limiting our understanding of the evolution of multispecies interactions. We manipulated aboveground herbivory on the legume Medicago lupulina in the field and quantified its effect on number of symbionts and the per-symbiont impact on plant performance in two belowground symbioses: mutualistic rhizobia bacteria (Ensifer meliloti) and parasitic root-knot nematodes (Meloidogyne hapla). We found that herbivores modified the number of rhizobia nodules, as well as the benefit per nodule. However, each effect was specific to a distinct herbivory regime: natural herbivory affected nodule number, whereas leafhoppers (Cicadellidae) weakened the per nodule benefit. We did not detect any effect of herbivory on nematode gall number or the cost of infection. Our data demonstrate that distinguishing between symbiont number from the fitness consequences of symbiosis is crucial to accurately infer how pairwise interactions will evolve in a community.

Evolution of the CBL and CIPK gene families in Medicago: genome-wide characterization, pervasive duplication, and expression pattern under salt and drought stress

BACKGROUND

Calcineurin B-like proteins (CBLs) are ubiquitous Ca²⁺ sensors that mediate plant responses to various stress and developmental processes by interacting with CBL-interacting protein kinases (CIPKs). CBLs and CIPKs play essential roles in acclimatization of crop plants. However, evolution of these two gene families in the genus Medicago is poorly understood.

RESULTS

A total of 68 CBL and 135 CIPK genes have been identified in five genomes from Medicago. Among these genomes, the gene number of CBLs and CIPKs shows no significant difference at the haploid genome level. Phylogenetic and comprehensive characteristic analyses reveal that CBLs and CIPKs are classified into four clades respectively, which is validated by distribution of conserved motifs. The synteny analysis indicates that the whole genome duplication events (WGDs) have contributed to the expansion of both families. Expression analysis demonstrates that two MsCBLs and three MsCIPKs are specifically expressed in roots, mature leaves, developing flowers and nitrogen fixing nodules of Medicago sativa spp. sativa, the widely grown tetraploid species. In particular, the expression of these five genes was highly up-regulated in roots when exposed to salt and drought stress, indicating crucial roles in stress responses.

CONCLUSIONS

Our study leads to a comprehensive understanding of evolution of CBL and CIPK gene families in Medicago, but also provides a rich resource to further address the functions of CBL-CIPK complexes in cultivated species and their closely related wild relatives.

MtING2 encodes an ING domain PHD finger protein which affects Medicago growth, flowering, global patterns of H3K4me3, and gene expression

Flowering of the reference legume Medicago truncatula is promoted by winter cold (vernalization) followed by long-day photoperiods (VLD) similar to winter annual Arabidopsis. However, Medicago lacks FLC and CO, key regulators of Arabidopsis VLD flowering. Most plants have two INHIBITOR OF GROWTH (ING) genes (ING1 and ING2), encoding proteins with an ING domain with two anti-parallel alpha-helices and a plant homeodomain (PHD) finger, but their genetic role has not been previously described. In Medicago, Mting1 gene-edited mutants developed and flowered normally, but an Mting2-1 Tnt1 insertion mutant and gene-edited Mting2 mutants had developmental abnormalities including delayed flowering particularly in VLD, compact architecture, abnormal leaves with extra leaflets but no trichomes, and smaller seeds and barrels. Mting2 mutants had reduced expression of activators of flowering, including the FT-like gene MtFTa1, and increased expression of the candidate repressor MtTFL1c, consistent with the delayed flowering of the mutant. MtING2 overexpression complemented Mting2-1, but did not accelerate flowering in wild type. The MtING2 PHD finger bound H3K4me2/3 peptides weakly in vitro, but analysis of gene-edited mutants indicated that it was dispensable to MtING2 function in wild-type plants. RNA sequencing experiments indicated that >7000 genes are mis-expressed in the Mting2-1 mutant, consistent with its strong mutant phenotypes. Interestingly, ChIP-seq analysis identified >5000 novel H3K4me3 locations in the genome of Mting2-1 mutants compared to wild type R108. Overall, our mutant study has uncovered an important physiological role of a plant ING2 gene in development, flowering, and gene expression, which likely involves an epigenetic mechanism.

Combining GWAS and population genomic analyses to characterize coevolution in a legume-rhizobia symbiosis

The mutualism between legumes and rhizobia is clearly the product of past coevolution. However, the nature of ongoing evolution between these partners is less clear. To characterize the nature of recent coevolution between legumes and rhizobia, we used population genomic analysis to characterize selection on functionally annotated symbiosis genes as well as on symbiosis gene candidates identified through a two-species association analysis. For the association analysis, we inoculated each of 202 accessions of the legume host Medicago truncatula with a community of 88 Sinorhizobia (Ensifer) meliloti strains. Multistrain inoculation, which better reflects the ecological reality of rhizobial selection in nature than single-strain inoculation, allows strains to compete for nodulation opportunities and host resources and for hosts to preferentially form nodules and provide resources to some strains. We found extensive host by symbiont, that is, genotype-by-genotype, effects on rhizobial fitness and some annotated rhizobial genes bear signatures of recent positive selection. However, neither genes responsible for this variation nor annotated host symbiosis genes are enriched for signatures of either positive or balancing selection. This result suggests that stabilizing selection dominates selection acting on symbiotic traits and that variation in these traits is under mutation-selection balance. Consistent with the lack of positive selection acting on host genes, we found that among-host variation in growth was similar whether plants were grown with rhizobia or N-fertilizer, suggesting that the symbiosis may not be a major driver of variation in plant growth in multistrain contexts.

Priority effects alter the colonization success of a host-associated parasite and mutualist

Priority effects shape the assembly of free-living communities and host-associated communities. However, the current literature does not fully incorporate two features of host-symbiont interactions-correlated host responses to multiple symbionts and ontogenetic changes in host responses to symbionts-leading to an incomplete picture of the role of priority effects in host-associated communities. We factorially manipulated the inoculation timing of two plant symbionts (mutualistic rhizobia bacteria and parasitic root-knot nematodes) and tested how host age at arrival, arrival order, and arrival synchrony affected symbiont colonization success in the model legume Medicago truncatula. We found that host age, arrival order, and arrival synchrony significantly affected colonization of one or both symbionts. Host age at arrival only affected nematodes but not rhizobia: younger plants were more heavily infected than older plants. By contrast, arrival order only affected rhizobia but not nematodes: plants formed more rhizobia nodules when rhizobia arrived before nematodes. Finally, synchronous arrival decreased colonization both symbionts, an effect that depended on host age. Our results demonstrate that priority effects compromise the host's ability to control colonization by two major symbionts, and suggest that the role of correlated host responses and host ontogeny in the assembly of host-associated communities deserve further attention.

Amino Acid Polymorphisms in the VHIID Conserved Motif of Nodulation Signaling Pathways 2 Distinctly Modulate Symbiotic Signaling and Nodule Morphogenesis in Medicago truncatula

Legumes establish an endosymbiotic association with nitrogen-fixing soil bacteria. Following the mutual recognition of the symbiotic partner, the infection process is controlled by the induction of the signaling pathway and subsequent activation of symbiosis-related host genes. One of the protein complexes regulating nitrogen-fixing root nodule symbiosis is formed by GRAS domain regulatory proteins Nodulation Signaling Pathways 1 and 2 (NSP1 and NSP2) that control the expression of several early nodulation genes. Here, we report on a novel point mutant allele (nsp2-6) affecting the function of the NSP2 gene and compared the mutant with the formerly identified nsp2-3 mutant. Both mutants carry a single amino acid substitution in the VHIID motif of the NSP2 protein. We found that the two mutant alleles show dissimilar root hair response to bacterial infection. Although the nsp2-3 mutant developed aberrant infection threads, rhizobia were able to colonize nodule cells in this mutant. The encoded NSP2 proteins of the nsp2-3 and the novel nsp2 mutants interact with NSP1 diversely and, as a consequence, the activation of early nodulin genes and nodule organogenesis are arrested in the new nsp2 allele. The novel mutant with amino acid substitution D244H in NSP2 shows similar defects in symbiotic responses as a formerly identified nsp2-2 mutant carrying a deletion in the NSP2 gene. Additionally, we found that rhizobial strains induce delayed nodule formation on the roots of the ns2-3 weak allele. Our study highlights the importance of a conserved Asp residue in the VHIID motif of NSP2 that is required for the formation of a functional NSP1-NSP2 signaling module. Furthermore, our results imply the involvement of NSP2 during differentiation of symbiotic nodule cells.

Ophryocystis sitonae sp. nov., (Neogregarinida: Ophryocystidae) parasitizing Sitona humeralis Stephens, 1831 (Coleoptera: Curculionidae)

Sitona humeralis Stephens 1831 (Coleoptera: Curculionidae) is an important pest of the Medicago and Vicia species in Turkey, and this study was conducted the determine the natural pathogens of this beetle. In the present study, a new neogregarine was observed in Malpighian tubules of the S. humeralis, collected from Ordu (Turkey) on the wild Medicago species. The yellowish oocysts were the most notable feature of the current neogregarine. The Giemsa-stained mature oocysts were fusiform shaped and measured 8.7 ± 0.7 (7.12-11.11) μm in length and 4.1 ± 0.3 (3.05-5.01) μm in width. The smooth oocyst wall was relatively thin (175-230 nm), and polar plugs were non-evident (weight = 380 nm, height = 500 nm). The 18S rDNA gene of the current neogregarine was sequenced and compared with fifteen sequences from GenBank. Morphological, ultrastructural, and molecular features indicate that the described neogregarine in S. humeralis differed from the all known Ophryocystis species and named here Ophryocystis sitonae sp. nov.

Evolutionary and expression dynamics of LRR-RLKs and functional establishment of KLAVIER homolog in shoot mediated regulation of AON in chickpea symbiosis

Chickpea shoot exogenously treated with cytokinin showed stunted phenotype of root, shoot and significantly reduced nodule numbers. Genome-wide identification of LRR-RLKs in chickpea and Medicago resulted in 200 and 371 genes respectively. Gene duplication analysis revealed that LRR-RLKs family expanded through segmental duplications in chickpea and tandem duplications in Medicago. Expression profiling of LRR-RLKs revealed their involvement in cytokinin signaling and plant organ development. Overexpression of KLAVIER ortholog of chickpea, Ca_LRR-RLK147, in roots revealed its localization in the membrane but showed no effect on root nodulation despite increased cle peptide levels. Two findings (i) drastic effect on nodule number by exogenous cytokinin treatment to only shoot and restoration to normal nodulation by treatment to both root and shoot tissue and (ii) no effect on nodule number by overexpression of Ca_LRR-RLK147 establishes the fact that despite presence of cle peptides in root, the function of Ca_LRR-RLK147 was shoot mediated during AON.

QTL dissection and mining of candidate genes for Ascochyta fabae and Orobanche crenata resistance in faba bean (Vicia faba L.)

BACKGROUND

Ascochyta blight caused by Ascochyta fabae Speg. and broomrape (Orobanche crenata) are among the economically most significant pathogens of faba bean. Several QTLs conferring resistance against the two pathogens have been identified and validated in different genetic backgrounds. The aim of this study was to saturate the most stable QTLs for ascochyta and broomrape resistance in two Recombinant Inbred Line (RIL) populations, 29H x Vf136 and Vf6 x Vf136, to identify candidate genes conferring resistance against these two pathogens.

RESULTS

We exploited the synteny between faba bean and the model species Medicago truncatula by selecting a set of 219 genes encoding putative WRKY transcription factors and defense related proteins falling within the target QTL intervals, for genotyping and marker saturation in the two RIL populations. Seventy and 50 of the candidate genes could be mapped in 29H x Vf136 and Vf6 x Vf136, respectively. Besides the strong reduction of the QTL intervals, the mapping process allowed replacing previous dominant and pedigree-specific RAPD flanking markers with robust and transferrable SNP markers, revealing promising candidates for resistance against the two pathogens.

CONCLUSIONS

Although further efforts in association mapping and expression studies will be required to corroborate the candidate genes for resistance, the fine-mapping approach proposed here increases the genetic resolution of relevant QTL regions and paves the way for an efficient deployment of useful alleles for faba bean ascochyta and broomrape resistance through marker-assisted breeding.

MtExpress, a Comprehensive and Curated RNAseq-based Gene Expression Atlas for the Model Legume Medicago truncatula

Although RNA sequencing (RNAseq) has been becoming the main transcriptomic approach in the model legume Medicago truncatula, there is currently no genome-wide gene expression atlas covering the whole set of RNAseq data published for this species. Nowadays, such a tool is highly valuable to provide a global view of gene expression in a wide range of conditions and tissues/organs. Here, we present MtExpress, a gene expression atlas that compiles an exhaustive set of published M. truncatula RNAseq data (https://medicago.toulouse.inrae.fr/MtExpress). MtExpress makes use of recent releases of M. truncatula genome sequence and annotation, as well as up-to-date tools to perform mapping, quality control, statistical analysis and normalization of RNAseq data. MtExpress combines semi-automated pipelines with manual re-labeling and organization of samples to produce an attractive and user-friendly interface, fully integrated with other available Medicago genomic resources. Importantly, MtExpress is highly flexible, in terms of both queries, e.g. allowing searches with gene names and orthologous gene IDs from Arabidopsis and other legume species, and outputs, to customize visualization and redirect gene study to relevant Medicago webservers. Thanks to its semi-automated pipeline, MtExpress will be frequently updated to follow the rapid pace of M. truncatula RNAseq data publications, as well as the constant improvement of genome annotation. MtExpress also hosts legacy GeneChip expression data originally stored in the Medicago Gene Expression Atlas, as a very valuable and complementary resource.

Extensive genomic rearrangements mediated by repetitive sequences in plastomes of Medicago and its relatives

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.

First Report of Aphanomyces euteiches Race 1 and Race 2 Causing Aphanomyces Root Rot on Alfalfa (Medicago sativa) in South Dakota

Aphanomyces euteiches causes Aphanomyces root rot (ARR) in alfalfa (Medicago sativa), along with root rot on many other legumes, including pea, clover, and lentil (Malvick et al., 2009). In 2020, South Dakota (SD) planted the most acres of alfalfa in the United States, which demonstrates the importance of alfalfa to the state. Several SD growers reported alfalfa establishment problems likely to be associated with ARR. Soil samples were collected from 16 fields under commercial alfalfa production in Lake County, SD in June 2020. Composite soil samples based on 24 subsamples were collected in a W-shaped pattern at a depth of 15 cm. Collected soil was sieved, and 80 cm3 was placed in plastic pots (6 cm x 6 cm). Each pot was planted with 25 seeds, covered with an additional 15 cm3 soil, and placed in a growth chamber with a 16-hour photoperiod at temperatures of 24 and 19 ℃ (day and night). Alfalfa seedlings, including Saranac (susceptible to races R1 and R2), WAPH-1 (resistant only to R1), WAPH-5 (resistant to both R1 and R2), and Mustang 625 (resistant to both R1 and R2 and coated with mefenoxam) grew in collected soil for 7 days, followed by 4 days under flooded conditions. Trays were drained, and at 21 days after planting (DAP), roots were removed from soil, washed in distilled water, and rated to measure severity of disease symptoms (Samac et al., 2015). The average severity index (ASI) used a 1-5 disease severity scale, 5 being a dead plant and 1 being no symptoms present (http://www.naaic.org/stdtests/Aphano.html). Race was based on ASI where R1 included an ASI of ≥3 for Saranac and <3 for WAPH-1, and R2 included an ASI of >3.0 for Saranac and WAPH-1 and <3.0 for WAPH-5 (Malvick and Grau, 2001). Race-typing experiments were repeated twice with six replicate pots per alfalfa cultivar per experiment and determined the presence of both R1 and R2 in Lake County, SD. ASI values for Mustang 625 and WAPH-5 were similar across all fields evaluated, which indicates limited confounding effects of other root rotting pathogens. DNA was extracted from three symptomatic roots from each field and was PCR amplified using A. euteiches specific primers (Vandemark et al., 2002). A PCR product was observed in all 16 fields evaluated, and the absence of a product was observed when DNA was extracted from alfalfa roots grown in vermiculite. Following race-typing, infected alfalfa roots were surfaced sterilized and placed on Aphanomyces selective media consisting of mefenoxam and benomyl in cornmeal agar (CMA) (Pfender et al., 1984). Isolates were identified as A. euteiches based on hyphal morphology (Malvick and Grau, 2001). Alfalfa seedlings (Saranac) were grown in vermiculite under growth conditions used for the race-typing assay and inoculated 6 DAP with two isolates of A. euteiches. Inoculation was completed using half plates of one week old A. euteiches mycelium on CMA blended with one liter of water (Samac et al., 2015). At 35 DAP, control alfalfa seedlings inoculated with blended CMA and water remained asymptomatic, and alfalfa infected with A. euteiches displayed symptoms including honey-brown colored lesions. For confirmation of Koch's postulates, DNA from three re-infected seedlings was again PCR amplified using A. euteiches specific primers and confirmed our previous work. This is the first report of either R1 or R2 of A. euteiches causing ARR on alfalfa in SD. To avoid future yield loss, SD growers should consider planting available alfalfa cultivars that have resistance to both races of A. euteiches.

Three Common Symbiotic ABC Subfamily B Transporters in Medicago truncatula Are Regulated by a NIN-Independent Branch of the Symbiosis Signaling Pathway

Several ATP-binding cassette (ABC) transporters involved in the arbuscular mycorrhizal symbiosis and nodulation have been identified. We describe three previously unreported ABC subfamily B transporters, named AMN1, AMN2, and AMN3 (ABCB for mycorrhization and nodulation), that are expressed early during infection by rhizobia and arbuscular mycorrhizal fungi. These ABCB transporters are strongly expressed in symbiotically infected tissues, including in root-hair cells with rhizobial infection threads and arbusculated cells. During nodulation, the expression of these genes is highly induced by rhizobia and purified Nod factors and is dependent on DMI3 but is not dependent on other known major regulators of infection, such as NIN, NSP1, or NSP2. During mycorrhization their expression is dependent on DMI3 and RAM1 but not on NSP1 and NSP2. Therefore, they may be commonly regulated through a distinct branch of the common symbiotic pathway. Mutants with exonic Tnt1-transposon insertions were isolated for all three genes. None of the single or double mutants showed any differences in colonization by either rhizobia or mycorrhizal fungi, but the triple amn1 amn2 amn3 mutant showed an increase in nodule number. Further studies are needed to identify potential substrates of these transporters and understand their roles in these beneficial symbioses.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Cleaning the Medicago Microarray Database to Improve Gene Function Analysis

Transcriptomics studies have been facilitated by the development of microarray and RNA-Seq technologies, with thousands of expression datasets available for many species. However, the quality of data can be highly variable, making the combined analysis of different datasets difficult and unreliable. Most of the microarray data for Medicago truncatula, the barrel medic, have been stored and made publicly accessible on the web database Medicago truncatula Gene Expression atlas (MtGEA). The aim of this work is to ameliorate the quality of the MtGEA database through a general method based on logical and statistical relationships among parameters and conditions. The initial 716 columns available in the dataset were reduced to 607 by evaluating the quality of data through the sum of the expression levels over the entire transcriptome probes and Pearson correlation among hybridizations. The reduced dataset shows great improvements in the consistency of the data, with a reduction in both false positives and false negatives resulting from Pearson correlation and GO enrichment analysis among genes. The approach we used is of general validity and our intent is to extend the analysis to other plant microarray databases.

Multiple Mutualism Effects generate synergistic selection and strengthen fitness alignment in the interaction between legumes, rhizobia and mycorrhizal fungi

Nearly all organisms participate in multiple mutualisms, and complementarity within these complex interactions can result in synergistic fitness effects. However, it remains largely untested how multiple mutualisms impact eco-evolutionary dynamics in interacting species. We tested how multiple microbial mutualists-N-fixing bacteria and mycorrrhizal fungi-affected selection and heritability of traits in their shared host plant (Medicago truncatula), as well as fitness alignment between partners. Our results demonstrate for the first time that multiple mutualisms synergistically affect the selection and heritability of host traits and enhance fitness alignment between mutualists. Specifically, we found interaction with multiple microbial symbionts doubled the strength of natural selection on a plant architectural trait, resulted in 2- to 3-fold higher heritability of plant reproductive success, and more than doubled fitness alignment between N-fixing bacteria and plants. These findings show synergism generated by multiple mutualisms extends to key components of microevolutionary change, emphasising the importance of multiple mutualism effects on evolutionary trajectories.

Advancing knowledge of the plant nuclear periphery and its application for crop science

In this review, we explore recent advances in knowledge of the structure and dynamics of the plant nuclear envelope. As a paradigm, we focused our attention on the Linker of Nucleoskeleton and Cytoskeleton (LINC) complex, a structurally conserved bridging complex comprising SUN domain proteins in the inner nuclear membrane and KASH domain proteins in the outer nuclear membrane. Studies have revealed that this bridging complex has multiple functions with structural roles in positioning the nucleus within the cell, conveying signals across the membrane and organizing chromatin in the 3D nuclear space with impact on gene transcription. We also provide an up-to-date survey in nuclear dynamics research achieved so far in the model plant Arabidopsis thaliana that highlights its potential impact on several key plant functions such as growth, seed maturation and germination, reproduction and response to biotic and abiotic stress. Finally, we bring evidences that most of the constituents of the LINC Complex and associated components are, with some specificities, conserved in monocot and dicot crop species and are displaying very similar functions to those described for Arabidopsis. This leads us to suggest that a better knowledge of this system and a better account of its potential applications will in the future enhance the resilience and productivity of crop plants.

The two-component system ActJK is involved in acid stress tolerance and symbiosis in Sinorhizobium meliloti

The nitrogen-fixing α-proteobacterium Sinorhizobium meliloti genome codifies at least 50 response regulator (RR) proteins mediating different and, in many cases, unknown processes. RR-mutant library screening allowed us to identify genes potentially implicated in survival to acid conditions. actJ mutation resulted in a strain with reduced growth rate under mildly acidic conditions as well as a lower capacity to tolerate a sudden shift to lethal acidic conditions compared with the parental strain. Mutation of the downstream gene actK, which encodes for a histidine kinase, showed a similar phenotype in acidic environments suggesting a functional two-component system. Interestingly, even though nodulation kinetics, quantity, and macroscopic morphology of Medicago sativa nodules were not affected in actJ and actK mutants, ActK was required to express the wild-type nitrogen fixation phenotype and ActJK was necessary for full bacteroid development and nodule occupancy. The actJK regulatory system presented here provides insights into an evolutionary process in rhizobium adaptation to acidic environments and suggests that actJK-controlled functions are crucial for optimal symbiosis development.

Phylogeny and Species Delimitation of Chinese Medicago (Leguminosae) and Its Relatives Based on Molecular and Morphological Evidence

Medicago and its relatives, Trigonella and Melilotus comprise the most important forage resources globally. The alfalfa selected from the wild relatives has been cultivated worldwide as the forage queen. In the Flora of China, 15 Medicago, eight Trigonella, and four Melilotus species are recorded, of which six Medicago and two Trigonella species are introduced. Although several studies have been conducted to investigate the phylogenetic relationship within the three genera, many Chinese naturally distributed or endemic species are not included in those studies. Therefore, the taxonomic identity and phylogenetic relationship of these species remains unclear. In this study, we collected samples representing 18 out of 19 Chinese naturally distributed species of these three genera and three introduced Medicago species, and applied an integrative approach by combining evidences from population-based morphological clusters and molecular data to investigate species boundaries. A total of 186 individuals selected from 156 populations and 454 individuals from 124 populations were collected for genetic and morphological analyses, respectively. We sequenced three commonly used DNA barcodes (trnH-psbA, trnK-matK, and ITS) and one nuclear marker (GA3ox1) for phylogenetic analyses. We found that 16 out of 21 species could be well delimited based on phylogenetic analyses and morphological clusters. Two Trigonella species may be merged as one species or treated as two subspecies, and Medicago falcata should be treated as a subspecies of the M. sativa complex. We further found that major incongruences between the chloroplast and nuclear trees mainly occurred among the deep diverging lineages, which may be resulted from hybridization, incomplete lineage sorting and/or sampling errors. Further studies involving a finer sampling of species associated with large scale genomic data should be employed to better understand the species delimitation of these three genera.

Medicago sativa and Medicago truncatula Show Contrasting Root Metabolic Responses to Drought

Drought is an environmental stressor that affects crop yield worldwide. Understanding plant physiological responses to stress conditions is needed to secure food in future climate conditions. In this study, we applied a combination of plant physiology and metabolomic techniques to understand plant responses to progressive water deficit focusing on the root system. We chose two legume plants with contrasting tolerance to drought, the widely cultivated alfalfa Medicago sativa (Ms) and the model legume Medicago truncatula (Mt) for comparative analysis. Ms taproot (tapR) and Mt fibrous root (fibR) biomass increased during drought, while a progressive decline in water content was observed in both species. Metabolomic analysis allowed the identification of key metabolites in the different tissues tested. Under drought, carbohydrates, abscisic acid, and proline predominantly accumulated in leaves and tapRs, whereas flavonoids increased in fibRs in both species. Raffinose-family related metabolites accumulated during drought. Along with an accumulation of root sucrose in plants subjected to drought, both species showed a decrease in sucrose synthase (SUS) activity related to a reduction in the transcript level of SUS1, the main SUS gene. This study highlights the relevance of root carbon metabolism during drought conditions and provides evidence on the specific accumulation of metabolites throughout the root system.

Medicago PHYA promotes flowering, primary stem elongation and expression of flowering time genes in long days

BACKGROUND

Flowering time is an important trait for productivity in legumes, which include many food and fodder plants. Medicago truncatula (Medicago) is a model temperate legume used to study flowering time pathways. Like Arabidopsis thaliana (Arabidopsis), its flowering is promoted by extended periods of cold (vernalization, V), followed by warm long day (LD) photoperiods. However, Arabidopsis flowering-time genes such as the FLOWERING LOCUS C (FLC)/ MADS AFFECTING FLOWERING (MAF) clade are missing and CONSTANS-LIKE (CO-LIKE) genes do not appear to have a role in Medicago or Pisum sativum (pea). Another photoperiodic regulator, the red/far red photoreceptor PHYTOCHROME A (PHYA), promotes Arabidopsis flowering by stabilizing the CO protein in LD. Interestingly, despite the absence of CO-LIKE function in pea, PsPHYA plays a key role in promoting LD photoperiodic flowering and plant architecture. Medicago has one homolog of PHYA, MtPHYA, but its function is not known.

RESULTS

Genetic analysis of two MtPHYA Tnt1 insertion mutant alleles indicates that MtPHYA has an important role in promoting Medicago flowering and primary stem elongation in VLD and LD and in perception of far-red wavelengths in seedlings. MtPHYA positively regulates the expression of MtE1-like (MtE1L), a homologue of an important legume-specific flowering time gene, E1 in soybean and other Medicago LD-regulated flowering-time gene homologues, including the three FLOWERING LOCUS T-LIKE (FT-LIKE) genes, MtFTa1, MtFTb1 and MtFTb2 and the two FRUITFULL-LIKE (FUL-LIKE) genes MtFULa and MtFULb. MtPHYA also modulates the expression of the circadian clock genes, GIGANTEA (GI) and TIMING OF CAB EXPRESSION 1a (TOC1a). Genetic analyses indicate that Mtphya-1 Mte1l double mutants flowered at the same time as the single mutants. However, Mtphya-1 Mtfta1 double mutants had a weak additive effect in delaying flowering and in reduction of primary axis lengths beyond what was conferred by either of the single mutants.

CONCLUSION

MtPHYA has an important role in LD photoperiodic control of flowering, plant architecture and seedling de-etiolation under far-red wavelengths in Medicago. It promotes the expression of LD-induced flowering time genes and modulates clock-related genes. In addition to MtFTa1, MtPHYA likely regulates other targets during LD floral induction in Medicago.

Plant defensin antibacterial mode of action against Pseudomonas species

BACKGROUND

Though many plant defensins exhibit antibacterial activity, little is known about their antibacterial mode of action (MOA). Antimicrobial peptides with a characterized MOA induce the expression of multiple bacterial outer membrane modifications, which are required for resistance to these membrane-targeting peptides. Mini-Tn5-lux mutant strains of Pseudomonas aeruginosa with Tn insertions disrupting outer membrane protective modifications were assessed for sensitivity against plant defensin peptides. These transcriptional lux reporter strains were also evaluated for lux gene expression in response to sublethal plant defensin exposure. Also, a plant pathogen, Pseudomonas syringae pv. syringae was modified through transposon mutagenesis to create mutants that are resistant to in vitro MtDef4 treatments.

RESULTS

Plant defensins displayed specific and potent antibacterial activity against strains of P. aeruginosa. A defensin from Medicago truncatula, MtDef4, induced dose-dependent gene expression of the aminoarabinose modification of LPS and surface polycation spermidine production operons. The ability for MtDef4 to damage bacterial outer membranes was also verified visually through fluorescent microscopy. Another defensin from M. truncatula, MtDef5, failed to induce lux gene expression and limited outer membrane damage was detected with fluorescent microscopy. The transposon insertion site on MtDef4 resistant P. syringae pv. syringae mutants was sequenced, and modifications of ribosomal genes were identified to contribute to enhanced resistance to plant defensin treatments.

CONCLUSIONS

MtDef4 damages the outer membrane similar to polymyxin B, which stimulates antimicrobial peptide resistance mechanisms to plant defensins. MtDef5, appears to have a different antibacterial MOA. Additionally, the MtDef4 antibacterial mode of action may also involve inhibition of translation.

Identification of nodulation-related genes in Medicago truncatula using genome-wide association studies and co-expression networks

Genome-wide association studies (GWAS) have proven to be a valuable approach for identifying genetic intervals associated with phenotypic variation in Medicago truncatula. These intervals can vary in size, depending on the historical local recombination. Typically, significant intervals span numerous gene models, limiting the ability to resolve high-confidence candidate genes underlying the trait of interest. Additional genomic data, including gene co-expression networks, can be combined with the genetic mapping information to successfully identify candidate genes. Co-expression network analysis provides information about the functional relationships of each gene through its similarity of expression patterns to other well-defined clusters of genes. In this study, we integrated data from GWAS and co-expression networks to pinpoint candidate genes that may be associated with nodule-related phenotypes in M. truncatula. We further investigated a subset of these genes and confirmed that several had existing evidence linking them nodulation, including MEDTR2G101090 (PEN3-like), a previously validated gene associated with nodule number.

Physical Dormancy Release in Medicago truncatula Seeds Is Related to Environmental Variations

Seed dormancy and timing of its release is an important developmental transition determining the survival of individuals, populations, and species in variable environments. Medicago truncatula was used as a model to study physical seed dormancy at the ecological and genetics level. The effect of alternating temperatures, as one of the causes releasing physical seed dormancy, was tested in 178 M. truncatula accessions over three years. Several coefficients of dormancy release were related to environmental variables. Dormancy varied greatly (4-100%) across accessions as well as year of experiment. We observed overall higher physical dormancy release under more alternating temperatures (35/15 °C) in comparison with less alternating ones (25/15 °C). Accessions from more arid climates released dormancy under higher experimental temperature alternations more than accessions originating from less arid environments. The plasticity of physical dormancy can probably distribute the germination through the year and act as a bet-hedging strategy in arid environments. On the other hand, a slight increase in physical dormancy was observed in accessions from environments with higher among-season temperature variation. Genome-wide association analysis identified 136 candidate genes related to secondary metabolite synthesis, hormone regulation, and modification of the cell wall. The activity of these genes might mediate seed coat permeability and, ultimately, imbibition and germination.

Activity of Saponins from Medicago Species against Phytoparasitic Nematodes

Content of bioactive saponins of Medicago species suggests that they may also exert, as previously demonstrated on M. sativa, nematicidal properties exploitable for the formulation of new products for sustainable phytoparasitic nematode management. This study was addressed to highlight the bioactivity of saponins from five different Medicago species still poorly known for their biological efficacy, i.e., M. heyniana, M. hybrida, M. lupulina, M. murex and M. truncatula, against the plant parasitic nematodes Meloidogyne incognita, Xiphinema index and Globodera rostochiensis. The bioactivity of the extracts from the five Medicago species was assessed by in vitro assays on the juveniles (J2) and eggs of M. incognita and G. rostochiensis and the adult females of X. index. The suppressiveness to M. incognita of soil treatments with the Medicago plant biomasses was also investigated in a tomato experiment. The nematicidal activity of the five Medicago species was reported and discussed in relation to their phytochemical profile.

PRPs localized to the middle lamellae are required for cortical tissue integrity in Medicago truncatula roots

A family of repetitive proline-rich proteins interact with acidic pectins and play distinct roles in legume root cell walls affecting cortical and vascular structure. A proline-rich protein (PRP) family, composed of tandemly repeated Pro-Hyp-Val-X-Lys pentapeptide motifs, is found primarily in the Leguminosae. Four distinct size classes within this family are encoded by seven tightly linked genes: MtPRP1, MtPRP2 and MtPRP3, and four nearly identical MtPRP4 genes. Promoter fusions to β-glucuronidase showed strong expression in the stele of hairy roots for all 4 PRP genes tested, with additional expression in the cortex for PRP1, PRP2 and PRP4. All except MtPRP4 are strongly expressed in non-tumorous roots, and secreted and ionically bound to root cell walls. These PRPs are absent from root epidermal cell walls, and PRP accumulation is highly localized within the walls of root cortical and vascular tissues. Within xylem tissue, PRPs are deposited in secondary thickenings where it is spatially exclusive to lignin. In newly differentiating xylem, PRPs are deposited in the regularly spaced paired-pits and pit membranes that hydraulically connect neighboring xylem elements. Hairpin-RNA knock-down constructs reducing PRP expression in Medicago truncatula hairy root tumors disrupted cortical and vascular patterning. Immunoblots showed that the knockdown tumors had potentially compensating increases in the non-targeted PRPs, all of which cross-react with the anti-PRP antibodies. However, PRP3 knockdown differed from knockdown of PRP1 and PRP2 in that it greatly reduced viability of hairy root tumors. We hypothesize that repetitive PRPs interact with acidic pectins to form block-copolymer gels that can play distinct roles in legume root cell walls.

Using clear plastic CD cases as low-cost mini-rhizotrons to phenotype root traits

PREMISE

We developed a novel low-cost method to visually phenotype belowground structures in the plant rhizosphere. We devised the method introduced here to address the difficulties encountered growing plants in seed germination pouches for long-term experiments and the high cost of other mini-rhizotron alternatives.

METHODS AND RESULTS

The method described here took inspiration from homemade ant farms commonly used as an educational tool in elementary schools. Using compact disc (CD) cases, we developed mini-rhizotrons for use in the field and laboratory using the burclover Medicago lupulina.

CONCLUSIONS

Our method combines the benefits of pots and germination pouches. In CD mini-rhizotrons, plants grew significantly larger than in germination pouches, and unlike pots, it is possible to measure roots without destructive sampling. Our protocol is a cheaper, widely available alternative to more destructive methods, which could facilitate the study of belowground phenotypes and processes by scientists with fewer resources.

SNARE Complexity in Arbuscular Mycorrhizal Symbiosis

How cells control the proper delivery of vesicles and their associated cargo to specific plasma membrane (PM) domains upon internal or external cues is a major question in plant cell biology. A widely held hypothesis is that expansion of plant exocytotic machinery components, such as SNARE proteins, has led to a diversification of exocytotic membrane trafficking pathways to function in specific biological processes. A key biological process that involves the creation of a specialized PM domain is the formation of a host-microbe interface (the peri-arbuscular membrane) in the symbiosis with arbuscular mycorrhizal fungi. We have previously shown that the ability to intracellularly host AM fungi correlates with the evolutionary expansion of both v- (VAMP721d/e) and t-SNARE (SYP132α) proteins, that are essential for arbuscule formation in Medicago truncatula. Here we studied to what extent the symbiotic SNAREs are different from their non-symbiotic family members and whether symbiotic SNAREs define a distinct symbiotic membrane trafficking pathway. We show that all tested SYP1 family proteins, and most of the non-symbiotic VAMP72 members, are able to complement the defect in arbuscule formation upon knock-down/-out of their symbiotic counterparts when expressed at sufficient levels. This functional redundancy is in line with the ability of all tested v- and t-SNARE combinations to form SNARE complexes. Interestingly, the symbiotic t-SNARE SYP132α appeared to occur less in complex with v-SNAREs compared to the non-symbiotic syntaxins in arbuscule-containing cells. This correlated with a preferential localization of SYP132α to functional branches of partially collapsing arbuscules, while non-symbiotic syntaxins accumulate at the degrading parts. Overexpression of VAMP721e caused a shift in SYP132α localization toward the degrading parts, suggesting an influence on its endocytic turn-over. These data indicate that the symbiotic SNAREs do not selectively interact to define a symbiotic vesicle trafficking pathway, but that symbiotic SNARE complexes are more rapidly disassembled resulting in a preferential localization of SYP132α at functional arbuscule branches.

Overexpression of Medicago MtCDFd1_1 Causes Delayed Flowering in Medicago via Repression of MtFTa1 but Not MtCO-Like Genes

Optimizing flowering time is crucial for maximizing crop productivity, but gaps remain in the knowledge of the mechanisms underpinning temperate legume flowering. Medicago, like winter annual Arabidopsis, accelerates flowering after exposure to extended cold (vernalization, V) followed by long-day (LD) photoperiods. In Arabidopsis, photoperiodic flowering is triggered through CO, a photoperiodic switch that directly activates the FT gene encoding a mobile florigen and potent activator of flowering. In Arabidopsis, several CYCLING DOF FACTORs (CDFs), including AtCDF1, act redundantly to repress CO and thus FT expression, until their removal in LD by a blue-light-induced F-BOX1/GIGANTEA (FKF1/GI) complex. Medicago possesses a homolog of FT, MtFTa1, which acts as a strong activator of flowering. However, the regulation of MtFTa1 does not appear to involve a CO-like gene. Nevertheless, work in pea suggests that CDFs may still regulate flowering time in temperate legumes. Here, we analyze the function of Medicago MtCDF genes with a focus on MtCDFd1_1 in flowering time and development. MtCDFd1_1 causes strong delays to flowering when overexpressed in Arabidopsis and shows a cyclical diurnal expression in Medicago with peak expression at dawn, consistent with AtCDF genes like AtCDF1. However, MtCDFd1_1 lacks predicted GI or FKF1 binding domains, indicating possible differences in its regulation from AtCDF1. In Arabidopsis, CDFs act in a redundant manner, and the same is likely true of temperate legumes as no flowering time phenotypes were observed when MtCDFd1_1 or other MtCDFs were knocked out in Medicago Tnt1 lines. Nevertheless, overexpression of MtCDFd1_1 in Medicago plants resulted in late flowering relative to wild type in inductive vernalized long-day (VLD) conditions, but not in vernalized short days (VSDs), rendering them day neutral. Expression of MtCO-like genes was not affected in the transgenic lines, but LD-induced genes MtFTa1, MtFTb1, MtFTb2, and MtSOC1a showed reduced expression. Plants carrying both the Mtfta1 mutation and 35S:MtCDFd1_1 flowered no later than the Mtfta1 plants. This indicates that 35S:MtCDFd1_1 likely influences flowering in VLD via repressive effects on MtFTa1 expression. Overall, our study implicates MtCDF genes in photoperiodic regulation in Medicago by working redundantly to repress FT-like genes, particularly MtFTa1, but in a CO-independent manner, indicating differences from the Arabidopsis model.

Lost and Found: Return of the Inverted Repeat in the Legume Clade Defined by Its Absence

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.

Regulation of Resource Partitioning Coordinates Nitrogen and Rhizobia Responses and Autoregulation of Nodulation in Medicago truncatula

Understanding how plants respond to nitrogen in their environment is crucial for determining how they use it and how the nitrogen use affects other processes related to plant growth and development. Under nitrogen limitation the activity and affinity of uptake systems is increased in roots, and lateral root formation is regulated in order to adapt to low nitrogen levels and scavenge from the soil. Plants in the legume family can form associations with rhizobial nitrogen-fixing bacteria, and this association is tightly regulated by nitrogen levels. The effect of nitrogen on nodulation has been extensively investigated, but the effects of nodulation on plant nitrogen responses remain largely unclear. In this study, we integrated molecular and phenotypic data in the legume Medicago truncatula and determined that genes controlling nitrogen influx are differently expressed depending on whether plants are mock or rhizobia inoculated. We found that a functional autoregulation of nodulation pathway is required for roots to perceive, take up, and mobilize nitrogen as well as for normal root development. Our results together revealed that autoregulation of nodulation, root development, and the location of nitrogen are processes balanced by the whole plant system as part of a resource-partitioning mechanism.

[Proteomic Profile of the Bacterium Sinorhizobium meliloti Depends on Its Life Form and Host Plant Species]

The importance of root nodule bacteria in biotechnology is determined by their distinctive feature: symbiotic nitrogen fixation resulting in the production of organic nitrogen-containing compounds. While interacting with host legume plants, the cells of these bacteria undergo global changes at all levels of expression of genetic information leading to the formation in root nodules of so-called bacteroids functioning as nitrogen fixation factories. The molecular mechanisms underlying plant-microbial symbiosis are actively investigated, and one of the most interesting and poorly studied aspects of this problem is the species-specificity of interaction between root nodule bacteria and host plants. In this work we have performed the proteomic analysis of the Sinorhizobium meliloti bacteroids isolated from two legume species: alfalfa (Medicago sativa L.) and yellow sweet clover (Melilotus officinalis L.). It has been shown that the S. meliloti bacteroids produce a lot of proteins (many of them associated with symbiosis) in a host-specific manner, i.e., only in certain host plant species. It has been demonstrated for the first time that the levels of expression in bacteroids of the genes encoding the ExoZ and MscL proteins responsible for the synthesis of surface lipopolysaccha-rides and formation of a large conductance mechanosensitive channel, respectively, depend on a host plant species that confirms the results of proteomic analysis. Overall, our data show that the regulation of bacteroid development by the host plant has species-specific features.

Genome-Wide Association Analyses in the Model Rhizobium Ensifer meliloti

Genome-wide association studies (GWAS) can identify genetic variants responsible for naturally occurring and quantitative phenotypic variation. Association studies therefore provide a powerful complement to approaches that rely on de novo mutations for characterizing gene function. Although bacteria should be amenable to GWAS, few GWAS have been conducted on bacteria, and the extent to which nonindependence among genomic variants (e.g., linkage disequilibrium [LD]) and the genetic architecture of phenotypic traits will affect GWAS performance is unclear. We apply association analyses to identify candidate genes underlying variation in 20 biochemical, growth, and symbiotic phenotypes among 153 strains of Ensifer meliloti For 11 traits, we find genotype-phenotype associations that are stronger than expected by chance, with the candidates in relatively small linkage groups, indicating that LD does not preclude resolving association candidates to relatively small genomic regions. The significant candidates show an enrichment for nucleotide polymorphisms (SNPs) over gene presence-absence variation (PAV), and for five traits, candidates are enriched in large linkage groups, a possible signature of epistasis. Many of the variants most strongly associated with symbiosis phenotypes were in genes previously identified as being involved in nitrogen fixation or nodulation. For other traits, apparently strong associations were not stronger than the range of associations detected in permuted data. In sum, our data show that GWAS in bacteria may be a powerful tool for characterizing genetic architecture and identifying genes responsible for phenotypic variation. However, careful evaluation of candidates is necessary to avoid false signals of association.IMPORTANCE Genome-wide association analyses are a powerful approach for identifying gene function. These analyses are becoming commonplace in studies of humans, domesticated animals, and crop plants but have rarely been conducted in bacteria. We applied association analyses to 20 traits measured in Ensifer meliloti, an agriculturally and ecologically important bacterium because it fixes nitrogen when in symbiosis with leguminous plants. We identified candidate alleles and gene presence-absence variants underlying variation in symbiosis traits, antibiotic resistance, and use of various carbon sources; some of these candidates are in genes previously known to affect these traits whereas others were in genes that have not been well characterized. Our results point to the potential power of association analyses in bacteria, but also to the need to carefully evaluate the potential for false associations.

Effect of Acyl Activating Enzyme (AAE) 3 on the growth and development of Medicago truncatula

The Acyl-Activating Enzyme (AAE) 3 gene encodes an oxalyl-CoA synthetase that catalyzes the conversion of oxalate to oxalyl-CoA in a CoA and ATP-dependent manner. Although the biochemical activity of AAE3 has been established, its biological role in plant growth and development remains unclear. To advance our understanding of the role of AAE3 in plant growth and development, we report here the characterization of two Medicago truncatula AAE3 (Mtaae3) mutants. Characterization of a Mtaae3 RNAi mutant revealed an accumulation of calcium oxalate crystals and increased seed permeability. These phenotypes were also exhibited in the Arabidopsis aae3 (Ataae3) mutants. Unlike the Ataae3 mutants, the Mtaae3 RNAi mutant did not show a reduction in vegetative growth, decreased seed germination, or increased seed calcium concentration. In an effort to clarify these phenotypic differences, a Mtaae3 Tnt1 mutant was identified and characterized. This Mtaae3 Tnt1 mutant displayed reduced vegetative growth, decreased seed germination, and increased seed calcium concentration as well as an accumulation of calcium oxalate crystals and increased seed permeability as found in Ataae3. Overall, the results presented here show the importance of AAE3 in the growth and development of plants. In addition, this study highlights the ability to separate specific growth and development phenotypes based on the level of AAE3 gene expression.

Pollen limitation and resource limitation affect the reproductive success of Medicago sativa L

Background

A large proportion of the flowers and ovules of plants do not develop into fruits and seeds. Plant reproduction may be limited because of pollen limitation and resource limitation. Medicago sativa L. is an ecologically important species in northwest China. We conducted a pollen supplementation experiment to determine the degree of pollen limitation in this species and detect the possible effects of resource allocation on pollen supplementation. We crossed two factors, pollen level (natural condition and hand pollinated) and resource level (control, water added, and fertilizer added), to estimate the effects of pollen addition and resource limitation on the opening of flowers and seed set. We also analyzed the floral characters, visitation frequency of pollinators and pollinator activity to estimate the effect of pollinators on the reproduction of M. sativa.

Results

Our results indicated that addition of pollen to some flowers did not divert resources from other flowers and that the addition of pollen boosted the seed set per flower, with no effect on flower number. The primary effect of resource limitation was on the number of flowers produced; however, there was no significant effect on seed set per flower. These findings showed that pollen limitation was an important limiting factor for seed set. In addition, Andrena lebedevi Popov was identified as the most effective pollinator, and pollinator visiting and activity affected reproduction success in M. sativa.

Conclusions

We found outcrossing was dominant in the breeding system and insect pollination played an important role in outcrossing. These findings have identified the dominant factor influencing seed set of M. sativa. This study aspires to contribute to a better understanding of pollen limitation, resource limitation and reproductive success.

Functional Genomics and Flowering Time in Medicago truncatula: An Overview

Flowering time is an important trait that influences adaptation and yield in many crop legumes. Both the inherent earliness of flowering and the degree to which it is responsive to environmental factors determine both the eco-geographic range across which crops can be successfully grown and the seasonal cycles most suitable for production. This chapter will provide a brief review of studies investigating the genetic control of flowering time in Medicago truncatula.

CRISPR/Cas9 and TALENs generate heritable mutations for genes involved in small RNA processing of Glycine max and Medicago truncatula

Processing of double-stranded RNA precursors into small RNAs is an essential regulator of gene expression in plant development and stress response. Small RNA processing requires the combined activity of a functionally diverse group of molecular components. However, in most of the plant species, there are insufficient mutant resources to functionally characterize each encoding gene. Here, mutations in loci encoding protein machinery involved in small RNA processing in soya bean and Medicago truncatula were generated using the CRISPR/Cas9 and TAL-effector nuclease (TALEN) mutagenesis platforms. An efficient CRISPR/Cas9 reagent was used to create a bi-allelic double mutant for the two soya bean paralogous Double-stranded RNA-binding2 (GmDrb2a and GmDrb2b) genes. These mutations, along with a CRISPR/Cas9-generated mutation of the M. truncatula Hua enhancer1 (MtHen1) gene, were determined to be germ-line transmissible. Furthermore, TALENs were used to generate a mutation within the soya bean Dicer-like2 gene. CRISPR/Cas9 mutagenesis of the soya bean Dicer-like3 gene and the GmHen1a gene was observed in the T0 generation, but these mutations failed to transmit to the T1 generation. The irregular transmission of induced mutations and the corresponding transgenes was investigated by whole-genome sequencing to reveal a spectrum of non-germ-line-targeted mutations and multiple transgene insertion events. Finally, a suite of combinatorial mutant plants were generated by combining the previously reported Gmdcl1a, Gmdcl1b and Gmdcl4b mutants with the Gmdrb2ab double mutant. Altogether, this study demonstrates the synergistic use of different genome engineering platforms to generate a collection of useful mutant plant lines for future study of small RNA processing in legume crops.