Genetic variation and genetic complexity of nodule occupancy in soybean inoculated with USDA110 and USDA123 rhizobium strains.

BACKGROUND: Symbiotic nitrogen fixation differs among Bradyrhizobium japonicum strains. Soybean inoculated with USDA123 has a lower yield than strains known to have high nitrogen fixation efficiency, such as USDA110. In the main soybean-producing area in the Midwest of the United States, USDA123 has a high nodule incidence in field-grown soybean and is competitive but inefficient in nitrogen fixation. In this study, a high-throughput system was developed to characterize nodule number among 1,321 Glycine max and 69 Glycine soja accessions single inoculated with USDA110 and USDA123. RESULTS: Seventy-three G. max accessions with significantly different nodule number of USDA110 and USDA123 were identified. After double inoculating 35 of the 73 accessions, it was observed that PI189939, PI317335, PI324187B, PI548461, PI562373, and PI628961 were occupied by USDA110 and double-strain nodules but not by USDA123 nodules alone. PI567624 was only occupied by USDA110 nodules, and PI507429 restricted all strains. Analysis showed that 35 loci were associated with nodule number in G. max when inoculated with strain USDA110 and 35 loci with USDA123. Twenty-three loci were identified in G. soja when inoculated with strain USDA110 and 34 with USDA123. Only four loci were common across two treatments, and each locus could only explain 0.8 to 1.5% of phenotypic variation. CONCLUSIONS: High-throughput phenotyping systems to characterize nodule number and occupancy were developed, and soybean germplasm restricting rhizobium strain USDA123 but preferring USDA110 was identified. The larger number of minor effects and a small few common loci controlling the nodule number indicated trait genetic complexity and strain-dependent nodulation restriction. The information from the present study will add to the development of cultivars that limit USDA123, thereby increasing nitrogen fixation efficiency and productivity.

The Genome of the Acid Soil-Adapted Strain Rhizobium favelukesii OR191 Encodes Determinants for Effective Symbiotic Interaction With Both an Inverted Repeat Lacking Clade and a Phaseoloid Legume Host.

Although Medicago sativa forms highly effective symbioses with the comparatively acid-sensitive genus Ensifer, its introduction into acid soils appears to have selected for symbiotic interactions with acid-tolerant R. favelukesii strains. Rhizobium favelukesii has the unusual ability of being able to nodulate and fix nitrogen, albeit sub-optimally, not only with M. sativa but also with the promiscuous host Phaseolus vulgaris. Here we describe the genome of R. favelukesii OR191 and genomic features important for the symbiotic interaction with both of these hosts. The OR191 draft genome contained acid adaptation loci, including the highly acid-inducible lpiA/acvB operon and olsC, required for production of lysine- and ornithine-containing membrane lipids, respectively. The olsC gene was also present in other acid-tolerant Rhizobium strains but absent from the more acid-sensitive Ensifer microsymbionts. The OR191 symbiotic genes were in general more closely related to those found in Medicago microsymbionts. OR191 contained the nodA, nodEF, nodHPQ, and nodL genes for synthesis of polyunsaturated, sulfated and acetylated Nod factors that are important for symbiosis with Medicago, but contained a truncated nodG, which may decrease nodulation efficiency with M. sativa. OR191 contained an E. meliloti type BacA, which has been shown to specifically protect Ensifer microsymbionts from Medicago nodule-specific cysteine-rich peptides. The nitrogen fixation genes nifQWZS were present in OR191 and P. vulgaris microsymbionts but absent from E. meliloti-Medicago microsymbionts. The ability of OR191 to nodulate and fix nitrogen symbiotically with P. vulgaris indicates that this host has less stringent requirements for nodulation than M. sativa but may need rhizobial strains that possess nifQWZS for N2-fixation to occur. OR191 possessed the exo genes required for the biosynthesis of succinoglycan, which is required for the Ensifer-Medicago symbiosis. However, 1H-NMR spectra revealed that, in the conditions tested, OR191 exopolysaccharide did not contain a succinyl substituent but instead contained a 3-hydroxybutyrate moiety, which may affect its symbiotic performance with Medicago hosts. These findings provide a foundation for the genetic basis of nodulation requirements and symbiotic effectiveness with different hosts.

Genetic Diversity and Phylogenetic Relationships of Annual and Perennial Glycine Species.

We have estimated the average genetic diversity of two Glycine annual and six perennial species based upon 76 orthologous gene sets and performed phylogenetic analysis, divergence analysis and tests for departure from neutrality of the eight species using 52 orthologous gene sets. In addition, 367 orthologous gene sets were used to estimate the relationships of 11 G. canescens accessions. Among the perennials, G. canescens showed the highest nucleotide diversity. The other perennials, except for G. tomentella, had higher nucleotide diversity than the two annuals. Phylogenetic analysis of the Glycine showed a similar genome grouping with the previous report except for G. cyrtoloba and G. stenophita which formed a sister clade in the study. Divergence analysis supported the phylogenetic relationships that G. falcata was the most divergent from G. max, followed by G. cyrtoloba, G. syndetika, G. tomentella D3, G. stenophita and G. canescens Most genic sequences were homogeneous in the levels of polymorphism and divergence between G. max and other Glycine species based on the HKA test, thus, Glycine perennials may have experienced a very similar evolution as inferred by trans-specific mutation analysis. The greater genetic diversity of most perennial Glycine species and their origins from the warmer and drier climates of Australia suggests the perennials maybe a potential source of heat and drought resistance that will be of value in the face of climate change.

A Proteomic Network for Symbiotic Nitrogen Fixation Efficiency in Bradyrhizobium elkanii.

Rhizobia colonize legumes and reduce N2 to NH3 in root nodules. The current model is that symbiotic rhizobia bacteroids avoid assimilating this NH3. Instead, host legume cells form glutamine from NH3, and the nitrogen is returned to the bacteroid as dicarboxylates, peptides, and amino acids. In soybean cells surrounding bacteroids, glutamine also is converted to ureides. One problem for soybean cultivation is inefficiency in symbiotic N2 fixation, the biochemical basis of which is unknown. Here, the proteomes of bacteroids of Bradyrhizobium elkanii USDA76 isolated from N2 fixation-efficient Peking and -inefficient Williams 82 soybean nodules were analyzed by mass spectrometry. Nearly half of the encoded bacterial proteins were quantified. Efficient bacteroids produced greater amounts of enzymes to form Nod factors and had increased amounts of signaling proteins, transporters, and enzymes needed to generate ATP to power nitrogenase and to acquire resources. Parallel investigation of nodule proteins revealed that Peking had no significantly greater accumulation of enzymes needed to assimilate NH3 than Williams 82. Instead, efficient bacteroids had increased amounts of enzymes to produce amino acids, including glutamine, and to form ureide precursors. These results support a model for efficient symbiotic N2 fixation in soybean where the bacteroid assimilates NH3 for itself.

High-quality draft genome sequence of Ensifer meliloti Mlalz-1, a microsymbiont of Medicago laciniata (L.) miller collected in Lanzarote, Canary Islands, Spain.

10.1601/nm.1335 Mlalz-1 (INSDC = ATZD00000000) is an aerobic, motile, Gram-negative, non-spore-forming rod that was isolated from an effective nitrogen-fixing nodule of Medicago laciniata (L.) Miller from a soil sample collected near the town of Guatiza on the island of Lanzarote, the Canary Islands, Spain. This strain nodulates and forms an effective symbiosis with the highly specific host M. laciniata. This rhizobial genome was sequenced as part of the DOE Joint Genome Institute 2010 Genomic Encyclopedia for Bacteria and Archaea-Root Nodule Bacteria (GEBA-RNB) sequencing project. Here the features of 10.1601/nm.1335 Mlalz-1 are described, together with high-quality permanent draft genome sequence information and annotation. The 6,664,116 bp high-quality draft genome is arranged in 99 scaffolds of 100 contigs, containing 6314 protein-coding genes and 74 RNA-only encoding genes. Strain Mlalz-1 is closely related to 10.1601/nm.1335 10.1601/strainfinder?urlappend=%3Fid%3DIAM+12611 T, 10.1601/nm.1334 A 321T and 10.1601/nm.17831 10.1601/strainfinder?urlappend=%3Fid%3DORS+1407 T, based on 16S rRNA gene sequences. gANI values of ≥98.1% support the classification of strain Mlalz-1 as 10.1601/nm.1335. Nodulation of M. laciniata requires a specific nodC allele, and the nodC gene of strain Mlalz-1 shares ≥98% sequence identity with nodC of M. laciniata-nodulating 10.1601/nm.1328 strains, but ≤93% with nodC of 10.1601/nm.1328 strains that nodulate other Medicago species. Strain Mlalz-1 is unique among sequenced 10.1601/nm.1335 strains in possessing genes encoding components of a T2SS and in having two versions of the adaptive acid tolerance response lpiA-acvB operon. In 10.1601/nm.1334 strain 10.1601/strainfinder?urlappend=%3Fid%3DWSM+419, lpiA is essential for enhancing survival in lethal acid conditions. The second copy of the lpiA-acvB operon of strain Mlalz-1 has highest sequence identity (> 96%) with that of 10.1601/nm.1334 strains, which suggests genetic recombination between strain Mlalz-1 and 10.1601/nm.1334 and the horizontal gene transfer of lpiA-acvB.

High-quality permanent draft genome sequence of the Bradyrhizobium elkanii type strain USDA 76T, isolated from Glycine max (L.) Merr.

Bradyrhizobium elkanii USDA 76T (INSCD = ARAG00000000), the type strain for Bradyrhizobium elkanii, is an aerobic, motile, Gram-negative, non-spore-forming rod that was isolated from an effective nitrogen-fixing root nodule of Glycine max (L. Merr) grown in the USA. Because of its significance as a microsymbiont of this economically important legume, B. elkanii USDA 76T was selected as part of the DOE Joint Genome Institute 2010 Genomic Encyclopedia for Bacteria and Archaea-Root Nodule Bacteria sequencing project. Here the symbiotic abilities of B. elkanii USDA 76T are described, together with its genome sequence information and annotation. The 9,484,767 bp high-quality draft genome is arranged in 2 scaffolds of 25 contigs, containing 9060 protein-coding genes and 91 RNA-only encoding genes. The B. elkanii USDA 76T genome contains a low GC content region with symbiotic nod and fix genes, indicating the presence of a symbiotic island integration. A comparison of five B. elkanii genomes that formed a clique revealed that 356 of the 9060 protein coding genes of USDA 76T were unique, including 22 genes of an intact resident prophage. A conserved set of 7556 genes were also identified for this species, including genes encoding a general secretion pathway as well as type II, III, IV and VI secretion system proteins. The type III secretion system has previously been characterized as a host determinant for Rj and/or rj soybean cultivars. Here we show that the USDA 76T genome contains genes encoding all the type III secretion system components, including a translocon complex protein NopX required for the introduction of effector proteins into host cells. While many bradyrhizobial strains are unable to nodulate the soybean cultivar Clark (rj1), USDA 76T was able to elicit nodules on Clark (rj1), although in reduced numbers, when plants were grown in Leonard jars containing sand or vermiculite. In these conditions, we postulate that the presence of NopX allows USDA 76T to introduce various effector molecules into this host to enable nodulation.

Biogeography of a Novel Ensifer meliloti Clade Associated with the Australian Legume Trigonella suavissima.

Here, we describe a novel clade within Ensifer meliloti and consider how geographic and ecological isolation contributed to the limited distribution of this group. Members of the genus Ensifer are best known for their ability to form nitrogen-fixing symbioses with forage legumes of three related genera, Medicago L., Melilotus Mill., and Trigonella L., which are members of the tribe Trifolieae. These legumes have a natural distribution extending from the Mediterranean Basin through western Asia, where there is an unsurpassed number of species belonging to these genera. Trigonella suavissima L. is unusual in that it is the only species in the tribe Trifolieae that is native to Australia. We compared the genetic diversity and taxonomic placement of rhizobia nodulating T. suavissima with those of members of an Ensifer reference collection. Our goal was to determine if the T. suavissima rhizobial strains, like their plant host, are naturally limited to the Australian continent. We used multilocus sequence analysis to estimate the genetic relatedness of 56 T. suavissima symbionts to 28 Ensifer reference strains. Sequence data were partitioned according to the replicons in which the loci are located. The results were used to construct replicon-specific phylogenetic trees. In both the chromosomal and chromid trees, the Australian strains formed a distinct clade within E. meliloti The strains also shared few alleles with Ensifer reference strains from other continents. Carbon source utilization assays revealed that the strains are also unusual in their ability to utilize 2-oxoglutarate as a sole carbon source. A strategy was outlined for locating similar strains elsewhere.IMPORTANCE In this study, we employed a biogeographical approach to investigate the origins of a symbiotic relationship between an Australian legume and its nitrogen-fixing rhizobia. The question of the ancestral origins of these symbionts is based on the observation that the legume host is not closely related to other native Australian legumes. Previous research has shown that the legume host Trigonella suavissima is instead closely related to legumes native to the Mediterranean Basin and western Asia, suggesting that it may have been introduced in Australia from those regions. This led to the question of whether its rhizobia may have been introduced as well. In this study, we were unable to find persuasive evidence supporting this hypothesis. Instead, our results suggest either that the T. suavissima rhizobia are native to Australia or that our methods for locating their close relatives elsewhere are inadequate. A strategy to investigate the latter alternative is proposed.

Development and application of a multilocus sequence analysis method for the identification of genotypes within genus Bradyrhizobium and for establishing nodule occupancy of soybean (Glycine max L. Merr).

A multilocus sequence typing (MLST) method based on allelic variation of seven chromosomal loci was developed for characterizing genotypes (GT) within the genus Bradyrhizobium. With the method, 29 distinct multilocus GT were identified among 190 culture collection soybean strains. The occupancy of 347 nodules taken from uninoculated field-grown soybean plants also was determined. The bacteroid GT were either the same as or were closely related to GT identified among strains in the culture collection. Double-nodule occupancy estimates of 2.9% were much lower than values published based on serology. Of the 347 nodules examined, 337 and 10 were occupied by Bradyrhizobium japonicum and B. elkanii, respectively. The collection strains within the species B. japonicum and B. elkaniialso were compared with Bradyrhizobium cultures from other legumes. In many cases, the observed GT varied more according to their geographic origin than by their trap hosts of isolation. In other cases, there were no apparent relationships with either the legume or geographic source. The MLST method that was developed should be a useful tool in determining the influence of geographic location, temperature, season, soil type, and host plant cultivar on the distribution of GT of Bradyrhizobium spp.

Application of multilocus sequence typing to study the genetic structure of megaplasmids in medicago-nodulating rhizobia.

A multilocus sequence typing (MLST) analysis was used to examine the genetic structure and diversity within the two large extrachromosomal replicons in Medicago-nodulating rhizobia (Sinorhizobium meliloti and Sinorhizobium medicae). The allelic diversity within these replicons was high compared to the reported diversity within the corresponding chromosomes of the same strains (P. van Berkum et al., J. Bacteriol. 188:5570-5577, 2006). Also, there was strong localized linkage disequilibrium (LD) between certain pSymA loci: e.g., nodC and nifD. Although both of these observations could be explained by positive (or diversifying) selection by plant hosts, results of tests for positive selection did not provide consistent support for this hypothesis. The strong LD observed between the nodC and nifD genes could also be explained by their close proximity on the pSymA replicon. Evidence was obtained that some nodC alleles had a history of intragenic recombination, while other alleles of this locus had a history of intergenic recombination. Both types of recombination were associated with a decline in symbiotic competence with Medicago sativa as the host plant. The combined observations of LD between the nodC and nifD genes and intragenic recombination within one of these loci indicate that the symbiotic gene region on the pSymA plasmid has evolved as a clonal segment, which has been laterally transferred within the natural populations.

Chromosomal and symbiotic relationships of rhizobia nodulating Medicago truncatula and M. laciniata.

Multilocus sequence typing (MLST) is a sequence-based method used to characterize bacterial genomes. This method was used to examine the genetic structure of Medicago-nodulating rhizobia at the Amra site, which is located in an arid region of Tunisia. Here the annual medics Medicago laciniata and M. truncatula are part of the natural flora. The goal of this study was to identify whether distinct chromosomal groups of rhizobia nodulate M. laciniata because of its restricted requirement for specific rhizobia. The MLST analysis involved determination of sequence variation in 10 chromosomal loci of 74 isolates each of M. laciniata and M. truncatula. M. truncatula was used as a control trap host, because unlike M. laciniata, it has relatively unrestricted rhizobial requirements. Allelic diversity among the plasmid nodC alleles in the isolates was also determined. The 148 isolates were placed into 26 chromosomal sequence types (STs), only 3 of which had been identified previously. The rhizobia of M. laciniata were shown to be part of the general Medicago-nodulating population in the soil because 99.95% of the isolates had chromosomal genotypes similar to those recovered from M. truncatula. However, the isolates recovered from M. laciniata were less diverse than those recovered from M. truncatula, and they also harbored an unusual nodC allele. This could perhaps be best explained by horizontal transfer of the different nodC alleles among members of the Medicago-nodulating rhizobial population at the field site. Evidence indicating a history of lateral transfer of rhizobial symbiotic genes across distinct chromosomal backgrounds is provided.

Multilocus sequence typing as an approach for population analysis of Medicago-nodulating rhizobia.

Multilocus sequence typing (MLST), a sequence-based method to characterize bacterial genomes, was used to examine the genetic structure in a large collection of Medicago-nodulating rhizobial strains. This is the first study where MLST has been applied in conjunction with eBURST analysis to determine the population genetic structure of nonpathogenic bacteria recovered from the soil environment. Sequence variation was determined in 10 chromosomal loci of 231 strains that predominantly originated from southwest Asia. Genetic diversity for each locus ranged from 0.351 to 0.819, and the strains examined were allocated to 91 different allelic profiles or sequence types (STs). The genus Medicago is nodulated by at least two groups of rhizobia with divergent chromosomes that have been classified as Sinorhizobium meliloti and Sinorhizobium medicae. Evidence was obtained that the degree of genetic exchange among the chromosomes across these groups is limited. The symbiosis with Medicago polymorpha of nine strains placed in one of these groups, previously identified as S. medicae, ranged from ineffective to fully effective, indicating that there was no strong relationship between symbiotic phenotype and chromosomal genotype.