Saturday, 12 November 2016

Genetic and Genomic Diversity Studies of Acacia Symbionts in Senegal Reveal New Species of Mesorhizobium with a Putative Geographical Pattern

Published Date

Author 



Abstract

Acacia senegal (L) Willd. and Acacia seyal Del. are highly nitrogen-fixing and moderately salt tolerant species. In this study we focused on the genetic and genomic diversity of Acaciamesorhizobia symbionts from diverse origins in Senegal and investigated possible correlations between the genetic diversity of the strains, their soil of origin, and their tolerance to salinity. We first performed a multi-locus sequence analysis on five markers gene fragments on a collection of 47 mesorhizobia strains of A. senegal and A. seyal from 8 localities. Most of the strains (60%) clustered with the M. plurifarium type strain ORS 1032T, while the others form four new clades (MSP1 to MSP4). We sequenced and assembled seven draft genomes: four in the M. plurifarium clade (ORS3356, ORS3365, STM8773 and ORS1032T), one in MSP1 (STM8789), MSP2 (ORS3359) and MSP3 (ORS3324). The average nucleotide identities between these genomes together with the MLSA analysis reveal three new species of Mesorhizobium. A great variability of salt tolerance was found among the strains with a lack of correlation between the genetic diversity of mesorhizobia, their salt tolerance and the soils samples characteristics. A putative geographical pattern of A. senegal symbionts between the dryland north part and the center of Senegal was found, reflecting adaptations to specific local conditions such as the water regime. However, the presence of salt does not seem to be an important structuring factor of Mesorhizobium species.

Introduction

Soil salinization, which might occur naturally or as a consequence of mismanaged irrigation, is a major and growing environmental problem, especially in arid and semi-arid areas of the world. Approximately 800 million hectares of land worldwide is affected by salt [1]. In Senegal, 6% of lands, mainly in coastal areas are affected by the phenomenon of salinization [2].
Salinization is a desertification factor causing degradation of biological, chemical and physical properties of soils [3]. A consequence of this degradation of soil properties is a decrease in fertility, which leads to a reduction in crop yields, land abandonment, and the loss of natural vegetation replaced by huge expanses of saline areas commonly called “tanne” in Senegal. Thus, the supply of services from natural forests becomes insufficient to meet demand, while poverty increases in rural areas.
Reclamation of large areas of saline land through the world seems difficult due to economic and climatic constraints. However, salty soil reclamation could be done using bioremediation method through planting halophyte plants which can absorb salt from soil and enhance the bioproductivity [4,5]. In salted areas of Senegal, despite the disappearance of vegetation islands of legume species with moderate salt tolerance such as AcaciaProsopis and Sesbaniaspecies can be found.
The Acacia genus offers several species adapted to degraded environments, particularly in salt affected areas. A. seyal Del. and A. senegal (L.) Willd. are indigenous woody legumes with important socio-economic roles. First they are widely used in reforestation processes [6] but also for producing gum arabic which is a very important source of income in the Sahel [7,8]. The ability of these species to establish in poor and degraded soils might be due to their aptitude to contract associations with nitrogen-fixing bacteria (rhizobia) and arbuscular mycorrhizal fungi that occur naturally in their rhizosphere [9]. Indeed, symbioses with microorganisms are powerful factors of plant adaptation to adverse environmental conditions, including the lack of major nutrients (nitrogen, phosphorus), biotic (pathogens, phytophagous) and abiotic (drought, salinity) stresses (for review [10]).
Rhizobial populations are known to vary in their tolerance to major environmental factors [11]. It has been reported that salt stress decreases legume growth and nitrogen fixation activity of nodules (for review see [10]). However, inoculation with salt-tolerant strains of rhizobia can enhance the nodulation and nitrogen-fixing ability of the leguminous plants growing under saline conditions, for example in Acacia [1214]. Furthermore, the ability of legume hosts to grow and survive in saline soils was also shown to improve when they were inoculated with salt-tolerant strains of rhizobia [1517].
Strains of the genus Mesorhizobium can establish nitrogen-fixing symbiosis with legume species from temperate, tropical, sub-tropical and arctic areas [18], or associate endophytically with legume plants [19]. Strains isolated from root nodules of the non-legume genus Parasponia [20] and different tropical tree legumes, such as AcaciaLeucaenaProsopisChamaecrista and Lotus in West Africa (Senegal), East Africa (Kenya, Sudan), South America (Mexico, Brazil) and Europe (the Canary island), were described as Mesorhizobium plurifarium[2130]. The large distribution suggests their adaptation to several eco-climatic conditions [31,32]. In a phenotypic comparative study, the M. plurifarium type strain was found more tolerant to heat and salt than other type strains of Mesorhizobium species [31].
Studies based on the sequencing of the 16S rRNA gene and the 16S–23S intergenic spacer of rhizobial strains associated to A. seyal from diverse agro-ecological zones in Senegal reported the predominance of Mesorhizobium genomic groups closely related to Mesorhizobium plurifarium with putatively several new species that remain to be defined [22,33]. Inferring diversity in the Mesorhizobium genus proved to be difficult when using only the 16S rRNA ribosomal taxonomic marker due to its high conservation across Mesorhizobium species [34]. The use of alternative phylogenetic markers is thus very important for species definition within this genus, since some species (e.g. M. mediterraneum and M. temperatum or M. metalliduransand M. gobiense) cannot be distinguished by their 16S rRNA sequence alone [32,35]. A Multi Locus Sequence Analysis (MLSA) with different core genes has been used previously for phylogenetic purposes in the Mesorhizobium genus [3540].
The average nucleotide identity (ANI) of whole genomes has been recently proposed as an alternative to DNA-DNA hybridizations (DDH) to infer bacterial species affiliation with values of ANI >95% on 69% of conserved DNA matching with the 70% species cut-off of DDH usually kept in taxonomic studies of bacteria [41,42]. The rapid development of bacterial genome sequencing at low cost coupled with comparative genomics software development (using either Blast or Mummer algorithms) as jSpecies [43] or MUMi [44] give opportunities to use such correlations to infer rapidly the species of a given strain.
It is often difficult to correlate the genetic diversity of rhizobia and their tolerance to several stresses as salinity. Several studies have underlined the lack of correlation between the sampling sites characteristics, the genetic diversity of rhizobia and their tolerance to salinity [33,45,46], while others authors could link the origin of the soil with the salt tolerance of rhizobia [47].
In this study we analysed the genetic and genomic diversity of A. senegal and A. seyalmesorhizobia symbionts from diverse origins in Senegal and investigated possible correlations between the genetic diversity of the strains, their soil of origin (being under salt-stress or not), and their tolerance to salinity. We first studied at a fine scale the genetic diversity of a collection of mesorhizobia (using Multi Locus Sequence Analysis and genomic fingerprints), inferred their species affiliation using draft genome sequencing and ANI values, and then compared the diversity patterns with salt-tolerance phenotypes and the geographical origin of isolates.

Materials and Methods

Bacterial culture and maintenance

The strains used in this study are listed in Table 1. They originate either from previous studies (See reference in Table 1) or were isolated for this study from nodules on the roots of Acacia senegal or A. seyal growing in pots on soils collected from the field (rhizospheric soil of Acacia). A total of 8 locations (under salt-stress or not, see Fig. 1 and Table 1 for electrical conductivity of soils and gps coordinates), and 36 strains from A. senegal and 11 from A. seyal, were studied. No specific permissions were required for the sampling locations. All strains were kept in 20% (v/v) glycerol at -80°C and cultured either in TY [48] or YEM [49] media in a shaking incubator at 28°C.
thumbnail
Table 1. Bacterial collection of Acacia symbionts used in this study.
thumbnail
Fig 1. Localization of sampling sites in Senegal.
Sampling sites are indicated in blue (D = Dahra, K = Kamb, B = Bambey, F = Foundiougne, Ng = Ngane, V = Vélor, Nd = Ndiafate, No = Nonane, G = Goudiry). The red color indicates salt affected soils (Electroconductivity >4000 μS cm-1). Background colors indicate the different climatic zones of Senegal: light yellow = Sahelian zone (annual rainfall 250–500 mm), Dark yellow = Sudano-Sahelian zone (annual rainfall 500–900 mm), light green = Sudanean zone (annual rainfall 900–1100 mm), dark green = Guinean zone (annual rainfall > 1100 mm).

Phenotypic tests

Tolerance to sodium chloride (NaCl) of rhizobia was tested in 96 well microplates (Nunc Microwell) in broth TY medium [48]. Microplates containing medium supplemented with increasing amounts of NaCl (0 to 600 mM) were inoculated with pure rhizobial culture suspensions (inoculation with 10 μl at optical density (OD) of 1 in a final volume of 200 μl per well, to reach OD 0.05) and incubated on a rotary shaker (160 rpm) at 28°C. Growth was monitored during 72 h by measuring the OD at 600 nm using a microplate spectrophotometer (TECAN-Infinite M200).

Molecular methods

DNA extraction, PCR and sequencing. DNA extraction was performed using a K proteinase lysis protocol as previously described [50]. All PCR amplifications were performed with Go-Taq polymerase (Promega) following manufacturer instructions. The primers used for PCR and sequencing are described in Table 2. The 16S ribosomal DNA (rDNA) was amplified using the universal eubacterial primers FGPS6 and FGPS1509 [51]. The 16S rDNA amplification was carried out as previously described [52]. Fragments of the house-keeping genes atpD, dnaJ, gyrB, and recA were amplified as described before [35,39]. The glnA gene was amplified using either GSI3–58F and GSI2–1143R [53] primers or glnA572F and glnA1143R [54]. PCR products were purified and sequenced by Genoscreen Inc.
thumbnail
Table 2. Primers used in this study for PCR and Sanger sequencing.
Draft genome sequencing and assembly. For whole genome sequencing, strains were grown in 50 ml of broth Yeast-mannitol medium and DNA isolation was performed using a CTAB (Cetyl trimethyl ammonium bromide) bacterial genomic DNA isolation method. DNA quality and quantity was evaluated on a Nanodrop spectrophotometer. The draft genome of 7 strains was produced using Illumina Hiseq technology at Montpellier Genomix platform (MGX), with paired-end sequencing on 700 bp fragments, read length of 100 bp, at 2000X coverage (1/3 Illumina lane per bacteria, around 60 million reads per strain). Sequences were assembled on CLC Genomic workbench v5, and contigs were filtered on size (>500 bp) and reads coverage (minimum of 200X, average cover at 2000X). Genomes were automatically annotated using the Microscope platform [55]. The comparative genomic study of the 7 strains is part of a separate article.
Rep PCR amplification. Repetitive extragenic palindromic PCR (Rep-PCR) genomic fingerprinting was generated with primers REP1R and REP21, as previously reported [56]. PCR mix and conditions were as described by Mishra et al [57]. Rep PCR amplification product were electrophoresed in a 1% agarose gel in TAE 1X buffer and 0.5 μg of ethidium bromide per ml and photographed under UV light.

Phylogenetic analyses

Gene fragments sequences were corrected with Chromas Pro v1.33 software (Technelysium) and aligned using either ClustalX [58] or muscle as implemented in MEGA6 [59]. Alignments were corrected manually under Genedoc software [60] when necessary, and recombination in the datasets was evaluated using the Recombination Detection Program (RDP) v4.35 [61]. Recombination was inferred as true when at least two programs of RDP (RDP, Geneconv, Bootscan, MaxChi, Chimaera, Siscan or 3Seq) could detect the same event. Single marker phylogenies were built with MEGA6 using either Neighbor-Joining (with Kimura 2 distance correction method [62] or Maximum likelihood analyses with 1000 bootstrap replicates. A Bayesian phylogenetic tree was built from the concatenate of all 5 gene fragment alignments, using a Markov chain Monte-Carlo (MCMC) analysis. The priors used for the MCMC analysis were based on a GTR+I+G model with 6 types of substitutions, with parameters estimated by maximum likelihood with Modeltest 3.6 [63].

Average nucleotide identities of whole genomes

Whole genome comparisons using Average Nucleotide Identities (ANI) between genome sequences were produced using jSpecies [43] that use both Blast and Mummer alignments to evaluate whole genome homologies. A blast approach on 1000 bp windows was preferred due to the draft status of genomes. Available genomes of Mesorhizobium were included in the analysis (see Table 3 for Accession numbers). Cut-offs for species delineation were 95% ANI on 69% of conserved DNA according to Goris et al. [41].
thumbnail
Table 3. Genome characteristics of strains sequenced as part of this study as well as reference genomes used for average nucleotide identities calculations (ANI).

Statistical analysis

Basic statistics as well as multiple correspondence analyses (MCA) were conducted in R v3.1.0 software using Ade4 and FactoMineR modules. Quantitative variables (soil pH, soil electro-conductivity and strain tolerance to salinity were transformed into qualitative variables (using different classes). We tested for significant differentiation of populations (in terms of species proportion) between the two sahelian and sudanian climatic zones, and among the three groups of salt tolerance that contained at least three individuals (200, 300 and 400 mM). Unbiased P values were estimated with an exact G tests implemented in Genepop [64].

Accession numbers

The sequences determined in this study have been deposited in the GenBank database and accession numbers are indicated in S1 Table: 16S rRNA (KJ609577-KJ609603); atpD (KJ648182-KJ648227); dnaJ (KJ648228-KJ648270); glnA (KJ883298-KJ883340); gyrB (KJ648271-KJ648315); rec(KJ609604-KJ609649). Genome assemblies produced in this study have been deposited in the European Nucleotide Archive under Project numbers PRJEB6721 to PRJEB726 and accession numbers ERP006359 to ERP006364.

Results

Multi-Locus Sequence Analysis of the Mesorhizobium collection isolated from two Acacia species from salt-contrasted soils

We first sequenced 1000 bp of the 16S rRNA marker (encompassing its variable part) in all strains (except those already sequenced) and produced a 16S rRNA phylogeny presented in Fig. 2. The 16S-derived phylogeny was not well resolved within the Mesorhizobium plurifariumclade, with poor bootstrap values at mostly all branch nodes, and long branches for several strains indicating nucleotide divergence.
thumbnail
Fig 2. Phylogeny of the 16S rRNA marker of the Mesorhizobium collection, built by neighbor joining from a distance matrix corrected by the Kimura-2 method.
The scale bar indicates the number of substitutions per site. Numbers at tree nodes indicate % of bootstrap replicates.
In order to better resolve the genetic diversity of the collection, we produced a multi locus sequence analysis on five house-keeping gene fragments (recA, gyrB, glnA, dnaJ, atpD). Gene markers were chosen according to their use and performance in previous Mesorhizobiumdiversity studies. The vast majority of gene fragments could be amplified and sequenced in almost all strains using methods and primer sets listed in Mat&Methods and Table 2, respectively. Aligned partitions of each gene were built, including reference strains (either type strains of Mesorhizobium species and available genomes). The list of reference strains used is presented in S1 Table. Recombination within datasets was assessed to remove sequences displaying horizontal gene transfer (see Mat&Methods for details). Recombination was detected in gyrB of ORS3369, glnA of ORS3404, and dnaJ of ORS3448. These sequences were removed from the datasets to avoid conflicting phylogenetic signals.
The alignments of all genes were concatenated to produce a global alignment of 2637 bp, with atpD (1–421 bp), dnaJ (422–1195 bp), glnA (1196–1606 bp), gyrB (1607–2237 bp) and recA (2238–2637 bp). A Bayesian phylogeny was then built using the prior (estimated by ML) and run parameters as shown on Fig. 3A, and a consensus tree was built. Bootstraps from 1000 replicates obtained from another analysis by maximum likelihood (the ML model being the same than used for the priors of the Bayesian study on all markers) were added to the tree nodes on Fig. 3A. The phylogenetic tree obtained was much more resolved than the 16S rRNA tree, with Mesorhizobium strains from A. senegal and A. seyal being splitted in five clades supported by high posterior probabilities and bootstraps from the Bayesian and ML analyses, respectively. These clades were named M. plurifarium (MP) for clade I, and MSP1 to MSP4 for the others clades as they did not include any known type strain of Mesorhizobium species. The first clade includes the highest number of strains (29) from the collection together with the type strain of Mesorhizobium plurifarium ORS 1032T. This clade includes both strains isolated from A. seyal (9 strains) and A. senegal (20 strains), unlike other clades which contain only strains isolated from A. senegal (MSP1, MSP4) or A. seyal (MSP2, MSP3). The clade MSP1 includes 14 strains, while MSP2 to MSP4 were rare as they comprise one to two strains.
thumbnail
Fig 3. Bayesian phylogeny of the 5 markers alignment (atpD-dnaJ-glnA-gyrB-recA) (A) and rep-PCR fingerprints (B) of the Mesorhizobium collection.
The priors (estimated by ML) and run parameters for the Bayesian phylogeny are indicated on the left side of the tree. The tree shown is a 50% majrule consensus of all trees produced by the 4 Markov chains. Posterior probabilities (upper number) and bootstraps from 1000 replicates from a ML phylogeny (below number) are indicated at each tree nodes (/ indicate that there is no bootstrap available because the tree node was not common between the Bayesian tree and the ML tree). The scale bar indicates the number of substitutions per site. For the rep-PCR fingerprints (B), see the Mat&Methods section. Each profile was numbered on the right side of the gel. Rep-PCR profiles were not obtained for ORS3400, ORS3598 and ORS1032.

Genomic diversity of bacteria assessed by Rep-PCR fingerprints

Molecular typing of Mesorhizobium strains by Rep-PCR amplification generated multiple amplification products ranging in size from 200 bp to 5000 bp (Fig. 3B). The results show a high intraspecific variability in the Mesorhizobium clades, with very few identical genome fingerprints among strains. This result mirrors the high number of sampling sites and species diversity and the almost absence of clonal isolates in our collection.

Genome sequencing of representative strains & Average Nucleotide identities

Draft genomes of strains were produced in order to evaluate the genomic diversity of M. plurifarium strains nodulating A. seyal and A. senegal as well as to determine if several species were present in the collection. Whole genome of four strains of M. plurifarium (including the type strain ORS1032T) and three strains of Mesorhizobium sp. belonging to clade MSP1 to MSP3 were sequenced. No genome was chosen in clade MSP4 as our preliminary analyses did not detect this new clade, and our criteria included also the salt tolerance and strains in this group did not exhibit interesting phenotypes. Information on the sequencing method is given in the Mat&Method section, and descriptive information on the genomes can be found in Table 3.
Average nucleotide identities (ANI) were calculated to evaluate species affiliation of strains which genomes were sequenced. According to Richter and Rossello-Mora [43] there is a correlation between the percentage of average nucleotide identity (ANI %) and the percentage of DNA-DNA hybridization, which is a major criteria in bacterial species delineation. The Table 4contain the matrix of nucleotide identities between whole genomes calculated under jSpecies [43] by the blast method on windows of 1000 bp in size. We colored in gray the ANI% value of strains which genomes identities exhibit the criteria of belonging to the same species as defined by Goris et al [41]. The latter article reported that when strains share an ANI % > 95% on more than 69% of conserved DNA, they would belong to the same species. According to these criteria, ORS3356, ORS3365 and DJ20 share more than 95% ANI with the type strain of M. plurifarium ORS1032T, and are classified in clade MP in the MLSA tree. The strains STM8789, ORS3359 and ORS3324, respectively belonging to MSP1, MSP2 and MSP3, did not share these criteria with any Mesorhizobium species included in the analysis.
thumbnail
Table 4. Average nucleotide identities (ANI) between Mesorhizobium genomes
% of ANI were calculated using blastn on 1000 bp windows of the genomes, with jSpecies (see Math&Methods). In bold are indicated ANI matching with the species affiliation cut-off: >95%, with the % of conserved DNA indicated between parentheses (>69%), as defined by Goris et al. and Konstantinidis & Tiedje [41,42].

Tolerance to salinity of the collection

Salt tolerance test of all strains (in triplicates) were conducted in microplate in TY broth medium supplemented with 0 to 600 mM of NaCl. Our aim was to establish putative correlations between the salt tolerance phenotype, the species affiliation and/or geographical origin of the isolates. We considered that tolerance of a given strain to a particular concentration of salt was true when the bacteria could still grow at 50% compared to its growth in the same medium without salt (considered as the 100% of growth). As the collection was quite high in number with tests in 7 different concentrations (0, 100, 200, 300, 400, 500 mM) and 3 OD sampling times (24, 41 and 68 hours), we chose to represent the data of the OD growth of the strain at 41h post inoculation of the medium, using a radar graphic, shown in Fig. 4. The same radar at 24h and 68h is presented in S1 Fig., and histograms of growth% compared to the unsalted control with standard deviation are given for every strain in S2 Fig. For M. plurifarium (MP) and MSP1 species, half of the strains were highly impacted in their growth at 300 mM (growth <50% compared to control) and most could not grow (OD not exceeding 0.3) at 400 mM of salt in the medium. Some strains exhibited a better tolerance to salinity and could still grow up to 400 mM (ORS1032, Dj20, Sd11, Sod14, Sod10 and Sod15). Some strains could even grow at 500 mM of salt: Dj16 and Dj20 for MP, Sod10 for MSP1. The few strains of MSP2, MSP3 and MSP4 were poorly tolerant to salt (even at 200 mM), and strains of MSP4 did not even grew well at low salt concentrations, probably due to the growth conditions that might be unadapted to this species.
thumbnail
Fig 4. Radar of Mesorhizobium strain optical density when grown in TY broth medium in 200 ul in a 96 well microplate with different concentrations of NaCl (100 to 500 mM).
The data presented here are at 41h post-inoculation with 0.05 OD as a starting point. Radars at 24h and 68h are presented in S1 Fig., and histograms of growth for each strain at 41h are presented in S2 Fig.

Correspondence Analyses between salt responses, genetic typing and the origin of strains

We performed a multivariate correspondence analysis (MCA) in order to study the distribution of our strains according to the following parameters: tolerance of strains, species, plant host, and pedoclimatic conditions of the sampling sites (soil pH and electro-conductivity (EC), geographical origin, eco-geographical and climatic zones). Quantitative variables were transformed in qualitative variable using different classes with biological relevance (pH, tolerance to salinity, EC). The results are presented as a MCA plot in S3 Fig. (A) and factors from the MCA are represented individually in S3 Fig. (B). As observed, there was no correlation between the distribution of strains and their salt tolerance, their species affiliation, and/or their geographical or soil of origin (being salted or not). We tested for significant differentiation of populations (in terms of species proportion) between the two sSahelian and Sudano-sahelian climatic zones, and among the three groups of salt tolerance that contained at least three individuals (200, 300 and 400 mM). There was no significant difference among the three groups of salt tolerance (all pairs of population, p values from 0.45 to 1). Conversely, the two climatic zones (marked by different annual rainfall amounts) were significantly differentiated (p = 0), reflecting a different composition of species between them. Indeed the MP, MSP2 and MSP3 species were only detected in the Sudano-Sahelien zone, while MSP1 was detected in both Sudano-Sahelian and Sahelian zones, and MSP4 was only found in the Sahelian zone. However the species MSP2, MSP3 and MSP4 contain very few strains and thus we cannot conclude about their geographical pattern (though we can conclude that these species are quite rare compared to MSP1 and MP). For MP and MSP1 the geographical pattern seems to be clear and thus related to the eco-climatic zones (rainfall) rather than the salinity of the soils.

Discussion

Mesorhizobium strains of Acacia seyal and A. senegal belong mainly to M. plurifarium but also to at least three new species

The MLSA phylogeny based on 5 housekeeping genes showed a higher species diversity of Mesorhizobium strains nodulating A. seyal and A. senegal than previously expected. Most of the strains (60%) clustered with the M. plurifarium type strain ORS 1032T, as previously reported for some of them [2225,27]. The rest of the strains were distributed in four new clades, which happen to be new species for at least 3 of them (MSP1 to MSP3) according to the average nucleotide identities of genomes of representative strains in each clade. Unfortunately we did not sequence a representative strain of MSP4 as we did not initially expect this clade in the first results of the single markers phylogenies. The ANI has been proposed as an alternative to DNA-DNA hybridizations (DDH) to infer bacterial species affiliation ([41,42]. If we apply the cut-offs of species delineation as previously published (ANI>95% on 69% of conserved DNA), then the MLSA phylogeny of Acacia mesorhizobia fits perfectly with the ANI-based species affiliation. The use of ANI seems to be a good species assessment in the Mesorhizobium genus, as there was also no conflict of ANI values between known Mesorhizobium species (Table 4, all known species comparison gave ANI <95%). The ANI between reference strains also confirms the belonging of the strain MAFF303099 to a separate species from M. loti (89.35% ANI with M. loti USDA3471T) as suggested by Turner et al. [65] and Wang et al. [66]. Our study also confirms the good performance of the 5 markers used (atpD, dnaJ, gyrB, recA and glnA) in Mesorhizobium diversity studies as previously reported for rhizobia [28,3537,39,53,67].
Our results thus suggest three new species of Mesorhizobium nodulating both Acacia species. Two of the new clades were anticipated to be separate species from MP as they clustered separately in the IGS spacer phylogeny (as ORS3359 and ORS3324) of the Diouf et al. [22] study, though not in their 16S rDNA phylogeny (this study, [22], [24]). It is interesting to note that MSP1 was only detected in A. senegal nodules, while MP was found in the nodules of both Acacia species. A recent analysis of nodulation genes clustered Mesorhizobium strains (MP, MSP2, MSP3) according to their host of origin: strains from the MSP1 clade (all from A. senegal) were clustered with M. plurifarium strains from A. senegal in their nodA, nodC and nifH phylogenies (Bakhoum et al., Microbial Ecology, in press). In the case of A. seyal strains, these are all grouped together in the nodA and nodC phylogenies [33]. The only exception of this latter study was the nodC of ORS3324 that grouped with the nodC of the Ensifer arboristype strain. As we sequenced the full genome of ORS3324, we analyzed its nodC in the genome data and found out that it grouped together with the other A. seyal mesorhizobia (ORS3359 and ORS3324 sharing 100% nucleotide on nodA and nodC), thus the nodC fragment published in Diouf et al. [33] shall be considered as an error.
The nodulation ability of strains was assessed on A. seyal and A. senegal and all strains were able to nodulate both species, whatever their species affiliation or their nodA allele. Such ability is correlated with previous articles showing that strains of the large M. plurifarium clade (also described as Cluster U in previous studies) share similar nodulation host range on AcaciaProsopis and Leucaena species [25,68], and thus their host range shall not explain their geographical distribution.

The tolerance to NaCl is highly variable among M. plurifarium and MSP1 species—whatever their geographical and pedoclimatic origin

A great variability of response to salt of Acacia mesorhizobia from MP and MSP1 was found. Conversely, strains of MSP2 to MSP4 did not exhibit high salt tolerance, but given the low number of strains in each of these species it is difficult to conclude to a species effect. Taken altogether, the mesorhizobia strains were mostly salt tolerant at 200 mM (i.e. with a growth at least of 50% compared to the control without salt), as 95% and 43% of strains could still grow at 200 mM and 300 mM of NaCl, respectively. For the species MP, 100% of strains could grow at 200 mM, while 48% and 14% of strains tolerate 300 mM and 400 mM of NaCl, respectively, and one (STM8773) tolerate 500 mM. It has been previously shown that the M. plurifarium type strain ORS1032T was more tolerant to heat and salt than several other type strains of Mesorhizobium species [31]. In our study, ORS1032T tolerate 400 mM of salt, together with others strains. The MP species could thus be more tolerant than others Mesorhizobiumspecies, but a larger study including many strains from the different species is required to infer this question. Another species detected in this study, MSP1, seems well adapted to salt tolerance. The ability of some symbiotically effective strains to tolerate high salinity is promising with regard to improving host plant reestablishment in salt affected soils [46,69].
We investigated possible links between the genetic diversity of mesorhizobia, their salt tolerance and the soils samples characteristics. However, we found a total lack of correlation between these parameters. The most salt-tolerant strains (up to 500 mM, as STM8773 and STM8789) were isolated from non-saline areas. These results are similar to previous studies of Thrall et al. [46] who found also salt-tolerant strains (up to 400 mM) from nodules of Acacia decurrens in non-saline sites. Several studies have underlined the lack of correlation between the sampling sites characteristics, the genetic diversity of rhizobia and the tolerance to salinity of rhizobial isolates ([33,45,46]. In Senegal, Diouf et al [22] found a weak influence of soil characteristics (pH and salt) on the distribution of rhizobial populations of A. seyal in the Groundnut Basin, as observed in our study at a larger scale. A possible explanation for this lack of correlation could be that the soil is not homogeneous as shown by van Asten et al. [70] who found large differences in salinity and alkalinity levels at short distances in salt affected areas. Such conditions would allow the persistence of various phenotypic traits among rhizobial populations. Rhizobia might also be protected in the nodules, explaining why non-tolerant strains can be detected in highly saline soils. The accumulation of certain compounds, including osmoprotectants (as trehalose and Poly β-hydroxybutyrate) and compatible solutes, may also increase the osmotic tolerance of rhizobia [7173].
On the other hand, we found a putative geographical pattern of A. senegal symbionts between the dryland north part and the center of Senegal. The MSP1 strains were found both in the center and north part of the country while MSP4 and M. plurifarium species were found only in the north and the center part of the country, respectively. Such specific distribution of Mesorhizobium species has also been observed on Caragana spp. symbionts in three eco-regions of China [74]. The geographical pattern of symbionts observed in our study could be linked to the annual rainfall characteristics of each site, as observed for Vigna unguiculatasymbionts in Senegal [75]. Indeed, the M. plurifarium strains were only found in the center of the Sudano-Sahelian zone where annual rainfall is between 500 and 900 mm, while no strain of this clade was found in the Kamb and Dahra sites, located in the arid regions of the sahelian zone where annual rainfall varied between 250 and 500 mm (Fig. 1). This result implies that M. plurifarium strains would not be well adapted to dryland conditions. On the other hand, the MSP1 strains from the Kamb soil were shown to be tolerant to water stress when testing their growth with different concentrations of polyethylene glycol [24]. Wade et al [75] found also an eco-geographical diversity of cowpea bradyrhizobia in Senegal marked by the dominance of two genetic types depending on annual rainfall. These authors concluded a possible role of the water regime and the pH in shaping the cowpea bradyrhizobia genotypic distribution, noticing that strains isolated from the northern region were, generally, more adapted to water stress and slightly alkaline soils. The observed geographical distribution of Acacia mesorhizobia could thus reflect particular adaptations of each species to specific local conditions as the water regime, but the presence of salt does not seem to be an important structuring factor of Mesorhizobiumspecies.

Supporting Information

S1_Table.docx
Acacia mesorhizobia from Senegal and type strain of rhizobial species with corresponding GenBank records (accesssion numbers or Gene ID).

S1 Table. Genbank accession numbers of all sequences from this study.

Acacia mesorhizobia from Senegal and type strain of rhizobial species with corresponding GenBank records (accesssion numbers or Gene ID).
doi:10.1371/journal.pone.0117667.s001
(DOCX)

S1 Fig. Radar of Mesorhizobium strain optical density when grown in TY broth medium with different concentrations of NaCl (100 to 500 mM).

The data presented here are at 24h, 41h and 68h post-inoculation with 0.05 OD as a starting point. Each color indicates a concentration of NaCl.
doi:10.1371/journal.pone.0117667.s002
(DOCX)

S2 Fig. Histogram of growth for each strain when grown in TY broth medium with different concentrations of NaCl (100 to 500 mM), at 41h post-inoculation.

A) Optical density, B) % of growth compared to the same strain growth in TY without salt (B). Standard deviations were calculated from triplicates wells. The arrows indicate the names of the best tolerant strains of Mesorhizobium in our bacterial growing conditions.
doi:10.1371/journal.pone.0117667.s003
(DOCX)

S3 Fig. Multivariate correspondence analysis of strains, geographical and soil factors.

A) Multivariate correspondence analysis factor map for Mesorhizobium strains (in blue) with factors: strain tolerance (T100 to T500 corresponding to tolerance to 100 to 500 mM of NaCl, in green), pH of soil of origin (in medium blue, defined into 5 classes: 4, 5, 6, 6.5, 7), species assignment (in black, MP for M. plurifarium, MSP1 to MSP4 as defined in the study), name of site of origin (in red), presence of salt in soil of origin (in pale blue, S = Salted (Electroconductivity>4000 μS.cm), MS = moderately salted (400<EC<4000), NS = not salted (EC<400)), plant host (in orange), Ecogeographic zones (in brown), and climatic zone (in pink with Sah for Sahelian and Sud for sudano-sahelian zone). The MCA analysis was performed with MineFactorR under R software. B) Distribution of each factor in the MCA analysis. For Species: MP means M. plurifarium, MSP1 to MSP4 correspond to the new species. For NaCl tolerance: the number indicates the concentration in mM in the medium at which the bacteria could still grow >50% compared to the negative control without salt. For soil pH: exact values were transformed into classes by agglomerating close values to the value of 4, 5, 6, 6.5, and 7. For soil salt content, NS means not salted (Electroconductivity <400 uS.cm), MS moderately salted (400<EC<4000 uS.cm), and S salted (EC >4000 uS.cm, as defined by F.A.O.) For Clim (climatic zone): Sah = sahelian zone as defined in the map in Fig. 1. For Ecogeo (Ecogeographical region). BassinArachidier means groundnut basin. For Geo: Geographical origin, indicate the locality of origin of the strains.
doi:10.1371/journal.pone.0117667.s004
(DOCX)

Acknowledgments

We thank Pierre Tisseyre and Mathieu Ndigue Faye for their help in microbiology procedures and strain collection, and the LABGeM and the National Infrastructure « France Genomique » for Genome annotation and data incorporation in the Microscope platform.

Author Contributions

Conceived and designed the experiments: DD MN LM. Performed the experiments: FD AK ALQ NB DF HP MD. Analyzed the data: FD DD LM. Contributed reagents/materials/analysis tools: DD DF HP IN LM. Wrote the paper: FD DD LM.

References

  1. 1.Munns R (2002) Comparative physiology of salt and water stress. Plant Cell Environ 25: 239–250. pmid:11841667 doi: 10.1046/j.0016-8025.2001.00808.x 
  2. 2.Barbiero L, Mohamedou A, Laperrousaz C, Furian S, Cunnac S (2004) Polyphasic origin of salinity in the Senegal delta and middle valley. Catena 58: 101–124. doi: 10.1016/j.catena.2004.03.003 
  3. 3.Qadir M, Schubert S (2002) Degradation processes and nutrient constraints in sodic soils. Land Degradation & Development. 294: 275–294. doi: 10.1126/science.1255641. pmid:25593191 
  4. 4.De Lacerda CF, Cambraia J, Oliva MA, Ruiz HA (2005) Changes in growth and in solute concentrations in sorghum leaves and roots during salt stress recovery. Environ Exp Bot 54: 69–76. doi: 10.1016/j.envexpbot.2004.06.004 
  5. 5.Singh G (2009) Salinity-related desertification and management strategies: Indian experience. Land Degradation & Development. 385: 367–385. doi: 10.1126/science.1255641. pmid:25593191 
  6. 6.Midgley JJ, Bond WJ (2001) A synthesis of the demography of African acacias. J Trop Ecol 17: 871–886. doi: 10.1017/s026646740100164x 
  7. 7.Motlagh S, Ravines P, Karamallah KA, Ma Q (2006) The analysis of Acacia gums using electrophoresis. Food Hydrocoll 20: 848–854. doi: 10.1016/j.foodhyd.2005.08.007 
  8. 8.Al-Assaf S, Phillips GO, Aoki H, Sasaki Y (2007) Characterization and properties of Acacia senegal (L.) Willd. var. senegal with enhanced properties (Acacia (sen) SUPER GUMTM): Part 1—Controlled maturation of Acacia senegal var. senegal to increase viscoelasticity, produce a hydrogel form and convert a poor into a good emulsifier. Food Hydrocoll 21: 319–328. doi: 10.1016/j.foodhyd.2006.04.011 
  9. 9.Dommergues YR (1995) Nitrogen fixation by trees in relation to soil nitrogen economy. Fertil Res 42: 215–230. doi: 10.1007/bf00750516 
  10. 10.Zahran HH (1999) Rhizobium-legume symbiosis and nitrogen fixation under severe conditions and in an arid climate. Microbiol Mol Biol Rev 63: 968–89. pmid:10585971 
  11. 11.Wei G-H, Yang X-Y, Zhang Z-X, Yang Y-Z, Lindström K (2008) Strain Mesorhizobiumsp. CCNWGX035: A Stress-Tolerant Isolate from Glycyrrhiza glabra Displaying a Wide Host Range of Nodulation. Pedosphere 18: 102–112. doi: 10.1016/s1002-0160(07)60108-8 
  12. 12.Diouf D, Duponnois R, Tidiane Ba A, Neyra M, Lesueur D (2005) Symbiosis of Acacia auriculiformis and Acacia mangium with mycorrhizal fungi and Bradyrhizobium spp. improves salt tolerance in greenhouse conditions. Funct Plant Biol 32: 1143. doi: 10.1071/fp04069 
  13. 13.Coba de la Peña T, Pueyo JJ (2011) Legumes in the reclamation of marginal soils, from cultivar and inoculant selection to transgenic approaches. Agron Sustain Dev 32: 65–91. doi: 10.1007/s13593-011-0024-2 
  14. 14.Lal B, Khanna S (1994) Selection of salt-tolerant Rhizobium isolates of Acacia nilotica. World J Microbiol Biotechnol 10: 637–639. doi: 10.1007/BF00327949. pmid:24421185 
  15. 15.Thrall PH, Millsdom DA, Jeavons AC, Waayers M, Harvey GR, et al. (2005) Seed inoculation with effective root-nodule bacteria enhances revegetation success. J Appl Ecol 42: 740–751. doi: 10.1111/j.1365-2664.2005.01058.x 
  16. 16.Zou N, Dart PJ, Marcar NE (1995) Interaction of salinity and rhizobial strain on growth and N fixation by Acacia ampliceps. Soil Biol Biochem 27: 409–413. doi: 10.1016/0038-0717(95)98611-q 
  17. 17.Shamseldin A, Werner D (2005) High salt and high pH tolerance of new isolated Rhizobium etli strains from Egyptian soils. Curr Microbiol 50: 11–16. pmid:15702257 doi: 10.1007/s00284-004-4391-7 
  18. 18.Chen W-M, Zhu W-F, Bontemps C, Young JPW, Wei G-H (2011) Mesorhizobium camelthorni sp. nov., isolated from Alhagi sparsifolia. Int J Syst Evol Microbiol 61: 574–579. doi: 10.1099/ijs.0.022947-0. pmid:20382787 
  19. 19.Wei G-H, Yang X-Y, Zhang J-W, Gao J-M, Ma Y-Q, et al. (2007) Rhizobialide: a new stearolactone produced by Mesorhizobium sp. CCNWGX022, a rhizobial endophyte from Glycyrrhiza uralensis. Chem Biodivers 4: 893–898. pmid:17510984 doi: 10.1002/chin.200738201 
  20. 20.Op den Camp RHM, Polone E, Fedorova E, Roelofsen W, Squartini A, et al. (2012) Nonlegume Parasponia andersonii deploys a broad rhizobium host range strategy resulting in largely variable symbiotic effectiveness. Mol Plant Microbe Interact 25: 954–963. doi: 10.1094/MPMI-11-11-0304. pmid:22668002 
  21. 21.Wang ET, Kan FL, Tan ZY, Toledo I, Chen WX, et al. (2003) Diverse Mesorhizobium plurifarium populations native to Mexican soils. Arch Microbiol 180: 444–454. pmid:14576977 doi: 10.1007/s00203-003-0610-z 
  22. 22.Diouf D, Samba-Mbaye R, Lesueur D, Ba AT, Dreyfus B, et al. (2007) Genetic diversity of Acacia seyal Del. rhizobial populations indigenous to Senegalese soils in relation to salinity and pH of the sampling sites. Microb Ecol 54: 553–566. pmid:17406772 doi: 10.1007/s00248-007-9243-0 
  23. 23.Bakhoum N, Ndoye F, Kane A, Assigbetse K, Fall D, et al. (2012) Impact of rhizobial inoculation on Acacia senegal (L.) Willd. growth in greenhouse and soil functioning in relation to seed provenance and soil origin. World J Microbiol Biotechnol 28: 2567–2579. doi: 10.1007/s11274-012-1066-6. pmid:22806163 
  24. 24.Fall D, Diouf D, Ourarhi M, Faye A, Abdelmounen H, et al. (2008) Phenotypic and genotypic characteristics of Acacia senegal (L.) Willd. root-nodulating bacteria isolated from soils in the dryland part of Senegal. Lett Appl Microbiol 47: 85–97. doi: 10.1111/j.1472-765X.2008.02389.x. pmid:18565139 
  25. 25.De Lajudie P, Willems A, Nick G, Moreira F, Molouba F, et al. (1998) Characterization of tropical tree rhizobia and description of Mesorhizobium plurifarium sp. nov. Int J Syst Bacteriol 48 Pt 2: 369–382. pmid:9731275 doi: 10.1099/00207713-48-2-369 
  26. 26.Sarr A, Neyra M, Houeibib MAO, Ndoye I, Oihabi A, et al. (2005) Rhizobial populations in soils from natural Acacia senegal and Acacia nilotica forests in Mauritania and the Senegal river valley. Microb Ecol 50: 152–162. pmid:16184338 doi: 10.1007/s00248-004-0077-8 
  27. 27.Ba S, Willems A, de Lajudie P, Roche P, Jeder H, et al. (2002) Symbiotic and taxonomic diversity of rhizobia isolated from Acacia tortilis subsp. raddiana in Africa. Syst Appl Microbiol 25: 130–145. pmid:12086180 doi: 10.1078/0723-2020-00091 
  28. 28.Lorite MJ, Donate-Correa J, del Arco-Aguilar M, Pérez Galdona R, Sanjuán J, et al. (2010) Lotus endemic to the Canary Islands are nodulated by diverse and novel rhizobial species and symbiotypes. Syst Appl Microbiol 33: 282–290. doi: 10.1016/j.syapm.2010.03.006. pmid:20447791 
  29. 29.Hoque MS, Broadhurst LM, Thrall PH (2011) Genetic characterization of root-nodule bacteria associated with Acacia salicina and A. stenophylla (Mimosaceae) across south-eastern Australia. Int J Syst Evol Microbiol 61: 299–309. doi: 10.1099/ijs.0.021014-0. pmid:20228207 
  30. 30.Odee DW, Haukka K, McInroy SG, Sprent JI, Sutherland JM, et al. (2002) Genetic and symbiotic characterization of rhizobia isolated from tree and herbaceous legumes grown in soils from ecologically diverse sites in Kenya. Soil Biol Biochem 34: 801–811. doi: 10.1016/s0038-0717(02)00009-3 
  31. 31.Laranjo M, Oliveira S (2011) Tolerance of Mesorhizobium type strains to different environmental stresses. Antonie Van Leeuwenhoek 99: 651–662. doi: 10.1007/s10482-010-9539-9. pmid:21152981 
  32. 32.Laranjo M, Alexandre A, Oliveira S (2014) Legume growth-promoting rhizobia: an overview on the Mesorhizobium genus. Microbiol Res 169: 2–17. doi: 10.1016/j.micres.2013.09.012. pmid:24157054 
  33. 33.Diouf D, Fall D, Chaintreuil C, Ba AT, Dreyfus B, et al. (2010) Phylogenetic analyses of symbiotic genes and characterization of functional traits of Mesorhizobium spp. strains associated with the promiscuous species Acacia seyal Del. J Appl Microbiol 108: 818–830. doi: 10.1111/j.1365-2672.2009.04500.x. pmid:19735328 
  34. 34.Nandasena KG, O’Hara GW, Tiwari RP, Willems A, Howieson JG (2009) Mesorhizobium australicum sp. nov. and Mesorhizobium opportunistum sp. nov., isolated from Biserrula pelecinus L. in Australia. Int J Syst Evol Microbiol 59: 2140–2147. doi: 10.1099/ijs.0.005728-0. pmid:19605705 
  35. 35.Gaunt MW, Turner SL, Rigottier-Gois L, Lloyd-Macgilp SA, Young JP (2001) Phylogenies of atpD and recA support the small subunit rRNA-based classification of rhizobia. Int J Syst Evol Microbiol 51: 2037–2048. pmid:11760945 doi: 10.1099/00207713-51-6-2037 
  36. 36.Weir BS, Turner SJ, Silvester WB, Park D-C, Young JM (2004) Unexpectedly diverse Mesorhizobium strains and Rhizobium leguminosarum nodulate native legume genera of New Zealand, while introduced legume weeds are nodulated by Bradyrhizobiumspecies. Appl Environ Microbiol 70: 5980–5987. pmid:15466541 doi: 10.1128/aem.70.10.5980-5987.2004 
  37. 37.Vinuesa P, Silva C, Lorite MJ, Izaguirre-Mayoral ML, Bedmar EJ, et al. (2005) Molecular systematics of rhizobia based on maximum likelihood and Bayesian phylogenies inferred from rrs, atpD, recA and nifH sequences, and their use in the classification of Sesbania microsymbionts from Venezuelan wetlands. Syst Appl Microbiol 28: 702–716. pmid:16261860 doi: 10.1016/j.syapm.2005.05.007 
  38. 38.Laranjo M, Alexandre A, Rivas R, Velázquez E, Young JPW, et al. (2008) Chickpea rhizobia symbiosis genes are highly conserved across multiple Mesorhizobium species. FEMS Microbiol Ecol 66: 391–400. doi: 10.1111/j.1574-6941.2008.00584.x. pmid:18795953 
  39. 39.Alexandre A, Laranjo M, Young JPW, Oliveira S (2008) dnaJ is a useful phylogenetic marker for alphaproteobacteria. Int J Syst Evol Microbiol 58: 2839–2849. doi: 10.1099/ijs.0.2008/001636-0. pmid:19060069 
  40. 40.Stepkowski T, Czaplińska M, Miedzinska K, Moulin L (2003) The variable part of the dnaK gene as an alternative marker for phylogenetic studies of rhizobia and related alpha Proteobacteria. Syst Appl Microbiol 26: 483–494. pmid:14666974 doi: 10.1078/072320203770865765 
  41. 41.Goris J, Konstantinidis KT, Klappenbach JA, Coenye T, Vandamme P, et al. (2007) DNA-DNA hybridization values and their relationship to whole-genome sequence similarities. Int J Syst Evol Microbiol 57: 81–91. pmid:17220447 doi: 10.1099/ijs.0.64483-0 
  42. 42.Konstantinidis KT, Tiedje JM (2007) Prokaryotic taxonomy and phylogeny in the genomic era: advancements and challenges ahead. Curr Opin Microbiol 10: 504–509. pmid:17923431 doi: 10.1016/j.mib.2007.08.006 
  43. 43.Richter M, Rosselló-Móra R (2009) Shifting the genomic gold standard for the prokaryotic species definition. Proc Natl Acad Sci U S A 106: 19126–19131. doi: 10.1073/pnas.0906412106. pmid:19855009 
  44. 44.Deloger M, El Karoui M, Petit M-A (2009) A genomic distance based on MUM indicates discontinuity between most bacterial species and genera. J Bacteriol 191: 91–99. doi: 10.1128/JB.01202-08. pmid:18978054 
  45. 45.Elboutahiri N, Thami-Alami I, Udupa SM (2010) Phenotypic and genetic diversity in Sinorhizobium meliloti and S. medicae from drought and salt affected regions of Morocco. BMC Microbiol 10: 15. doi: 10.1186/1471-2180-10-15. pmid:20089174 
  46. 46.Thrall PH, Broadhurst LM, Hoque MS, Bagnall DJ (2009) Diversity and salt tolerance of native Acacia rhizobia isolated from saline and non-saline soils. Austral Ecol 34: 950–963. doi: 10.1111/j.1442-9993.2009.01998.x 
  47. 47.Lorite MJ, Muñoz S, Olivares J, Soto MJ, Sanjuán J (2010) Characterization of strains unlike Mesorhizobium loti that nodulate lotus spp. in saline soils of Granada, Spain. Appl Environ Microbiol 76: 4019–4026. doi: 10.1128/AEM.02555-09. pmid:20435777 
  48. 48.Sambrook J, Fritsch EF, Manniatis T (1989) Molecular cloning: a laboratory manual. New York, NY: Cold Spring Harb. Lab. 
  49. 49.Vincent J (1970) A manual for the pratical study of root-nodule bacteria. I.B.P. Han. Ltd, Oxford: Blackwell Scientific Publications. 
  50. 50.Moulin L, Béna G, Boivin-Masson C, Stępkowski T (2004) Phylogenetic analyses of symbiotic nodulation genes support vertical and lateral gene co-transfer within the Bradyrhizobium genus. Mol Phylogenet Evol 30: 720–732. pmid:15012950 doi: 10.1016/s1055-7903(03)00255-0 
  51. 51.Normand P, Cournoyer B, Simonet P, Nazaret S (1992) Analysis of a ribosomal RNA operon in the actinomycete Frankia. Gene 111: 119–124. pmid:1372279 doi: 10.1016/0378-1119(92)90612-s 
  52. 52.Moulin L, Munive A, Dreyfus B, Boivin-Masson C (2001) Nodulation of legumes by members of the beta-subclass of Proteobacteria. Nature 411: 948–950. pmid:11418858 doi: 10.1038/35082070 
  53. 53.Turner SL, Young JP (2000) The glutamine synthetases of rhizobia: phylogenetics and evolutionary implications. Mol Biol Evol 17: 309–319. pmid:10677854 doi: 10.1093/oxfordjournals.molbev.a026311 
  54. 54.Martens M, Delaere M, Coopman R, De Vos P, Gillis M, et al. (2007) Multilocus sequence analysis of Ensifer and related taxa. Int J Syst Evol Microbiol 57: 489–503. pmid:17329774 doi: 10.1099/ijs.0.64344-0 
  55. 55.Vallenet D, Engelen S, Mornico D, Cruveiller S, Fleury L, et al. (2009) MicroScope: a platform for microbial genome annotation and comparative genomics. Database (Oxford) 2009: bap021. doi: 10.1093/database/bap021 
  56. 56.Versalovic J, Koeuth T, Lupski JR (1991) Distribution of repetitive DNA sequences in eubacteria and application to fingerprinting of bacterial genomes. Nucleic Acids Res 19: 6823–6831. pmid:1762913 doi: 10.1093/nar/19.24.6823 
  57. 57.Mishra RPN, Tisseyre P, Melkonian R, Chaintreuil C, Miché L, et al. (2012) Genetic diversity of Mimosa pudica rhizobial symbionts in soils of French Guiana: Investigating the origin and diversity of Burkholderia phymatum and other beta-rhizobia. FEMS Microbiol Ecol 79: 487–503. doi: 10.1111/j.1574-6941.2011.01235.x. pmid:22093060 
  58. 58.Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG (1997) The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 25: 4876–4882. pmid:9396791 doi: 10.1093/nar/25.24.4876 
  59. 59.Tamura K, Stecher G, Peterson D, Filipski A, Kumar S (2013) MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol 30: 2725–2729. doi: 10.1093/molbev/mst197. pmid:24132122 
  60. 60.Nicholas KB, Nicholas HB, Dierfield DW (1997) GeneDoc: Analysis and visualization of genetic variation. EMBNEWNEWS 4: 14. 
  61. 61.Martin DP, Lemey P, Lott M, Moulton V, Posada D, et al. (2010) RDP3: a flexible and fast computer program for analyzing recombination. Bioinformatics 26: 2462–2463. doi: 10.1093/bioinformatics/btq467. pmid:20798170 
  62. 62.Saitou N, Nei M (1987) The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 4: 406–425. pmid:3447015 
  63. 63.Posada D, Buckley TR (2004) Model selection and model averaging in phylogenetics: advantages of akaike information criterion and bayesian approaches over likelihood ratio tests. Syst Biol 53: 793–808. pmid:15545256 doi: 10.1080/10635150490522304 
  64. 64.Raymond M, Rousset F (1995) GENEPOP(Version1.2):Population Genetics Software for Exact Tests and Ecumenicism. J Hered 86: 248–249. 
  65. 65.Turner SL, Knight KAL, Young JPW (2002) Identification and analysis of rhizobial plasmid origins of transfer. FEMS Microbiol Ecol 42: 227–234. doi: 10.1111/j.1574-6941.2002.tb01012.x. pmid:19709282 
  66. 66.Wang S, Hao B, Li J, Gu H, Peng J, et al. (2014) Whole-genome sequencing of Mesorhizobium huakuii 7653R provides molecular insights into host specificity and symbiosis island dynamics. BMC Genomics 15: 440. doi: 10.1186/1471-2164-15-440. pmid:24906389 
  67. 67.Tan HW, Weir BS, Carter N, Heenan PB, Ridgway HJ, et al. (2012) Rhizobia with 16S rRNA and nifH similar to Mesorhizobium huakuii but Novel recA, glnII, nodA and nodC genes are symbionts of New Zealand Carmichaelinae. PLoS One 7: e47677. doi: 10.1371/journal.pone.0047677. pmid:23118889 
  68. 68.Lortet G, Méar N, Lorquin J, Dreyfus B, De Lajudie P, et al. (1996) Nod factor thin-layer chromatography profiling as a tool to characterize symbiotic specificity of rhizobial strains: application to Sinorhizobium saheliS. teranga, and Rhizobium sp. strains isolated from Acacia and Sesbania. Mol Plant Microbe Interact 9: 736–747. doi: 10.1094/mpmi-9-0736 
  69. 69.Thrall PH, Bever JD, Slattery JF (2008) Rhizobial mediation of Acacia adaptation to soil salinity : evidence of underlying trade-offs and tests of expected patterns. J Ecol 96: 746–755. doi: 10.1111/j.1365-2745.2008.01381.x 
  70. 70.Van Asten PJ, Barbiéro L, Wopereis MC, Maeght J, van der Zee SEAT (2003) Actual and potential salt-related soil degradation in an irrigated rice scheme in the Sahelian zone of Mauritania. Agric Water Manag 60: 13–32. doi: 10.1016/s0378-3774(02)00155-5 
  71. 71.Fernandez-Aunión C, Hamouda TB, Iglesias-Guerra F, Argandoña M, Reina-Bueno M, et al. (2010) Biosynthesis of compatible solutes in rhizobial strains isolated from Phaseolus vulgaris nodules in Tunisian fields. BMC Microbiol 10: 192. doi: 10.1186/1471-2180-10-192. pmid:20633304 
  72. 72.Vriezen JAC, de Bruijn FJ, Nüsslein K (2007) Responses of rhizobia to desiccation in relation to osmotic stress, oxygen, and temperature. Appl Environ Microbiol 73: 3451–3459. pmid:17400779 doi: 10.1128/aem.02991-06 
  73. 73.Arora NK, Singhal V, Maheshwari DK (2006) Salinity-induced accumulation of poly-β-hydroxybutyrate in rhizobia indicating its role in cell protection. World J Microbiol Biotechnol 22: 603–606. doi: 10.1007/s11274-005-9077-1 
  74. 74.Lu YL, Chen WF, Wang ET, Han LL, Zhang XX, et al. (2009) Mesorhizobium shangrilense sp. nov., isolated from root nodules of Caragana species. Int J Syst Evol Microbiol 59: 3012–3018. doi: 10.1099/ijs.0.007393-0. pmid:19643905 
  75. 75.Wade TK, Le Quéré A, Laguerre G, N’zoué A, Ndione J-A, et al. (2014) Eco-geographical diversity of cowpea bradyrhizobia in Senegal is marked by dominance of two genetic types. Syst Appl Microbiol 37: 129–139. doi: 10.1016/j.syapm.2013.10.002. pmid:24373721 

For further details log on website :
http://journals.plos.org/plosone/article?id=10.1371/journal.pone.0117667

No comments:

Post a Comment