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

Ex vivo lipidomics reveal monoacylglycerols as substrates for a fatty acid amide hydrolase in the legume Medicago truncatula

Fatty acid amide hydrolase (FAAH) is a conserved hydrolase in eukaryotes with promiscuous activity toward a range of acylamide substrates. The native substrate repertoire for FAAH has just begun to be explored in plant systems outside the model Arabidopsis thaliana. Here, we used ex vivo lipidomics to identify potential endogenous substrates for Medicago truncatula FAAH1 (MtFAAH1). We incubated recombinant MtFAAH1 with lipid mixtures extracted from M. truncatula and resolved their profiles via gas chromatography-mass spectrometry (GC-MS). Data revealed that besides N-acylethanolamines (NAEs), sn-1 or sn-2 isomers of monoacylglycerols (MAGs) were substrates for MtFAAH1. Combined with in vitro and computational approaches, our data support both amidase and esterase activities for MtFAAH1. MAG-mediated hydrolysis via MtFAAH1 may be linked to biological roles that are yet to be discovered.

The phosphate starvation response regulator PHR2 antagonizes arbuscule maintenance in Medicago

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.

Methylated chalcones are required for rhizobial nod gene induction in the Medicago truncatula rhizosphere

Legume nodulation requires the detection of flavonoids in the rhizosphere by rhizobia to activate their production of Nod factor countersignals. Here we investigated the flavonoids involved in nodulation of Medicago truncatula. We biochemically characterized five flavonoid-O-methyltransferases (OMTs) and a lux-based nod gene reporter was used to investigate the response of Sinorhizobium medicae NodD1 to various flavonoids. We found that chalcone-OMT 1 (ChOMT1) and ChOMT3, but not OMT2, 4, and 5, were able to produce 4,4'-dihydroxy-2'-methoxychalcone (DHMC). The bioreporter responded most strongly to DHMC, while isoflavones important for nodulation of soybean (Glycine max) showed no activity. Mutant analysis revealed that loss of ChOMT1 strongly reduced DHMC levels. Furthermore, chomt1 and omt2 showed strongly reduced bioreporter luminescence in their rhizospheres. In addition, loss of both ChOMT1 and ChOMT3 reduced nodulation, and this phenotype was strengthened by the further loss of OMT2. We conclude that: the loss of ChOMT1 greatly reduces root DHMC levels; ChOMT1 or OMT2 are important for nod gene activation in the rhizosphere; and ChOMT1/3 and OMT2 promote nodulation. Our findings suggest a degree of exclusivity in the flavonoids used for nodulation in M. truncatula compared to soybean, supporting a role for flavonoids in rhizobial host range.

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.

Light-sensitive short hypocotyl genes confer symbiotic nodule identity in the legume Medicago truncatula

Legumes produce specialized root nodules that are distinct from lateral roots in morphology and function, with nodules intracellularly hosting nitrogen-fixing bacteria. We have previously shown that a lateral root program underpins nodule initiation, but there must be additional developmental regulators that confer nodule identity. Here, we show two members of the LIGHT-SENSITIVE SHORT HYPOCOTYL (LSH) transcription factor family, predominantly known to define shoot meristem complexity and organ boundaries, function as regulators of nodule organ identity. In parallel to the root initiation program, LSH1/LSH2 recruit a program into the root cortex that mediates the divergence into nodules, in particular with cell divisions in the mid-cortex. This includes regulation of auxin and cytokinin, promotion of NODULE ROOT1/2 and Nuclear Factor YA1, and suppression of the lateral root program. A principal outcome of LSH1/LSH2 function is the production of cells able to accommodate nitrogen-fixing bacteria, a key feature unique to nodules.

Biochemical analysis of the TPS-a subfamily in Medicago truncatula

Terpenes are important mediators of plant chemical response to environmental cues. Here, we describe the genome-wide identification and biochemical characterization of TPS-a members in Medicago truncatula, a model legume crop. Genome mining identified thirty-nine full-length terpene synthases with a significant number predicted to produce monoterpenes and sesquiterpenes. Biochemical characterization of the TPS-a subfamily associated with sesquiterpene biosynthesis revealed such compounds, that exhibit substantial biological activity in other plants. Gene expression analysis using qPCR and the Medicago gene atlas illustrated distinct tissue and time-based variation in expression in leaves and roots. Together our work establishes the gene-to-metabolite relationships for sesquiterpene synthases in M. truncatula. Understanding the biosynthetic capacity is a foundational step to defining the ecological roles of this important family of compounds.

CRISPR/Cas9-Mediated Generation of Mutant Lines in Medicago truncatula Indicates a Symbiotic Role of MtLYK10 during Nodule Formation

CRISPR/Cas9 systems are commonly used for plant genome editing; however, the generation of homozygous mutant lines in Medicago truncatula remains challenging. Here, we present a CRISPR/Cas9-based protocol that allows the efficient generation of M. truncatula mutants. Gene editing was performed for the LysM receptor kinase gene MtLYK10 and two major facilitator superfamily transporter genes. The functionality of CRISPR/Cas9 vectors was tested in Nicotiana benthamiana leaves by editing a co-transformed GUSPlus gene. Transformed M. truncatula leaf explants were regenerated to whole plants at high efficiency (80%). An editing efficiency (frequency of mutations at a given target site) of up to 70% was reached in the regenerated plants. Plants with MtLYK10 knockout mutations were propagated, and three independent homozygous mutant lines were further characterized. No off-target mutations were identified in these lyk10 mutants. Finally, the lyk10 mutants and wild-type plants were compared with respect to the formation of root nodules induced by nitrogen-fixing Sinorhizobium meliloti bacteria. Nodule formation was considerably delayed in the three lyk10 mutant lines. Surprisingly, the size of the rare nodules in mutant plants was higher than in wild-type plants. In conclusion, the symbiotic characterization of lyk10 mutants generated with the developed CRISPR/Cas9 protocol indicated a role of MtLYK10 in nodule formation.

Expression and mutagenesis studies in the Medicago truncatula iron transporter MtVTL8 confirm its role in symbiotic nitrogen fixation and reveal amino acids essential for transport

The model legume Medicago truncatula establishes a symbiosis with soil bacteria (rhizobia) that carry out symbiotic nitrogen fixation (SNF) in plant root nodules. SNF requires the exchange of nutrients between the plant and rhizobia in the nodule that occurs across a plant-derived symbiosome membrane. One iron transporter, belonging to the Vacuolar iron Transporter-Like (VTL) family, MtVTL8, has been identified as essential for bacteria survival and therefore SNF. In this work we investigated the spatial expression of MtVTL8 in nodules and addressed whether it could be functionally interchangeable with a similar nodule-expressed iron transporter, MtVTL4. Using a structural model for MtVTL8 and the previously hypothesized mechanism for iron transport in a phylogenetically-related Vacuolar Iron Transporter (VIT), EgVIT1 with known crystal structure, we identified critical amino acids and obtained their mutants. Mutants were tested in planta for complementation of an SNF defective line and in an iron sensitive mutant yeast strain. An extended phylogenetic assessment of VTLs and VITs showed that amino acids critical for function are conserved differently in VTLs vs. VITs. Our studies showed that some amino acids are essential for iron transport leading us to suggest a model for MtVTL8 function, one that is different for other iron transporters (VITs) studied so far. This study extends the understanding of iron transport mechanisms in VTLs as well as those used in SNF.

MODMS: a multi-omics database for facilitating biological studies on alfalfa (Medicago sativa L.)

Alfalfa (Medicago sativa L.) is a globally important forage crop. It also serves as a vegetable and medicinal herb because of its excellent nutritional quality and significant economic value. Multi-omics data on alfalfa continue to accumulate owing to recent advances in high-throughput techniques, and integrating this information holds great potential for expediting genetic research and facilitating advances in alfalfa agronomic traits. Therefore, we developed a comprehensive database named MODMS (multi-omics database of M. sativa) that incorporates multiple reference genomes, annotations, comparative genomics, transcriptomes, high-quality genomic variants, proteomics, and metabolomics. This report describes our continuously evolving database, which provides researchers with several convenient tools and extensive omics data resources, facilitating the expansion of alfalfa research. Further details regarding the MODMS database are available at https://modms.lzu.edu.cn/.

MtWOX2 and MtWOX9-1 Effects on the Embryogenic Callus Transcriptome in Medicago truncatula

WOX family transcription factors are well-known regulators of plant development, controlling cell proliferation and differentiation in diverse organs and tissues. Several WOX genes have been shown to participate in regeneration processes which take place in plant cell cultures in vitro, but the effects of most of them on tissue culture development have not been discovered yet. In this study, we evaluated the effects of MtWOX2 gene overexpression on the embryogenic callus development and transcriptomic state in Medicago truncatula. According to our results, overexpression of MtWOX2 leads to an increase in callus weight. Furthermore, transcriptomic changes in MtWOX2 overexpressing calli are, to a large extent, opposite to the changes caused by overexpression of MtWOX9-1, a somatic embryogenesis stimulator. These results add new information about the mechanisms of interaction between different WOX genes and can be useful for the search of new regeneration regulators.

Genome-wide identification and characterization of filamentation temperature-sensitive H (FtsH) genes and expression analysis in response to multiple stresses in Medicago truncatula

BACKGROUND

Filamentation temperature-sensitive H (FtsH) is an AAA+ ATP-dependent protease that plays a vital role in plant environmental adaption and tolerance. However, little is known about the function of the FtsH gene family in the most important legume model plant, Medicago truncatula.

METHODS AND RESULTS

To identify and investigate the potential stress adaptation roles of FtsH gene family in M. truncatula, we conducted a series of genome-wide characterization and expression analyses. Totally, twenty MtFtsH genes were identified, which were unevenly distributed across eight chromosomes and classified into six evolution groups based on their phylogenetic relationships, with each group containing similar structures and motifs. Furthermore, MtFtsH genes exhibited a high degree of collinearity and homology with leguminous plants such as alfalfa and soybean. Multiple cis-elements in the upstream region of MtFtsH genes were also identified that responded to light, abiotic stress, and phytohormones. Public RNA-seq data indicated that MtFtsH genes were induced under both salt and drought stresses, and our transcript expression analysis showed that MtFtsH genes of MtFtsH1, MtFtsH2, MtFtsH4, MtFtsH9, and MtFtsH10 were up-regulated after ABA, H2O2, PEG, and NaCl treatments. These results suggest that MtFtsH genes may play a critical role in drought and high salt stress responses and the adaption processes of plants.

CONCLUSIONS

This study provides a systematic analysis of FtsH gene family in M. truncatula, serving as a valuable molecular theoretical basis for future functional investigations. Our findings also extend the pool of potential candidate genes for the genetic improvement of abiotic stress tolerance in legume crops.

Comparative RNA sequencing-based transcriptome profiling of ten grapevine rootstocks: shared and specific sets of genes respond to mycorrhizal symbiosis

Arbuscular mycorrhizal symbiosis improves water and nutrient uptake by plants and provides them other ecosystem services. Grapevine is one of the major crops in the world. Vitis vinifera scions generally are grafted onto a variety of rootstocks that confer different levels of resistance against different pests, tolerance to environmental stress, and influence the physiology of the scions. Arbuscular mycorrhizal fungi are involved in the root architecture and in the immune response to soil-borne pathogens. However, the fine-tuned regulation and the transcriptomic plasticity of rootstocks in response to mycorrhization are still unknown. We compared the responses of 10 different grapevine rootstocks to arbuscular mycorrhizal symbiosis (AMS) formed with Rhizophagus irregularis DAOM197198 using RNA sequencing-based transcriptome profiling. We have highlighted a few shared regulation mechanisms, but also specific rootstock responses to R. irregularis colonization. A set of 353 genes was regulated by AMS in all ten rootstocks. We also compared the expression level of this set of genes to more than 2000 transcriptome profiles from various grapevine varieties and tissues to identify a class of transcripts related to mycorrhizal associations in these 10 rootstocks. Then, we compared the response of the 351 genes upregulated by mycorrhiza in grapevine to their Medicago truncatula homologs in response to mycorrhizal colonization based on available transcriptomic studies. More than 97% of the 351 M. truncatula-homologous grapevine genes were expressed in at least one mycorrhizal transcriptomic study, and 64% in every single RNAseq dataset. At the intra-specific level, we described, for the first time, shared and specific grapevine rootstock genes in response to R. irregularis symbiosis. At the inter-specific level, we defined a shared subset of mycorrhiza-responsive genes.

A snapshot of the transcriptome of Medicago truncatula (Fabales: Fabaceae) shoots and roots in response to an arbuscular mycorrhizal fungus and the pea aphid (Acyrthosiphon pisum) (Hemiptera: Aphididae)

Plants simultaneously interact with belowground symbionts such as arbuscular mycorrhizal (AM) fungi and aboveground antagonists such as aphids. Generally, plants gain access to valuable resources including nutrients and water through the AM symbiosis and are more resistant to pests. Nevertheless, aphids' performance improves on mycorrhizal plants, and it remains unclear whether a more nutritious food source and/or attenuated defenses are the contributing factors. This study examined the shoot and root transcriptome of barrel medic (Medicago truncatula Gaertn.) plants highly colonized by the AM fungus Rhizophagus irregularis (Blaszk., Wubet, Renker, and Buscot) C. Walker and A. Schüßler (Glomerales: Glomeraceae) and exposed to 7 days of mixed age pea aphid (Acyrthosiphon pisum (Harris)) herbivory. The RNA-seq samples chosen for this study showed that aphids were heavier when fed mycorrhizal plants compared to nonmycorrhizal plants. We hypothesized that (i) insect-related plant defense pathways will be downregulated in shoots of mycorrhizal plants with aphids compared to nonmycorrhizal plants with aphids; (ii) pathways involved in nutrient acquisition, carbohydrate-related and amino acid transport will be upregulated in shoots of mycorrhizal plants with aphids compared to nonmycorrhizal plants with aphids; and (iii) roots of mycorrhizal plants with aphids will exhibit mycorrhiza-induced resistance. The transcriptome data revealed that the gene repertoire related to defenses, nutrient transport, and carbohydrates differs between nonmycorrhizal and mycorrhizal plants with aphids, which could explain the weight gain in aphids. We also identified novel candidate genes that are differentially expressed in nonmycorrhizal plants with aphids, thus setting the stage for future functional studies.

A Novel Role of Medicago truncatula KNAT3/4/5-like Class 2 KNOX Transcription Factors in Drought Stress Tolerance

Class 2 KNOX homeobox transcription factors (KNOX2) play a role in promoting cell differentiation in several plant developmental processes. In Arabidopsis, they antagonize the meristematic KNOX1 function during leaf development through the modulation of phytohormones. In Medicago truncatula, three KNOX2 genes belonging to the KNAT3/4/5-like subclass (Mt KNAT3/4/5-like or MtKNOX3-like) redundantly works upstream of a cytokinin-signaling module to control the symbiotic root nodule formation. Their possible role in the response to abiotic stress is as-of-yet unknown. We produced transgenic M. truncatula lines, in which the expression of four MtKNOX3-like genes was knocked down by RNA interference. When tested for response to water withdrawal in the soil, RNAi lines displayed a lower tolerance to drought conditions compared to the control lines, measured as increased leaf water loss, accelerated leaf wilting time, and faster chlorophyll loss. Reanalysis of a transcriptomic M. truncatula drought stress experiment via cluster analysis and gene co-expression networks pointed to a possible role of MtKNOX3-like transcription factors in repressing a proline dehydrogenase gene (MtPDH), specifically at 4 days after water withdrawal. Proline measurement and gene expression analysis of transgenic RNAi plants compared to the controls confirmed the role of KNOX3-like genes in inhibiting proline degradation through the regulation of the MtPDH gene.

Development of a knowledge graph framework to ease and empower translational approaches in plant research: a use-case on grain legumes

While the continuing decline in genotyping and sequencing costs has largely benefited plant research, some key species for meeting the challenges of agriculture remain mostly understudied. As a result, heterogeneous datasets for different traits are available for a significant number of these species. As gene structures and functions are to some extent conserved through evolution, comparative genomics can be used to transfer available knowledge from one species to another. However, such a translational research approach is complex due to the multiplicity of data sources and the non-harmonized description of the data. Here, we provide two pipelines, referred to as structural and functional pipelines, to create a framework for a NoSQL graph-database (Neo4j) to integrate and query heterogeneous data from multiple species. We call this framework Orthology-driven knowledge base framework for translational research (Ortho_KB). The structural pipeline builds bridges across species based on orthology. The functional pipeline integrates biological information, including QTL, and RNA-sequencing datasets, and uses the backbone from the structural pipeline to connect orthologs in the database. Queries can be written using the Neo4j Cypher language and can, for instance, lead to identify genes controlling a common trait across species. To explore the possibilities offered by such a framework, we populated Ortho_KB to obtain OrthoLegKB, an instance dedicated to legumes. The proposed model was evaluated by studying the conservation of a flowering-promoting gene. Through a series of queries, we have demonstrated that our knowledge graph base provides an intuitive and powerful platform to support research and development programmes.

The saponin bomb: a nucleolar-localized β-glucosidase hydrolyzes triterpene saponins in Medicago truncatula

Plants often protect themselves from their own bioactive defense metabolites by storing them in less active forms. Consequently, plants also need systems allowing correct spatiotemporal reactivation of such metabolites, for instance under pathogen or herbivore attack. Via co-expression analysis with public transcriptomes, we determined that the model legume Medicago truncatula has evolved a two-component system composed of a β-glucosidase, denominated G1, and triterpene saponins, which are physically separated from each other in intact cells. G1 expression is root-specific, stress-inducible, and coregulated with that of the genes encoding the triterpene saponin biosynthetic enzymes. However, the G1 protein is stored in the nucleolus and is released and united with its typically vacuolar-stored substrates only upon tissue damage, partly mediated by the surfactant action of the saponins themselves. Subsequently, enzymatic removal of carbohydrate groups from the saponins creates a pool of metabolites with an increased broad-spectrum antimicrobial activity. The evolution of this defense system benefited from both the intrinsic condensation abilities of the enzyme and the bioactivity properties of its substrates. We dub this two-component system the saponin bomb, in analogy with the mustard oil and cyanide bombs, commonly used to describe the renowned β-glucosidase-dependent defense systems for glucosinolates and cyanogenic glucosides.

Two legume fatty acid amide hydrolase isoforms with distinct preferences for microbial- and plant-derived acylamides

Fatty acid amide hydrolase (FAAH) is a widely conserved amidase in eukaryotes, perhaps best known for inactivating N-acylethanolamine lipid mediators. However, FAAH enzymes hydrolyze a wide range of acylamide substrates. Analysis of FAAHs from multiple angiosperm species revealed two conserved phylogenetic groups that differed in key conserved residues in the substrate binding pocket. While the foundation group of plant FAAHs, designated FAAH1, has been studied at the structural and functional level in Arabidopsis thaliana, nothing is known about FAAH2 members. Here, we combined computational and biochemical approaches to compare the structural and enzymatic properties of two FAAH isoforms in the legume Medicago truncatula designated MtFAAH1 and MtFAAH2a. Differences in structural and physicochemical properties of the substrate binding pockets, predicted from homology modeling, molecular docking, and molecular dynamic simulation experiments, suggested that these two FAAH isoforms would exhibit differences in their amidohydrolase activity profiles. Indeed, kinetic studies of purified, recombinant MtFAAHs indicated a reciprocal preference for acylamide substrates with MtFAAH1 more efficiently utilizing long-chain acylamides, and MtFAAH2a more efficiently hydrolyzing short-chain and aromatic acylamides. This first report of the enzymatic behavior of two phylogenetically distinct plant FAAHs will provide a foundation for further investigations regarding FAAH isoforms in legumes and other plant species.

Transcript profiling of hairy vetch (Vicia villosa Roth) identified interesting genes for seed dormancy

Hairy vetch, a diploid annual legume species, has a robust growth habit, high biomass yield, and winter hardy characteristics. Seed hardness is a major constraint for growing hairy vetch commercially. Hard seeded cultivars are valuable as forages, whereas soft seeded and shatter resistant cultivars have advantages for their use as a cover crop. Transcript analysis of hairy vetch was performed to understand the genetic mechanisms associated with important hairy vetch traits. RNA was extracted from leaves, flowers, immature pods, seed coats, and cotyledons of contrasting soft and hard seeded "AU Merit" plants. A range of 31.22-79.18 Gb RNA sequence data per tissue sample were generated with estimated coverage of 1040-2639×. RNA sequence assembly and mapping of the contigs against the Medicago truncatula (V4.0) genome identified 76,422 gene transcripts. A total of 24,254 transcripts were constitutively expressed in hairy vetch tissues. Key genes, such as KNOX4 (a class II KNOTTED-like homeobox KNOXII gene), qHs1 (endo-1,4-β-glucanase), GmHs1-1 (calcineurin-like metallophosphoesterase), chitinase, shatterproof 1 and 2 (SHP1, SHP2), shatter resistant 1-5 (SHAT1-5)(NAC transcription factor), PDH1 (prephenate dehydrogenase 1), and pectin methylesterases with a potential role in seed hardness and pod shattering, were further explored based on genes involved in seed hardness from other species to query the hairy vetch transcriptome data. Identification of interesting candidate genes in hairy vetch can facilitate the development of improved cultivars with desirable seed characteristics for use as a forage and as a cover crop.

Sucrose-induced auxin conjugate hydrolase restores symbiosis in a Medicago cytokinin perception mutant

Rhizobia-legume interactions recruit cytokinin for the induction of nodule primordia in the cortex. Cytokinin signaling regulates auxin transport and biosynthesis, causing local auxin accumulation, which triggers cortical cell division. Since sugar signaling can trigger auxin responses, we explored whether sugar treatments could rescue symbiosis in the Medicago truncatula cytokinin response 1 (cre1) mutant. Herein, we demonstrate that sucrose and its nonmetabolizable isomer turanose can trigger auxin response and recover functional symbiosis in cre1, indicating sucrose signaling to be necessary for the restoration of symbiosis. In both M. truncatula A17 (wild type) and cre1, sucrose signaling significantly upregulated IAA-Ala Resistant 3 (IAR33), encoding an auxin conjugate hydrolase, in rhizobia-infected as well as in uninfected roots. Knockdown of IAR33 (IAR33-KD) significantly reduced nodulation in A17, highlighting the importance of deconjugation-mediated auxin accumulation during nodule inception. In cre1, IAR33-KD restricted the sucrose-mediated restoration of functional symbiosis, suggesting that deconjugation-mediated auxin accumulation plays a key role in the absence of CRE1-mediated auxin biosynthesis and transport control. Overexpression of IAR33 also restored functional symbiosis in cre1, further suggesting that IAR33 mediates auxin accumulation in response to sucrose signaling. Since all the observed sucrose-mediated responses were common to A17 and cre1, deconjugation-mediated auxin response appeared to be independent of CRE1, which normally governs local auxin accumulation in the presence of rhizobia. We propose that sucrose-dependent restoration of symbiosis in cre1 occurs by the activation of IAR33-mediated auxin deconjugation.

Long-lasting impact of chitooligosaccharide application on strigolactone biosynthesis and fungal accommodation promotes arbuscular mycorrhiza in Medicago truncatula

The establishment of arbuscular mycorrhiza (AM) between plants and Glomeromycotina fungi is preceded by the exchange of chemical signals: fungal released Myc-factors, including chitooligosaccharides (CO) and lipo-chitooligosaccharides (LCO), activate plant symbiotic responses, while root-exuded strigolactones stimulate hyphal branching and boost CO release. Furthermore, fungal signaling reinforcement through CO application was shown to promote AM development in Medicago truncatula, but the cellular and molecular bases of this effect remained unclear. Here, we focused on long-term M. truncatula responses to CO treatment, demonstrating its impact on the transcriptome of both mycorrhizal and nonmycorrhizal roots over several weeks and providing an insight into the mechanistic bases of the CO-dependent promotion of AM colonization. CO treatment caused the long-lasting regulation of strigolactone biosynthesis and fungal accommodation-related genes. This was mirrored by an increase in root didehydro-orobanchol content, and the promotion of accommodation responses to AM fungi in root epidermal cells. Lastly, an advanced downregulation of AM symbiosis marker genes was observed at the latest time point in CO-treated plants, in line with an increased number of senescent arbuscules. Overall, CO treatment triggered molecular, metabolic, and cellular responses underpinning a protracted acceleration of AM development.

Symbiotic Nodule Development and Efficiency in the Medicago truncatula Mtefd-1 Mutant Is Highly Dependent on Sinorhizobium Strains

Symbiotic nitrogen fixation (SNF) can play a key role in agroecosystems to reduce the negative impact of nitrogen fertilizers. Its efficiency is strongly affected by the combination of bacterial and plant genotypes, but the mechanisms responsible for the differences in the efficiency of rhizobium strains are not well documented. In Medicago truncatula, SNF has been mostly studied using model systems, such as M. truncatula A17 in interaction with Sinorhizobium meliloti Sm2011. Here we analyzed both the wild-type (wt) A17 and the Mtefd-1 mutant in interaction with five S. meliloti and two Sinorhizobium medicae strains. ETHYLENE RESPONSE FACTOR REQUIRED FOR NODULE DIFFERENTIATION (MtEFD) encodes a transcription factor, which contributes to the control of nodule number and differentiation in M. truncatula. We found that, in contrast to Sm2011, four strains induce functional (Fix+) nodules in Mtefd-1, although less efficient for SNF than in wt A17. In contrast, the Mtefd-1 hypernodulation phenotype is not strain-dependent. We compared the plant nodule transcriptomes in response to SmBL225C, a highly efficient strain with A17, versus Sm2011, in wt and Mtefd-1 backgrounds. This revealed faster nodule development with SmBL225C and early nodule senescence with Sm2011. These RNA sequencing analyses allowed us to identify candidate plant factors that could drive the differential nodule phenotype. In conclusion, this work shows the value of having a set of rhizobium strains to fully evaluate the biological importance of a plant symbiotic gene.

Genome-Wide Identification, Characterization, and Expression Analysis of SPIRAL1 Family Genes in Legume Species

The SPIRAL1 (SPR1) gene family encodes microtubule-associated proteins that are essential for the anisotropic growth of plant cells and abiotic stress resistance. Currently, little is known about the characteristics and roles of the gene family outside of Arabidopsis thaliana. This study intended to investigate the SPR1 gene family in legumes. In contrast to that of A. thaliana, the gene family has undergone shrinking in the model legume species Medicago truncatula and Glycine max. While the orthologues of SPR1 were lost, very few SPR1-Like (SP1L) genes were identified given the genome size of the two species. Specifically, the M. truncatula and G. max genomes only harbor two MtSP1L and eight GmSP1L genes, respectively. Multiple sequence alignment showed that all these members contain conserved N- and C-terminal regions. Phylogenetic analysis clustered the legume SP1L proteins into three clades. The SP1L genes showed similar exon-intron organizations and similar architectures in their conserved motifs. Many essential cis-elements are present in the promoter regions of the MtSP1L and GmSP1L genes associated with growth and development, plant hormones, light, and stress. The expression analysis revealed that clade 1 and clade 2 SP1L genes have relatively high expression in all tested tissues in Medicago and soybean, suggesting their function in plant growth and development. MtSP1L-2, as well as clade 1 and clade 2 GmSP1L genes, display a light-dependent expression pattern. The SP1L genes in clade 2 (MtSP1L-2, GmSP1L-3, and GmSP1L-4) were significantly induced by sodium chloride treatment, suggesting a potential role in the salt-stress response. Our research provides essential information for the functional studies of SP1L genes in legume species in the future.

Widely conserved AHL transcription factors are essential for NCR gene expression and nodule development in Medicago

Symbiotic nitrogen fixation by Rhizobium bacteria in the cells of legume root nodules alleviates the need for nitrogen fertilizers. Nitrogen fixation requires the endosymbionts to differentiate into bacteroids which can be reversible or terminal. The latter is controlled by the plant, it is more beneficial and has evolved in multiple clades of the Leguminosae family. The plant effectors of terminal differentiation in inverted repeat-lacking clade legumes (IRLC) are nodule-specific cysteine-rich (NCR) peptides, which are absent in legumes such as soybean where there is no terminal differentiation of rhizobia. It was assumed that NCRs co-evolved with specific transcription factors, but our work demonstrates that expression of NCR genes does not require NCR-specific transcription factors. Introduction of the Medicago truncatula NCR169 gene under its own promoter into soybean roots resulted in its nodule-specific expression, leading to bacteroid changes associated with terminal differentiation. We identified two AT-Hook Motif Nuclear Localized (AHL) transcription factors from both M. truncatula and soybean nodules that bound to AT-rich sequences in the NCR169 promoter inducing its expression. Whereas mutation of NCR169 arrested bacteroid development at a late stage, the absence of MtAHL1 or MtAHL2 completely blocked bacteroid differentiation indicating that they also regulate other NCR genes required for the development of nitrogen-fixing nodules. Regulation of NCRs by orthologous transcription factors in non-IRLC legumes opens up the possibility of increasing the efficiency of nitrogen fixation in legumes lacking NCRs.

Stabilization of membrane topologies by proteinaceous remorin scaffolds

In plants, the topological organization of membranes has mainly been attributed to the cell wall and the cytoskeleton. Additionally, few proteins, such as plant-specific remorins have been shown to function as protein and lipid organizers. Root nodule symbiosis requires continuous membrane re-arrangements, with bacteria being finally released from infection threads into membrane-confined symbiosomes. We found that mutations in the symbiosis-specific SYMREM1 gene result in highly disorganized perimicrobial membranes. AlphaFold modelling and biochemical analyses reveal that SYMREM1 oligomerizes into antiparallel dimers and may form a higher-order membrane scaffolding structure. This was experimentally confirmed when expressing this and other remorins in wall-less protoplasts is sufficient where they significantly alter and stabilize de novo membrane topologies ranging from membrane blebs to long membrane tubes with a central actin filament. Reciprocally, mechanically induced membrane indentations were equally stabilized by SYMREM1. Taken together we describe a plant-specific mechanism that allows the stabilization of large-scale membrane conformations independent of the cell wall.

The putative transporter MtUMAMIT14 participates in nodule formation in Medicago truncatula

Transport systems are crucial in many plant processes, including plant-microbe interactions. Nodule formation and function in legumes involve the expression and regulation of multiple transport proteins, and many are still uncharacterized, particularly for nitrogen transport. Amino acids originating from the nitrogen-fixing process are an essential form of nitrogen for legumes. This work evaluates the role of MtN21 (henceforth MtUMAMIT14), a putative transport system from the MtN21/EamA-like/UMAMIT family, in nodule formation and nitrogen fixation in Medicago truncatula. To dissect this transporter's role, we assessed the expression of MtUMAMIT14 using GUS staining, localized the corresponding protein in M. truncatula root and tobacco leaf cells, and investigated two independent MtUMAMIT14 mutant lines. Our results indicate that MtUMAMIT14 is localized in endosomal structures and is expressed in both the infection zone and interzone of nodules. Comparison of mutant and wild-type M. truncatula indicates MtUMAMIT14, the expression of which is dependent on the presence of NIN, DNF1, and DNF2, plays a role in nodule formation and nitrogen-fixation. While the function of the transporter is still unclear, our results connect root nodule nitrogen fixation in legumes with the UMAMIT family.

Anthocyanin pigmentation as a quantitative visual marker for arbuscular mycorrhizal fungal colonization of Medicago truncatula roots

Visualization of root colonization by arbuscular mycorrhizal fungi (AMF) is the most elementary experiment in the field of mycorrhizal symbiosis. The most widely used approach for evaluating levels of AMF colonization is staining with trypan blue or ink, which is scored using the time-consuming grid intersection method. Here we demonstrate the use of an anthocyanin-based visual marker system for visualizing AMF colonization of Medicago truncatula roots. Expression of MtLAP1, a transcription factor which regulates the production of anthocyanins, from the AMF-induced Kunitz Protease Inhibitor 106 promoter, allowed the visualization of arbuscules in live plant tissues without microscopy or staining. This marker system allowed straightforward qualitative evaluation of the ram1, vpy and dmi3 AMF phenotypes using Agrobacterium rhizogenes hairy-root transformation. For the strigolactone biosynthesis mutant carotenoid cleavage dioxygenase 8a and a novel mutant scooby, which show quantitative AMF symbiotic phenotypes, the amount of anthocyanins in the roots estimated by spectrophotometry correlated very well with colonization levels estimated by staining and scoring using the grid intersection method. The LAP1-based marker system therefore provides a highly efficient approach for mutant screening and monitoring of AMF colonization in live tissues by eye, or for quantitative assessment using a simple and quick photometric assay.

A legume kinesin controls vacuole morphogenesis for rhizobia endosymbiosis

Symbioses between legumes and rhizobia require establishment of the plant-derived symbiosome membrane, which surrounds the rhizobia and accommodates the symbionts by providing an interface for nutrient and signal exchange. The host cytoskeleton and endomembrane trafficking systems play central roles in the formation of a functional symbiotic interface for rhizobia endosymbiosis; however, the underlying mechanisms remain largely unknown. Here we demonstrate that the nodulation-specific kinesin-like calmodulin-binding protein (nKCBP), a plant-specific microtubule-based kinesin motor, controls central vacuole morphogenesis in symbiotic cells in Medicago truncatula. Phylogenetic analysis further indicated that nKCBP duplication occurs solely in legumes of the clade that form symbiosomes. Knockout of nKCBP results in central vacuole deficiency, defective symbiosomes and abolished nitrogen fixation. nKCBP decorates linear particles along microtubules, and crosslinks microtubules with the actin cytoskeleton, to control central vacuole formation by modulating vacuolar vesicle fusion in symbiotic cells. Together, our findings reveal that rhizobia co-opted nKCBP to achieve symbiotic interface formation by regulating cytoskeletal assembly and central vacuole morphogenesis during nodule development.

Signaling by reactive molecules and antioxidants in legume nodules

Legume nodules are symbiotic structures formed as a result of the interaction with rhizobia. Nodules fix atmospheric nitrogen into ammonia that is assimilated by the plant and this process requires strict metabolic regulation and signaling. Reactive oxygen species (ROS) and reactive nitrogen species (RNS) are involved as signal molecules at all stages of symbiosis, from rhizobial infection to nodule senescence. Also, reactive sulfur species (RSS) are emerging as important signals for an efficient symbiosis. Homeostasis of reactive molecules is mainly accomplished by antioxidant enzymes and metabolites and is essential to allow redox signaling while preventing oxidative damage. Here, we examine the metabolic pathways of reactive molecules and antioxidants with an emphasis on their functions in signaling and protection of symbiosis. In addition to providing an update of recent findings while paying tribute to original studies, we identify several key questions. These include the need of new methodologies to detect and quantify ROS, RNS, and RSS, avoiding potential artifacts due to their short lifetimes and tissue manipulation; the regulation of redox-active proteins by post-translational modification; the production and exchange of reactive molecules in plastids, peroxisomes, nuclei, and bacteroids; and the unknown but expected crosstalk between ROS, RNS, and RSS in nodules.

MtNF-YC6 and MtNF-YC11 are involved in regulating the transcriptional program of arbuscular mycorrhizal symbiosis

Arbuscular mycorrhizal fungi are obligate symbionts that transfer mineral nutrients to host plants through arbuscules, a fungal structure specialized for exchange for photosynthetic products. MtNF-YC6 and MtNF-YC11, which encode the C subunits of nuclear factor Y (NF-Y) family in Medicago truncatula are induced specifically by arbuscular mycorrhizal symbiosis (AMS). A previous study showed that MtNF-YC6 and MtNF-YC11 are activated in cortical cells of mycorrhizal roots, but the gene functions were unknown. Herein, we identified both MtNF-YB17 and MtNF-YB12 as the interacting partners of MtNF-YC6 and MtNF-YC11 in yeast and plants. MtNF-YB17 was highly induced by AMS and activated in cortical cells only in mycorrhizal roots but MtNF-YB12 was not affected. The formation of B/C heterodimers led the protein complexes to transfer from the cytoplasm to the nucleus. Silencing MtNF-YC6 and C11 by RNA interference (RNAi) resulted in decreased colonization efficiency and arbuscule richness. Coincidently, genes associated with arbuscule development and degeneration in RNAi roots were also downregulated. In silico analysis showed CCAAT-binding motifs in the promoter regions of downregulated genes, further supporting the involvement of NF-Y complexes in transcriptional regulation of symbiosis. Taken together, this study identifies MtNF-YC6- or MtNF-YC11-containing protein complexes as novel transcriptional regulators of symbiotic program and provides a list of potential downstream target genes. These data will help to further dissect the AMS regulatory network.

DNA demethylation and hypermethylation are both required for late nodule development in Medicago

Plant epigenetic regulations are involved in transposable element silencing, developmental processes and responses to the environment^1-7. They often involve modifications of DNA methylation, particularly through the DEMETER (DME) demethylase family and RNA-dependent DNA methylation (RdDM)⁸. Root nodules host rhizobia that can fix atmospheric nitrogen for the plant's benefit in nitrogen-poor soils. The development of indeterminate nodules, as in Medicago truncatula, involves successive waves of gene activation^9-12, control of which raises interesting questions. Using laser capture microdissection (LCM) coupled to RNA-sequencing (SYMbiMICS data¹¹), we previously identified 4,309 genes (termed NDD) activated in the nodule differentiation and nitrogen fixation zones, 36% of which belong to co-regulated genomic regions dubbed symbiotic islands¹³. We found MtDME to be upregulated in the differentiation zone and required for nodule development, and we identified 474 differentially methylated regions hypomethylated in the nodule by analysing ~2% of the genome⁴. Here, we coupled LCM and whole-genome bisulfite sequencing for a comprehensive view of DNA methylation, integrated with gene expression at the tissue level. Furthermore, using CRISPR-Cas9 mutagenesis of MtDRM2, we showed the importance of RdDM for CHH hypermethylation and nodule development. We thus proposed a model of DNA methylation dynamics during nodule development.

The Hidden Potential of High-Throughput RNA-Seq Re-Analysis, a Case Study for DHDPS, Key Enzyme of the Aspartate-Derived Lysine Biosynthesis Pathway and Its Role in Abiotic and Biotic Stress Responses in Soybean

DHDPS is a key enzyme in the aspartate-derived lysine biosynthesis pathway and an evident object of study for biofortification strategies in plants. DHDPS isoforms with novel regulatory properties in Medicago truncatula were demonstrated earlier and hypothesized to be involved in abiotic and biotic stress responses. Here, we present a phylogenetic analysis of the DHPDS gene family in land plants which establishes the existence of a legume-specific class of DHDPS, termed DHDPS B-type, distinguishable from the DHDPS A-type commonly present in all land plants. The G. max genome comprises two A-type DHDPS genes (Gm.DHDPS-A1; Glyma.09G268200, Gm.DHDPS-A2; Glyma.18G221700) and one B-type (Gm.DHDPS-B; Glyma.03G022300). To further investigate the expression pattern of the G. max DHDPS isozymes in different plant tissues and under various stress conditions, 461 RNA-seq experiments were exploited and re-analyzed covering two expression atlases, 13 abiotic and 5 biotic stress studies. Gm.DHDPS-B is seen almost exclusively expressed in roots and nodules in addition to old cotyledons or senescent leaves while both DHDPS A-types are expressed constitutively in all tissues analyzed with the highest expression in mature seeds. Furthermore, Gm.DHDPS-B expression is significantly upregulated in some but not all stress responses including salt stress, flooding, ethylene or infection with Phytophthora sojae and coincides with downregulation of DHDPS A-types. In conclusion, we demonstrate the potential of an in-depth RNA-seq re-analysis for the guidance of future experiments and to expand on current knowledge.

Regulation of the Zinc Deficiency Response in the Legume Model Medicago truncatula

The zinc deficiency response in Arabidopsis thaliana is regulated by F-group basic region leucine-zipper (F-bZIP) transcription factors, and there is evidence of evolutionary conservation of this regulatory network in land plants. Fundamental knowledge on the zinc homeostasis regulation in crop species will contribute to improving their zinc nutritional value. Legumes are protein-rich crops, used worldwide as part of traditional diets and as animal forage, being therefore a good target for micronutrient biofortification. Here, we identified F-bZIP transcription factors in representative legume species and functionally characterized the two F-bZIPs from Medicago truncatula. Results indicate that MtFbZIP1 is the functional homolog of A. thaliana bZIP19 and bZIP23, while MtFbZIP2 does not play a role in the zinc deficiency response. Additionally, analysis of M. truncatula genes from the Zrt/Irt-like protein (ZIP) family of zinc transporters or encoding nicotianamine synthase enzymes that produce the zinc ligand nicotianamine, support the conservation of the F-bZIP-regulated zinc deficiency response in M. truncatula. Phylogenetic analysis of F-bZIP homologs enriched in legume species reinforces the branching into two groups, with MtFbZIP1 and MtFbZIP2 mapping in Groups 1 and 2, respectively. This phylogeny combined with the functional characterization of MtFbZIPs supports the suggested conservation of the zinc deficiency response associated with Group 1 F-bZIPs, and the more variable evolutionary paths associated with Group 2. Overall, we provide novel insight on the mechanisms of response to zinc deficiency in M. truncatula, which contributes to developing strategies for improving zinc content in legume crops.

MtEFD and MtEFD2: Two transcription factors with distinct neofunctionalization in symbiotic nodule development

Rhizobium-legume nitrogen-fixing symbiosis involves the formation of a specific organ, the root nodule, which provides bacteria with the proper cellular environment for atmospheric nitrogen fixation. Coordinated differentiation of plant and bacterial cells is an essential step of nodule development, for which few transcriptional regulators have been characterized. Medicago truncatula ETHYLENE RESPONSE FACTOR REQUIRED FOR NODULE DIFFERENTIATION (MtEFD) encodes an APETALA2/ETHYLENE RESPONSIVE FACTOR (ERF) transcription factor, the mutation of which leads to both hypernodulation and severe defects in nodule development. MtEFD positively controls a negative regulator of cytokinin signaling, the RESPONSE REGULATOR 4 (MtRR4) gene. Here we showed that that the Mtefd-1 mutation affects both plant and bacterial endoreduplication in nodules, as well as the expression of hundreds of genes in young and mature nodules, upstream of known regulators of symbiotic differentiation. MtRR4 expressed with the MtEFD promoter complemented Mtefd-1 hypernodulation but not the nodule differentiation phenotype. Unexpectedly, a nonlegume homolog of MtEFD, AtERF003 in Arabidopsis (Arabidopsis thaliana), could efficiently complement both phenotypes of Mtefd-1, in contrast to the MtEFD paralog MtEFD2 expressed in the root and nodule meristematic zone. A domain swap experiment showed that MtEFD2 differs from MtEFD by its C-terminal fraction outside the DNA binding domain. Furthermore, clustered regularly interspaced short palindromic repeats-CRISPR associated protein 9 (CRISPR-Cas9) mutagenesis of MtEFD2 led to a reduction in the number of nodules formed in Mtefd-1, with downregulation of a set of genes, including notably NUCLEAR FACTOR-YA1 (MtNF-YA1) and MtNF-YB16, which are essential for nodule meristem establishment. We, therefore, conclude that nitrogen-fixing symbiosis recruited two proteins originally expressed in roots, MtEFD and MtEFD2, with distinct functions and neofunctionalization processes for each of them.

Alfalfa (Medicago sativa L.) pho2 mutant plants hyperaccumulate phosphate

In this article, we describe a set of novel alfalfa (Medicago sativa L.) plants that hyper-accumulate Phosphate ion (Pi) at levels 3- to 6-fold higher than wild-type. This alfalfa germplasm will have practical applications reclaiming Pi from contaminated or enriched soil or be used in conservation buffer strips to protect waterways from Pi run-off. Hyper-accumulating alfalfa plants were generated by targeted mutagenesis of PHOSPHATE2 (PHO2) using newly created CRISPR/Cas9 reagents and an improved mutant screening strategy. PHO2 encodes a ubiquitin conjugating E2 enzyme (UBC24) previously characterized in Arabidopsis thaliana, Medicago truncatula, and Oryza sativa. Mutations of PHO2 disrupt Pi homeostasis resulting in Pi hyper-accumulation. Successful CRISPR/Cas9 editing of PHO2 demonstrates that this is an efficient mutagenesis tool in alfalfa despite its complex autotetraploid genome structure. Arabidopsis and M. truncatula ortholog genes were used to identify PHO2 haplotypes in outcrossing tetraploid M. sativa with the aim of generating heritable mutations in both PHO2-like genes (PHO2-B and PHO2-C). After delivery of the reagent and regeneration from transformed leaf explants, plants with mutations in all haplotypes of PHO2-B and PHO2-C were identified. These plants were evaluated for morphology, Pi accumulation, heritable transmission of targeted mutations, segregation of mutant haplotypes and removal of T-DNA(s). The Agrobacterium-mediated transformation assay and gene editing reagents reported here were also evaluated for further optimization for future alfalfa functional genomic studies.

Functional Characterization of MtrGSTF7, a Glutathione S-Transferase Essential for Anthocyanin Accumulation in Medicago truncatula

Flavonoids are essential compounds widespread in plants and exert many functions such as defence, definition of organ colour and protection against stresses. In Medicago truncatula, flavonoid biosynthesis and accumulation is finely regulated in terms of tissue specificity and induction by external factors, such as cold and other stresses. Among flavonoids, anthocyanin precursors are synthesised in the cytoplasm, transported to the tonoplast, then imported into the vacuole for further modifications and storage. In the present work, we functionally characterised MtrGSTF7, a phi-class glutathione S-transferase involved in anthocyanin transport to the tonoplast. The mtrgstf7 mutant completely lost the ability to accumulate anthocyanins in leaves both under control and anthocyanin inductive conditions. On the contrary, this mutant showed an increase in the levels of soluble proanthocyanidins (Pas) in their seeds with respect to the wild type. By complementation and expression data analysis, we showed that, differently from A. thaliana and similarly to V. vinifera, transport of anthocyanin and proanthocyanidins is likely carried out by different GSTs belonging to the phi-class. Such functional diversification likely results from the plant need to finely tune the accumulation of diverse classes of flavonoids according to the target organs and developmental stages.

Multidimensional Gene Regulatory Landscape of Motor Organ Pulvinus in the Model Legume Medicago truncatula

Nyctinastic leaf movement of Fabaceae is driven by the tiny motor organ pulvinus located at the base of the leaf or leaflet. Despite the increased understanding of the essential role of ELONGATED PETIOLULE1 (ELP1)/PETIOLE LIKE PULVINUS (PLP) orthologs in determining pulvinus identity in legumes, key regulatory components and molecular mechanisms underlying this movement remain largely unclear. Here, we used WT pulvinus and the equivalent tissue in the elp1 mutant to carry out transcriptome and proteome experiments. The omics data indicated that there are multiple cell biological processes altered at the gene expression and protein abundance level during the pulvinus development. In addition, comparative analysis of different leaf tissues provided clues to illuminate the possible common primordium between pulvinus and petiole, as well as the function of ELP1. Furthermore, the auxin pathway, cell wall composition and chloroplast distribution were altered in elp1 mutants, verifying their important roles in pulvinus development. This study provides a comprehensive insight into the motor organ of the model legume Medicago truncatula and further supplies a rich dataset to facilitate the identification of novel players involved in nyctinastic movement.