Epigenetic control during root development and symbiosis.

The roots of plants play multiples functions that are essential for growth and development, including anchoring to the soil and water and nutrient acquisition. These underground organs exhibit the plasticity to modify their root system architecture in response to environmental cues allowing adaptation to change in water and nutrient availability. In addition, roots enter in mutualistic interactions with soil microorganisms, e.g. the root nodule symbiosis established between a limited group of plants and nitrogen fixing soil bacteria and the arbuscular mycorrhiza symbiosis involving most land plants and fungi of the Glomeromycetes phylum. In the past 20 years, genetic approaches allowed the identification and functional characterization of genes required for the specific programs of root development, root nodule and arbuscular mycorrhiza symbioses. These genetic studies provided evidence that the program of root nodule symbiosis recruited components of the arbuscular mycorrhiza symbiosis and the root developmental programs. The execution of these programs is strongly influenced by epigenetic changes -DNA methylation and histone post-translational modifications- that alter chromatin conformation modifying the expression of key genes. In this review, we summarize recent advances that highlighted how DNA methylation and histone post-translational modifications, as well as chromatin remodeling factors and long non-coding RNAs, shape the root system architecture and allow the successful establishment of both root nodule and arbuscular mycorrhiza symbioses. We anticipate that the analysis of dynamic epigenetic changes and chromatin 3D structure in specific single-cells or tissue types of root organs will illuminate our understanding of how root developmental and symbiotic programs are orchestrated, opening exciting questions and new perspectives to modulate agronomical and ecological traits linked to nutrient acquisition.

A lateral organ boundaries domain transcription factor acts downstream of the auxin response factor 2 to control nodulation and root architecture in Medicago truncatula.

Legume plants develop two types of root postembryonic organs, lateral roots and symbiotic nodules, using shared regulatory components. The module composed by the microRNA390, the Trans-Acting SIRNA3 (TAS3) RNA and the Auxin Response Factors (ARF)2, ARF3, and ARF4 (miR390/TAS3/ARFs) mediates the control of both lateral roots and symbiotic nodules in legumes. Here, a transcriptomic approach identified a member of the Lateral Organ Boundaries Domain (LBD) family of transcription factors in Medicago truncatula, designated MtLBD17/29a, which is regulated by the miR390/TAS3/ARFs module. ChIP-PCR experiments evidenced that MtARF2 binds to an Auxin Response Element present in the MtLBD17/29a promoter. MtLBD17/29a is expressed in root meristems, lateral root primordia, and noninfected cells of symbiotic nodules. Knockdown of MtLBD17/29a reduced the length of primary and lateral roots and enhanced lateral root formation, whereas overexpression of MtLBD17/29a produced the opposite phenotype. Interestingly, both knockdown and overexpression of MtLBD17/29a reduced nodule number and infection events and impaired the induction of the symbiotic genes Nodulation Signaling Pathway (NSP) 1 and 2. Our results demonstrate that MtLBD17/29a is regulated by the miR390/TAS3/ARFs module and a direct target of MtARF2, revealing a new lateral root regulatory hub recruited by legumes to act in the root nodule symbiotic program.

Nuclear and cytoplasmic lncRNAs in root tips of the model legume Medicago truncatula under control and submergence.

In this study, we aimed to identify long noncoding RNAs (lncRNAs) in root tips of the model legume Medicago truncatula using previously generated nuclear, total polyA, ribosome-associated polyA, and Riboseq RNA datasets, which might shed light on their localization and potential regulatory roles. RNA-seq data were mapped to the version 5 of the M. truncatula A17 genome and analyzed to identify genome annotated lncRNAs and putative new root tip (NRT) lncRNAs. lncRNAs were classified according to their genomic location relative to chromatin accessible regions, protein-coding genes and transposable elements (TE), finding differences between annotated lncRNAs and NRT lncRNAs, both in their genomic position as well as in the type of TEs in their vicinity. We investigated their response to submergence and found a set of regulated lncRNAs that were preferentially upregulated in the nucleus, some of which were located nearby genes of the conserved submergence upregulated gene families, and chromatin accessible regions suggesting a potential regulatory role. Finally, the accumulation of lncRNAs under submergence was validated by reverse transcription quantitative polymerase chain reaction on nuclear RNA, providing additional evidence of their localization, which could ultimately be required for their function.

NIPK, a protein pseudokinase that interacts with the C subunit of the transcription factor NF-Y, is involved in rhizobial infection and nodule organogenesis.

Heterotrimeric Nuclear Factor Y (NF-Y) transcription factors are key regulators of the symbiotic program that controls rhizobial infection and nodule organogenesis. Using a yeast two-hybrid screening, we identified a putative protein kinase of Phaseolus vulgaris that interacts with the C subunit of the NF-Y complex. Physical interaction between NF-YC1 Interacting Protein Kinase (NIPK) and NF-YC1 occurs in the cytoplasm and the plasma membrane. Only one of the three canonical amino acids predicted to be required for catalytic activity is conserved in NIPK and its putative homologs from lycophytes to angiosperms, indicating that NIPK is an evolutionary conserved pseudokinase. Post-transcriptional silencing on NIPK affected infection and nodule organogenesis, suggesting NIPK is a positive regulator of the NF-Y transcriptional complex. In addition, NIPK is required for activation of cell cycle genes and early symbiotic genes in response to rhizobia, including NF-YA1 and NF-YC1. However, strain preference in co-inoculation experiments was not affected by NIPK silencing, suggesting that some functions of the NF-Y complex are independent of NIPK. Our work adds a new component associated with the NF-Y transcriptional regulators in the context of nitrogen-fixing symbiosis.

Identification of conserved and new miRNAs that affect nodulation and strain selectivity in the Phaseolus vulgaris-Rhizobium etli symbiosis through differential analysis of host small RNAs.

Phaseolus vulgaris plants from the Mesoamerican centre of genetic diversification establish a preferential and more efficient root nodule symbiosis with sympatric Rhizobium etli strains. This is mediated by changes in host gene expression, which might occur either at the transcriptional or at the post-transcriptional level. However, the implication of small RNA (sRNA)-mediated control of gene expression in strain selectivity has remained elusive. sRNA sequencing was used to identify host microRNAs (miRNAs) differentially regulated in roots at an early stage of the symbiotic interaction, which were further characterized by applying a reverse genetic approach. In silico analysis identified known and new miRNAs that accumulated to a greater extent in the preferential and more efficient interaction. One of them, designated as Pvu-miR5924, participates in the mechanisms that determine the selection of R. etli strains that will colonize the nodules. In addition, the functional analysis of Pvu-miR390b verified that this miRNA is a negative modulator of nodule formation and bacterial infection. This study not only extended the list of miRNAs identified in P. vulgaris but also enabled the identification of miRNAs that play relevant functions in nodule formation, rhizobial infection and the selection of the rhizobial strains that will occupy the nodule.

Transcriptomic analysis of Mesoamerican and Andean Phaseolus vulgaris accessions revealed mRNAs and lncRNAs associated with strain selectivity during symbiosis.

Legume plants establish a nitrogen-fixing symbiosis with soil bacteria known as rhizobia. Compatibility between legumes and rhizobia is determined at species-specific level, but variations in the outcome of the symbiotic process are also influenced by the capacity of the plant to discriminate and select specific strains that are better partners. We compared the transcriptional response of two genetically diverse accessions of Phaseolus vulgaris from Mesoamerica and South Andes to Rhizobium etli strains that exhibit variable degrees of symbiotic affinities. Our results indicate that the plant genotype is the major determinant of the transcriptional reprogramming occurring in roots at early stages of the symbiotic interaction. Differentially expressed genes (DEGs) regulated in the Mesoamerican and the Andean accessions in response to specific strains are different, but they belong to the same functional categories. The common and strain-specific transcriptional responses to rhizobia involve distinct transcription factors and cis-elements present in the promoters of DEGs in each accession, showing that diversification and domestication of common bean at different geographic regions influenced the evolution of symbiosis differently in each genetic pool. Quantitative PCR analysis validated our transcriptional datasets, which constitute a valuable source of coding and non-coding candidate genes to further unravel the molecular determinants governing the mechanisms by which plants select bacterial strains that produce a better symbiotic outcome.

Translational regulation in pathogenic and beneficial plant-microbe interactions.

Plants are surrounded by a vast diversity of microorganisms. Limiting pathogenic microorganisms is crucial for plant survival. On the other hand, the interaction of plants with beneficial microorganisms promotes their growth or allows them to overcome nutrient deficiencies. Balancing the number and nature of these interactions is crucial for plant growth and development, and thus, for crop productivity in agriculture. Plants use sophisticated mechanisms to recognize pathogenic and beneficial microorganisms and genetic programs related to immunity or symbiosis. Although most research has focused on characterizing changes in the transcriptome during plant-microbe interactions, the application of techniques such as Translating Ribosome Affinity Purification (TRAP) and Ribosome profiling allowed examining the dynamic association of RNAs to the translational machinery, highlighting the importance of the translational level of control of gene expression in both pathogenic and beneficial interactions. These studies revealed that the transcriptional and the translational responses are not always correlated, and that translational control operates at cell-specific level. In addition, translational control is governed by cis-elements present in the 5'mRNA leader of regulated mRNAs, e.g. upstream open reading frames (uORFs) and sequence-specific motifs. In this review, we summarize and discuss the recent advances made in the field of translational control during pathogenic and beneficial plant-microbe interactions.

Auxin Response Factor 2 (ARF2), ARF3, and ARF4 Mediate Both Lateral Root and Nitrogen Fixing Nodule Development in Medicago truncatula.

Auxin Response Factors (ARFs) constitute a large family of transcription factors that mediate auxin-regulated developmental programs in plants. ARF2, ARF3, and ARF4 are post-transcriptionally regulated by the microRNA390 (miR390)/trans-acting small interference RNA 3 (TAS3) module through the action of TAS3-derived trans - acting small interfering RNAs (ta-siRNA). We have previously reported that constitutive activation of the miR390/TAS3 pathway promotes elongation of lateral roots but impairs nodule organogenesis and infection by rhizobia during the nitrogen-fixing symbiosis established between Medicago truncatula and its partner Sinorhizobium meliloti. However, the involvement of the targets of the miR390/TAS3 pathway, i.e., MtARF2, MtARF3, MtARF4a, and MtARF4b, in root development and establishment of the nitrogen-fixing symbiosis remained unexplored. Here, promoter:reporter fusions showed that expression of both MtARF3 and MtARF4a was associated with lateral root development; however, only the MtARF4a promoter was active in developing nodules. In addition, up-regulation of MtARF2, MtARF3, and MtARF4a/b in response to rhizobia depends on Nod Factor perception. We provide evidence that simultaneous knockdown of MtARF2, MtARF3, MtARF4a, and MtARF4b or mutation in MtARF4a impaired nodule formation, and reduced initiation and progression of infection events. Silencing of MtARF2, MtARF3, MtARF4a, and MtARF4b altered mRNA levels of the early nodulation gene nodulation signaling pathway 2 (MtNSP2). In addition, roots with reduced levels of MtARF2, MtARF3, MtARF4a, and MtARF4b, as well as arf4a mutant plants exhibited altered root architecture, causing a reduction in primary and lateral root length, but increasing lateral root density. Taken together, our results suggest that these ARF members are common key players of the morphogenetic programs that control root development and the formation of nitrogen-fixing nodules.

TRAP-SEQ of Eukaryotic Translatomes Applied to the Detection of Polysome-Associated Long Noncoding RNAs.

Translating ribosome affinity purification (TRAP) technology allows the isolation of polysomal complexes and the RNAs associated with at least one 80S ribosome. TRAP consists of the stabilization and affinity purification of polysomes containing a tagged version of a ribosomal protein. Quantitative assessment of the TRAP RNA is achieved by direct sequencing (TRAP-SEQ), which provides accurate quantitation of ribosome-associated RNAs, including long noncoding RNAs (lncRNAs). Here we present an updated procedure for TRAP-SEQ, as well as a primary analysis guide for identification of ribosome-associated lncRNAs. This methodology enables the study of dynamic association of lncRNAs by assessing rapid changes in their transcript levels in polysomes at organ or cell-type level, during development, or in response to endogenous or exogenous stimuli.

To keep or not to keep: mRNA stability and translatability in root nodule symbiosis.

Post-transcriptional control of gene expression allows plants to rapidly adapt to changes in their environment. Under low nitrogen conditions, legume plants engage into a symbiosis with soil bacteria that results in the formation of root nodules, where bacteria are allocated and fix atmospheric nitrogen for the plant's benefit. Recent studies highlighted the importance of small RNA-mediated mechanisms in the control of bacterial infection, nodule organogenesis, and the long-distance signaling that balances plant growth and nodulation. Examples of such mechanisms are shoot-to-root mobile microRNAs and small RNA fragments derived from degradation of bacterial transfer RNAs that repress complementary mRNAs in the host plant. Mechanisms of selective mRNA translation also contribute to rapidly modulate the expression of nodulation genes in a cell-specific manner during symbiosis. Here, the most recent advances made on the regulation of mRNA stability and translatability, and the emerging roles of long non-coding RNAs in symbiosis are summarized.

Reprogramming of Root Cells during Nitrogen-Fixing Symbiosis Involves Dynamic Polysome Association of Coding and Noncoding RNAs.

Translational control is a widespread mechanism that allows the cell to rapidly modulate gene expression in order to provide flexibility and adaptability to eukaryotic organisms. We applied translating ribosome affinity purification combined with RNA sequencing to characterize translational regulation of mRNAs at early stages of the nitrogen-fixing symbiosis established between Medicago truncatula and Sinorhizobium meliloti Our analysis revealed a poor correlation between transcriptional and translational changes and identified hundreds of regulated protein-coding and long noncoding RNAs (lncRNAs), some of which are regulated in specific cell types. We demonstrated that a short variant of the lncRNA Trans-acting small interference RNA3 (TAS3) increased its association to the translational machinery in response to rhizobia. Functional analysis revealed that this short variant of TAS3 might act as a target mimic that captures microRNA390, contributing to reduce trans acting small interference Auxin Response Factor production and modulating nodule formation and rhizobial infection. The analysis of alternative transcript variants identified a translationally upregulated mRNA encoding subunit 3 of the SUPERKILLER complex (SKI3), which participates in mRNA decay. Knockdown of SKI3 decreased nodule initiation and development, as well as the survival of bacteria within nodules. Our results highlight the importance of translational control and mRNA decay pathways for the successful establishment of the nitrogen-fixing symbiosis.

The PvNF-YA1 and PvNF-YB7 Subunits of the Heterotrimeric NF-Y Transcription Factor Influence Strain Preference in the Phaseolus vulgaris-Rhizobium etli Symbiosis.

Transcription factors of the Nuclear Factor Y (NF-Y) family play essential functions in plant development and plasticity, including the formation of lateral root organs such as lateral root and symbiotic nodules. NF-Ys mediate transcriptional responses by acting as heterotrimers composed of three subunits, NF-YA, NF-YB, and NF-YC, which in plants are encoded by relatively large gene families. We have previously shown that, in the Phaseolus vulgaris × Rhizobium etli interaction, the PvNF-YC1 subunit is involved not only in the formation of symbiotic nodules, but also in the preference exhibited by the plant for rhizobial strains that are more efficient and competitive in nodule formation. PvNF-YC1 forms a heterotrimer with the PvNF-YA1 and PvNF-YB7 subunits. Here, we used promoter:reporter fusions to show that both PvNF-YA1 and PvNF-YB7 are expressed in symbiotic nodules. In addition, we report that knock-down of PvNF-YA1 and its close paralog PvNF-YA9 abolished nodule formation by either high or low efficient strains and arrested rhizobial infection. On the other hand, knock-down of PvNF-YB7 only affected the symbiotic outcome of the high efficient interaction, suggesting that other symbiotic NF-YB subunits might be involved in the more general mechanisms of nodule formation. More important, we present functional evidence supporting that both PvNF-YA1 and PvNF-YB7 are part of the mechanisms that allow P. vulgaris plants to discriminate and select those bacterial strains that perform better in nodule formation, most likely by acting in the same heterotrimeric complex that PvNF-YC1.

Small GTPases in plant biotic interactions.

The superfamily of small monomeric GTPases originated in a common ancestor of eukaryotic multicellular organisms and, since then, it has evolved independently in each lineage to cope with the environmental challenges imposed by their different life styles. Members of the small GTPase family function in the control of vesicle trafficking, cytoskeleton rearrangements and signaling during crucial biological processes, such as cell growth and responses to environmental cues. In this review, we discuss the emerging roles of these small GTPases in the pathogenic and symbiotic interactions established by plants with microorganisms present in their nearest environment, in which membrane trafficking is crucial along the different steps of the interaction, from recognition and signal transduction to nutrient exchange.

Compatibility between Legumes and Rhizobia for the Establishment of a Successful Nitrogen-Fixing Symbiosis.

The root nodule symbiosis established between legumes and rhizobia is an exquisite biological interaction responsible for fixing a significant amount of nitrogen in terrestrial ecosystems. The success of this interaction depends on the recognition of the right partner by the plant within the richest microbial ecosystems on Earth, the soil. Recent metagenomic studies of the soil biome have revealed its complexity, which includes microorganisms that affect plant fitness and growth in a beneficial, harmful, or neutral manner. In this complex scenario, understanding the molecular mechanisms by which legumes recognize and discriminate rhizobia from pathogens, but also between distinct rhizobia species and strains that differ in their symbiotic performance, is a considerable challenge. In this work, we will review how plants are able to recognize and select symbiotic partners from a vast diversity of surrounding bacteria. We will also analyze recent advances that contribute to understand changes in plant gene expression associated with the outcome of the symbiotic interaction. These aspects of nitrogen-fixing symbiosis should contribute to translate the knowledge generated in basic laboratory research into biotechnological advances to improve the efficiency of the nitrogen-fixing symbiosis in agronomic systems.

Comparative phylogenetic and expression analysis of small GTPases families in legume and non-legume plants.

BACKGROUND: Small monomeric GTPases act as molecular switches in several processes that involve polar cell growth, participating mainly in vesicle trafficking and cytoskeleton rearrangements. This gene superfamily has largely expanded in plants through evolution as compared with other Kingdoms, leading to the suggestion that members of each subfamily might have acquired new functions associated to plant-specific processes. Legume plants engage in a nitrogen-fixing symbiotic interaction with rhizobia in a process that involves polar growth processes associated with the infection throughout the root hair. To get insight into the evolution of small GTPases associated with this process, we use a comparative genomic approach to establish differences in the Ras GTPase superfamily between legume and non-legume plants. RESULTS: Phylogenetic analyses did not show clear differences in the organization of the different subfamilies of small GTPases between plants that engage or not in nodule symbiosis. Protein alignments revealed a strong conservation at the sequence level of small GTPases previously linked to nodulation by functional genetics. Interestingly, one Rab and three Rop proteins showed conserved amino acid substitutions in legumes, but these changes do not alter the predicted conformational structure of these proteins. Although the steady-state levels of most small GTPases do not change in response to rhizobia, we identified a subset of Rab, Rop and Arf genes whose transcript levels are modulated during the symbiotic interaction, including their spatial distribution along the indeterminate nodule. CONCLUSIONS: This study provides a comprehensive study of the small GTPase superfamily in several plant species. The genetic program associated to root nodule symbiosis includes small GTPases to fulfill specific functions during infection and formation of the symbiosomes. These GTPases seems to have been recruited from members that were already present in common ancestors with plants as distant as monocots since we failed to detect asymmetric evolution in any of the subfamily trees. Expression analyses identified a number of legume members that can have undergone neo- or sub-functionalization associated to the spatio-temporal transcriptional control during the onset of the symbiotic interaction.

The MicroRNA390/TAS3 Pathway Mediates Symbiotic Nodulation and Lateral Root Growth.

Legume roots form two types of postembryonic organs, lateral roots and symbiotic nodules. Nodule formation is the result of the interaction of legumes with rhizobia and requires the mitotic activation and differentiation of root cells as well as an independent, but coordinated, program that allows infection by rhizobia. MicroRNA390 (miR390) is an evolutionarily conserved microRNA that targets the Trans-Acting Short Interference RNA3 (TAS3) transcript. Cleavage of TAS3 by ARGONAUTE7 results in the production of trans-acting small interference RNAs, which target mRNAs encoding AUXIN RESPONSE FACTOR2 (ARF2), ARF3, and ARF4. Here, we show that activation of the miR390/TAS3 regulatory module by overexpression of miR390 in Medicago truncatula promotes lateral root growth but prevents nodule organogenesis, rhizobial infection, and the induction of two key nodulation genes, Nodulation Signaling Pathway1 (NSP1) and NSP2 Accordingly, inactivation of the miR390/TAS3 module, either by expression of a miR390 target mimicry construct or mutations in ARGONAUTE7, enhances nodulation and rhizobial infection, alters the spatial distribution of the nodules, and increases the percentage of nodules with multiple meristems. Our results revealed a key role of the miR390/TAS3 pathway in legumes as a modulator of lateral root organs, playing opposite roles in lateral root and nodule development.

The monomeric GTPase RabA2 is required for progression and maintenance of membrane integrity of infection threads during root nodule symbiosis.

KEY MESSAGE: Progression of the infection canal that conducts rhizobia to the nodule primordium requires a functional Rab GTPase located in Golgi/trans-Golgi that also participate in root hair polar growth. Common bean (Phaseolus vulgaris) symbiotically associates with its partner Rhizobium etli, resulting in the formation of root nitrogen-fixing nodules. Compatible bacteria can reach cortical cells in a tightly regulated infection process, in which the specific recognition of signal molecules is a key step to select the symbiotic partner. In this work, we show that RabA2, a monomeric GTPase from common bean, is required for the progression of the infection canal, referred to as the infection thread (IT), toward the cortical cells. Expression of miss-regulated mutant variants of RabA2 resulted in an increased number of abortive infection events, including bursting of ITs and a reduction in the number of nodules. Nodules formed in these plants were small and contained infected cells with disrupted symbiosome membranes, indicating either early senescence of these cells or defects in the formation of the symbiosome membrane during bacterial release. RabA2 localized to mobile vesicles around the IT, but mutations that affect GTP hydrolysis or GTP/GDP exchange modified this localization. Colocalization of RabA2 with ArfA1 and a Golgi marker indicates that RabA2 localizes in Golgi stacks and the trans-Golgi network. Our results suggest that RabA2 is part of the vesicle transport events required to maintain the integrity of the membrane during IT progression.

Plant NF-Y transcription factors: Key players in plant-microbe interactions, root development and adaptation to stress.

NF-Ys are heterotrimeric transcription factors composed by the NF-YA, NF-YB and NF-YC subunits. In plants, NF-Y subunits are encoded by multigene families whose members show structural and functional diversifications. An increasing number of NF-Y genes has been shown to play key roles during different stages of root nodule and arbuscular mycorrhizal symbiosis, as well as during the interaction of plants with pathogenic microorganisms. Individual members of the NF-YA and NF-YB families have also been implicated in the development of primary and lateral roots. In addition, different members of the NF-YA and NF-YB gene families from mono- and di-cotyledonous plants have been involved in plant responses to water and nutrient scarcity. This review presents the most relevant and striking results concerning these NF-Y subunits. A phylogenetic analysis of the functionally characterized NF-Y genes revealed that, across plant species, NF-Y proteins functioning in the same biological process tend to belong to common phylogenetic groups. Finally, we discuss the forthcoming challenges of plant NF-Y research, including the detailed dissection of expression patterns, the elucidation of functional specificities as well as the characterization of the potential NF-Y-mediated epigenetic mechanisms by which they control the expression of their target genes. This article is part of a Special Issue entitled: Nuclear Factor Y in Development and Disease, edited by Prof. Roberto Mantovani.