Comparison of 432 Pseudomonas strains through integration of genomic, functional, metabolic and expression data

An ecological-evolutionary perspective on the genomic diversity and habitat preferences of the Acidobacteriota

Members of the phylum Acidobacteriota inhabit a wide range of ecosystems including soils. We analysed the global patterns of distribution and habitat preferences of various Acidobacteriota lineages across major ecosystems (soil, engineered, host-associated, marine, non-marine saline and alkaline and terrestrial non-soil ecosystems) in 248 559 publicly available metagenomic datasets. Classes Terriglobia, Vicinamibacteria, Blastocatellia and Thermoanaerobaculia were highly ubiquitous and showed a clear preference to soil over non-soil habitats, while classes Aminicenantia and Holophagae showed preferences to non-soil habitats. However, while specific preferences were observed, most Acidobacteriota lineages were habitat generalists rather than specialists, with genomic and/or metagenomic fragments recovered from soil and non-soil habitats at various levels of taxonomic resolution. Comparative analysis of 1930 genomes strongly indicates that phylogenetic affiliation plays a more important role than the habitat from which the genome was recovered in shaping the genomic characteristics and metabolic capacities of the Acidobacteriota. The observed lack of strong habitat specialization and habitat-transition-driven lineage evolution in the Acidobacteriota suggest ready cross-colonization between soil and non-soil habitats. We posit that such capacity is key to the successful establishment of Acidobacteriota as a major component in soil microbiomes post-ecosystem disturbance events or during pedogenesis.

Shotgun proteomic analyses of Pseudomonas species isolated from fish products

Numerous Pseudomonas species can infect aquatic animals, such as farmed rainbow trout, sea trout, sea bass, and sea bream, by causing disease or stress reactions. In aquaculture facilities, a number of Pseudomonas species have been isolated and identified as the main pathogens. The present study describes the characterization of 18 Pseudomonas strains, isolated from fish products using shotgun proteomics. The bacterial proteomes obtained were further analyzed to identify the main functional pathway proteins involved. In addition, this study revealed the presence of 1015 non-redundant peptides related to virulence factors. An additional 25 species-specific peptides were identified as putative Pseudomonas spp. biomarkers. The results constitute the largest dataset, described thus far for the rapid identification and characterization of Pseudomonas species present in edible fish; furthermore, these data can provide the basis for further research into the development of new therapies against these harmful pathogens.

An ecological-evolutionary perspective on the genomic diversity and habitat preferences of the Acidobacteriota

Members of the phylum Acidobacteriota inhabit a wide range of ecosystems including soils. We analyzed the global patterns of distribution and habitat preferences of various Acidobacteriota lineages across major ecosystems (soil, engineered, host-associated, marine, non-marine saline and alkaline, and terrestrial non-soil ecosystem) in 248,559 publicly available metagenomic datasets. Classes Terriglobia, Vicinamibacteria, Blastocatellia, and Thermoanaerobaculia were highly ubiquitous and showed clear preference to soil over non-soil habitats, class Polarisedimenticolia showed comparable ubiquity and preference between soil and non-soil habitats, while classes Aminicenantia and Holophagae showed preferences to non-soil habitats. However, while specific preferences were observed, most Acidobacteriota lineages were habitat generalists rather than specialists, with genomic and/or metagenomic fragments recovered from soil and non-soil habitats at various levels of taxonomic resolution. Comparative analysis of 1930 genomes strongly indicates that phylogenetic affiliation plays a more important role than the habitat from which the genome was recovered in shaping the genomic characteristics and metabolic capacities of the Acidobacteriota. The observed lack of strong habitat specialization and habitat transition driven lineage evolution in the Acidobacteriota suggest ready cross colonization between soil and non-soil habitats. We posit that such capacity is key to the successful establishment of Acidobacteriota as a major component in soil microbiomes post ecosystem disturbance events or during pedogenesis.

Comparative genomics reveals distinct diversification patterns among LysR-type transcriptional regulators in the ESKAPE pathogen Pseudomonas aeruginosa

Pseudomonas aeruginosa, a harmful nosocomial pathogen associated with cystic fibrosis and burn wounds, encodes for a large number of LysR-type transcriptional regulator proteins. To understand how and why LTTR proteins evolved with such frequency and to establish whether any relationships exist within the distribution we set out to identify the patterns underpinning LTTR distribution in P. aeruginosa and to uncover cluster-based relationships within the pangenome. Comparative genomic studies revealed that in the JGI IMG database alone ~86 000 LTTRs are present across the sequenced genomes (n=699). They are widely distributed across the species, with core LTTRs present in >93 % of the genomes and accessory LTTRs present in <7 %. Analysis showed that subsets of core LTTRs can be classified as either variable (typically specific to P. aeruginosa) or conserved (and found to be distributed in other Pseudomonas species). Extending the analysis to the more extensive Pseudomonas database, PA14 rooted analysis confirmed the diversification patterns and revealed PqsR, the receptor for the Pseudomonas quinolone signal (PQS) and 2-heptyl-4-quinolone (HHQ) quorum-sensing signals, to be amongst the most variable in the dataset. Successful complementation of the PAO1 pqsR - mutant using representative variant pqsR sequences suggests a degree of structural promiscuity within the most variable of LTTRs, several of which play a prominent role in signalling and communication. These findings provide a new insight into the diversification of LTTR proteins within the P. aeruginosa species and suggests a functional significance to the cluster, conservation and distribution patterns identified.

Classification of the plant-associated lifestyle of Pseudomonas strains using genome properties and machine learning

The rhizosphere, the region of soil surrounding roots of plants, is colonized by a unique population of Plant Growth Promoting Rhizobacteria (PGPR). Many important PGPR as well as plant pathogens belong to the genus Pseudomonas. There is, however, uncertainty on the divide between beneficial and pathogenic strains as previously thought to be signifying genomic features have limited power to separate these strains. Here we used the Genome properties (GP) common biological pathways annotation system and Machine Learning (ML) to establish the relationship between the genome wide GP composition and the plant-associated lifestyle of 91 Pseudomonas strains isolated from the rhizosphere and the phyllosphere representing both plant-associated phenotypes. GP enrichment analysis, Random Forest model fitting and feature selection revealed 28 discriminating features. A test set of 75 new strains confirmed the importance of the selected features for classification. The results suggest that GP annotations provide a promising computational tool to better classify the plant-associated lifestyle.

Cold-adaptive traits identified by comparative genomic analysis of a lipase-producing Pseudomonas sp. HS6 isolated from snow-covered soil of Sikkim Himalaya and molecular simulation of lipase for wide substrate specificity

The genomic analysis of industrially important bacteria can help in understanding their capability to withstand extreme environments and shed light on their metabolic capabilities. The whole genome of a previously reported broad temperature active lipase-producing Pseudomonas sp. HS6, isolated from snow-covered soil of the Sikkim Himalayan Region, was analyzed to understand the capability of the bacterium to withstand cold temperatures and study its lipolytic nature. Pseudomonas sp. HS6 was found to be psychrotolerant with an optimal growth temperature ranging between 25 and 30 °C, with the ability to grow at 5 °C. The genome harbours various cold-adaptation genes, such as cold-shock proteins, fatty acid alteration, and cold stress-tolerance genes, supporting the psychrotolerant nature of the organism. The comparative analysis of Pseudomonas sp. HS6 genome showed the presence of amino acid substitutions in genes that favor efficient functioning and flexibility at cold temperatures. Genome mining revealed the presence of four triacylglycerol lipases, among which the putative lipase 3 was highly similar to the broad temperature-active lipase purified and characterized in our previous study. In silico studies of putative lipase 3 revealed broad substrate specificity with partial and no inhibition of the enzyme activity in the presence of PMSF and orlistat. The presence of genes associated with cold adaptations and true lipases with activity at broad temperature and substrate specificity in the genome of Pseudomonas sp. HS6 makes this bacterium a suitable candidate for industrial applications.

Immuno-Informatics Based Peptides: An Approach for Vaccine Development Against Outer Membrane Proteins of Pseudomonas Genus

Pseudomonas genus is among the top nosocomial pathogens known to date. Being highly opportunistic, members of pseudomonas genus are most commonly connected with nosocomial infections of urinary tract and ventilator-associated pneumonia. Nevertheless, vaccine development for this pathogenic genus is slow because of no information regarding immunity correlated functional mechanism. In this present work, an immunoinformatics pipeline is used for vaccine development based on epitope-based peptide design, which can result in crucial immune response against outer membrane proteins of pseudomonas genus. A total of 127 outer membrane proteins were analysed, studied and out of them three sequences were obtained to be the producer of non-allergic, highly antigenic T-cell and B-cell epitopes which show good binding affinity towards class II HLA molecules. After performing rigorous screening utilizing docking, simulation, modelling techniques, we had one nonameric peptide (WLLATGIFL)as a good vaccine candidate. The predicted epitopes needs to be further validated for its apt use as vaccine. This work paves a new way with extensive therapeutic application against Pseudomonas genus and their associated diseases.

Pangenome inventory of Burkholderia sensu lato, Burkholderia sensu stricto, and the Burkholderia cepacia complex reveals the uniqueness of Burkholderia catarinensis

Here the pangenome analysis of Burkholderia sensu lato (s.l.) was performed for the first time, together with an updated analysis of the pangenome of Burkholderia sensu stricto, and Burkholderia cepacia complex (Bcc) focusing on the Bcc B. catarinensis specific features of its re-sequenced genome. The pangenome of Burkholderia s.l., Burkholderia s.s., and of the Bcc are open, composed of more than 96% of accessory genes, and more than 62% of unknown genes. Functional annotations showed that secondary metabolism genes belong to the variable portion of genomes, which might explain their production of several compounds with varied bioactivities. Taken together, this work shows the great variability and uniqueness of these genomes and reveals an underexplored unknown potential in poorly characterized genes. Regarding B. catarinensis 89^T, its genome harbors genes related to hydrolases production and plant growth promotion. This draft genome will be valuable for further investigation of its biotechnological potentials.

A pan-genome method to determine core regions of the Bacillus subtilis and Escherichia coli genomes

Background: Synthetic engineering of bacteria to produce industrial products is a burgeoning field of research and application. In order to optimize genome design, designers need to understand which genes are essential, which are optimal for growth, and locations in the genome that will be tolerated by the organism when inserting engineered cassettes.

Methods: We present a pan-genome based method for the identification of core regions in a genome that are strongly conserved at the species level.

Results: We show that the core regions determined by our method contain all or almost all essential genes. This demonstrates the accuracy of our method as essential genes should be core genes. We show that we outperform previous methods by this measure. We also explain why there are exceptions to this rule for our method.

Conclusions: We assert that synthetic engineers should avoid deleting or inserting into these core regions unless they understand and are manipulating the function of the genes in that region. Similarly, if the designer wishes to streamline the genome, non-core regions and in particular low penetrance genes would be good targets for deletion. Care should be taken to remove entire cassettes with similar penetrance of the genes within cassettes as they may harbor toxin/antitoxin genes which need to be removed in tandem. The bioinformatic approach introduced here saves considerable time and effort relative to knockout studies on single isolates of a given species and captures a broad understanding of the conservation of genes that are core to a species.

Dynamics of the compartmentalized Streptomyces chromosome during metabolic differentiation

Bacteria of the genus Streptomyces are prolific producers of specialized metabolites, including antibiotics. The linear chromosome includes a central region harboring core genes, as well as extremities enriched in specialized metabolite biosynthetic gene clusters. Here, we show that chromosome structure in Streptomyces ambofaciens correlates with genetic compartmentalization during exponential phase. Conserved, large and highly transcribed genes form boundaries that segment the central part of the chromosome into domains, whereas the terminal ends tend to be transcriptionally quiescent compartments with different structural features. The onset of metabolic differentiation is accompanied by a rearrangement of chromosome architecture, from a rather 'open' to a 'closed' conformation, in which highly expressed specialized metabolite biosynthetic genes form new boundaries. Thus, our results indicate that the linear chromosome of S. ambofaciens is partitioned into structurally distinct entities, suggesting a link between chromosome folding, gene expression and genome evolution.

Comparative genomic insights into the epidemiology and virulence of plant pathogenic pseudomonads from Turkey

Pseudomonas is a highly diverse genus that includes species that cause disease in both plants and animals. Recently, pathogenic pseudomonads from the Pseudomonas syringae and Pseudomonas fluorescens species complexes have caused significant outbreaks in several agronomically important crops in Turkey, including tomato, citrus, artichoke and melon. We characterized 169 pathogenic Pseudomonas strains associated with recent outbreaks in Turkey via multilocus sequence analysis and whole-genome sequencing, then used comparative and evolutionary genomics to characterize putative virulence mechanisms. Most of the isolates are closely related to other plant pathogens distributed among the primary phylogroups of P. syringae, although there are significant numbers of P. fluorescens isolates, which is a species better known as a rhizosphere-inhabiting plant-growth promoter. We found that all 39 citrus blast pathogens cluster in P. syringae phylogroup 2, although strains isolated from the same host do not cluster monophyletically, with lemon, mandarin orange and sweet orange isolates all being intermixed throughout the phylogroup. In contrast, 20 tomato pith pathogens are found in two independent lineages: one in the P. syringae secondary phylogroups, and the other from the P. fluorescens species complex. These divergent pith necrosis strains lack characteristic virulence factors like the canonical tripartite type III secretion system, large effector repertoires and the ability to synthesize multiple bacterial phytotoxins, suggesting they have alternative molecular mechanisms to cause disease. These findings highlight the complex nature of host specificity among plant pathogenic pseudomonads.

A pan-genome method to determine core regions of the Bacillus subtilis and Escherichia coli genomes

Background: Synthetic engineering of bacteria to produce industrial products is a burgeoning field of research and application. In order to optimize genome design, designers need to understand which genes are essential, which are optimal for growth, and locations in the genome that will be tolerated by the organism when inserting engineered cassettes. Methods: We present a pan-genome based method for the identification of core regions in a genome that are strongly conserved at the species level. Results: We show that the core regions determined by our method contain all or almost all essential genes. This demonstrates the accuracy of our method as essential genes should be core genes. We show that we outperform previous methods by this measure. We also explain why there are exceptions to this rule for our method. Conclusions: We assert that synthetic engineers should avoid deleting or inserting into these core regions unless they understand and are manipulating the function of the genes in that region. Similarly, if the designer wishes to streamline the genome, non-core regions and in particular low penetrance genes would be good targets for deletion. Care should be taken to remove entire cassettes with similar penetrance of the genes within cassettes as they may harbor toxin/antitoxin genes which need to be removed in tandem. The bioinformatic approach introduced here saves considerable time and effort relative to knockout studies on single isolates of a given species and captures a broad understanding of the conservation of genes that are core to a species.

A navigation guide of synthetic biology tools for Pseudomonas putida

Pseudomonas putida is a microbial chassis of huge potential for industrial and environmental biotechnology, owing to its remarkable metabolic versatility and ability to sustain difficult redox reactions and operational stresses, among other attractive characteristics. A wealth of genetic and in silico tools have been developed to enable the unravelling of its physiology and improvement of its performance. However, the rise of this microbe as a promising platform for biotechnological applications has resulted in diversification of tools and methods rather than standardization and convergence. As a consequence, multiple tools for the same purpose have been generated, whilst most of them have not been embraced by the scientific community, which has led to compartmentalization and inefficient use of resources. Inspired by this and by the substantial increase in popularity of P. putida, we aim herein to bring together and assess all currently available (wet and dry) synthetic biology tools specific for this microbe, focusing on the last 5 years. We provide information on the principles, functionality, advantages and limitations, with special focus on their use in metabolic engineering. Additionally, we compare the tool portfolio for P. putida with those for other bacterial chassis and discuss potential future directions for tool development. Therefore, this review is intended as a reference guide for experts and new 'users' of this promising chassis.

Phenotype and multi-omics comparison of Staphylococcus and Streptococcus uncovers pathogenic traits and predicts zoonotic potential

BACKGROUND

Staphylococcus and Streptococcus species can cause many different diseases, ranging from mild skin infections to life-threatening necrotizing fasciitis. Both genera consist of commensal species that colonize the skin and nose of humans and animals, and of which some can display a pathogenic phenotype.

RESULTS

We compared 235 Staphylococcus and 315 Streptococcus genomes based on their protein domain content. We show the relationships between protein persistence and essentiality by integrating essentiality predictions from two metabolic models and essentiality measurements from six large-scale transposon mutagenesis experiments. We identified clusters of strains within species based on proteins associated to similar biological processes. We built Random Forest classifiers that predicted the zoonotic potential. Furthermore, we identified shared attributes between of Staphylococcus aureus and Streptococcus pyogenes that allow them to cause necrotizing fasciitis.

CONCLUSIONS

Differences observed in clustering of strains based on functional groups of proteins correlate with phenotypes such as host tropism, capability to infect multiple hosts and drug resistance. Our method provides a solid basis towards large-scale prediction of phenotypes based on genomic information.

A metabolic and physiological design study of Pseudomonas putida KT2440 capable of anaerobic respiration

BACKGROUND

Pseudomonas putida KT2440 is a metabolically versatile, HV1-certified, genetically accessible, and thus interesting microbial chassis for biotechnological applications. However, its obligate aerobic nature hampers production of oxygen sensitive products and drives up costs in large scale fermentation. The inability to perform anaerobic fermentation has been attributed to insufficient ATP production and an inability to produce pyrimidines under these conditions. Addressing these bottlenecks enabled growth under micro-oxic conditions but does not lead to growth or survival under anoxic conditions.

RESULTS

Here, a data-driven approach was used to develop a rational design for a P. putida KT2440 derivative strain capable of anaerobic respiration. To come to the design, data derived from a genome comparison of 1628 Pseudomonas strains was combined with genome-scale metabolic modelling simulations and a transcriptome dataset of 47 samples representing 14 environmental conditions from the facultative anaerobe Pseudomonas aeruginosa.

CONCLUSIONS

The results indicate that the implementation of anaerobic respiration in P. putida KT2440 would require at least 49 additional genes of known function, at least 8 genes encoding proteins of unknown function, and 3 externally added vitamins.

Metabolic reconstruction of Pseudomonas chlororaphis ATCC 9446 to understand its metabolic potential as a phenazine-1-carboxamide-producing strain

Pseudomonas chlororaphis is a plant-associated bacterium with reported antagonistic activity against different organisms and plant growth-promoting properties. P. chlororaphis possesses exciting biotechnological features shared with another Pseudomonas with a nonpathogenic phenotype. Part of the antagonistic role of P. chlororaphis is due to its production of a wide variety of phenazines. To expand the knowledge of the metabolic traits of this organism, we constructed the first experimentally validated genome-scale model of P. chlororaphis ATCC 9446, containing 1267 genes and 2289 reactions, and analyzed strategies to maximize its potential for the production of phenazine-1-carboxamide (PCN). The resulting model also describes the capability of P. chlororaphis to carry out the denitrification process and its ability to consume sucrose (Scr), trehalose, mannose, and galactose as carbon sources. Additionally, metabolic network analysis suggested fatty acids as the best carbon source for PCN production. Moreover, the optimization of PCN production was performed with glucose and glycerol. The optimal PCN production phenotype requires an increased carbon flux in TCA and glutamine synthesis. Our simulations highlight the intrinsic H₂O₂ flux associated with PCN production, which may generate cellular stress in an overproducing strain. These results suggest that an improved antioxidative strategy could lead to optimal performance of phenazine-producing strains of P. chlororaphis. KEY POINTS : • This is the first publication of a metabolic model for a strain of P. chlororaphis. • Genome-scale model is worthy tool to increase the knowledge of a non model organism. • Fluxes simulations indicate a possible effect of H₂O₂ on phenazines production. • P. chlororaphis can be a suitable model for a wide variety of compounds.

The Great ESKAPE: Exploring the Crossroads of Bile and Antibiotic Resistance in Bacterial Pathogens

Throughout the course of infection, many pathogens encounter bactericidal conditions that threaten the viability of the bacteria and impede the establishment of infection. Bile is one of the most innately bactericidal compounds present in humans, functioning to reduce the bacterial burden in the gastrointestinal tract while also aiding in digestion. It is becoming increasingly apparent that pathogens successfully resist the bactericidal conditions of bile, including bacteria that do not normally cause gastrointestinal infections. This review highlights the ability of Enterococcus, Staphylococcus, Klebsiella, Acinetobacter, Pseudomonas, Enterobacter (ESKAPE), and other enteric pathogens to resist bile and how these interactions can impact the sensitivity of bacteria to various antimicrobial agents. Given that pathogen exposure to bile is an essential component to gastrointestinal transit that cannot be avoided, understanding how bile resistance mechanisms align with antimicrobial resistance is vital to our ability to develop new, successful therapeutics in an age of widespread and increasing antimicrobial resistance.

Combined Transcriptome Sequencing of Mycoplasma hyopneumoniae and Infected Pig Lung Tissue Reveals Up-Regulation of Bacterial F1-Like ATPase and Down-Regulation of the P102 Cilium Adhesin in vivo

Mycoplasma hyopneumoniae (M. hyopneumoniae) causes enzootic pneumonia in pigs but it is still largely unknown which host-pathogen interactions enable persistent infection and cause disease. In this study, we analyzed the host and bacterial transcriptomes during infection using RNA sequencing. Comparison of the transcriptome of lung lesion tissue from infected pigs with lung tissue from non-infected animals, identified 424 differentially expressed genes (FDR < 0.01 and fold change > 1.5LOG2). These genes were part of the following major pathways of the immune system: interleukin signaling (type 4, 10, 13, and 18), regulation of Toll-like receptors by endogenous ligand and activation of C3 and C5 in the complement system. Besides analyzing the lung transcriptome, a sampling protocol was developed to obtain enough bacterial mRNA from infected lung tissue for RNA sequencing. This was done by flushing infected lobes in the lung, and subsequently enriching for bacterial RNA. On average, 2.2 million bacterial reads were obtained per biological replicate to analyze the bacterial in vivo transcriptome. We compared the in vivo bacterial transcriptome with the transcriptome of bacteria grown in vitro and identified 22 up-regulated and 30 down-regulated genes (FDR < 0.01 and fold change > 2LOG2). Six out of seven genes in the operon encoding the mycoplasma specific F1-like ATPase (MHP_RS02445-MHP_RS02475) and all genes in the operon MHP_RS01965-MHP_RS01990 with functions related to nucleotide metabolism, spermidine transport and glycerol-3-phoshate transport were up-regulated in vivo. Down-regulated in vivo were genes related to glycerol uptake, cilium adhesion (P102), cell division and myo-inositol metabolism. In addition to providing a novel method to isolate bacterial mRNA from infected lung, this study provided insights into changes in gene expression during infection, which could help development of novel treatment strategies against enzootic pneumonia caused by M. hyopneumoniae.

Comparative Analysis of the Core Proteomes among the Pseudomonas Major Evolutionary Groups Reveals Species-Specific Adaptations for Pseudomonas aeruginosa and Pseudomonas chlororaphis

The Pseudomonas genus includes many species living in diverse environments and hosts. It is important to understand which are the major evolutionary groups and what are the genomic/proteomic components they have in common or are unique. Towards this goal, we analyzed 494 complete Pseudomonas proteomes and identified 297 core-orthologues. The subsequent phylogenomic analysis revealed two well-defined species (Pseudomonas aeruginosa and Pseudomonas chlororaphis) and four wider phylogenetic groups (Pseudomonas fluorescens, Pseudomonas stutzeri, Pseudomonas syringae, Pseudomonas putida) with a sufficient number of proteomes. As expected, the genus-level core proteome was highly enriched for proteins involved in metabolism, translation, and transcription. In addition, between 39–70% of the core proteins in each group had a significant presence in each of all the other groups. Group-specific core proteins were also identified, with P. aeruginosa having the highest number of these and P. fluorescens having none. We identified several P. aeruginosa-specific core proteins (such as CntL, CntM, PlcB, Acp1, MucE, SrfA, Tse1, Tsi2, Tse3, and EsrC) that are known to play an important role in its pathogenicity. Finally, a holin family bacteriocin and a mitomycin-like biosynthetic protein were found to be core-specific for P. cholororaphis and we hypothesize that these proteins may confer a competitive advantage against other root-colonizers.

Presence of Legionella spp. in cooling towers: the role of microbial diversity, Pseudomonas, and continuous chlorine application

Legionnaires' disease (LD) is a severe pneumonia caused by several species of the genus Legionella, most frequently by Legionella pneumophila. Cooling towers are the most common source for large community-associated outbreaks. Colonization, survival, and proliferation of L. pneumophila in cooling towers are necessary for outbreaks to occur. These steps are affected by the chemical and physical parameters of the cooling tower environment. We hypothesize that the bacterial community residing in the cooling tower could also affect the presence of L. pneumophila. A 16S rRNA gene targeted amplicon sequencing approach was used to study the bacterial community of cooling towers and its relationship with the Legionella spp. and L. pneumophila communities. The results indicated that the water source shaped the bacterial community of cooling towers. Several taxa were enriched and positively correlated with Legionella spp. and L. pneumophila. In contrast, Pseudomonas showed a strong negative correlation with Legionella spp. and several other genera. Most importantly, continuous chlorine application reduced microbial diversity and promoted the presence of Pseudomonas creating a non-permissive environment for Legionella spp. This suggests that disinfection strategies as well as the resident microbial population influences the ability of Legionella spp. to colonize cooling towers.

Genome-guided analysis allows the identification of novel physiological traits in Trichococcus species

BACKGROUND

The genus Trichococcus currently contains nine species: T. flocculiformis, T. pasteurii, T. palustris, T. collinsii, T. patagoniensis, T. ilyis, T. paludicola, T. alkaliphilus, and T. shcherbakoviae. In general, Trichococcus species can degrade a wide range of carbohydrates. However, only T. pasteurii and a non-characterized strain of Trichococcus, strain ES5, have the capacity of converting glycerol to mainly 1,3-propanediol. Comparative genomic analysis of Trichococcus species provides the opportunity to further explore the physiological potential and uncover novel properties of this genus.

RESULTS

In this study, a genotype-phenotype comparative analysis of Trichococcus strains was performed. The genome of Trichococcus strain ES5 was sequenced and included in the comparison with the other nine type strains. Genes encoding functions related to e.g. the utilization of different carbon sources (glycerol, arabinan and alginate), antibiotic resistance, tolerance to low temperature and osmoregulation could be identified in all the sequences analysed. T. pasteurii and Trichococcus strain ES5 contain a operon with genes encoding necessary enzymes for 1,3-PDO production from glycerol. All the analysed genomes comprise genes encoding for cold shock domains, but only five of the Trichococcus species can grow at 0 °C. Protein domains associated to osmoregulation mechanisms are encoded in the genomes of all Trichococcus species, except in T. palustris, which had a lower resistance to salinity than the other nine studied Trichococcus strains.

CONCLUSIONS

Genome analysis and comparison of ten Trichococcus strains allowed the identification of physiological traits related to substrate utilization and environmental stress resistance (e.g. to cold and salinity). Some substrates were used by single species, e.g. alginate by T. collinsii and arabinan by T. alkaliphilus. Strain ES5 may represent a subspecies of Trichococcus flocculiformis and contrary to the type strain (DSM 2094T), is able to grow on glycerol with the production of 1,3-propanediol.

The Empusa code generator and its application to GBOL, an extendable ontology for genome annotation

The RDF data model facilitates integration of diverse data available in structured and semi-structured formats. To obtain a coherent RDF graph the chosen ontology must be consistently applied. However, addition of new diverse data causes the ontology to evolve, which could lead to accumulation of unintended erroneous composites. Thus, there is a need for a gate keeping system that compares the intended content described in the ontology with the actual content of the resource. The Empusa code generator facilitates creation of composite RDF resources from disparate sources. Empusa can convert a schema into an associated application programming interface (API), that can be used to perform data consistency checks and generates Markdown documentation to make persistent URLs resolvable. Using Empusa consistency is ensured within and between the ontology and the content of the resource. As an illustration of the potential of Empusa, we present the Genome Biology Ontology Language (GBOL). GBOL uses and extends current ontologies to provide a formal representation of genomic entities, along with their properties, relations and provenance.

In silico-guided engineering of Pseudomonas putida towards growth under micro-oxic conditions

BACKGROUND

Pseudomonas putida is a metabolically versatile, genetically accessible, and stress-robust species with outstanding potential to be used as a workhorse for industrial applications. While industry recognises the importance of robustness under micro-oxic conditions for a stable production process, the obligate aerobic nature of P. putida, attributed to its inability to produce sufficient ATP and maintain its redox balance without molecular oxygen, severely limits its use for biotechnology applications.

RESULTS

Here, a combination of genome-scale metabolic modelling and comparative genomics is used to pinpoint essential [Formula: see text]-dependent processes. These explain the inability of the strain to grow under anoxic conditions: a deficient ATP generation and an inability to synthesize essential metabolites. Based on this, several P. putida recombinant strains were constructed harbouring acetate kinase from Escherichia coli for ATP production, and a class I dihydroorotate dehydrogenase and a class III anaerobic ribonucleotide triphosphate reductase from Lactobacillus lactis for the synthesis of essential metabolites. Initial computational designs were fine-tuned by means of adaptive laboratory evolution.

CONCLUSIONS

We demonstrated the value of combining in silico approaches, experimental validation and adaptive laboratory evolution for microbial design by making the strictly aerobic Pseudomonas putida able to grow under micro-oxic conditions.

Risk-Based Bioengineering Strategies for Reliable Bacterial Vaccine Production

Design of a reliable process for bacterial antigen production requires understanding of and control over critical process parameters. Current methods for process design use extensive screening experiments for determining ranges of critical process parameters yet fail to give clear insights into how they influence antigen potency. To address this gap, we propose to apply constraint-based, genome-scale metabolic models to reduce the need of experimental screening for strain selection and to optimize strains based on model driven iterative Design-Build-Test-Learn (DBTL) cycles. Application of these systematic methods has not only increased the understanding of how metabolic network properties influence antigen potency, but also allows identification of novel critical process parameters that need to be controlled to achieve high process reliability.

SAPP: functional genome annotation and analysis through a semantic framework using FAIR principles

Summary

To unlock the full potential of genome data and to enhance data interoperability and reusability of genome annotations we have developed SAPP, a Semantic Annotation Platform with Provenance. SAPP is designed as an infrastructure supporting FAIR de novo computational genomics but can also be used to process and analyze existing genome annotations. SAPP automatically predicts, tracks and stores structural and functional annotations and associated dataset- and element-wise provenance in a Linked Data format, thereby enabling information mining and retrieval with Semantic Web technologies. This greatly reduces the administrative burden of handling multiple analysis tools and versions thereof and facilitates multi-level large scale comparative analysis.

Availability and implementation

SAPP is written in JAVA and freely available at https://gitlab.com/sapp and runs on Unix-like operating systems. The documentation, examples and a tutorial are available at https://sapp.gitlab.io.

Microbial strains isolated from CO2-venting Kolumbo submarine volcano show enhanced co-tolerance to acidity and antibiotics

As ocean acidification intensifies, there is growing global concern about the impacts that future pH levels are likely to have on marine life and ecosystems. By analogy, a steep decrease of seawater pH with depth is encountered inside the Kolumbo submarine volcano (northeast Santorini) as a result of natural CO₂ venting, making this system ideal for ocean acidification research. Here, we investigated whether the increase of acidity towards deeper layers of Kolumbo crater had any effect on relevant phenotypic traits of bacterial isolates. A total of 31 Pseudomonas strains were isolated from both surface- (SSL) and deep-seawater layers (DSL), with the latter presenting a significantly higher acid tolerance. In particular, the DSL strains were able to cope with H⁺ levels that were 18 times higher. Similarly, the DSL isolates exhibited a significantly higher tolerance than SSL strains against six commonly used antibiotics and As(III). More importantly, a significant positive correlation was revealed between antibiotics and acid tolerance across the entire set of SSL and DSL isolates. Our findings imply that Pseudomonas species with higher resilience to antibiotics could be favored by the prospect of acidifying oceans. Further studies are required to determine if this feature is universal across marine bacteria and to assess potential ecological impacts.

CiliateGEM: an open-project and a tool for predictions of ciliate metabolic variations and experimental condition design

BACKGROUND

The study of cell metabolism is becoming central in several fields such as biotechnology, evolution/adaptation and human disease investigations. Here we present CiliateGEM, the first metabolic network reconstruction draft of the freshwater ciliate Tetrahymena thermophila. We also provide the tools and resources to simulate different growth conditions and to predict metabolic variations. CiliateGEM can be extended to other ciliates in order to set up a meta-model, i.e. a metabolic network reconstruction valid for all ciliates. Ciliates are complex unicellular eukaryotes of presumably monophyletic origin, with a phylogenetic position that is equal from plants and animals. These cells represent a new concept of unicellular system with a high degree of species, population biodiversity and cell complexity. Ciliates perform in a single cell all the functions of a pluricellular organism, including locomotion, feeding, digestion, and sexual processes.

RESULTS

After generating the model, we performed an in-silico simulation with the presence and absence of glucose. The lack of this nutrient caused a 32.1% reduction rate in biomass synthesis. Despite the glucose starvation, the growth did not stop due to the use of alternative carbon sources such as amino acids.

CONCLUSIONS

The future models obtained from CiliateGEM may represent a new approach to describe the metabolism of ciliates. This tool will be a useful resource for the ciliate research community in order to extend these species as model organisms in different research fields. An improved understanding of ciliate metabolism could be relevant to elucidate the basis of biological phenomena like genotype-phenotype relationships, population genetics, and cilia-related disease mechanisms.

Comparative Genomics Highlights Symbiotic Capacities and High Metabolic Flexibility of the Marine Genus Pseudovibrio

Pseudovibrio is a marine bacterial genus members of which are predominantly isolated from sessile marine animals, and particularly sponges. It has been hypothesized that Pseudovibrio spp. form mutualistic relationships with their hosts. Here, we studied Pseudovibrio phylogeny and genetic adaptations that may play a role in host colonization by comparative genomics of 31 Pseudovibrio strains, including 25 sponge isolates. All genomes were highly similar in terms of encoded core metabolic pathways, albeit with substantial differences in overall gene content. Based on gene composition, Pseudovibrio spp. clustered by geographic region, indicating geographic speciation. Furthermore, the fact that isolates from the Mediterranean Sea clustered by sponge species suggested host-specific adaptation or colonization. Genome analyses suggest that Pseudovibrio hongkongensis UST20140214-015BT is only distantly related to other Pseudovibrio spp., thereby challenging its status as typical Pseudovibrio member. All Pseudovibrio genomes were found to encode numerous proteins with SEL1 and tetratricopeptide repeats, which have been suggested to play a role in host colonization. For evasion of the host immune system, Pseudovibrio spp. may depend on type III, IV, and VI secretion systems that can inject effector molecules into eukaryotic cells. Furthermore, Pseudovibrio genomes carry on average seven secondary metabolite biosynthesis clusters, reinforcing the role of Pseudovibrio spp. as potential producers of novel bioactive compounds. Tropodithietic acid, bacteriocin, and terpene biosynthesis clusters were highly conserved within the genus, suggesting an essential role in survival, for example through growth inhibition of bacterial competitors. Taken together, these results support the hypothesis that Pseudovibrio spp. have mutualistic relations with sponges.

Fermentation products in the cystic fibrosis airways induce aggregation and dormancy-associated expression profiles in a CF clinical isolate of Pseudomonas aeruginosa

Pseudomonas aeruginosa is a well-known dominant opportunistic pathogen in cystic fibrosis (CF) with a wide range of metabolic capacities. However, P. aeruginosa does not colonize the airways alone, and benefits from the metabolic products of neighboring cells-especially volatile molecules that can travel between different parts of the airways easily. Here, we present a study that investigates the metabolic, gene expression profiles and phenotypic responses of a P. aeruginosa clinical isolate to fermentation products lactic acid and 2,3-butanediol, metabolites that are produced by facultative anaerobic members of the CF polymicrobial community and potential biomarkers of disease progression. Although previous studies have successfully investigated the metabolic and transcriptional profiles of P. aeruginosa, most have used common lab reference strains that may differ in important ways from clinical isolates. Using transcriptomics and metabolomics with gas chromatography time of flight mass spectrometry, we observe that fermentation products induce pyocyanin production along with the expression of genes involved in P. aeruginosa amino acid utilization, dormancy and aggregative or biofilm modes of growth. These findings have important implications for how interactions within the diverse CF microbial community influence microbial physiology, with potential clinical consequences.

In silico engineering of Pseudomonas metabolism reveals new biomarkers for increased biosurfactant production

BACKGROUND

Rhamnolipids, biosurfactants with a wide range of biomedical applications, are amphiphilic molecules produced on the surfaces of or excreted extracellularly by bacteria including Pseudomonas aeruginosa. However, Pseudomonas putida is a non-pathogenic model organism with greater metabolic versatility and potential for industrial applications.

METHODS

We investigate in silico the metabolic capabilities of P. putida for rhamnolipids biosynthesis using statistical, metabolic and synthetic engineering approaches after introducing key genes (RhlA and RhlB) from P. aeruginosa into a genome-scale model of P. putida. This pipeline combines machine learning methods with multi-omic modelling, and drives the engineered P. putida model toward an optimal production and export of rhamnolipids out of the membrane.

RESULTS

We identify a substantial increase in synthesis of rhamnolipids by the engineered model compared to the control model. We apply statistical and machine learning techniques on the metabolic reaction rates to identify distinct features on the structure of the variables and individual components driving the variation of growth and rhamnolipids production. We finally provide a computational framework for integrating multi-omics data and identifying latent pathways and genes for the production of rhamnolipids in P. putida.

CONCLUSIONS

We anticipate that our results will provide a versatile methodology for integrating multi-omics data for topological and functional analysis of P. putida toward maximization of biosurfactant production.

Metabolic modeling of energy balances in Mycoplasma hyopneumoniae shows that pyruvate addition increases growth rate

Mycoplasma hyopneumoniae is cultured on large-scale to produce antigen for inactivated whole-cell vaccines against respiratory disease in pigs. However, the fastidious nutrient requirements of this minimal bacterium and the low growth rate make it challenging to reach sufficient biomass yield for antigen production. In this study, we sequenced the genome of M. hyopneumoniae strain 11 and constructed a high quality constraint-based genome-scale metabolic model of 284 chemical reactions and 298 metabolites. We validated the model with time-series data of duplicate fermentation cultures to aim for an integrated model describing the dynamic profiles measured in fermentations. The model predicted that 84% of cellular energy in a standard M. hyopneumoniae cultivation was used for non-growth associated maintenance and only 16% of cellular energy was used for growth and growth associated maintenance. Following a cycle of model-driven experimentation in dedicated fermentation experiments, we were able to increase the fraction of cellular energy used for growth through pyruvate addition to the medium. This increase in turn led to an increase in growth rate and a 2.3 times increase in the total biomass concentration reached after 3-4 days of fermentation, enhancing the productivity of the overall process. The model presented provides a solid basis to understand and further improve M. hyopneumoniae fermentation processes. Biotechnol. Bioeng. 2017;114: 2339-2347. © 2017 Wiley Periodicals, Inc.

Draft Genome Sequence of Actinomyces succiniciruminis Strain Am4 T , Isolated from Cow Rumen Fluid

Actinomyces succiniciruminis strain Am4 T , isolated from cow rumen fluid, can metabolize a range of substrates including complex carbohydrates to organic acids. Here, we report a 3.33-Mbp draft genome of Actinomyces succiniciruminis .

Concurrent Haloalkanoate Degradation and Chlorate Reduction by Pseudomonas chloritidismutans AW-1T

Haloalkanoates are environmental pollutants that can be degraded aerobically by microorganisms producing hydrolytic dehalogenases. However, there is a lack of information about the anaerobic degradation of haloalkanoates. Genome analysis of Pseudomonas chloritidismutans AW-1^T, a facultative anaerobic chlorate-reducing bacterium, showed the presence of two putative haloacid dehalogenase genes, the l-DEX gene and dehI, encoding an l-2-haloacid dehalogenase (l-DEX) and a halocarboxylic acid dehydrogenase (DehI), respectively. Hence, we studied the concurrent degradation of haloalkanoates and chlorate as a yet-unexplored trait of strain AW-1^T The deduced amino acid sequences of l-DEX and DehI revealed 33 to 37% and 26 to 86% identities with biochemically/structurally characterized l-DEX and the d- and dl-2-haloacid dehalogenase enzymes, respectively. Physiological experiments confirmed that strain AW-1^T can grow on chloroacetate, bromoacetate, and both l- and d-α-halogenated propionates with chlorate as an electron acceptor. Interestingly, growth and haloalkanoate degradation were generally faster with chlorate as an electron acceptor than with oxygen as an electron acceptor. In line with this, analyses of l-DEX and DehI dehalogenase activities using cell-free extract (CFE) of strain AW-1^T grown on dl-2-chloropropionate under chlorate-reducing conditions showed up to 3.5-fold higher dehalogenase activity than the CFE obtained from AW-1^T cells grown on dl-2-chloropropionate under aerobic conditions. Reverse transcription-quantitative PCR showed that the l-DEX gene was expressed constitutively independently of the electron donor (haloalkanoates or acetate) or acceptor (chlorate or oxygen), whereas the expression of dehI was induced by haloalkanoates. Concurrent degradation of organic and inorganic halogenated compounds by strain AW-1^T represents a unique metabolic capacity in a single bacterium, providing a new piece of the puzzle of the microbial halogen cycle.IMPORTANCE Halogenated organic and inorganic compounds are important environmental pollutants that have carcinogenic and genotoxic effects on both animals and humans. Previous research studied the degradation of organic and inorganic halogenated compounds separately but not concurrently. This study shows concurrent degradation of halogenated alkanoates and chlorate as an electron donor and acceptor, respectively, coupled to growth in a single bacterium, Pseudomonas chloritidismutans AW-1^T Hence, besides biogenesis of molecular oxygen from chlorate reduction enabling a distinctive placement of strain AW-1^T between aerobic and anaerobic microorganisms, we can now add another unique metabolic potential of this bacterium to the roster. The degradation of different halogenated compounds under anoxic conditions by a single bacterium is also of interest for the natural halogen cycle in different aquatic and terrestrial ecosystems where ample natural production of halogenated compounds has been documented.

Persistence of Functional Protein Domains in Mycoplasma Species and their Role in Host Specificity and Synthetic Minimal Life

Mycoplasmas are the smallest self-replicating organisms and obligate parasites of a specific vertebrate host. An in-depth analysis of the functional capabilities of mycoplasma species is fundamental to understand how some of simplest forms of life on Earth succeeded in subverting complex hosts with highly sophisticated immune systems. In this study we present a genome-scale comparison, focused on identification of functional protein domains, of 80 publically available mycoplasma genomes which were consistently re-annotated using a standardized annotation pipeline embedded in a semantic framework to keep track of the data provenance. We examined the pan- and core-domainome and studied predicted functional capability in relation to host specificity and phylogenetic distance. We show that the pan- and core-domainome of mycoplasma species is closed. A comparison with the proteome of the “minimal” synthetic bacterium JCVI-Syn3.0 allowed us to classify domains and proteins essential for minimal life. Many of those essential protein domains, essential Domains of Unknown Function (DUFs) and essential hypothetical proteins are not persistent across mycoplasma genomes suggesting that mycoplasma species support alternative domain configurations that bypass their essentiality. Based on the protein domain composition, we could separate mycoplasma species infecting blood and tissue. For selected genomes of tissue infecting mycoplasmas, we could also predict whether the host is ruminant, pig or human. Functionally closely related mycoplasma species, which have a highly similar protein domain repertoire, but different hosts could not be separated. This study provides a concise overview of the functional capabilities of mycoplasma species, which can be used as a basis to further understand host-pathogen interaction or to design synthetic minimal life.