Refining the transcriptional landscapes for distinct clades of virulent phages infecting Pseudomonas aeruginosa.

The introduction of high-throughput sequencing has resulted in a surge of available bacteriophage genomes, unveiling their tremendous genomic diversity. However, our current understanding of the complex transcriptional mechanisms that dictate their gene expression during infection is limited to a handful of model phages. Here, we applied ONT-cappable-seq to reveal the transcriptional architecture of six different clades of virulent phages infecting Pseudomonas aeruginosa. This long-read microbial transcriptomics approach is tailored to globally map transcription start and termination sites, transcription units, and putative RNA-based regulators on dense phage genomes. Specifically, the full-length transcriptomes of LUZ19, LUZ24, 14-1, YuA, PAK_P3, and giant phage phiKZ during early, middle, and late infection were collectively charted. Beyond pinpointing traditional promoter and terminator elements and transcription units, these transcriptional profiles provide insights in transcriptional attenuation and splicing events and allow straightforward validation of Group I intron activity. In addition, ONT-cappable-seq data can guide genome-wide discovery of novel regulatory element candidates, including noncoding RNAs and riboswitches. This work substantially expands the number of annotated phage-encoded transcriptional elements identified to date, shedding light on the intricate and diverse gene expression regulation mechanisms in Pseudomonas phages, which can ultimately be sourced as tools for biotechnological applications in phage and bacterial engineering.

Posttranslational modifications in bacteria during phage infection.

During phage infection, both virus and bacteria attempt to gain and/or maintain control over critical bacterial functions, through a plethora of strategies. These strategies include posttranslational modifications (PTMs, including phosphorylation, ribosylation, and acetylation), as rapid and dynamic regulators of protein behavior. However, to date, knowledge on the topic remains scarce and fragmented, while a more systematic investigation lies within reach. The release of AlphaFold, which advances PTM enzyme discovery and functional elucidation, and the increasing inclusivity and scale of mass spectrometry applications to new PTM types, could significantly accelerate research in the field. In this review, we highlight the current knowledge on PTMs during phage infection, and conceive a possible pipeline for future research, following an enzyme-target-function scheme.

FLAMS: Find Lysine Acylations and other Modification Sites.

SUMMARY: Today, hundreds of post-translational modification (PTM) sites are routinely identified at once, but the comparison of new experimental datasets to already existing ones is hampered by the current inability to search most PTM databases at the protein residue level. We present FLAMS (Find Lysine Acylations and other Modification Sites), a Python3-based command line and web-tool that enables researchers to compare their PTM sites to the contents of the CPLM, the largest dedicated protein lysine modification database, and dbPTM, the most comprehensive general PTM database, at the residue level. FLAMS can be integrated into PTM analysis pipelines, allowing researchers to quickly assess the novelty and conservation of PTM sites across species in newly generated datasets, aiding in the functional assessment of sites and the prioritization of sites for further experimental characterization. AVAILABILITY AND IMPLEMENTATION: FLAMS is implemented in Python3, and freely available under an MIT license. It can be found as a command line tool at https://github.com/hannelorelongin/FLAMS, pip and conda; and as a web service at https://www.biw.kuleuven.be/m2s/cmpg/research/CSB/tools/flams/.

Molecular Evolution of Aryl Hydrocarbon Receptor Signaling Pathway Genes.

The Aryl hydrocarbon receptor is an ancient transcriptional factor originally discovered as a sensor of dioxin. In addition to its function as a receptor of environmental toxicants, it plays an important role in development. Although a significant amount of research has been carried out to understand the AHR signal transduction pathway and its involvement in species' susceptibility to environmental toxicants, none of them to date has comprehensively studied its evolutionary origins. Studying the evolutionary origins of molecules can inform ancestral relationships of genes. The vertebrate genome has been shaped by two rounds of whole-genome duplications (WGD) at the base of vertebrate evolution approximately 600 million years ago, followed by lineage-specific gene losses, which often complicate the assignment of orthology. It is crucial to understand the evolutionary origins of this transcription factor and its partners, to distinguish orthologs from ancient non-orthologous homologs. In this study, we have investigated the evolutionary origins of proteins involved in the AHR pathway. Our results provide evidence of gene loss and duplications, crucial for understanding the functional connectivity of humans and model species. Multiple studies have shown that 2R-ohnologs (genes and proteins that have survived from the 2R-WGD) are enriched in signaling components relevant to developmental disorders and cancer. Our findings provide a link between the AHR pathway's evolutionary trajectory and its potential mechanistic involvement in pathogenesis.

Genomic diversity and metabolic potential of marine Pseudomonadaceae.

Recent changes in the taxonomy of the Pseudomonadaceae family have led to the delineation of three new genera (Atopomonas, Halopseudomonas and Stutzerimonas). However, the genus Pseudomonas remains the most densely populated and displays a broad genetic diversity. Pseudomonas are able to produce a wide variety of secondary metabolites which drives important ecological functions and have a great impact in sustaining their lifestyles. While soilborne Pseudomonas are constantly examined, we currently lack studies aiming to explore the genetic diversity and metabolic potential of marine Pseudomonas spp. In this study, 23 Pseudomonas strains were co-isolated with Vibrio strains from three marine microalgal cultures and rpoD-based phylogeny allowed their assignment to the Pseudomonas oleovorans group (Pseudomonas chengduensis, Pseudomonas toyotomiensis and one new species). We combined whole genome sequencing on three selected strains with an inventory of marine Pseudomonas genomes to assess their phylogenetic assignations and explore their metabolic potential. Our results revealed that most strains are incorrectly assigned at the species level and half of them do not belong to the genus Pseudomonas but instead to the genera Halopseudomonas or Stutzerimonas. We highlight the presence of 26 new species (Halopseudomonas (n = 5), Stutzerimonas (n = 7) and Pseudomonas (n = 14)) and describe one new species, Pseudomonas chaetocerotis sp. nov. (type strain 536T = LMG 31766T = DSM 111343T). We used genome mining to identify numerous BGCs coding for the production of diverse known metabolites (i.e., osmoprotectants, photoprotectants, quorum sensing molecules, siderophores, cyclic lipopeptides) but also unknown metabolites (e.g., ARE, hybrid ARE-DAR, siderophores, orphan NRPS gene clusters) awaiting chemical characterization. Finally, this study underlines that marine environments host a huge diversity of Pseudomonadaceae that can drive the discovery of new secondary metabolites.

DARTpaths, an in silico platform to investigate molecular mechanisms of compounds.

SUMMARY: Xpaths is a collection of algorithms that allow for the prediction of compound-induced molecular mechanisms of action by integrating phenotypic endpoints of different species; and proposes follow-up tests for model organisms to validate these pathway predictions. The Xpaths algorithms are applied to predict developmental and reproductive toxicity (DART) and implemented into an in silico platform, called DARTpaths. AVAILABILITY AND IMPLEMENTATION: All code is available on GitHub https://github.com/Xpaths/dartpaths-app under Apache license 2.0, detailed overview with demo is available at https://www.vivaltes.com/dartpaths/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Combination of pre-adapted bacteriophage therapy and antibiotics for treatment of fracture-related infection due to pandrug-resistant Klebsiella pneumoniae.

A 30-year-old bombing victim with a fracture-related pandrug-resistant Klebsiella pneumoniae infection after long-term (>700 days) antibiotic therapy is treated with a pre-adapted bacteriophage along with meropenem and colistin, followed by ceftazidime/avibactam. This salvage therapy results in objective clinical, microbiological and radiological improvement of the patient's wounds and overall condition. In support, the bacteriophage and antibiotic combination is highly effective against the patient's K. pneumoniae strain in vitro, in 7-day mature biofilms and in suspensions.

Genome-driven elucidation of phage-host interplay and impact of phage resistance evolution on bacterial fitness.

When considering the interactions between bacteriophages and their host, the issue of phage-resistance emergence is a key element in understanding the ecological impact of phages on the bacterial population. It is also an essential parameter for the implementation of phage therapy to combat antibiotic-resistant pathogens. This study investigates the phenotypic and genetic responses of five Pseudomonas aeruginosa strains (PAO1, A5803, AA43, CHA, and PAK) to the infection by seven phages with distinct evolutionary backgrounds and recognised receptors (LPS/T4P). Emerging phage-insensitivity was generally accompanied by self and cross-resistance mechanisms. Significant differences were observed between the reference PAO1 responses compared to other clinical representatives. LPS-dependent phage infections in clinical strains selected for mutations in the "global regulatory" and "other" genes, rather than in the LPS-synthesis clusters detected in PAO1 clones. Reduced fitness, as proxied by the growth rate, was correlated with large deletion (20-500 kbp) and phage carrier state. Multi-phage resistance was significantly correlated with a reduced growth rate but only in the PAO1 population. In addition, we observed that the presence of prophages decreased the lytic phage maintenance seemingly protecting the host against carrier state and occasional lytic phage propagation, thus preventing a significant reduction in bacterial growth rate.

Digital phagograms: predicting phage infectivity through a multilayer machine learning approach.

Machine learning has been broadly implemented to investigate biological systems. In this regard, the field of phage biology has embraced machine learning to elucidate and predict phage-host interactions, based on receptor-binding proteins, (anti-)defense systems, prophage detection, and life cycle recognition. Here, we highlight the enormous potential of integrating information from omics data with insights from systems biology to better understand phage-host interactions. We conceptualize and discuss the potential of a multilayer model that mirrors the phage infection process, integrating adsorption, bacterial pan-immune components and hijacking of the bacterial metabolism to predict phage infectivity. In the future, this model can offer insights into the underlying mechanisms of the infection process, and digital phagograms can support phage cocktail design and phage engineering.

Shopping for phages? Unpacking design rules for therapeutic phage cocktails.

In bacteriophage therapy, the combination of different phages into a single cocktail is of critical importance to overcome the narrow host range of single phage isolates. Today, the design of therapeutic cocktails is often akin to a black box and relies largely on intuition and (pre-)availability of isolates in local collections. Here we show that straightforward host range analysis can disclose design rules and we propose to apply/translate a data mining approach, historically applied in the field of marketing ('shopping cart analysis') to explore patterns in phage combinations. The technique is broadly applicable to host range datasets and can serve in combination with other molecular-based approaches to propose rationales for phage cocktail design.

Towards a reporting guideline for developmental and reproductive toxicology testing in C. elegans and other nematodes.

Implementation of reliable methodologies allowing Reduction, Refinement, and Replacement (3Rs) of animal testing is a process that takes several decades and is still not complete. Reliable methods are essential for regulatory hazard assessment of chemicals where differences in test protocol can influence the test outcomes and thus affect the confidence in the predictive value of the organisms used as an alternative for mammals. Although test guidelines are common for mammalian studies, they are scarce for non-vertebrate organisms that would allow for the 3Rs of animal testing. Here, we present a set of 30 reporting criteria as the basis for such a guideline for Developmental and Reproductive Toxicology (DART) testing in the nematode Caenorhabditis elegans. Small organisms like C. elegans are upcoming in new approach methodologies for hazard assessment; thus, reliable and robust test protocols are urgently needed. A literature assessment of the fulfilment of the reporting criteria demonstrates that although studies describe methodological details, essential information such as compound purity and lot/batch number or type of container is often not reported. The formulated set of reporting criteria for C. elegans testing can be used by (i) researchers to describe essential experimental details (ii) data scientists that aggregate information to assess data quality and include data in aggregated databases (iii) regulators to assess study data for inclusion in regulatory hazard assessment of chemicals.

Unraveling Protein Interactions between the Temperate Virus Bam35 and Its Bacillus Host Using an Integrative Yeast Two Hybrid-High Throughput Sequencing Approach.

Bacillus virus Bam35 is the model Betatectivirus and member of the family Tectiviridae, which is composed of tailless, icosahedral, and membrane-containing bacteriophages. Interest in these viruses has greatly increased in recent years as they are thought to be an evolutionary link between diverse groups of prokaryotic and eukaryotic viruses. Additionally, betatectiviruses infect bacteria of the Bacillus cereus group, which are known for their applications in industry and notorious since it contains many pathogens. Here, we present the first protein-protein interactions (PPIs) network for a tectivirus-host system by studying the Bam35-Bacillus thuringiensis model using a novel approach that integrates the traditional yeast two-hybrid system and high-throughput sequencing (Y2H-HTS). We generated and thoroughly analyzed a genomic library of Bam35's host B. thuringiensis HER1410 and screened interactions with all the viral proteins using different combinations of bait-prey couples. Initial analysis of the raw data enabled the identification of over 4000 candidate interactions, which were sequentially filtered to produce 182 high-confidence interactions that were defined as part of the core virus-host interactome. Overall, host metabolism proteins and peptidases were particularly enriched within the detected interactions, distinguishing this host-phage system from the other reported host-phage PPIs. Our approach also suggested biological roles for several Bam35 proteins of unknown function, including the membrane structural protein P25, which may be a viral hub with a role in host membrane modification during viral particle morphogenesis. This work resulted in a better understanding of the Bam35-B. thuringiensis interaction at the molecular level and holds great potential for the generalization of the Y2H-HTS approach for other virus-host models.

The complete genome of 2,6-dichlorobenzamide (BAM) degrader Aminobacter sp. MSH1 suggests a polyploid chromosome, phylogenetic reassignment, and functions of plasmids.

Aminobacter sp. MSH1 (CIP 110285) can use the pesticide dichlobenil and its recalcitrant transformation product, 2,6-dichlorobenzamide (BAM), as sole source of carbon, nitrogen, and energy. The concentration of BAM in groundwater often exceeds the threshold limit for drinking water, requiring additional treatment in drinking water treatment plants or closure of the affected abstraction wells. Biological treatment with MSH1 is considered a potential sustainable alternative to remediate BAM-contamination in drinking water production. We present the complete genome of MSH1, which was determined independently in two institutes at Aarhus University and KU Leuven. Divergences were observed between the two genomes, i.e. one of them lacked four plasmids compared to the other. Besides the circular chromosome and the two previously described plasmids involved in BAM catabolism, pBAM1 and pBAM2, the genome of MSH1 contained two megaplasmids and three smaller plasmids. The MSH1 substrain from KU Leuven showed a reduced genome lacking a megaplasmid and three smaller plasmids and was designated substrain MK1, whereas the Aarhus variant with all plasmids was designated substrain DK1. A plasmid stability experiment indicate that substrain DK1 may have a polyploid chromosome when growing in R2B medium with more chromosomes than plasmids per cell. Finally, strain MSH1 is reassigned as Aminobacter niigataensis MSH1.

Bacteriophage Therapy for Difficult-to-Treat Infections: The Implementation of a Multidisciplinary Phage Task Force (The PHAGEFORCE Study Protocol).

In times where only a few novel antibiotics are to be expected, antimicrobial resistance remains an expanding global health threat. In case of chronic infections caused by therapy-resistant pathogens, physicians have limited therapeutic options, which are often associated with detrimental consequences for the patient. This has resulted in a renewed interest in alternative strategies, such as bacteriophage (phage) therapy. However, there are still important hurdles that currently impede the more widespread implementation of phage therapy in clinical practice. First, the limited number of good-quality case series and clinical trials have failed to show the optimal application protocol in terms of route of administration, frequency of administration, treatment duration and phage titer. Second, there is limited information on the systemic effects of phage therapy. Finally, in the past, phage therapy has been applied intuitively in terms of the selection of phages and their combination as parts of phage cocktails. This has led to an enormous heterogeneity in previously published studies, resulting in a lack of reliable safety and efficacy data for phage therapy. We hereby present a study protocol that addresses these scientific hurdles using a multidisciplinary approach, bringing together the experience of clinical, pharmaceutical and molecular microbiology experts.

The Ever-Expanding Pseudomonas Genus: Description of 43 New Species and Partition of the Pseudomonas putida Group.

The genus Pseudomonas hosts an extensive genetic diversity and is one of the largest genera among Gram-negative bacteria. Type strains of Pseudomonas are well known to represent only a small fraction of this diversity and the number of available Pseudomonas genome sequences is increasing rapidly. Consequently, new Pseudomonas species are regularly reported and the number of species within the genus is constantly evolving. In this study, whole genome sequencing enabled us to define 43 new Pseudomonas species and provide an update of the Pseudomonas evolutionary and taxonomic relationships. Phylogenies based on the rpoD gene and whole genome sequences, including, respectively, 316 and 313 type strains of Pseudomonas, revealed sixteen groups of Pseudomonas and, together with the distribution of cyclic lipopeptide biosynthesis gene clusters, enabled the partitioning of the P. putida group into fifteen subgroups. Pairwise average nucleotide identities were calculated between type strains and a selection of 60 genomes of non-type strains of Pseudomonas. Forty-one strains were incorrectly assigned at the species level and among these, 19 strains were shown to represent an additional 13 new Pseudomonas species that remain to be formally classified. This work pinpoints the importance of correct taxonomic assignment and phylogenetic classification in order to perform integrative studies linking genetic diversity, lifestyle, and metabolic potential of Pseudomonas spp.

Integrated Multi-Omics Analysis of Mechanisms Underlying Yeast Ethanol Tolerance.

For yeast cells, tolerance to high levels of ethanol is vital both in their natural environment and in industrially relevant conditions. We recently genotyped experimentally evolved yeast strains adapted to high levels of ethanol and identified mutations linked to ethanol tolerance. In this study, by integrating genomic sequencing data with quantitative proteomics profiles from six evolved strains (data set identifier PXD006631) and construction of protein interaction networks, we elucidate exactly how the genotype and phenotype are related at the molecular level. Our multi-omics approach points to the rewiring of numerous metabolic pathways affected by genomic and proteomic level changes, from energy-producing and lipid pathways to differential regulation of transposons and proteins involved in cell cycle progression. One of the key differences is found in the energy-producing metabolism, where the ancestral yeast strain responds to ethanol by switching to respiration and employing the mitochondrial electron transport chain. In contrast, the ethanol-adapted strains appear to have returned back to energy production mainly via glycolysis and ethanol fermentation, as supported by genomic and proteomic level changes. This work is relevant for synthetic biology where systems need to function under stressful conditions, as well as for industry and in cancer biology, where it is important to understand how the genotype relates to the phenotype.

Molecular dynamics shows complex interplay and long-range effects of post-translational modifications in yeast protein interactions.

Post-translational modifications (PTMs) play a vital, yet often overlooked role in the living cells through modulation of protein properties, such as localization and affinity towards their interactors, thereby enabling quick adaptation to changing environmental conditions. We have previously benchmarked a computational framework for the prediction of PTMs' effects on the stability of protein-protein interactions, which has molecular dynamics simulations followed by free energy calculations at its core. In the present work, we apply this framework to publicly available data on Saccharomyces cerevisiae protein structures and PTM sites, identified in both normal and stress conditions. We predict proteome-wide effects of acetylations and phosphorylations on protein-protein interactions and find that acetylations more frequently have locally stabilizing roles in protein interactions, while the opposite is true for phosphorylations. However, the overall impact of PTMs on protein-protein interactions is more complex than a simple sum of local changes caused by the introduction of PTMs and adds to our understanding of PTM cross-talk. We further use the obtained data to calculate the conformational changes brought about by PTMs. Finally, conservation of the analyzed PTM residues in orthologues shows that some predictions for yeast proteins will be mirrored to other organisms, including human. This work, therefore, contributes to our overall understanding of the modulation of the cellular protein interaction networks in yeast and beyond.

Genomics of an endemic cystic fibrosis Burkholderia multivorans strain reveals low within-patient evolution but high between-patient diversity.

Burkholderia multivorans is a member of the Burkholderia cepacia complex (Bcc), notorious for its pathogenicity in persons with cystic fibrosis. Epidemiological surveillance suggests that patients predominantly acquire B. multivorans from environmental sources, with rare cases of patient-to-patient transmission. Here we report on the genomic analysis of thirteen isolates from an endemic B. multivorans strain infecting four cystic fibrosis patients treated in different pediatric cystic fibrosis centers in Belgium, with no evidence of cross-infection. All isolates share an identical sequence type (ST-742) but whole genome analysis shows that they exhibit peculiar patterns of genomic diversity between patients. By combining short and long reads sequencing technologies, we highlight key differences in terms of small nucleotide polymorphisms indicative of low rates of adaptive evolution within patient, and well-defined, hundred kbps-long segments of high enrichment in mutations between patients. In addition, we observed large structural genomic variations amongst the isolates which revealed different plasmid contents, active roles for transposase IS3 and IS5 in the deactivation of genes, and mobile prophage elements. Our study shows limited within-patient B. multivorans evolution and high between-patient strain diversity, indicating that an environmental microdiverse reservoir must be present for this endemic strain, in which active diversification is taking place. Furthermore, our analysis also reveals a set of 30 parallel adaptations across multiple patients, indicating that the specific genomic background of a given strain may dictate the route of adaptation within the cystic fibrosis lung.

Reliable Identification of Environmental Pseudomonas Isolates Using the rpoD Gene.

The taxonomic affiliation of Pseudomonas isolates is currently assessed by using the 16S rRNA gene, MultiLocus Sequence Analysis (MLSA), or whole genome sequencing. Therefore, microbiologists are facing an arduous choice, either using the universal marker, knowing that these affiliations could be inaccurate, or engaging in more laborious and costly approaches. The rpoD gene, like the 16S rRNA gene, is included in most MLSA procedures and has already been suggested for the rapid identification of certain groups of Pseudomonas. However, a comprehensive overview of the rpoD-based phylogenetic relationships within the Pseudomonas genus is lacking. In this study, we present the rpoD-based phylogeny of 217 type strains of Pseudomonas and defined a cutoff value of 98% nucleotide identity to differentiate strains at the species level. To validate this approach, we sequenced the rpoD of 145 environmental isolates and complemented this analysis with whole genome sequencing. The rpoD sequence allowed us to accurately assign Pseudomonas isolates to 20 known species and represents an excellent first diagnostic tool to identify new Pseudomonas species. Finally, rpoD amplicon sequencing appears as a reliable and low-cost alternative, particularly in the case of large environmental studies with hundreds or thousands of isolates.

Completed Genomic Sequence of Bacillus thuringiensis HER1410 Reveals a Cry-Containing Chromosome, Two Megaplasmids, and an Integrative Plasmidial Prophage.

Bacillus thuringiensis is the most used biopesticide in agriculture. Its entomopathogenic capacity stems from the possession of plasmid-borne insecticidal crystal genes (cry), traditionally used as discriminant taxonomic feature for that species. As such, crystal and plasmid identification are key to the characterization of this species. To date, about 600 B. thuringiensis genomes have been reported, but less than 5% have been completed, while the other draft genomes are incomplete, hindering full plasmid delineation. Here we present the complete genome of Bacillus thuringiensis HER1410, a strain closely related to B. thuringiensis entomocidus and a known host for a variety of Bacillus phages. The combination of short and long-read techniques allowed fully resolving the genome and delineation of three plasmids. This enabled the accurate detection of an unusual location of a unique cry gene, cry1Ba4, located in a genomic island near the chromosome replication origin. Two megaplasmids, pLUSID1 and pLUSID2 could be delineated: pLUSID1 (368 kb), a likely conjugative plasmid involved in virulence, and pLUSID2 (156 kb) potentially related to the sporulation process. A smaller plasmidial prophage pLUSID3, with a dual lifestyle whose integration within the chromosome causes the disruption of a flagellar key component. Finally, phylogenetic analysis placed this strain within a clade comprising members from the B. thuringiensis serovar thuringiensis and other serovars and with B. cereus s. s in agreement with the intermingled taxonomy of B. cereus sensu lato group.