Genomics Reveals a Unique Clone of Burkholderia cenocepacia Harboring an Actively Excising Novel Genomic Island

Impact of template denaturation prior to whole genome amplification on gene detection in high GC-content species, Burkholderia mallei and B. pseudomallei

OBJECTIVE

In this study, we sought to determine the types and prevalence of antimicrobial resistance determinants (ARDs) in Burkholderia spp. strains using the Antimicrobial Resistance Determinant Microarray (ARDM).

RESULTS

Whole genome amplicons from 22 B. mallei (BM) and 37 B. pseudomallei (BP) isolates were tested for > 500 ARDs using ARDM v.3.1. ARDM detected the following Burkholderia spp.-derived genes, aac(6), blaBP/MBL-3, blaABPS, penA-BP, and qacE, in both BM and BP while blaBP/MBL-1, macB, blaOXA-42/43 and penA-BC were observed in BP only. The method of denaturing template for whole genome amplification greatly affected the numbers and types of genes detected by the ARDM. BlaTEM was detected in nearly a third of BM and BP amplicons derived from thermally, but not chemically denatured templates. BlaTEM results were confirmed by PCR, with 81% concordance between methods. Sequences from 414-nt PCR amplicons (13 preparations) were 100% identical to the Klebsiella pneumoniae reference gene. Although blaTEM sequences have been observed in B. glumae, B. cepacia, and other undefined Burkholderia strains, this is the first report of such sequences in BM/BP/B. thailandensis (BT) clade. These results highlight the importance of sample preparation in achieving adequate genome coverage in methods requiring untargeted amplification before analysis.

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.

Genomic features, antimicrobial susceptibility, and epidemiological insights into Burkholderia cenocepacia clonal complex 31 isolates from bloodstream infections in India

Introduction

Burkholderia cepacia complex (Bcc) clonal complex (CC) 31, the predominant lineage causing devastating outbreaks globally, has been a growing concern of infections in non-cystic fibrosis (NCF) patients in India. B. cenocepacia is very challenging to treat owing to its virulence determinants and antibiotic resistance. Improving the management of these infections requires a better knowledge of their resistance patterns and mechanisms.

Methods

Whole-genome sequences of 35 CC31 isolates obtained from patient samples, were analyzed against available 210 CC31 genomes in the NCBI database to glean details of resistance, virulence, mobile elements, and phylogenetic markers to study genomic diversity and evolution of CC31 lineage in India.

Results

Genomic analysis revealed that 35 isolates belonging to CC31 were categorized into 11 sequence types (ST), of which five STs were reported exclusively from India. Phylogenetic analysis classified 245 CC31 isolates into eight distinct clades (I-VIII) and unveiled that NCF isolates are evolving independently from the global cystic fibrosis (CF) isolates forming a distinct clade. The detection rate of seven classes of antibiotic-related genes in 35 isolates was 35 (100%) for tetracyclines, aminoglycosides, and fluoroquinolones; 26 (74.2%) for sulphonamides and phenicols; 7 (20%) for beta-lactamases; and 1 (2.8%) for trimethoprim resistance genes. Additionally, 3 (8.5%) NCF isolates were resistant to disinfecting agents and antiseptics. Antimicrobial susceptibility testing revealed that majority of NCF isolates were resistant to chloramphenicol (77%) and levofloxacin (34%). NCF isolates have a comparable number of virulence genes to CF isolates. A well-studied pathogenicity island of B. cenocepacia, GI11 is present in ST628 and ST709 isolates from the Indian Bcc population. In contrast, genomic island GI15 (highly similar to the island found in B. pseudomallei strain EY1) is exclusively reported in ST839 and ST824 isolates from two different locations in India. Horizontal acquisition of lytic phage ST79 of pathogenic B. pseudomallei is demonstrated in ST628 isolates Bcc1463, Bcc29163, and BccR4654 amongst CC31 lineage.

Discussion

The study reveals a high diversity of CC31 lineages among B. cenocepacia isolates from India. The extensive information from this study will facilitate the development of rapid diagnostic and novel therapeutic approaches to manage B. cenocepacia infections.

Microbial paracetamol degradation involves a high diversity of novel amidase enzyme candidates

Pharmaceuticals are relatively new to nature and often not completely removed in wastewater treatment plants (WWTPs). Consequently, these micropollutants end up in water bodies all around the world posing a great environmental risk. One exception to this recalcitrant conversion is paracetamol, whose full degradation has been linked to several microorganisms. However, the genes and corresponding proteins involved in microbial paracetamol degradation are still elusive. In order to improve our knowledge of the microbial paracetamol degradation pathway, we inoculated a bioreactor with sludge of a hospital WWTP (Pharmafilter, Delft, NL) and fed it with paracetamol as the sole carbon source. Paracetamol was fully degraded without any lag phase and the enriched microbial community was investigated by metagenomic and metatranscriptomic analyses, which demonstrated that the microbial community was very diverse. Dilution and plating on paracetamol-amended agar plates yielded two Pseudomonas sp. isolates: a fast-growing Pseudomonas sp. that degraded 200 mg/L of paracetamol in approximately 10 h while excreting 4-aminophenol, and a slow-growing Pseudomonas sp. that degraded paracetamol without obvious intermediates in more than 90 days. Each Pseudomonas sp. contained a different highly-expressed amidase (31% identity to each other). These amidase genes were not detected in the bioreactor metagenome suggesting that other as-yet uncharacterized amidases may be responsible for the first biodegradation step of paracetamol. Uncharacterized deaminase genes and genes encoding dioxygenase enzymes involved in the catabolism of aromatic compounds and amino acids were the most likely candidates responsible for the degradation of paracetamol intermediates based on their high expression levels in the bioreactor metagenome and the Pseudomonas spp. genomes. Furthermore, cross-feeding between different community members might have occurred to efficiently degrade paracetamol and its intermediates in the bioreactor. This study increases our knowledge about the ongoing microbial evolution towards biodegradation of pharmaceuticals and points to a large diversity of (amidase) enzymes that are likely involved in paracetamol metabolism in WWTPs.

Virulence Potential and Treatment Options of Multidrug-Resistant (MDR) Acinetobacter baumannii

Acinetobacter baumannii is an opportunistic pathogen which is undoubtedly known for a high rate of morbidity and mortality in hospital-acquired infections. A. baumannii causes life-threatening infections, including; ventilator-associated pneumonia (VAP), meningitis, bacteremia, and wound and urinary tract infections (UTI). In 2017, the World Health Organization listed A. baumannii as a priority-1 pathogen. The prevalence of A. baumannii infections and outbreaks emphasizes the direct need for the use of effective therapeutic agents for treating such infections. Available antimicrobials, such as; carbapenems, tigecycline, and colistins have insufficient effectiveness due to the appearance of multidrug-resistant strains, accentuating the need for alternative and novel therapeutic remedies. To understand and overcome this menace, the knowledge of recent discoveries on the virulence factors of A. baumannii is needed. Herein, we summarized the role of various virulence factors, including; outer membrane proteins, efflux pumps, biofilm, penicillin-binding proteins, and siderophores/iron acquisition systems. We reviewed the recent scientific literature on different A. baumannii virulence factors and the effective antimicrobial agents for the treatment and management of bacterial infections.

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.

Genomic analyses of Burkholderia cenocepacia reveal multiple species with differential host-adaptation to plants and humans

Background

Burkholderia cenocepacia is a human opportunistic pathogen causing devastating symptoms in patients suffering from immunodeficiency and cystic fibrosis. Out of the 303 B. cenocepacia strains with available genomes, the large majority were isolated from a clinical context. However, several isolates originate from other environmental sources ranging from aerosols to plant endosphere. Plants can represent reservoirs for human infections as some pathogens can survive and sometimes proliferate in the rhizosphere. We therefore investigated if B. cenocepacia had the same potential.

Results

We selected genome sequences from 31 different strains, representative of the diversity of ecological niches of B. cenocepacia, and conducted comparative genomic analyses in the aim of finding specific niche or host-related genetic determinants. Phylogenetic analyses and whole genome average nucleotide identity suggest that strains, registered as B. cenocepacia, belong to at least two different species. Core-genome analyses show that the clade enriched in environmental isolates lacks multiple key virulence factors, which are conserved in the sister clade where most clinical isolates fall, including the highly virulent ET12 lineage. Similarly, several plant associated genes display an opposite distribution between the two clades. Finally, we suggest that B. cenocepacia underwent a host jump from plants/environment to animals, as supported by the phylogenetic analysis. We eventually propose a name for the new species that lacks several genetic traits involved in human virulence.

Conclusion

Regardless of the method used, our studies resulted in a disunited perspective of the B. cenocepacia species. Strains currently affiliated to this taxon belong to at least two distinct species, one having lost several determining animal virulence factors.

Whole-Genome Comparisons Among the Genus Shewanella Reveal the Enrichment of Genes Encoding Ankyrin-Repeats Containing Proteins in Sponge-Associated Bacteria

The bacterial members of the genus Shewanella are widely distributed and inhabit both freshwater and marine environments. Some members of Shewanella have gained considerable attention due to its ability to survive in redox-stratified environments. However, a gap of knowledge exists on the key genomic features of the sponge-associated Shewanella sp. involving the successful host-bacteria interaction, as sponge-symbiotic Shewanella are largely underrepresented in the public repositories. With the aim of identifying the genomic signatures of sponge-Shewanella association, we generated a high-quality genome data of a sponge-associated, Shewanella sp. OPT22, isolated from the intertidal marine sponge Ophlitaspongia papilla and performed comprehensive comparative analyses of 68 genome strains of the genus Shewanella including two previously reported genomes of sponge-associated bacteria, Shewanella spongiae KCTC 22492 and Shewanella sp. Alg231_23. The 16S rRNA-based phylogenetic reconstruction showed the well-supported affiliation of OPT22 and KCTC 22492 with previously reported sponge-associated bacteria, affirming the “sponge-specific” nature of these two bacterial strains isolated from different marine sponge species from the Atlantic and Pacific (East Sea) Oceans, respectively. The genome comparison of the 68 strains of Shewanella inhabiting different habitats revealed the unusual/previously unreported abundance of genes encoding for ankyrin-repeat containing proteins (ANKs) in the genomes of the two sponge-associated strains, OPT22 (ANKs; n = 45) and KCTC 22492 (ANKs; n = 52), which might be involved in sponge-Shewanella interactions. Focused analyses detected the syntenic organization of the gene cluster encoding major secretion system (type III/IV/VI) components and the presence of effector homologs in OPT22 and KCTC 22492 that seem to play a role in the virulence of the sponge bacteria. The genomic island (GI) of Shewanella sp. OPT22 was identified to localize a gene cluster encoding T4SS components and ANK (n = 1), whereas S. spongiae KCTC 22492 harbored a total of seven ANKs within multiple GIs. GIs may play a pivotal role in the dissemination of symbioses-related genes (ANKs) through the horizontal gene transfer, contributing to the diversification and adaptation of sponge-associated Shewanella. Overall, the genome analyses of Shewanella isolates from marine sponges revealed genomic repertoires that might be involved in establishing successful symbiotic relationships with the sponge hosts.

Insights into the genome diversity and virulence of two clinical isolates of Burkholderia cenocepacia

PURPOSE

Burkholderia cenocepacia is among the most common members of the Burkholderia cepacia complex (Bcc) isolated from patients with cystic fibrosis (CF). The factors triggering the high rates of morbidity and mortality in CF patients are not well elucidated. In this study, we aim to highlight the genome diversity of two clinical isolates of B. cenocepacia through comparative genome analysis.

METHODOLOGY

The repertoire of virulence factors and resistance genes compared to reference strains J2315 and K56-2 was elucidated. The isolates were screened for the presence of phages and insertion sequences. Two methods were combined to obtain an accurate prediction of genomic islands (GIs): the cumulative GC profile and the IslandViewer web tool. To study evolutionary relatedness, whole genome-based single-nucleotide polymorphism (wgSNP) analysis was also performed with 43 publically available strains of the Bcc of various sequence types.Results/Key findings. Genome-based species identification of the two isolates BC-AUH and BC-BMEH confirmed the species as B. cenocepacia. Both belonged to ST-602, a double-locus variant of ST-32 (CC31), genomovar IIIA, and carried a large number of antibiotic resistance genes. Eighteen GIs were predicted in BC-AUH and BC-BMEH, occupying 9.3 and 6.1 % of the respective genomes. Comparison to J2315 revealed 89 and 85 genes unique to BC-BMEH and BC-AUH, respectively. Additionally, 1823 intergenic SNPs were detected between BC-BMEH and BC-AUH.

CONCLUSION

This study mapped existing genetic variations in B. cenocepacia associated with notorious outcomes in CF patients, and the data obtained provide comprehensive, genome-inferred insights and multifactorial examination of an important human pathogen.

Burkholderia cenocepacia Infections in Cystic Fibrosis Patients: Drug Resistance and Therapeutic Approaches

Burkholderia cenocepacia is an opportunistic pathogen particularly dangerous for cystic fibrosis (CF) patients. It can cause a severe decline in CF lung function possibly developing into a life-threatening systemic infection known as cepacia syndrome. Antibiotic resistance and presence of numerous virulence determinants in the genome make B. cenocepacia extremely difficult to treat. Better understanding of its resistance profiles and mechanisms is crucial to improve management of these infections. Here, we present the clinical distribution of B. cenocepacia described in the last 6 years and methods for identification and classification of epidemic strains. We also detail new antibiotics, clinical trials, and alternative approaches reported in the literature in the last 5 years to tackle B. cenocepacia resistance issue. All together these findings point out the urgent need of new and alternative therapies to improve CF patients’ life expectancy.

PacBio But Not Illumina Technology Can Achieve Fast, Accurate and Complete Closure of the High GC, Complex Burkholderia pseudomallei Two-Chromosome Genome

Although PacBio third-generation sequencers have improved the read lengths of genome sequencing which facilitates the assembly of complete genomes, no study has reported success in using PacBio data alone to completely sequence a two-chromosome bacterial genome from a single library in a single run. Previous studies using earlier versions of sequencing chemistries have at most been able to finish bacterial genomes containing only one chromosome with de novo assembly. In this study, we compared the robustness of PacBio RS II, using one SMRT cell and the latest P6-C4 chemistry, with Illumina HiSeq 1500 in sequencing the genome of Burkholderia pseudomallei, a bacterium which contains two large circular chromosomes, very high G+C content of 68–69%, highly repetitive regions and substantial genomic diversity, and represents one of the largest and most complex bacterial genomes sequenced, using a reference genome generated by hybrid assembly using PacBio and Illumina datasets with subsequent manual validation. Results showed that PacBio data with de novo assembly, but not Illumina, was able to completely sequence the B. pseudomallei genome without any gaps or mis-assemblies. The two large contigs of the PacBio assembly aligned unambiguously to the reference genome, sharing >99.9% nucleotide identities. Conversely, Illumina data assembled using three different assemblers resulted in fragmented assemblies (201–366 contigs), sharing only 92.2–100% and 92.0–100% nucleotide identities to chromosomes I and II reference sequences, respectively, with no indication that the B. pseudomallei genome consisted of two chromosomes with four copies of ribosomal operons. Among all assemblies, the PacBio assembly recovered the highest number of core and virulence proteins, and housekeeping genes based on whole-genome multilocus sequence typing (wgMLST). Most notably, assembly solely based on PacBio outperformed even hybrid assembly using both PacBio and Illumina datasets. Hybrid approach generated only 74 contigs, while the PacBio data alone with de novo assembly achieved complete closure of the two-chromosome B. pseudomallei genome without additional costly bench work and further sequencing. PacBio RS II using P6-C4 chemistry is highly robust and cost-effective and should be the platform of choice in sequencing bacterial genomes, particularly for those that are well-known to be difficult-to-sequence.