Tomasz Stępkowski

Phylogenetic analyses of symbiotic nodulation genes support vertical and lateral gene co-transfer within the Bradyrhizobium genus

Symbiotic nitrogen fixing bacteria-known as rhizobia-harbour a set of nodulation (nod) genes that control the synthesis of modified lipo-chitooligosaccharides, called Nod factors that are required for legume nodulation. The nodA gene, which is essential for symbiosis, is responsible for the attachment of the fatty acid group to the oligosaccharide backbone. The nodZ, nolL, and noeI genes are involved in specific modifications of Nod factors common to bradyrhizobia, i.e., the transfer of a fucosyl group on the Nod factor core, fucose acetylation and fucose methylation, respectively. PCR amplification, sequencing and phylogenetic analysis of nodA gene sequences from a collection of diverse Bradyrhizobium strains revealed the monophyletic character with the possible exception of photosynthetic Bradyrhizobium, despite high sequence diversity. The distribution of the nodZ, nolL, and noeI genes in the studied strains, as assessed by gene amplification, hybridization or sequencing, was found to correlate with the nodA tree topology. Moreover, the nodA, nodZ, and noeI phylogenies were largely congruent, but did not closely follow the taxonomy of the strains shown by the housekeeping 16S rRNA and dnaK genes. Additionally, the distribution of nodZ, noeI, and nolL genes suggested that their presence may be related to the requirements of their legume hosts. These data indicated that the spread and maintenance of nodulation genes within the Bradyrhizobium genus occurred through vertical transmission, although lateral gene transfer also played a significant role.

Diversification of Lupine Bradyrhizobium Strains: Evidence from Nodulation Gene Trees

ABSTRACT Bradyrhizobium strains isolated in Europe from Genisteae and serradella legumes form a distinct lineage, designated clade II, on nodulation gene trees. Clade II bradyrhizobia appear to prevail also in the soils of Western Australia and South Africa following probably accidental introduction with seeds of their lupine and serradella hosts. Given this potential for dispersal, we investigated Bradyrhizobium isolates originating from a range of native New World lupines, based on phylogenetic analyses of nodulation (nodA, nodZ, noeI) and housekeeping (atpD, dnaK, glnII, recA) genes. The housekeeping gene trees revealed considerable diversity among lupine bradyrhizobia, with most isolates placed in the Bradyrhizobium japonicum lineage, while some European strains were closely related to Bradyrhizobium canariense. The nodA gene tree resolved seven strongly supported groups (clades I to VII) that correlated with strain geographical origins and to some extent with major Lupinus clades. All European strains were placed in clade II, whereas only a minority of New World strains was placed in this clade. This work, as well as our previous studies, suggests that clade II diversified predominately in the Old World, possibly in the Mediterranean. Most New World isolates formed subclade III.2, nested in a large “pantropical” clade III, which appears to be New World in origin, although it also includes strains originating from nonlupine legumes. Trees generated using nodZ and noeI gene sequences accorded well with the nodA tree, but evidence is presented that the noeI gene may not be required for nodulation of lupine and that loss of this gene is occurring.

The Variable Part of the dnaK Gene as an Alternative Marker for Phylogenetic Studies of Rhizobia and Related Alpha Proteobacteria

DnaK is the 70 kDa chaperone that prevents protein aggregation and supports the refolding of damaged proteins. Due to sequence conservation and its ubiquity this chaperone has been widely used in phylogenetic studies. In this study, we applied the less conserved part that encodes the so-called alpha-subdomain of the substrate-binding domain of DnaK for phylogenetic analysis of rhizobia and related non-symbiotic alpha-Proteobacteria. A single 330 bp DNA fragment was routinely amplified from DNA templates isolated from the species of the genera, Azorhizobium, Bradyrhizobium, Mesorhizobium, Rhizobium and Sinorhizobium, but also from some non-symbiotic alpha Proteobacteria such as Blastochloris, Chelatobacter and Chelatococcus. Phylogenetic analyses revealed high congruence between dnaK sequences and 16S rDNA trees, but they were not identical. In contrast, the partition homogeneity tests revealed that dnaK sequence data could be combined with other housekeeping genes such as recA, atpD or glnA. The dnaK trees exhibited good resolution in the cases of the genera Mesorhizobium, Sinorhizobium and Rhizobium, even better than usually shown by 16S rDNA phylogeny. The dnaK phylogeny supported the close phylogenetic relationship of Rhizobium galegae and Agrobacterium tumefaciens (R. radiobacter) C58, which together formed a separate branch within the fast-growing rhizobia, albeit closer to the genus Sinorhizobium. The Rhizobium and Sinorhizobium genera carried an insertion composed of two amino acids, which additionally supported the phylogenetic affinity of these two genera, as well as their distinctness from the Mesorhizobium genus. Consistently with the phylogeny shown by 16S-23S rDNA intergenic region sequences, the dnaK trees divided the genus Bradyrhizobium into three main lineages, corresponding to B. japonicum, B. elkanii, and photosynthetic Bradyrhizobium strains that infect Aeschynomene plants. Our results suggest that the 330 bp dnaK sequences could be used as an additional taxonomic marker for rhizobia and related species (alternatively to the 16S rRNA gene phylogeny).

Cowpea and peanut in southern Africa are nodulated by diverse Bradyrhizobium strains harboring nodulation genes that belong to the large pantropical clade common in Africa.

Cowpea (Vigna unguiculata) and peanut (Arachis hypogaea) in southern Africa are nodulated by a genetically diverse group of Bradyrhizobium strains. To determine the identity of these bacteria, a collection of 22 isolates originating from the root nodules of both hosts in Botswana and South Africa was investigated using the combined sequences for the core genome genes rrs, recA, and glnII. These data separated the majority of the isolates into one of three unique lineages that most likely represent novel Bradyrhizobium species. Some isolates were also conspecific with B. yuanmingense and with B. elkanii, although none grouped with B. japonicum, B. canariense or B. liaoningense. To study the evolution of nodulation genes in these bacteria, the common nodulation gene, nodA, and host-specific nodulation genes, nodZ, noeE, and noeI, were analyzed. The nodA phylogeny showed that the cowpea and peanut Bradyrhizobium isolates represent various locally adapted groups or ecotypes that form part of Clade III of the seven known BradyrhizobiumnodA clades. This large and highly diverse clade comprises all strains from sub-Saharan Africa, as well as some originating from the Americas, Australia, Indonesia, China and Japan. Some similar groupings were supported by the other nodulation genes, although the overall phylogenies for the nodulation genes were incongruent with that inferred from the core genome genes, suggesting that horizontal gene transfer significantly influences the evolution of cowpea and peanut root-nodule bacteria. Furthermore, identification of the nodZ, noeI, and noeE genes in the isolates tested indicates that African Bradyrhizobium species may produce highly decorated nodulation factors, which potentially represent an important adaptation enabling nodulation of a great variety of legumes inhabiting the African continent.

Phylogeny of nodulation genes and symbiotic properties of Genista tinctoria bradyrhizobia

Pairwise comparisons of Genista tinctoria (dyer’s weed) rhizobium nodA, nodC, and nodZ gene sequences to those available in databanks revealed their highest sequence identities to nodulation loci of Bradyrhizobium sp. (Lupinus) strains and rhizobia from other genistoid legumes. On phylogenetic trees, genistoid microsymbionts were grouped together in monophyletic clusters, which suggested that their nodulation genes evolved from a common ancestor. G. tinctoria nodulators formed symbioses not only with the native host, but also with other plants of Genisteae tribe such as: Lupinus luteus, Sarothamnus scoparius, and Chamaecytisus ratisbonensis, and they were classified as the genistoid cross-inoculation group. The dyer’s weed root nodules were designated as indeterminate with apical meristem consisting of infected and uninfected cells.

Bradyrhizobium canariense and Bradyrhizobium japonicum are the two dominant rhizobium species in root nodules of lupin and serradella plants growing in Europe.

Forty three Bradyrhizobium strains isolated in Poland from root nodules of lupin species (Lupinus albus, L. angustifolius and L. luteus), and pink serradella (Ornithopus sativus) were examined based on phylogenetic analyses of three housekeeping (atpD, glnII and recA) and nodulation (nodA) gene sequences. Additionally, seven strains originating from root-nodules of yellow serradella (O. compressus) from Asinara Island (Italy) were included in this study. Phylogenetic trees revealed that 15 serradella strains, including all yellow serradella isolates, and six lupin strains grouped in Bradyrhizobium canariense (BC) clade, whereas eight strains from pink serradella and 15 lupin strains were assigned to Bradyrhizobium japonicum (BJ1). Apparently, these species are the two dominant groups in soils of central Europe, in the nodules of lupin and serradella plants. Only three strains belonged to other chromosomal lineages: one formed a cluster that was sister to B. canariense, one strain grouped outside the branch formed by B. japonicum super-group, and one strain occupied a distant position in the genus Bradyrhizobium, clustering with strains of the Rhodopseudomonas genus. All strains in nodulation nodA gene tree grouped in a cluster referred to as Clade II, which is in line with earlier data on this clade dominance among Bradyrhizobium strains in Europe. The nodA tree revealed four well-supported subgroups within Clade II (II.1-II.4). Interestingly, all B. canariense strains clustered in subgroup II.1 whereas B. japonicum strains dominated subgroups II.2-II.4.

Distinct Bradyrhizbium communities nodulate legumes native to temperate and tropical monsoon Australia

Geographic isolation and growing climate aridity played major roles in the evolution of Australian legumes. It is likely that these two factors also impacted on the evolution of their root-nodule bacteria. To investigate this issue, we applied a multilocus sequence analysis (MLSA) approach to examine Bradyrhizobium isolates originating from temperate areas of Western Australia (WA) and the tropical-monsoon area of the Northern Territory (NT). The isolates were mostly collected from the nodules of legumes belonging to tribes, genera and species endemic or native to Australia. Phylogenetic analyses of sequences for the housekeeping atpD, dnaK, glnII, gyrB, recA and 16S rRNA genes and nodulation nodA gene revealed that most isolates belonged to groups that are distinct from non-Australian Bradyrhizobium isolates, which is in line with earlier studies based on 16S rRNA gene sequence analyses. Phylogenetic analysis of the nodA data allowed identification of five major Clades among the WA and NT isolates. All WA isolates grouped in a subgroup I.1 of Clade I with strains originating from temperate eastern Australia. In contrast, the NT isolates formed part of Clades I (subgroup I.2), III (subgroup III.3), IV, V and X. Of these nodA clades, Clade I, Clade IV, Clade X presumably have an Australian origin. Overall, these data demonstrate that the impact of geographic isolation of the Australian landmass is manifested by the presence of numerous unique clusters in housekeeping and nodulation gene trees. In addition, the intrinsic climate characteristics of temperate WA and tropical-monsoon NT were responsible for the formation of distinct legume communities selecting for unrelated Bradyrhizobium groups.

Low sequence similarity and gene content of symbiotic clusters of Bradyrhizobium sp. WM9 (Lupinus) indicate early divergence of lupin lineage in the genus Bradyrhizobium

Two sequenced nodulation regions of lupin Bradyrhizobium sp. WM9 carried the majority of genes involved in the Nod factor production. The nod region I harbored: nolA, nodD, nodA, nodB, nodC, nodS, nodI, nodJ, nolO, nodZ, fixR, nifA, fixA, nodM, nolK and noeL. This gene arrangement resembled that found in the nodulation region of Bradyrhizobium japonicum USDA110, however strain WM9 harbored only one nodD gene copy, while the nodM, nolK and noeL genes had no counterparts in the 410 kb symbiotic region of strain USDA110. Region II harbored nolL and nodW, but lacked an nodV gene. Both regions carried ORFs that lacked similarity to the published USDA110 sequences, though they had homologues in symbiotic regions of Rhizobium etli, Sinorhizobium sp. NGR234 and Mesorhizobium loti. These differences in gene content, as well as a low average sequence identity (70%) of symbiotic genes with respect to B. japonicum USDA110 were in contrast with the phylogenetic relationship of USDA110 and WM9 revealed by the analysis of 16S rDNA and dnaK sequences. This most likely reflected an early divergence of symbiotic loci, and possible co-speciation with distinct legumes. During this process the loss of a noeI gene and the acquisition of a nolL gene could be regarded as an adaptation towards these legumes that responded to Nod factors carrying 4-O-acetylfucose rather than 2-O-methylfucose. This explained various responses of lupins and serradella plants to infection by mutants in nodZ and nolL genes, knowing that serradella is a stringent legume while lupins are more promiscuous legumes.

The root-nodule symbiosis between Sarothamnus scoparius L. and its microsymbionts

When nitrogen fixing root nodules are formed, Sarothamnus scoparius (broom) is inoculated with its microsymbionts. Nodules studied under light and electron microscopy exhibited typical indeterminate nodule histology with apical, persistent meristem, age gradient of nodule tissues, and open vascular bundles, and also with some particular features such as: the presence of mitotic activity in the infected meristematic cells, lack of infection threads, distribution of bacteria by process of host cell division, and occurrence of a large bacteroid zone only with infected cells. The results of cross-inoculation tests have shown a broad host range for S. scoparius microsymbionts including not only the native host but also species such as: Lupinus luteus, Ornithopus sativa, Lotus corniculatus, Genista tinctoria, Chamaecitisus ratisbonensis, Macroptilium atropurpureum, and Phaseolus vulgaris. In addition, our data established a close symbiotic relationship of S. scoparius nodulators to Bradyrhizobium sp. (Lupinus) by comparison of the partial sequence of nodC gene of the strain CYT7, specific for the broom, to those from Bradyrhizobium sp. (Lupinus) strain D1 and others available in the public databases.

Crystal Structures of NodS N-Methyltransferase from Bradyrhizobium japonicum in Ligand-Free Form and as SAH Complex.

NodS is an S-adenosyl-L-methionine (SAM)-dependent N-methyltransferase that is involved in the biosynthesis of Nod factor (NF) in rhizobia, which are bacterial symbionts of legume plants. NF is a modified chitooligosaccharide (COS) signal molecule that is recognized by the legume host, where it initiates symbiotic processes leading to atmospheric nitrogen fixation. We report the crystal structure of recombinant NodS protein from Bradyrhizobium japonicum, which infects lupine and serradella legumes. Two crystal forms--ligand-free NodS and NodS in complex with S-adenosyl-L-homocysteine, which is a by-product of the methylation reaction--were obtained, and their structures were refined to resolutions of 2.43 Å and 1.85 Å, respectively. Although the overall fold (consisting of a seven-stranded β-sheet flanked by layers of helices) is similar to those of other SAM-dependent methyltransferases, NodS has specific features reflecting the unique character of its oligosaccharide substrate. In particular, the N-terminal helix and its connecting loop get ordered upon SAM binding, thereby closing the methyl donor cavity and shaping a long surface canyon that is clearly the binding site for the acceptor molecule. Comparison of the two structural forms of NodS suggests that there are also other conformational changes taking place upon the binding of the donor substrate. As an enzyme that methylates a COS substrate, NodS is the first example among all SAM-dependent methyltransferases to have its three-dimensional structure elucidated. Gaining insight about how NodS binds its donor and acceptor substrates helps to better understand the mechanism of NodS activity and the basis of its functional difference in various rhizobia.