Continuing evolution of Burkholderia mallei through genome reduction and large-scale rearrangements

Analysis of In-Patient Evolution of Escherichia coli Reveals Potential Links to Relapse of Bone and Joint Infections

Bone and joint infections (BJIs) are difficult to treat and affect a growing number of patients, in which relapses are observed in 10-20% of case. These relapses, which call for prolonged antibiotic treatment and increase resistance emergence risk, may originate from ill-understood adaptation of the pathogen to the host. Here, we investigated 3 pairs of Escherichia coli strains from BJI cases and their relapses to unravel adaptations within patients. Whole-genome comparison presented evidence for positive selection and phenotypic characterization showed that biofilm formation remained unchanged, contrary to what is usually described in such cases. Although virulence was not modified, we identified the loss of 2 virulence factors contributing to immune system evasion in one of the studied strains. Other strategies, including global growth optimization and colicin production, likely allowed the strains to outcompete competitors. This work highlights the variety of strategies allowing in-patient adaptation in BJIs.

Isolation of Burkholderia pseudomallei from a goat in New Caledonia: implications for animal and human health monitoring and serological tool comparison

BACKGROUND

Melioidosis is a serious bacterial infection caused by Burkholderia pseudomallei, a gram-negative bacterium commonly found in soil and water. It can affect both humans and animals, and is endemic in regions such as Southeast Asia and Northern Australia. In recent years, there have been reports of an emergence of human melioidosis in other areas, including New Caledonia.

RESULTS

During standard laboratory analysis in New Caledonia in 2021, a strain of B. pseudomallei was isolated from a goat. The strain was characterized using both MLST and WGS techniques and was found to cluster with previously described local human strains from the area. In parallel, several serological tests (CFT, ELISA, Luminex (Hcp1, GroEL, BPSS1840), arrays assay and a latex agglutination test) were performed on animals from the farm where the goat originated, and/or from three other neighboring farms. Using two commercial ELISA kits, seropositive animals were found only on the farm where the infected goat originated and tests based on recombinant proteins confirmed the usefulness of the Hcp1 protein for the diagnosis of melioidosis in animals.

CONCLUSIONS

Despite the regular reports of human cases, this is the first confirmed case of melioidosis in an animal in New Caledonia. These results confirm the presence of the bacterium in the region and highlight the importance of vigilance for both animal and human health. It is critical that all health partners, including breeders, veterinarians, and biologists, work together to monitor and prevent the spread of the disease.

Burkholderia pseudomallei Complex Subunit and Glycoconjugate Vaccines and Their Potential to Elicit Cross-Protection to Burkholderia cepacia Complex

Burkholderia are a group of Gram-negative bacteria that can cause a variety of diseases in at-risk populations. B. pseudomallei and B. mallei, the etiological agents of melioidosis and glanders, respectively, are the two clinically relevant members of the B. pseudomallei complex (Bpc). The development of vaccines against Bpc species has been accelerated in recent years, resulting in numerous promising subunits and glycoconjugate vaccines incorporating a variety of antigens. However, a second group of pathogenic Burkholderia species exists known as the Burkholderia cepacia complex (Bcc), a group of opportunistic bacteria which tend to affect individuals with weakened immunity or cystic fibrosis. To date, there have been few attempts to develop vaccines to Bcc species. Therefore, the primary goal of this review is to provide a broad overview of the various subunit antigens that have been tested in Bpc species, their protective efficacy, study limitations, and known or suspected mechanisms of protection. Then, we assess the reviewed Bpc antigens for their amino acid sequence conservation to homologous proteins found in Bcc species. We propose that protective Bpc antigens with a high degree of Bpc-to-Bcc sequence conservation could serve as components of a pan-Burkholderia vaccine capable of protecting against both disease-causing groups.

Live tracking of a plant pathogen outbreak reveals rapid and successive, multidecade plasmid reduction

Quickly understanding the genomic changes that lead to pathogen emergence is necessary to launch mitigation efforts and reduce harm. In this study, we tracked in real time a 2022 bacterial plant disease outbreak in U.S. geraniums (Pelargonium × hortorum) caused by Xhp2022, a novel lineage of Xanthomonas hortorum. Genomes from 31 Xhp2022 isolates from seven states showed limited chromosomal variation and all contained a single plasmid (p93). Time tree and single nucleotide polymorphism whole-genome analysis estimated that Xhp2022 emerged within the last decade. The phylogenomic analysis determined that p93 resulted from the cointegration of three plasmids (p31, p45, and p66) found sporadically across isolates from previous outbreaks. Although p93 had a 49 kb nucleotide reduction, it retained putative fitness genes, which became predominant in the 2022 outbreak. Overall, we demonstrated, through rapid whole-genome sequencing and analysis, a recent, traceable event of genome reduction for niche adaptation typically observed over millennia in obligate and fastidious pathogens.IMPORTANCEThe geranium industry, valued at $4 million annually, faces an ongoing Xanthomonas hortorum pv. pelargonii (Xhp) pathogen outbreak. To track and describe the outbreak, we compared the genome structure across historical and globally distributed isolates. Our research revealed Xhp population has not had chromosome rearrangements since 1974 and has three distinct plasmids. In 2012, we found all three plasmids in individual Xhp isolates. However, in 2022, the three plasmids co-integrated into one plasmid named p93. p93 retained putative fitness genes but lost extraneous genomic material. Our findings show that the 2022 strain group of the bacterial plant pathogen Xanthomonas hortorum underwent a plasmid reduction. We also observed several Xanthomonas species from different years, hosts, and continents have similar plasmids to p93, possibly due to shared agricultural settings. We noticed parallels between genome efficiency and reduction that we see across millennia with obligate parasites with increased niche specificity.

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.

Recombination in Bacterial Genomes: Evolutionary Trends

Bacterial organisms have undergone homologous recombination (HR) and horizontal gene transfer (HGT) multiple times during their history. These processes could increase fitness to new environments, cause specialization, the emergence of new species, and changes in virulence. Therefore, comprehensive knowledge of the impact and intensity of genetic exchanges and the location of recombination hotspots on the genome is necessary for understanding the dynamics of adaptation to various conditions. To this end, we aimed to characterize the functional impact and genomic context of computationally detected recombination events by analyzing genomic studies of any bacterial species, for which events have been detected in the last 30 years. Genomic loci where the transfer of DNA was detected pertained to mobile genetic elements (MGEs) housing genes that code for proteins engaged in distinct cellular processes, such as secretion systems, toxins, infection effectors, biosynthesis enzymes, etc. We found that all inferences fall into three main lifestyle categories, namely, ecological diversification, pathogenesis, and symbiosis. The latter primarily exhibits ancestral events, thus, possibly indicating that adaptation appears to be governed by similar recombination-dependent mechanisms.

Pan-genome analysis of the Burkholderia gladioli PV. Cocovenenans reveal the extent of variation in the toxigenic gene cluster

Burkholderia gladioli has been reported as the pathogen responsible for cases of foodborne illness in many countries. The poisonous bongkrekic acid (BA) produced by B. gladioli was linked to a gene cluster absent in non-pathogenic strains. The whole genome sequence of eight bacteria strains, which were screened from the collected 175 raw food and environmental samples, were assembled and analyzed to detect a significant association of 19 protein-coding genes with the pathogenic status. Except for the common BA synthesis-related gene, several other genes, including the toxin-antitoxin genes, were also absent in the non-pathogenic strains. The bacteria strains with the BA gene cluster were found to form a single cluster in the analysis of all B. gladioli genome assemblies for the variants in the gene cluster. Divergence of this cluster was detected in the analysis for both the flanking sequences and those of the whole genome level, which indicates its complex origin. Genome recombination was found to cause a precise sequence deletion in the gene cluster region, which was found to be predominant in the non-pathogenic strains indicating the possible effect of horizontal gene transfer. Our study provided new information and resources for understanding the evolution and divergence of the B. gladioli species.

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.

An IS-mediated, RecA-dependent, bet-hedging strategy in Burkholderia thailandensis

Adaptation to fluctuating environmental conditions is difficult to achieve. Phase variation mechanisms can overcome this difficulty by altering genomic architecture in a subset of individuals, creating a phenotypically heterogeneous population with subpopulations optimized to persist when conditions change, or are encountered, suddenly. We have identified a phase variation system in Burkholderia thailandensis that generates a genotypically and phenotypically heterogeneous population. Genetic analyses revealed that RecA-mediated homologous recombination between a pair of insertion sequence (IS) 2-like elements duplicates a 208.6 kb region of DNA that contains 157 coding sequences. RecA-mediated homologous recombination also resolves merodiploids, and hence copy number of the region is varied and dynamic within populations. We showed that the presence of two or more copies of the region is advantageous for growth in a biofilm, and a single copy is advantageous during planktonic growth. While IS elements are well known to contribute to evolution through gene inactivation, polar effects on downstream genes, and altering genomic architecture, we believe that this system represents a rare example of IS element-mediated evolution in which the IS elements provide homologous sequences for amplification of a chromosomal region that provides a selective advantage under specific growth conditions, thereby expanding the lifestyle repertoire of the species.

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.

Reduced and Nonreduced Genomes in Paraburkholderia Symbionts of Social Amoebas

The social amoeba Dictyostelium discoideum is a predatory soil protist frequently used for studying host-pathogen interactions. A subset of D. discoideum strains isolated from soil persistently carry symbiotic Paraburkholderia, recently formally described as P. agricolaris, P. bonniea, and P. hayleyella. The three facultative symbiont species of D. discoideum present a unique opportunity to study a naturally occurring symbiosis in a laboratory model protist. There is a large difference in genome size between P. agricolaris (8.7 million base pairs [Mbp]) versus P. hayleyella and P. bonniea (4.1 Mbp). We took a comparative genomics approach and compared the three genomes of D. discoideum symbionts to 12 additional Paraburkholderia genomes to test for genome evolution patterns that frequently accompany host adaptation. Overall, P. agricolaris is difficult to distinguish from other Paraburkholderia based on its genome size and content, but the reduced genomes of P. bonniea and P. hayleyella display characteristics indicative of genome streamlining rather than deterioration during adaptation to their protist hosts. In addition, D. discoideum-symbiont genomes have increased secretion system and motility genes that may mediate interactions with their host. Specifically, adjacent BurBor-like type 3 and T6SS-5-like type 6 secretion system operons shared among all three D. discoideum-symbiont genomes may be important for host interaction. Horizontal transfer of these secretion system operons within the amoeba host environment may have contributed to the unique ability of these symbionts to establish and maintain a symbiotic relationship with D. discoideum. IMPORTANCE Protists are a diverse group of typically single cell eukaryotes. Bacteria and archaea that form long-term symbiotic relationships with protists may evolve in additional ways than those in relationships with multicellular eukaryotes such as plants, animals, or fungi. Social amoebas are a predatory soil protist sometimes found with symbiotic bacteria living inside their cells. They present a unique opportunity to explore a naturally occurring symbiosis in a protist frequently used for studying host-pathogen interactions. We show that one amoeba-symbiont species is similar to other related bacteria in genome size and content, while the two reduced-genome-symbiont species show characteristics of genome streamlining rather than deterioration during adaptation to their host. We also identify sets of genes present in all three amoeba-symbiont genomes that are potentially used for host-symbiont interactions. Because the amoeba symbionts are distantly related, the amoeba host environment may be where these genes were shared among symbionts.

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

Burkholderia mallei is the etiological agent of glanders, a highly contagious and often fatal disease in equids. Due to the high genetic clonality of B. mallei, high-resolution typing assays are necessary to differentiate between individual strains. Here we report on the development and validation of a robust and reproducible core genome-based Multi Locus Sequence Typing Assay (cgMLST) for B. mallei, which is based on 3328 gene targets and enables high-resolution typing at the strain level. The assay was validated using a set of 120 B. mallei genomes from public databases and 23 newly sequenced outbreak strains from in-house strain collections. In this cgMLST analysis, strains from different geographic regions were clearly distinguished by at least 70 allele differences, allowing spatial clustering while closely related and epidemiologically related strains were separated by only zero to three alleles. Neither the different sequencing technologies nor the assembly strategies had an influence on the cgMLST results. The developed cgMLST is highly robust, reproducible and can be used for outbreak investigations, source tracking and molecular characterization of new B. mallei isolates.

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

Burkholderia (B.) mallei is a host-adapted equine pathogen that causes glanders, a re-emerging zoonotic disease, which is endemic in Pakistan and other developing countries and seriously impacts the global equine movement. Due to globalization, the geographical restriction of diseases vanishes and the lack of awareness of and experience with eradicated diseases in industrialized countries also promotes the re-introduction of infections in these regions. Owing to the high equine population, the Pakistani province Punjab is a potential hotspot where several glanders outbreaks have been seen over last two decades. For determining the genomic diversity of B. mallei in this and other equine-populated prefectures, the genomes of 19 B. mallei strains isolated between 1999 and 2020 in different locations were sequenced and their genotypes were determined. Particularly, for genetically highly homogenous pathogens like B. mallei genotyping techniques require a high discriminatory power for enabling differentiation on the strain level. Thus, core-genome single nucleotide polymorphism (cgSNP) analysis was applied for distinguishing the highly similar strains. Furthermore, a whole-genome sequence-based core genome multi locus sequence typing (cgMLST) scheme, specific to B. mallei, was developed and additionally applied to the data. It was found that B. mallei genotypes in Pakistan persisted over time and space and genotype clusters preferred connection with a time point rather than the place of isolation, probably due to frequent equine movement, which promotes the spread of glanders. The cgMLST approach proved to work in accord with SNP typing and may help to investigate future glanders outbreaks.

Assessment of the control measures of the category A diseases of Animal Health Law: Burkholderia mallei (Glanders)

EFSA received a mandate from the European Commission to assess the effectiveness of some of the control measures against diseases included in the Category A list according to Regulation (EU) 2016/429 on transmissible animal diseases ('Animal Health Law'). This opinion belongs to a series of opinions where these control measures will be assessed, with this opinion covering the assessment of control measures for glanders. In this opinion, EFSA and the AHAW Panel of experts review the effectiveness of: (i) clinical and laboratory sampling procedures, (ii) monitoring period and (iii) the minimum radius of the protection and surveillance zone, and the minimum length of time the measures should be applied in these zones. The general methodology used for this series of opinions has been published elsewhere. Considering the epidemiology and distribution of glanders, it was foreseen that three different situations could lead to a suspicion of the disease. Sampling procedures were defined for each of the three different suspicion types, which can also be applied in most of the other scenarios assessed. The monitoring period (6 months) was assessed as effective in all scenarios. The AHAW Panel of experts considered the minimum radius and duration of the existing protection and surveillance zone, set at the establishment level, effective. Recommendations provided for each of the scenarios assessed aim to support the European Commission in the drafting of further pieces of legislation, as well as for plausible ad hoc requests in relation to glanders.

Molecular epidemiology of Burkholderia mallei isolates from India (2015-2016): New SNP markers for strain tracing

Glanders, caused by a bacterium called B. mallei, is primarily an infectious horse and human disease. Although its incidence is rare in developed countries, it is nonetheless prevalent in several geographical areas of the world. There is a lack of cost-effective, rapid and specific molecular typing tools for epidemiological tracing of glanders cases. We previously reported an SNP-based typing method that categorizes global B. mallei strains into three lineages (L1 to L3), as well as additional branches, sub-branches and groups. However, further discrimination of the Indian and Pakistani isolates within the L2B2sB2 sub-branch was not possible due to the lack of sufficient epidemiological markers. In this study, 10 B. mallei strains isolated from four states in India during 2015-2016 were whole genome sequenced; SNP analysis further confirmed their position in the L2B2sB2 branch. To better track the strains, four new markers targeting Indian or Pakistani strains, and specifically targeting sub-groups within the Indian strains, were identified. The new SNP markers were tested and validated on the 10 Indian isolates included in this study as well as on 6 contemporary B. mallei Pakistani strains. These rapid and discriminating typing tools will contribute to the epidemiological monitoring of B. mallei infections, particularly in South Asia and the Middle East, endemic regions of the disease.

The Burkholderia pseudomallei hmqA-G Locus Mediates Competitive Fitness against Environmental Gram-Positive Bacteria

Burkholderia pseudomallei is an opportunistic pathogen that is responsible for the disease melioidosis in humans and animals. The microbe is a tier 1 select agent because it is highly infectious by the aerosol route, it is inherently resistant to multiple antibiotics, and no licensed vaccine currently exists. Naturally acquired infections result from contact with contaminated soil or water sources in regions of endemicity. There have been few reports investigating the molecular mechanism(s) utilized by B. pseudomallei to survive and persist in ecological niches harboring microbial competitors. Here, we report the isolation of Gram-positive bacteria from multiple environmental sources and show that ∼45% of these isolates are inhibited by B. pseudomallei in head-to-head competition assays. Two competition-deficient B. pseudomallei transposon mutants were identified that contained insertion mutations in the hmqA-G operon. This large biosynthetic gene cluster encodes the enzymes that produce a family of secondary metabolites called 4-hydroxy-3-methyl-2-alkylquinolines (HMAQs). Liquid chromatography and mass spectrometry conducted on filter-sterilized culture supernatants revealed five HMAQs and N-oxide derivatives that were produced by the parental strain but were absent in an isogenic hmqD deletion mutant. The results demonstrate that B. pseudomallei inhibits the growth of environmental Gram-positive bacteria in a contact-independent manner via the production of HMAQs by the hmqA-G operon. IMPORTANCE Burkholderia pseudomallei naturally resides in water, soil, and the rhizosphere and its success as an opportunistic pathogen is dependent on the ability to persist in these harsh habitats long enough to come into contact with a susceptible host. In addition to adapting to limiting nutrients and diverse chemical and physical challenges, B. pseudomallei also has to interact with a variety of microbial competitors. Our research shows that one of the ways in which B. pseudomallei competes with Gram-positive environmental bacteria is by exporting a diverse array of closely related antimicrobial secondary metabolites.

The Peptidoglycan-associated lipoprotein Pal contributes to the virulence of Burkholderia mallei and provides protection against lethal aerosol challenge

Burkholderia mallei

is a highly pathogenic bacterium that causes the fatal zoonosis glanders. The organism specifies multiple membrane proteins, which represent prime targets for the development of countermeasures given their location at the host-pathogen interface. We investigated one of these proteins, Pal, and discovered that it is involved in the ability of B. mallei to resist complement-mediated killing and replicate inside host cells in vitro, is expressed in vivo and induces antibodies during the course of infection, and contributes to virulence in a mouse model of aerosol infection. A mutant in the pal gene of the B. mallei wild-type strain ATCC 23344 was found to be especially attenuated, as BALB/c mice challenged with the equivalent of 5,350 LD₅₀ completely cleared infection. Based on these findings, we tested the hypothesis that a vaccine containing the Pal protein elicits protective immunity against aerosol challenge. To achieve this, the pal gene was cloned in the vaccine vector Parainfluenza Virus 5 (PIV5) and mice immunized with the virus were infected with a lethal dose of B. mallei. These experiments revealed that a single dose of PIV5 expressing Pal provided 80% survival over a period of 40 days post-challenge. In contrast, only 10% of mice vaccinated with a PIV5 control virus construct survived infection. Taken together, our data establish that the Peptidoglycan-associated lipoprotein Pal is a critical virulence determinant of B. mallei and effective target for developing a glanders vaccine.

Comparative Genomic Analysis of the Foodborne Pathogen Burkholderia gladioli pv. cocovenenans Harboring a Bongkrekic Acid Biosynthesis Gene Cluster

The environmental bacterium Burkholderia gladioli pv. cocovenenans (B. cocovenenans) has been linked to fatal food poisoning cases in Asia and Africa. Bongkrekic acid (BA), a mitochondrial toxin produced by B. cocovenenans, is thought to be responsible for these outbreaks. While there are over 80 species in the Burkholderia genus, B. cocovenenans is the only pathovar capable of producing BA and causing human death. However, the genomic features of B. gladioli and the evolution of the BA biosynthesis gene cluster, bon, in B. cocovenenans remain elusive. In this study, 239 whole genome sequences (WGSs) of B. gladioli, isolated from 12 countries collected over 100 years, were used to analyze the intra-species genomic diversity and phylogenetic relationships of B. gladioli and to explore the origin and evolution of the bon gene cluster. Our results showed that the genome-wide average nucleotide identity (ANI) values were above 97.29% for pairs of B. gladioli genomes. Thirty-six of the 239 (15.06%) B. gladioli genomes, isolated from corn, rice, fruits, soil, and patients from Asia, Europe, North America, and South America, contained the bon gene cluster and formed three clades within the phylogenetic tree. Pan- and core-genome analysis suggested that the BA biosynthesis genes were recently acquired. Comparative genome analysis of the bon gene cluster showed that complex recombination events contributed to this toxin biosynthesis gene cluster's evolution and formation. This study suggests that a better understanding of the genomic diversity and evolution of this lethal foodborne pathovar will potentially contribute to B. cocovenenans food poisoning outbreak prevention.

High Rates of Genome Rearrangements and Pathogenicity of Shigella spp

Shigella are pathogens originating within the Escherichia lineage but frequently classified as a separate genus. Shigella genomes contain numerous insertion sequences (ISs) that lead to pseudogenisation of affected genes and an increase of non-homologous recombination. Here, we study 414 genomes of E. coli and Shigella strains to assess the contribution of genomic rearrangements to Shigella evolution. We found that Shigella experienced exceptionally high rates of intragenomic rearrangements and had a decreased rate of homologous recombination compared to pathogenic and non-pathogenic E. coli. The high rearrangement rate resulted in independent disruption of syntenic regions and parallel rearrangements in different Shigella lineages. Specifically, we identified two types of chromosomally encoded E3 ubiquitin-protein ligases acquired independently by all Shigella strains that also showed a high level of sequence conservation in the promoter and further in the 5′-intergenic region. In the only available enteroinvasive E. coli (EIEC) strain, which is a pathogenic E. coli with a phenotype intermediate between Shigella and non-pathogenic E. coli, we found a rate of genome rearrangements comparable to those in other E. coli and no functional copies of the two Shigella-specific E3 ubiquitin ligases. These data indicate that the accumulation of ISs influenced many aspects of genome evolution and played an important role in the evolution of intracellular pathogens. Our research demonstrates the power of comparative genomics-based on synteny block composition and an important role of non-coding regions in the evolution of genomic islands.

Burkholderia ubonensis High-Level Tetracycline Resistance Is Due to Efflux Pump Synergy Involving a Novel TetA(64) Resistance Determinant

Burkholderia ubonensis, a nonpathogenic soil bacterium belonging to the Burkholderia cepacia complex (Bcc), is highly resistant to some clinically significant antibiotics. The concern is that B. ubonensis may serve as a resistance reservoir for Bcc or B. pseudomallei complex (Bpc) organisms that are opportunistic human pathogens. Using a B. ubonensis strain highly resistant to tetracycline (MIC, ≥256 µg/ml), we identified and characterized tetA(64) that encodes a novel tetracycline-specific efflux pump of the major facilitator superfamily. TetA(64) and associated TetR(64) regulator expression are induced by tetracyclines. Although TetA(64) is the primary tetracycline and doxycycline resistance determinant, maximum tetracycline and doxycycline resistance requires synergy between TetA(64) and the nonspecific AmrAB-OprA resistance nodulation cell division efflux pump. TetA(64) does not efflux minocycline, tigecycline, and eravacycline. Comprehensive screening of genome sequences showed that TetA(64) is unequally distributed in the Bcc and absent from the Bpc. It is present in some major cystic fibrosis pathogens, like Burkholderia cenocepacia, but absent from others like Burkholderia multivorans The tetR(64)-tetA(64) genes are located in a region of chromosome 1 that is highly conserved in Burkholderia sp. Because there is no evidence for transposition, the tetR(64)-tetA(64) genes may have been acquired by homologous recombination after horizontal gene transfer. Although Burkholderia species contain a resident multicomponent efflux pump that allows them to respond to tetracyclines up to a certain concentration, the acquisition of the single-component TetA(64) by some species likely provides the synergy that these bacteria need to defend against high tetracycline concentrations in niche environments.

An Evolutionary Arms Race Between Burkholderia pseudomallei and Host Immune System: What Do We Know?

A better understanding of co-evolution between pathogens and hosts holds promise for better prevention and control strategies. This review will explore the interactions between Burkholderia pseudomallei, an environmental and opportunistic pathogen, and the human host immune system. B. pseudomallei causes "Melioidosis," a rapidly fatal tropical infectious disease predicted to affect 165,000 cases annually worldwide, of which 89,000 are fatal. Genetic heterogeneities were reported in both B. pseudomallei and human host population, some of which may, at least in part, contribute to inter-individual differences in disease susceptibility. Here, we review (i) a multi-host-pathogen characteristic of the interaction; (ii) selection pressures acting on B. pseudomallei and human genomes with the former being driven by bacterial adaptation across ranges of ecological niches while the latter are driven by human encounter of broad ranges of pathogens; (iii) the mechanisms that generate genetic diversity in bacterial and host population particularly in sequences encoding proteins functioning in host-pathogen interaction; (iv) reported genetic and structural variations of proteins or molecules observed in B. pseudomallei-human host interactions and their implications in infection outcomes. Together, these predict bacterial and host evolutionary trajectory which continues to generate genetic diversity in bacterium and operates host immune selection at the molecular level.

Burkholderia collagen-like protein 8, Bucl8, is a unique outer membrane component of a putative tetrapartite efflux pump in Burkholderia pseudomallei and Burkholderia mallei

Bacterial efflux pumps are an important pathogenicity trait because they extrude a variety of xenobiotics. Our laboratory previously identified in silico Burkholderia collagen-like protein 8 (Bucl8) in the hazardous pathogens Burkholderia pseudomallei and Burkholderia mallei. We hypothesize that Bucl8, which contains two predicted tandem outer membrane efflux pump domains, is a component of a putative efflux pump. Unique to Bucl8, as compared to other outer membrane proteins, is the presence of an extended extracellular region containing a collagen-like (CL) domain and a non-collagenous C-terminus (Ct). Molecular modeling and circular dichroism spectroscopy with a recombinant protein, corresponding to this extracellular CL-Ct portion of Bucl8, demonstrated that it adopts a collagen triple helix, whereas functional assays screening for Bucl8 ligands identified binding to fibrinogen. Bioinformatic analysis of the bucl8 gene locus revealed it resembles a classical efflux-pump operon. The bucl8 gene is co-localized with downstream fusCDE genes encoding fusaric acid (FA) resistance, and with an upstream gene, designated as fusR, encoding a LysR-type transcriptional regulator. Using reverse transcriptase (RT)-qPCR, we defined the boundaries and transcriptional organization of the fusR-bucl8-fusCDE operon. We found exogenous FA induced bucl8 transcription over 80-fold in B. pseudomallei, while deletion of the entire bucl8 locus decreased the minimum inhibitory concentration of FA 4-fold in its isogenic mutant. We furthermore showed that the putative Bucl8-associated pump expressed in the heterologous Escherichia coli host confers FA resistance. On the contrary, the Bucl8-associated pump did not confer resistance to a panel of clinically-relevant antimicrobials in Burkholderia and E. coli. We finally demonstrated that deletion of the bucl8-locus drastically affects the growth of the mutant in L-broth. We determined that Bucl8 is a component of a novel tetrapartite efflux pump, which confers FA resistance, fibrinogen binding, and optimal growth.

A genetic variant of Burkholderia mallei detected in Kuwait: Consequences for the PCR diagnosis of glanders

Glanders is a contagious zoonotic disease caused by Burkholderia mallei. Following the detection of glanders positive horses using the OIE complement fixation test, the tissues of two horses were analysed by PCR. While PCR systems targeting the Burkholderia pseudomallei complex gave positive signals, the species-specific PCR systems targeting B. mallei (fliP-IS407A) and B. pseudomallei (orf11)-the OIE recommended targets-resulted in negative signals. However, the presence of B. mallei in these tissues was confirmed with a recently described B. mallei-specific real-time PCR system and genotyping with MLST- and SNP-based methods, performed on the most positive tissue, identified a genotype closely related to B. mallei strains recently isolated in the Middle East. This study leads to recommendations regarding the use of PCR systems for the molecular diagnosis of glanders, especially in regions where the circulating B. mallei strains have not yet been fully genetically characterized.

Genomic loss in environmental and isogenic morphotype isolates of Burkholderia pseudomallei is associated with intracellular survival and plaque-forming efficiency

BACKGROUND

Burkholderia pseudomallei is an environmental bacterium that causes melioidosis. A facultative intracellular pathogen, B. pseudomallei can induce multinucleated giant cells (MNGCs) leading to plaque formation in vitro. B. pseudomallei can switch colony morphotypes under stress conditions. In addition, different isolates have been reported to have varying virulence in vivo, but genomic evolution and the relationship with plaque formation is poorly understood.

METHODOLOGY/PRINCIPLE FINDINGS

To gain insights into genetic underpinnings of virulence of B. pseudomallei, we screened plaque formation of 52 clinical isolates and 11 environmental isolates as well as 4 isogenic morphotype isolates of B. pseudomallei strains K96243 (types II and III) and 153 (types II and III) from Thailand in A549 and HeLa cells. All isolates except one environmental strain (A4) and K96243 morphotype II were able to induce plaque formation in both cell lines. Intracellular growth assay and confocal microscopy analyses demonstrated that the two plaque-forming-defective isolates were also impaired in intracellular replication, actin polymerization and MNGC formation in infected cells. Whole genome sequencing analysis and PCR revealed that both isolates had a large genomic loss on the same region in chromosome 2, which included Bim cluster, T3SS-3 and T6SS-5 genes.

CONCLUSIONS/SIGNIFICANCE

Our plaque screening and genomic studies revealed evidence of impairment in plaque formation in environmental isolates of B. pseudomallei that is associated with large genomic loss of genes important for intracellular multiplication and MNGC formation. These findings suggest that the genomic and phenotypic differences of environmental isolates may be associated with clinical infection.

A New Phylogenomic Approach For Quantifying Horizontal Gene Transfer Trends in Prokaryotes

It is well established nowadays that among prokaryotes, various families of orthologous genes exhibit conflicting evolutionary history. A prime factor for this conflict is horizontal gene transfer (HGT) - the transfer of genetic material not via vertical descent. Thus, the prevalence of HGT is challenging the meaningfulness of the classical Tree of Life concept. Here we present a comprehensive study of HGT representing the entire prokaryotic world. We mainly rely on a novel analytic approach for analyzing an aggregate of gene histories, by means of the quartet plurality distribution (QPD) that we develop. Through the analysis of real and simulated data, QPD is used to reveal evidence of a barrier against HGT, separating the archaea from the bacteria and making HGT between the two domains, in general, quite rare. In contrast, bacteria's confined HGT is substantially more frequent than archaea's. Our approach also reveals that despite intensive HGT, a strong tree-like signal can be extracted, corroborating several previous works. Thus, QPD, which enables one to analytically combine information from an aggregate of gene trees, can be used for understanding patterns and rates of HGT in prokaryotes, as well as for validating or refuting models of horizontal genetic transfers and evolution in general.

Development of a microsphere-based immunoassay for the serological detection of glanders in equids

Burkholderia mallei is the etiologic agent of glanders, an infectious disease of solipeds, with renewed scientific interest due to its increasing incidence in different parts of the world. More rapid, sensitive and specific assays are required by laboratories for confirmatory testing of this disease. A microsphere-based immunoassay consisting of beads coated with B. mallei recombinant proteins (BimA, GroEL, Hcp1, and TssB) has been developed for the serological diagnosis of glanders. The proteins' performance was compared with the OIE reference complement fixation test (CFT) and an indirect enzyme-linked immunosorbent assay (iELISA) on a large panel of sera comprised of uninfected horses (n=198) and clinically confirmed cases of glanders from India and Pakistan (n=99). Using Receiver Operating Characteristics (ROC) analysis and adjusting the cutoff levels, Hcp1 (Se=100%, Sp=99.5%) and GroEL (Se= 97%, Sp=99.5%) antigens exhibited the best specificity and sensitivity. Neither Hcp1 and GroEL proteins, nor iELISA reacted with doubtful and positive CFT samples from glanders free countries which further confirmed the false positive reactions seen in CFT.

Pathogen to commensal? Longitudinal within-host population dynamics, evolution, and adaptation during a chronic >16-year Burkholderia pseudomallei infection

Although acute melioidosis is the most common outcome of Burkholderia pseudomallei infection, we have documented a case, P314, where disease severity lessened with time, and the pathogen evolved towards a commensal relationship with the host. In the current study, we used whole-genome sequencing to monitor this long-term symbiotic relationship to better understand B. pseudomallei persistence in P314’s sputum despite intensive initial therapeutic regimens. We collected and sequenced 118 B. pseudomallei isolates from P314’s airways over a >16-year period, and also sampled the patient’s home environment, recovering six closely related B. pseudomallei isolates from the household water system. Using comparative genomics, we identified 126 SNPs in the core genome of the 124 isolates or 162 SNPs/indels when the accessory genome was included. The core SNPs were used to construct a phylogenetic tree, which demonstrated a close relationship between environmental and clinical isolates and detailed within-host evolutionary patterns. The phylogeny had little homoplasy, consistent with a strictly clonal mode of genetic inheritance. Repeated sampling revealed evidence of genetic diversification, but frequent extinctions left only one successful lineage through the first four years and two lineages after that. Overall, the evolution of this population is nonadaptive and best explained by genetic drift. However, some genetic and phenotypic changes are consistent with in situ adaptation. Using a mouse model, P314 isolates caused greatly reduced morbidity and mortality compared to the environmental isolates. Additionally, potentially adaptive phenotypes emerged and included differences in the O-antigen, capsular polysaccharide, motility, and colony morphology. The >13-year co-existence of two long-lived lineages presents interesting hypotheses that can be tested in future studies to provide additional insights into selective pressures, niche differentiation, and microbial adaptation. This unusual melioidosis case presents a rare example of the evolutionary progression towards commensalism by a highly virulent pathogen within a single human host.

Pathogens frequently jump between different hosts, and associated adaptation may lead to the emergence of new infectious agents. Such host-jumping evolution is witnessed through endpoint analyses but these cannot capture genetic changes in lineages that have gone extinct. In this study, we have identified and monitored an example of the evolution of a bacterium often deadly to its mammalian host, in an unprecedented case whereby disease lessened through time and the pathogen became a part of the commensal human flora. We used genomic analyses to characterize more than 16 years of this evolutionary process and the stepwise mutations that control pathogen interactions with the patient. Soon after infection, mutational changes occurred that allowed the bacterium to remain in the airways without causing disease. This shift towards avirulence was determined based on clinical data and virulence testing in an animal model. In addition, mutations occurred that contributed to the persistence of the bacteria in the patient's lungs. Finally, we found evidence for the evolutionary emergence and persistence of two distinct lineages of the bacterium over the last 13 years, presenting interesting questions about niche utilization. Bacteria are ubiquitous in the human body and almost all are beneficial or benign. In this study, we document the evolutionary conversion of a normally deadly bacterium towards a commensal.

Satellites in the prokaryote world

Background

Satellites or tandem repeats are very abundant in many eukaryotic genomes. Occasionally they have been reported to be present in some prokaryotes, but to our knowledge there is no general comparative study on their occurrence. For this reason we present here an overview of the distribution and properties of satellites in a set of representative species. Our results provide novel insights into the evolutionary relationship between eukaryotes, Archaea and Bacteria.

Results

We have searched all possible satellites present in the NCBI reference group of genomes in Archaea (142 species) and in Bacteria (119 species), detecting 2735 satellites in Archaea and 1067 in Bacteria. We have found that the distribution of satellites is very variable in different organisms. The archaeal Methanosarcina class stands out for the large amount of satellites in their genomes. Satellites from a few species have similar characteristics to those in eukaryotes, but most species have very few satellites: only 21 species in Archaea and 18 in Bacteria have more than 4 satellites/Mb. The distribution of satellites in these species is reminiscent of what is found in eukaryotes, but we find two significant differences: most satellites have a short length and many of them correspond to segments of genes coding for amino acid repeats. Transposition of non-coding satellites throughout the genome occurs rarely: only in the bacteria Leptospira interrogans and the archaea Methanocella conradii we have detected satellite families of transposed satellites with long repeats.

Conclusions

Our results demonstrate that the presence of satellites in the genome is not an exclusive feature of eukaryotes. We have described a few prokaryotes which do contain satellites. We present a discussion on their eventual evolutionary significance.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1504-2) contains supplementary material, which is available to authorized users.

An in situ high-throughput screen identifies inhibitors of intracellular Burkholderia pseudomallei with therapeutic efficacy

Burkholderia pseudomallei, the etiologic agent of melioidosis, is an environmental organism that inhabits tropical soils and kills an estimated 90,000 people each year. Caused by an intracellular and often drug-resistant pathogen, melioidosis is notoriously difficult to treat, with mortality rates approaching 50% in some settings despite appropriate diagnosis and clinical management. Using a high-throughput, cell-based phenotypic screen we have discovered 2 antibiotic candidates with improved in vivo efficacy compared to the current standard of care: a fluoroquinolone analog, burkfloxacin, and an FDA-approved antifungal drug, flucytosine. As a widely used antifungal with a well-known safety profile, the potential to repurpose flucytosine for treating melioidosis may represent a rapid route to clinical translation.

Burkholderia pseudomallei (Bp) and Burkholderia mallei (Bm) are Tier-1 Select Agents that cause melioidosis and glanders, respectively. These are highly lethal human infections with limited therapeutic options. Intercellular spread is a hallmark of Burkholderia pathogenesis, and its prominent ties to virulence make it an attractive therapeutic target. We developed a high-throughput cell-based phenotypic assay and screened ∼220,000 small molecules for their ability to disrupt intercellular spread by Burkholderia thailandensis, a closely related BSL-2 surrogate. We identified 268 hits, and cross-species validation found 32 hits that also disrupt intercellular spread by Bp and/or Bm. Among these were a fluoroquinolone analog, which we named burkfloxacin (BFX), which potently inhibits growth of intracellular Burkholderia, and flucytosine (5-FC), an FDA-approved antifungal drug. We found that 5-FC blocks the intracellular life cycle at the point of type VI secretion system 5 (T6SS-5)-mediated cell–cell spread. Bacterial conversion of 5-FC to 5-fluorouracil and subsequently to fluorouridine monophosphate is required for potent and selective activity against intracellular Burkholderia. In a murine model of fulminant respiratory melioidosis, treatment with BFX or 5-FC was significantly more effective than ceftazidime, the current antibiotic of choice, for improving survival and decreasing bacterial counts in major organs. Our results demonstrate the utility of cell-based phenotypic screening for Select Agent drug discovery and warrant the advancement of BFX and 5-FC as candidate therapeutics for melioidosis in humans.

A general protein O-glycosylation machinery conserved in Burkholderia species improves bacterial fitness and elicits glycan immunogenicity in humans

The Burkholderia genus encompasses many Gram-negative bacteria living in the rhizosphere. Some Burkholderia species can cause life-threatening human infections, highlighting the need for clinical interventions targeting specific lipopolysaccharide proteins. Burkholderia cenocepacia O-linked protein glycosylation has been reported, but the chemical structure of the O-glycan and the machinery required for its biosynthesis are unknown and could reveal potential therapeutic targets. Here, using bioinformatics approaches, gene-knockout mutants, purified recombinant proteins, LC-MS-based analyses of O-glycans, and NMR-based structural analyses, we identified a B. cenocepacia O-glycosylation (ogc) gene cluster necessary for synthesis, assembly, and membrane translocation of a lipid-linked O-glycan, as well as its structure, which consists of a β-Gal-(1,3)-α-GalNAc-(1,3)-β-GalNAc trisaccharide. We demonstrate that the ogc cluster is conserved in the Burkholderia genus, and we confirm the production of glycoproteins with similar glycans in the Burkholderia species: B. thailandensis, B. gladioli, and B. pseudomallei Furthermore, we show that absence of protein O-glycosylation severely affects bacterial fitness and accelerates bacterial clearance in a Galleria mellonella larva infection model. Finally, our experiments revealed that patients infected with B. cenocepacia, Burkholderia multivorans, B. pseudomallei, or Burkholderia mallei develop O-glycan-specific antibodies. Together, these results highlight the importance of general protein O-glycosylation in the biology of the Burkholderia genus and its potential as a target for inhibition or immunotherapy approaches to control Burkholderia infections.

On the Impact of the Pangenome and Annotation Discrepancies While Building Protein Sequence Databases for Bacteria Proteogenomics

In proteomics, peptide information within mass spectrometry (MS) data from a specific organism sample is routinely matched against a protein sequence database that best represent such organism. However, if the species/strain in the sample is unknown or genetically poorly characterized, it becomes challenging to determine a database which can represent such sample. Building customized protein sequence databases merging multiple strains for a given species has become a strategy to overcome such restrictions. However, as more genetic information is publicly available and interesting genetic features such as the existence of pan- and core genes within a species are revealed, we questioned how efficient such merging strategies are to report relevant information. To test this assumption, we constructed databases containing conserved and unique sequences for 10 different species. Features that are relevant for probabilistic-based protein identification by proteomics were then monitored. As expected, increase in database complexity correlates with pangenomic complexity. However, Mycobacterium tuberculosis and Bordetella pertussis generated very complex databases even having low pangenomic complexity. We further tested database performance by using MS data from eight clinical strains from M. tuberculosis, and from two published datasets from Staphylococcus aureus. We show that by using an approach where database size is controlled by removing repeated identical tryptic sequences across strains/species, computational time can be reduced drastically as database complexity increases.

The autotransporter protein BatA is a protective antigen against lethal aerosol infection with Burkholderia mallei and Burkholderia pseudomallei

BACKGROUND

Burkholderia mallei and Burkholderia pseudomallei are the causative agents of glanders and melioidosis, respectively. There is no vaccine to protect against these highly-pathogenic and intrinsically antibiotic-resistant bacteria, and there is concern regarding their use as biological warfare agents. For these reasons, B. mallei and B. pseudomallei are classified as Tier 1 organisms by the U.S. Federal Select Agent Program and the availability of effective countermeasures represents a critical unmet need.

METHODS

Vaccines (subunit and vectored) containing the surface-exposed passenger domain of the conserved Burkholderia autotransporter protein BatA were administered to BALB/c mice and the vaccinated animals were challenged with lethal doses of wild-type B. mallei and B. pseudomallei strains via the aerosol route. Mice were monitored for signs of illness for a period of up to 40 days post-challenge and tissues from surviving animals were analyzed for bacterial burden at study end-points.

RESULTS

A single dose of recombinant Parainfluenza Virus 5 (PIV5) expressing BatA provided 74% and 60% survival in mice infected with B. mallei and B. pseudomallei, respectively. Vaccination with PIV5-BatA also resulted in complete bacterial clearance from the lungs and spleen of 78% and 44% of animals surviving lethal challenge with B. pseudomallei, respectively. In contrast, all control animals vaccinated with a PIV5 construct expressing an irrelevant antigen and infected with B. pseudomallei were colonized in those tissues.

CONCLUSION

Our study indicates that the autotransporter BatA is a valuable target for developing countermeasures against B. mallei and B. pseudomallei and demonstrates the utility of the PIV5 viral vaccine delivery platform to elicit cross-protective immunity against the organisms.

Micro-evolution of three Streptococcus species: selection, antigenic variation, and horizontal gene inflow

Background

The genus Streptococcus comprises pathogens that strongly influence the health of humans and animals. Genome sequencing of multiple Streptococcus strains demonstrated high variability in gene content and order even in closely related strains of the same species and created a newly emerged object for genomic analysis, the pan-genome. Here we analysed the genome evolution of 25 strains of Streptococcus suis, 50 strains of Streptococcus pyogenes and 28 strains of Streptococcus pneumoniae.

Results

Fractions of the pan-genome, unique, periphery, and universal genes differ in size, functional composition, the level of nucleotide substitutions, and predisposition to horizontal gene transfer and genomic rearrangements. The density of substitutions in intergenic regions appears to be correlated with selection acting on adjacent genes, implying that more conserved genes tend to have more conserved regulatory regions.

The total pan-genome of the genus is open, but only due to strain-specific genes, whereas other pan-genome fractions reach saturation. We have identified the set of genes with phylogenies inconsistent with species and non-conserved location in the chromosome; these genes are rare in at least one species and have likely experienced recent horizontal transfer between species. The strain-specific fraction is enriched with mobile elements and hypothetical proteins, but also contains a number of candidate virulence-related genes, so it may have a strong impact on adaptability and pathogenicity.

Mapping the rearrangements to the phylogenetic tree revealed large parallel inversions in all species. A parallel inversion of length 15 kB with breakpoints formed by genes encoding surface antigen proteins PhtD and PhtB in S. pneumoniae leads to replacement of gene fragments that likely indicates the action of an antigen variation mechanism.

Conclusions

Members of genus Streptococcus have a highly dynamic, open pan-genome, that potentially confers them with the ability to adapt to changing environmental conditions, i.e. antibiotic resistance or transmission between different hosts. Hence, integrated analysis of all aspects of genome evolution is important for the identification of potential pathogens and design of drugs and vaccines.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1403-6) contains supplementary material, which is available to authorized users.

The Burkholderia Type VI Secretion System 5: Composition, Regulation and Role in Virulence

The soil saprophyte and Tier I select agent Burkholderia pseudomallei can cause rapidly fatal infections in humans and animals. The capability of switching to an intracellular life cycle during infection appears to be a decisive trait of B. pseudomallei for causing disease. B. pseudomallei harbors multiple type VI secretion systems (T6SSs) orthologs of which are present in the surrogate organism Burkholderia thailandensis. Upon host cell entry and vacuolar escape into the cytoplasm, B. pseudomallei and B. thailandensis manipulate host cells by utilizing the T6SS-5 (also termed T6SS1) to form multinucleated giant cells for intercellular spread. Disruption of the T6SS-5 in B. thailandensis causes a drastic attenuation of virulence in wildtype but not in mice lacking the central innate immune adapter protein MyD88. This result suggests that the T6SS-5 is deployed by the bacteria to overcome innate immune responses. However, important questions in this field remain unsolved including the mechanism underlying T6SS-5 activity and its physiological role during infection. In this review, we summarize the current knowledge on the components and regulation of the T6SS-5 as well as its role in virulence in mammalian hosts.

Antibodies Are Major Drivers of Protection against Lethal Aerosol Infection with Highly Pathogenic Burkholderia spp

Burkholderia pseudomallei and Burkholderia mallei are the causative agents of melioidosis and glanders, respectively. There is no vaccine to protect against these highly pathogenic bacteria, and there is concern regarding their emergence as global public health (B. pseudomallei) and biosecurity (B. mallei) threats. In this issue of mSphere, an article by Khakhum and colleagues (N. Khakhum, P. Bharaj, J. N. Myers, D. Tapia, et al., mSphere 4:e00570-18, 2019, https://doi.org/10.1128/mSphere.00570-18) describes a novel vaccination platform with excellent potential for cross-protection against both Burkholderia species. The report also highlights the importance of antibodies in immunity against these facultative intracellular organisms.

Genome rearrangements and selection in multi-chromosome bacteria Burkholderia spp

BACKGROUND

The genus Burkholderia consists of species that occupy remarkably diverse ecological niches. Its best known members are important pathogens, B. mallei and B. pseudomallei, which cause glanders and melioidosis, respectively. Burkholderia genomes are unusual due to their multichromosomal organization, generally comprised of 2-3 chromosomes.

RESULTS

We performed integrated genomic analysis of 127 Burkholderia strains. The pan-genome is open with the saturation to be reached between 86,000 and 88,000 genes. The reconstructed rearrangements indicate a strong avoidance of intra-replichore inversions that is likely caused by selection against the transfer of large groups of genes between the leading and the lagging strands. Translocated genes also tend to retain their position in the leading or the lagging strand, and this selection is stronger for large syntenies. Integrated reconstruction of chromosome rearrangements in the context of strains phylogeny reveals parallel rearrangements that may indicate inversion-based phase variation and integration of new genomic islands. In particular, we detected parallel inversions in the second chromosomes of B. pseudomallei with breakpoints formed by genes encoding membrane components of multidrug resistance complex, that may be linked to a phase variation mechanism. Two genomic islands, spreading horizontally between chromosomes, were detected in the B. cepacia group.

CONCLUSIONS

This study demonstrates the power of integrated analysis of pan-genomes, chromosome rearrangements, and selection regimes. Non-random inversion patterns indicate selective pressure, inversions are particularly frequent in a recent pathogen B. mallei, and, together with periods of positive selection at other branches, may indicate adaptation to new niches. One such adaptation could be a possible phase variation mechanism in B. pseudomallei.

MarR Family Transcription Factors from Burkholderia Species: Hidden Clues to Control of Virulence-Associated Genes

Species within the genus Burkholderia exhibit remarkable phenotypic diversity. Genomic plasticity, including genome reduction and horizontal gene transfer, has been correlated with virulence traits in several species. However, the conservation of virulence genes in species otherwise considered to have limited potential for infection suggests that phenotypic diversity may not be explained solely on the basis of genetic diversity. Instead, differential organization and control of gene regulatory networks may underlie many phenotypic differences. In this review, we evaluate how regulation of gene expression by members of the multiple antibiotic resistance regulator (MarR) family of transcription factors may contribute to shaping the physiological diversity of Burkholderia species, with a focus on the clinically relevant human pathogens. All Burkholderia species encode a relatively large number of MarR proteins, a feature common to bacteria that must respond to environmental changes such as those associated with host invasion. However, evolution of gene regulatory networks has likely resulted in orthologous transcription factors controlling disparate sets of genes. Adaptation to, and survival in, diverse habitats, including a human or plant host, is key to the success of Burkholderia species as (opportunistic) pathogens, and recent reports suggest that control of virulence-associated genes by MarR proteins features prominently among the survival strategies employed by these species. We suggest that identification of MarR regulons will contribute significantly to clarification of virulence determinants and phenotypic diversity.

The heritable natural competency trait of Burkholderia pseudomallei in other Burkholderia species through comE and crp

Natural competency requires uptake of exogenous DNA from the environment and the integration of that DNA into recipient bacteria can be used for DNA-repair or genetic diversification. The Burkholderia genus is unique in that only some of the species and strains are naturally competent. We identified and characterized two genes, comE and crp, from naturally competent B. pseudomallei 1026b that play a role in DNA uptake and catabolism. Single-copies of rhamnose-inducible comE and crp genes were integrated into a Tn7 attachment-site in non-naturally competent Burkholderia including pathogens B. pseudomallei K96243, B. cenocepacia K56-2, and B. mallei ATCC23344. Strains expressing comE or crp were assayed for their ability to uptake and catabolize DNA. ComE and Crp allowed non-naturally competent Burkholderia species to catabolize DNA, uptake exogenous gfp DNA and express GFP. Furthermore, we used synthetic comE and crp to expand the utility of the λ-red recombineering system for genetic manipulation of non-competent Burkholderia species. A newly constructed vector, pKaKa4, was used to mutate the aspartate semialdehyde dehydrogenase (asd) gene in four B. mallei strains, leading to the complete attenuation of these tier-1 select-agents. These strains have been excluded from select-agent regulations and will be of great interest to the field.

High-resolution melting PCR analysis for rapid genotyping of Burkholderia mallei

Burkholderia (B.) mallei is the causative agent of glanders. A previous work conducted on single-nucleotide polymorphisms (SNP) extracted from the whole genome sequences of 45 B. mallei isolates identified 3 lineages for this species. In this study, we designed a high-resolution melting (HRM) method for the screening of 15 phylogenetically informative SNPs within the genome of B. mallei that subtype the species into 3 lineages and 12 branches/sub-branches/groups. The present results demonstrate that SNP-based genotyping represent an interesting approach for the molecular epidemiology analysis of B. mallei.

Antibodies against In Vivo-Expressed Antigens Are Sufficient To Protect against Lethal Aerosol Infection with Burkholderia mallei and Burkholderia pseudomallei

Burkholderia mallei, a facultative intracellular bacterium and tier 1 biothreat, causes the fatal zoonotic disease glanders. The organism possesses multiple genes encoding autotransporter proteins, which represent important virulence factors and targets for developing countermeasures in pathogenic Gram-negative bacteria. In the present study, we investigated one of these autotransporters, BatA, and demonstrate that it displays lipolytic activity, aids in intracellular survival, is expressed in vivo, elicits production of antibodies during infection, and contributes to pathogenicity in a mouse aerosol challenge model. A mutation in the batA gene of wild-type strain ATCC 23344 was found to be particularly attenuating, as BALB/c mice infected with the equivalent of 80 median lethal doses cleared the organism. This finding prompted us to test the hypothesis that vaccination with the batA mutant strain elicits protective immunity against subsequent infection with wild-type bacteria. We discovered that not only does vaccination provide high levels of protection against lethal aerosol challenge with B. mallei ATCC 23344, it also protects against infection with multiple isolates of the closely related organism and causative agent of melioidosis, Burkholderia pseudomallei Passive-transfer experiments also revealed that the protective immunity afforded by vaccination with the batA mutant strain is predominantly mediated by IgG antibodies binding to antigens expressed exclusively in vivo Collectively, our data demonstrate that BatA is a target for developing medical countermeasures and that vaccination with a mutant lacking expression of the protein provides a platform to gain insights regarding mechanisms of protective immunity against B. mallei and B. pseudomallei, including antigen discovery.

Genome-scale analysis of the genes that contribute to Burkholderia pseudomallei biofilm formation identifies a crucial exopolysaccharide biosynthesis gene cluster

Burkholderia pseudomallei, the causative agent of melioidosis, is an important public health threat due to limited therapeutic options for treatment. Efforts to improve therapeutics for B. pseudomallei infections are dependent on the need to understand the role of B. pseudomallei biofilm formation and its contribution to antibiotic tolerance and persistence as these are bacterial traits that prevent effective therapy. In order to reveal the genes that regulate and/or contribute to B. pseudomallei 1026b biofilm formation, we screened a sequence defined two-allele transposon library and identified 118 transposon insertion mutants that were deficient in biofilm formation. These mutants include transposon insertions in genes predicted to encode flagella, fimbriae, transcriptional regulators, polysaccharides, and hypothetical proteins. Polysaccharides are key constituents of biofilms and B. pseudomallei has the capacity to produce a diversity of polysaccharides, thus there is a critical need to link these biosynthetic genes with the polysaccharides they produce to better understand their biological role during infection. An allelic exchange deletion mutant of the entire B. pseudomallei biofilm-associated exopolysaccharide biosynthetic cluster was decreased in biofilm formation and produced a smooth colony morphology suggestive of the loss of exopolysaccharide production. Conversely, deletion of the previously defined capsule I polysaccharide biosynthesis gene cluster increased biofilm formation. Bioinformatics analyses combined with immunoblot analysis and glycosyl composition studies of the partially purified exopolysaccharide indicate that the biofilm-associated exopolysaccharide is neither cepacian nor the previously described acidic exopolysaccharide. The biofilm-associated exopolysaccharide described here is also specific to the B. pseudomallei complex of bacteria. Since this novel exopolysaccharide biosynthesis cluster is retained in B. mallei, it is predicted to have a role in colonization and infection of the host. These findings will facilitate further advances in understanding the pathogenesis of B. pseudomallei and improve diagnostics and therapeutic treatment strategies.

Within-Host Evolution of Burkholderia pseudomallei during Chronic Infection of Seven Australasian Cystic Fibrosis Patients

Cystic fibrosis (CF) is a genetic disorder characterized by progressive lung function decline. CF patients are at an increased risk of respiratory infections, including those by the environmental bacterium Burkholderia pseudomallei, the causative agent of melioidosis. Here, we compared the genomes of B. pseudomallei isolates collected between ~4 and 55 months apart from seven chronically infected CF patients. Overall, the B. pseudomallei strains showed evolutionary patterns similar to those of other chronic infections, including emergence of antibiotic resistance, genome reduction, and deleterious mutations in genes involved in virulence, metabolism, environmental survival, and cell wall components. We documented the first reported B. pseudomallei hypermutators, which were likely caused by defective MutS. Further, our study identified both known and novel molecular mechanisms conferring resistance to three of the five clinically important antibiotics for melioidosis treatment. Our report highlights the exquisite adaptability of microorganisms to long-term persistence in their environment and the ongoing challenges of antibiotic treatment in eradicating pathogens in the CF lung. Convergent evolution with other CF pathogens hints at a degree of predictability in bacterial evolution in the CF lung and potential targeted eradication of chronic CF infections in the future.IMPORTANCEBurkholderia pseudomallei, the causative agent of melioidosis, is an environmental opportunistic bacterium that typically infects immunocompromised people and those with certain risk factors such as cystic fibrosis (CF). Patients with CF tend to develop chronic melioidosis infections, for reasons that are not well understood. This report is the first to describe B. pseudomallei evolution within the CF lung during chronic infection. We show that the pathways by which B. pseudomallei adapts to the CF lung are similar to those seen in better-studied CF pathogens such as Pseudomonas aeruginosa, Staphylococcus aureus, and Burkholderia cepacia complex species. Adaptations include the accumulation of antibiotic resistance, loss of nonessential genes, metabolic alterations, and virulence factor attenuation. Known and novel mechanisms of resistance to three of the five antibiotics used in melioidosis treatment were identified. Similar pathways of evolution in CF pathogens, including B. pseudomallei, provide exciting avenues for more-targeted treatment of chronic, recalcitrant infections.

DBSecSys 2.0: a database of Burkholderia mallei and Burkholderia pseudomallei secretion systems

Background

Burkholderia mallei and B. pseudomallei are the causative agents of glanders and melioidosis, respectively, diseases with high morbidity and mortality rates. B. mallei and B. pseudomallei are closely related genetically; B. mallei evolved from an ancestral strain of B. pseudomallei by genome reduction and adaptation to an obligate intracellular lifestyle. Although these two bacteria cause different diseases, they share multiple virulence factors, including bacterial secretion systems, which represent key components of bacterial pathogenicity. Despite recent progress, the secretion system proteins for B. mallei and B. pseudomallei, their pathogenic mechanisms of action, and host factors are not well characterized.

Results

We previously developed a manually curated database, DBSecSys, of bacterial secretion system proteins for B. mallei. Here, we report an expansion of the database with corresponding information about B. pseudomallei. DBSecSys 2.0 contains comprehensive literature-based and computationally derived information about B. mallei ATCC 23344 and literature-based and computationally derived information about B. pseudomallei K96243. The database contains updated information for 163 B. mallei proteins from the previous database and 61 additional B. mallei proteins, and new information for 281 B. pseudomallei proteins associated with 5 secretion systems, their 1,633 human- and murine-interacting targets, and 2,400 host-B. mallei interactions and 2,286 host-B. pseudomallei interactions. The database also includes information about 13 pathogenic mechanisms of action for B. mallei and B. pseudomallei secretion system proteins inferred from the available literature or computationally. Additionally, DBSecSys 2.0 provides details about 82 virulence attenuation experiments for 52 B. mallei secretion system proteins and 98 virulence attenuation experiments for 61 B. pseudomallei secretion system proteins. We updated the Web interface and data access layer to speed-up users’ search of detailed information for orthologous proteins related to secretion systems of the two pathogens.

Conclusions

The updates of DBSecSys 2.0 provide unique capabilities to access comprehensive information about secretion systems of B. mallei and B. pseudomallei. They enable studies and comparisons of corresponding proteins of these two closely related pathogens and their host-interacting partners.

The database is available at http://dbsecsys.bhsai.org.

The Effects of Signal Erosion and Core Genome Reduction on the Identification of Diagnostic Markers

Whole-genome sequence (WGS) data are commonly used to design diagnostic targets for the identification of bacterial pathogens. To do this effectively, genomics databases must be comprehensive to identify the strict core genome that is specific to the target pathogen. As additional genomes are analyzed, the core genome size is reduced and there is erosion of the target-specific regions due to commonality with related species, potentially resulting in the identification of false positives and/or false negatives.

A comparative analysis of 1,130 Burkholderia genomes identified unique markers for many named species, including the human pathogens B. pseudomallei and B. mallei. Due to core genome reduction and signature erosion, only 38 targets specific to B. pseudomallei/mallei were identified. By using only public genomes, a larger number of markers were identified, due to undersampling, and this larger number represents the potential for false positives. This analysis has implications for the design of diagnostics for other species where the genomic space of the target and/or closely related species is not well defined.