Effect of Non-Lethal Selection on Spontaneous Revertants of Frameshift Mutations: The Escherichia coli hisF Case

Microorganisms possess the potential to adapt to fluctuations in environmental parameters, and their evolution is driven by the continuous generation of mutations. The reversion of auxotrophic mutations has been widely studied; however, little is known about the reversion of frameshift mutations resulting in amino acid auxotrophy and on the structure and functioning of the protein encoded by the revertant mutated gene. The aims of this work were to analyze the appearance of reverse mutations over time and under different selective pressures and to investigate revertant enzymes' three-dimensional structures and their correlation with a different growth ability. Escherichia coli FB182 strain, carrying the hisF892 single nucleotide deletion resulting in histidine auxotrophy, was subjected to different selective pressures, and revertant mutants were isolated and characterized. The obtained results allowed us to identify different indels of different lengths located in different positions in the hisF gene, and relations with the incubation time and the selective pressure applied were observed. Moreover, the structure of the different mutant proteins was consistent with the respective revertant ability to grow in absence of histidine, highlighting a correlation between the mutations and the catalytic activity of the mutated HisF enzyme.

Comparison of Highly and Weakly Virulent Dickeya solani Strains, With a View on the Pangenome and Panregulon of This Species.

Bacteria belonging to the genera Dickeya and Pectobacterium are responsible for significant economic losses in a wide variety of crops and ornamentals. During last years, increasing losses in potato production have been attributed to the appearance of Dickeya solani. The D. solani strains investigated so far share genetic homogeneity, although different virulence levels were observed among strains of various origins. The purpose of this study was to investigate the genetic traits possibly related to the diverse virulence levels by means of comparative genomics. First, we developed a new genome assembly pipeline which allowed us to complete the D. solani genomes. Four de novo sequenced and ten publicly available genomes were used to identify the structure of the D. solani pangenome, in which 74.8 and 25.2% of genes were grouped into the core and dispensable genome, respectively. For D. solani panregulon analysis, we performed a binding site prediction for four transcription factors, namely CRP, KdgR, PecS and Fur, to detect the regulons of these virulence regulators. Most of the D. solani potential virulence factors were predicted to belong to the accessory regulons of CRP, KdgR, and PecS. Thus, some differences in gene expression could exist between D. solani strains. The comparison between a highly and a low virulent strain, IFB0099 and IFB0223, respectively, disclosed only small differences between their genomes but significant differences in the production of virulence factors like pectinases, cellulases and proteases, and in their mobility. The D. solani strains also diverge in the number and size of prophages present in their genomes. Another relevant difference is the disruption of the adhesin gene fhaB2 in the highly virulent strain. Strain IFB0223, which has a complete adhesin gene, is less mobile and less aggressive than IFB0099. This suggests that in this case, mobility rather than adherence is needed in order to trigger disease symptoms. This study highlights the utility of comparative genomics in predicting D. solani traits involved in the aggressiveness of this emerging plant pathogen.

Creation and Characterization of a Genomically Hybrid Strain in the Nitrogen-Fixing Symbiotic Bacterium Sinorhizobium meliloti.

Many bacteria, often associated with eukaryotic hosts and of relevance for biotechnological applications, harbor a multipartite genome composed of more than one replicon. Biotechnologically relevant phenotypes are often encoded by genes residing on the secondary replicons. A synthetic biology approach to developing enhanced strains for biotechnological purposes could therefore involve merging pieces or entire replicons from multiple strains into a single genome. Here we report the creation of a genomic hybrid strain in a model multipartite genome species, the plant-symbiotic bacterium Sinorhizobium meliloti. We term this strain as cis-hybrid, since it is produced by genomic material coming from the same species' pangenome. In particular, we moved the secondary replicon pSymA (accounting for nearly 20% of total genome content) from a donor S. meliloti strain to an acceptor strain. The cis-hybrid strain was screened for a panel of complex phenotypes (carbon/nitrogen utilization phenotypes, intra- and extracellular metabolomes, symbiosis, and various microbiological tests). Additionally, metabolic network reconstruction and constraint-based modeling were employed for in silico prediction of metabolic flux reorganization. Phenotypes of the cis-hybrid strain were in good agreement with those of both parental strains. Interestingly, the symbiotic phenotype showed a marked cultivar-specific improvement with the cis-hybrid strains compared to both parental strains. These results provide a proof-of-principle for the feasibility of genome-wide replicon-based remodelling of bacterial strains for improved biotechnological applications in precision agriculture.

Trade, Diplomacy, and Warfare: The Quest for Elite Rhizobia Inoculant Strains.

Rhizobia form symbiotic nitrogen-fixing nodules on leguminous plants, which provides an important source of fixed nitrogen input into the soil ecosystem. The improvement of symbiotic nitrogen fixation is one of the main challenges facing agriculture research. Doing so will reduce the usage of chemical nitrogen fertilizer, contributing to the development of sustainable agriculture practices to deal with the increasing global human population. Sociomicrobiological studies of rhizobia have become a model for the study of the evolution of mutualistic interactions. The exploitation of the wide range of social interactions rhizobia establish among themselves, with the soil and root microbiota, and with the host plant, could constitute a great advantage in the development of a new generation of highly effective rhizobia inoculants. Here, we provide a brief overview of the current knowledge on three main aspects of rhizobia interaction: trade of fixed nitrogen with the plant; diplomacy in terms of communication and possible synergistic effects; and warfare, as antagonism and plant control over symbiosis. Then, we propose new areas of investigation and the selection of strains based on the combination of the genetic determinants for the relevant rhizobia symbiotic behavioral phenotypes.

Role and Regulation of ACC Deaminase Gene in Sinorhizobium meliloti: Is It a Symbiotic, Rhizospheric or Endophytic Gene?

Plant-associated bacteria exhibit a number of different strategies and specific genes allow bacteria to communicate and metabolically interact with plant tissues. Among the genes found in the genomes of plant-associated bacteria, the gene encoding the enzyme 1-aminocyclopropane-1-carboxylate (ACC) deaminase (acdS) is one of the most diffused. This gene is supposed to be involved in the cleaving of plant-produced ACC, the precursor of the plant stress-hormone ethylene toning down the plant response to infection. However, few reports are present on the actual role in rhizobia, one of the most investigated groups of plant-associated bacteria. In particular, still unclear is the origin and the role of acdS in symbiotic competitiveness and on the selective benefit it may confer to plant symbiotic rhizobia. Here we present a phylogenetic and functional analysis of acdS orthologs in the rhizobium model-species Sinorhizobium meliloti. Results showed that acdS orthologs present in S. meliloti pangenome have polyphyletic origin and likely spread through horizontal gene transfer, mediated by mobile genetic elements. When acdS ortholog from AK83 strain was cloned and assayed in S. meliloti 1021 (lacking acdS), no modulation of plant ethylene levels was detected, as well as no increase in fitness for nodule occupancy was found in the acdS-derivative strain compared to the parental one. Surprisingly, AcdS was shown to confer the ability to utilize formamide and some dipeptides as sole nitrogen source. Finally, acdS was shown to be negatively regulated by a putative leucine-responsive regulator (LrpL) located upstream to acdS sequence (acdR). acdS expression was induced by root exudates of both legumes and non-leguminous plants. We conclude that acdS in S. meliloti is not directly related to symbiotic interaction, but it could likely be involved in the rhizospheric colonization or in the endophytic behavior.

Metabolic modelling reveals the specialization of secondary replicons for niche adaptation in Sinorhizobium meliloti.

The genome of about 10% of bacterial species is divided among two or more large chromosome-sized replicons. The contribution of each replicon to the microbial life cycle (for example, environmental adaptations and/or niche switching) remains unclear. Here we report a genome-scale metabolic model of the legume symbiont Sinorhizobium meliloti that is integrated with carbon utilization data for 1,500 genes with 192 carbon substrates. Growth of S. meliloti is modelled in three ecological niches (bulk soil, rhizosphere and nodule) with a focus on the role of each of its three replicons. We observe clear metabolic differences during growth in the tested ecological niches and an overall reprogramming following niche switching. In silico examination of the inferred fitness of gene deletion mutants suggests that secondary replicons evolved to fulfil a specialized function, particularly host-associated niche adaptation. Thus, genes on secondary replicons might potentially be manipulated to promote or suppress host interactions for biotechnological purposes.

Mixed Nodule Infection in Sinorhizobium meliloti-Medicago sativa Symbiosis Suggest the Presence of Cheating Behavior.

In the symbiosis between rhizobia and legumes, host plants can form symbiotic root nodules with multiple rhizobial strains, potentially showing different symbiotic performances in nitrogen fixation. Here, we investigated the presence of mixed nodules, containing rhizobia with different degrees of mutualisms, and evaluate their relative fitness in the Sinorhizobium meliloti-Medicago sativa model symbiosis. We used three S. meliloti strains, the mutualist strains Rm1021 and BL225C and the non-mutualist AK83. We performed competition experiments involving both in vitro and in vivo symbiotic assays with M. sativa host plants. We show the occurrence of a high number (from 27 to 100%) of mixed nodules with no negative effect on both nitrogen fixation and plant growth. The estimation of the relative fitness as non-mutualist/mutualist ratios in single nodules shows that in some nodules the non-mutualist strain efficiently colonized root nodules along with the mutualist ones. In conclusion, we can support the hypothesis that in S. meliloti-M. sativa symbiosis mixed nodules are formed and allow non-mutualist or less-mutualist bacterial partners to be less or not sanctioned by the host plant, hence allowing a potential form of cheating behavior to be present in the nitrogen fixing symbiosis.

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.

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.

Evaluation of the Performances of Ribosomal Database Project (RDP) Classifier for Taxonomic Assignment of 16S rRNA Metabarcoding Sequences Generated from Illumina-Solexa NGS.

Here we report a benchmark of the effect of bootstrap cut-off values of the RDP Classifier tool in terms of data retention along the different taxonomic ranks by using Illumina reads. Results provide guidelines for planning sequencing depths and selection of bootstrap cut-off in taxonomic assignments.

Exploring the dynamics of bacterial community composition in soil: the pan-bacteriome approach.

We performed a longitudinal study (repeated observations of the same sample over time) to investigate both the composition and structure of temporal changes of bacterial community composition in soil mesocosms, subjected to three different treatments (water and 5 or 25 mg kg(-1) of dried soil Cd(2+)). By analogy with the pan genome concept, we identified a core bacteriome and an accessory bacteriome. Resident taxa were assigned to the core bacteriome, while occasional taxa were assigned to the accessory bacteriome. Core and accessory bacteriome represented roughly 35 and 50 % of the taxa detected, respectively, and were characterized by different taxonomic signatures from phylum to genus level while 15 % of the taxa were found to be unique to a particular sample. In particular, the core bacteriome was characterized by higher abundance of members of Planctomycetes, Actinobacteria, Verrucomicrobia and Acidobacteria, while the accessory bacteriome included more members of Firmicutes, Clamydiae and Proteobacteria, suggesting potentially different responses to environmental changes of members from these phyla. We conclude that the pan-bacteriome model may be a useful approach to gain insight for modeling bacterial community structure and inferring different abilities of bacteria taxa.

Mapping contigs using CONTIGuator.

Obtaining bacterial genomic sequences has become a routine task in today's biology. The emergence of the comparative genomics approach has led to an increasing number of bacterial species having more than one strain sequenced, thus facilitating the annotation process. On the other hand, many genomic sequences are now left in the "draft" status, as a series of contigs, mainly for the labor-intensive finishing task. As a result, many genomic analyses are incomplete (e.g., in their annotation) or impossible to be performed (e.g., structural genomics analysis). Many approaches have been recently developed to facilitate the finishing process or at least to produce higher quality scaffolds; taking advantage of the comparative genomics paradigm, closely related genomes are used to align the contigs and determine their relative orientation, thus facilitating the finishing process, but also producing higher quality scaffolds. In this chapter we present the use of the CONTIGuator algorithm, which aligns the contigs from a draft genome to a closely related closed genome and resolves their relative orientation based on this alignment, producing a scaffold and a series of PCR primer pairs for the finishing process. The CONTIGuator algorithm is also capable of handling multipartite genomes (i.e., genomes having chromosomes and other plasmids), univocally mapping contigs to the most similar replicon. The program also produces a series of contig maps that allow to perform structural genomics analysis on the draft genome. The functionalities of the web interface, as well as the command line version, are presented.

Biogeography of Sinorhizobium meliloti nodulating alfalfa in different Croatian regions.

Sinorhizobium meliloti is a nitrogen-fixing rhizobium symbiont of legumes, widespread in many temperate environments the high genetic diversity of which enables it to thrive as a symbiont of host legumes and free-living in soil. Soil type, together with geographic differences and host plant genotype, seem to be prominent factors in shaping rhizobial genetic diversity. While a large body of research supports the idea that the genetic structure of free-living microbial taxa exhibits a clear biogeographic pattern, few investigations have been performed on the biogeographic pattern of S. meliloti genotypes in a restricted geographic range. In the present study, a collection of 128 S. meliloti isolates from three different regions in Croatia was investigated to analyze the relationship between genetic diversity, geographic distribution, soil features and isolate phenotypes by using amplified fragment length polymorphism (AFLP) as a genome-wide scanning method. Results obtained led to the conclusion that the genotypes of isolates cluster according to the region of origin and that the differentiation of S. meliloti populations can be mainly ascribed to geographic isolation following an isolation-by-distance model, with a strong distance-decay relationship of genetic similarity with distance, in which local soil conditions are not the major component influencing the isolate phenotypes or their genomic differentiation.

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.

Replicon-dependent bacterial genome evolution: the case of Sinorhizobium meliloti.

Many bacterial species, such as the alphaproteobacterium Sinorhizobium meliloti, are characterized by open pangenomes and contain multipartite genomes consisting of a chromosome and other large-sized replicons, such as chromids, megaplasmids, and plasmids. The evolutionary forces in both functional and structural aspects that shape the pangenome of species with multipartite genomes are still poorly understood. Therefore, we sequenced the genomes of 10 new S. meliloti strains, analyzed with four publicly available additional genomic sequences. Results indicated that the three main replicons present in these strains (a chromosome, a chromid, and a megaplasmid) partly show replicon-specific behaviors related to strain differentiation. In particular, the pSymB chromid was shown to be a hot spot for positively selected genes, and, unexpectedly, genes resident in the pSymB chromid were also found to be more widespread in distant taxa than those located in the other replicons. Moreover, through the exploitation of a DNA proximity network, a series of conserved "DNA backbones" were found to shape the evolution of the genome structure, with the rest of the genome experiencing rearrangements. The presented data allow depicting a scenario where the pSymB chromid has a distinctive role in intraspecies differentiation and in evolution through positive selection, whereas the pSymA megaplasmid mostly contributes to structural fluidity and to the emergence of new functions, indicating a specific evolutionary role for each replicon in the pangenome evolution.

Permanent draft genome sequences of the symbiotic nitrogen fixing Ensifer meliloti strains BO21CC and AK58.

Ensifer (syn. Sinorhizobium) meliloti is an important symbiotic bacterial species that fixes nitrogen. Strains BO21CC and AK58 were previously investigated for their substrate utilization and their plant-growth promoting abilities showing interesting features. Here, we describe the complete genome sequence and annotation of these strains. BO21CC and AK58 genomes are 6,985,065 and 6,974,333 bp long with 6,746 and 6,992 genes predicted, respectively.

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.

Spread of the group II intron RmInt1 and its insertion sequence target sites in the plant endosymbiont Sinorhizobium meliloti.

RmInt1 is a mobile group II intron from Sinorhizobium meliloti that is exceptionally abundant in this bacterial species. We compared the presence of RmInt1 and two of its insertion sequence homing sites (ISRm2011-2 and ISRm10-2) in two phylogenetic clusters (I and II) identified by AFLP analysis in a collection of S. meliloti field isolates from Italy. Both clusters contained several copies of the ISRm2011-2 element, which is present at high copy number in almost all S. meliloti isolates. By contrast, isolates from cluster I harbored no copies of ISRm10-2 and only a truncated copy of RmInt1, despite the absence of constraints on intron mobility in this genetic background, whereas cluster II strains harbored several copies of this intron. The absence of ISRm10-2 from one of the strains of this cluster suggests that this element was acquired more recently than the other two elements. Furthermore, studies of insertional polymorphisms in cluster II strains revealed the acquisition of ISRm10-2 and subsequent retrohoming of RmInt1 to this homing site. These results highlight the role of intron homing sites (ISs) in facilitating intron dispersal and the dynamic and ongoing nature of the spread of the group II intron RmInt1 in S. meliloti.

CONTIGuator: a bacterial genomes finishing tool for structural insights on draft genomes.

Recent developments in sequencing technologies have given the opportunity to sequence many bacterial genomes with limited cost and labor, compared to previous techniques. However, a limiting step of genome sequencing is the finishing process, needed to infer the relative position of each contig and close sequencing gaps. An additional degree of complexity is given by bacterial species harboring more than one replicon, which are not contemplated by the currently available programs. The availability of a large number of bacterial genomes allows geneticists to use complete genomes (possibly from the same species) as templates for contigs mapping.Here we present CONTIGuator, a software tool for contigs mapping over a reference genome which allows the visualization of a map of contigs, underlining loss and/or gain of genetic elements and permitting to finish multipartite genomes. The functionality of CONTIGuator was tested using four genomes, demonstrating its improved performances compared to currently available programs.Our approach appears efficient, with a clear visualization, allowing the user to perform comparative structural genomics analysis on draft genomes. CONTIGuator is a Python script for Linux environments and can be used on normal desktop machines and can be downloaded from http://contiguator.sourceforge.net.