Sequencing-Based Genotyping of Pakistani Burkholderia mallei Strains: A Useful Way for Investigating Glanders Outbreaks

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.

Tracking the footsteps of Burkholderia mallei: determination of the molecular differences and potential resistance genes

Background/aim

Chemical biological radiological nuclear threats are at an important point in the agenda of world health today, as they can cause mass deaths. B. mallei attracts attention as a potential biological warfare agent due to its features such as multidrug resistance, a rapid transmission mechanism via aerosol, the absence of a complete treatment protocol for the infection it causes, and the absence of an approved vaccine for protection against the bacteria. B. mallei suspect samples must be studied by experienced personnel in biosafety level III laboratories. B mallei is a difficult and troublesome pathogen to diagnose and many unknowns about B. mallei today. Therefore, the aim of the study was to determine the molecular differences and potential resistance genes of B mallei strains.

Materials and methods

Determination of the molecular differences and potential resistance genes of B mallei strains with new bioinformatics approaches by comparatively examining the data of 29 B mallei strains, 10 of which were isolated from Türkiye, on the genome list of the National Biotechnology Information Center (NCBI).

Results

According to the genome annotations of the origins, the origin containing the highest number of CDS which is 5172 was found as the 11th strain obtained in Türkiye in 1949. The origin with the highest number of pseudogenes was determined as 23,344 (China 7) origin. Two hundred and eighty-five pseudogenes found in this strain were obtained from a knee effusion in Myanmar. According to chromosome 2 data, B. mallei strain was determined as the most similar strain to ATCC 23344, line 11 with NCTC 10229 strain, and SAVP1 strain was determined as the least similar strain. When the antimicrobial resistance gene markers of the isolates included in the study were examined, amrA and amrB, qacG ade, Burkholderia pseudomallei Omp38 were found to be carrying.

Conclusion

In terms of public health, it was thought that the data obtained as a result of our study about B mallei, which is defined as a biological weapon, is very valuable for creating treatment protocols to be applied to possible epidemics in the future. In addition, the available genetic epidemiological data of these strains belonging to a category that is dangerous to work with in a laboratory environment were reviewed.

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.

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.

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.

Comparison of Illumina and Oxford Nanopore Technology for genome analysis of Francisella tularensis, Bacillus anthracis, and Brucella suis

BACKGROUND

Bacterial epidemiology needs to understand the spread and dissemination of strains in a One Health context. This is important for highly pathogenic bacteria such as Bacillus anthracis, Brucella species, and Francisella tularensis. Whole genome sequencing (WGS) has paved the way for genetic marker detection and high-resolution genotyping. While such tasks are established for Illumina short-read sequencing, Oxford Nanopore Technology (ONT) long-read sequencing has yet to be evaluated for such highly pathogenic bacteria with little genomic variations between strains. In this study, three independent sequencing runs were performed using Illumina, ONT flow cell version 9.4.1, and 10.4 for six strains of each of Ba. anthracis, Br. suis and F. tularensis. Data from ONT sequencing alone, Illumina sequencing alone and two hybrid assembly approaches were compared.

RESULTS

As previously shown, ONT produces ultra-long reads, while Illumina produces short reads with higher sequencing accuracy. Flow cell version 10.4 improved sequencing accuracy over version 9.4.1. The correct (sub-)species were inferred from all tested technologies, individually. Moreover, the sets of genetic markers for virulence, were almost identical for the respective species. The long reads of ONT allowed to assemble not only chromosomes of all species to near closure, but also virulence plasmids of Ba. anthracis. Assemblies based on nanopore data alone, Illumina data alone, and both hybrid assemblies correctly detected canonical (sub-)clades for Ba. anthracis and F. tularensis as well as multilocus sequence types for Br. suis. For F. tularensis, high-resolution genotyping using core-genome MLST (cgMLST) and core-genome Single-Nucleotide-Polymorphism (cgSNP) typing produced highly comparable results between data from Illumina and both ONT flow cell versions. For Ba. anthracis, only data from flow cell version 10.4 produced similar results to Illumina for both high-resolution typing methods. However, for Br. suis, high-resolution genotyping yielded larger differences comparing Illumina data to data from both ONT flow cell versions.

CONCLUSIONS

In summary, combining data from ONT and Illumina for high-resolution genotyping might be feasible for F. tularensis and Ba. anthracis, but not yet for Br. suis. The ongoing improvement of nanopore technology and subsequent data analysis may facilitate high-resolution genotyping for all bacteria with highly stable genomes in future.

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.