Unzipped genome assemblies of polyploid root-knot nematodes reveal unusual and clade-specific telomeric repeats.

Using long-read sequencing, we assembled and unzipped the polyploid genomes of Meloidogyne incognita, M. javanica and M. arenaria, three of the most devastating plant-parasitic nematodes. We found the canonical nematode telomeric repeat to be missing in these and other Meloidogyne genomes. In addition, we find no evidence for the enzyme telomerase or for orthologs of C. elegans telomere-associated proteins, suggesting alternative lengthening of telomeres. Instead, analyzing our assembled genomes, we identify species-specific composite repeats enriched mostly at one extremity of contigs. These repeats are G-rich, oriented, and transcribed, similarly to canonical telomeric repeats. We confirm them as telomeric using fluorescent in situ hybridization. These repeats are mostly found at one single end of chromosomes in these species. The discovery of unusual and specific complex telomeric repeats opens a plethora of perspectives and highlights the evolutionary diversity of telomeres despite their central roles in senescence, aging, and chromosome integrity.

AUXIN RESPONSIVE FACTOR8 regulates development of the feeding site induced by root-knot nematodes in tomato.

Root-knot nematodes (RKN) from the genus Meloidogyne induce the dedifferentiation of root vascular cells into giant multinucleate feeding cells. These feeding cells result from an extensive reprogramming of gene expression, and auxin is known to be a key player in their development. However, little is known about how the auxin signal is transmitted during giant cell development. Integrative analyses combining transcriptome and small non-coding RNA datasets with the specific sequencing of cleaved transcripts identified genes targeted by miRNAs in tomato (Solanum lycopersicum) galls. The two auxin-responsive transcription factors ARF8A and ARF8B, and their miRNA167 regulators, were identified as robust gene-miRNA pair candidates to be involved in the tomato response to M. incognita. Spatiotemporal expression analysis using promoter-ß-glucuronidase (GUS) fusions showed the up-regulation of ARF8A and ARF8B in RKN-induced feeding cells and surrounding cells. The generation and phenotyping of CRISPR (clustered regularly interspaced palindromic repeats) mutants demonstrated the role of ARF8A and ARF8B in giant cell development and allowed the characterization of their downstream regulated genes.

Copper microRNAs modulate the formation of giant feeding cells induced by the root knot nematode Meloidogyne incognita in Arabidopsis thaliana.

Root-knot nematodes (RKNs) are root endoparasites that induce the dedifferentiation of a few root cells and the reprogramming of their gene expression to generate giant hypermetabolic feeding cells. We identified two microRNA families, miR408 and miR398, as upregulated in Arabidopsis thaliana and Solanum lycopersicum roots infected by RKNs. In plants, the expression of these two conserved microRNA families is known to be activated by the SPL7 transcription factor in response to copper starvation. By combining functional approaches, we deciphered the network involving these microRNAs, their regulator and their targets. MIR408 expression was located within nematode-induced feeding cells like its regulator SPL7 and was regulated by copper. Moreover, infection assays with mir408 and spl7 knockout mutants or lines expressing targets rendered resistant to cleavage by miR398 demonstrated the essential role of the SPL7/MIR408/MIR398 module in the formation of giant feeding cells. Our findings reveal how perturbation of plant copper homeostasis, via the SPL7/MIR408/MIR398 module, modulates the development of nematode-induced feeding cells.

Root-associated microbiota drive phytoremediation strategies to lead of Sonchus Asper (L.) Hill as revealed by intercropping-induced modifications of the rhizosphere microbiome.

Intercropping or assistant endophytes promote phytoremediation capacities of hyperaccumulators and enhance their tolerance to heavy metal (HM) stress. Findings from a previous study showed that intercropping the hyperaccumulator Sonchus asper (L.) Hill grown in HM-contaminated soils with maize improved the remediating properties and indicated an excluder-to-hyperaccumulator switched mode of action towards lead. In the current study, RNA-Seq analysis was conducted on Sonchus roots grown under intercropping or monoculture systems to explore the molecular events underlying this shift in lead sequestering strategy. The findings showed that intercropping only slightly affects S. asper transcriptome but significantly affects expression of root-associated microbial genomes. Further, intercropping triggers significant reshaping of endophytic communities associated with a 'root-to-shoot' transition of lead sequestration and improved phytoremediation capacities of S. asper. These findings indicate that accumulator activities of a weed are partially attributed to the root-associated microbiota, and a complex network of plant-microbe-plant interactions shapes the phytoremediation potential of S. asper. Analysis showed that intercropping may significantly change the structure of root-associated communities resulting in novel remediation properties, thus providing a basis for improving phytoremediation practices to restore contaminated soils.

Chromatin Landscape Dynamics in the Early Development of the Plant Parasitic Nematode Meloidogyne incognita.

In model organisms, epigenome dynamics underlies a plethora of biological processes. The role of epigenetic modifications in development and parasitism in nematode pests remains unknown. The root-knot nematode Meloidogyne incognita adapts rapidly to unfavorable conditions, despite its asexual reproduction. However, the mechanisms underlying this remarkable plasticity and their potential impact on gene expression remain unknown. This study provides the first insight into contribution of epigenetic mechanisms to this plasticity, by studying histone modifications in M. incognita. The distribution of five histone modifications revealed the existence of strong epigenetic signatures, similar to those found in the model nematode Caenorhabditis elegans. We investigated their impact on chromatin structure and their distribution relative to transposable elements (TE) loci. We assessed the influence of the chromatin landscape on gene expression at two developmental stages: eggs, and pre-parasitic juveniles. H3K4me3 histone modification was strongly correlated with high levels of expression for protein-coding genes implicated in stage-specific processes during M. incognita development. We provided new insights in the dynamic regulation of parasitism genes kept under histone modifications silencing. In this pioneering study, we establish a comprehensive framework for the importance of epigenetic mechanisms in the regulation of the genome expression and its stability in plant-parasitic nematodes.

Movements of transposable elements contribute to the genomic plasticity and species diversification in an asexually reproducing nematode pest.

Despite reproducing without sexual recombination, Meloidogyne incognita is an adaptive and versatile phytoparasitic nematode. This species displays a global distribution, can parasitize a large range of plants, and can overcome plant resistance in a few generations. The mechanisms underlying this adaptability remain poorly known. At the whole-genome level, only a few single nucleotide variations have been observed across different geographical isolates with distinct ranges of compatible hosts. Exploring other factors possibly involved in genomic plasticity is thus important. Transposable elements (TEs), by their repetitive nature and mobility, can passively and actively impact the genome dynamics. This is particularly expected in polyploid hybrid genomes such as the one of M. incognita. Here, we have annotated the TE content of M. incognita, analyzed the statistical properties of this TE landscape, and used whole-genome pool-seq data to estimate the mobility of these TEs across twelve geographical isolates, presenting variations in ranges of compatible host plants. DNA transposons are more abundant than retrotransposons, and the high similarity of TE copies to their consensus sequences suggests they have been at least recently active. We have identified loci in the genome where the frequencies of presence of a TE showed substantial variations across the different isolates. Overall, variations in TE frequencies across isolates followed their phylogenetic divergence, suggesting TEs participate in the species diversification. Compared with the M. incognita reference genome, we detected isolate and lineage-specific de novo insertion of some TEs, including within genic regions or in the upstream regulatory regions. We validated by PCR the insertion of some of these TEs inside genic regions, confirming TE movements have possible functional impacts. Overall, we show DNA transposons can drive genomic plasticity in M. incognita and their role in genome evolution of other parthenogenetic animal deserves further investigation.

Genome Expression Dynamics Reveal the Parasitism Regulatory Landscape of the Root-Knot Nematode Meloidogyne incognita and a Promoter Motif Associated with Effector Genes.

Root-knot nematodes (genus Meloidogyne) are the major contributor to crop losses caused by nematodes. These nematodes secrete effector proteins into the plant, derived from two sets of pharyngeal gland cells, to manipulate host physiology and immunity. Successful completion of the life cycle, involving successive molts from egg to adult, covers morphologically and functionally distinct stages and will require precise control of gene expression, including effector genes. The details of how root-knot nematodes regulate transcription remain sparse. Here, we report a life stage-specific transcriptome of Meloidogyne incognita. Combined with an available annotated genome, we explore the spatio-temporal regulation of gene expression. We reveal gene expression clusters and predicted functions that accompany the major developmental transitions. Focusing on effectors, we identify a putative cis-regulatory motif associated with expression in the dorsal glands, providing an insight into effector regulation. We combine the presence of this motif with several other criteria to predict a novel set of putative dorsal gland effectors. Finally, we show this motif, and thereby its utility, is broadly conserved across the Meloidogyne genus, and we name it Mel-DOG. Taken together, we provide the first genome-wide analysis of spatio-temporal gene expression in a root-knot nematode and identify a new set of candidate effector genes that will guide future functional analyses.

The molecular signatures of compatible and incompatible pollination in Arabidopsis.

BACKGROUND: Fertilization in flowering plants depends on the early contact and acceptance of pollen grains by the receptive papilla cells of the stigma. Deciphering the specific transcriptomic response of both pollen and stigmatic cells during their interaction constitutes an important challenge to better our understanding of this cell recognition event. RESULTS: Here we describe a transcriptomic analysis based on single nucleotide polymorphisms (SNPs) present in two Arabidopsis thaliana accessions, one used as female and the other as male. This strategy allowed us to distinguish 80% of transcripts according to their parental origins. We also developed a tool which predicts male/female specific expression for genes without SNP. We report an unanticipated transcriptional activity triggered in stigma upon incompatible pollination and show that following compatible interaction, components of the pattern-triggered immunity (PTI) pathway are induced on the female side. CONCLUSIONS: Our work unveils the molecular signatures of compatible and incompatible pollinations both at the male and female side. We provide invaluable resource and tools to identify potential new molecular players involved in pollen-stigma interaction.

Author Correction: Genome assembly and annotation of Meloidogyne enterolobii, an emerging parthenogenetic root-knot nematode.

A Correction to this paper has been published: https://doi.org/10.1038/s41597-020-00747-0.

Genome assembly and annotation of Meloidogyne enterolobii, an emerging parthenogenetic root-knot nematode.

Root-knot nematodes (genus Meloidogyne) are plant parasites causing huge economic loss in the agricultural industry and affecting severely numerous developing countries. Control methods against these plant pests are sparse, the preferred one being the deployment of plant cultivars bearing resistance genes against Meloidogyne species. However, M. enterolobii is not controlled by the resistance genes deployed in the crop plants cultivated in Europe. The recent identification of this species in Europe is thus a major concern. Here, we sequenced the genome of M. enterolobii using short and long-read technologies. The genome assembly spans 240 Mbp with contig N50 size of 143 kbp, enabling high-quality annotations of 59,773 coding genes, 4,068 non-coding genes, and 10,944 transposable elements (spanning 8.7% of the genome). We validated the genome size by flow cytometry and the structure, quality and completeness by bioinformatics metrics. This ensemble of resources will fuel future projects aiming at pinpointing the genome singularities, the origin, diversity, and adaptive potential of this emerging plant pest.

Characterization of Two Satellite DNA Families in the Genome of the Oomycete Plant Pathogen Phytophthora parasitica.

Satellite DNA is a class of repetitive sequences that are organized in long arrays of tandemly repeated units in most eukaryotes. Long considered as selfish DNA, satellite sequences are now proposed to contribute to genome integrity. Despite their potential impact on the architecture and evolution of the genome, satellite DNAs have not been investigated in oomycetes due to the paucity of genomic data and the difficulty of assembling highly conserved satellite arrays. Yet gaining knowledge on the structure and evolution of genomes of oomycete pathogens is crucial to understanding the mechanisms underlying adaptation to their environment and to proposing efficient disease control strategies. A de novo assembly of the genome of Phytophthora parasitica, an important oomycete plant pathogen, led to the identification of several families of tandemly repeated sequences varying in size, copy number, and sequence conservation. Among them, two abundant families, designated as PpSat1 and PpSat2, displayed typical features of satellite DNA and were collectively designated as PpSat. These two satellite families differ by their length, sequence, organization, genomic environment, and evolutionary dynamics. PpSat1, but not PpSat2, presented homologs among oomycetes. This observation, as well as the characterization of transcripts of PpSat families, suggested that these satellite DNA families likely play a conserved role within this important group of pathogens.

Gene copy number variations as signatures of adaptive evolution in the parthenogenetic, plant-parasitic nematode Meloidogyne incognita.

Adaptation to changing environmental conditions represents a challenge to parthenogenetic organisms, and until now, how phenotypic variants are generated in clones in response to the selection pressure of their environment remains poorly known. The obligatory parthenogenetic root-knot nematode species Meloidogyne incognita has a worldwide distribution and is the most devastating plant-parasitic nematode. Despite its asexual reproduction, this species exhibits an unexpected capacity of adaptation to environmental constraints, for example, resistant hosts. Here, we used a genomewide comparative hybridization strategy to evaluate variations in gene copy numbers between genotypes of M. incognita resulting from two parallel experimental evolution assays on a susceptible vs. resistant host plant. We detected gene copy number variations (CNVs) associated with the ability of the nematodes to overcome resistance of the host plant, and this genetic variation may reflect an adaptive response to host resistance in this parthenogenetic species. The CNV distribution throughout the nematode genome is not random and suggests the occurrence of genomic regions more prone to undergo duplications and losses in response to the selection pressure of the host resistance. Furthermore, our analysis revealed an outstanding level of gene loss events in nematode genotypes that have overcome the resistance. Overall, our results support the view that gene loss could be a common class of adaptive genetic mechanism in response to a challenging new biotic environment in clonal animals.

Characterization of siRNAs clusters in Arabidopsis thaliana galls induced by the root-knot nematode Meloidogyne incognita.

BACKGROUND: Root-knot nematodes (RKN), genus Meloidogyne, are plant parasitic worms that have the ability to transform root vascular cylinder cells into hypertrophied, multinucleate and metabolically over-active feeding cells. Redifferentiation into feeding cells is the result of a massive transcriptional reprogramming of root cells targeted by RKN. Since RKN are able to induce similar feeding cells in roots of thousands of plant species, these worms are thought to manipulate essential and conserved plant molecular pathways. RESULTS: Small non-coding RNAs of uninfected roots and infected root galls induced by M. incognita from Arabidopsis thaliana were sequenced by high throughput sequencing. SiRNA populations were analysed by using the Shortstack algorithm. We identified siRNA clusters that are differentially expressed in infected roots and evidenced an over-representation of the 23-24 nt siRNAs in infected tissue. This size corresponds to heterochromatic siRNAs (hc-siRNAs) which are known to regulate expression of transposons and genes at the transcriptional level, mainly by inducing DNA methylation. CONCLUSIONS: Correlation of siRNA clusters expression profile with transcriptomic data identified several protein coding genes that are candidates to be regulated by siRNAs at the transcriptional level by RNA directed DNA methylation (RdDM) pathway either directly or indirectly via silencing of neighbouring transposable elements.

Genome-wide expert annotation of the epigenetic machinery of the plant-parasitic nematodes Meloidogyne spp., with a focus on the asexually reproducing species.

BACKGROUND: The renewed interest in epigenetics has led to the understanding that both the environment and individual lifestyle can directly interact with the epigenome to influence its dynamics. Epigenetic phenomena are mediated by DNA methylation, stable chromatin modifications and non-coding RNA-associated gene silencing involving specific proteins called epigenetic factors. Multiple organisms, ranging from plants to yeast and mammals, have been used as model systems to study epigenetics. The interactions between parasites and their hosts are models of choice to study these mechanisms because the selective pressures are strong and the evolution is fast. The asexually reproducing root-knot nematodes (RKN) offer different advantages to study the processes and mechanisms involved in epigenetic regulation. RKN genomes sequencing and annotation have identified numerous genes, however, which of those are involved in the adaption to an environment and potentially relevant to the evolution of plant-parasitism is yet to be discovered. RESULTS: Here, we used a functional comparative annotation strategy combining orthology data, mining of curated genomics as well as protein domain databases and phylogenetic reconstructions. Overall, we show that (i) neither RKN, nor the model nematode Caenorhabditis elegans possess any DNA methyltransferases (DNMT) (ii) RKN do not possess the complete machinery for DNA methylation on the 6th position of adenine (6mA) (iii) histone (de)acetylation and (de)methylation pathways are conserved between C. elegans and RKN, and the corresponding genes are amplified in asexually reproducing RKN (iv) some specific non-coding RNA families found in plant-parasitic nematodes are dissimilar from those in C. elegans. In the asexually reproducing RKN Meloidogyne incognita, expression data from various developmental stages supported the putative role of these proteins in epigenetic regulations. CONCLUSIONS: Our results refine previous predictions on the epigenetic machinery of model species and constitute the most comprehensive description of epigenetic factors relevant to the plant-parasitic lifestyle and/or asexual mode of reproduction of RKN. Providing an atlas of epigenetic factors in RKN is an informative resource that will enable researchers to explore their potential role in adaptation of these parasites to their environment.

A root-knot nematode small glycine and cysteine-rich secreted effector, MiSGCR1, is involved in plant parasitism.

Root-knot nematodes, Meloidogyne spp., are obligate endoparasites that maintain a biotrophic relationship with their hosts. They infect roots as microscopic vermiform second-stage juveniles, and establish specialized feeding structures called 'giant-cells', from which they withdraw water and nutrients. The nematode effector proteins secreted in planta are key elements in the molecular dialogue of parasitism. Here, we compared Illumina RNA-seq transcriptomes for M. incognita obtained at various points in the lifecycle, and identified 31 genes more strongly expressed in parasitic stages than in preparasitic juveniles. We then selected candidate effectors for functional characterization. Quantitative real-time PCR and in situ hybridizations showed that the validated differentially expressed genes are predominantly specifically expressed in oesophageal glands of the nematode. We also soaked the nematodes in siRNA to silence these genes and to determine their role in pathogenicity. The silencing of the dorsal gland specific-Minc18876 and its paralogues resulted in a significant, reproducible decrease in the number of mature females with egg masses, demonstrating a potentially important role for the small glycine- and cysteine-rich effector MiSGCR1 in early stages of plant-nematode interaction. Finally, we report that MiSGCR1 suppresses plant cell death induced by bacterial or oomycete triggers of plant defense.

Transcriptome Profiling of Neurosensory Perception Genes in Wing Tissue of Two Evolutionary Distant Insect Orders: Diptera (Drosophila melanogaster) and Hemiptera (Acyrthosiphon pisum).

The neurogenesis and neuronal functions in insect wing have been understudied mainly due to technical hindrances that have prevented electrophysiology studies for decades. The reason is that the nano-architecture of the wing chemosensory bristles hampers the receptors accessibility of odorants/tastants to receptors in fixed setup, whereas in nature, the wing flapping mixes these molecules in bristle lymph. In this report, we analyzed the transcriptome of the wing tissue of two species phylogenetically strongly divergent: Drosophila melanogaster a generic model for diptera order (complete metamorphosis) and the aphid acyrthosiphon pisum, representative of hemiptera order (incomplete metamorphosis) for which a conditional winged/wingless polyphenism is under control of population density and resources. The transcriptome shows that extensive gene networks involved in chemosensory perception are active in adult wing for both species. Surprisingly, the specific transcripts of genes that are commonly found in eye were present in Drosophila wing but not in aphid. The analysis reveals that in the aphid conditional wing, expressed genes show strong similarities with those in the gut epithelia. This suggests that the epithelial cell layer between the cuticle sheets is persistent at least in young aphid adult, whereas it disappears after emergence in Drosophila. Despite marked differences between the two transcriptomes, the results highlight the probable universalism of wing chemosensory function in the holometabolous and hemimetabolous orders of winged insects.

The Transcriptomes of Xiphinema index and Longidorus elongatus Suggest Independent Acquisition of Some Plant Parasitism Genes by Horizontal Gene Transfer in Early-Branching Nematodes.

Nematodes have evolved the ability to parasitize plants on at least four independent occasions, with plant parasites present in Clades 1, 2, 10 and 12 of the phylum. In the case of Clades 10 and 12, horizontal gene transfer of plant cell wall degrading enzymes from bacteria and fungi has been implicated in the evolution of plant parasitism. We have used ribonucleic acid sequencing (RNAseq) to generate reference transcriptomes for two economically important nematode species, Xiphinema index and Longidorus elongatus, representative of two genera within the early-branching Clade 2 of the phylum Nematoda. We used a transcriptome-wide analysis to identify putative horizontal gene transfer events. This represents the first in-depth transcriptome analysis from any plant-parasitic nematode of this clade. For each species, we assembled ~30 million Illumina reads into a reference transcriptome. We identified 62 and 104 transcripts, from X. index and L. elongatus, respectively, that were putatively acquired via horizontal gene transfer. By cross-referencing horizontal gene transfer prediction with a phylum-wide analysis of Pfam domains, we identified Clade 2-specific events. Of these, a GH12 cellulase from X. index was analysed phylogenetically and biochemically, revealing a likely bacterial origin and canonical enzymatic function. Horizontal gene transfer was previously shown to be a phenomenon that has contributed to the evolution of plant parasitism among nematodes. Our findings underline the importance and the extensiveness of this phenomenon in the evolution of plant-parasitic life styles in this speciose and widespread animal phylum.

Characterization of microRNAs from Arabidopsis galls highlights a role for miR159 in the plant response to the root-knot nematode Meloidogyne incognita.

Root knot nematodes (RKN) are root parasites that induce the genetic reprogramming of vascular cells into giant feeding cells and the development of root galls. MicroRNAs (miRNAs) regulate gene expression during development and plant responses to various stresses. Disruption of post-transcriptional gene silencing in Arabidopsis ago1 or ago2 mutants decrease the infection rate of RKN suggesting a role for this mechanism in the plant-nematode interaction. By sequencing small RNAs from uninfected Arabidopsis roots and from galls 7 and 14 d post infection with Meloidogyne incognita, we identified 24 miRNAs differentially expressed in gall as putative regulators of gall development. Moreover, strong activity within galls was detected for five miRNA promoters. Analyses of nematode development in an Arabidopsis miR159abc mutant had a lower susceptibility to RKN, suggesting a role for the miR159 family in the plant response to M. incognita. Localization of mature miR159 within the giant and surrounding cells suggested a role in giant cell and gall. Finally, overexpression of miR159 in galls at 14 d post inoculation was associated with the repression of the miR159 target MYB33 which expression is restricted to the early stages of infection. Overall, these results implicate the miR159 in plant responses to RKN.

Hybridization and polyploidy enable genomic plasticity without sex in the most devastating plant-parasitic nematodes.

Root-knot nematodes (genus Meloidogyne) exhibit a diversity of reproductive modes ranging from obligatory sexual to fully asexual reproduction. Intriguingly, the most widespread and devastating species to global agriculture are those that reproduce asexually, without meiosis. To disentangle this surprising parasitic success despite the absence of sex and genetic exchanges, we have sequenced and assembled the genomes of three obligatory ameiotic and asexual Meloidogyne. We have compared them to those of relatives able to perform meiosis and sexual reproduction. We show that the genomes of ameiotic asexual Meloidogyne are large, polyploid and made of duplicated regions with a high within-species average nucleotide divergence of ~8%. Phylogenomic analysis of the genes present in these duplicated regions suggests that they originated from multiple hybridization events and are thus homoeologs. We found that up to 22% of homoeologous gene pairs were under positive selection and these genes covered a wide spectrum of predicted functional categories. To biologically assess functional divergence, we compared expression patterns of homoeologous gene pairs across developmental life stages using an RNAseq approach in the most economically important asexually-reproducing nematode. We showed that >60% of homoeologous gene pairs display diverged expression patterns. These results suggest a substantial functional impact of the genome structure. Contrasting with high within-species nuclear genome divergence, mitochondrial genome divergence between the three ameiotic asexuals was very low, signifying that these putative hybrids share a recent common maternal ancestor. Transposable elements (TE) cover a ~1.7 times higher proportion of the genomes of the ameiotic asexual Meloidogyne compared to the sexual relative and might also participate in their plasticity. The intriguing parasitic success of asexually-reproducing Meloidogyne species could be partly explained by their TE-rich composite genomes, resulting from allopolyploidization events, and promoting plasticity and functional divergence between gene copies in the absence of sex and meiosis.

Nod Factor Effects on Root Hair-Specific Transcriptome of Medicago truncatula: Focus on Plasma Membrane Transport Systems and Reactive Oxygen Species Networks.

Root hairs are involved in water and nutrient uptake, and thereby in plant autotrophy. In legumes, they also play a crucial role in establishment of rhizobial symbiosis. To obtain a holistic view of Medicago truncatula genes expressed in root hairs and of their regulation during the first hours of the engagement in rhizobial symbiotic interaction, a high throughput RNA sequencing on isolated root hairs from roots challenged or not with lipochitooligosaccharides Nod factors (NF) for 4 or 20 h was carried out. This provided a repertoire of genes displaying expression in root hairs, responding or not to NF, and specific or not to legumes. In analyzing the transcriptome dataset, special attention was paid to pumps, transporters, or channels active at the plasma membrane, to other proteins likely to play a role in nutrient ion uptake, NF electrical and calcium signaling, control of the redox status or the dynamic reprogramming of root hair transcriptome induced by NF treatment, and to the identification of papilionoid legume-specific genes expressed in root hairs. About 10% of the root hair expressed genes were significantly up- or down-regulated by NF treatment, suggesting their involvement in remodeling plant functions to allow establishment of the symbiotic relationship. For instance, NF-induced changes in expression of genes encoding plasma membrane transport systems or disease response proteins indicate that root hairs reduce their involvement in nutrient ion absorption and adapt their immune system in order to engage in the symbiotic interaction. It also appears that the redox status of root hair cells is tuned in response to NF perception. In addition, 1176 genes that could be considered as "papilionoid legume-specific" were identified in the M. truncatula root hair transcriptome, from which 141 were found to possess an ortholog in every of the six legume genomes that we considered, suggesting their involvement in essential functions specific to legumes. This transcriptome provides a valuable resource to investigate root hair biology in legumes and the roles that these cells play in rhizobial symbiosis establishment. These results could also contribute to the long-term objective of transferring this symbiotic capacity to non-legume plants.