Comparative genomics of Burkholderia multivorans, a ubiquitous pathogen with a highly conserved genomic structure

Effects of simulated low-temperature thermal remediation on the microbial community of a tropical creosote contaminated soil

In the search for more sustainable remediation strategies for PAH-contaminated soils, an integrated application of thermal remediation and bioremediation (TEB) may allow the use of less impacting temperatures by associating heating to biological degradation. However, the influence of heating on soil microbiota remains poorly understood, especially in soils from tropical regions. This work investigated the effects of low-temperature heating on creosote-contaminated soil bacteria. We used culture-dependent and 16 S rRNA sequencing methods to compare the microbial community of soil samples heated to 60 and 100 oC for 1 h in microcosms. Heating to 60 °C reduced the density of cultivable heterotrophic bacteria compared to control soil (p < 0.05), and exposure to 100 °C inactivated the viable heterotrophic community. Burkholderia-Caballeronia-Paraburkholderia (BCP) group and Sphingobium were the predominant genera. Temperature and incubation time affected the Bray-Curtis dissimilarity index (p < 0.05). At 60 °C and 30 days incubation, the relative abundance of Sphingobium decreased and BCP increased dominance. The network of heated soil after 30 days of incubation showed fewer nodes and edges but maintained its density and complexity. Both main genera are associated with PAH degradation, suggesting functional redundancy and a likely potential of soil microbiota to maintain biodegradation ability after exposure to higher temperatures. We concluded that TEB can be considered as a potential strategy to bioremediate creosote-contaminated soils, allowing biodegradation in temperature ranges where thermal remediation does not completely remove contaminants. However, we recommend further research to determine degradation rates with this technology.

Recombineering enables genome mining of novel siderophores in a non-model Burkholderiales strain

Iron is essential for bacterial survival, and most bacteria capture iron by producing siderophores. Burkholderiales bacteria produce various types of bioactive secondary metabolites, such as ornibactin and malleobactin siderophores. In this study, the genome analysis of Burkholderiales genomes showed a putative novel siderophore gene cluster crb, which is highly similar to the ornibactin and malleobactin gene clusters but does not have pvdF, a gene encoding a formyltransferase for N-δ‑hydroxy-ornithine formylation. Establishing the bacteriophage recombinase Redγ-Redαβ7029 mediated genome editing system in a non-model Burkholderiales strain Paraburkholderia caribensis CICC 10960 allowed the rapid identification of the products of crb gene cluster, caribactins A-F (1-6). Caribactins contain a special amino acid residue N-δ‑hydroxy-N-δ-acetylornithine (haOrn), which differs from the counterpart N-δ‑hydroxy-N-δ-formylornithine (hfOrn) in ornibactin and malleobactin, owing to the absence of pvdF. Gene inactivation showed that the acetylation of hOrn is catalyzed by CrbK, whose homologs probably not be involved in the biosynthesis of ornibactin and malleobactin, showing possible evolutionary clues of these siderophore biosynthetic pathways from different genera. Caribactins promote biofilm production and enhance swarming and swimming abilities, suggesting that they may play crucial roles in biofilm formation. This study also revealed that recombineering has the capability to mine novel secondary metabolites from non-model Burkholderiales species.

Identification of two distinct phylogenomic lineages and model strains for the understudied cystic fibrosis lung pathogen Burkholderia multivorans

Burkholderia multivorans is the dominant Burkholderia pathogen recovered from lung infection in people with cystic fibrosis. However, as an understudied pathogen there are knowledge gaps in relation to its population biology, phenotypic traits and useful model strains. A phylogenomic study of B. multivorans was undertaken using a total of 283 genomes, of which 73 were sequenced and 49 phenotypically characterized as part of this study. Average nucleotide identity analysis (ANI) and phylogenetic alignment of core genes demonstrated that the B. multivorans population separated into two distinct evolutionary clades, defined as lineage 1 (n=58 genomes) and lineage 2 (n=221 genomes). To examine the population biology of B. multivorans, a representative subgroup of 77 B. multivorans genomes (28 from the reference databases and the 49 novel short-read genome sequences) were selected based on multilocus sequence typing (MLST), isolation source and phylogenetic placement criteria. Comparative genomics was used to identify B. multivorans lineage-specific genes - ghrB_1 in lineage 1 and glnM_2 in lineage 2 - and diagnostic PCRs targeting them were successfully developed. Phenotypic analysis of 49 representative B. multivorans strains showed considerable inter-strain variance, but the majority of the isolates tested were motile and capable of biofilm formation. A striking absence of B. multivorans protease activity in vitro was observed, but no lineage-specific phenotypic differences were demonstrated. Using phylogenomic and phenotypic criteria, three model B. multivorans CF strains were identified, BCC0084 (lineage 1), BCC1272 (lineage 2a) and BCC0033 lineage 2b, and their complete genome sequences determined. B. multivorans CF strains BCC0033 and BCC0084, and the environmental reference strain, ATCC 17616, were all capable of short-term survival within a murine lung infection model. By mapping the population biology, identifying lineage-specific PCRs and model strains, we provide much needed baseline resources for future studies of B. multivorans.

Comprehensive genome analysis of Burkholderia contaminans SK875, a quorum-sensing strain isolated from the swine

The Burkholderia cepacia complex (BCC) is a Gram-negative bacterial, including Burkholderia contaminans species. Although the plain Burkholderia is pervasive from taxonomic and genetic perspectives, a common characteristic is that they may use the quorum-sensing (QS) system. In our previous study, we generated the complete genome sequence of Burkholderia contaminans SK875 isolated from the respiratory tract. To our knowledge, this is the first study to report functional genomic features of B. contaminans SK875 for understanding the pathogenic characteristics. In addition, comparative genomic analysis for five B. contaminans genomes was performed to provide comprehensive information on the disease potential of B. contaminans species. Analysis of average nucleotide identity (ANI) showed that the genome has high similarity (> 96%) with other B. contaminans strains. Five B. contaminans genomes yielded a pangenome of 8832 coding genes, a core genome of 5452 genes, the accessory genome of 2128 genes, and a unique genome of 1252 genes. The 186 genes were specific to B. contaminans SK875, including toxin higB-2, oxygen-dependent choline dehydrogenase, and hypothetical proteins. Genotypic analysis of the antimicrobial resistance of B. contaminans SK875 verified resistance to tetracycline, fluoroquinolone, and aminoglycoside. Compared with the virulence factor database, we identified 79 promising virulence genes such as adhesion system, invasions, antiphagocytic, and secretion systems. Moreover, 45 genes of 57 QS-related genes that were identified in B. contaminans SK875 indicated high sequence homology with other B. contaminans strains. Our results will help to gain insight into virulence, antibiotic resistance, and quorum sensing for B. contaminans species.

Examining the genomic features of human and plant-associated Burkholderia strains

Humans and plants have evolved in the near omnipresence of a microbial milieu, and the factors that govern host-microbe interactions continue to require scientific exploration. To better understand if and to what degree patterns between microbial genomic features and host association (i.e., human and plant) exist, I analyzed the genomes of select Burkholderia strains-a bacterial genus comprised of both human and plant-associated strains-that were isolated from either humans or plants. To this end, I uncovered host-specific, genomic patterns related to metabolic pathway potentials in addition to convergent features that may be related to pathogenic overlap between hosts. Together, these findings detail the genomic associations of human and plant-associated Burkholderia strains and provide a framework for future investigations that seek to link host-host transmission potentials.

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

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

Burkholderia: An Untapped but Promising Bacterial Genus for the Conversion of Aromatic Compounds

Burkholderia, a bacterial genus comprising more than 120 species, is typically reported to inhabit soil and water environments. These Gram-negative bacteria harbor a variety of aromatic catabolic pathways and are thus potential organisms for bioremediation of sites contaminated with aromatic pollutants. However, there are still substantial gaps in our knowledge of these catabolic processes that must be filled before these pathways and organisms can be harnessed for biotechnological applications. This review presents recent discoveries on the catabolism of monoaromatic and polycyclic aromatic hydrocarbons, as well as of heterocyclic compounds, by a diversity of Burkholderia strains. We also present a perspective on the beneficial features of Burkholderia spp. and future directions for their potential utilization in the bioremediation and bioconversion of aromatic compounds.

Versatile and Dynamic Symbioses Between Insects and Burkholderia Bacteria

Symbiotic associations with microorganisms represent major sources of ecological and evolutionary innovations in insects. Multiple insect taxa engage in symbioses with bacteria of the genus Burkholderia, a diverse group that is widespread across different environments and whose members can be mutualistic or pathogenic to plants, fungi, and animals. Burkholderia symbionts provide nutritional benefits and resistance against insecticides to stinkbugs, defend Lagria beetle eggs against pathogenic fungi, and may be involved in nitrogen metabolism in ants. In contrast to many other insect symbioses, the known associations with Burkholderia are characterized by environmental symbiont acquisition or mixed-mode transmission, resulting in interesting ecological and evolutionary dynamics of symbiont strain composition. Insect-Burkholderia symbioses present valuable model systems from which to derive insights into general principles governing symbiotic interactions because they are often experimentally and genetically tractable and span a large fraction of the diversity of functions, localizations, and transmission routes represented in insect symbioses.

Burkholderia cepacia Complex Species Differ in the Frequency of Variation of the Lipopolysaccharide O-Antigen Expression During Cystic Fibrosis Chronic Respiratory Infection

Burkholderia cepacia complex (Bcc) bacteria can adapt to the lung environment of cystic fibrosis (CF) patients resulting in the emergence of a very difficult to eradicate heterogeneous population leading to chronic infections associated with rapid lung function loss and increased mortality. Among the important phenotypic modifications is the variation of the lipopolysaccharide (LPS) structure at level of the O-antigen (OAg) presence, influencing adherence, colonization and the ability to evade the host defense mechanisms. The present study was performed to understand whether the loss of OAg expression during CF infection can be considered a general phenomenon in different Bcc species favoring its chronicity. In fact, it is still not clear why different Bcc species/strains differ in their ability to persist in the CF lung and pathogenic potential. The systematic two-decade-retrospective-longitudinal-screening conducted covered 357 isolates retrieved from 19 chronically infected patients receiving care at a central hospital in Lisbon. The study involved 21 Bcc strains of six/seven Bcc species/lineages, frequently or rarely isolated from CF patients worldwide. Different strains/clonal variants obtained during infection gave rise to characteristic OAg-banding patterns. The two most prevalent and feared species, B. cenocepacia and B. multivorans, showed a tendency to lose the OAg along chronic infection. B. cenocepacia recA lineage IIIA strains known to lead to particularly destructive infections exhibit the most frequent OAg loss, compared with lineage IIIB. The switch frequency increased with the duration of infection and the level of lung function deterioration. For the first time, it is shown that the rarely found B. cepacia and B. contaminans, whose representation in the cohort of patients examined is abnormally high, keep the OAg even during 10- or 15-year infections. Data from co-infections with different Bcc species reinforced these conclusions. Concerning the two other rarely found species examined, B. stabilis exhibited a stable OAg expression phenotype over the infection period while for the single clone of the more distantly related B. dolosa species, the OAg-chain was absent from the beginning of the 5.5-year infection until the patient dead. This work reinforces the relevance attributed to the OAg-expression switch suggesting marked differences in the various Bcc species.

Comparative Genomics of Environmental and Clinical Burkholderia cenocepacia Strains Closely Related to the Highly Transmissible Epidemic ET12 Lineage

The Burkholderia cenocepacia epidemic ET12 lineage belongs to the genomovar IIIA including the reference strain J2315, a highly transmissible epidemic B. cenocepacia lineage. Members of this lineage are able to cause lung infections in immunocompromised and cystic fibrosis patients. In this study, we describe the genome of F01, an environmental B. cenocepacia strain isolated from soil in Burkina Faso that is, to our knowledge, the most closely related strain to this epidemic lineage. A comparative genomic analysis was performed on this new isolate, in association with five clinical and one environmental B. cenocepacia strains whose genomes were previously sequenced. Antibiotic resistances, virulence phenotype, and genomic contents were compared and discussed with an emphasis on virulent and antibiotic determinants. Surprisingly, no significant differences in antibiotic resistance and virulence were found between clinical and environmental strains, while the most important genomic differences were related to the number of prophages identified in their genomes. The ET12 lineage strains showed a noticeable greater number of prophages (partial or full-length), especially compared to the phylogenetically related environmental F01 strain (i.e., 5–6 and 3 prophages, respectively). Data obtained suggest possible involvements of prophages in the clinical success of opportunistic pathogens.

Phylogeographic, genomic, and meropenem susceptibility analysis of Burkholderia ubonensis

The bacterium Burkholderia ubonensis is commonly co-isolated from environmental specimens harbouring the melioidosis pathogen, Burkholderia pseudomallei. B. ubonensis has been reported in northern Australia and Thailand but not North America, suggesting similar geographic distribution to B. pseudomallei. Unlike most other Burkholderia cepacia complex (Bcc) species, B. ubonensis is considered non-pathogenic, although its virulence potential has not been tested. Antibiotic resistance in B. ubonensis, particularly towards drugs used to treat the most severe B. pseudomallei infections, has also been poorly characterised. This study examined the population biology of B. ubonensis, and includes the first reported isolates from the Caribbean. Phylogenomic analysis of 264 B. ubonensis genomes identified distinct clades that corresponded with geographic origin, similar to B. pseudomallei. A small proportion (4%) of strains lacked the 920kb chromosome III replicon, with discordance of presence/absence amongst genetically highly related strains, demonstrating that the third chromosome of B. ubonensis, like other Bcc species, probably encodes for a nonessential pC3 megaplasmid. Multilocus sequence typing using the B. pseudomallei scheme revealed that one-third of strains lack the "housekeeping" narK locus. In comparison, all strains could be genotyped using the Bcc scheme. Several strains possessed high-level meropenem resistance (≥32 μg/mL), a concern due to potential transmission of this phenotype to B. pseudomallei. In silico analysis uncovered a high degree of heterogeneity among the lipopolysaccharide O-antigen cluster loci, with at least 35 different variants identified. Finally, we show that Asian B. ubonensis isolate RF23-BP41 is avirulent in the BALB/c mouse model via a subcutaneous route of infection. Our results provide several new insights into the biology of this understudied species.