Integrated multi-omics and genetic analyses reveal molecular determinants underlying Arabidopsis snap33 mutant phenotype.

The secretory pathway is essential for plant immunity, delivering diverse antimicrobial molecules into the extracellular space. Arabidopsis thaliana soluble N-ethylmaleimide-sensitive-factor attachment protein receptor SNAP33 is a key actor of this process. The snap33 mutant displays dwarfism and necrotic lesions, however the molecular determinants of its macroscopic phenotypes remain elusive. Here, we isolated several new snap33 mutants that exhibited constitutive cell death and H2O2 accumulation, further defining snap33 as an autoimmune mutant. We then carried out quantitative transcriptomic and proteomic analyses showing that numerous defense transcripts and proteins were up-regulated in the snap33 mutant, among which genes/proteins involved in defense hormone, pattern-triggered immunity, and nucleotide-binding domain leucine-rich-repeat receptor signaling. qRT-PCR analyses and hormone dosages supported these results. Furthermore, genetic analyses elucidated the diverse contributions of the main defense hormones and some nucleotide-binding domain leucine-rich-repeat receptor signaling actors in the establishment of the snap33 phenotype, emphasizing the preponderant role of salicylic acid over other defense phytohormones. Moreover, the accumulation of pattern-triggered immunity and nucleotide-binding domain leucine-rich-repeat receptor signaling proteins in the snap33 mutant was confirmed by immunoblotting analyses and further shown to be salicylic acid-dependent. Collectively, this study unveiled molecular determinants underlying the Arabidopsis snap33 mutant phenotype and brought new insights into autoimmunity signaling.

Impact of high atmospheric carbon dioxide on the biotic stress response of the model cereal species Brachypodium distachyon.

Losses due to disease and climate change are among the most important issues currently facing crop production. It is therefore important to establish the impact of climate change, and particularly of high carbon dioxide (hCO2), on plant immunity in cereals, which provide 60% of human calories. The aim of this study was to determine if hCO2 impacts Brachypodium distachyon immunity, a model plant for temperate cereals. Plants were grown in air (430 ppm CO2) and at two high CO2 conditions, one that is relevant to projections within the coming century (1000 ppm) and a concentration sufficient to saturate photosynthesis (3000 ppm). The following measurements were performed: phenotyping and growth, salicylic acid contents, pathogen resistance tests, and RNAseq analysis of the transcriptome. Improved shoot development was observed at both 1000 and 3000 ppm. A transcriptomic analysis pointed to an increase in primary metabolism capacity under hCO2. Alongside this effect, up-regulation of genes associated with secondary metabolism was also observed. This effect was especially evident for the terpenoid and phenylpropanoid pathways, and was accompanied by enhanced expression of immunity-related genes and accumulation of salicylic acid. Pathogen tests using the fungus Magnaporthe oryzae revealed that hCO2 had a complex effect, with enhanced susceptibility to infection but no increase in fungal development. The study reveals that immunity in B. distachyon is modulated by growth at hCO2 and allows identification of pathways that might play a role in this effect.

Conditional generation of free radicals by selective activation of alkoxyamines: towards more effective and less toxic targeting of brain tumors.

Brain tumors are an important cause of suffering and death. Glioblastoma are the most frequent primary tumors of the central nervous system in adults. They are associated with a very poor prognosis, since only 10% of GBM patients survive 5 years after diagnosis. Medulloblastoma are the most frequent brain malignancies in childhood; they affect the cerebellum in children under 10 years of age in 75% of cases. The current multimodal treatment comes at the expense of serious and often long-lasting side effects. Herein, we propose the synthesis of a library of novel alkoxyamines as anticancer drug candidates. The most efficient molecule, ALK4, was selected based on its ability to inhibit both survival and migration of GBM and MB cells in 2D cultures and in 3D tumor spheroids. A fluorescent derivative was used to show the early cytosolic accumulation of ALK4 in tumor cells. Spontaneous homolysis of ALK4 led to the release of alkyl radicals, which triggered the generation of reactive oxygen species, fragmentation of the mitochondrial network and ultimately apoptosis. To control its homolytic process, the selected alkoxyamine was bioconjugated to a peptide selectively recognized by matrix metalloproteases. This bioconjugate, named ALK4-MMPp, successfully inhibited survival, proliferation, and invasion of GBM and MB tumor micromasses. We further developed innovative brain and cerebellum organotypic models to monitor treatment response over time. It confirmed that ALK4-MMPp significantly impaired tumor progression, while no significant damage was observed on normal brain tissue. Lastly, we showed that ALK4-MMPp was well-tolerated in vivo by zebrafish embryos. This study provides a new strategy to control the activation of alkoxyamines, and revealed the bioconjugate ALK4-MMPp bioconjugate as a good anticancer drug candidate.

The Arabidopsis transcription factor NLP2 regulates early nitrate responses and integrates nitrate assimilation with energy and carbon skeleton supply.

Nitrate signaling improves plant growth under limited nitrate availability and, hence, optimal resource use for crop production. Whereas several transcriptional regulators of nitrate signaling have been identified, including the Arabidopsis thaliana transcription factor NIN-LIKE PROTEIN7 (NLP7), additional regulators are expected to fine-tune this pivotal physiological response. Here, we characterized Arabidopsis NLP2 as a top-tier transcriptional regulator of the early nitrate response gene regulatory network. NLP2 interacts with NLP7 in vivo and shares key molecular features such as nitrate-dependent nuclear localization, DNA-binding motif, and some target genes with NLP7. Genetic, genomic, and metabolic approaches revealed a specific role for NLP2 in the nitrate-dependent regulation of carbon and energy-related processes that likely influence plant growth under distinct nitrogen environments. Our findings highlight the complementarity and specificity of NLP2 and NLP7 in orchestrating a multitiered nitrate regulatory network that links nitrate assimilation with carbon and energy metabolism for efficient nitrogen use and biomass production.

Drought specifically downregulates mineral nutrition: Plant ionomic content and associated gene expression.

One of the main limiting factors of plant yield is drought, and while the physiological responses to this environmental stress have been broadly described, research addressing its impact on mineral nutrition is scarce. Brassica napus and Triticum aestivum were subjected to moderate or severe water deficit, and their responses to drought were assessed by functional ionomic analysis, and derived calculation of the net uptake of 20 nutrients. While the uptake of most mineral nutrients decreased, Fe, Zn, Mn, and Mo uptake were impacted earlier and at a larger scale than most physiological parameters assessed (growth, ABA concentration, gas exchanges and photosynthetic activity). Additionally, in B. napus, the patterns of 183 differentially expressed genes in leaves related to the ionome (known ionomic genes, KIGs) or assumed to be involved in transport of a given nutrient were analyzed. This revealed three patterns of gene expression under drought consisting of up (transport of Cl and Co), down (transport of N, P, B, Mo, and Ni), or mixed levels (transport of S, Mg, K, Zn, Fe, Cu, or Mn) of regulation. The three patterns of gene regulations are discussed in relation to specific gene functions, changes of leaf ionomic composition and with consideration of the crosstalks that have been established between elements. It is suggested that the observed reduction in Fe uptake occurred via a specific response to drought, leading indirectly to reduced uptake of Zn and Mn, and these may be taken up by common transporters encoded by genes that were downregulated.

Analyzing Multifactorial RNA-Seq Experiments with DicoExpress.

The proper use of statistical modeling in NGS data analysis requires an advanced level of expertise. There has recently been a growing consensus on using generalized linear models for differential analysis of RNA-Seq data and the advantage of mixture models to perform co-expression analysis. To offer a managed setting to use these modeling approaches, we developed DiCoExpress that provides a standardized R pipeline to perform an RNA-Seq analysis. Without any particular knowledge in statistics or R programming, beginners can perform a complete RNA-Seq analysis from quality controls to co-expression through differential analysis based on contrasts inside a generalized linear model. An enrichment analysis is proposed both on the lists of differentially expressed genes, and the co-expressed gene clusters. This video tutorial is conceived as a step-by-step protocol to help users take full advantage of DiCoExpress and its potential in empowering the biological interpretation of an RNA-Seq experiment.

Genome-Wide Transcriptomic Analysis of the Effects of Infection with the Hemibiotrophic Fungus Colletotrichum lindemuthianum on Common Bean.

Bean anthracnose caused by the hemibiotrophic fungus Colletotrichum lindemuthianum is one of the most important diseases of common bean (Phaseolus vulgaris) in the world. In the present study, the whole transcriptome of common bean infected with C. lindemuthianum during compatible and incompatible interactions was characterized at 48 and 72 hpi, corresponding to the biotrophy phase of the infection cycle. Our results highlight the prominent role of pathogenesis-related (PR) genes from the PR10/Bet vI family as well as a complex interplay of different plant hormone pathways including Ethylene, Salicylic acid (SA) and Jasmonic acid pathways. Gene Ontology enrichment analysis reveals that infected common bean seedlings responded by down-regulation of photosynthesis, ubiquitination-mediated proteolysis and cell wall modifications. In infected common bean, SA biosynthesis seems to be based on the PAL pathway instead of the ICS pathway, contrarily to what is described in Arabidopsis. Interestingly, ~30 NLR were up-regulated in both contexts. Overall, our results suggest that the difference between the compatible and incompatible reaction is more a question of timing and strength, than a massive difference in differentially expressed genes between these two contexts. Finally, we used RT-qPCR to validate the expression patterns of several genes, and the results showed an excellent agreement with deep sequencing.

Corrigendum: Genotypic Variation of Nitrogen Use Efficiency and Amino Acid Metabolism in Barley.

[This corrects the article DOI: 10.3389/fpls.2021.807798.].

Comparative Omics Analysis of Brassica napus Roots Subjected to Six Individual Macronutrient Deprivations Reveals Deficiency-Specific Genes and Metabolomic Profiles.

The early and specific diagnosis of a macronutrient deficiency is challenging when seeking to better manage fertilizer inputs in the context of sustainable agriculture. Consequently, this study explored the potential for transcriptomic and metabolomic analysis of Brassica napus roots to characterize the effects of six individual macronutrient deprivations (N, Mg, P, S, K, and Ca). Our results showed that before any visual phenotypic response, all macronutrient deprivations led to a large modulation of the transcriptome and metabolome involved in various metabolic pathways, and some were common to all macronutrient deprivations. Significantly, comparative transcriptomic analysis allowed the definition of a subset of 3282, 2011, 6325, 1384, 439, and 5157 differentially expressed genes (DEGs) specific to N, Mg, P, S, K, and Ca deprivations, respectively. Surprisingly, gene ontology term enrichment analysis performed on this subset of specific DEGs highlighted biological processes that are common to a number of these macronutrient deprivations, illustrating the complexity of nutrient interactions. In addition, a set of 38 biochemical compounds that discriminated the macronutrient deprivations was identified using a metabolic approach. The opportunity to use these specific DEGs and/or biochemical compounds as potential molecular indicators to diagnose macronutrient deficiency is discussed.

Genome-Wide Identification of Key Components of RNA Silencing in Two Phaseolus vulgaris Genotypes of Contrasting Origin and Their Expression Analyses in Response to Fungal Infection.

RNA silencing serves key roles in a multitude of cellular processes, including development, stress responses, metabolism, and maintenance of genome integrity. Dicer, Argonaute (AGO), double-stranded RNA binding (DRB) proteins, RNA-dependent RNA polymerase (RDR), and DNA-dependent RNA polymerases known as Pol IV and Pol V form core components to trigger RNA silencing. Common bean (Phaseolus vulgaris) is an important staple crop worldwide. In this study, we aimed to unravel the components of the RNA-guided silencing pathway in this non-model plant, taking advantage of the availability of two genome assemblies of Andean and Meso-American origin. We identified six PvDCLs, thirteen PvAGOs, 10 PvDRBs, 5 PvRDRs, in both genotypes, suggesting no recent gene amplification or deletion after the gene pool separation. In addition, we identified one PvNRPD1 and one PvNRPE1 encoding the largest subunits of Pol IV and Pol V, respectively. These genes were categorized into subgroups based on phylogenetic analyses. Comprehensive analyses of gene structure, genomic localization, and similarity among these genes were performed. Their expression patterns were investigated by means of expression models in different organs using online data and quantitative RT-PCR after pathogen infection. Several of the candidate genes were up-regulated after infection with the fungus Colletotrichum lindemuthianum.

Genotypic Variation of Nitrogen Use Efficiency and Amino Acid Metabolism in Barley.

Owing to the large genetic diversity of barley and its resilience under harsh environments, this crop is of great value for agroecological transition and the need for reduction of nitrogen (N) fertilizers inputs. In the present work, we investigated the diversity of a North African barley genotype collection in terms of growth under limiting N (LN) or ample N (HN) supply and in terms of physiological traits including amino acid content in young seedlings. We identified a Moroccan variety, Laanaceur, accumulating five times more lysine in its leaves than the others under both N nutritional regimes. Physiological characterization of the barley collection showed the genetic diversity of barley adaptation strategies to LN and highlighted a genotype x environment interaction. In all genotypes, N limitation resulted in global biomass reduction, an increase in C concentration, and a higher resource allocation to the roots, indicating that this organ undergoes important adaptive metabolic activity. The most important diversity concerned leaf nitrogen use efficiency (LNUE), root nitrogen use efficiency (RNUE), root nitrogen uptake efficiency (RNUpE), and leaf nitrogen uptake efficiency (LNUpE). Using LNUE as a target trait reflecting barley capacity to deal with N limitation, this trait was positively correlated with plant nitrogen uptake efficiency (PNUpE) and RNUpE. Based on the LNUE trait, we determined three classes showing high, moderate, or low tolerance to N limitation. The transcriptomic approach showed that signaling, ionic transport, immunity, and stress response were the major functions affected by N supply. A candidate gene encoding the HvNRT2.10 transporter was commonly up-regulated under LN in the three barley genotypes investigated. Genes encoding key enzymes required for lysine biosynthesis in plants, dihydrodipicolinate synthase (DHPS) and the catabolic enzyme, the bifunctional Lys-ketoglutarate reductase/saccharopine dehydrogenase are up-regulated in Laanaceur and likely account for a hyperaccumulation of lysine in this genotype. Our work provides key physiological markers of North African barley response to low N availability in the early developmental stages.

Retrograde signalling from the mitochondria to the nucleus translates the positive effect of ethylene on dormancy breaking of Arabidopsis thaliana seeds.

Ethylene and reactive oxygen species (ROS) regulate seed dormancy alleviation, but the molecular basis of their action and crosstalk remains largely unknown. Here we studied the mechanism of Arabidopsis seed dormancy release by ethylene using cell imaging, and genetic and transcriptomics approaches, in order to tackle its possible interaction with ROS homeostasis. We found that the effect of ethylene on seed germination required ROS production by the mitochondrial electron transport chain. Seed response to ethylene involved a mitochondrial retrograde response (MRR) through nuclear ROS production and upregulation of the MRR components AOX1a and ANAC013, but also required the activation of the ethylene canonical pathway. Together our data allowed deciphering of the mode of action of ethylene on seed germination and the associated dynamics of ROS production. Our findings highlight the occurrence of retrograde signalling in seed germination.

Gene expression predictions and networks in natural populations supports the omnigenic theory.

BACKGROUND: Recent literature on the differential role of genes within networks distinguishes core from peripheral genes. If previous works have shown contrasting features between them, whether such categorization matters for phenotype prediction remains to be studied. RESULTS: We measured 17 phenotypic traits for 241 cloned genotypes from a Populus nigra collection, covering growth, phenology, chemical and physical properties. We also sequenced RNA for each genotype and built co-expression networks to define core and peripheral genes. We found that cores were more differentiated between populations than peripherals while being less variable, suggesting that they have been constrained through potentially divergent selection. We also showed that while cores were overrepresented in a subset of genes statistically selected for their capacity to predict the phenotypes (by Boruta algorithm), they did not systematically predict better than peripherals or even random genes. CONCLUSION: Our work is the first attempt to assess the importance of co-expression network connectivity in phenotype prediction. While highly connected core genes appear to be important, they do not bear enough information to systematically predict better quantitative traits than other gene sets.

DiCoExpress: a tool to process multifactorial RNAseq experiments from quality controls to co-expression analysis through differential analysis based on contrasts inside GLM models.

BACKGROUND: RNAseq is nowadays the method of choice for transcriptome analysis. In the last decades, a high number of statistical methods, and associated bioinformatics tools, for RNAseq analysis were developed. More recently, statistical studies realised neutral comparison studies using benchmark datasets, shedding light on the most appropriate approaches for RNAseq data analysis. RESULTS: DiCoExpress is a script-based tool implemented in R that includes methods chosen based on their performance in neutral comparisons studies. DiCoExpress uses pre-existing R packages including FactoMineR, edgeR and coseq, to perform quality control, differential, and co-expression analysis of RNAseq data. Users can perform the full analysis, providing a mapped read expression data file and a file containing the information on the experimental design. Following the quality control step, the user can move on to the differential expression analysis performed using generalized linear models thanks to the automated contrast writing function. A co-expression analysis is implemented using the coseq package. Lists of differentially expressed genes and identified co-expression clusters are automatically analyzed for enrichment of annotations provided by the user. We used DiCoExpress to analyze a publicly available RNAseq dataset on the transcriptional response of Brassica napus L. to silicon treatment in plant roots and mature leaves. This dataset, including two biological factors and three replicates for each condition, allowed us to demonstrate in a tutorial all the features of DiCoExpress. CONCLUSIONS: DiCoExpress is an R script-based tool allowing users to perform a full RNAseq analysis from quality controls to co-expression analysis through differential analysis based on contrasts inside generalized linear models. DiCoExpress focuses on the statistical modelling of gene expression according to the experimental design and facilitates the data analysis leading the biological interpretation of the results.

Mutations of the AtYAK1 Kinase Suppress TOR Deficiency in Arabidopsis.

The target of rapamycin (TOR) kinase is a conserved energy sensor that regulates growth in response to environmental cues. However, little is known about the TOR signaling pathway in plants. We used Arabidopsis lines affected in the lethal with SEC13 protein 8 (LST8-1) gene, a core element of the TOR complex, to search for suppressor mutations. Two suppressor lines with improved growth were isolated that carried mutations in the Yet Another Kinase 1 (AtYAK1) gene encoding a member of the dual-specificity tyrosine phosphorylation-regulated kinase (DYRK) family. Atyak1 mutations partly rescued the developmental defects of lst8-1-1 mutants and conferred resistance to the TOR inhibitor AZD-8055. Moreover, atyak1 mutations suppressed the transcriptomic and metabolic perturbations as well as the abscisic acid (ABA) hypersensitivity of the lst8-1-1 mutants. AtYAK1 interacted with the regulatory-associated protein of TOR (RAPTOR), a component of the TOR complex, and was phosphorylated by TOR. Thus, our findings reveal that AtYAK1 is a TOR effector that probably needs to be switched off to activate plant growth.

Genetic diversity of rhizobia associated with root nodules of white lupin (Lupinus albus L.) in Tunisian calcareous soils.

With a view to introducing white lupin (Lupinus albus L.) for cultivation in Tunisian calcareous soils, compatible indigenous rhizobia for nitrogen-fixing symbiosis were investigated and characterized. Two L. albus varieties, Mekna and Lumen, were used to trap rhizobia in soil samples collected from 56 sites with high active lime contents (0-49%). Nodulation occurred in only 15 soils. The local variety, Mekna, developed significantly more root nodules and had a trapping capacity in more soils than the imported variety Lumen. A phylogenetic analysis based on the partial 16S-23S ribosomal RNA internal transcribed spacer region (ITS) and multi-locus sequence analysis (MLSA) of three chromosomal housekeeping genes, recA, atpD and dnaK, showed that strains were affiliated to Agrobacterium, Rhizobium, and Neorhizobium, with large internal diversity, including separate lineages. Infectivity tests highlighted some nodulation specificity at the plant variety level, since the strains originating from Mekna could only nodulate this variety, while strains trapped in Lumen could nodulate both varieties. When inoculated, almost all strains resulted in a significant increase in plant shoot dry weight on L. albus. Although Agrobacterium sp. strains isolated from L. albus could nodulate and had a plant growth promoting effect, no nodA and nodC genes could be amplified. This is discussed together with the absence of bradyrhizobia and the general infrequency of L. albus-nodulating rhizobia in Tunisian soils. The adapted and efficient rhizobial strains reported here were promising candidates for inoculant development and represent a contribution towards successful cultivation of L. albus in Tunisia, especially the most promising Mekna variety.

The rhizosphere of the halophytic grass Sporobolus robustus Kunth hosts rhizobium genospecies that are efficient on Prosopis juliflora (Sw.) DC and Vachellia seyal (Del.) P.J.H. Hurter seedlings.

The aim of this study was to survey the abundance and genetic diversity of legume-nodulating rhizobia (LNR) in the rhizosphere of a salt-tolerant grass, Sporobolus robustus Kunth, in the dry and rainy seasons along a salinity gradient, and to test their effectiveness on Prosopis juliflora (SW.) DC and Vachellia seyal (Del.) P.J.H. Hurter seedlings. The results showed a significant decrease in LNR population density and diversity in response to salinity, particularly during the dry season. A phylogenetic analysis of the 16S-23S rRNA ITS region clustered the 232 rhizobium isolates into three genera and 12 distinct representative genotypes: Mesorhizobium (8 genotypes), Ensifer (2 genotypes) and Rhizobium (2 genotypes). Of these genotypes, 2 were only found in the dry season, 4 exclusively in the rainy season and 6 were found in both seasons. Isolates of the Mesorhizobium and Ensifer genera were more abundant than those of Rhizobium, with 55%, 44% and 1% of the total strains, respectively. The abundance of the Mesorhizobium isolates appeared to increase in the dry season, suggesting that they were more adapted to environmental aridity than Ensifer genospecies. Conversely, Ensifer genospecies were more tolerant of high salinity levels than the other genospecies. However, Ensifer genospeciesproved to be the most efficient strains on P. juliflora and V. seyal seedlings. We concluded that S. robustus hosts efficient rhizobium strains in its rhizosphere, suggesting its ability to act as a nurse plant to facilitate seedling recruitment of P. juliflora and V. seyal in saline soils.

Phylogenetic diversity of Bradyrhizobium strains nodulating Calicotome spinosa in the Northeast of Algeria.

Fifty-two slow-growing strains were isolated from root nodules of Calicotome spinosa grown in the Northeast of Algeria and grouped in 24 rep-PCR clusters. One representative strain for each profile was further phylogenetically characterized. The nearly complete 16S rRNA gene sequence indicated that all strains were affiliated to Bradyrhizobium. Multi-Locus Sequence Analysis (MLSA) of the atpD, glnII and recA genes and of the 16S-23S rRNA internal transcribed spacer (ITS) showed that these strains formed four divergent clusters: one close to Bradyrhizobium canariense and Bradyrhizobium lupini and three others separate from all the described species, representing three putative new Bradyrhizobium species. A phylogenetic analysis based on the nodC gene sequence affiliated the strains to either of the two symbiovars, genistearum or retamae.

Combining laser-assisted microdissection (LAM) and RNA-seq allows to perform a comprehensive transcriptomic analysis of epidermal cells of Arabidopsis embryo.

BACKGROUND: Genome-wide characterization of tissue- or cell-specific gene expression is a recurrent bottleneck in biology. We have developed a sensitive approach based on ultra-low RNA sequencing coupled to laser assisted microdissection for analyzing different tissues of the small Arabidopsis embryo. METHODS AND RESULTS: We first characterized the number of genes detected according to the quantity of tissue yield and total RNA extracted. Our results revealed that as low as 0.02 mm2 of tissue and 50 pg of total RNA can be used without compromising the number of genes detected. The optimised protocol was used to compare the epidermal versus mesophyll cell transcriptomes of cotyledons at the torpedo-shaped stage of embryo development. The approach was validated by the recovery of well-known epidermal genes such AtML1 or AtPDF2 and genes involved in flavonoid and cuticular waxes pathways. Moreover, the interest and sensitivity of this approach were highlighted by the characterization of several transcription factors preferentially expressed in epidermal cells. CONCLUSION: This technical advance unlocks some current limitations of transcriptomic analyses and allows to investigate further and efficiently new biological questions for which only a very small amounts of cells need to be isolated. For instance, it paves the way to increasing the spatial accuracy of regulatory networks in developing small embryo of Arabidopsis or other plant tissues.

Vicia faba L. in the Bejaia region of Algeria is nodulated by Rhizobium leguminosarum sv. viciae, Rhizobium laguerreae and two new genospecies.

Fifty-eight rhizobial strains were isolated from root nodules of Vicia faba cv. Equina and Vicia faba cv. Minor by the host-trapping method in soils collected from eleven sites in Bejaia, Eastern Algeria. Eleven genotypic groups were distinguished based on the combined PCR/RFLP of 16S rRNA, 16S-23S rRNA intergenic spacer and symbiotic (nodC and nodD-F) genes and further confirmed by multilocus sequence analysis (MLSA) of three housekeeping genes (recA, atpD and rpoB), the 16S rRNA gene and the nodulation genes nodC and nodD. Of the 11 genotypes, 5 were dominant and 2 were the most represented. Most of the strains shared high nodD gene sequence similarity with Rhizobium leguminosarum sv. viciae; their nodC sequences were similar to both Rhizobium leguminosarum and Rhizobium laguerreae. Sequence analyses of the 16S-23S rRNA intergenic spacer showed that all the new strains were phylogenetically related to those described from Vicia sativa and V. faba in several African, European, American and Asian countries, with which they form a group related to Rhizobium leguminosarum. Phylogenetic analysis based on MLSA of 16S rRNA, recA, atpD and rpoB genes allowed the affiliations of strain AM11R to Rhizobium leguminosarum sv. viciae and of strains EB1 and ES8 to Rhizobium laguerreae. In addition, two separate clades with <97% similarity may represent two novel genospecies within the genus Rhizobium.