Interest of bacterial pangenome analyses in clinical microbiology

Unlocking the Mycobacteroides abscessus pan-genome using computational tools: insights into evolutionary dynamics and lifestyle

Mycobacteroides abscessus is a non-tuberculous mycobacteria implicated in causing lung infections. It is difficult to control owing to resistance to antibiotics and disinfectants. This work was aimed at comprehending: the pan-genome architecture, evolutionary dynamics, and functionalities of pan-genome components linked to COGs and KEGG. Around 2802 core genes were present in each strain of the M. abscessus genome. The number of accessory genes ranged from 1615 to 2481. The open pan-genome of M. abscessus was attributed to the accessory genes underlining its adaptability in the host. Phylogenetic analysis revealed cluster-based relationships and highlighted factors shaping variability and adaptive capabilities. Transcription, metabolism, and pathogenic genes were vital for M. abscessus lifestyle. The accessory genes contributed to the diverse metabolic capability. The incidence of a significant portion of secondary metabolite biosynthesis genes provided insights for investigating their biosynthetic gene clusters. Additionally, a high proportion of xenobiotic biodegradation genes highlighted potential metabolic capabilities. In silico screening identified a potential vaccine candidate among hypothetical proteins in COGs. Functional analysis of M. abscessus pan-genome components unveiled factors associated with virulence, pathogenicity, infection establishment, persistence, and resistance. Notable amongst them were: MMPL family transporters, PE-PPE domain-containing proteins, TetR family transcriptional regulators, ABC transporters, Type-I, II, III, VII secretion proteins, DUF domain-containing proteins, cytochrome P450, VapC family toxin, virulence factor Mce family protein, type II toxin-antitoxin system. Overall, these results enhanced understanding of the metabolism, host-pathogen dynamics, pathogenic lifestyle, and adaptations. This will facilitate further investigations for combating infections and designing suitable therapies.

Unveiling intraspecific diversity and evolutionary dynamics of the foodborne pathogen Bacillus paranthracis through high-quality pan-genome analysis

Understanding the evolutionary dynamics of foodborne pathogens throughout host-associated habitats is of utmost importance. Bacterial pan-genomes, as dynamic entities, are strongly influenced by ecological lifestyles. As a phenotypically diverse species in the Bacillus cereus group, Bacillus paranthracis is recognized as an emerging foodborne pathogen and a probiotic simultaneously. This poorly understood species is a suitable study model for adaptive pan-genome evolution. In this study, we determined the biogeographic distribution, abundance, genetic diversity, and genotypic profiles of key genetic elements of B. paranthracis. Metagenomic read recruitment analyses demonstrated that B. paranthracis members are globally distributed and abundant in host-associated habitats. A high-quality pan-genome of B. paranthracis was subsequently constructed to analyze the evolutionary dynamics involved in ecological adaptation comprehensively. The open pan-genome indicated a flexible gene repertoire with extensive genetic diversity. Significant divergences in the phylogenetic relationships, functional enrichment, and degree of selective pressure between the different components demonstrated different evolutionary dynamics between the core and accessory genomes driven by ecological forces. Purifying selection and gene loss are the main signatures of evolutionary dynamics in B. paranthracis pan-genome. The plasticity of the accessory genome is characterized by horizontal gene transfer (HGT), massive gene losses, and weak purifying or positive selection, which might contribute to niche-specific adaptation. In contrast, although the core genome dominantly undergoes purifying selection, its association with HGT and positively selected mutations indicates its potential role in ecological diversification. Furthermore, host fitness-related dynamics are characterized by the loss of secondary metabolite biosynthesis gene clusters (BGCs) and CAZyme-encoding genes and the acquisition of antimicrobial resistance (AMR) and virulence genes via HGT. This study offers a case study of pan-genome evolution to investigate the ecological adaptations reflected by biogeographical characteristics, thereby advancing the understanding of intraspecific diversity and evolutionary dynamics of foodborne pathogens.

Unveiling a Comprehensive Multi-epitope Subunit Vaccine Strategy Against Salmonella subsp. enterica: Bridging Core, Subtractive Proteomics, and Immunoinformatics

Salmonella subsp. enterica (SE) presents a significant global health challenge in both developed and developing countries. Current SE vaccines have limitations, targeting specific strains and demonstrating moderate efficacy in adults, while also being unsuitable for young children and often unaffordable in regions with lower income levels where the disease is prevalent. To address these challenges, this study employed a computational approach integrating core proteomics, subtractive proteomics, and immunoinformatics to develop a universal SE vaccine and identify potential drug targets. Analysis of the core proteome of 185 SE strains revealed 1964 conserved proteins. Subtractive proteomics identified 9 proteins as potential vaccine candidates and 41 as novel drug targets. Using reverse vaccinology-based immunoinformatics, four multi-epitope-based subunit vaccine constructs (MESVCs) were designed, aiming to stimulate cytotoxic T lymphocyte, helper T lymphocyte, and linear B lymphocyte responses. These constructs underwent comprehensive evaluations for antigenicity, immunogenicity, toxicity, hydropathicity, and physicochemical properties. Predictive modeling, refinement, and validation were conducted to determine the secondary and tertiary structures of the SE-MESVCs, followed by docking studies with MHC-I, MHC-II, and TLR4 receptors. Molecular docking assessments showed favorable binding with all three receptors, with SE-MESVC-4 exhibiting the most promising binding energy. Molecular dynamics simulations confirmed the binding affinity and stability of SE-MESVC-4 with the TLR4/MD2 complex. Additionally, codon optimization and in silico cloning verified the efficient translation and successful expression of SE-MESVC-4 in Escherichia coli (E. coli) str. K12. Subsequent in silico immune simulation evaluated the efficacy of SE-MESVC-4 in triggering an effective immune response. These results suggest that SE-MESVC-4 may induce both humoral and cellular immune responses, making it a potential candidate for an effective SE vaccine. However, further experimental investigations are necessary to validate the immunogenicity and efficacy of SE-MESVC-4, bringing us closer to effectively combating SE infections.

RIBAP: a comprehensive bacterial core genome annotation pipeline for pangenome calculation beyond the species level

Microbial pangenome analysis identifies present or absent genes in prokaryotic genomes. However, current tools are limited when analyzing species with higher sequence diversity or higher taxonomic orders such as genera or families. The Roary ILP Bacterial core Annotation Pipeline (RIBAP) uses an integer linear programming approach to refine gene clusters predicted by Roary for identifying core genes. RIBAP successfully handles the complexity and diversity of Chlamydia, Klebsiella, Brucella, and Enterococcus genomes, outperforming other established and recent pangenome tools for identifying all-encompassing core genes at the genus level. RIBAP is a freely available Nextflow pipeline at github.com/hoelzer-lab/ribap and zenodo.org/doi/10.5281/zenodo.10890871.

PanDelos-frags: A methodology for discovering pangenomic content of incomplete microbial assemblies

Pangenomics was originally defined as the problem of comparing the composition of genes into gene families within a set of bacterial isolates belonging to the same species. The problem requires the calculation of sequence homology among such genes. When combined with metagenomics, namely for human microbiome composition analysis, gene-oriented pangenome detection becomes a promising method to decipher ecosystem functions and population-level evolution. Established computational tools are able to investigate the genetic content of isolates for which a complete genomic sequence is available. However, there is a plethora of incomplete genomes that are available on public resources, which only a few tools may analyze. Incomplete means that the process for reconstructing their genomic sequence is not complete, and only fragments of their sequence are currently available. However, the information contained in these fragments may play an essential role in the analyses. Here, we present PanDelos-frags, a computational tool which exploits and extends previous results in analyzing complete genomes. It provides a new methodology for inferring missing genetic information and thus for managing incomplete genomes. PanDelos-frags outperforms state-of-the-art approaches in reconstructing gene families in synthetic benchmarks and in a real use case of metagenomics. PanDelos-frags is publicly available at https://github.com/InfOmics/PanDelos-frags.

Diversity and prevalence of type VI secretion system effectors in clinical Pseudomonas aeruginosa isolates

Pseudomonas aeruginosa is an opportunistic pathogen and a major driver of morbidity and mortality in people with Cystic Fibrosis (CF). The Type VI secretion system (T6SS) is a molecular nanomachine that translocates effectors across the bacterial membrane into target cells or the extracellular environment enabling intermicrobial interaction. P. aeruginosa encodes three T6SS clusters, the H1-, H2- and H3-T6SS, and numerous orphan islands. Genetic diversity of T6SS-associated effectors in P. aeruginosa has been noted in reference strains but has yet to be explored in clinical isolates. Here, we perform a comprehensive bioinformatic analysis of the pangenome and T6SS effector genes in 52 high-quality clinical P. aeruginosa genomes isolated from CF patients and housed in the Personalised Approach to P. aeruginosa strain repository. We confirm that the clinical CF isolate pangenome is open and principally made up of accessory and unique genes that may provide strain-specific advantages. We observed genetic variability in some effector/immunity encoding genes and show that several well-characterised vgrG and PAAR islands are absent from numerous isolates. Our analysis shows clear evidence of disruption to T6SS genomic loci through transposon, prophage, and mobile genetic element insertions. We identified an orphan vgrG island in P. aeruginosa strain PAK and five clinical isolates using in silico analysis which we denote vgrG7, predicting a gene within this cluster to encode a Tle2 lipase family effector. Close comparison of T6SS loci in clinical isolates compared to reference P. aeruginosa strain PAO1 revealed the presence of genes encoding eight new T6SS effectors with the following putative functions: cytidine deaminase, lipase, metallopeptidase, NADase, and pyocin. Finally, the prevalence of characterised and putative T6SS effectors were assessed in 532 publicly available P. aeruginosa genomes, which suggests the existence of accessory effectors. Our in silico study of the P. aeruginosa T6SS exposes a level of genetic diversity at T6SS genomic loci not seen to date within P. aeruginosa, particularly in CF isolates. As understanding the effector repertoire is key to identifying the targets of T6SSs and its efficacy, this comprehensive analysis provides a path for future experimental characterisation of these mediators of intermicrobial competition and host manipulation.

Recent Advances in Genomics-Based Approaches for the Development of Intracellular Bacterial Pathogen Vaccines

Infectious diseases continue to be a leading cause of morbidity and mortality worldwide. The majority of infectious diseases are caused by intracellular pathogenic bacteria (IPB). Historically, conventional vaccination drives have helped control the pathogenesis of intracellular bacteria and the emergence of antimicrobial resistance, saving millions of lives. However, in light of various limitations, many diseases that involve IPB still do not have adequate vaccines. In response to increasing demand for novel vaccine development strategies, a new area of vaccine research emerged following the advent of genomics technology, which changed the paradigm of vaccine development by utilizing the complete genomic data of microorganisms against them. It became possible to identify genes related to disease virulence, genetic patterns linked to disease virulence, as well as the genetic components that supported immunity and favorable vaccine responses. Complete genomic databases, and advancements in transcriptomics, metabolomics, structural genomics, proteomics, immunomics, pan-genomics, synthetic genomics, and population biology have allowed researchers to identify potential vaccine candidates and predict their effects in patients. New vaccines have been created against diseases for which previously there were no vaccines available, and existing vaccines have been improved. This review highlights the key issues and explores the evolution of vaccines. The increasing volume of IPB genomic data, and their application in novel genome-based techniques for vaccine development, were also examined, along with their characteristics, and the opportunities and obstacles involved. Critically, the application of genomics technology has helped researchers rapidly select and evaluate candidate antigens. Novel vaccines capable of addressing the limitations associated with conventional vaccines have been developed and pressing healthcare issues are being addressed.

Critical assessment of pan-genomic analysis of metagenome-assembled genomes

Pan-genome analyses of metagenome-assembled genomes (MAGs) may suffer from the known issues with MAGs: fragmentation, incompleteness and contamination. Here, we conducted a critical assessment of pan-genomics of MAGs, by comparing pan-genome analysis results of complete bacterial genomes and simulated MAGs. We found that incompleteness led to significant core gene (CG) loss. The CG loss remained when using different pan-genome analysis tools (Roary, BPGA, Anvi'o) and when using a mixture of MAGs and complete genomes. Contamination had little effect on core genome size (except for Roary due to in its gene clustering issue) but had major influence on accessory genomes. Importantly, the CG loss was partially alleviated by lowering the CG threshold and using gene prediction algorithms that consider fragmented genes, but to a less degree when incompleteness was higher than 5%. The CG loss also led to incorrect pan-genome functional predictions and inaccurate phylogenetic trees. Our main findings were supported by a study of real MAG-isolate genome data. We conclude that lowering CG threshold and predicting genes in metagenome mode (as Anvi'o does with Prodigal) are necessary in pan-genome analysis of MAGs. Development of new pan-genome analysis tools specifically for MAGs are needed in future studies.

The potential impact of a probiotic: Akkermansia muciniphila in the regulation of blood pressure—the current facts and evidence

Akkermansia muciniphila (A. muciniphila) is present in the human gut microbiota from infancy and gradually increases in adulthood. The potential impact of the abundance of A. muciniphila has been studied in major cardiovascular diseases including elevated blood pressure or hypertension (HTN). HTN is a major factor in premature death worldwide, and approximately 1.28 billion adults aged 30–79 years have hypertension. A. muciniphila is being considered a next-generation probiotic and though numerous studies had highlighted the positive role of A. muciniphila in lowering/controlling the HTN, however, few studies had highlighted the negative impact of increased abundance of A. muciniphila in the management of HTN. Thus, in the review, we aimed to discuss the current facts, evidence, and controversy about the role of A. muciniphila in the pathophysiology of HTN and its potential effect on HTN management/regulation, which could be beneficial in identifying the drug target for the management of HTN.

Whole-genome analyses reveal a novel prophage and cgSNPs-derived sublineages of Brachyspira hyodysenteriae ST196

Background

Brachyspira (B.) hyodysenteriae is a fastidious anaerobe spirochete that can cause swine dysentery, a severe mucohaemorragic colitis that affects pig production and animal welfare worldwide. In Switzerland, the population of B. hyodysenteriae is characterized by the predominance of macrolide-lincosamide-resistant B. hyodysenteriae isolates of sequence type (ST) ST196, prompting us to obtain deeper insights into the genomic structure and variability of ST196 using pangenome and whole genome variant analyses.

Results

The draft genome of 14 B. hyodysenteriae isolates of ST196, sampled during a 7-year period from geographically distant pig herds, was obtained by whole-genome sequencing (WGS) and compared to the complete genome of the B. hyodysenteriae isolate Bh743-7 of ST196 used as reference. Variability results revealed the existence of 30 to 52 single nucleotide polymorphisms (SNPs), resulting in eight sublineages of ST196. The pangenome analysis led to the identification of a novel prophage, pphBhCH20, of the Siphoviridae family in a single isolate of ST196, which suggests that horizontal gene transfer events may drive changes in genomic structure.

Conclusions

This study contributes to the catalogue of publicly available genomes and provides relevant bioinformatic tools and information for further comparative genomic analyses for B. hyodysenteriae. It reveals that Swiss B. hyodysenteriae isolates of the same ST may have evolved independently over time by point mutations and acquisition of larger genetic elements. In line with this, the third type of mobile genetic element described so far in B. hyodysenteriae, the novel prophage pphBhCH20, has been identified in a single isolate of B. hyodysenteriae of ST196.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08347-5.