Insertion sequence evolutionary patterns highlight convergent genetic inactivations and recent genomic island acquisitions among epidemic Burkholderia cenocepacia

Clonal diversity, antimicrobial resistance, and genome features among nonfermenting gram-negative bacteria isolated from patients with cystic fibrosis in Russia

Nonfermenting gram-negative (NFGN) bacteria were isolated from cystic fibrosis (CF) patients and subjected to susceptibility testing and whole-genome sequencing. Among 170 enrolled CF patients, 112 (65.9%) were colonized with at least 1 key NFGN species. The species-specific infection rate was highest for Pseudomonas aeruginosa (40.6%) followed by Stenotrophomonas maltophilia (14.1%), Achromobacter spp. (9.4%), and Burkholderia cepacia complex (Bcc, 8.2%) demonstrating a significant age-dependent increase for P. aeruginosa and Achromobacter spp., but not for S. maltophilia or Bcc. P. aeruginosa sequence types (STs) related to high-risk epidemic and global CF clones were carried by 12 (7.1%) and 13 (7.6%) patients, respectively. In total, 47% NFGN isolates, predominantly P. aeruginosa, harbored at least 1 plasmid-borne resistance gene; 5 ST235 isolates carried blaVIM2. Pathogenicity island-borne virulence genes were harbored by 9% NFGN isolates. These findings in conjunction with frequent early colonization by Bcc raised serious concerns regarding infection control in Russian CF centers.

Genomic analysis of Burkholderia cenocepacia isolated from a liver abscess in a patient with cystic fibrosis

Burkholderia cenocepacia complex is associated with high transmissibility, virulence, and poor prognosis in cystic fibrosis (CF) patients. However, extrapulmonary infections are rare. We investigated the genome of a B. cenocepacia IIIA isolated from a liver abscess in a Brazilian CF patient and compared it to strain J2315. The whole genome was sequenced, and contigs were annotated by Rapid Annotation using Subsystem Technology. The Pathosystems Resource Integration Center was used to map antimicrobial and virulence genes. The genomic island (GIs) analysis was performed using two prediction methods, and the presence of putative plasmids and insertion sequences (ISs) was investigated. The isolate was confirmed as B. cenocepacia IIIA to ST-28 (ET12 lineage). A total of 64 genes for antimicrobial resistance and 47 genes related to virulence were identified. Among the virulence factors, there was a predominance of factors related to the invasion mechanism, to the flagellar biosynthesis protein, and to the RNA polymerase sigma factor for flagellar operon (cdpA). Two IS families (IS3 and IS5) and only one plasmid were found. On average 56 GIs were predicted by at least one of the methods applied. Comparative analysis showed resistance mechanisms and virulence factors revealing invasive determinants used by B. cenocepacia IIIA (ET12) in the process of disease spread to other infection sites (extrapulmonary) of highly virulent strains in CF patients.

Comparative Evolutionary Patterns of Burkholderia cenocepacia and B. multivorans During Chronic Co-infection of a Cystic Fibrosis Patient Lung

During chronic respiratory infections of cystic fibrosis (CF) patients, bacteria adaptively evolve in response to the nutritional and immune environment as well as influence other infecting microbes. The present study was designed to gain insights into the genetic mechanisms underlying adaptation and diversification by the two most prevalent pathogenic species of the Burkholderia cepacia complex (Bcc), B. cenocepacia and B. multivorans. Herein, we study the evolution of both of these species during coinfection of a CF patient for 4.4 years using genome sequences of 9 B. multivorans and 11 B. cenocepacia. This co-infection spanned at least 3 years following initial infection by B. multivorans and ultimately ended in the patient's death by cepacia syndrome. Both species acquired several mutations with accumulation rates of 2.08 (B. cenocepacia) and 2.27 (B. multivorans) SNPs/year. Many of the mutated genes are associated with oxidative stress response, transition metal metabolism, defense mechanisms against antibiotics, and other metabolic alterations consistent with the idea that positive selection might be driven by the action of the host immune system, antibiotic therapy and low oxygen and iron concentrations. Two orthologous genes shared by B. cenocepacia and B. multivorans were found to be under strong selection and accumulated mutations associated with lineage diversification. One gene encodes a nucleotide sugar dehydratase involved in lipopolysaccharide O-antigen (OAg) biosynthesis (wbiI). The other gene encodes a putative two-component regulatory sensor kinase protein required to sense and adapt to oxidative- and heavy metal- inducing stresses. This study contributes to understanding of shared and species-specific evolutionary patterns of B. cenocepacia and B. multivorans evolving in the same CF lung environment.

Comparative analysis of the Burkholderia cenocepacia K56-2 essential genome reveals cell envelope functions that are uniquely required for survival in species of the genus Burkholderia

Burkholderia cenocepacia K56-2 belongs to the Burkholderia cepacia complex, a group of Gram-negative opportunistic pathogens that have large and dynamic genomes. In this work, we identified the essential genome of B. cenocepacia K56-2 using high-density transposon mutagenesis and insertion site sequencing (Tn-seq circle). We constructed a library of one million transposon mutants and identified the transposon insertions at an average of one insertion per 27 bp. The probability of gene essentiality was determined by comparing of the insertion density per gene with the variance of neutral datasets generated by Monte Carlo simulations. Five hundred and eight genes were not significantly disrupted, suggesting that these genes are essential for survival in rich, undefined medium. Comparison of the B. cenocepacia K56-2 essential genome with that of the closely related B. cenocepacia J2315 revealed partial overlapping, suggesting that some essential genes are strain-specific. Furthermore, 158 essential genes were conserved in B. cenocepacia and two species belonging to the Burkholderia pseudomallei complex, B. pseudomallei K96243 and Burkholderia thailandensis E264. Porins, including OpcC, a lysophospholipid transporter, LplT, and a protein involved in the modification of lipid A with aminoarabinose were found to be essential in Burkholderia genomes but not in other bacterial essential genomes identified so far. Our results highlight the existence of cell envelope processes that are uniquely essential in species of the genus Burkholderia for which the essential genomes have been identified by Tn-seq.

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.

Phylogeographic diversity and mosaicism of the Helicobacter pylori tfs integrative and conjugative elements

Background

The genome of the gastric pathogen Helicobacter pylori is characterised by considerable variation of both gene sequence and content, much of which is contained within three large genomic islands comprising the cag pathogenicity island (cagPAI) and two mobile integrative and conjugative elements (ICEs) termed tfs3 and tfs4. All three islands are implicated as virulence factors, although whereas the cagPAI is well characterised, understanding of how the tfs elements influence H. pylori interactions with different human hosts is significantly confounded by limited definition of their distribution, diversity and structural representation in the global H. pylori population.

Results

To gain a global perspective of tfs ICE population dynamics we established a bioinformatics workflow to extract and precisely define the full tfs pan-gene content contained within a global collection of 221 draft and complete H. pylori genome sequences. Complete (ca. 35-55kbp) and remnant tfs ICE clusters were reconstructed from a dataset comprising > 12,000 genes, from which orthologous gene complements and distinct alleles descriptive of different tfs ICE types were defined and classified in comparative analyses. The genetic variation within defined ICE modular segments was subsequently used to provide a complete description of tfs ICE diversity and a comprehensive assessment of their phylogeographic context. Our further examination of the apparent ICE modular types identified an ancient and complex history of ICE residence, mobility and interaction within particular H. pylori phylogeographic lineages and further, provided evidence of both contemporary inter-lineage and inter-species ICE transfer and displacement.

Conclusions

Our collective results establish a clear view of tfs ICE diversity and phylogeographic representation in the global H. pylori population, and provide a robust contextual framework for elucidating the functional role of the tfs ICEs particularly as it relates to the risk of gastric disease associated with different tfs ICE genotypes.

Electronic supplementary material

The online version of this article (10.1186/s13100-018-0109-4) contains supplementary material, which is available to authorized users.

What matters in chronic Burkholderia cenocepacia infection in cystic fibrosis: Insights from comparative genomics

Burkholderia cenocepacia causes severe pulmonary infections in cystic fibrosis (CF) patients. Since the bacterium is virtually untreatable by antibiotics, chronic infections persist for years and might develop into fatal septic pneumonia (cepacia syndrome, CS). To devise new strategies to combat chronic B. cenocepacia infections, it is essential to obtain comprehensive knowledge about their pathogenesis. We conducted a comparative genomic analysis of 32 Czech isolates of epidemic clone B. cenocepacia ST32 isolated from various stages of chronic infection in 8 CF patients. High numbers of large-scale deletions were found to occur during chronic infection, affecting preferentially genomic islands and nonessential replicons. Recombination between insertion sequences (IS) was inferred as the mechanism behind deletion formation; the most numerous IS group was specific for the ST32 clone and has undergone transposition burst since its divergence. Genes functionally related to transition metal metabolism were identified as hotspots for deletions and IS insertions. This functional category was also represented among genes where nonsynonymous point mutations and indels occurred parallelly among patients. Another category exhibiting parallel mutations was oxidative stress protection; mutations in catalase KatG resulted in impaired detoxification of hydrogen peroxide. Deep sequencing revealed substantial polymorphism in genes of both categories within the sputum B. cenocepacia ST32 populations, indicating extensive adaptive evolution. Neither oxidative stress response nor transition metal metabolism genes were previously reported to undergo parallel evolution during chronic CF infection. Mutations in katG and copper metabolism genes were overrepresented in patients where chronic infection developed into CS. Among professional phagocytes, macrophages use both hydrogen peroxide and copper for their bactericidal activity; our results thus tentatively point to macrophages as suspects in pathogenesis towards the fatal CS.

The large Burkholderia cenocepacia populations which persist in cystic fibrosis lungs during many years of chronic infections have an inherent potential for adaptive evolution. The results provided by comparative genomics are key in understanding the processes involved. Mutational events which have taken place allow us to deductively reconstruct the history of chronic infection and to identify driving forces acting upon the bacteria. Beyond the conventional point mutation analysis of next generation sequencing data, we observed interesting phenomena such as large deletions and transposable element movement which represent another facet of adaptive evolution of B. cenocepacia during chronic infection. We also found, unexpectedly, that adaptive evolution in B. cenocepacia strain ST32 affects a set of genes conspicuously different from related species B. dolosa; these appear to be linked to host immune response. Our study provides clues to the complex puzzle of chronic B. cenocepacia infection establishment, persistence and outcome in cystic fibrosis.

Identification and analysis of genomic islands in Burkholderia cenocepacia AU 1054 with emphasis on pathogenicity islands

BACKGROUND

Genomic islands (GIs) are genomic regions that reveal evidence of horizontal DNA transfer. They can code for many functions and may augment a bacterium's adaptation to its host or environment. GIs have been identified in strain J2315 of Burkholderia cenocepacia, whereas in strain AU 1054 there has been no published works on such regions according to our text mining and keyword search in Medline.

RESULTS

In this study, we identified 21 GIs in AU 1054 by combining two computational tools. Feature analyses suggested that the predictions are highly reliable and hence illustrated the advantage of joint predictions by two independent methods. Based on putative virulence factors, four GIs were further identified as pathogenicity islands (PAIs). Through experiments of gene deletion mutants in live bacteria, two putative PAIs were confirmed, and the virulence factors involved were identified as lipA and copR. The importance of the genes lipA (from PAI 1) and copR (from PAI 2) for bacterial invasion and replication indicates that they are required for the invasive properties of B. cenocepacia and may function as virulence determinants for bacterial pathogenesis and host infection.

CONCLUSIONS

This approach of in silico prediction of GIs and subsequent identification of potential virulence factors in the putative island regions with final validation using wet experiments could be used as an effective strategy to rapidly discover novel virulence factors in other bacterial species and strains.

Structure and fate of a Pseudomonas aeruginosa population originating from a combined sewer and colonizing a wastewater treatment lagoon

The efficacy of a wastewater treatment lagoon (WWTL) at preventing the spread of Pseudomonas aeruginosa into natural aquatic habitats was investigated. A WWTL and its connected combined sewer and brook were exhaustively sampled. Physico-chemical analyses showed a stratification of the first pond according to pH, temperature and oxygen content. The P. aeruginosa counts partially matched this stratification with higher values among the bottom anaerobic waters of the first half of this pond. Genotyping of 494 WWTL P. aeruginosa strains was performed and led to the definition of 85 lineages. Dominant lineages were observed, with some being found all over the WWTL including the connected brook. IS5 was used as an indicator of genomic changes, and 1 to 12 elements were detected among 16 % of the strains. IS-driven lasR (genetic regulator) disruptions were detected among nine strains that were not part of the dominant lineages. These insertional mutants did not show significant elastase activities but showed better growth than the PAO1 reference strain in WWTL waters. Differences in growth patterns were related to a better survival of these mutants at an alkaline pH and a better ability at using some C-sources such as alanine. The opportunistic colonization of a WWTL by P. aeruginosa can involve several metabolic strategies which appeared lineage specific. Some clones appeared more successful than others at disseminating from a combined sewer toward the overflow of a WWTL.

Niches, population structure and genome reduction in Ochrobactrum intermedium: clues to technology-driven emergence of pathogens

Ochrobactrum intermedium is considered as an emerging human environmental opportunistic pathogen with mild virulence. The distribution of isolates and sequences described in literature and databases showed frequent association with human beings and polluted environments. As population structures are related to bacterial lifestyles, we investigated by multi-locus approach the genetic structure of a population of 65 isolates representative of the known natural distribution of O. intermedium. The population was further surveyed for genome dynamics using pulsed-field gel electrophoresis and genomics. The population displayed a clonal epidemic structure with events of recombination that occurred mainly in clonal complexes. Concerning biogeography, clones were shared by human and environments and were both cosmopolitan and local. The main cosmopolitan clone was genetically and genomically stable, and grouped isolates that all harbored an atypical insertion in the rrs. Ubiquitism and stability of this major clone suggested a clonal succes in a particular niche. Events of genomic reduction were detected in the population and the deleted genomic content was described for one isolate. O. intermedium displayed allopatric characters associated to a tendancy of genome reduction suggesting a specialization process. Considering its relatedness with Brucella, this specialization might be a commitment toward pathogenic life-style that could be driven by technological selective pressure related medical and industrial technologies.

The Nocardia cyriacigeorgica GUH-2 genome shows ongoing adaptation of an environmental Actinobacteria to a pathogen's lifestyle

Background

Nocardia cyriacigeorgica is recognized as one of the most prevalent etiological agents of human nocardiosis. Human exposure to these Actinobacteria stems from direct contact with contaminated environmental matrices. The full genome sequence of N. cyriacigeorgica strain GUH-2 was studied to infer major trends in its evolution, including the acquisition of novel genetic elements that could explain its ability to thrive in multiple habitats.

Results

N. cyriacigeorgica strain GUH-2 genome size is 6.19 Mb-long, 82.7% of its CDS have homologs in at least another actinobacterial genome, and 74.5% of these are found in N. farcinica. Among N. cyriacigeorgica specific CDS, some are likely implicated in niche specialization such as those involved in denitrification and RuBisCO production, and are found in regions of genomic plasticity (RGP). Overall, 22 RGP were identified in this genome, representing 11.4% of its content. Some of these RGP encode a recombinase and IS elements which are indicative of genomic instability. CDS playing part in virulence were identified in this genome such as those involved in mammalian cell entry or encoding a superoxide dismutase. CDS encoding non ribosomal peptide synthetases (NRPS) and polyketide synthases (PKS) were identified, with some being likely involved in the synthesis of siderophores and toxins. COG analyses showed this genome to have an organization similar to environmental Actinobacteria.

Conclusion

N. cyriacigeorgica GUH-2 genome shows features suggesting a diversification from an ancestral saprophytic state. GUH-2 ability at acquiring foreign DNA was found significant and to have led to functional changes likely beneficial for its environmental cycle and opportunistic colonization of a human host.