A Putative Type III Secretion System Effector Encoded by the MA20_12780 Gene in Bradyrhizobium japonicum Is-34 Causes Incompatibility with Rj4 Genotype Soybeans

Comparative genomic analysis of Bradyrhizobium strains with natural variability in the efficiency of nitrogen fixation, competitiveness, and adaptation to stressful edaphoclimatic conditions

Bradyrhizobium is known for fixing atmospheric nitrogen in symbiosis with agronomically important crops. This study focused on two groups of strains, each containing eight natural variants of the parental strains, Bradyrhizobium japonicum SEMIA 586 (=CNPSo 17) or Bradyrhizobium diazoefficiens SEMIA 566 (=CNPSo 10). CNPSo 17 and CNPSo 10 were used as commercial inoculants for soybean crops in Brazil at the beginning of the crop expansion in the southern region in the 1960s-1970s. Variants derived from these parental strains were obtained in the late 1980s through a strain selection program aimed at identifying elite strains adapted to a new cropping frontier in the central-western Cerrado region, with a higher capacity of biological nitrogen fixation (BNF) and competitiveness. Here, we aimed to detect genetic variations possibly related to BNF, competitiveness for nodule occupancy, and adaptation to the stressful conditions of the Brazilian Cerrado soils. High-quality genome assemblies were produced for all strains. The core genome phylogeny revealed that strains of each group are closely related, as confirmed by high average nucleotide identity values. However, variants accumulated divergences resulting from horizontal gene transfer, genomic rearrangements, and nucleotide polymorphisms. The B. japonicum group presented a larger pangenome and a higher number of nucleotide polymorphisms than the B. diazoefficiens group, possibly due to its longer adaptation time to the Cerrado soil. Interestingly, five strains of the B. japonicum group carry two plasmids. The genetic variability found in both groups is discussed considering the observed differences in their BNF capacity, competitiveness for nodule occupancy, and environmental adaptation.IMPORTANCEToday, Brazil is a global leader in the study and use of biological nitrogen fixation with soybean crops. As Brazilian soils are naturally void of soybean-compatible bradyrhizobia, strain selection programs were established, starting with foreign isolates. Selection searched for adaptation to the local edaphoclimatic conditions, higher efficiency of nitrogen fixation, and strong competitiveness for nodule occupancy. We analyzed the genomes of two parental strains of Bradyrhizobium japonicum and Bradyrhizobium diazoefficiens and eight variant strains derived from each parental strain. We detected two plasmids in five strains and several genetic differences that might be related to adaptation to the stressful conditions of the soils of the Brazilian Cerrado biome. We also detected genetic variations in specific regions that may impact symbiotic nitrogen fixation. Our analysis contributes to new insights into the evolution of Bradyrhizobium, and some of the identified differences may be applied as genetic markers to assist strain selection programs.

Broad-spectrum ubiquitin/ubiquitin-like deconjugation activity of the rhizobial effector NopD from Bradyrhizobium (sp. XS1150)

The post-translational modification of proteins by ubiquitin-like modifiers (UbLs), such as SUMO, ubiquitin, and Nedd8, regulates a vast array of cellular processes. Dedicated UbL deconjugating proteases families reverse these modifications. During bacterial infection, effector proteins, including deconjugating proteases, are released to disrupt host cell defenses and promote bacterial survival. NopD, an effector protein from rhizobia involved in legume nodule symbiosis, exhibits deSUMOylation activity and, unexpectedly, also deubiquitination and deNeddylation activities. Here, we present two crystal structures of Bradyrhizobium (sp. XS1150) NopD complexed with either Arabidopsis SUMO2 or ubiquitin at 1.50 Å and 1.94 Å resolution, respectively. Despite their low sequence similarity, SUMO and ubiquitin bind to a similar NopD interface, employing a unique loop insertion in the NopD sequence. In vitro binding and activity assays reveal specific residues that distinguish between deubiquitination and deSUMOylation. These unique multifaceted deconjugating activities against SUMO, ubiquitin, and Nedd8 exemplify an optimized bacterial protease that disrupts distinct UbL post-translational modifications during host cell infection.

The Rhizobial Type 3 Secretion System: The Dr. Jekyll and Mr. Hyde in the Rhizobium–Legume Symbiosis

Rhizobia are soil bacteria that can establish a symbiotic association with legumes. As a result, plant nodules are formed on the roots of the host plants where rhizobia differentiate to bacteroids capable of fixing atmospheric nitrogen into ammonia. This ammonia is transferred to the plant in exchange of a carbon source and an appropriate environment for bacterial survival. This process is subjected to a tight regulation with several checkpoints to allow the progression of the infection or its restriction. The type 3 secretion system (T3SS) is a secretory system that injects proteins, called effectors (T3E), directly into the cytoplasm of the host cell, altering host pathways or suppressing host defense responses. This secretion system is not present in all rhizobia but its role in symbiosis is crucial for some symbiotic associations, showing two possible faces as Dr. Jekyll and Mr. Hyde: it can be completely necessary for the formation of nodules, or it can block nodulation in different legume species/cultivars. In this review, we compile all the information currently available about the effects of different rhizobial effectors on plant symbiotic phenotypes. These phenotypes are diverse and highlight the importance of the T3SS in certain rhizobium–legume symbioses.

Reference nodule transcriptomes for Melilotus officinalis and Medicago sativa cv. Algonquin

Host/symbiont compatibility is a hallmark of the symbiotic nitrogen-fixing interaction between rhizobia and legumes, mediated in part by plant-produced nodule-specific cysteine-rich (NCR) peptides and the bacterial BacA membrane protein that can act as a NCR peptide transporter. In addition, the genetic and metabolic properties supporting symbiotic nitrogen fixation often differ between compatible partners, including those sharing a common partner, highlighting the need for multiple study systems. Here, we report high-quality nodule transcriptome assemblies for Medicago sativa cv. Algonquin and Melilotus officinalis, two legumes able to form compatible symbioses with Sinorhizobium meliloti. The compressed M. sativa and M. officinalis assemblies consisted of 79,978 and 64,593 contigs, respectively, of which 33,341 and 28,278 were assigned putative annotations, respectively. As expected, the two transcriptomes showed broad similarity at a global level. We were particularly interested in the NCR peptide profiles of these plants, as these peptides drive bacterial differentiation during the symbiosis. A total of 412 and 308 NCR peptides were predicted from the M. sativa and M. officinalis transcriptomes, respectively, with approximately 9% of the transcriptome of both species consisting of NCR transcripts. Notably, transcripts encoding highly cationic NCR peptides (isoelectric point > 9.5), which are known to have antimicrobial properties, were ∼2-fold more abundant in M. sativa than in M. officinalis, and ∼27-fold more abundant when considering only NCR peptides in the six-cysteine class. We hypothesize that the difference in abundance of highly cationic NCR peptides explains our previous observation that some rhizobial bacA alleles which can support symbiosis with M. officinalis are unable to support symbiosis with M. sativa.

Varietas Delectat: Exploring Natural Variations in Nitrogen-Fixing Symbiosis Research

The nitrogen-fixing symbiosis between leguminous plants and soil bacteria collectively called rhizobia plays an important role in the global nitrogen cycle and is an essential component of sustainable agriculture. Genetic determinants directing the development and functioning of the interaction have been identified with the help of a very limited number of model plants and bacterial strains. Most of the information obtained from the study of model systems could be validated on crop plants and their partners. The investigation of soybean cultivars and different rhizobia, however, has revealed the existence of ineffective interactions between otherwise effective partners that resemble gene-for-gene interactions described for pathogenic systems. Since then, incompatible interactions between natural isolates of model plants, called ecotypes, and different bacterial partner strains have been reported. Moreover, diverse phenotypes of both bacterial mutants on different host plants and plant mutants with different bacterial strains have been described. Identification of the genetic factors behind the phenotypic differences did already and will reveal novel functions of known genes/proteins, the role of certain proteins in some interactions, and the fine regulation of the steps during nodule development.

The Type III Effectome of the Symbiotic Bradyrhizobium vignae Strain ORS3257

Many Bradyrhizobium strains are able to establish a Nod factor-independent symbiosis with the leguminous plant Aeschynomene indica by the use of a type III secretion system (T3SS). Recently, an important advance in the understanding of the molecular factors supporting this symbiosis has been achieved by the in silico identification and functional characterization of 27 putative T3SS effectors (T3Es) of Bradyrhizobium vignae ORS3257. In the present study, we experimentally extend this catalog of T3Es by using a multi-omics approach. Transcriptome analysis under non-inducing and inducing conditions in the ORS3257 wild-type strain and the ttsI mutant revealed that the expression of 18 out of the 27 putative effectors previously identified, is under the control of TtsI, the global transcriptional regulator of T3SS and T3Es. Quantitative shotgun proteome analysis of culture supernatant in the wild type and T3SS mutant strains confirmed that 15 of the previously determined candidate T3Es are secreted by the T3SS. Moreover, the combined approaches identified nine additional putative T3Es and one of them was experimentally validated as a novel effector. Our study underscores the power of combined proteome and transcriptome analyses to complement in silico predictions and produce nearly complete effector catalogs. The establishment of the ORS3257 effectome will form the basis for a full appraisal of the symbiotic properties of this strain during its interaction with various host legumes via different processes.

One door closes, another opens: when nodulation impairment with natural hosts extends rhizobial host-range

The rhizobium-legume symbiosis is the best-understood plant-microbe association. The high degree of specificity observed in this relationship is supported by a complex exchange of signals between the two components of the symbiosis. Findings reported in last years indicate that multiple molecular mechanisms, such as the production of a particular set of nodulation factors at a very specific concentration or a suitable arsenal of effectors secreted through the type III secretion system, have been adjusted during evolution to ensure and optimize the recognition of specific rhizobial strains by its legume host. Qualitative or quantitative changes in the production of these symbiotic molecular determinants are detrimental for nodulation with its natural host but, in some cases, can also result beneficial for the rhizobium since it extends the nodulation host-range to other legumes. Potential repercussion of the extension in the nodulation host-range of rhizobia is discussed.

Determinants of Host Range Specificity in Legume-Rhizobia Symbiosis

Leguminous plants possess the almost unique ability to enter symbiosis with soil-resident, nitrogen fixing bacteria called rhizobia. During this symbiosis, the bacteria physically colonize specialized organs on the roots of the host plant called nodules, where they reduce atmospheric nitrogen into forms that can be assimilated by the host plant and receive photosynthates in return. In order for nodule development to occur, there is extensive chemical cross-talk between both parties during the formative stages of the symbiosis. The vast majority of the legume family are capable of forming root nodules and typically rhizobia are only able to fix nitrogen within the context of this symbiotic association. However, many legume species only enter productive symbiosis with a few, or even single rhizobial species or strains, and vice-versa. Permitting symbiosis with only rhizobial strains that will be able to fix nitrogen with high efficiency is a crucial strategy for the host plant to prevent cheating by rhizobia. This selectivity is enforced at all stages of the symbiosis, with partner choice beginning during the initial communication between the plant and rhizobia. However, it can also be influenced even once nitrogen-fixing nodules have developed on the root. This review sets out current knowledge about the molecular mechanisms employed by both parties to influence host range during legume-rhizobia symbiosis.

NopD of Bradyrhizobium sp. XS1150 Possesses SUMO Protease Activity

Effectors secreted by the type III protein secretion system (T3SS) of rhizobia are host-specific determinants of the nodule symbiosis. Here, we have characterized NopD, a putative type III effector of Bradyrhizobium sp. XS1150. NopD was found to possess a functional N-terminal secretion signal sequence that could replace that of the NopL effector secreted by Sinorhizobium sp. NGR234. Recombinant NopD and the C-terminal domain of NopD alone can process small ubiquitin-related modifier (SUMO) proteins and cleave SUMO-conjugated proteins. Activity was abolished in a NopD variant with a cysteine-to-alanine substitution in the catalytic core (NopD-C₉₇₂A). NopD recognizes specific plant SUMO proteins (AtSUMO1 and AtSUMO2 of Arabidopsis thaliana; GmSUMO of Glycine max; PvSUMO of Phaseolus vulgaris). Subcellular localization analysis with A. thaliana protoplasts showed that NopD accumulates in nuclear bodies. NopD, but not NopD-C₉₇₂A, induces cell death when expressed in Nicotiana tabacum. Likewise, inoculation tests with constructed mutant strains of XS1150 indicated that nodulation of Tephrosia vogelii is negatively affected by the protease activity of NopD. In conclusion, our findings show that NopD is a symbiosis-related protein that can process specific SUMO proteins and desumoylate SUMO-conjugated proteins.

Mutants Disrupted in the Type III Secretion System of Bradyrhizobium elkanii BLY3-8 Overcame Nodulation Restriction by Rj3-genotype Soybean

Bradyrhizobium elkanii BLY3-8 does not form nodules on the roots of Rj3-genotype soybean (cultivar D-51). This is a cultivar-specific nodulation restriction. The genes A6X20_40975 and A6X20_41030 in strain BLY3-8 were predicted to encode the transcriptional activator and apparatus of the type III secretion system (T3SS) (the proteins TtsI and RhcJ), respectively. Mutants disrupted in these genes overcame the nodulation restriction. These results suggest that an effector injected via T3SS into Rj3-genotype soybean is involved in nodulation restriction by Rj3-genotype soybean.

Legume Nodules: Massive Infection in the Absence of Defense Induction

Plants of the legume family host massive intracellular bacterial populations in the tissues of specialized organs, the nodules. In these organs, the bacteria, named rhizobia, can fix atmospheric nitrogen and transfer it to the plant. This special metabolic skill provides to the legumes an advantage when they grow on nitrogen-scarce substrates. While packed with rhizobia, the nodule cells remain alive, metabolically active, and do not develop defense reactions. Here, we review our knowledge on the control of plant immunity during the rhizobia-legume symbiosis. We present the results of an evolutionary process that selected both divergence of microbial-associated molecular motifs and active suppressors of immunity on the rhizobial side and, on the legume side, active mechanisms that contribute to suppression of immunity.

Genetic and Molecular Mechanisms Underlying Symbiotic Specificity in Legume-Rhizobium Interactions

Legumes are able to form a symbiotic relationship with nitrogen-fixing soil bacteria called rhizobia. The result of this symbiosis is to form nodules on the plant root, within which the bacteria can convert atmospheric nitrogen into ammonia that can be used by the plant. Establishment of a successful symbiosis requires the two symbiotic partners to be compatible with each other throughout the process of symbiotic development. However, incompatibility frequently occurs, such that a bacterial strain is unable to nodulate a particular host plant or forms nodules that are incapable of fixing nitrogen. Genetic and molecular mechanisms that regulate symbiotic specificity are diverse, involving a wide range of host and bacterial genes/signals with various modes of action. In this review, we will provide an update on our current knowledge of how the recognition specificity has evolved in the context of symbiosis signaling and plant immunity.

Adaptive evolution of rhizobial symbiotic compatibility mediated by co-evolved insertion sequences

Mutualism between bacteria and eukaryotes has essential roles in the history of life, but the evolution of their compatibility is poorly understood. Here we show that different Sinorhizobium strains can form either nitrogen-fixing nodules or uninfected pseudonodules on certain cultivated soybeans, while being all effective microsymbionts of some wild soybeans. However, a few well-infected nodules can be found on a commercial soybean using inocula containing a mixed pool of Tn5 insertion mutants derived from an incompatible strain. Reverse genetics and genome sequencing of compatible mutants demonstrated that inactivation of T3SS (type three secretion system) accounted for this phenotypic change. These mutations in the T3SS gene cluster were dominated by parallel transpositions of insertion sequences (ISs) other than the introduced Tn5. This genetic and phenotypic change can also be achieved in an experimental evolution scenario on a laboratory time scale using incompatible wild-type strains as inocula. The ISs acting in the adaptive evolution of Sinorhizobium strains exhibit broader phyletic and replicon distributions than other ISs, and prefer target sequences of low GC% content, a characteristic feature of symbiosis plasmid where T3SS genes are located. These findings suggest an important role of co-evolved ISs in the adaptive evolution of rhizobial compatibility.

Type 3 Secretion System (T3SS) of Bradyrhizobium sp. DOA9 and Its Roles in Legume Symbiosis and Rice Endophytic Association

The Bradyrhizobium sp. DOA9 strain isolated from a paddy field has the ability to nodulate a wide spectrum of legumes. Unlike other bradyrhizobia, this strain has a symbiotic plasmid harboring nod, nif, and type 3 secretion system (T3SS) genes. This T3SS cluster contains all the genes necessary for the formation of the secretory apparatus and the transcriptional activator (TtsI), which is preceded by a nod-box motif. An in silico search predicted 14 effectors putatively translocated by this T3SS machinery. In this study, we explored the role of the T3SS in the symbiotic performance of DOA9 by evaluating the ability of a T3SS mutant (ΩrhcN) to nodulate legumes belonging to Dalbergioid, Millettioid, and Genistoid tribes. Among the nine species tested, four (Arachis hypogea, Vigna radiata, Crotalaria juncea, and Macroptilium atropurpureum) responded positively to the rhcN mutation (ranging from suppression of plant defense reactions, an increase in the number of nodules and a dramatic improvement in nodule development and infection), one (Stylosanthes hamata) responded negatively (fewer nodules and less nitrogen fixation) and four species (Aeschynomene americana, Aeschynomene afraspera, Indigofera tinctoria, and Desmodium tortuosum) displayed no phenotype. We also tested the role of the T3SS in the ability of the DOA9 strain to endophytically colonize rice roots, but detected no effect of the T3SS mutation, in contrast to what was previously reported in the Bradyrhizobium SUTN9-2 strain. Taken together, these data indicate that DOA9 contains a functional T3SS that interferes with the ability of the strain to interact symbiotically with legumes but not with rice.

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

Bradyrhizobium elkanii USDA 76^T (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 76^T 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 76^T 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 76^T 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 76^T 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 76^T 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 76^T 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 76^T to introduce various effector molecules into this host to enable nodulation.

The Soybean Rfg1 Gene Restricts Nodulation by Sinorhizobium fredii USDA193

Sinorhizobium fredii is a fast-growing rhizobial species that can establish a nitrogen-fixing symbiosis with a wide range of legume species including soybeans (Glycine max). In soybeans, this interaction shows a high level of specificity such that particular S. fredii strains nodulate only a limited set of plant genotypes. Here we report the identification of a dominant gene in soybeans that restricts nodulation with S. fredii USDA193. Genetic mapping in an F2 population revealed co-segregation of the underlying locus with the previously cloned Rfg1 gene. The Rfg1 allele encodes a member of the Toll-interleukin receptor/nucleotide-binding site/leucine-rich repeat class of plant resistance proteins that restricts nodulation by S. fredii strains USDA257 and USDA205, and an allelic variant of this gene also restricts nodulation by Bradyrhizobium japonicum USDA122. By means of complementation tests and CRISPR/Cas9-mediated gene knockouts, we demonstrate that the Rfg1 allele also is responsible for resistance to nodulation by S. fredii USDA193. Therefore, the Rfg1 allele likely provides broad-spectrum resistance to nodulation by many S. fredii and B. japonicum strains in soybeans.