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

Abstract

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.

The Ever-Expanding Pseudomonas Genus: Description of 43 New Species and Partition of the Pseudomonas putida Group

“…We applied the same methodology to an expanded set of Pseudomonas isolates and identified 17 additional new species. To provide high-quality genomes for the type strains of these 43 new species, we combined Illumina and Nanopore sequencing [ 32 ]. An overview of the different sequencing methodologies used for the entire set of strains is shown in Table S2 .…”

Genome-driven elucidation of phage-host interplay and impact of phage resistance evolution on bacterial fitness

“…The low frequency of mutations in LPS-related genes, detected in our study, indicates the prevention of specific features crucial for clinically evolved strain. This phenomenon could be further investigated in a wider context of bacterial virulence and by complementing the list of genetic changes with those from larger structural variation analyses (transposon activity, prophage mobility, mobile islands) based on long reads sequencing data [48].…”

Bacteriophage Therapy for Difficult-to-Treat Infections: The Implementation of a Multidisciplinary Phage Task Force (The PHAGEFORCE Study Protocol)

“…Antibiograms will be compared to the baseline antibiogram to determine if the susceptibility to antibiotics has changed and look for synergistic interactions [ 17 ]. Furthermore, bacterial isolates will be fully characterized by whole genome sequencing, using a combination of Illumina and Nanopore sequencing [ 18 ] to obtain optimal genome assemblies. This will be crucial to improve the functional annotation of the genomes.…”