Hundreds of antimicrobial peptides create a selective barrier for insect gut symbionts.

The spatial organization of gut microbiota is crucial for the functioning of the gut ecosystem, although the mechanisms that organize gut bacterial communities in microhabitats are only partially understood. The gut of the insect Riptortus pedestris has a characteristic microbiota biogeography with a multispecies community in the anterior midgut and a monospecific bacterial population in the posterior midgut. We show that the posterior midgut region produces massively hundreds of specific antimicrobial peptides (AMPs), the Crypt-specific Cysteine-Rich peptides (CCRs) that have membrane-damaging antimicrobial activity against diverse bacteria but posterior midgut symbionts have elevated resistance. We determined by transposon-sequencing the genetic repertoire in the symbiont Caballeronia insecticola to manage CCR stress, identifying different independent pathways, including AMP-resistance pathways unrelated to known membrane homeostasis functions as well as cell envelope functions. Mutants in the corresponding genes have reduced capacity to colonize the posterior midgut, demonstrating that CCRs create a selective barrier and resistance is crucial in gut symbionts. Moreover, once established in the gut, the bacteria differentiate into a CCR-sensitive state, suggesting a second function of the CCR peptide arsenal in protecting the gut epithelia or mediating metabolic exchanges between the host and the gut symbionts. Our study highlights the evolution of an extreme diverse AMP family that likely contributes to establish and control the gut microbiota.

The noncoding RNA CcnA modulates the master cell cycle regulators CtrA and GcrA in Caulobacter crescentus.

Bacteria are powerful models for understanding how cells divide and accomplish global regulatory programs. In Caulobacter crescentus, a cascade of essential master regulators supervises the correct and sequential activation of DNA replication, cell division, and development of different cell types. Among them, the response regulator CtrA plays a crucial role coordinating all those functions. Here, for the first time, we describe the role of a novel factor named CcnA (cell cycle noncoding RNA A), a cell cycle-regulated noncoding RNA (ncRNA) located at the origin of replication, presumably activated by CtrA, and responsible for the accumulation of CtrA itself. In addition, CcnA may be also involved in the inhibition of translation of the S-phase regulator, GcrA, by interacting with its 5' untranslated region (5' UTR). Performing in vitro experiments and mutagenesis, we propose a mechanism of action of CcnA based on liberation (ctrA) or sequestration (gcrA) of their ribosome-binding site (RBS). Finally, its role may be conserved in other alphaproteobacterial species, such as Sinorhizobium meliloti, representing indeed a potentially conserved process modulating cell cycle in Caulobacterales and Rhizobiales.

Allosteric regulation of histidine kinases by their cognate response regulator determines cell fate.

The two-component phosphorylation network is of critical importance for bacterial growth and physiology. Here, we address plasticity and interconnection of distinct signal transduction pathways within this network. In Caulobacter crescentus antagonistic activities of the PleC phosphatase and DivJ kinase localized at opposite cell poles control the phosphorylation state and subcellular localization of the cell fate determinator protein DivK. We show that DivK functions as an allosteric regulator that switches PleC from a phosphatase into an autokinase state and thereby mediates a cyclic di-GMP-dependent morphogenetic program. Through allosteric activation of the DivJ autokinase, DivK also stimulates its own phosphorylation and polar localization. These data suggest that DivK is the central effector of an integrated circuit that operates via spatially organized feedback loops to control asymmetry and cell fate determination in C. crescentus. Thus, single domain response regulators can facilitate crosstalk, feedback control, and long-range communication among members of the two-component network.

OrpR is a σ54 -dependent activator using an iron-sulfur cluster for redox sensing in Desulfovibrio vulgaris Hildenborough.

Enhancer binding proteins (EBPs) are key players of σ54 -regulation that control transcription in response to environmental signals. In the anaerobic microorganism Desulfovibrio vulgaris Hildenborough (DvH), orp operons have been previously shown to be coregulated by σ54 -RNA polymerase, the integration host factor IHF and a cognate EBP, OrpR. In this study, ChIP-seq experiments indicated that the OrpR regulon consists of only the two divergent orp operons. In vivo data revealed that (i) OrpR is absolutely required for orp operons transcription, (ii) under anaerobic conditions, OrpR binds on the two dedicated DNA binding sites and leads to high expression levels of the orp operons, (iii) increasing the redox potential of the medium leads to a drastic down-regulation of the orp operons expression. Moreover, combining functional and biophysical studies on the anaerobically purified OrpR leads us to propose that OrpR senses redox potential variations via a redox-sensitive [4Fe-4S]2+ cluster in the sensory PAS domain. Overall, the study herein presents the first characterization of a new Fe-S redox regulator belonging to the σ54 -dependent transcriptional regulator family probably advantageously selected by cells adapted to the anaerobic lifestyle to monitor redox stress conditions.

Methylation-dependent transcriptional regulation of crescentin gene (creS) by GcrA in Caulobacter crescentus.

In Caulobacter crescentus the combined action of chromosome replication and the expression of DNA methyl-transferase CcrM at the end of S-phase maintains a cyclic alternation between a full- to hemi-methylated chromosome. This transition of the chromosomal methylation pattern affects the DNA-binding properties of the transcription factor GcrA that controls the several key cell cycle functions. However, the molecular mechanism by which GcrA and methylation are linked to transcription is not fully elucidated yet. Using a combination of cell biology, genetics, and in vitro analysis, we deciphered how GcrA integrates the methylation pattern of several S-phase expressed genes to their transcriptional output. We demonstrated in vitro that transcription of ctrA from the P1 promoter in its hemi-methylated state is activated by GcrA, while in its fully methylated state GcrA had no effect. Further, GcrA and methylation together influence a peculiar distribution of creS transcripts, encoding for crescentin, the protein responsible for the characteristic shape of Caulobacter cells. This gene is duplicated at the onset of chromosome replication and the two hemi-methylated copies are spatially segregated. Our results indicated that GcrA transcribed only the copy where coding strand is methylated. In vitro transcription assay further substantiated this finding. As several of the cell cycle-regulated genes are also under the influence of methylation and GcrA-dependent transcriptional regulation, this could be a mechanism responsible for maintaining the gene transcription dosage during the S-phase.

A new factor controlling cell envelope integrity in Alphaproteobacteria in the context of cell cycle, stress response and infection.

The bacterial envelope is a remarkable and complex compartment of the prokaryotic cell in which many essential functions take place. The article by Herrou and collaborators (Herrou et al., in press), by a clever combination of structural analysis, genetics and functional characterization in free-living bacterial cells and during infection in animal models, elucidates a new factor, named EipA, that plays a major role in Brucella spp envelope biogenesis and cell division. The authors demonstrate a genetic connection between eipA and lipopolysaccharide synthesis, specifically genes involved in the synthesis of the O-antigen portion of lipopolysaccharide (LPS). Beyond its crucial role in Brucella physiology, the conservation of EipA in the class Alphaproteobacteria urges microbiologists to pursue future investigation of its homologs in other species belonging to this important group of bacteria.

The nucleoid as a scaffold for the assembly of bacterial signaling complexes.

The FrzCD chemoreceptor from the gliding bacterium Myxococcus xanthus forms cytoplasmic clusters that occupy a large central region of the cell body also occupied by the nucleoid. In this work, we show that FrzCD directly binds to the nucleoid with its N-terminal positively charged tail and recruits active signaling complexes at this location. The FrzCD binding to the nucleoid occur in a DNA-sequence independent manner and leads to the formation of multiple distributed clusters that explore constrained areas. This organization might be required for cooperative interactions between clustered receptors as observed in membrane-bound chemosensory arrays.

CtrA controls cell division and outer membrane composition of the pathogen Brucella abortus.

Brucella abortus is a pathogen infecting cattle, able to survive, traffic, and proliferate inside host cells. It belongs to the Alphaproteobacteria, a phylogenetic group comprising bacteria with free living, symbiotic, and pathogenic lifestyles. An essential regulator of cell cycle progression named CtrA was described in the model bacterium Caulobacter crescentus. This regulator is conserved in many alphaproteobacteria, but the evolution of its regulon remains elusive. Here we identified promoters that are CtrA targets using ChIP-seq and we found that CtrA binds to promoters of genes involved in cell cycle progression, in addition to numerous genes encoding outer membrane components involved in export of membrane proteins and synthesis of lipopolysaccharide. Analysis of a conditional B. abortus ctrA loss of function mutant confirmed that CtrA controls cell division. Impairment of cell division generates elongated and branched morphologies, that are also detectable inside HeLa cells. Surprisingly, abnormal bacteria are able to traffic to the endoplasmic reticulum, the usual replication niche of B. abortus in host cells. We also found that CtrA depletion affected outer membrane composition, in particular the abundance and spatial distribution of Omp25. Control of the B. abortus envelope composition by CtrA indicates the plasticity of the CtrA regulon along evolution.

Cell Cycle Control by the Master Regulator CtrA in Sinorhizobium meliloti.

In all domains of life, proper regulation of the cell cycle is critical to coordinate genome replication, segregation and cell division. In some groups of bacteria, e.g. Alphaproteobacteria, tight regulation of the cell cycle is also necessary for the morphological and functional differentiation of cells. Sinorhizobium meliloti is an alphaproteobacterium that forms an economically and ecologically important nitrogen-fixing symbiosis with specific legume hosts. During this symbiosis S. meliloti undergoes an elaborate cellular differentiation within host root cells. The differentiation of S. meliloti results in massive amplification of the genome, cell branching and/or elongation, and loss of reproductive capacity. In Caulobacter crescentus, cellular differentiation is tightly linked to the cell cycle via the activity of the master regulator CtrA, and recent research in S. meliloti suggests that CtrA might also be key to cellular differentiation during symbiosis. However, the regulatory circuit driving cell cycle progression in S. meliloti is not well characterized in both the free-living and symbiotic state. Here, we investigated the regulation and function of CtrA in S. meliloti. We demonstrated that depletion of CtrA cause cell elongation, branching and genome amplification, similar to that observed in nitrogen-fixing bacteroids. We also showed that the cell cycle regulated proteolytic degradation of CtrA is essential in S. meliloti, suggesting a possible mechanism of CtrA depletion in differentiated bacteroids. Using a combination of ChIP-Seq and gene expression microarray analysis we found that although S. meliloti CtrA regulates similar processes as C. crescentus CtrA, it does so through different target genes. For example, our data suggest that CtrA does not control the expression of the Fts complex to control the timing of cell division during the cell cycle, but instead it negatively regulates the septum-inhibiting Min system. Our findings provide valuable insight into how highly conserved genetic networks can evolve, possibly to fit the diverse lifestyles of different bacteria.

Evolution of Intra-specific Regulatory Networks in a Multipartite Bacterial Genome.

Reconstruction of the regulatory network is an important step in understanding how organisms control the expression of gene products and therefore phenotypes. Recent studies have pointed out the importance of regulatory network plasticity in bacterial adaptation and evolution. The evolution of such networks within and outside the species boundary is however still obscure. Sinorhizobium meliloti is an ideal species for such study, having three large replicons, many genomes available and a significant knowledge of its transcription factors (TF). Each replicon has a specific functional and evolutionary mark; which might also emerge from the analysis of their regulatory signatures. Here we have studied the plasticity of the regulatory network within and outside the S. meliloti species, looking for the presence of 41 TFs binding motifs in 51 strains and 5 related rhizobial species. We have detected a preference of several TFs for one of the three replicons, and the function of regulated genes was found to be in accordance with the overall replicon functional signature: house-keeping functions for the chromosome, metabolism for the chromid, symbiosis for the megaplasmid. This therefore suggests a replicon-specific wiring of the regulatory network in the S. meliloti species. At the same time a significant part of the predicted regulatory network is shared between the chromosome and the chromid, thus adding an additional layer by which the chromid integrates itself in the core genome. Furthermore, the regulatory network distance was found to be correlated with both promoter regions and accessory genome evolution inside the species, indicating that both pangenome compartments are involved in the regulatory network evolution. We also observed that genes which are not included in the species regulatory network are more likely to belong to the accessory genome, indicating that regulatory interactions should also be considered to predict gene conservation in bacterial pangenomes.

DNA binding of the cell cycle transcriptional regulator GcrA depends on N6-adenosine methylation in Caulobacter crescentus and other Alphaproteobacteria.

Several regulators are involved in the control of cell cycle progression in the bacterial model system Caulobacter crescentus, which divides asymmetrically into a vegetative G1-phase (swarmer) cell and a replicative S-phase (stalked) cell. Here we report a novel functional interaction between the enigmatic cell cycle regulator GcrA and the N6-adenosine methyltransferase CcrM, both highly conserved proteins among Alphaproteobacteria, that are activated early and at the end of S-phase, respectively. As no direct biochemical and regulatory relationship between GcrA and CcrM were known, we used a combination of ChIP (chromatin-immunoprecipitation), biochemical and biophysical experimentation, and genetics to show that GcrA is a dimeric DNA-binding protein that preferentially targets promoters harbouring CcrM methylation sites. After tracing CcrM-dependent N6-methyl-adenosine promoter marks at a genome-wide scale, we show that these marks recruit GcrA in vitro and in vivo. Moreover, we found that, in the presence of a methylated target, GcrA recruits the RNA polymerase to the promoter, consistent with its role in transcriptional activation. Since methylation-dependent DNA binding is also observed with GcrA orthologs from other Alphaproteobacteria, we conclude that GcrA is the founding member of a new and conserved class of transcriptional regulators that function as molecular effectors of a methylation-dependent (non-heritable) epigenetic switch that regulates gene expression during the cell cycle.

A Salmonella type three secretion effector/chaperone complex adopts a hexameric ring-like structure.

Many bacterial pathogens use type three secretion systems (T3SS) to inject virulence factors, named effectors, directly into the cytoplasm of target eukaryotic cells. Most of the T3SS components are conserved among plant and animal pathogens, suggesting a common mechanism of recognition and secretion of effectors. However, no common motif has yet been identified for effectors allowing T3SS recognition. In this work, we performed a biochemical and structural characterization of the Salmonella SopB/SigE chaperone/effector complex by small-angle X-ray scattering (SAXS). Our results showed that the SopB/SigE complex is assembled in dynamic homohexameric-ring-shaped structures with an internal tunnel. In this ring, the chaperone maintains a disordered N-terminal end of SopB molecules, in a good position to be reached and processed by the T3SS. This ring dimensionally fits the ring-organized molecules of the injectisome, including ATPase hexameric rings; this organization suggests that this structural feature is important for ATPase recognition by T3SS. Our work constitutes the first evidence of the oligomerization of an effector, analogous to the organization of the secretion machinery, obtained in solution. As effectors share neither sequence nor structural identity, the quaternary oligomeric structure could constitute a strategy evolved to promote the specificity and efficiency of T3SS recognition.

DuctApe: a suite for the analysis and correlation of genomic and OmniLog™ Phenotype Microarray data.

Addressing the functionality of genomes is one of the most important and challenging tasks of today's biology. In particular the ability to link genotypes to corresponding phenotypes is of interest in the reconstruction and biotechnological manipulation of metabolic pathways. Over the last years, the OmniLog™ Phenotype Microarray (PM) technology has been used to address many specific issues related to the metabolic functionality of microorganisms. However, computational tools that could directly link PM data with the gene(s) of interest followed by the extraction of information on gene-phenotype correlation are still missing. Here we present DuctApe, a suite that allows the analysis of both genomic sequences and PM data, to find metabolic differences among PM experiments and to correlate them with KEGG pathways and gene presence/absence patterns. As example, an application of the program to four bacterial datasets is presented. The source code and tutorials are available at http://combogenomics.github.io/DuctApe/.

The DivJ, CbrA and PleC system controls DivK phosphorylation and symbiosis in Sinorhizobium meliloti.

Sinorhizobium meliloti is a soil bacterium that invades the root nodules it induces on Medicago sativa, whereupon it undergoes an alteration of its cell cycle and differentiates into nitrogen-fixing, elongated and polyploid bacteroid with higher membrane permeability. In Caulobacter crescentus, a related alphaproteobacterium, the principal cell cycle regulator, CtrA, is inhibited by the phosphorylated response regulator DivK. The phosphorylation of DivK depends on the histidine kinase DivJ, while PleC is the principal phosphatase for DivK. Despite the importance of the DivJ in C. crescentus, the mechanistic role of this kinase has never been elucidated in other Alphaproteobacteria. We show here that the histidine kinases DivJ together with CbrA and PleC participate in a complex phosphorylation system of the essential response regulator DivK in S. meliloti. In particular, DivJ and CbrA are involved in DivK phosphorylation and in turn CtrA inactivation, thereby controlling correct cell cycle progression and the integrity of the cell envelope. In contrast, the essential PleC presumably acts as a phosphatase of DivK. Interestingly, we found that a DivJ mutant is able to elicit nodules and enter plant cells, but fails to establish an effective symbiosis suggesting that proper envelope and/or low CtrA levels are required for symbiosis.

Analysis of the CtrA pathway in Magnetospirillum reveals an ancestral role in motility in alphaproteobacteria.

Developmental events across the prokaryotic life cycle are highly regulated at the transcriptional and posttranslational levels. Key elements of a few regulatory networks are conserved among phylogenetic groups of bacteria, although the features controlled by these conserved systems are as diverse as the organisms encoding them. In this work, we probed the role of the CtrA regulatory network, conserved throughout the Alphaproteobacteria, in the magnetotactic bacterium Magnetospirillum magneticum strain AMB-1, which possesses unique intracellular organization and compartmentalization. While we have shown that CtrA in AMB-1 is not essential for viability, it is required for motility, and its putative phosphorylation state dictates the ability of CtrA to activate the flagellar biosynthesis gene cascade. Gene expression analysis of strains expressing active and inactive CtrA alleles points to the composition of the extended CtrA regulon, including both direct and indirect targets. These results, combined with a bioinformatic study of the AMB-1 genome, enabled the prediction of an AMB-1-specific CtrA binding site. Further, phylogenetic studies comparing CtrA sequences from alphaproteobacteria in which the role of CtrA has been experimentally examined reveal an ancestral role of CtrA in the regulation of motility and suggest that its essential functions in other alphaproteobacteria were acquired subsequently.

Exploring the plant-associated bacterial communities in Medicago sativa L.

BACKGROUND: Plant-associated bacterial communities caught the attention of several investigators which study the relationships between plants and soil and the potential application of selected bacterial species in crop improvement and protection. Medicago sativa L. is a legume crop of high economic importance as forage in temperate areas and one of the most popular model plants for investigations on the symbiosis with nitrogen fixing rhizobia (mainly belonging to the alphaproteobacterial species Sinorhizobium meliloti). However, despite its importance, no studies have been carried out looking at the total bacterial community associated with the plant. In this work we explored for the first time the total bacterial community associated with M. sativa plants grown in mesocosms conditions, looking at a wide taxonomic spectrum, from the class to the single species (S. meliloti) level. RESULTS: Results, obtained by using Terminal-Restriction Fragment Length Polymorphism (T-RFLP) analysis, quantitative PCR and sequencing of 16 S rRNA gene libraries, showed a high taxonomic diversity as well as a dominance by members of the class Alphaproteobacteria in plant tissues. Within Alphaproteobacteria the families Sphingomonadaceae and Methylobacteriaceae were abundant inside plant tissues, while soil Alphaproteobacteria were represented by the families of Hyphomicrobiaceae, Methylocystaceae, Bradyirhizobiaceae and Caulobacteraceae. At the single species level, we were able to detect the presence of S. meliloti populations in aerial tissues, nodules and soil. An analysis of population diversity on nodules and soil showed a relatively low sharing of haplotypes (30-40%) between the two environments and between replicate mesocosms, suggesting drift as main force shaping S. meliloti population at least in this system. CONCLUSIONS: In this work we shed some light on the bacterial communities associated with M. sativa plants, showing that Alphaproteobacteria may constitute an important part of biodiversity in this system, which includes also the well known symbiont S. meliloti. Interestingly, this last species was also found in plant aerial part, by applying cultivation-independent protocols, and a genetic diversity analysis suggested that population structure could be strongly influenced by random drift.

Structural insights into ChpT, an essential dimeric histidine phosphotransferase regulating the cell cycle in Caulobacter crescentus.

Two-component and phosphorelay signal-transduction proteins are crucial for bacterial cell-cycle regulation in Caulobacter crescentus. ChpT is an essential histidine phosphotransferase that controls the activity of the master cell-cycle regulator CtrA by phosphorylation. Here, the 2.2 Å resolution crystal structure of ChpT is reported. ChpT is a homodimer and adopts the domain architecture of the intracellular part of class I histidine kinases. Each subunit consists of two distinct domains: an N-terminal helical hairpin domain and a C-terminal α/ß domain. The two N-terminal domains are adjacent within the dimer, forming a four-helix bundle. The ChpT C-terminal domain adopts an atypical Bergerat ATP-binding fold.

Regulation of the bacterial cell cycle by an integrated genetic circuit.

How bacteria regulate cell cycle progression at a molecular level is a fundamental but poorly understood problem. In Caulobacter crescentus, two-component signal transduction proteins are crucial for cell cycle regulation, but the connectivity of regulators involved has remained elusive and key factors are unidentified. Here we identify ChpT, an essential histidine phosphotransferase that controls the activity of CtrA, the master cell cycle regulator. We show that the essential histidine kinase CckA initiates two phosphorelays, each requiring ChpT, which lead to the phosphorylation and stabilization of CtrA. Downregulation of CckA activity therefore results in the dephosphorylation and degradation of CtrA, which in turn allow the initiation of DNA replication. Furthermore, we show that CtrA triggers its own destruction by promoting cell division and inducing synthesis of the essential regulator DivK, which feeds back to downregulate CckA immediately before S phase. Our results define a single integrated circuit whose components and connectivity can account for the cell cycle oscillations of CtrA in Caulobacter.

DNA Methylation in Ensifer Species during Free-Living Growth and during Nitrogen-Fixing Symbiosis with Medicago spp.

Methylation of specific DNA sequences is ubiquitous in bacteria and has known roles in immunity and regulation of cellular processes, such as the cell cycle. Here, we explored DNA methylation in bacteria of the genus Ensifer, including its potential role in regulating terminal differentiation during nitrogen-fixing symbiosis with legumes. Using single-molecule real-time sequencing, six genome-wide methylated motifs were identified across four Ensifer strains, five of which were strain-specific. Only the GANTC motif, recognized by the cell cycle-regulated CcrM methyltransferase, was methylated in all strains. In actively dividing cell cultures, methylation of GANTC motifs increased progressively from the ori to ter regions in each replicon, in agreement with a cell cycle-dependent regulation of CcrM. In contrast, there was near full genome-wide GANTC methylation in the early stage of symbiotic differentiation. This was followed by a moderate decrease in the overall extent of methylation and a progressive decrease in chromosomal GANTC methylation from the ori to ter regions in later stages of differentiation. Based on these observations, we suggest that CcrM activity is dysregulated and constitutive during terminal differentiation, which we hypothesize is a driving factor for endoreduplication of terminally differentiated bacteroids. IMPORTANCE Nitrogen fixation by rhizobia in symbiosis with legumes is economically and ecologically important. The symbiosis can involve a complex bacterial transformation-terminal differentiation-that includes major shifts in the transcriptome and cell cycle. Epigenetic regulation is an important regulatory mechanism in diverse bacteria; however, the roles of DNA methylation in rhizobia and symbiotic nitrogen fixation have been poorly investigated. We show that aside from cell cycle regulation, DNA methyltransferases are unlikely to have conserved roles in the biology of bacteria of the genus Ensifer. However, we present evidence consistent with an interpretation that the cell cycle methyltransferase CcrM is dysregulated during symbiosis, which we hypothesize may be a key factor driving the cell cycle switch in terminal differentiation required for effective symbioses.

Exploring the symbiotic pangenome of the nitrogen-fixing bacterium Sinorhizobium meliloti.

BACKGROUND: Sinorhizobium meliloti is a model system for the studies of symbiotic nitrogen fixation. An extensive polymorphism at the genetic and phenotypic level is present in natural populations of this species, especially in relation with symbiotic promotion of plant growth. AK83 and BL225C are two nodule-isolated strains with diverse symbiotic phenotypes; BL225C is more efficient in promoting growth of the Medicago sativa plants than strain AK83. In order to investigate the genetic determinants of the phenotypic diversification of S. meliloti strains AK83 and BL225C, we sequenced the complete genomes for these two strains. RESULTS: With sizes of 7.14 Mbp and 6.97 Mbp, respectively, the genomes of AK83 and BL225C are larger than the laboratory strain Rm1021. The core genome of Rm1021, AK83, BL225C strains included 5124 orthologous groups, while the accessory genome was composed by 2700 orthologous groups. While Rm1021 and BL225C have only three replicons (Chromosome, pSymA and pSymB), AK83 has also two plasmids, 260 and 70 Kbp long. We found 65 interesting orthologous groups of genes that were present only in the accessory genome, consequently responsible for phenotypic diversity and putatively involved in plant-bacterium interaction. Notably, the symbiosis inefficient AK83 lacked several genes required for microaerophilic growth inside nodules, while several genes for accessory functions related to competition, plant invasion and bacteroid tropism were identified only in AK83 and BL225C strains. Presence and extent of polymorphism in regulons of transcription factors involved in symbiotic interaction were also analyzed. Our results indicate that regulons are flexible, with a large number of accessory genes, suggesting that regulons polymorphism could also be a key determinant in the variability of symbiotic performances among the analyzed strains. CONCLUSIONS: In conclusions, the extended comparative genomics approach revealed a variable subset of genes and regulons that may contribute to the symbiotic diversity.