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

Phage treatment of Pseudomonas aeruginosa yields a phage-resistant population with different susceptibility to innate immune responses and mild effects on metabolic profiles

In this study, we have investigated innate immune activation capacity and metabolic features of a population of P. aeruginosa PAO1 phage-resistant mutants with diverse genetic modification (large genomic deletions and point mutations) arising after exposure to phages targetting lipopolysaccharide (LPS) or Type-4 pili (T4P). Deletions led to the loss of genes involved in LPS synthesis, cell envelope permeability, efflux systems, biofilm production, oxidative stress tolerance, and DNA repair. Loss of LPS O antigen resulted in bacterial sensitivity to serum complement and stimulation of inflammatory cascades but did not cause increased phagocytosis, while T4P phage-resistant mutants were more effectively phagocytized than LPS-defective mutants. Changes in the utilization of different carbon, nitrogen, sulphur, and phosphorus sources were identified, especially in mutants where the two phage DNA persisted in the bacterial population (pseudolysogeny). However, the metabolic changes did not directly correlate with single-gene mutations or the large gene deletions, suggesting they reflect adaptive changes to the gene modifications that arise during the selection of resistant mutants. In contrast, phage-resistant mutants were susceptible to humoral innate immune responses, suggesting that phage resistance may be a beneficial outcome of phage therapy.

Proteogenomic Characterization of Pseudomonas veronii SM-20 Growing on Phenanthrene as Only Carbon and Energy Source

In this study, we conducted an extensive investigation of the biodegradation capabilities and stress response of the newly isolated strain Pseudomonas veronii SM-20 in order, to assess its potential for bioremediation of sites contaminated with polycyclic aromatic hydrocarbons (PAHs). Initially, phenotype microarray technology demonstrated the strain's proficiency in utilizing various carbon sources and its resistance to certain stressors. Genomic analysis has identified numerous genes involved in aromatic hydrocarbon metabolism. Biodegradation assay analyzed the depletion of phenanthrene (PHE) when it was added as a sole carbon and energy source. We found that P. veronii strain SM-20 degraded approximately 25% of PHE over a 30-day period, starting with an initial concentration of 600 µg/mL, while being utilized for growth. The degradation process involved PHE oxidation to an unstable arene oxide and 9,10-phenanthrenequinone, followed by ring-cleavage. Comparative proteomics provided a comprehensive understanding of how the entire proteome responded to PHE exposure, revealing the strain's adaptation in terms of aromatic metabolism, surface properties, and defense mechanism. In conclusion, our findings shed light on the promising attributes of P. veronii SM-20 and offer valuable insights for the use of P. veronii species in environmental restoration efforts targeting PAH-impacted sites.

Genomic, Molecular, and Phenotypic Characterization of Arthrobacter sp. OVS8, an Endophytic Bacterium Isolated from and Contributing to the Bioactive Compound Content of the Essential Oil of the Medicinal Plant Origanum vulgare L

Medicinal plants play an important role in the discovery of new bioactive compounds with antimicrobial activity, thanks to their pharmacological properties. However, members of their microbiota can also synthesize bioactive molecules. Among these, strains belonging to the genera Arthrobacter are commonly found associated with the plant's microenvironments, showing plant growth-promoting (PGP) activity and bioremediation properties. However, their role as antimicrobial secondary metabolite producers has not been fully explored. The aim of this work was to characterize the Arthrobacter sp. OVS8 endophytic strain, isolated from the medicinal plant Origanum vulgare L., from molecular and phenotypic viewpoints to evaluate its adaptation and influence on the plant internal microenvironments and its potential as a producer of antibacterial volatile molecules (VOCs). Results obtained from the phenotypic and genomic characterization highlight its ability to produce volatile antimicrobials effective against multidrug-resistant (MDR) human pathogens and its putative PGP role as a producer of siderophores and degrader of organic and inorganic pollutants. The outcomes presented in this work identify Arthrobacter sp. OVS8 as an excellent starting point toward the exploitation of bacterial endophytes as antibiotics sources.

Biogeographic patterns in populations of marine Pseudoalteromonas atlantica isolates

Intra-specific genomic diversity is well documented in microbes. The question, however, remains whether natural selection or neutral evolution is the major contributor to this diversity. We undertook this study to estimate genomic diversity in Pseudoalteromonas atlantica populations and whether the diversity, if present, could be attributed to environmental factors or distance effects. We isolated and sequenced twenty-three strains of P. atlantica from three geographically distant deep marine basins and performed comparative genomic analyses to study the genomic diversity of populations among these basins. Average nucleotide identity followed a strictly geographical pattern. In two out of three locations, the strains within the location exhibited >99.5% identity, whereas, among locations, the strains showed <98.11% identity. Phylogenetic and pan-genome analysis also reflected the biogeographical separation of the strains. Strains from the same location shared many accessory genes and clustered closely on the phylogenetic tree. Phenotypic diversity between populations was studied in ten out of twenty-three strains testing carbon and nitrogen source utilization and osmotolerance. A genetic basis for phenotypic diversity could be established in most cases but was apparently not influenced by local environmental conditions. Our study suggests that neutral evolution may have a substantial role in the biodiversity of P. atlantica.

High-quality genome-scale metabolic network reconstruction of probiotic bacterium Escherichia coli Nissle 1917

BACKGROUND

Escherichia coli Nissle 1917 (EcN) is a probiotic bacterium used to treat various gastrointestinal diseases. EcN is increasingly being used as a chassis for the engineering of advanced microbiome therapeutics. To aid in future engineering efforts, our aim was to construct an updated metabolic model of EcN with extended secondary metabolite representation.

RESULTS

An updated high-quality genome-scale metabolic model of EcN, iHM1533, was developed based on comparison with 55 E. coli/Shigella reference GEMs and manual curation, including expanded secondary metabolite pathways (enterobactin, salmochelins, aerobactin, yersiniabactin, and colibactin). The model was validated and improved using phenotype microarray data, resulting in an 82.3% accuracy in predicting growth phenotypes on various nutrition sources. Flux variability analysis with previously published ¹³C fluxomics data validated prediction of the internal central carbon fluxes. A standardised test suite called Memote assessed the quality of iHM1533 to have an overall score of 89%. The model was applied by using constraint-based flux analysis to predict targets for optimisation of secondary metabolite production. Modelling predicted design targets from across amino acid metabolism, carbon metabolism, and other subsystems that are common or unique for influencing the production of various secondary metabolites.

CONCLUSION

iHM1533 represents a well-annotated metabolic model of EcN with extended secondary metabolite representation. Phenotype characterisation and the iHM1533 model provide a better understanding of the metabolic capabilities of EcN and will help future metabolic engineering efforts.

Deciphering the global roles of Cold shock proteins in Listeria monocytogenes nutrient metabolism and stress tolerance

Listeria monocytogenes (Lm) accounts for serious public health and food safety problems owing to its stress resilience and pathogenicity. Based on their regulatory involvement in global gene expression events, cold-shock domain family proteins (Csps) are crucial in expression of various stress fitness and virulence phenotypes in bacteria. Lm possesses three Csps (CspA, CspB, and CspD) whose regulatory roles in the context of the genetic diversity of this bacterium are not yet fully understood. We examined the impacts of Csps deficiency on Lm nutrient metabolism and stress tolerance using a set of csp deletion mutants generated in different genetic backgrounds. Phenotype microarrays (PM) analysis showed that the absence of Csps in ∆cspABD reduces carbon (C-) source utilization capacity and increases Lm sensitivity to osmotic, pH, various chemical, and antimicrobial stress conditions. Single and double csp deletion mutants in different Lm genetic backgrounds were used to further dissect the roles of individual Csps in these phenotypes. Selected PM-based observations were further corroborated through targeted phenotypic assays, confirming that Csps are crucial in Lm for optimal utilization of various C-sources including rhamnose and glucose as well as tolerance against NaCl, β-phenyethylamine (PEA), and food relevant detergent stress conditions. Strain and genetic lineage background-based differences, division of labour, epistasis, and functional redundancies among the Csps were uncovered with respect to their roles in various processes including C-source utilization, cold, and PEA stress resistance. Finally, targeted transcriptome analysis was performed, revealing the activation of csp gene expression under defined stress conditions and the impact of Csps on expression regulation of selected rhamnose utilization genes. Overall, our study shows that Csps play important roles in nutrient utilization and stress responses in Lm strains, contributing to traits that are central to the public health and food safety impacts of this pathogen.

Genome-Scale Metabolic Modelling of Lifestyle Changes in Rhizobium leguminosarum

Biological nitrogen fixation in rhizobium-legume symbioses is of major importance for sustainable agricultural practices. To establish a mutualistic relationship with their plant host, rhizobia transition from free-living bacteria in soil to growth down infection threads inside plant roots and finally differentiate into nitrogen-fixing bacteroids. We reconstructed a genome-scale metabolic model for Rhizobium leguminosarum and integrated the model with transcriptome, proteome, metabolome, and gene essentiality data to investigate nutrient uptake and metabolic fluxes characteristic of these different lifestyles. Synthesis of leucine, polyphosphate, and AICAR is predicted to be important in the rhizosphere, while myo-inositol catabolism is active in undifferentiated nodule bacteria in agreement with experimental evidence. The model indicates that bacteroids utilize xylose and glycolate in addition to dicarboxylates, which could explain previously described gene expression patterns. Histidine is predicted to be actively synthesized in bacteroids, consistent with transcriptome and proteome data for several rhizobial species. These results provide the basis for targeted experimental investigation of metabolic processes specific to the different stages of the rhizobium-legume symbioses. IMPORTANCE Rhizobia are soil bacteria that induce nodule formation on plant roots and differentiate into nitrogen-fixing bacteroids. A detailed understanding of this complex symbiosis is essential for advancing ongoing efforts to engineer novel symbioses with cereal crops for sustainable agriculture. Here, we reconstruct and validate a genome-scale metabolic model for Rhizobium leguminosarum bv. viciae 3841. By integrating the model with various experimental data sets specific to different stages of symbiosis formation, we elucidate the metabolic characteristics of rhizosphere bacteria, undifferentiated bacteria inside root nodules, and nitrogen-fixing bacteroids. Our model predicts metabolic flux patterns for these three distinct lifestyles, thus providing a framework for the interpretation of genome-scale experimental data sets and identifying targets for future experimental studies.

Mesorhizobium comanense sp. nov., isolated from groundwater

Strain 3P27G6^T was isolated from an artesian well connected to the thermal water basin of Comano Terme, Province of Trento, Italy. In phylogenetic analyses based on multilocus sequence analysis, strain 3P27G6^T clustered together with Mesorhizobium australicum WSM2073^T. Genome sequencing produced a 99.51 % complete genome, with a length of 7 363 057 bp and G+C content of 63.53 mol%, containing 6897 coding sequences, 55 tRNA and three rRNA. Average nucleotide identity analysis revealed that all distances calculated between strain 3P27G6^T and the other Mesorhizobium genomes were below 0.9, indicating that strain 3P27G6^T represents a new species. Therefore, we propose the name Mesorhizobium comanense sp. nov. with the type strain 3P27G6^T (=DSM 110654^T=CECT 30067^T). Strain 3P27G6^T is a Gram-negative, rod-shaped, aerobic bacterium. Growth condition, antibiotic susceptibility, metabolic and fatty acid-methyl esters profiles of the strain were determined. Only few nodulation and nitrogen fixation genes were found in the genome, suggesting that this strain may not be specialized in nodulation or in nitrogen fixation.

Characterization and engineering of Streptomyces griseofuscus DSM 40191 as a potential host for heterologous expression of biosynthetic gene clusters

Streptomyces griseofuscus DSM 40191 is a fast growing Streptomyces strain that remains largely underexplored as a heterologous host. Here, we report the genome mining of S. griseofuscus, followed by the detailed exploration of its phenotype, including the production of native secondary metabolites and ability to utilise carbon, nitrogen, sulphur and phosphorus sources. Furthermore, several routes for genetic engineering of S. griseofuscus were explored, including use of GusA-based vectors, CRISPR-Cas9 and CRISPR-cBEST-mediated knockouts. Two out of the three native plasmids were cured using CRISPR-Cas9 technology, leading to the generation of strain S. griseofuscus DEL1. DEL1 was further modified by the full deletion of a pentamycin BGC and an unknown NRPS BGC, leading to the generation of strain DEL2, lacking approx. 500 kbp of the genome, which corresponds to a 5.19% genome reduction. DEL2 can be characterized by faster growth and inability to produce three main native metabolites: lankacidin, lankamycin, pentamycin and their derivatives. To test the ability of DEL2 to heterologously produce secondary metabolites, the actinorhodin BGC was used. We were able to observe a formation of a blue halo, indicating a potential production of actinorhodin by both DEL2 and a wild type.

Analysis of the Oxidative Stress Regulon Identifies soxS as a Genetic Target for Resistance Reversal in Multidrug-Resistant Klebsiella pneumoniae

In bacteria, the defense system deployed to counter oxidative stress is orchestrated by three transcriptional factors, SoxS, SoxR, and OxyR. Although the regulon that these factors control is known in many bacteria, similar data are not available for Klebsiella pneumoniae. To address this data gap, oxidative stress was artificially induced in K. pneumoniae MGH78578 using paraquat and the corresponding oxidative stress regulon recorded using transcriptome sequencing (RNA-seq). The soxS gene was significantly induced during oxidative stress, and a knockout mutant was constructed to explore its functionality. The wild type and mutant were grown in the presence of paraquat and subjected to RNA-seq to elucidate the soxS regulon in K. pneumoniae MGH78578. Genes that are commonly regulated both in the oxidative stress and soxS regulons were identified and denoted as the oxidative SoxS regulon; these included a group of genes specifically regulated by SoxS. Efflux pump-encoding genes and global regulators were identified as part of this regulon. Consequently, the isogenic soxS mutant was found to exhibit a reduction in the minimum bactericidal concentration against tetracycline compared to that of the wild type. Impaired efflux activity, allowing tetracycline to be accumulated in the cytoplasm to bactericidal levels, was further evaluated using a tetraphenylphosphonium (TPP⁺) accumulation assay. The soxS mutant was also susceptible to tetracycline in vivo in a zebrafish embryo model. We conclude that the soxS gene could be considered a genetic target against which an inhibitor could be developed and used in combinatorial therapy to combat infections associated with multidrug-resistant K. pneumoniae.

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.

Microbial high throughput phenomics: The potential of an irreplaceable omics

The phenotype-genotype landscape is a projection coming from detailed phenotypic and genotypic data under environmental pressure. Although phenome of microbes or microbial consortia mirrors the functional expression of a genome or set of genomes, metabolic traits rely on the phenotype. Phenomics has the potential to revolution functional genomics. In this review, we discuss why and how phenomics was developed. We described how phenomics may extend our understanding of the assembly of microbial consortia and their functionality, and then we outlined the novel applications within the study of phenomes using Omnilog platform together with a revision of its current application to study lactic acid bacteria (LAB) metabolic traits during food processing. LAB were proposed as a suitable model system to analyze and discuss the implementation and exploitation of this emerging omics approach. We introduced the 'phenotype switching', as a new phenotype microarray approach to get insights in bacterial physiology. An overview of methodologies and tools to manage and analyze the generated data was provided. Finally, pro and cons of pipelines developed so far, including the most innovative ones were critically analyzed. We propose an R pipeline, recently deposited, which allows to automatically analyze Omnilog data integrating the latest approaches and implementing the new concepts described here.

Evolution of Listeria monocytogenes During a Persistent Human Prosthetic Hip Joint Infection

Listeria monocytogenes associated prosthetic joint infections (PJI) are a rare but increasing clinical problem of listeriosis. We characterized two isolates of the same L. monocytogenes strain isolated within five years of each other from a recurrent human prosthetic joint infection. The two isolates although clonally identical were phenotypically distinct confirming that the original infection strain had evolved within the human host PJI environment giving rise to a phenotypically distinct variant. The recurrent PJI isolate displayed various phenotypic differences compared to the parental original PJI isolate including diminished growth and carbon source metabolism, as well as altered morphology and increased stress sensitivity. The PJI isolates were both diminished in virulence due to an identical truncation mutation in the major virulence regulator PrfA. Genome wide sequence comparison provided conclusive evidence that the two isolates were identical clonal descendants of the same L. monocytogenes strain that had evolved through acquisition of various single nucleotide polymorphisms (SNPs) as well as insertion and deletion events (InDels) during a persistent human PJI. Acquired genetic changes included a specific mutation causing premature stop codon (PMSC) and truncation of RNAse J1 protein. Based on analysis of this naturally truncated as well as other complete RNAse J1 deletion mutants we show that the long-term survival of this specific L. monocytogenes strain within the prosthetic joint might in part be explained by the rnjA PMSC mutation that diminishes virulence and activation of the host immune system in a zebrafish embryo localized infection model. Overall our analysis of this special natural case provides insights into random mutation events and molecular mechanisms that might be associated with the adaptation and short-term evolution of this specific L. monocytogenes strain within a persistent human PJI environment.

The metabolic shift in highly and weakly virulent Dickeya solani strains is more affected by temperature than by mutations in genes encoding global virulence regulators

Global warming may shortly increase the risk of disease development on plants. Significant differences in the metabolic activity screened with Phenotype Microarray at 22°C and 28°C were observed between D. solani strains with high and low virulence level. Highly virulent D. solani was characterized by a higher number of metabolized compounds and a faster metabolism and was more tolerant to non-favorable pH and osmolarity. Metabolic phenotyping showed for the first time that the mutation in pecT gene, which encodes a global repressor of virulence, affects several pathways of the basic cell metabolism. PecT mutants had a higher maceration capacity of potato tissue and showed a higher pectinolytic activity than the wild-type strains. On the contrary, mutation in expI gene, which encoded the signaling molecules synthase crucial for quorum sensing, had an insignificant effect on the cell metabolism, although it slightly reduced the potato tissue maceration. The ability to utilize most of the tested compounds was higher at 28°C, while the survival at non-favorable pH and osmolarity was higher at 22°C. These results proved that the temperature of incubation had the most significant impact on the D. solani metabolic profiles.

Genomic and phenomic analysis of a marine bacterium, Photobacterium marinum J15

Photobacterium species are widely distributed in the marine environment. The overall metabolism of this genus remains largely unknown. In order to improve our knowledge on this bacterium, the relationship between the genome and phenome of the Photobacterium isolate was analyzed. The cream colored, Gram-negative, rod-shaped and motile bacterial strain, J15, was isolated from marine water of Tanjung Pelepas, Johor, Malaysia. The 5,684,538 bp genome of strain J15 comprised 3 contigs (2 chromosomes and 1 plasmid) with G + C content of 46.39 % and contained 4924 protein-coding genes including 180 tRNAs and 40 rRNAs. The phenotypic microarray (PM) as analyzed using BIOLOG showed the utilization of; i) 93 of the 190 carbon sources tested, where 61 compounds were used efficiently; ii) 41 of the 95 nitrogen sources tested, where 22 compounds were used efficiently; and iii) 3 of the 94 phosphorous and sulphur sources tested. Furthermore, high tolerance to osmotic stress, basic pH and toxic compounds as well as resistance to antibiotics of strain J15 were determined by BIOLOG PM. The ANI and kSNP analyses revealed that strain J15 to be the same species with Photobacterium marinum AK15 with ANI value of 96.93 % and bootstrapping value of 100 in kSNP. Based on the ANI and kSNP analyses, strain J15 was identified as P. marinum J15.

Variable Carbon Source Utilization, Stress Resistance, and Virulence Profiles Among Listeria monocytogenes Strains Responsible for Listeriosis Outbreaks in Switzerland

A combination of phenotype microarrays, targeted stress resistance and virulence assays and comparative genome analysis was used to compare a set of Listeria monocytogenes strains including those involved in previous Swiss foodborne listeriosis outbreaks. Despite being highly syntenic in gene content these strains showed significant phenotypic variation in utilization of different carbon (C)-sources as well as in resistance of osmotic and pH stress conditions that are relevant to host and food associated environments. An outbreak strain from the 2005 Swiss Tomme cheese listeriosis outbreak (Lm3163) showed the highest versatility in C-sources utilized whereas the strain responsible for the 1983 to 1987 Vacherin Montd'or cheese listeriosis outbreak (LL195) showed the highest tolerance to both osmotic and pH stress conditions among the examined strains. Inclusion of L-norvaline led to enhanced resistance of acidic stress in all the examined strains and there were strain-strain-specific differences observed in the ability of other amino acids and urea to enhance acid stress resistance in L. monocytogenes. A strain dependent inhibition pattern was also observed upon inclusion of β-phenylethylamine under alkaline stress conditions. In targeted phenotypic analysis the strain-specific differences in salt stress tolerance uncovered in phenotypic microarrays were corroborated and variations in host cell invasion and virulence among the examined strains were also revealed. Outbreak associated strains representing lineage I serotype 4b showed superior pathogenicity in a zebrafish infection model whilst Lm3163 a lineage II serotype 1/2a outbreak strain demonstrated the highest cellular invasion capacity amongst the tested strains. A genome wide sequence comparison of the strains only revealed few genetic differences between the strains suggesting that variations in gene regulation and expression are largely responsible for the phenotypic differences revealed among the examined strains. Our results have generated data that provides a potential basis for the future design of improved Listeria specific media to enhance routine detection and isolation of this pathogen as well as provide knowledge for developing novel methods for its control in food.

Metabolic Modeling of Pectobacterium parmentieri SCC3193 Provides Insights into Metabolic Pathways of Plant Pathogenic Bacteria

Understanding plant⁻microbe interactions is crucial for improving plants' productivity and protection. Constraint-based metabolic modeling is one of the possible ways to investigate the bacterial adaptation to different ecological niches and may give insights into the metabolic versatility of plant pathogenic bacteria. We reconstructed a raw metabolic model of the emerging plant pathogenic bacterium Pectobacterium parmentieri SCC3193 with the use of KBase. The model was curated by using inParanoind and phenotypic data generated with the use of the OmniLog system. Metabolic modeling was performed through COBRApy Toolbox v. 0.10.1. The curated metabolic model of P. parmentieri SCC3193 is highly reliable, as in silico obtained results overlapped up to 91% with experimental data on carbon utilization phenotypes. By mean of flux balance analysis (FBA), we predicted the metabolic adaptation of P. parmentieri SCC3193 to two different ecological niches, relevant for the persistence and plant colonization by this bacterium: soil and the rhizosphere. We performed in silico gene deletions to predict the set of essential core genes for this bacterium to grow in such environments. We anticipate that our metabolic model will be a valuable element for defining a set of metabolic targets to control infection and spreading of this plant pathogen.

O-Antigen-Dependent Colicin Insensitivity of Uropathogenic Escherichia coli

The outer membrane of Gram-negative bacteria presents a significant barrier for molecules entering the cell. Nevertheless, colicins, which are antimicrobial proteins secreted by Escherichia coli, can target other E. coli cells by binding to cell surface receptor proteins and activating their import, resulting in cell death. Previous studies have documented high rates of nonspecific resistance (insensitivity) of various E. coli strains toward colicins that is independent of colicin-specific immunity and is instead associated with lipopolysaccharide (LPS) in the outer membrane. This observation poses a contradiction: why do E. coli strains have colicin-expressing plasmids, which are energetically costly to retain, if cells around them are likely to be naturally insensitive to the colicin they produce? Here, using a combination of transposon sequencing and phenotypic microarrays, we show that colicin insensitivity of uropathogenic E. coli sequence type 131 (ST131) is dependent on the production of its O-antigen but that minor changes in growth conditions render the organism sensitive toward colicins. The reintroduction of O-antigen into E. coli K-12 demonstrated that it is the density of O-antigen that is the dominant factor governing colicin insensitivity. We also show, by microscopy of fluorescently labelled colicins, that growth conditions affect the degree of occlusion by O-antigen of outer membrane receptors but not the clustered organization of receptors. The result of our study demonstrate that environmental conditions play a critical role in sensitizing E. coli toward colicins and that O-antigen in LPS is central to this role.IMPORTANCEEscherichia coli infections can be a major health burden, especially with the organism becoming increasingly resistant to "last-resort" antibiotics such as carbapenems. Although colicins are potent narrow-spectrum antimicrobials with potential as future antibiotics, high levels of naturally occurring colicin insensitivity have been documented which could limit their efficacy. We identify O-antigen-dependent colicin insensitivity in a clinically relevant uropathogenic E. coli strain and show that this insensitivity can be circumvented by minor changes to growth conditions. The results of our study suggest that colicin insensitivity among E. coli organisms has been greatly overestimated, and as a consequence, colicins could in fact be effective species-specific antimicrobials targeting pathogenic E. coli such as uropathogenic E. coli (UPEC).

Exposure to Sub-inhibitory Concentrations of the Chemosensitizer 1-(1-Naphthylmethyl)-Piperazine Creates Membrane Destabilization in Multi-Drug Resistant Klebsiella pneumoniae

Antimicrobial efflux is one of the important mechanisms causing multi-drug resistance (MDR) in bacteria. Chemosensitizers like 1-(1-naphthylmethyl)-piperazine (NMP) can inhibit an efflux pump and therefore can overcome MDR. However, secondary effects of NMP other than efflux pump inhibition are rarely investigated. Here, using phenotypic assays, phenotypic microarray and transcriptomic assays we show that NMP creates membrane destabilization in MDR Klebsiella pneumoniae MGH 78578 strain. The NMP mediated membrane destabilization activity was measured using β-lactamase activity, membrane potential alteration studies, and transmission electron microscopy assays. Results from both β-lactamase and membrane potential alteration studies shows that both outer and inner membranes are destabilized in NMP exposed K. pneumoniae MGH 78578 cells. Phenotypic Microarray and RNA-seq were further used to elucidate the metabolic and transcriptional signals underpinning membrane destabilization. Membrane destabilization happens as early as 15 min post-NMP treatment. Our RNA-seq data shows that many genes involved in envelope stress response were differentially regulated in the NMP treated cells. Up-regulation of genes encoding the envelope stress response and repair systems show the distortion in membrane homeostasis during survival in an environment containing sub-inhibitory concentration of NMP. In addition, the lsr operon encoding the production of autoinducer-2 responsible for biofilm production was down-regulated resulting in reduced biofilm formation in NMP treated cells, a phenotype confirmed by crystal violet-based assays. We postulate that the early membrane disruption leads to destabilization of inner membrane potential, impairing ATP production and consequently resulting in efflux pump inhibition.

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.

Galleria mellonella as an insect model for P. destructans, the cause of White-nose Syndrome in bats

Pseudogymnoascus destructans is the fungal pathogen responsible for White-nose Syndrome (WNS), a disease that has killed millions of bats in North America over the last decade. A major obstacle to research on P. destructans has been the lack of a tractable infection model for monitoring virulence. Here, we establish a high-throughput model of infection using larvae of Galleria mellonella, an invertebrate used to study host-pathogen interactions for a wide range of microbial species. We demonstrate that P. destructans can kill G. mellonella larvae in an inoculum-dependent manner when infected larvae are housed at 13°C or 18°C. Larval killing is an active process, as heat-killed P. destructans spores caused significantly decreased levels of larval death compared to live spores. We also show that fungal spores that were germinated prior to inoculation were able to kill larvae 3–4 times faster than non-germinated spores. Lastly, we identified chemical inhibitors of P. destructans and used G. mellonella to evaluate these inhibitors for their ability to reduce virulence. We demonstrate that amphotericin B can effectively block larval killing by P. destructans and thereby establish that this infection model can be used to screen biocontrol agents against this fungal pathogen.

Paracetamol - toxicity and microbial utilization. Pseudomonas moorei KB4 as a case study for exploring degradation pathway

Paracetamol, a widely used analgesic and antipyretic drug, is currently one of the most emerging pollutants worldwide. Besides its wide prevalence in the literature only several bacterial strains able to degrade this compound have been described. In this study, we isolated six new bacterial strains able to remove paracetamol. The isolated strains were identified as the members of Pseudomonas, Bacillus, Acinetobacter and Sphingomonas genera and characterized phenotypically and biochemically using standard methods. From the isolated strains, Pseudomonas moorei KB4 was able to utilize 50 mg L^-1 of paracetamol. As the main degradation products, p-aminophenol and hydroquinone were identified. Based on the measurements of specific activity of acyl amidohydrolase, deaminase and hydroquinone 1,2-dioxygenase and the results of liquid chromatography analyses, we proposed a mechanism of paracetamol degradation by KB4 strain under co-metabolic conditions with glucose. Additionally, toxicity bioassays and the influence of various environmental factors, including pH, temperature, heavy metals at no-observed-effective-concentrations, and the presence of aromatic compounds on the efficiency and mechanism of paracetamol degradation by KB4 strain were determined. This comprehensive study about paracetamol biodegradation will be helpful in designing a treatment systems of wastewaters contaminated with paracetamol.

Growth Score: a single metric to define growth in 96-well phenotype assays

High-throughput phenotype assays are a cornerstone of systems biology as they allow direct measurements of mutations, genes, strains, or even different genera. High-throughput methods also require data analytic methods that reduce complex time-series data to a single numeric evaluation. Here, we present the Growth Score, an improvement on the previous Growth Level formula. There is strong correlation between Growth Score and Growth Level, but the new Growth Score contains only essential growth curve properties while the formula of the previous Growth Level was convoluted and not easily interpretable. Several programs can be used to estimate the parameters required to calculate the Growth Score metric, including our PMAnalyzer pipeline.

Integrated phenotypic-genotypic approach to understand the influence of ultrasound on metabolic response of Lactobacillus sakei

The lethal effects of soundwaves on a range of microorganisms have been known for almost a century whereas, the use of ultrasound to promote or control their activity is much more recent. Moreover, the fundamental molecular mechanism influencing the behaviour of microorganisms subjected to ultrasonic waves is not well established. In this study, we investigated the influence of ultrasonic frequencies of 20, 45, 130 and 950 kHz on growth kinetics of Lactobacillus sakei. A significant increase in the growth rate of L. sakei was observed following ultrasound treatment at 20 kHz despite the treatment yielding a significant reduction of ca. 3 log cfu/mL in cells count. Scanning electron microscopy showed that ultrasound caused significant changes on the cell surface of L. sakei culture with the formation of pores "sonoporation". Phenotypic microarrays showed that all ultrasound treated L. sakei after exposure to various carbon, nitrogen, phosphorus and sulphur sources had significant variations in nutrient utilisation. Integration of this phenotypic data with the genome of L. sakei revealed that various metabolic pathways were being influenced by the ultrasound treatments. Results presented in this study showed that the physiological response of L. sakei in response to US is frequency dependent and that it can influence metabolic pathways. Hence, ultrasound treatments can be employed to modulate microbial activity for specialised applications.

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.

Ecology of cold environments: new insights of bacterial metabolic adaptation through an integrated genomic-phenomic approach

Cold environments dominate Earth's biosphere, hosting complex microbial communities with the ability to thrive at low temperatures. However, the underlying molecular mechanisms and the metabolic pathways involved in bacterial cold-adaptation mechanisms are still not fully understood. Herein, we assessed the metabolic features of the Antarctic bacterium Pseudoalteromonas haloplanktis TAC125 (PhTAC125), a model organism for cold-adaptation, at both 4 °C and 15 °C, by integrating genomic and phenomic (high-throughput phenotyping) data and comparing the obtained results to the taxonomically related Antarctic bacterium Pseudoalteromonas sp. TB41 (PspTB41). Although the genome size of PspTB41 is considerably larger than PhTAC125, the higher number of genes did not reflect any higher metabolic versatility at 4 °C as compared to PhTAC125. Remarkably, protein S-thiolation regulated by glutathione and glutathionylspermidine appeared to be a new possible mechanism for cold adaptation in PhTAC125. More in general, this study represents an example of how 'multi-omic' information might potentially contribute in filling the gap between genotypic and phenotypic features related to cold-adaptation mechanisms in bacteria.

Two stable variants of Burkholderia pseudomallei strain MSHR5848 express broadly divergent in vitro phenotypes associated with their virulence differences

Burkholderia pseudomallei (Bp), the agent of melioidosis, causes disease ranging from acute and rapidly fatal to protracted and chronic. Bp is highly infectious by aerosol, can cause severe disease with nonspecific symptoms, and is naturally resistant to multiple antibiotics. However, no vaccine exists. Unlike many Bp strains, which exhibit random variability in traits such as colony morphology, Bp strain MSHR5848 exhibited two distinct and relatively stable colony morphologies on sheep blood agar plates: a smooth, glossy, pale yellow colony and a flat, rough, white colony. Passage of the two variants, designated "Smooth" and "Rough", under standard laboratory conditions produced cultures composed of > 99.9% of the single corresponding type; however, both could switch to the other type at different frequencies when incubated in certain nutritionally stringent or stressful growth conditions. These MSHR5848 derivatives were extensively characterized to identify variant-associated differences. Microscopic and colony morphology differences on six differential media were observed and only the Rough variant metabolized sugars in selective agar. Antimicrobial susceptibilities and lipopolysaccharide (LPS) features were characterized and phenotype microarray profiles revealed distinct metabolic and susceptibility disparities between the variants. Results using the phenotype microarray system narrowed the 1,920 substrates to a subset which differentiated the two variants. Smooth grew more rapidly in vitro than Rough, yet the latter exhibited a nearly 10-fold lower lethal dose for mice than Smooth. Finally, the Smooth variant was phagocytosed and replicated to a greater extent and was more cytotoxic than Rough in macrophages. In contrast, multiple locus sequence type (MLST) analysis, ribotyping, and whole genome sequence analysis demonstrated the variants' genetic conservation; only a single consistent genetic difference between the two was identified for further study. These distinct differences shown by two variants of a Bp strain will be leveraged to better understand the mechanism of Bp phenotypic variability and to possibly identify in vitro markers of infection.

The pangenome of (Antarctic) Pseudoalteromonas bacteria: evolutionary and functional insights

BACKGROUND

Pseudoalteromonas is a genus of ubiquitous marine bacteria used as model organisms to study the biological mechanisms involved in the adaptation to cold conditions. A remarkable feature shared by these bacteria is their ability to produce secondary metabolites with a strong antimicrobial and antitumor activity. Despite their biotechnological relevance, representatives of this genus are still lacking (with few exceptions) an extensive genomic characterization, including features involved in the evolution of secondary metabolites production. Indeed, biotechnological applications would greatly benefit from such analysis.

RESULTS

Here, we analyzed the genomes of 38 strains belonging to different Pseudoalteromonas species and isolated from diverse ecological niches, including extreme ones (i.e. Antarctica). These sequences were used to reconstruct the largest Pseudoalteromonas pangenome computed so far, including also the two main groups of Pseudoalteromonas strains (pigmented and not pigmented strains). The downstream analyses were conducted to describe the genomic diversity, both at genus and group levels. This allowed highlighting a remarkable genomic heterogeneity, even for closely related strains. We drafted all the main evolutionary steps that led to the current structure and gene content of Pseudoalteromonas representatives. These, most likely, included an extensive genome reduction and a strong contribution of Horizontal Gene Transfer (HGT), which affected biotechnologically relevant gene sets and occurred in a strain-specific fashion. Furthermore, this study also identified the genomic determinants related to some of the most interesting features of the Pseudoalteromonas representatives, such as the production of secondary metabolites, the adaptation to cold temperatures and the resistance to abiotic compounds.

CONCLUSIONS

This study poses the bases for a comprehensive understanding of the evolutionary trajectories followed in time by this peculiar bacterial genus and for a focused exploitation of their biotechnological potential.

Phenotypic Profiling Reveals that Candida albicans Opaque Cells Represent a Metabolically Specialized Cell State Compared to Default White Cells

The white-opaque switch is a bistable, epigenetic transition affecting multiple traits in Candida albicans including mating, immunogenicity, and niche specificity. To compare how the two cell states respond to external cues, we examined the fitness, phenotypic switching, and filamentation properties of white cells and opaque cells under 1,440 different conditions at 25°C and 37°C. We demonstrate that white and opaque cells display striking differences in their integration of metabolic and thermal cues, so that the two states exhibit optimal fitness under distinct conditions. White cells were fitter than opaque cells under a wide range of environmental conditions, including growth at various pHs and in the presence of chemical stresses or antifungal drugs. This difference was exacerbated at 37°C, consistent with white cells being the default state of C. albicans in the mammalian host. In contrast, opaque cells showed greater fitness than white cells under select nutritional conditions, including growth on diverse peptides at 25°C. We further demonstrate that filamentation is significantly rewired between the two states, with white and opaque cells undergoing filamentous growth in response to distinct external cues. Genetic analysis was used to identify signaling pathways impacting the white-opaque transition both in vitro and in a murine model of commensal colonization, and three sugar sensing pathways are revealed as regulators of the switch. Together, these findings establish that white and opaque cells are programmed for differential integration of metabolic and thermal cues and that opaque cells represent a more metabolically specialized cell state than the default white state.

IMPORTANCE

Epigenetic transitions are an important mechanism by which microbes adapt to external stimuli. For Candida albicans, such transitions are crucial for adaptation to complex, fluctuating environments, and therefore contribute to its success as a human pathogen. The white-opaque switch modulates multiple C. albicans attributes, from sexual competency to niche specificity. Here, we demonstrate that metabolic circuits are extensively rewired between white and opaque states, so that the two cell types exhibit optimal fitness under different nutritional conditions and at different temperatures. We thereby establish that epigenetic events can profoundly alter the metabolism of fungal cells. We also demonstrate that epigenetic switching regulates filamentation and biofilm formation, two phenotypes closely associated with pathogenesis. These experiments reveal that white cells, considered the most clinically relevant form of C. albicans, are a "general-purpose" state suited to many environments, whereas opaque cells appear to represent a more metabolically specialized form of the species.

Phenotypic and genomic characterization of the antimicrobial producer Rheinheimera sp. EpRS3 isolated from the medicinal plant Echinacea purpurea: insights into its biotechnological relevance

In recent years, there has been increasing interest in plant microbiota; however, despite medicinal plant relevance, very little is known about their highly complex endophytic communities. In this work, we report on the genomic and phenotypic characterization of the antimicrobial compound producer Rheinheimera sp. EpRS3, a bacterial strain isolated from the rhizospheric soil of the medicinal plant Echinacea purpurea. In particular, EpRS3 is able to inhibit growth of different bacterial pathogens (Bcc, Acinetobacter baumannii, and Klebsiella pneumoniae) which might be related to the presence of gene clusters involved in the biosynthesis of different types of secondary metabolites. The outcomes presented in this work highlight the fact that the strain possesses huge biotechnological potential; indeed, it also shows antimicrobial effects upon well-described multidrug-resistant (MDR) human pathogens, and it affects plant root elongation and morphology, mimicking indole acetic acid (IAA) action.

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.

The in vitro and in vivo effects of constitutive light expression on a bioluminescent strain of the mouse enteropathogen Citrobacter rodentium

Bioluminescent reporter genes, such as those from fireflies and bacteria, let researchers use light production as a non-invasive and non-destructive surrogate measure of microbial numbers in a wide variety of environments. As bioluminescence needs microbial metabolites, tagging microorganisms with luciferases means only live metabolically active cells are detected. Despite the wide use of bioluminescent reporter genes, very little is known about the impact of continuous (also called constitutive) light expression on tagged bacteria. We have previously made a bioluminescent strain of Citrobacter rodentium, a bacterium which infects laboratory mice in a similar way to how enteropathogenic Escherichia coli (EPEC) and enterohaemorrhagic E. coli (EHEC) infect humans. In this study, we compared the growth of the bioluminescent C. rodentium strain ICC180 with its non-bioluminescent parent (strain ICC169) in a wide variety of environments. To understand more about the metabolic burden of expressing light, we also compared the growth profiles of the two strains under approximately 2,000 different conditions. We found that constitutive light expression in ICC180 was near-neutral in almost every non-toxic environment tested. However, we also found that the non-bioluminescent parent strain has a competitive advantage over ICC180 during infection of adult mice, although this was not enough for ICC180 to be completely outcompeted. In conclusion, our data suggest that constitutive light expression is not metabolically costly to C. rodentium and supports the view that bioluminescent versions of microbes can be used as a substitute for their non-bioluminescent parents to study bacterial behaviour in a wide variety of environments.

Comparative genome-scale modelling of Staphylococcus aureus strains identifies strain-specific metabolic capabilities linked to pathogenicity

Staphylococcus aureus is a preeminent bacterial pathogen capable of colonizing diverse ecological niches within its human host. We describe here the pangenome of S. aureus based on analysis of genome sequences from 64 strains of S. aureus spanning a range of ecological niches, host types, and antibiotic resistance profiles. Based on this set, S. aureus is expected to have an open pangenome composed of 7,411 genes and a core genome composed of 1,441 genes. Metabolism was highly conserved in this core genome; however, differences were identified in amino acid and nucleotide biosynthesis pathways between the strains. Genome-scale models (GEMs) of metabolism were constructed for the 64 strains of S. aureus These GEMs enabled a systems approach to characterizing the core metabolic and panmetabolic capabilities of the S. aureus species. All models were predicted to be auxotrophic for the vitamins niacin (vitamin B3) and thiamin (vitamin B1), whereas strain-specific auxotrophies were predicted for riboflavin (vitamin B2), guanosine, leucine, methionine, and cysteine, among others. GEMs were used to systematically analyze growth capabilities in more than 300 different growth-supporting environments. The results identified metabolic capabilities linked to pathogenic traits and virulence acquisitions. Such traits can be used to differentiate strains responsible for mild vs. severe infections and preference for hosts (e.g., animals vs. humans). Genome-scale analysis of multiple strains of a species can thus be used to identify metabolic determinants of virulence and increase our understanding of why certain strains of this deadly pathogen have spread rapidly throughout the world.

Phenotype MicroArray™ system in the study of fungal functional diversity and catabolic versatility

Fungi cover a range of important ecological functions associated with nutrient and carbon cycling in leaf litter and soil. As a result, research on existing relationships between fungal functional diversity, decomposition rates and competition is of key interest. Indeed, availability of nutrients in soil is largely the consequence of organic matter degradation dynamics. The Biolog^® Phenotype MicroArrays™ (PM) system allows for the testing of fungi against many different carbon sources at any one time. The use and potential of the PM system as a tool for studying niche overlap and catabolic versatility of saprotrophic fungi is discussed here, and examples of its application are provided.

Draft genome sequence and overview of the purple non sulfur bacterium Rhodopseudomonas palustris 42OL

Rhodopseudomonas palustris strain 42OL was isolated in 1973 from a sugar refinery waste treatment pond. The strain has been prevalently used for hydrogen production processes using a wide variety of waste-derived substrates, and cultured both indoors and outdoors, either freely suspended or immobilized. R. palustris 42OL was suitable for many other applications and capable of growing in very different culturing conditions, revealing a wide metabolic versatility. The analysis of the genome sequence allowed to identify the metabolic pathways for hydrogen and poly-β-hydroxy-butyrate production, and confirmed the ability of using a wide range of organic acids as substrates.

Genomic and phenotypic characterization of the species Acinetobacter venetianus

Crude oil is a complex mixture of hydrocarbons and other organic compounds that can produce serious environmental problems and whose removal is highly demanding in terms of human and technological resources. The potential use of microbes as bioremediation agents is one of the most promising fields in this area. Members of the species Acinetobacter venetianus have been previously characterized for their capability to degrade n-alkanes and thus may represent interesting model systems to implement this process. Although a preliminary experimental characterization of the overall hydrocarbon degradation capability has been performed for five of them, to date, the genetic/genomic features underlying such molecular processes have not been identified. Here we have integrated genomic and phenotypic information for six A. venetianus strains, i.e. VE-C3, RAG-1^T, LUH 13518, LUH 7437, LUH 5627 and LUH 8758. Besides providing a thorough description of the A. venetianus species, these data were exploited to infer the genetic features (presence/absence patterns of genes) and the short-term evolutionary events possibly responsible for the variability in n-alkane degradation efficiency of these strains, including the mechanisms of interaction with the fuel droplet and the subsequent catabolism of this pollutant.

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.

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.

The influence of transcriptional regulatory networks on the evolution of bacterial pangenomes has not yet been elucidated, even though the role of transcriptional regulation is widely recognized. Using the model symbiont Sinorhizobium meliloti we have predicted the regulatory targets of 41 transcription factors in 51 strains and 5 other rhizobial species, showing a correlation between regulon diversity and pangenome evolution, through upstream sequence diversity and accessory genome composition. We have also shown that genes not wired to the regulatory network are more likely to belong to the accessory genome, thus suggesting that inclusion in the regulatory circuits may be an indicator of gene conservation. We have also highlighted a series of transcription factors that preferentially regulate genes belonging to one of the three replicons of this species, indicating the presence of replicon-specific regulatory modules, with peculiar functional signatures. At the same time the chromid shares a significant part of the regulatory network with the chromosome, indicating an additional way by which this replicon integrates itself in the pangenome.

Defining orthologs and pangenome size metrics

Since the advent of ultra-massive sequencing techniques, the consequent drop-off in both price and time required made feasible the sequencing of increasingly more genomes from microbes belonging to the same taxonomic unit. Eventually, this led to the concept of pangenome, that is, the entire set of genes present in a group of representatives of the same genus/species, which, in turn, can be divided into core genome, defined as the set of those genes present in all the genomes under study, and a dispensable genome, the set of genes possessed only by one or a subset of organism. When analyzing a pangenome, an interesting point is to measure its size, thus estimating the gene repertoire of a given taxonomic group. This is usually performed counting the novel genes added to the overall pangenome when new genomes are sequenced and annotated. A pangenome can be also classified as open or close: in an open pangenome its size increases indefinitely when adding new genomes; thus sequencing additional strains will likely yield novel genes. Conversely, in a close pangenome, adding new genomes will not lead to the discovery of new coding capabilities. A central point in pangenomics is the definition of homology relationships between genes belonging to different genomes. This may turn into the search of those genes with similar sequences between different organisms (and including both paralogous and orthologous genes). In this chapter, methods for finding groups of orthologs between genomes and for estimating the pangenome size are discussed. Also, working codes to address these tasks are provided.

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.

Microbial bioinformatics for food safety and production

In the production of fermented foods, microbes play an important role. Optimization of fermentation processes or starter culture production traditionally was a trial-and-error approach inspired by expert knowledge of the fermentation process. Current developments in high-throughput 'omics' technologies allow developing more rational approaches to improve fermentation processes both from the food functionality as well as from the food safety perspective. Here, the authors thematically review typical bioinformatics techniques and approaches to improve various aspects of the microbial production of fermented food products and food safety.

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.

High-throughput phenomics

Standard protocols are available in order to apply Phenotype MicroArray (PM) technology to characterize different groups of microorganisms. Nevertheless, there is the need to pay attention to several crucial steps in order to obtain high-quality and reproducible data from PM, such as the choice of the Dye mix, the type and concentration of the carbon source in metabolic experiments, the use of a buffered medium. A systematic research of auxotrophies in strains to be tested should be carefully evaluated before starting with PM experiments. Detailed protocols to obtain defined and reproducible phenotypic profiles for bacteria and yeasts are shown. Moreover, the innovative software opm R packages and DuctApe suite for the analysis of kinetic data produced by PM and panphenome description are reported.