Metabolic capacity of Sinorhizobium (Ensifer) meliloti strains as determined by phenotype MicroArray analysis

Metabolically versatile psychrotolerant Antarctic bacterium Pseudomonas sp. ANT_H12B is an efficient producer of siderophores and accompanying metabolites (SAM) useful for agricultural purposes

BACKGROUND

Bacterial siderophores are chelating compounds with the potential of application in agriculture, due to their plant growth-promoting (PGP) properties, however, high production and purification costs are limiting factors for their wider application. Cost-efficiency of the production could be increased by omitting purification processes, especially since siderophores accompanying metabolites (SAM) often also possess PGP traits. In this study, the metabolism versatility of Pseudomonas sp. ANT_H12B was used for the optimization of siderophores production and the potential of these metabolites and SAM was characterized in the context of PGP properties.

RESULTS

The metabolic diversity of ANT_H12B was examined through genomic analysis and phenotype microarrays. The strain was found to be able to use numerous C, N, P, and S sources, which allowed for the design of novel media suitable for efficient production of siderophores in the form of pyoverdine (223.50-512.60 μM). Moreover, depending on the culture medium, the pH of the siderophores and SAM solutions varied from acidic (pH < 5) to alkaline (pH > 8). In a germination test, siderophores and SAM were shown to have a positive effect on plants, with a significant increase in germination percentage observed in beetroot, pea, and tobacco. The PGP potential of SAM was further elucidated through GC/MS analysis, which revealed other compounds with PGP potential, such as indolic acetic acids, organic acids, fatty acids, sugars and alcohols. These compounds not only improved seed germination but could also potentially be beneficial for plant fitness and soil quality.

CONCLUSIONS

Pseudomonas sp. ANT_H12B was presented as an efficient producer of siderophores and SAM which exhibit PGP potential. It was also shown that omitting downstream processes could not only limit the costs of siderophores production but also improve their agricultural potential.

McpT, a Broad-Range Carboxylate Chemoreceptor in Sinorhizobium meliloti

Chemoreceptors enable the legume symbiont Sinorhizobium meliloti to detect and respond to specific chemicals released from their host plant alfalfa, which allows the establishment of a nitrogen-fixing symbiosis. The periplasmic region (PR) of transmembrane chemoreceptors act as the sensory input module for chemotaxis systems via binding of specific ligands, either directly or indirectly. S. meliloti has six transmembrane and two cytosolic chemoreceptors. However, the function of only three of the transmembrane receptors have been characterized so far, with McpU, McpV, and McpX serving as general amino acid, short-chain carboxylate, and quaternary ammonium compound sensors, respectively. In the present study, we analyzed the S. meliloti chemoreceptor McpT. High-throughput differential scanning fluorimetry assays, using Biolog phenotype microarray plates, identified 15 potential ligands for McpT^PR, with the majority classified as mono-, di-, and tricarboxylates. S. meliloti exhibited positive chemotaxis toward seven selected carboxylates, namely, α-ketobutyrate, citrate, glyoxylate, malate, malonate, oxalate, and succinate. These carboxylates were detected in seed exudates of the alfalfa host. Deletion of mcpT resulted in a significant decrease of chemotaxis to all carboxylates except for citrate. Isothermal titration calorimetry revealed that McpT^PR bound preferentially to the monocarboxylate glyoxylate and with lower affinity to the dicarboxylates malate, malonate, and oxalate. However, no direct binding was detected for the remaining three carboxylates that elicited an McpT-dependent chemotaxis response. Taken together, these results demonstrate that McpT is a broad-range carboxylate chemoreceptor that mediates chemotactic response via direct ligand binding and an indirect mechanism that needs to be identified. IMPORTANCE Nitrate pollution is one of the most widespread and challenging environmental problems that is mainly caused by the agricultural overapplication of nitrogen fertilizers. Biological nitrogen fixation by the endosymbiont Sinorhizobium meliloti enhances the growth of its host Medicago sativa (alfalfa), which also efficiently supplies the soil with nitrogen. Establishment of the S. meliloti-alfalfa symbiosis relies on the early exchange and recognition of chemical signals. The present study contributes to the disclosure of this complex molecular dialogue by investigating the underlying mechanisms of carboxylate sensing in S. meliloti. Understanding individual steps that govern the S. meliloti-alfalfa molecular cross talk helps in the development of efficient, commercial bacterial inoculants that promote the growth of alfalfa, which is the most cultivated forage legume in the world, and improves soil fertility.

Nonadditive Transcriptomic Signatures of Genotype-by-Genotype Interactions during the Initiation of Plant-Rhizobium Symbiosis

Rhizobia are ecologically important, facultative plant-symbiotic microbes. In nature, there is a large variability in the association of rhizobial strains and host plants of the same species. Here, we evaluated whether plant and rhizobial genotypes influence the initial transcriptional response of rhizobium following perception of a host plant. RNA sequencing of the model rhizobium Sinorhizobium meliloti exposed to root exudates or luteolin (an inducer of nod genes, involved in the early steps of symbiotic interaction) was performed on a combination of three S. meliloti strains and three alfalfa varieties as host plants. The response to root exudates involved hundreds of changes in the rhizobium transcriptome. Of the differentially expressed genes, 35% were influenced by the strain genotype, 16% were influenced by the plant genotype, and 29% were influenced by strain-by-host plant genotype interactions. We also examined the response of a hybrid S. meliloti strain in which the symbiotic megaplasmid (∼20% of the genome) was mobilized between two of the above-mentioned strains. Dozens of genes were upregulated in the hybrid strain, indicative of nonadditive variation in the transcriptome. In conclusion, this study demonstrated that transcriptional responses of rhizobia upon perception of legumes are influenced by the genotypes of both symbiotic partners and their interaction, suggesting a wide spectrum of genetic determinants involved in the phenotypic variation of plant-rhizobium symbiosis.IMPORTANCE A sustainable way for meeting the need of an increased global food demand should be based on a holobiont perspective, viewing crop plants as intimately associated with their microbiome, which helps improve plant nutrition, tolerance to pests, and adverse climate conditions. However, the genetic repertoire needed for efficient association with plants by the microbial symbionts is still poorly understood. The rhizobia are an exemplary model of facultative plant symbiotic microbes. Here, we evaluated whether genotype-by-genotype interactions could be identified in the initial transcriptional response of rhizobium perception of a host plant. We performed an RNA sequencing study to analyze the transcriptomes of different rhizobial strains elicited by root exudates of three alfalfa varieties as a proxy of an early step of the symbiotic interaction. The results indicated strain- and plant variety-dependent variability in the observed transcriptional changes, providing fundamentally novel insights into the genetic basis of rhizobium-plant interactions. Our results provide genetic insights and perspective to aid in the exploitation of natural rhizobium variation for improvement of legume growth in agricultural ecosystems.

Nonadditive Transcriptomic Signatures of Genotype-by-Genotype Interactions during the Initiation of Plant-Rhizobium Symbiosis

Rhizobia are ecologically important, facultative plant-symbiotic microbes. In nature, there is a large variability in the association of rhizobial strains and host plants of the same species. Here, we evaluated whether plant and rhizobial genotypes influence the initial transcriptional response of rhizobium following perception of a host plant. RNA sequencing of the model rhizobium Sinorhizobium meliloti exposed to root exudates or luteolin (an inducer of nod genes, involved in the early steps of symbiotic interaction) was performed on a combination of three S. meliloti strains and three alfalfa varieties as host plants. The response to root exudates involved hundreds of changes in the rhizobium transcriptome. Of the differentially expressed genes, 35% were influenced by the strain genotype, 16% were influenced by the plant genotype, and 29% were influenced by strain-by-host plant genotype interactions. We also examined the response of a hybrid S. meliloti strain in which the symbiotic megaplasmid (∼20% of the genome) was mobilized between two of the above-mentioned strains. Dozens of genes were upregulated in the hybrid strain, indicative of nonadditive variation in the transcriptome. In conclusion, this study demonstrated that transcriptional responses of rhizobia upon perception of legumes are influenced by the genotypes of both symbiotic partners and their interaction, suggesting a wide spectrum of genetic determinants involved in the phenotypic variation of plant-rhizobium symbiosis.IMPORTANCE A sustainable way for meeting the need of an increased global food demand should be based on a holobiont perspective, viewing crop plants as intimately associated with their microbiome, which helps improve plant nutrition, tolerance to pests, and adverse climate conditions. However, the genetic repertoire needed for efficient association with plants by the microbial symbionts is still poorly understood. The rhizobia are an exemplary model of facultative plant symbiotic microbes. Here, we evaluated whether genotype-by-genotype interactions could be identified in the initial transcriptional response of rhizobium perception of a host plant. We performed an RNA sequencing study to analyze the transcriptomes of different rhizobial strains elicited by root exudates of three alfalfa varieties as a proxy of an early step of the symbiotic interaction. The results indicated strain- and plant variety-dependent variability in the observed transcriptional changes, providing fundamentally novel insights into the genetic basis of rhizobium-plant interactions. Our results provide genetic insights and perspective to aid in the exploitation of natural rhizobium variation for improvement of legume growth in agricultural ecosystems.

Symbiotic and Nonsymbiotic Members of the Genus Ensifer (syn. Sinorhizobium) Are Separated into Two Clades Based on Comparative Genomics and High-Throughput Phenotyping

Rhizobium–legume symbioses serve as paradigmatic examples for the study of mutualism evolution. The genus Ensifer (syn. Sinorhizobium) contains diverse plant-associated bacteria, a subset of which can fix nitrogen in symbiosis with legumes. To gain insights into the evolution of symbiotic nitrogen fixation (SNF), and interkingdom mutualisms more generally, we performed extensive phenotypic, genomic, and phylogenetic analyses of the genus Ensifer. The data suggest that SNF likely emerged several times within the genus Ensifer through independent horizontal gene transfer events. Yet, the majority (105 of 106) of the Ensifer strains with the nodABC and nifHDK nodulation and nitrogen fixation genes were found within a single, monophyletic clade. Comparative genomics highlighted several differences between the “symbiotic” and “nonsymbiotic” clades, including divergences in their pangenome content. Additionally, strains of the symbiotic clade carried 325 fewer genes, on average, and appeared to have fewer rRNA operons than strains of the nonsymbiotic clade. Initial characterization of a subset of ten Ensifer strains identified several putative phenotypic differences between the clades. Tested strains of the nonsymbiotic clade could catabolize 25% more carbon sources, on average, than strains of the symbiotic clade, and they were better able to grow in LB medium and tolerate alkaline conditions. On the other hand, the tested strains of the symbiotic clade were better able to tolerate heat stress and acidic conditions. We suggest that these data support the division of the genus Ensifer into two main subgroups, as well as the hypothesis that pre-existing genetic features are required to facilitate the evolution of SNF in bacteria.

Genome-scale metabolic reconstruction of the symbiosis between a leguminous plant and a nitrogen-fixing bacterium

The mutualistic association between leguminous plants and endosymbiotic rhizobial bacteria is a paradigmatic example of a symbiosis driven by metabolic exchanges. Here, we report the reconstruction and modelling of a genome-scale metabolic network of Medicago truncatula (plant) nodulated by Sinorhizobium meliloti (bacterium). The reconstructed nodule tissue contains five spatially distinct developmental zones and encompasses the metabolism of both the plant and the bacterium. Flux balance analysis (FBA) suggests that the metabolic costs associated with symbiotic nitrogen fixation are primarily related to supporting nitrogenase activity, and increasing N2-fixation efficiency is associated with diminishing returns in terms of plant growth. Our analyses support that differentiating bacteroids have access to sugars as major carbon sources, ammonium is the main nitrogen export product of N2-fixing bacteria, and N2 fixation depends on proton transfer from the plant cytoplasm to the bacteria through acidification of the peribacteroid space. We expect that our model, called 'Virtual Nodule Environment' (ViNE), will contribute to a better understanding of the functioning of legume nodules, and may guide experimental studies and engineering of symbiotic nitrogen fixation.

Diversity and phenotypic analyses of salt- and heat-tolerant wild bean Phaseolus filiformis rhizobia native of a sand beach in Baja California and description of Ensifer aridi sp. nov

In northern Mexico, aridity, salinity and high temperatures limit areas that can be cultivated. To investigate the nature of nitrogen-fixing symbionts of Phaseolus filiformis, an adapted wild bean species native to this region, their phylogenies were inferred by MLSA. Most rhizobia recovered belong to the proposed new species Ensifer aridi. Phylogenetic analyses of nodC and nifH show that Mexican isolates carry symbiotic genes acquired through horizontal gene transfer that are divergent from those previously characterized among bean symbionts. These strains are salt tolerant, able to grow in alkaline conditions, high temperatures, and capable of utilizing a wide range of carbohydrates and organic acids as carbon sources for growth. This study improves the knowledge on diversity, geographic distribution and evolution of bean-nodulating rhizobia in Mexico and further enlarges the spectrum of microsymbiont with which Phaseolus species can interact with, including cultivated bean varieties, notably under stressed environments. Here, the species Ensifer aridi sp. nov. is proposed as strain type of the Moroccan isolate LMR001^T (= LMG 31426^T; = HAMBI 3707^T) recovered from desert sand dune.

Municipal solid wastes as a resource for environmental recovery: Impact of water treatment residuals and compost on the microbial and biochemical features of As and trace metal-polluted soils

In this study we evaluated the microbiological and biochemical impact of iron-based water treatment residuals (Fe-WTRs) and municipal solid waste compost (MSWC), alone and combined, on three different soils co-contaminated with arsenic (As) and trace-metals (TM), i.e. Pb, Cu and Zn. Overall, all the amendments considered significantly increased the abundance of culturable heterotrophic bacteria, with MSWC showing the greatest impact across all soils (up to a 24% increase). In most of treated soils this was accompanied by a significant reduction of both the (culturable) fungal/bacterial ratio, and the proportion of culturable As(V)- and As(III)-resistant bacteria with respect to total bacterial population. The catabolic potential and versatility of the resident microbial communities (assessed by community level physiological profile) was highly soil-dependent and substantial increases of both parameters were observed in the amended soils with the higher total As concentration (from approx. 749 to 22,600 mg kg^-1). Moreover, both carbon source utilisation profile and 16S rRNA soil metagenome sequencing indicated a significant impact of MSWC and Fe-WTRs on the structure and diversity of soil microbial communities, with Proteobacteria, Actinobacteria and Firmicutes being the most affected taxa. The assessment of selected soil enzyme activities (dehydrogenase, urease and β-glucosidase) indicated an increase of metabolic functioning especially in soils treated with MSWC (e.g. dehydrogenase activity increased up to 19.5-fold in the most contaminated soil treated with MSWC). Finally, the microbial and biochemical features of treated (and untreated) contaminated soils (i.e. total bacterial counts, catabolic potential and versatility and soil enzyme activities) were highly correlated with the concentrations of labile As and TM in these latter soils and supported a clear role of the tested amendments (especially MSWC) as As- and TM-immobilising agents.

Sinorhizobium meliloti Chemoreceptor McpV Senses Short-Chain Carboxylates via Direct Binding

Sinorhizobium meliloti is a soil-dwelling endosymbiont of alfalfa that has eight chemoreceptors to sense environmental stimuli during its free-living state. The functions of two receptors have been characterized, with McpU and McpX serving as general amino acid and quaternary ammonium compound sensors, respectively. Both receptors use a dual Cache (calcium channels and chemotaxis receptors) domain for ligand binding. We identified that the ligand-binding periplasmic region (PR) of McpV contains a single Cache domain. Homology modeling revealed that McpV^PR is structurally similar to a sensor domain of a chemoreceptor with unknown function from Anaeromyxobacter dehalogenans, which crystallized with acetate in its binding pocket. We therefore assayed McpV for carboxylate binding and S. meliloti for carboxylate sensing. Differential scanning fluorimetry identified 10 potential ligands for McpV^PR Nine of these are monocarboxylates with chain lengths between two and four carbons. We selected seven compounds for capillary assay analysis, which established positive chemotaxis of the S. meliloti wild type, with concentrations of peak attraction at 1 mM for acetate, propionate, pyruvate, and glycolate, and at 100 mM for formate and acetoacetate. Deletion of mcpV or mutation of residues essential for ligand coordination abolished positive chemotaxis to carboxylates. Using microcalorimetry, we determined that dissociation constants of the seven ligands with McpV^PR were in the micromolar range. An McpV^PR variant with a mutation in the ligand coordination site displayed no binding to isobutyrate or propionate. Of all the carboxylates tested as attractants, only glycolate was detected in alfalfa seed exudates. This work examines the relevance of carboxylates and their sensor to the rhizobium-legume interaction.IMPORTANCE Legumes share a unique association with certain soil-dwelling bacteria known broadly as rhizobia. Through concerted interorganismal communication, a legume allows intracellular infection by its cognate rhizobial species. The plant then forms an organ, the root nodule, dedicated to housing and supplying fixed carbon and nutrients to the bacteria. In return, the engulfed rhizobia, differentiated into bacteroids, fix atmospheric N₂ into ammonium for the plant host. This interplay is of great benefit to the cultivation of legumes, such as alfalfa and soybeans, and is initiated by chemotaxis to the host plant. This study on carboxylate chemotaxis contributes to the understanding of rhizobial survival and competition in the rhizosphere and aids the development of commercial inoculants.

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.

Techniques for Large-Scale Bacterial Genome Manipulation and Characterization of the Mutants with Respect to In Silico Metabolic Reconstructions

The rate at which all genes within a bacterial genome can be identified far exceeds the ability to characterize these genes. To assist in associating genes with cellular functions, a large-scale bacterial genome deletion approach can be employed to rapidly screen tens to thousands of genes for desired phenotypes. Here, we provide a detailed protocol for the generation of deletions of large segments of bacterial genomes that relies on the activity of a site-specific recombinase. In this procedure, two recombinase recognition target sequences are introduced into known positions of a bacterial genome through single cross-over plasmid integration. Subsequent expression of the site-specific recombinase mediates recombination between the two target sequences, resulting in the excision of the intervening region and its loss from the genome. We further illustrate how this deletion system can be readily adapted to function as a large-scale in vivo cloning procedure, in which the region excised from the genome is captured as a replicative plasmid. We next provide a procedure for the metabolic analysis of bacterial large-scale genome deletion mutants using the Biolog Phenotype MicroArray™ system. Finally, a pipeline is described, and a sample Matlab script is provided, for the integration of the obtained data with a draft metabolic reconstruction for the refinement of the reactions and gene-protein-reaction relationships in a metabolic reconstruction.

A Key Regulator of the Glycolytic and Gluconeogenic Central Metabolic Pathways in Sinorhizobium meliloti

The order Rhizobiales contains numerous agriculturally, biotechnologically, and medically important bacteria, including the rhizobia, and the genera Agrobacterium, Brucella, and Methylobacterium, among others. These organisms tend to be metabolically versatile, but there has been relatively little investigation into the regulation of their central carbon metabolic pathways. Here, RNA-sequencing and promoter fusion data are presented to show that the PckR protein is a key regulator of central carbon metabolism in Sinorhizobium meliloti; during growth with gluconeogenic substrates, PckR represses expression of the complete Entner-Doudoroff glycolytic pathway and induces expression of the pckA and fbaB gluconeogenic genes. Electrophoretic mobility shift assays indicate that PckR binds an imperfect palindromic sequence that overlaps the promoter or transcriptional start site in the negatively regulated promoters, or is present in tandem upstream the promoter motifs in the positively regulated promoters. Genetic and in vitro electrophoretic mobility shift assay experiments suggest that elevated concentrations of a PckR effector ligand results in the dissociation of PckR from its target binding site, and evidence is presented that suggests phosphoenolpyruvate may function as the effector. Characterization of missense pckR alleles identified three conserved residues important for increasing the affinity of PckR for its cognate effector molecule. Bioinformatics analyses illustrates that PckR is limited to a narrow phylogenetic range consisting of the Rhizobiaceae, Phyllobacteriaceae, Brucellaceae, and Bartonellaceae families. These data provide novel insights into the regulation of the core carbon metabolic pathways of this pertinent group of α-proteobacteria.

Insights from the complete genome sequence of Clostridium tyrobutyricum provide a platform for biotechnological and industrial applications

Genetic research enables the evolution of novel biochemical reactions for the production of valuable chemicals from environmentally-friendly raw materials. However, the choice of appropriate microorganisms to support these reactions, which must have strong robustness and be capable of a significant product output, is a major difficulty. In the present study, the complete genome of the Clostridium tyrobutyricum strain CCTCC W428, a hydrogen- and butyric acid-producing bacterium with increased oxidative tolerance was analyzed. A total length of 3,011,209 bp of the C. tyrobutyricum genome with a GC content of 31.04% was assembled, and 3038 genes were discovered. Furthermore, a comparative clustering of proteins from C. tyrobutyricum CCTCC W428, C. acetobutylicum ATCC 824, and C. butyricum KNU-L09 was conducted. The results of genomic analysis indicate that butyric acid is produced by CCTCC W428 from butyryl-CoA through acetate reassimilation via CoA transferase, instead of the well-established phosphotransbutyrylase-butyrate kinase pathway. In addition, we identified ten proteins putatively involved in hydrogen production and 21 proteins associated with CRISPR systems, together with 358 ORFs related to ABC transporters and transcriptional regulators. Enzymes, such as oxidoreductases, HNH endonucleases, and catalase, were also found in this species. The genome sequence illustrates that C. tyrobutyricum has several desirable traits, and is expected to be suitable as a platform for the high-level production of bulk chemicals as well as bioenergy.

Bacterial Root Microbiome of Plants Growing in Oil Sands Reclamation Covers

Oil sands mining in northern Alberta impacts a large footprint, but the industry is committed to reclaim all disturbed land to an ecologically healthy state in response to environmental regulations. However, these newly reconstructed landscapes may be limited by several factors that include low soil nutrient levels and reduced microbial activity. Rhizosphere microorganisms colonize plant roots providing hosts with nutrients, stimulating growth, suppressing disease and increasing tolerance to abiotic stress. High-throughput sequencing techniques can be used to provide a detailed characterization of microbial community structure. This study used 16S rRNA amplicon sequencing to characterize the bacterial root microbiome associated with annual barley (Hordeum vulgare) and sweet clover (Melilotus albus) growing in an oil sands reclamation area. Our results indicate that Proteobacteria dominated the endosphere, whereas other phyla such as Acidobacteria and Gemmatimonadetes were restricted to the rhizosphere, suggesting that plants have the ability to select for certain soil bacterial consortia. The bacterial community in the endosphere compartments were less rich and diverse compared to the rhizosphere. Furthermore, it was apparent that sweet clover plants were more selective, as the community exhibited a lower richness and diversity compared to barley. Members of the family Rhizobiaceae, such as Sinorhizobium and Rhizobium were mainly associated with clover, whereas Acholeplasma (wall-less bacteria transmitted by insects) was unique to barley. Genera from the Enterobacteriaceae family, such as Yersinia and Lentzea were also mostly detected in barley, while other genera such Pseudomonas and Pantoea were able to successfully colonize both plants. Endophytic bacterial profiles varied within the same plant species at different sampling locations; however, these differences were driven by factors other than slope positions or cover management. Our results suggest that bacterial endophytic communities of plants growing in land reclamation systems are a subset of the rhizosphere community and selection is driven by plant factors.

Using phenotype microarrays in the assessment of the antibiotic susceptibility profile of bacteria isolated from wastewater in on-site treatment facilities

The scope of the study was to apply Phenotype Biolog MicroArray (PM) technology to test the antibiotic sensitivity of the bacterial strains isolated from on-site wastewater treatment facilities. In the first step of the study, the percentage values of resistant bacteria from total heterotrophic bacteria growing on solid media supplemented with various antibiotics were determined. In the untreated wastewater, the average shares of kanamycin-, streptomycin-, and tetracycline-resistant bacteria were 53, 56, and 42%, respectively. Meanwhile, the shares of kanamycin-, streptomycin-, and tetracycline-resistant bacteria in the treated wastewater were 39, 33, and 29%, respectively. To evaluate the antibiotic susceptibility of the bacteria present in the wastewater, using the phenotype microarrays (PMs), the most common isolates from the treated wastewater were chosen: Serratia marcescens ss marcescens, Pseudomonas fluorescens, Stenotrophomonas maltophilia, Stenotrophomonas rhizophila, Microbacterium flavescens, Alcaligenes faecalis ss faecalis, Flavobacterium hydatis, Variovorax paradoxus, Acinetobacter johnsonii, and Aeromonas bestiarum. The strains were classified as multi-antibiotic-resistant bacteria. Most of them were resistant to more than 30 antibiotics from various chemical classes. Phenotype microarrays could be successfully used as an additional tool for evaluation of the multi-antibiotic resistance of environmental bacteria and in preliminary determination of the range of inhibition concentration.

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.

The genome of some bacteria consists of two or more chromosomes or replicons. Here, diCenzo et al. integrate genome-scale metabolic modelling and growth data from a collection of mutants of the plant symbiont Sinorhizobium meliloti to estimate the fitness contribution of each replicon in three environments.

Tailoring of global transcription sigma D factor by random mutagenesis to improve Escherichia coli tolerance towards low-pHs

Bioconversion processes of organic acid or acid hydrolysis of raw material for microbial metabolism often suffer limitations as a result of microbial sensitivity in low-pH conditions. We adopted a three-step method called RAndom Insertional-deletional Strand Exchange mutagenesis (RAISE) to engineer the components of global regulator Sigma D factor (RpoD) of Escherichia coli to improve its acid tolerance. The best strain Mutant VII was identified from random mutagenesis libraries based on the growth performance, which exhibited much higher growth rate than the control (0.22h(-1) vs. 0.15h(-1)) at pH as low as 3.17. Combined transcriptome and phenome analysis of E. coli was carried out to better understand the global effects of RpoD on the regulatory networks. Our analysis showed that 95 (2.1%) of all E. coli genes were induced and 178 (4.0%) genes were repressed, including those for trehalose biosynthesis, nucleotides biosynthesis, carbon metabolism, amino acid utilization, except for acid resistance. Also regulated were the master regulators (ArcA, EvgA, H-NS and RpoS) and gene/operon-specific transcription factors (GadX, GadW, AppY, YdeO, KdgR). These results demonstrated that RpoD acts as global regulator in the growth phase of E. coli and consequently improves acid tolerances.

Divergent genes in potential inoculant Sinorhizobium strains are related to DNA replication, recombination, and repair

To serve as inoculants of legumes, nitrogen-fixing rhizobium strains should be competitive and tolerant of diverse environments. We hybridized the genomes of symbiotically efficient and salt tolerant Sinorhizobium inoculant strains onto the Sinorhizobium meliloti Rm1021 microarray. The number of variable genes, that is, divergent or putatively multiplied genes, ranged from 503 to 1556 for S. meliloti AK23, S. meliloti STM 1064 and S. arboris HAMBI 1552. The numbers of divergent genes affiliated with the symbiosis plasmid pSymA and related to DNA replication, recombination and repair were significantly higher than expected. The variation was mainly in the accessory genome, implying that it was important in shaping the adaptability of the strains.

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.

Genome and Phenotype Microarray Analyses of Rhodococcus sp. BCP1 and Rhodococcus opacus R7: Genetic Determinants and Metabolic Abilities with Environmental Relevance

In this paper comparative genome and phenotype microarray analyses of Rhodococcus sp. BCP1 and Rhodococcus opacus R7 were performed. Rhodococcus sp. BCP1 was selected for its ability to grow on short-chain n-alkanes and R. opacus R7 was isolated for its ability to grow on naphthalene and on o-xylene. Results of genome comparison, including BCP1, R7, along with other Rhodococcus reference strains, showed that at least 30% of the genome of each strain presented unique sequences and only 50% of the predicted proteome was shared. To associate genomic features with metabolic capabilities of BCP1 and R7 strains, hundreds of different growth conditions were tested through Phenotype Microarray, by using Biolog plates and plates manually prepared with additional xenobiotic compounds. Around one-third of the surveyed carbon sources was utilized by both strains although R7 generally showed higher metabolic activity values compared to BCP1. Moreover, R7 showed broader range of nitrogen and sulphur sources. Phenotype Microarray data were combined with genomic analysis to genetically support the metabolic features of the two strains. The genome analysis allowed to identify some gene clusters involved in the metabolism of the main tested xenobiotic compounds. Results show that R7 contains multiple genes for the degradation of a large set of aromatic and PAHs compounds, while a lower variability in terms of genes predicted to be involved in aromatic degradation was found in BCP1. This genetic feature can be related to the strong genetic pressure exerted by the two different environment from which the two strains were isolated. According to this, in the BCP1 genome the smo gene cluster involved in the short-chain n-alkanes degradation, is included in one of the unique regions and it is not conserved in the Rhodococcus strains compared in this work. Data obtained underline the great potential of these two Rhodococcus spp. strains for biodegradation and environmental decontamination processes.

Genomic characterization of Sinorhizobium meliloti AK21, a wild isolate from the Aral Sea Region

The symbiotic, nitrogen-fixing bacterium Sinorhizobium meliloti has been widely studied due to its ability to improve crop yields through direct interactions with leguminous plants. S. meliloti AK21 is a wild type strain that forms nodules on Medicago plants in saline and drought conditions in the Aral Sea Region. The aim of this work was to establish the genetic similarities and differences between S. meliloti AK21 and the reference strain S. meliloti 1021. Comparative genome hybridization with the model reference strain S. meliloti 1021 yielded 365 variable genes, grouped into 11 regions in the three main replicons in S. meliloti AK21. The most extensive regions of variability were found in the symbiotic plasmid pSymA, which also contained the largest number of orthologous and polymorphic sequences identified by suppression subtractive hybridization. This procedure identified a large number of divergent sequences and others without homology in the databases, the further investigation of which could provide new insight into the alternative metabolic pathways present in S. meliloti AK21. We identified a plasmid replication module from the repABC replicon family, together with plasmid mobilization-related genes (traG and a VirB9-like protein), which suggest that this indigenous isolate harbors an accessory plasmid. Furthermore, the transcriptomic profiles reflected differences in gene content and regulation between S. meliloti AK21 and S. meliloti 1021 (ExpR and PhoB regulons), but provided evidence for an as yet unknown, alternative mechanism involving activation of the cbb3 terminal oxidase. Finally, phenotypic microarrays characterization revealed a greater versatility of substrate use and chemical degradation than for S. meliloti 1021.

The integrated microbial genome resource of analysis

Integrated Microbial Genomes and Metagenomes (IMG) is a biocomputational system that allows to provide information and support for annotation and comparative analysis of microbial genomes and metagenomes. IMG has been developed by the US Department of Energy (DOE)-Joint Genome Institute (JGI). IMG platform contains both draft and complete genomes, sequenced by Joint Genome Institute and other public and available genomes. Genomes of strains belonging to Archaea, Bacteria, and Eukarya domains are present as well as those of viruses and plasmids. Here, we provide some essential features of IMG system and case study for pangenome analysis.

From pangenome to panphenome and back

The ability to relate genomic differences in bacterial species to their variability in expressed phenotypes is one of the most challenging tasks in today's biology. Such task is of paramount importance towards the understanding of biotechnologically relevant pathways and possibly for their manipulation. Fundamental prerequisites are the genome-wide reconstruction of metabolic pathways and a comprehensive measurement of cellular phenotypes. Cellular pathways can be reliably reconstructed using the KEGG database, while the OmniLog™ Phenotype Microarray (PM) technology may be used to measure nearly 2,000 growth conditions over time. However, few computational tools that can directly link PM data with the gene(s) of interest followed by the extraction of information on gene-phenotype correlation are available. In this chapter the use of the DuctApe software suite is presented, which allows the joint analysis of bacterial genomic and phenomic data, highlighting those pathways and reactions most probably associated with phenotypic variability. A case study on four Sinorhizobium meliloti strains is presented; more example datasets are available online.

Examination of prokaryotic multipartite genome evolution through experimental genome reduction

Many bacteria carry two or more chromosome-like replicons. This occurs in pathogens such as Vibrio cholerea and Brucella abortis as well as in many N₂-fixing plant symbionts including all isolates of the alfalfa root-nodule bacteria Sinorhizobium meliloti. Understanding the evolution and role of this multipartite genome organization will provide significant insight into these important organisms; yet this knowledge remains incomplete, in part, because technical challenges of large-scale genome manipulations have limited experimental analyses. The distinct evolutionary histories and characteristics of the three replicons that constitute the S. meliloti genome (the chromosome (3.65 Mb), pSymA megaplasmid (1.35 Mb), and pSymB chromid (1.68 Mb)) makes this a good model to examine this topic. We transferred essential genes from pSymB into the chromosome, and constructed strains that lack pSymB as well as both pSymA and pSymB. This is the largest reduction (45.4%, 3.04 megabases, 2866 genes) of a prokaryotic genome to date and the first removal of an essential chromid. Strikingly, strains lacking pSymA and pSymB (ΔpSymAB) lost the ability to utilize 55 of 74 carbon sources and various sources of nitrogen, phosphorous and sulfur, yet the ΔpSymAB strain grew well in minimal salts media and in sterile soil. This suggests that the core chromosome is sufficient for growth in a bulk soil environment and that the pSymA and pSymB replicons carry genes with more specialized functions such as growth in the rhizosphere and interaction with the plant. These experimental data support a generalized evolutionary model, in which non-chromosomal replicons primarily carry genes with more specialized functions. These large secondary replicons increase the organism's niche range, which offsets their metabolic burden on the cell (e.g. pSymA). Subsequent co-evolution with the chromosome then leads to the formation of a chromid through the acquisition of functions core to all niches (e.g. pSymB).

Rhizobia are free-living bacteria of agricultural and environmental importance that form root-nodules on leguminous plants and provide these plants with fixed nitrogen. Many of the rhizobia have a multipartite genome, as do several plant and animal pathogens. All isolates of the alfalfa symbiont, Sinorhizobium meliloti, carry three large replicons, the chromosome (∼3.7 Mb), pSymA megaplasmid (∼1.4 Mb), and pSymB chromid (∼1.7 Mb). To gain insight into the role and evolutionary history of these replicons, we have ‘reversed evolution’ by constructing a S. meliloti strain consisting solely of the chromosome and lacking the pSymB chromid and pSymA megaplasmid. As the resulting strain was viable, we could perform a detailed phenotypic analysis and these data provided significant insight into the biology and metabolism of S. meliloti. The data lend direct experimental evidence in understanding the evolution and role of the multipartite genome. Specifically the large secondary replicons increase the organism's niche range, and this advantage offsets the metabolic burden of these replicons on the cell. Additionally, the single-chromosome strain offers a useful platform to facilitate future forward genetic approaches to understanding and manipulating the symbiosis and plant-microbe interactions.

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.

Cell growth inhibition upon deletion of four toxin-antitoxin loci from the megaplasmids of Sinorhizobium meliloti

Toxin and antitoxin (TA) gene pairs are addiction systems that are present in many microbial genomes. Sinorhizobium meliloti is an N2-fixing bacterial symbiont of alfalfa and other leguminous plants, and its genome consists of three large replicons, a circular chromosome (3.7 Mb) and the megaplasmids pSymA (1.4 Mb) and pSymB (1.7 Mb). S. meliloti carries 211 predicted type II TA genes, each encoding a toxin or an antitoxin. We constructed defined deletion strains that collectively removed the entire pSymA and pSymB megaplasmids except for their oriV regions. Of approximately 100 TA genes on pSymA and pSymB, we identified four whose loss was associated with cell death or stasis unless copies of the genes were supplied in trans. Orthologs of three of these loci have been characterized in other organisms (relB/E [sma0471/sma0473], Fic [DOC] [sma2105], and VapC [PIN] [orf2230/sma2231]), and this report contains the first experimental proof that RES/Xre (smb21127/smb21128) loci can function as a TA system. Transcriptome sequencing (RNA-seq) analysis did not reveal transcriptional differences between the TA systems to account for why deletion of the four "active" systems resulted in cell toxicity. These data suggest that severe cell growth phenotypes result from the loss of a few TA systems and that loss of most TA systems may result in more subtle phenotypes. These four TA systems do not appear to play a direct role in the S. meliloti-alfalfa symbiosis, as strains lacking these TA systems had a symbiotic N2 fixation phenotype that was indistinguishable from the wild type.

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/.

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.

Polyphasic characterization of rhizobia isolated from Leucaena leucocephala from Panxi, China

Leucaena leucocephala was introduced into Panxi, Sichuan, China, in the 1980s and 1990s for afforestation and preventing water loss and soil erosion in this area. The co-introduction of rhizobial symbionts of introduced plants has drawn attention since they may influence local soil communities. We studied the phylogenetic position of the L. leucocephala isolates and assessed if the rhizobia were introduced together with the host to Panxi, Sichuan, China. The glnII and atpD genes of fifteen representative isolates were sequenced and analyzed, and applied multilocus sequence analyses in which the housekeeping genes recA, glnII and atpD were included. Furthermore, we estimated the within species diversity directly with 23S rDNA and IGS RFLP and indirectly through phenotypic analysis of forty L. leucocephala isolates. The isolates represented seven species and 38 diversified strains in the genera Ensifer, Mesorhizobium, Bradyrhizobium and Rhizobium. The within species diversity of the Ensifer isolates was large, proposing a potential to occupy novel niches. There was not conclusive evidence to show that any of the strains would have been co-introduced with L. leucocephala. On the contrary, we came to a conclusion that the possible introduction should not be inferred from sequence data alone.

Phenotype microarray screening of carbon sources used by Vibrio cholerae identifies genes regulated by the cAMP receptor protein

The cyclic AMP receptor protein (CRP) regulates genes involved in carbon source metabolism, iron uptake, and virulence in bacteria. Identifying the carbon sources utilized by bacteria that are regulated by CRP will help elucidate the CRP regulation cascade and associated responses to environmental stimuli. CRP-dependent regulation of carbon source metabolism in Vibrio cholerae is not thoroughly understood. To identify the candidate carbon sources utilized by V. cholerae that are affected by CRP, we used high-throughput screening to compare the metabolic differences between wild-type and CRP mutant strains of V. cholerae O1 El Tor. Phenotype microarray was used for primary screening of the wild-type and mutant strains, followed by minimal media growth assays and quantitative RT-PCR to validate the candidate carbon sources. In total, 24 carbon sources were subject to CRP regulation, 11 of which have not been previously reported in bacteria. The genes known to be involved in the metabolism of 4 of the carbon sources identified were verified by quantitative RT-PCR. In addition, gel shift experiments showed that CRP bound directly to VCA0053 and VC0391 promoters. Overall, this comprehensive analysis of CRP-mediated catabolite control in V. cholerae has identified new candidate carbon sources for in-depth experimental studies.

Explaining microbial phenotypes on a genomic scale: GWAS for microbes

There is an increasing availability of complete or draft genome sequences for microbial organisms. These data form a potentially valuable resource for genotype-phenotype association and gene function prediction, provided that phenotypes are consistently annotated for all the sequenced strains. In this review, we address the requirements for successful gene-trait matching. We outline a basic protocol for microbial functional genomics, including genome assembly, annotation of genotypes (including single nucleotide polymorphisms, orthologous groups and prophages), data pre-processing, genotype-phenotype association, visualization and interpretation of results. The methodologies for association described herein can be applied to other data types, opening up possibilities to analyze transcriptome-phenotype associations, and correlate microbial population structure or activity, as measured by metagenomics, to environmental parameters.

The success of acinetobacter species; genetic, metabolic and virulence attributes

An understanding of why certain Acinetobacter species are more successful in causing nosocomial infections, transmission and epidemic spread in healthcare institutions compared with other species is lacking. We used genomic, phenotypic and virulence studies to identify differences between Acinetobacter species. Fourteen strains representing nine species were examined. Genomic analysis of six strains showed that the A. baumannii core genome contains many genes important for diverse metabolism and survival in the host. Most of the A. baumannii core genes were also present in one or more of the less clinically successful species. In contrast, when the accessory genome of an individual A. baumannii strain was compared to a strain of a less successful species (A. calcoaceticus RUH2202), many operons with putative virulence function were found to be present only in the A. baumannii strain, including the csu operon, the acinetobactin chromosomal cluster, and bacterial defence mechanisms. Phenotype microarray analysis showed that compared to A. calcoaceticus (RUH2202), A. baumannii ATCC 19606^T was able to utilise nitrogen sources more effectively and was more tolerant to pH, osmotic and antimicrobial stress. Virulence differences were also observed, with A. baumannii ATCC 19606^T, A. pittii SH024, and A. nosocomialis RUH2624 persisting and forming larger biofilms on human skin than A. calcoaceticus. A. baumannii ATCC 19606^T and A. pittii SH024 were also able to survive in a murine thigh infection model, whereas the other two species were eradicated. The current study provides important insights into the elucidation of differences in clinical relevance among Acinetobacter species.

Biolog phenotype microarrays

Phenotype microarrays nicely complement traditional genomic, transcriptomic, and proteomic analysis by offering opportunities for researchers to ground microbial systems analysis and modeling in a broad yet quantitative assessment of the organism's physiological response to different metabolites and environments. Biolog phenotype assays achieve this by coupling tetrazolium dyes with minimally defined nutrients to measure the impact of hundreds of carbon, nitrogen, phosphorous, and sulfur sources on redox reactions that result from compound-induced effects on the electron transport chain. Over the years, we have used Biolog's reproducible and highly sensitive assays to distinguish closely related bacterial isolates, to understand their metabolic differences, and to model their metabolic behavior using flux balance analysis. This chapter describes Biolog phenotype microarray system components, reagents, and methods, particularly as they apply to bacterial identification, characterization, and metabolic analysis.

Cryptic prophages help bacteria cope with adverse environments

Phages are the most abundant entity in the biosphere and outnumber bacteria by a factor of 10. Phage DNA may also constitute 20% of bacterial genomes; however, its role is ill defined. Here, we explore the impact of cryptic prophages on cell physiology by precisely deleting all nine prophage elements (166 kbp) using Escherichia coli. We find that cryptic prophages contribute significantly to resistance to sub-lethal concentrations of quinolone and β-lactam antibiotics primarily through proteins that inhibit cell division (for example, KilR of rac and DicB of Qin). Moreover, the prophages are beneficial for withstanding osmotic, oxidative and acid stresses, for increasing growth, and for influencing biofilm formation. Prophage CPS-53 proteins YfdK, YfdO and YfdS enhanced resistance to oxidative stress, prophages e14, CPS-53 and CP4-57 increased resistance to acid, and e14 and rac proteins increased early biofilm formation. Therefore, cryptic prophages provide multiple benefits to the host for surviving adverse environmental conditions.

Up to 20% of bacterial genomes are made up of cryptic prophages, but their function is relatively unknown. In this study, the authors demonstrate that prophages influence the response of the host cell to stress and provide a competitive growth advantage in the presence of antibiotics.

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.

Role of the Sinorhizobium meliloti global regulator Hfq in gene regulation and symbiosis

The RNA-binding protein Hfq is a global regulator which controls diverse cellular processes in bacteria. To begin understanding the role of Hfq in the Sinorhizobium meliloti-Medicago truncatula nitrogen-fixing symbiosis, we defined free-living and symbiotic phenotypes of an hfq mutant. Over 500 transcripts were differentially accumulated in the hfq mutant of S. meliloti Rm1021 when grown in a shaking culture. Consistent with transcriptome-wide changes, the hfq mutant displayed dramatic alterations in metabolism of nitrogen-containing compounds, even though its carbon source utilization profiles were nearly identical to the wild type. The hfq mutant had reduced motility and was impaired for growth at alkaline pH. A deletion of hfq resulted in a reduced symbiotic efficiency, although the mutant was still able to initiate nodule development and differentiate into bacteroids.