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

Antibacterial Activity of a Natural Clay Mineral against Burkholderia cepacia Complex and Other Bacterial Pathogens Isolated from People with Cystic Fibrosis

There is an impending crisis in healthcare brought about by a new era of untreatable infections caused by bacteria resistant to all available antibiotics. Thus, there is an urgent need to identify novel antimicrobial agents to counter the continuing threat posed by formerly treatable infections. We previously reported that a natural mineral clay known as Kisameet clay (KC) is a potent inhibitor of the organisms responsible for acute infections. Chronic bacterial infections present another major challenge to treatment by antimicrobials, due to their prolonged nature, which results in repeated exposure to antibiotics and a constant selection for antimicrobial resistance. A prime example is bacteria belonging to the Burkholderia cepacia complex (Bcc), which particularly causes some of the most serious chronic lung infections in patients with cystic fibrosis (CF) associated with unpredictable clinical outcomes, poor prognosis, and high mortality rates. Eradication of these organisms from CF patients with limited effective antimicrobial options is a major challenge. Novel therapeutic approaches are urgently required. Here, we report the in vitro antibacterial activity of KC aqueous suspensions (1-10% w/v) and its aqueous extract (L100) against a collection of extensively and multi-drug resistant clinical isolates of Bcc, Pseudomonas aeruginosa, and Stenotrophomonas maltophilia isolated from patients with CF. These findings present a potential novel therapy for further investigation in the clinic.

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.

Adaptation and Evolution of Pathogens in the Cystic Fibrosis Lung

As opposed to acute respiratory infections, the persistent bacterial infections of the lung that characterize cystic fibrosis (CF) provide ample time for bacteria to evolve and adapt. The process of adaptation is recorded in mutations that accumulate over time in the genomes of the infecting bacteria. Some of these mutations lead to obvious phenotypic differences such as antibiotic resistance or the well-known mucoid phenotype of Pseudomonas aeruginosa. Other mutations may be just as important but harder to detect such as increased mutation rates, cell surface changes, and shifts in metabolism and nutrient acquisition. Remarkably, many of the adaptations occur again and again in different patients, signaling that bacteria are adapting to solve specific challenges in the CF respiratory tract. This parallel evolution even extends across distinct bacterial species. This review addresses the bacterial systems that are known to change in long-term CF infections with a special emphasis on cross-species comparisons. Consideration is given to how adaptation may impact health in CF, and the possible evolutionary mechanisms that lead to the repeated parallel adaptations.

Bacterial evolution during human infection: Adapt and live or adapt and die

Microbes are constantly evolving. Laboratory studies of bacterial evolution increase our understanding of evolutionary dynamics, identify adaptive changes, and answer important questions that impact human health. During bacterial infections in humans, however, the evolutionary parameters acting on infecting populations are likely to be much more complex than those that can be tested in the laboratory. Nonetheless, human infections can be thought of as naturally occurring in vivo bacterial evolution experiments, which can teach us about antibiotic resistance, pathogenesis, and transmission. Here, we review recent advances in the study of within-host bacterial evolution during human infection and discuss practical considerations for conducting such studies. We focus on 2 possible outcomes for de novo adaptive mutations, which we have termed "adapt-and-live" and "adapt-and-die." In the adapt-and-live scenario, a mutation is long lived, enabling its transmission on to other individuals, or the establishment of chronic infection. In the adapt-and-die scenario, a mutation is rapidly extinguished, either because it carries a substantial fitness cost, it arises within tissues that block transmission to new hosts, it is outcompeted by more fit clones, or the infection resolves. Adapt-and-die mutations can provide rich information about selection pressures in vivo, yet they can easily elude detection because they are short lived, may be more difficult to sample, or could be maladaptive in the long term. Understanding how bacteria adapt under each of these scenarios can reveal new insights about the basic biology of pathogenic microbes and could aid in the design of new translational approaches to combat bacterial infections.

In vivo evolution of an emerging zoonotic bacterial pathogen in an immunocompromised human host

Zoonotic transfer of animal pathogens to human hosts can generate novel agents, but the genetic events following such host jumps are not well studied. Here we characterize the mechanisms driving adaptive evolution of the emerging zoonotic pathogen Bordetella hinzii in a patient with interleukin-12 receptor β1 deficiency. Genomic sequencing of 24 B. hinzii isolates cultured from blood and stool over 45 months revealed a clonal lineage that had undergone extensive within-host genetic and phenotypic diversification. Twenty of 24 isolates shared an E9G substitution in the DNA polymerase III ε-subunit active site, resulting in a proofreading deficiency. Within this proofreading-deficient clade, multiple lineages with mutations in DNA repair genes and altered mutational spectra emerged and dominated clinical cultures for more than 12 months. Multiple enzymes of the tricarboxylic acid cycle and gluconeogenesis pathways were repeatedly mutated, suggesting rapid metabolic adaptation to the human environment. Furthermore, an excess of G:C > T:A transversions suggested that oxidative stress shaped genetic diversification during adaptation. We propose that inactivation of DNA proofreading activity in combination with prolonged, but sub-lethal, oxidative attack resulting from the underlying host immunodeficiency facilitated rapid genomic adaptation. These findings suggest a fundamental role for host immune phenotype in shaping pathogen evolution following zoonotic infection.

Bordetella hinzii is an emerging pathogen with zoonotic risk to humans, known to be able to cause respiratory tract infection, bacteremia and endocarditis. Here, applying whole genome sequencing to bacterial isolates, the authors characterize the mechanisms driving adaptive evolution in B. hinzii in a patient with interleukin-12 receptor β1 deficiency, suggesting a role for host immune phenotype in shaping within-host pathogen evolution following zoonotic infection.

Evaluation of Reproductive Tract Cytokines in Post-partum Beef Cows Relating to Reproductive Microbiota and Fertility Outcomes

The activity of the immune system in the reproductive tract has been proven to be crucial in the response to uterine diseases, normal reproductive functions, and tolerance to the allogeneic fetus during pregnancy. The objectives of the current study were to (1) evaluate uterine and vaginal cytokine concentrations in postpartum cows undergoing estrus synchronization followed by timed artificial insemination (TAI) and (2) correlate bacterial communities with cytokine concentrations. Postpartum Angus cows (n = 20) were subjected to a 7-Day Co-Synch protocol with pre-synchronization beginning 21 days prior (d −21) to TAI (d 0). Uterine and vaginal flushes were collected on d −21 and −2. Pregnancy was determined by transrectal ultrasound on d 30. Cytokines include interleukin (IL)-1b, IL-6, IL-10, transforming growth factor beta (TGF-β), and immunoglobin A (IgA) and concentrations were determined by commercial ELISA kits. No differences by day or pregnancy status in cytokine concentrations were detected in vaginal samples. No differences by day or pregnancy status in IgA, IL-10, or IL-1b concentrations were detected in uterine samples. Overall TGF-β concentrations in the uterus were greater in resulting pregnant than non-pregnant cows (44.0 ± 13.4 pg/mL vs. 14.7 ± 4.9 pg/mL; P = 0.047). Uterine TGF-β was correlated with the relative abundance of genera Treponema (r = −0.668; P = 0.049) in resulting non-pregnant cows on d −21 and with the relative abundance of genera Ureaplasma (r = 0.901; P = 0.0004) in resulting pregnant cows on d −2. In resulting pregnant animals, a tendency for a strong correlation was detected between d −2 progesterone concentrations and uterine TGF-β concentrations (r = 0.591, P = 0.07). Overall IL-6 concentrations in the uterus were greater in resulting non-pregnant than pregnant cows (198.7 ± 21.8 pg/mL vs. 144.3 ± 16.1 pg/mL; P = 0.045). A correlation was also detected between uterine IL-6 concentrations and the relative abundance of genera Butyrivibrio (r = 0.742; P = 0.022) in resulting non-pregnant cows on d −21. These results suggest possible relationships between different bacterial communities and cytokine concentrations within the uterus of beef cattle prior to TAI that may ultimately affect fertility outcomes.

Pan-Genome Analysis Reveals Host-Specific Functional Divergences in Burkholderia gladioli

Burkholderia gladioli has high versatility and adaptability to various ecological niches. Here, we constructed a pan-genome using 14 genome sequences of B. gladioli, which originate from different niches, including gladiolus, rice, humans, and nature. Functional roles of core and niche-associated genomes were investigated by pathway enrichment analyses. Consequently, we inferred the uniquely important role of niche-associated genomes in (1) selenium availability during competition with gladiolus host; (2) aromatic compound degradation in seed-borne and crude oil-accumulated environments, and (3) stress-induced DNA repair system/recombination in the cystic fibrosis-niche. We also identified the conservation of the rhizomide biosynthetic gene cluster in all the B. gladioli strains and the concentrated distribution of this cluster in human isolates. It was confirmed the absence of complete CRISPR/Cas system in both plant and human pathogenic B. gladioli and the presence of the system in B. gladioli living in nature, possibly reflecting the inverse relationship between CRISPR/Cas system and virulence.

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.

Bacterial defenses against a natural antibiotic promote collateral resilience to clinical antibiotics

Bacterial opportunistic human pathogens frequently exhibit intrinsic antibiotic tolerance and resistance, resulting in infections that can be nearly impossible to eradicate. We asked whether this recalcitrance could be driven by these organisms' evolutionary history as environmental microbes that engage in chemical warfare. Using Pseudomonas aeruginosa as a model, we demonstrate that the self-produced antibiotic pyocyanin (PYO) activates defenses that confer collateral tolerance specifically to structurally similar synthetic clinical antibiotics. Non-PYO-producing opportunistic pathogens, such as members of the Burkholderia cepacia complex, likewise display elevated antibiotic tolerance when cocultured with PYO-producing strains. Furthermore, by widening the population bottleneck that occurs during antibiotic selection and promoting the establishment of a more diverse range of mutant lineages, PYO increases apparent rates of mutation to antibiotic resistance to a degree that can rival clinically relevant hypermutator strains. Together, these results reveal an overlooked mechanism by which opportunistic pathogens that produce natural toxins can dramatically modulate the efficacy of clinical antibiotics and the evolution of antibiotic resistance, both for themselves and other members of clinically relevant polymicrobial communities.

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.

One gene, multiple ecological strategies: A biofilm regulator is a capacitor for sustainable diversity

Many organisms, including bacteria, live in fluctuating environments that require attachment and dispersal. These lifestyle decisions require processing of multiple external signals by several genetic pathways, but how they are integrated is largely unknown. We conducted multiple evolution experiments totaling >20,000 generations with Burkholderia cenocepacia populations grown in a model of the biofilm life cycle and identified parallel mutations in one gene, rpfR, that is a conserved central regulator. Because RpfR has multiple sensor and catalytic domains, different mutations can produce different ecological strategies that can coexist and even increase net growth. This study demonstrates that a single gene may coordinate complex life histories in biofilm-dwelling bacteria and that selection in defined environments can reshape niche breadth by single mutations.

Many bacteria cycle between sessile and motile forms in which they must sense and respond to internal and external signals to coordinate appropriate physiology. Maintaining fitness requires genetic networks that have been honed in variable environments to integrate these signals. The identity of the major regulators and how their control mechanisms evolved remain largely unknown in most organisms. During four different evolution experiments with the opportunist betaproteobacterium Burkholderia cenocepacia in a biofilm model, mutations were most frequently selected in the conserved gene rpfR. RpfR uniquely integrates two major signaling systems—quorum sensing and the motile–sessile switch mediated by cyclic-di-GMP—by two domains that sense, respond to, and control the synthesis of the autoinducer cis-2-dodecenoic acid (BDSF). The BDSF response in turn regulates the activity of diguanylate cyclase and phosphodiesterase domains acting on cyclic-di-GMP. Parallel adaptive substitutions evolved in each of these domains to produce unique life history strategies by regulating cyclic-di-GMP levels, global transcriptional responses, biofilm production, and polysaccharide composition. These phenotypes translated into distinct ecology and biofilm structures that enabled mutants to coexist and produce more biomass than expected from their constituents grown alone. This study shows that when bacterial populations are selected in environments challenging the limits of their plasticity, the evolved mutations not only alter genes at the nexus of signaling networks but also reveal the scope of their regulatory functions.

Impact of insertion sequences on convergent evolution of Shigella species

Shigella species are specialised lineages of Escherichia coli that have converged to become human-adapted and cause dysentery by invading human gut epithelial cells. Most studies of Shigella evolution have been restricted to comparisons of single representatives of each species; and population genomic studies of individual Shigella species have focused on genomic variation caused by single nucleotide variants and ignored the contribution of insertion sequences (IS) which are highly prevalent in Shigella genomes. Here, we investigate the distribution and evolutionary dynamics of IS within populations of Shigella dysenteriae Sd1, Shigella sonnei and Shigella flexneri. We find that five IS (IS1, IS2, IS4, IS600 and IS911) have undergone expansion in all Shigella species, creating substantial strain-to-strain variation within each population and contributing to convergent patterns of functional gene loss within and between species. We find that IS expansion and genome degradation are most advanced in S. dysenteriae and least advanced in S. sonnei; and using genome-scale models of metabolism we show that Shigella species display convergent loss of core E. coli metabolic capabilities, with S. sonnei and S. flexneri following a similar trajectory of metabolic streamlining to that of S. dysenteriae. This study highlights the importance of IS to the evolution of Shigella and provides a framework for the investigation of IS dynamics and metabolic reduction in other bacterial species.

Phenotypic traits of Burkholderia spp. associated with ecological adaptation and plant-host interaction

Burkholderia species have different lifestyles establishing mutualist or pathogenic associations with plants and animals. Changes in the ecological behavior of these bacteria may depend on genetic variations in response to niche adaptation. Here, we studied 15 Burkholderia strains isolated from different environments with respect to genetic and phenotypic traits. By Multilocus Sequence Analysis (MLSA) these isolates fell into 6 distinct groups. MLSA clusters did not correlate with strain antibiotic sensitivity, but with the bacterial ability to produce antimicrobial compounds and control orchid necrosis. Further, the B. seminalis strain TC3.4.2R3, a mutualistic bacterium, was inoculated into orchid plants and the interaction with the host was evaluated by analyzing the plant response and the bacterial oxidative stress response in planta. TC3.4.2R3 responded to plant colonization by increasing its own growth rate and by differential gene regulation upon oxidative stress caused by the plant, while reducing the plant's membrane lipid peroxidation. The bacterial responses to oxidative stress were recapitulated by bacterial exposure to the herbicide paraquat. We suggest that the ability of Burkholderia species to successfully establish in the rhizosphere correlates with genetic variation, whereas traits associated with antibiotic resistance are more likely to be categorized as strain specific.

Variation of Burkholderia cenocepacia cell wall morphology and mechanical properties during cystic fibrosis lung infection, assessed by atomic force microscopy

The influence that Burkholderia cenocepacia adaptive evolution during long-term infection in cystic fibrosis (CF) patients has on cell wall morphology and mechanical properties is poorly understood despite their crucial role in cell physiology, persistent infection and pathogenesis. Cell wall morphology and physical properties of three B. cenocepacia isolates collected from a CF patient over a period of 3.5 years were compared using atomic force microscopy (AFM). These serial clonal variants include the first isolate retrieved from the patient and two late isolates obtained after three years of infection and before the patient’s death with cepacia syndrome. A consistent and progressive decrease of cell height and a cell shape evolution during infection, from the typical rods to morphology closer to cocci, were observed. The images of cells grown in biofilms showed an identical cell size reduction pattern. Additionally, the apparent elasticity modulus significantly decreases from the early isolate to the last clonal variant retrieved from the patient but the intermediary highly antibiotic resistant clonal isolate showed the highest elasticity values. Concerning the adhesion of bacteria surface to the AFM tip, the first isolate was found to adhere better than the late isolates whose lipopolysaccharide (LPS) structure loss the O-antigen (OAg) during CF infection. The OAg is known to influence Gram-negative bacteria adhesion and be an important factor in B. cenocepacia adaptation to chronic infection. Results reinforce the concept of the occurrence of phenotypic heterogeneity and adaptive evolution, also at the level of cell size, form, envelope topography and physical properties during long-term infection.

The Effect of 2-Thiocyanatopyridine Derivative 11026103 on Burkholderia Cenocepacia: Resistance Mechanisms and Systemic Impact

Bacteria of the Burkholderia cepacia complex (Bcc) are associated with significant decline of lung functions in cystic fibrosis patients. Bcc infections are virtually impossible to eradicate due to their irresponsiveness to antibiotics. The 2-thiocyanatopyridine derivative 11026103 is a novel, synthetic compound active against Burkholderia cenocepacia. To characterize mechanisms of resistance to 11026103, B. cenocepacia was subjected to chemical mutagenesis, followed by whole genome sequencing. Parallel mutations in resistant isolates were localized in a regulatory protein of the efflux system Resistance-Nodulation-Division (RND)-9 (BCAM1948), RNA polymerase sigma factor (BCAL2462) and its cognate putative anti-sigma factor (BCAL2461). Transcriptomic analysis identified positive regulation of a major facilitator superfamily (MFS) efflux system BCAL1510-1512 by BCAL2462. Artificial overexpression of both efflux systems increased resistance to the compound. The effect of 11026103 on B. cenocepacia was analyzed by RNA-Seq and a competitive fitness assay utilizing an essential gene knockdown mutant library. 11026103 exerted a pleiotropic effect on transcription including profound downregulation of cluster of orthologous groups (COG) category “Translation, ribosomal structure, and biogenesis”. The competitive fitness assay identified many genes which modulated susceptibility to 11026103. In summary, 11026103 exerts a pleiotropic cellular response in B. cenocepacia which can be prevented by efflux system-mediated export.

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 and metagenomics of the metallomes

Recent achievements and advances in comparative genomic and metagenomic analyses of trace metals were reviewed and discussed.

A genome-wide association analysis reveals a potential role for recombination in the evolution of antimicrobial resistance in Burkholderia multivorans

Cystic fibrosis (CF) lung infections caused by members of the Burkholderia cepacia complex, such as Burkholderia multivorans, are associated with high rates of mortality and morbidity. We performed a population genomics study of 111 B. multivorans sputum isolates from one CF patient through three stages of infection including an early incident isolate, deep sampling of a one-year period of chronic infection occurring weeks before a lung transplant, and deep sampling of a post-transplant infection. We reconstructed the evolutionary history of the population and used a lineage-controlled genome-wide association study (GWAS) approach to identify genetic variants associated with antibiotic resistance. We found the incident isolate was basally related to the rest of the strains and more susceptible to antibiotics from three classes (β-lactams, aminoglycosides, quinolones). The chronic infection isolates diversified into multiple, distinct genetic lineages and showed reduced antimicrobial susceptibility to the same antibiotics. The post-transplant reinfection isolates derived from the same source as the incident isolate and were genetically distinct from the chronic isolates. They also had a level of susceptibility in between that of the incident and chronic isolates. We identified numerous examples of potential parallel pathoadaptation, in which multiple mutations were found in the same locus or even codon. The set of parallel pathoadaptive loci was enriched for functions associated with virulence and resistance. Our GWAS analysis identified statistical associations between a polymorphism in the ampD locus with resistance to β-lactams, and polymorphisms in an araC transcriptional regulator and an outer membrane porin with resistance to both aminoglycosides and quinolones. Additionally, these three loci were independently mutated four, three and two times, respectively, providing further support for parallel pathoadaptation. Finally, we identified a minimum of 14 recombination events, and observed that loci carrying putative parallel pathoadaptations and polymorphisms statistically associated with β-lactam resistance were over-represented in these recombinogenic regions.

Cystic fibrosis (CF) is the most common lethal genetic disorder affecting individuals of European descent. Most CF patients die at a young age due to chronic lung infections. Among the organisms involved in these infections are bacteria from the Burkholderia cepacia complex (BCC), which are strongly associated with poor clinical prognosis. This study examines how the most prevalent BCC species among CF patients, B. multivorans, evolves within a single CF patient by studying the first B. multivorans isolate recovered from the patient, one hundred isolates recovered over a one year period during the chronic infection phase, and an additional ten isolates recovered after the reinfection of the transplanted lungs. We found that B. multivorans diversify phenotypically and genetically within the CF lung over the course of the infection, and evolves into a complex population during the chronic infection phase. We found that isolates collected from the post-transplant reinfection were more closely related to descendants of the original isolate rather than those recovered in the chronic infection. We identify genetic variants statistically associated with resistance to the antibiotics, and showed that some of these variants were found in regions that show patterns of recombination (genetic exchange) between strains. We also found that genes which were mutated multiple times during overall infection were more likely to be found in regions showing signals consistent with recombination. The presence of multiple independent mutations in a gene is a very strong signal that the gene helps bacteria adapt to their environment. Overall, this study provides insight into how pathogens adapt to the host during long-term infections, specific genes associated with antibiotic resistance, and the origin of new and recurrent infections.