The association between Dioscorea sansibarensis and Orrella dioscoreae as a model for hereditary leaf symbiosis.

Hereditary, or vertically-transmitted, symbioses affect a large number of animal species and some plants. The precise mechanisms underlying transmission of functions of these associations are often difficult to describe, due to the difficulty in separating the symbiotic partners. This is especially the case for plant-bacteria hereditary symbioses, which lack experimentally tractable model systems. Here, we demonstrate the potential of the leaf symbiosis between the wild yam Dioscorea sansibarensis and the bacterium Orrella dioscoreae (O. dioscoreae) as a model system for hereditary symbiosis. O. dioscoreae is easy to grow and genetically manipulate, which is unusual for hereditary symbionts. These properties allowed us to design an effective antimicrobial treatment to rid plants of bacteria and generate whole aposymbiotic plants, which can later be re-inoculated with bacterial cultures. Aposymbiotic plants did not differ morphologically from symbiotic plants and the leaf forerunner tip containing the symbiotic glands formed normally even in the absence of bacteria, but microscopic differences between symbiotic and aposymbiotic glands highlight the influence of bacteria on the development of trichomes and secretion of mucilage. This is to our knowledge the first leaf symbiosis where both host and symbiont can be grown separately and where the symbiont can be genetically altered and reintroduced to the host.

Structures and roles of BcsD and partner scaffold proteins in proteobacterial cellulose secretion.

Cellulose is the world's most abundant biopolymer, and similar to its role as a cell wall component in plants, it is a prevalent constituent of the extracellular matrix in bacterial biofilms. Although bacterial cellulose (BC) was first described in the 19th century, it was only recently revealed that it is produced by several distinct types of Bcs secretion systems that feature multiple accessory subunits in addition to a catalytic BcsAB synthase tandem. We recently showed that crystalline cellulose secretion in the Gluconacetobacter genus (α-Proteobacteria) is driven by a supramolecular BcsH-BcsD scaffold-the "cortical belt"-which stabilizes the synthase nanoarrays through an unexpected inside-out mechanism for secretion system assembly. Interestingly, while bcsH is specific for Gluconacetobacter, bcsD homologs are widespread in Proteobacteria. Here, we examine BcsD homologs and their gene neighborhoods from several plant-colonizing ß- and γ-Proteobacteria proposed to secrete a variety of non-crystalline and/or chemically modified cellulosic polymers. We provide structural and mechanistic evidence that through different quaternary structure assemblies BcsD acts with proline-rich BcsH, BcsP, or BcsO partners across the proteobacterial clade to form synthase-interacting intracellular scaffolds that, in turn, determine the biofilm strength and architecture in species with strikingly different physiology and secreted biopolymers.

Motility-Independent Vertical Transmission of Bacteria in Leaf Symbiosis.

Hereditary symbioses have the potential to drive transgenerational effects, yet the mechanisms responsible for transmission of heritable plant symbionts are still poorly understood. The leaf symbiosis between Dioscorea sansibarensis and the bacterium Orrella dioscoreae offers an appealing model system to study how heritable bacteria are transmitted to the next generation. Here, we demonstrate that inoculation of apical buds with a bacterial suspension is sufficient to colonize newly formed leaves and propagules, and to ensure transmission to the next plant generation. Flagellar motility is not required for movement inside the plant but is important for the colonization of new hosts. Further, tissue-specific regulation of putative symbiotic functions highlights the presence of two distinct subpopulations of bacteria in the leaf gland and at the shoot meristem. We propose that bacteria in the leaf gland dedicate resources to symbiotic functions, while dividing bacteria in the shoot tip ensure successful colonization of meristematic tissue, glands, and propagules. Compartmentalization of intrahost populations together with tissue-specific regulation may serve as a robust mechanism for the maintenance of mutualism in leaf symbiosis. IMPORTANCE Hereditary symbioses with bacteria are common in the animal kingdom, but relatively unexplored in plants. Several plant species form associations with bacteria in their leaves, which is called leaf symbiosis. These associations are highly specific, but the mechanisms responsible for symbiont transmission are poorly understood. Using the association between the yam species Dioscorea sansibarensis and Orrella dioscoreae as a model leaf symbiosis, we show that bacteria are distributed to specific leaf structures via association with shoot meristems. Flagellar motility is required for initial infection but does not contribute to spread within host tissue. We also provide evidence that bacterial subpopulations at the meristem or in the symbiotic leaf gland differentially express key symbiotic genes. We argue that this separation of functional symbiont populations, coupled with tight control over bacterial infection and transmission, explain the evolutionary robustness of leaf symbiosis. These findings may provide insights into how plants may recruit and maintain beneficial symbionts at the leaf surface.

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.

µLAS technology for DNA isolation coupled to Cas9-assisted targeting for sequencing and assembly of a 30 kb region in plant genome.

Cas9-assisted targeting of DNA fragments in complex genomes is viewed as an essential strategy to obtain high-quality and continuous sequence data. However, the purity of target loci selected by pulsed-field gel electrophoresis (PFGE) has so far been insufficient to assemble the sequence in one contig. Here, we describe the µLAS technology to capture and purify high molecular weight DNA. First, the technology is optimized to perform high sensitivity DNA profiling with a limit of detection of 20 fg/µl for 50 kb fragments and an analytical time of 50 min. Then, µLAS is operated to isolate a 31.5 kb locus cleaved by Cas9 in the genome of the plant Medicago truncatula. Target purification is validated on a Bacterial Artificial Chromosome plasmid, and subsequently carried out in whole genome with µLAS, PFGE or by combining these techniques. PacBio sequencing shows an enrichment factor of the target sequence of 84 with PFGE alone versus 892 by association of PFGE with µLAS. These performances allow us to sequence and assemble one contig of 29 441 bp with 99% sequence identity to the reference sequence.

Whole-genome landscape of Medicago truncatula symbiotic genes.

Advances in deciphering the functional architecture of eukaryotic genomes have been facilitated by recent breakthroughs in sequencing technologies, enabling a more comprehensive representation of genes and repeat elements in genome sequence assemblies, as well as more sensitive and tissue-specific analyses of gene expression. Here we show that PacBio sequencing has led to a substantially improved genome assembly of Medicago truncatula A17, a legume model species notable for endosymbiosis studies1, and has enabled the identification of genome rearrangements between genotypes at a near-base-pair resolution. Annotation of the new M. truncatula genome sequence has allowed for a thorough analysis of transposable elements and their dynamics, as well as the identification of new players involved in symbiotic nodule development, in particular 1,037 upregulated long non-coding RNAs (lncRNAs). We have also discovered that a substantial proportion (~35% and 38%, respectively) of the genes upregulated in nodules or expressed in the nodule differentiation zone colocalize in genomic clusters (270 and 211, respectively), here termed symbiotic islands. These islands contain numerous expressed lncRNA genes and display differentially both DNA methylation and histone marks. Epigenetic regulations and lncRNAs are therefore attractive candidate elements for the orchestration of symbiotic gene expression in the M. truncatula genome.

Laser Capture Micro-Dissection Coupled to RNA Sequencing: A Powerful Approach Applied to the Model Legume Medicago truncatula in Interaction with Sinorhizobium meliloti.

Understanding the development of multicellular organisms requires the identification of regulators, notably transcription factors, and specific transcript populations associated with tissue differentiation. Laser capture microdissection (LCM) is one of the techniques that enable the analysis of distinct tissues or cells within an organ. Coupling this technique with RNA sequencing (RNAseq) makes it extremely powerful to obtain a genome-wide and dynamic view of gene expression. Moreover, RNA sequencing allows two or potentially more interacting organisms to be analyzed simultaneously. In this chapter, a LCM-RNAseq protocol optimized for root and symbiotic root nodule analysis is presented, using the model legume Medicago truncatula (in interaction with Sinorhizobium meliloti in the nodule samples). This includes the description of procedures for plant material fixation, embedding, and micro-dissection; it is followed by a presentation of techniques for RNA extraction and amplification, adapted for the simultaneous analysis of plant and bacterial cells in interaction or, more generally, polyadenylated and non-polyadenylated RNAs. Finally, step-by-step statistical analyses of RNAseq data are described. Those are critical for quality assessment of the whole procedure and for the identification of differentially expressed genes.

Reprogramming of DNA methylation is critical for nodule development in Medicago truncatula.

The legume-Rhizobium symbiosis leads to the formation of a new organ, the root nodule, involving coordinated and massive induction of specific genes. Several genes controlling DNA methylation are spatially regulated within the Medicago truncatula nodule, notably the demethylase gene, DEMETER (DME), which is mostly expressed in the differentiation zone. Here, we show that MtDME is essential for nodule development and regulates the expression of 1,425 genes, some of which are critical for plant and bacterial cell differentiation. Bisulphite sequencing coupled to genomic capture enabled the identification of 474 regions that are differentially methylated during nodule development, including nodule-specific cysteine-rich peptide genes. Decreasing DME expression by RNA interference led to hypermethylation and concomitant downregulation of 400 genes, most of them associated with nodule differentiation. Massive reprogramming of gene expression through DNA demethylation is a new epigenetic mechanism controlling a key stage of indeterminate nodule organogenesis during symbiotic interactions.

A Laser Dissection-RNAseq Analysis Highlights the Activation of Cytokinin Pathways by Nod Factors in the Medicago truncatula Root Epidermis.

Nod factors (NFs) are lipochitooligosaccharidic signal molecules produced by rhizobia, which play a key role in the rhizobium-legume symbiotic interaction. In this study, we analyzed the gene expression reprogramming induced by purified NF (4 and 24 h of treatment) in the root epidermis of the model legume Medicago truncatula Tissue-specific transcriptome analysis was achieved by laser-capture microdissection coupled to high-depth RNA sequencing. The expression of 17,191 genes was detected in the epidermis, among which 1,070 were found to be regulated by NF addition, including previously characterized NF-induced marker genes. Many genes exhibited strong levels of transcriptional activation, sometimes only transiently at 4 h, indicating highly dynamic regulation. Expression reprogramming affected a variety of cellular processes, including perception, signaling, regulation of gene expression, as well as cell wall, cytoskeleton, transport, metabolism, and defense, with numerous NF-induced genes never identified before. Strikingly, early epidermal activation of cytokinin (CK) pathways was indicated, based on the induction of CK metabolic and signaling genes, including the CRE1 receptor essential to promote nodulation. These transcriptional activations were independently validated using promoter:ß-glucuronidase fusions with the MtCRE1 CK receptor gene and a CK response reporter (TWO COMPONENT SIGNALING SENSOR NEW). A CK pretreatment reduced the NF induction of the EARLY NODULIN11 (ENOD11) symbiotic marker, while a CK-degrading enzyme (CYTOKININ OXIDASE/DEHYDROGENASE3) ectopically expressed in the root epidermis led to increased NF induction of ENOD11 and nodulation. Therefore, CK may play both positive and negative roles in M. truncatula nodulation.

Full-length de novo assembly of RNA-seq data in pea (Pisum sativum L.) provides a gene expression atlas and gives insights into root nodulation in this species.

Next-generation sequencing technologies allow an almost exhaustive survey of the transcriptome, even in species with no available genome sequence. To produce a Unigene set representing most of the expressed genes of pea, 20 cDNA libraries produced from various plant tissues harvested at various developmental stages from plants grown under contrasting nitrogen conditions were sequenced. Around one billion reads and 100 Gb of sequence were de novo assembled. Following several steps of redundancy reduction, 46 099 contigs with N50 length of 1667 nt were identified. These constitute the 'Caméor' Unigene set. The high depth of sequencing allowed identification of rare transcripts and detected expression for approximately 80% of contigs in each library. The Unigene set is now available online (http://bios.dijon.inra.fr/FATAL/cgi/pscam.cgi), allowing (i) searches for pea orthologs of candidate genes based on gene sequences from other species, or based on annotation, (ii) determination of transcript expression patterns using various metrics, (iii) identification of uncharacterized genes with interesting patterns of expression, and (iv) comparison of gene ontology pathways between tissues. This resource has allowed identification of the pea orthologs of major nodulation genes characterized in recent years in model species, as a major step towards deciphering unresolved pea nodulation phenotypes. In addition to a remarkable conservation of the early transcriptome nodulation apparatus between pea and Medicago truncatula, some specific features were highlighted. The resource provides a reference for the pea exome, and will facilitate transcriptome and proteome approaches as well as SNP discovery in pea.

An integrated analysis of plant and bacterial gene expression in symbiotic root nodules using laser-capture microdissection coupled to RNA sequencing.

Rhizobium-induced root nodules are specialized organs for symbiotic nitrogen fixation. Indeterminate-type nodules are formed from an apical meristem and exhibit a spatial zonation which corresponds to successive developmental stages. To get a dynamic and integrated view of plant and bacterial gene expression associated with nodule development, we used a sensitive and comprehensive approach based upon oriented high-depth RNA sequencing coupled to laser microdissection of nodule regions. This study, focused on the association between the model legume Medicago truncatula and its symbiont Sinorhizobium meliloti, led to the production of 942 million sequencing read pairs that were unambiguously mapped on plant and bacterial genomes. Bioinformatic and statistical analyses enabled in-depth comparison, at a whole-genome level, of gene expression in specific nodule zones. Previously characterized symbiotic genes displayed the expected spatial pattern of expression, thus validating the robustness of our approach. We illustrate the use of this resource by examining gene expression associated with three essential elements of nodule development, namely meristem activity, cell differentiation and selected signaling processes related to bacterial Nod factors and redox status. We found that transcription factor genes essential for the control of the root apical meristem were also expressed in the nodule meristem, while the plant mRNAs most enriched in nodules compared with roots were mostly associated with zones comprising both plant and bacterial partners. The data, accessible on a dedicated website, represent a rich resource for microbiologists and plant biologists to address a variety of questions of both fundamental and applied interest.

The symbiotic transcription factor MtEFD and cytokinins are positively acting in the Medicago truncatula and Ralstonia solanacearum pathogenic interaction.

• A plant-microbe dual biological system was set up involving the model legume Medicago truncatula and two bacteria, the soil-borne root pathogen Ralstonia solanacearum and the beneficial symbiont Sinorhizobium meliloti. • Comparison of transcriptomes under symbiotic and pathogenic conditions highlighted the transcription factor MtEFD (Ethylene response Factor required for nodule Differentiation) as being upregulated in both interactions, together with a set of cytokinin-related transcripts involved in metabolism, signaling and response. MtRR4 (Response Regulator), a cytokinin primary response gene negatively regulating cytokinin signaling and known as a target of MtEFD in nodulation processes, was retrieved in this set of transcripts. • Refined studies of MtEFD and MtRR4 expression during M. truncatula and R. solanacearum interaction indicated differential kinetics of induction and requirement of central regulators of bacterial pathogenicity, HrpG and HrpB. Similar to MtRR4, MtEFD upregulation during the pathogenic interaction was dependent on cytokinin perception mediated by the MtCRE1 (Cytokinin REsponse 1) receptor. • The use of M. truncatula efd-1 and cre1-1 mutants evidenced MtEFD and cytokinin perception as positive factors for bacterial wilt development. These factors therefore play an important role in both root nodulation and root disease development.

The CCAAT box-binding transcription factor NF-YA1 controls rhizobial infection.

Symbiosis between legume plants and soil rhizobia culminates in the formation of a novel root organ, the 'nodule', containing bacteria differentiated as facultative nitrogen-fixing organelles. MtNF-YA1 is a Medicago truncatula CCAAT box-binding transcription factor (TF), formerly called HAP2-1, highly expressed in mature nodules and required for nodule meristem function and persistence. Here a role for MtNF-YA1 during early nodule development is demonstrated. Detailed expression analysis based on RNA sequencing, quantitiative real-time PCR (qRT-PCR), as well as promoter-ß-glucuronidase (GUS) fusions reveal that MtNF-YA1 is first induced at the onset of symbiotic development during preparation for, and initiation and progression of, symbiotic infection. Moreover, using a new knock-out mutant, Mtnf-ya1-1, it is shown that MtNF-YA1 controls infection thread (IT) progression from initial root infection through colonization of nodule tissues. Extensive confocal and electronic microscopic observations suggest that the bulbous and erratic IT growth phenotypes observed in Mtnf-ya1-1 could be a consequence of the fact that walls of ITs in this mutant are thinner and less coherent than in the wild type. It is proposed that MtNF-YA1 controls rhizobial infection progression by regulating the formation and the wall of ITs.

MtbHLH1, a bHLH transcription factor involved in Medicago truncatula nodule vascular patterning and nodule to plant metabolic exchanges.

This study aimed at defining the role of a basic helix-loop-helix (bHLH) transcription factor gene from Medicago truncatula, MtbHLH1, whose expression is upregulated during the development of root nodules produced upon infection by rhizobia bacteria. We used MtbHLH1 promoter::GUS fusions and quantitative reverse-transcription polymerase chain reaction analyses to finely characterize the MtbHLH1 expression pattern. We altered MtbHLH1 function by expressing a dominantly repressed construct (CRES-T approach) and looked for possible MtbHLH1 target genes by transcriptomics. We found that MtbHLH1 is expressed in nodule primordia cells derived from pericycle divisions, in nodule vascular bundles (VBs) and in uninfected cells of the nitrogen (N) fixation zone. MtbHLH1 is also expressed in root tips, lateral root primordia cells and root VBs, and induced upon auxin treatment. Altering MtbHLH1 function led to an unusual phenotype, with a modified patterning of nodule VB development and a reduced growth of aerial parts of the plant, even though the nodules were able to fix atmospheric N. Several putative MtbHLH1 regulated genes were identified, including an asparagine synthase and a LOB (lateral organ boundary) transcription factor. Our results suggest that the MtbHLH1 gene is involved in the control of nodule vasculature patterning and nutrient exchanges between nodules and roots.

Transcription reprogramming during root nodule development in Medicago truncatula.

Many genes which are associated with root nodule development and activity in the model legume Medicago truncatula have been described. However information on precise stages of activation of these genes and their corresponding transcriptional regulators is often lacking. Whether these regulators are shared with other plant developmental programs also remains an open question. Here detailed microarray analyses have been used to study the transcriptome of root nodules induced by either wild type or mutant strains of Sinorhizobium meliloti. In this way we have defined eight major activation patterns in nodules and identified associated potential regulatory genes. We have shown that transcription reprogramming during consecutive stages of nodule differentiation occurs in four major phases, respectively associated with (i) early signalling events and/or bacterial infection; plant cell differentiation that is either (ii) independent or (iii) dependent on bacteroid differentiation; (iv) nitrogen fixation. Differential expression of several genes involved in cytokinin biosynthesis was observed in early symbiotic nodule zones, suggesting that cytokinin levels are actively controlled in this region. Taking advantage of databases recently developed for M. truncatula, we identified a small subset of gene expression regulators that were exclusively or predominantly expressed in nodules, whereas most other regulators were also activated under other conditions, and notably in response to abiotic or biotic stresses. We found evidence suggesting the activation of the jasmonate pathway in both wild type and mutant nodules, thus raising questions about the role of jasmonate during nodule development. Finally, quantitative RT-PCR was used to analyse the expression of a series of nodule regulator and marker genes at early symbiotic stages in roots and allowed us to distinguish several early stages of gene expression activation or repression.

The Medicago truncatula E3 ubiquitin ligase PUB1 interacts with the LYK3 symbiotic receptor and negatively regulates infection and nodulation.

LYK3 is a lysin motif receptor-like kinase of Medicago truncatula, which is essential for the establishment of the nitrogen-fixing, root nodule symbiosis with Sinorhizobium meliloti. LYK3 is a putative receptor of S. meliloti Nod factor signals, but little is known of how it is regulated and how it transduces these symbiotic signals. In a screen for LYK3-interacting proteins, we identified M. truncatula Plant U-box protein 1 (PUB1) as an interactor of the kinase domain. In planta, both proteins are localized and interact in the plasma membrane. In M. truncatula, PUB1 is expressed specifically in symbiotic conditions, is induced by Nod factors, and shows an overlapping expression pattern with LYK3 during nodulation. Biochemical studies show that PUB1 has a U-box-dependent E3 ubiquitin ligase activity and is phosphorylated by the LYK3 kinase domain. Overexpression and RNA interference studies in M. truncatula show that PUB1 is a negative regulator of the LYK3 signaling pathway leading to infection and nodulation and is important for the discrimination of rhizobia strains producing variant Nod factors. The potential role of PUB E3 ubiquitin ligases in controlling plant-microbe interactions and development through interacting with receptor-like kinases is discussed.

A remorin protein interacts with symbiotic receptors and regulates bacterial infection.

Remorin proteins have been hypothesized to play important roles during cellular signal transduction processes. Induction of some members of this multigene family has been reported during biotic interactions. However, no roles during host-bacteria interactions have been assigned to remorin proteins until now. We used root nodule symbiosis between Medicago truncatula and Sinorhizobium meliloti to study the roles of a remorin that is specifically induced during nodulation. Here we show that this oligomeric remorin protein attaches to the host plasma membrane surrounding the bacteria and controls infection and release of rhizobia into the host cytoplasm. It interacts with the core set of symbiotic receptors that are essential for perception of bacterial signaling molecules, and thus might represent a plant-specific scaffolding protein.

Analysis and modeling of the integrative response of Medicago truncatula to nitrogen constraints.

An integrative biology approach was conducted in Medicago truncatula for: (i) unraveling the coordinated regulation of NO3-, NH4+ and N(2) acquisition by legumes to fulfill the plant N demand; and (ii) modeling the emerging properties occurring at the whole plant level. Upon localized addition of a high level of mineral N, the three N acquisition pathways displayed similar systemic feedback repression to adjust N acquisition capacities to the plant N status. Genes associated to these responses were in contrast rather specific to the N source. Following an N deficit, NO3- fed plants maintained efficiently their N status through rapid functional and developmental up regulations while N(2) fed plants responded by long term plasticity of nodule development. Regulatory genes associated with various symbiotic stages were further identified. An ecophysiological model simulating relations between leaf area and roots N retrieval was developed and now furnishes an analysis grid to characterize a spontaneous or induced genetic variability for plant N nutrition.

EFD Is an ERF transcription factor involved in the control of nodule number and differentiation in Medicago truncatula.

Mechanisms regulating legume root nodule development are still poorly understood, and very few regulatory genes have been cloned and characterized. Here, we describe EFD (for ethylene response factor required for nodule differentiation), a gene that is upregulated during nodulation in Medicago truncatula. The EFD transcription factor belongs to the ethylene response factor (ERF) group V, which contains ERN1, 2, and 3, three ERFs involved in Nod factor signaling. The role of EFD in the regulation of nodulation was examined through the characterization of a null deletion mutant (efd-1), RNA interference, and overexpression studies. These studies revealed that EFD is a negative regulator of root nodulation and infection by Rhizobium and that EFD is required for the formation of functional nitrogen-fixing nodules. EFD appears to be involved in the plant and bacteroid differentiation processes taking place beneath the nodule meristem. We also showed that EFD activated Mt RR4, a cytokinin primary response gene that encodes a type-A response regulator. We propose that EFD induction of Mt RR4 leads to the inhibition of cytokinin signaling, with two consequences: the suppression of new nodule initiation and the activation of differentiation as cells leave the nodule meristem. Our work thus reveals a key regulator linking early and late stages of nodulation and suggests that the regulation of the cytokinin pathway is important both for nodule initiation and development.

A gene expression atlas of the model legume Medicago truncatula.

Legumes played central roles in the development of agriculture and civilization, and today account for approximately one-third of the world's primary crop production. Unfortunately, most cultivated legumes are poor model systems for genomic research. Therefore, Medicago truncatula, which has a relatively small diploid genome, has been adopted as a model species for legume genomics. To enhance its value as a model, we have generated a gene expression atlas that provides a global view of gene expression in all major organ systems of this species, with special emphasis on nodule and seed development. The atlas reveals massive differences in gene expression between organs that are accompanied by changes in the expression of key regulatory genes, such as transcription factor genes, which presumably orchestrate genetic reprogramming during development and differentiation. Interestingly, many legume-specific genes are preferentially expressed in nitrogen-fixing nodules, indicating that evolution endowed them with special roles in this unique and important organ. Comparative transcriptome analysis of Medicago versus Arabidopsis revealed significant divergence in developmental expression profiles of orthologous genes, which indicates that phylogenetic analysis alone is insufficient to predict the function of orthologs in different species. The data presented here represent an unparalleled resource for legume functional genomics, which will accelerate discoveries in legume biology.