Genetic diversity and spatial distribution of Burkholderia mallei by core genome-based multilocus sequence typing analysis

Virulome and phylogenomic profiling of a novel Burkholderia pseudomallei strain from an Indian clinical isolate

Highly pathogenic Burkholderia pseudomallei is the causative agent of melioidosis, a neglected tropical disease endemic in Southeast Asian tropical region. This bacterium encompasses diverse virulence factors which further undergo dynamic gene-expression flux as it transits through distinct environmental niches within the host which may lead to manifestation of differential clinical symptoms. B. pseudomallei, is classified as a Tier 1 select agent in the United States and regarded as a risk group 3 organism in India with the potential to be used as bioweapon. Considering these facts, it is vital to uncover both physiological and genetic heterogeneity of B. pseudomallei, particularly to identify any novel virulence factors that may contribute to pathogenicity. B. pseudomallei strain CM000113 was isolated from a clinical case in India, characterized it for its physiological, biochemical, and prominently genetic traits through WGS. It has a type 2 morphotype with faster doubling time and high biofilm producing capacity as compared to Pseudomonas aeruginosa. The genome size is 7.3 Mbp and it is phylogenetically close to B. pseudomallei strain Mahidol 1106a and Burkholderia mallei Turkey 2. We observed genetic heterogeneity, as key virulence factors that were identified shows sequence dissimilarity with reference strains. Additionally, presence of genomic islands, harbouring two virulence factors, GmhA and GmhB2, associated with pathogenesis indicates possibility of horizontal gene transfer. These results emphasize the need for an extensive study focusing the genome of B. pseudomallei and its associated heterogeneity, to identify molecular biomarkers aiding to develop point-of-care diagnostic kits for early diagnosis of melioidosis.

BakRep - a searchable large-scale web repository for bacterial genomes, characterizations and metadata

Bacteria are fascinating research objects in many disciplines for countless reasons, and whole-genome sequencing (WGS) has become the paramount methodology to advance our microbiological understanding. Meanwhile, access to cost-effective sequencing platforms has accelerated bacterial WGS to unprecedented levels, introducing new challenges in terms of data accessibility, computational demands, heterogeneity of analysis workflows and, thus, ultimately its scientific usability. To this end, a previous study released a uniformly processed set of 661 405 bacterial genome assemblies obtained from the European Nucleotide Archive as of November 2018. Building on these accomplishments, we conducted further genome-based analyses like taxonomic classification, multilocus sequence typing and annotation of all genomes. Here, we present BakRep, a searchable large-scale web repository of these genomes enriched with consistent genome characterizations and original metadata. The platform provides a flexible search engine combining taxonomic, genomic and metadata information, as well as interactive elements to visualize genomic features. Furthermore, all results can be downloaded for offline analyses via an accompanying command line tool. The web repository is accessible via https://bakrep.computational.bio.

First molecular characterization of Burkholderia mallei strains isolated from horses in Mongolia

Glanders, a highly contagious and often fatal disease affecting equids, is caused by Burkholderia mallei. Although sporadic cases of equine glanders have recently been documented in Mongolia, genome sequencing and molecular studies of the bacteria within this region are lacking. This study provided the first molecular characterization of B. mallei isolated from four native Mongolian horses from two different provinces in 2019 and 2022 by applying whole-genome sequencing with two SNP types (previously developed genotyping with 15 SNP markers that provide global coverage of the B. mallei population and the core genome coding SNP typing developed in this study). The Mongolian isolates were located within the L3B1 cluster, which was previously associated with the V-120 strain from Russia. Within the L3B1 cluster shared by neighboring countries, they were in a unique subbranch. In this study, specific SNP markers unique to the Mongolian strains were identified to track these strains using a high-resolution melting analysis (HRMA). This study revealed the unique phylogenetic background of Mongolian strains isolated from the eastern part of Mongolia. HRMA specific to the Mongolian subbranch may contribute to the molecular epidemiological monitoring of glanders in Mongolia and surrounding countries.

Development and implementation of a core genome multilocus sequence typing scheme for Haemophilus influenzae

Haemophilus influenzae is part of the human nasopharyngeal microbiota and a pathogen causing invasive disease. The extensive genetic diversity observed in H. influenzae necessitates discriminatory analytical approaches to evaluate its population structure. This study developed a core genome multilocus sequence typing (cgMLST) scheme for H. influenzae using pangenome analysis tools and validated the cgMLST scheme using datasets consisting of complete reference genomes (N = 14) and high-quality draft H. influenzae genomes (N = 2297). The draft genome dataset was divided into a development dataset (N = 921) and a validation dataset (N = 1376). The development dataset was used to identify potential core genes, and the validation dataset was used to refine the final core gene list to ensure the reliability of the proposed cgMLST scheme. Functional classifications were made for all the resulting core genes. Phylogenetic analyses were performed using both allelic profiles and nucleotide sequence alignments of the core genome to test congruence, as assessed by Spearman's correlation and ordinary least square linear regression tests. Preliminary analyses using the development dataset identified 1067 core genes, which were refined to 1037 with the validation dataset. More than 70% of core genes were predicted to encode proteins essential for metabolism or genetic information processing. Phylogenetic and statistical analyses indicated that the core genome allelic profile accurately represented phylogenetic relatedness among the isolates (R 2 = 0.945). We used this cgMLST scheme to define a high-resolution population structure for H. influenzae, which enhances the genomic analysis of this clinically relevant human pathogen.

Ralstonia pickettii bloodstream infections with potential genomic link to internationally distributed contaminated saline solution, Germany, October 2023

Ralstonia pickettii is a Gram-negative rod which may cause invasive infections when they contaminate liquid medical products. After R. pickettii was detected in blood cultures and a stem cell product from three patients in a tertiary care hospital in Germany, whole genome sequencing of these three isolates and two water isolates from the environment was performed. Core genome multilocus sequence typing analysis showed that the three patient isolates were closely related and there was a large distance to the environmental isolates. In a genomic comparison, the patients' isolates were distantly related to an R. pickettii strain from a cluster in Australia suspected to be caused by contaminated saline produced in India, while all liquid medical products with a link to all patients were produced in Europe or the United States. Our data point towards an ongoing risk by an unknown common source that could be traced back to medical products contaminated with R. pickettii and potentially distributed worldwide. Investigating invasive R. pickettii infections, identifying and testing medical products administered to the patients and timely whole genome sequencing may help identify the exact source of this potentially global outbreak.

Influence of genomic variations on glanders serodiagnostic antigens using integrative genomic and transcriptomic approaches

Glanders is a highly contagious and life-threatening zoonotic disease caused by Burkholderia mallei (B. mallei). Without an effective vaccine or treatment, early diagnosis has been regarded as the most effective method to prevent glanders transmission. Currently, the diagnosis of glanders is heavily reliant on serological tests. However, given that markedly different host immune responses can be elicited by genetically different strains of the same bacterial species, infection by B. mallei, whose genome is unstable and plastic, may result in various immune responses. This variability can make the serodiagnosis of glanders challenging. Therefore, there is a need for a comprehensive understanding and assessment of how B. mallei genomic variations impact the appropriateness of specific target antigens for glanders serodiagnosis. In this study, we investigated how genomic variations in the B. mallei genome affect gene content (gene presence/absence) and expression, with a special focus on antigens used or potentially used in serodiagnosis. In all the genome sequences of B. mallei isolates available in NCBI's RefSeq database (accessed in July 2023) and in-house sequenced samples, extensive small and large variations were observed when compared to the type strain ATCC 23344. Further pan-genome analysis of those assemblies revealed variations of gene content among all available genomes of B. mallei. Specifically, differences in gene content ranging from 31 to 715 genes with an average of 334 gene presence-absence variations were found in strains with complete or chromosome-level genome assemblies, using the ATCC 23344 strain as a reference. The affected genes included some encoded proteins used as serodiagnostic antigens, which were lost due mainly to structural variations. Additionally, a transcriptomic analysis was performed using the type strain ATCC 23344 and strain Zagreb which has been widely utilized to produce glanders antigens. In total, 388 significant differentially expressed genes were identified between these two strains, including genes related to bacterial pathogenesis and virulence, some of which were associated with genomic variations, particularly structural variations. To our knowledge, this is the first comprehensive study to uncover the impacts of genetic variations of B. mallei on its gene content and expression. These differences would have significant impacts on host innate and adaptive immunity, including antibody production, during infection. This study provides novel insights into B. mallei genetic variants, knowledge which will help to improve glanders serodiagnosis.

Virulence and resistance patterns of Vibrio cholerae non-O1/non-O139 acquired in Germany and other European countries

Global warming has caused an increase in the emergence of Vibrio species in marine and estuarine environments as well as fresh water bodies. Over the past decades, antimicrobial resistance (AMR) has evolved among Vibrio species toward various antibiotics commonly used for the treatment of Vibrio infections. In this study, we assessed virulence and resistance patterns of Vibrio cholerae non-O1/non-O139 strains derived from Germany and other European countries. A total of 63 clinical and 24 environmental Vibrio cholerae non-O1/non-O139 strains, collected between 2011 and 2021, were analyzed. In silico antibiotic resistances were compared with resistance phenotypes according to EUCAST breakpoints. Additionally, genetic relatedness between isolates was assessed by two cgMLST schemes (SeqSphere +, pubMLST). Both cgMLST schemes yielded similar results, indicating high genetic diversity among V. cholerae non-O1/non-O139 isolates. Some isolates were found to be genetically closely related (allelic distance < 20), which suggests an epidemiological link. Thirty-seven virulence genes (VGs) were identified among 87 V. cholerae non-O1/non-O139 isolates, which resulted in 38 virulence profiles (VPs). VPs were similar between clinical and environmental isolates, with the exception of one clinical isolate that displayed a higher abundance of VGs. Also, a cluster of 11 environmental isolates was identified to have the lowest number of VGs. Among all strains, the predominant virulence factors were quorum sensing protein (luxS), repeats-in-toxins (rtxC/rtxD), hemolysin (hlyA) and different type VI secretion systems (T6SS) genes. The genotypic profiles revealed antibiotic resistance genes (ARGs) associated with resistance to beta-lactams, quinolones, macrolides, tetracycline, antifolate, aminoglycosides, fosfomycin, phenicols and sulfonamide. Carbapenemase gene VCC-1 was detected in 10 meropenem-resistant V. cholerae non-O1/non-O139 isolates derived from surface water in Germany. The proportion of resistance among V. cholerae non-O1/non-O139 species isolates against first line treatment (3rd generation cephalosporin, tetracycline and fluoroquinolone) was low. Empirical treatment would likely have been effective for all of the clinical V. cholerae non-O1/non-O139 isolates examined. Nevertheless, carbapenem-resistant isolates have been present in fresh water in Germany and might represent a reservoir for ARGs. Monitoring antimicrobial resistance is crucial for public health authorities to minimize the risks for the human population.

Assessment of the Virulence of the Burkholderia mallei Strain BAC 86/19 in BALB/c Mice

Burkholderia mallei is an aerobic, Gram-negative, non-motile bacillus. As an obligate mammalian pathogen, it primarily affects solipeds. Although rarely transmitted to humans, the disease it causes, glanders, is classified as a zoonosis. The bacterium was officially eradicated in Brazil in 1969; however, it reemerged after three decades. This study aims to assess the virulence of a specific B. mallei strain, isolated in Brazil, in BALB/c mice through intranasal infection. The strain, B. mallei BAC 86/19, was obtained from the tracheal secretion of a young mare displaying positive serology but no clinical signs of glanders. Post-mortem examinations revealed macroscopic lesions consistent with the disease, however. In mice, the LD50 was determined to be approximately 1.59 × 105 colony-forming units (CFU)/animal. Mice exposed to either 0.1 × LD50 or 1 × LD50 displayed transient weight loss, which resolved after three or five days, respectively. B. mallei persisted within the liver and lung for five days post-infection and in the spleen for seven days. These findings underscore the detectable virulence of the Brazilian B. mallei BAC 86/19 strain in mice, which are relatively resilient hosts. This research points to the importance of the continued investigation of the virulence mechanisms and potential countermeasures associated with B. mallei infections, including their Brazilian isolates.

Glanders Diagnosis in an Asymptomatic Mare from Brazil: Insights from Serology, Microbiological Culture, Mass Spectrometry, and Genome Sequencing

This manuscript elucidates the occurrence of glanders in an asymptomatic mare from Brazil presenting positive Burkholderia mallei antibody titers. The diagnosis was established through a multi-pronged approach encompassing microbiological culture, mass spectrometry, and genome sequencing. The outbreak occurred in 2019 in Tatuí, São Paulo, Brazil, and the infected mare, despite displaying no clinical symptoms, had multiple miliary lesions in the liver, as well as intense catarrhal discharge in the trachea. Samples were collected from various organs and subjected to bacterial isolation, molecular detection, and identification. The strain was identified as B. mallei using PCR and confirmed by MALDI-TOF mass spectrometry. Whole-genome sequencing revealed a genome size of 5.51 Mb with a GC content of 65.8%, 5871 genes (including 4 rRNA and 53 tRNA genes), and 5583 coding DNA sequences (CDSs). Additionally, 227 predicted pseudogenes were detected. In silico analysis of different genomic loci that allow for differentiation with Burkholderia pseudomallei confirmed the identity of the isolate as B. mallei, in addition to the characteristic genome size. The BAC 86/19 strain was identified as lineage 3, sublineage 2, which includes other strains from Brazil, India, and Iran. The genome sequencing of this strain provides valuable information that can be used to better understand the pathogen and its epidemiology, as well as to develop diagnostic tools for glanders.

Molecular detection of Burkholderia mallei in different geographic regions of Brazil

Glanders is a contagious disease of equids caused by the Gram-negative bacterium Burkholderia mallei. In Brazil, the disease is considered to be reemerging and has been expanding, with records of equids with positive serology in most of the federative units. However, there are few reports describing the genotypic detection of the agent. This study demonstrated the detection of B. mallei by species-specific PCR directly from tissues or from bacterial cultures, followed by amplicon sequencing in equids (equines, mules, and asinines) with positive serology for glanders in all five geographic regions of Brazil. The molecular evidence of B. mallei infection in serologically positive equids in this study expands the possibility of strain isolation and the conduction of epidemiological characterizations based on molecular information. The microbiological detection of B. mallei in cultures from nasal and palate swabs, even in equids without clinical manifestations, raises the possibility of environmental elimination of the agent.

Expression of virulence and antimicrobial related proteins in Burkholderia mallei and Burkholderia pseudomallei

BACKGROUND

Burkholderia mallei and Burkholderia pseudomallei are both potential biological threat agents. Melioidosis caused by B. pseudomallei is endemic in Southeast Asia and Northern Australia, while glanders caused by B. mallei infections are rare. Here we studied the proteomes of different B. mallei and B. pseudomallei isolates to determine species specific characteristics.

METHODS

The expressed proteins of 5 B. mallei and 6 B. pseudomallei strains were characterized using liquid chromatography high-resolution tandem mass spectrometry (LC-HRMS/MS). Subsequently, expression of potential resistance and virulence related characteristics were analyzed and compared.

RESULTS

Proteome analysis can be used for the identification of B. mallei and B. pseudomallei. Both species were identified based on >60 discriminative peptides. Expression of proteins potentially involved in antimicrobial resistance, AmrAB-OprA, BpeAB-OprB, BpeEF-OprC, PenA as well as several other efflux pump related proteins and putative β-lactamases was demonstrated. Despite, the fact that efflux pump BpeAB-OprB was expressed in all isolates, no clear correlation with an antimicrobial phenotype and the efflux-pump could be established. Also consistent with the phenotypes, no amino acid mutations in PenA known to result in β-lactam resistance could be identified. In all studied isolates, the expression of virulence (related) factors Capsule-1 and T2SS was demonstrated. The expression of T6SS-1 was demonstrated in all 6 B. pseudomallei isolates and in 2 of the 5 B. mallei isolates. In all, except one B. pseudomallei isolate, poly-beta-1,6 N-acetyl-D-glucosamine export porin (Pga), important for biofilm formation, was detected, which were absent in the proteomes of B. mallei. Siderophores, iron binding proteins, malleobactin and malleilactone are possibly expressed in both species under standard laboratory growth conditions. Expression of multiple proteins from both the malleobactin and malleilactone polyketide synthase (PKS) and non-ribosomal peptide synthetase (NRPS) clusters was demonstrated in both species. All B. pseudomallei expressed at least seven of the nine proteins of the bactobolin synthase cluster (bactobolin, is a ribosome targeting antibiotic), while only in one B. mallei isolate expression of two proteins of this synthase cluster was identified.

CONCLUSIONS

Analyzing the expressed proteomes revealed differences between B. mallei and B. pseudomallei but also between isolates from the same species. Proteome analysis can be used not only to identify B. mallei and B. pseudomallei but also to characterize the presence of important factors that putatively contribute to the pathogenesis of B. mallei and B. pseudomallei.

Sequence-based detection and typing procedures for Burkholderia mallei: Assessment and prospects

Although glanders has been eradicated in most of the developed world, the disease still persists in various countries such as Brazil, India, Pakistan, Bangladesh, Nepal, Iran, Bahrain, UAE and Turkey. It is one of the notifiable diseases listed by the World Organization for Animal Health. Occurrence of glanders imposes restriction on equestrian events and restricts equine movement, thus causing economic losses to equine industry. The genetic diversity and global distribution of the causing agent, Burkholderia (B.) mallei, have not been assessed in detail and are complicated by the high clonality of this organism. Among the identification and typing methods, PCR-based methods for distinguishing B. mallei from its close relative B. pseudomallei as well as genotyping using tandem repeat regions (MLVA) are established. The advent and continuous advancement of the sequencing techniques and the reconstruction of closed genomes enable the development of genome guided epidemiological tools. For achieving a higher genomic resolution, genotyping methods based on whole genome sequencing data can be employed, like genome-wide single nucleotide polymorphisms. One of the limitations in obtaining complete genomic sequences for further molecular characterization of B. mallei is its high GC content. In this review, we aim to provide an overview of the widely used detection and typing methods for B. mallei and illustrate gaps that still require development. The genomic features of Burkholderia, their high homology and clonality will be first described from a comparative genomics perspective. Then, the commonly used molecular detection (PCR systems) and typing systems (e.g., multilocus sequence typing, variable number of tandem repeat analysis) will be presented and put in perspective with recently developed genomic methods. Also, the increasing availability of B. mallei genomic sequences and evolution of the sequencing methods offers exciting prospects for further refinement of B. mallei typing, that could overcome the difficulties presently encountered with this particular bacterium.