Phenotypic diversity and genotypic flexibility of Burkholderia cenocepacia during long-term chronic infection of cystic fibrosis lungs

Within-Host Evolution of Bacterial Pathogens in Acute and Chronic Infection

Bacterial pathogens undergo remarkable adaptive change in response to the selective forces they encounter during host colonization and infection. Studies performed over the past few decades have demonstrated that many general evolutionary processes can be discerned during the course of host adaptation, including genetic diversification of lineages, clonal succession events, convergent evolution, and balanced fitness trade-offs. In some cases, elevated mutation rates resulting from mismatch repair or proofreading deficiencies accelerate evolution, and active mobile genetic elements or phages may facilitate genome plasticity. The host immune response provides another critical component of the fitness landscapes guiding adaptation, and selection operating on pathogens at this level may lead to immune evasion and the establishment of chronic infection. This review summarizes recent advances in this field, with a special focus on different forms of bacterial genome plasticity in the context of infection, and considers clinical consequences of adaptive changes for the host.

Influence of genomic variations on glanders serodiagnostic antigens using integrative genomic and transcriptomic approaches

Glanders is a highly contagious and life-threatening zoonotic disease caused by Burkholderia mallei (B. mallei). Without an effective vaccine or treatment, early diagnosis has been regarded as the most effective method to prevent glanders transmission. Currently, the diagnosis of glanders is heavily reliant on serological tests. However, given that markedly different host immune responses can be elicited by genetically different strains of the same bacterial species, infection by B. mallei, whose genome is unstable and plastic, may result in various immune responses. This variability can make the serodiagnosis of glanders challenging. Therefore, there is a need for a comprehensive understanding and assessment of how B. mallei genomic variations impact the appropriateness of specific target antigens for glanders serodiagnosis. In this study, we investigated how genomic variations in the B. mallei genome affect gene content (gene presence/absence) and expression, with a special focus on antigens used or potentially used in serodiagnosis. In all the genome sequences of B. mallei isolates available in NCBI's RefSeq database (accessed in July 2023) and in-house sequenced samples, extensive small and large variations were observed when compared to the type strain ATCC 23344. Further pan-genome analysis of those assemblies revealed variations of gene content among all available genomes of B. mallei. Specifically, differences in gene content ranging from 31 to 715 genes with an average of 334 gene presence-absence variations were found in strains with complete or chromosome-level genome assemblies, using the ATCC 23344 strain as a reference. The affected genes included some encoded proteins used as serodiagnostic antigens, which were lost due mainly to structural variations. Additionally, a transcriptomic analysis was performed using the type strain ATCC 23344 and strain Zagreb which has been widely utilized to produce glanders antigens. In total, 388 significant differentially expressed genes were identified between these two strains, including genes related to bacterial pathogenesis and virulence, some of which were associated with genomic variations, particularly structural variations. To our knowledge, this is the first comprehensive study to uncover the impacts of genetic variations of B. mallei on its gene content and expression. These differences would have significant impacts on host innate and adaptive immunity, including antibody production, during infection. This study provides novel insights into B. mallei genetic variants, knowledge which will help to improve glanders serodiagnosis.

An IS-mediated, RecA-dependent, bet-hedging strategy in Burkholderia thailandensis

Adaptation to fluctuating environmental conditions is difficult to achieve. Phase variation mechanisms can overcome this difficulty by altering genomic architecture in a subset of individuals, creating a phenotypically heterogeneous population with subpopulations optimized to persist when conditions change, or are encountered, suddenly. We have identified a phase variation system in Burkholderia thailandensis that generates a genotypically and phenotypically heterogeneous population. Genetic analyses revealed that RecA-mediated homologous recombination between a pair of insertion sequence (IS) 2-like elements duplicates a 208.6 kb region of DNA that contains 157 coding sequences. RecA-mediated homologous recombination also resolves merodiploids, and hence copy number of the region is varied and dynamic within populations. We showed that the presence of two or more copies of the region is advantageous for growth in a biofilm, and a single copy is advantageous during planktonic growth. While IS elements are well known to contribute to evolution through gene inactivation, polar effects on downstream genes, and altering genomic architecture, we believe that this system represents a rare example of IS element-mediated evolution in which the IS elements provide homologous sequences for amplification of a chromosomal region that provides a selective advantage under specific growth conditions, thereby expanding the lifestyle repertoire of the 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.

Within-host evolution of bacterial pathogens during persistent infection of humans

Many bacterial pathogens can form persistent infections, providing an infectious reservoir, which allows for infection of new hosts. Currently, the molecular mechanisms and evolutionary dynamics driving persistence are still not well-understood. High-throughput sequencing methods have enabled the study of within-host evolution of persistent bacterial pathogens, revealing common trends among bacterial species in how they adapt to persist. We will focus on trends emerging from longitudinal human-cohort studies, including i) genome-size reduction, ii) metabolic adaptation to the host, iii) antimicrobial resistance, iv) changes in virulence and the bacterial cell surface, and v) hypermutation, and comment on where the field should focus going forward.

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.

Metabolomic profiling of Burkholderia cenocepacia in synthetic cystic fibrosis sputum medium reveals nutrient environment-specific production of virulence factors

Infections by Burkholderia cenocepacia lead to life-threatening disease in immunocompromised individuals, including those living with cystic fibrosis (CF). While genetic variation in various B. cenocepacia strains has been reported, it remains unclear how the chemical environment of CF lung influences the production of small molecule virulence factors by these strains. Here we compare metabolomes of three clinical B. cenocepacia strains in synthetic CF sputum medium (SCFM2) and in a routine laboratory medium (LB), in the presence and absence of the antibiotic trimethoprim. Using a mass spectrometry-based untargeted metabolomics approach, we identify several compound classes which are differentially produced in SCFM2 compared to LB media, including siderophores, antimicrobials, quorum sensing signals, and various lipids. Furthermore, we describe that specific metabolites are induced in the presence of the antibiotic trimethoprim only in SCFM2 when compared to LB. Herein, C13-acyl-homoserine lactone, a quorum sensing signal previously not known to be produced by B. cenocepacia as well as pyochelin-type siderophores were exclusively detected during growth in SCFM2 in the presence of trimethoprim. The comparative metabolomics approach described in this study provides insight into environment-dependent production of secondary metabolites by B. cenocepacia strains and suggests future work which could identify personalized strain-specific regulatory mechanisms involved in production of secondary metabolites. Investigations into whether antibiotics with different mechanisms of action induce similar metabolic alterations will inform development of combination treatments aimed at effective clearance of Burkholderia spp. pathogens.

The Burkholderia cenocepacia Type VI Secretion System Effector TecA Is a Virulence Factor in Mouse Models of Lung Infection

Burkholderia cenocepacia is a member of the Burkholderia cepacia complex (Bcc), a group of bacteria with members responsible for causing lung infections in cystic fibrosis (CF) patients. The most severe outcome of Bcc infection in CF patients is cepacia syndrome, a disease characterized by necrotizing pneumonia with bacteremia and sepsis. B. cenocepacia is strongly associated with cepacia syndrome, making it one of the most virulent members of the Bcc. Mechanisms underlying the pathogenesis of B. cenocepacia in lung infections and cepacia syndrome remain to be uncovered. B. cenocepacia is primarily an intracellular pathogen and encodes the type VI secretion system (T6SS) effector TecA, which is translocated into host phagocytes. TecA is a deamidase that inactivates multiple Rho GTPases, including RhoA. Inactivation of RhoA by TecA triggers assembly of the pyrin inflammasome, leading to secretion of proinflammatory cytokines, such as interleukin-1β, from macrophages. Previous work with the B. cenocepacia clinical isolate J2315 showed that TecA increases immunopathology during acute lung infection in C57BL/6 mice and suggested that this effector acts as a virulence factor by triggering assembly of the pyrin inflammasome. Here, we extend these results using a second B. cenocepacia clinical isolate, AU1054, to demonstrate that TecA exacerbates weight loss and lethality during lung infection in C57BL/6 mice and mice engineered to have a CF genotype. Unexpectedly, pyrin was dispensable for TecA virulence activity in both mouse infection models. Our findings establish that TecA is a B. cenocepacia virulence factor that exacerbates lung inflammation, weight loss, and lethality in mouse infection models.

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.

Evolution towards Virulence in a Burkholderia Two-Component System

Bacteria in the Burkholderia cepacia complex (BCC) are significant pathogens for people with cystic fibrosis (CF) and are often extensively antibiotic resistant. Here, we assess the impacts of clinically observed mutations in fixL, which encodes the sensor histidine kinase FixL. FixL along with FixJ compose a two-component system that regulates multiple phenotypes. Mutations in fixL across two species, B. dolosa and B. multivorans, have shown evidence of positive selection during chronic lung infection in CF. Herein, we find that BCC carrying the conserved, ancestral fixL sequence have lower survival in macrophages and in murine pneumonia models than mutants carrying evolved fixL sequences associated with clinical decline in CF patients. In vitro phosphotransfer experiments found that one evolved FixL protein, W439S, has a reduced ability to autophosphorylate and phosphorylate FixJ, while LacZ reporter experiments demonstrate that B. dolosa carrying evolved fixL alleles has reduced fix pathway activity. Interestingly, B. dolosa carrying evolved fixL alleles was less fit in a soil assay than those strains carrying the ancestral allele, demonstrating that increased survival of these variants in macrophages and the murine lung comes at a potential expense in their environmental reservoir. Thus, modulation of the two-component system encoded by fixLJ by point mutations is one mechanism that allows BCC to adapt to the host infection environment.

Beyond the pan-genome: current perspectives on the functional and practical outcomes of the distributed genome hypothesis

The principle of monoclonality with regard to bacterial infections was considered immutable prior to 30 years ago. This view, espoused by Koch for acute infections, has proven inadequate regarding chronic infections as persistence requires multiple forms of heterogeneity among the bacterial population. This understanding of bacterial plurality emerged from a synthesis of what-were-then novel technologies in molecular biology and imaging science. These technologies demonstrated that bacteria have complex life cycles, polymicrobial ecologies, and evolve in situ via the horizontal exchange of genic characters. Thus, there is an ongoing generation of diversity during infection that results in far more highly complex microbial communities than previously envisioned. This perspective is based on the fundamental tenet that the bacteria within an infecting population display genotypic diversity, including gene possession differences, which result from horizontal gene transfer mechanisms including transformation, conjugation, and transduction. This understanding is embodied in the concepts of the supragenome/pan-genome and the distributed genome hypothesis (DGH). These paradigms have fostered multiple researches in diverse areas of bacterial ecology including host-bacterial interactions covering the gamut of symbiotic relationships including mutualism, commensalism, and parasitism. With regard to the human host, within each of these symbiotic relationships all bacterial species possess attributes that contribute to colonization and persistence; those species/strains that are pathogenic also encode traits for invasion and metastases. Herein we provide an update on our understanding of bacterial plurality and discuss potential applications in diagnostics, therapeutics, and vaccinology based on perspectives provided by the DGH with regard to the evolution of pathogenicity.

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.

Epidemiology of E. coli in Cystic Fibrosis Airways Demonstrates the Capacity for Persistent Infection but Not Patient-Patient Transmission

Escherichia coli is frequently isolated from the respiratory secretions of cystic fibrosis (CF) patients yet is not considered a classical CF pathogen. Accordingly, little is known about the natural history of this organism in the CF airways, as well as the potential for patient-to-patient transmission. Patients attending the Calgary Adult CF Clinic (CACFC) between January 1983 and December 2016 with at least one E. coli-positive sputum culture were identified by retrospective review. Annual E. coli isolates from the CACFC biobank from each patient were typed by pulsed-field gel electrophoresis (PFGE) and isolates belonging to shared pulsotypes were sequenced. Single nucleotide polymorphism (SNP) and phylogenetic analysis were used to investigate the natural history of E. coli infection and identify potential transmission events. Forty-five patients with E. coli-positive sputum cultures were identified. Most patients had a single infection episode with a single pulsotype, while replacement of an initial pulsotype with a second was observed in three patients. Twenty-four had E. coli recovered from their sputum more than once and 18 patients had persistent infections (E. coli carriage >6 months with ≥3 positive cultures). Shared pulsotypes corresponded to known extraintestinal pathogenic E. coli strains: ST-131, ST-73, and ST-1193. Phylogenetic relationships and SNP distances among isolates within shared pulsotypes were consistent with independent acquisition of E. coli by individual patients. Most recent common ancestor date estimates of isolates between patients were inconsistent with patient-to-patient transmission. E. coli infection in CF is a dynamic process that appears to be characterized by independent acquisition within our patient population and carriage of unique sets of strains over time by individual patients.

Panel 3: Genomics, precision medicine and targeted therapies

OBJECTIVE

To review the most recent advances in human and bacterial genomics as applied to pathogenesis and clinical management of otitis media.

DATA SOURCES

PubMed articles published since the last meeting in June 2015 up to June 2019.

REVIEW METHODS

A panel of experts in human and bacterial genomics of otitis media was formed. Each panel member reviewed the literature in their respective fields and wrote draft reviews. The reviews were shared with all panel members, and a merged draft was created. The panel met at the 20th International Symposium on Recent Advances in Otitis Media in June 2019, discussed the review and refined the content. A final draft was made, circulated, and approved by the panel members.

CONCLUSION

Trans-disciplinary approaches applying pan-omic technologies to identify human susceptibility to otitis media and to understand microbial population dynamics, patho-adaptation and virulence mechanisms are crucial to the development of novel, personalized therapeutics and prevention strategies for otitis media.

IMPLICATIONS FOR PRACTICE

In the future otitis media prevention strategies may be augmented by mucosal immunization, combination vaccines targeting multiple pathogens, and modulation of the middle ear microbiome. Both treatment and vaccination may be tailored to an individual's otitis media phenotype as defined by molecular profiles obtained by using rapidly developing techniques in microbial and host genomics.

Genomic and phenotypic characterization of Burkholderia isolates from the potable water system of the International Space Station

The opportunistic pathogens Burkholderia cepacia and Burkholderia contaminans, both genomovars of the Burkholderia cepacia complex (BCC), are frequently cultured from the potable water dispenser (PWD) of the International Space Station (ISS). Here, we sequenced the genomes and conducted phenotypic assays to characterize these Burkholderia isolates. All recovered isolates of the two species fall within monophyletic clades based on phylogenomic trees of conserved single-copy core genes. Within species, the ISS-derived isolates all demonstrate greater than 99% average nucleotide identity (with 95-99% of genomes aligning) and share around 90% of the identified gene clusters from a pangenomic analysis-suggesting that the two groups are each composed of highly similar genomic lineages and their members may have all stemmed from the same two founding populations. The differences that can be observed between the recovered isolates at the pangenomic level are primarily located within putative plasmids. Phenotypically, macrophage intracellularization and lysis occurred at generally similar rates between all ISS-derived isolates, as well as with their respective type-terrestrial strain references. All ISS-derived isolates exhibited antibiotic sensitivity similar to that of the terrestrial reference strains, and minimal differences between isolates were observed. With a few exceptions, biofilm formation rates were generally consistent across each species. And lastly, though isolation date does not necessarily provide any insight into how long a given isolate had been aboard the ISS, none of the assayed physiology correlated with either date of isolation or distances based on nucleotide variation. Overall, we find that while the populations of Burkholderia present in the ISS PWS each maintain virulence, they are likely are not more virulent than those that might be encountered on planet and remain susceptible to clinically used antibiotics.

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.

Sex-Specific Linkages Between Taxonomic and Functional Profiles of Tick Gut Microbiomes

Ticks transmit the most diverse array of disease agents and harbor one of the most diverse microbial communities. Major progress has been made in the characterization of the taxonomic profiles of tick microbiota. However, the functional profiles of tick microbiome have been comparatively less studied. In this proof of concept we used state-of-the-art functional metagenomics analytical tools to explore previously reported datasets of bacteria found in male and female Ixodes ovatus, Ixodes persulcatus, and Amblyomma variegatum. Results showed that both taxonomic and functional profiles have differences between sexes of the same species. KEGG pathway analysis revealed that male and female of the same species had major differences in the abundance of genes involved in different metabolic pathways including vitamin B, amino acids, carbohydrates, nucleotides, and antibiotics among others. Partial reconstruction of metabolic pathways using KEGG enzymes suggests that tick microbiome form a complex metabolic network that may increase microbial community resilience and adaptability. Linkage analysis between taxonomic and functional profiles showed that among the KEGG enzymes with differential abundance in male and female ticks only 12% were present in single bacterial genera. The rest of these enzymes were found in more than two bacterial genera, and 27% of them were found in five up to ten bacterial genera. Comparison of bacterial genera contributing to the differences in the taxonomic and functional profiles of males and females revealed that while a small group of bacteria has a dual-role, most of the bacteria contribute only to functional or taxonomic differentiation between sexes. Results suggest that the different life styles of male and female ticks exert sex-specific evolutionary pressures that act independently on the phenomes (set of phenotypes) and genomes of bacteria in tick gut microbiota. We conclude that functional redundancy is a fundamental property of male and female tick microbiota and propose that functional metagenomics should be combined with taxonomic profiling of microbiota because both analyses are complementary.

Mucoid switch in Burkholderia cepacia complex bacteria: Triggers, molecular mechanisms and implications in pathogenesis

Bacteria produce a vast range of exopolysaccharides (EPSs) to thrive in diverse environmental niches and often display a mucoid phenotype in solid media. One such exopolysaccharide, cepacian, is produced by bacteria of the genus Burkholderia and is of interest due to its role in pathogenesis associated with lung infections in cystic fibrosis (CF) patients. Cepacian is a repeat-unit polymer that has been implicated in biofilm formation, immune system evasion, interaction with host cells, resistance against antimicrobials, and virulence. Its biosynthesis proceeds through the Wzy-dependent polymerization and secretion mechanism, which requires a multienzymatic complex. Key aspects of its structure, genetic organization, and the regulatory network involved in mucoid switch and regulation of cepacian biosynthesis at transcriptional and posttranscriptional levels are reviewed. It is also evaluated the importance of cepacian biosynthesis/regulation key players as evolutionary targets of selection and highlighted the complexity of the regulatory network, which allows cells to coordinate the expression of metabolic functions to the ones of the cell wall, in order to be successful in ever changing environments, including in the interaction with host cells.

Identification of Diverse Integron and Plasmid Structures Carrying a Novel Carbapenemase Among Pseudomonas Species

A novel carbapenem-hydrolyzing beta-lactamase, called IMP-63, was identified in three clonally distinct strains of Pseudomonas aeruginosa and two strains of Pseudomonas putida isolated within a 4 year timeframe in three French hospitals. The bla_IMP–63 gene that encodes this carbapenemase turned out to be located in the variable region of four integrons (In1297, In1574, In1573, and In1572) and to coexist with novel or rare gene cassettes (fosM, gcu170, gcuF1) and insertion elements (ISPsp7v, ISPa16v). All these integrons except one (In1574) were flanked by a copy of insertion sequence ISPa17 next to the orf6 putative gene, and were carried by non-conjugative plasmids (pNECK1, pROUSS1, pROUSS2, pROUE1). These plasmids exhibit unique modular structures and partial sequence homologies with plasmids previously identified in various non-fermenting environmental Gram-negative species. Lines of evidence suggest that ISPa17 promoted en bloc the transposition of IMP-63-encoding integrons on these different plasmids. As demonstrated by genotyping experiments, isolates of P. aeruginosa harboring the 28.9-kb plasmid pNECK1 and belonging to international “high-risk” clone ST308 were responsible for an outbreak in one hospital. Collectively, these data provide an insight into the complex and unpredictable routes of diffusion of some resistance determinants, here bla_IMP–63, among Pseudomonas species.

Comparative Genome Analysis of an Extensively Drug-Resistant Isolate of Avian Sequence Type 167 Escherichia coli Strain Sanji with Novel In Silico Serotype O89b:H9

E. coli strain Sanji is the first sequenced and analyzed genome of the recently emerged pathogenic XDR strains with sequence type ST167 and novel in silico serotype O89b:H9. Comparison of the genomes of Sanji with other ST167 strains revealed distinct sets of different plasmids, mobile IS elements, and antibiotic resistance genes in each genome, indicating that there exist multiple paths toward achieving XDR. The emergence of these pathogenic ST167 E. coli strains with diverse XDR capabilities highlights the difficulty of preventing or mitigating the development of XDR properties in bacteria and points to the importance of better understanding of the shared underlying virulence mechanisms and physiology of pathogenic bacteria.

Extensive drug resistance (XDR) is an escalating global problem. Escherichia coli strain Sanji was isolated from an outbreak of pheasant colibacillosis in Fujian province, China, in 2011. This strain has XDR properties, exhibiting sensitivity to carbapenems but no other classes of known antibiotics. Whole-genome sequencing revealed a total of 32 known antibiotic resistance genes, many associated with insertion sequence 26 (IS26) elements. These were found on the Sanji chromosome and 2 of its 6 plasmids, pSJ_255 and pSJ_82. The Sanji chromosome also harbors a type 2 secretion system (T2SS), a type 3 secretion system (T3SS), a type 6 secretion system (T6SS), and several putative prophages. Sanji and other ST167 strains have a previously uncharacterized O-antigen (O89b) that is most closely related to serotype O89 as determined on the basis of analysis of the wzm-wzt genes and in silico serotyping. This O89b-antigen gene cluster was also found in the genomes of a few other pathogenic sequence type 617 (ST617) and ST10 complex strains. A time-scaled phylogeny inferred from comparative single nucleotide variant analysis indicated that development of these O89b-containing lineages emerged about 30 years ago. Comparative sequence analysis revealed that the core genome of Sanji is nearly identical to that of several recently sequenced strains of pathogenic XDR E. coli belonging to the ST167 group. Comparison of the mobile elements among the different ST167 genomes revealed that each genome carries a distinct set of multidrug resistance genes on different types of plasmids, indicating that there are multiple paths toward the emergence of XDR in E. coli.

IMPORTANCEE. coli strain Sanji is the first sequenced and analyzed genome of the recently emerged pathogenic XDR strains with sequence type ST167 and novel in silico serotype O89b:H9. Comparison of the genomes of Sanji with other ST167 strains revealed distinct sets of different plasmids, mobile IS elements, and antibiotic resistance genes in each genome, indicating that there exist multiple paths toward achieving XDR. The emergence of these pathogenic ST167 E. coli strains with diverse XDR capabilities highlights the difficulty of preventing or mitigating the development of XDR properties in bacteria and points to the importance of better understanding of the shared underlying virulence mechanisms and physiology of pathogenic bacteria.

Author Video: An author video summary of this article is available.

Antagonistic Pleiotropy in the Bifunctional Surface Protein FadL (OmpP1) during Adaptation of Haemophilus influenzae to Chronic Lung Infection Associated with Chronic Obstructive Pulmonary Disease

Tracking bacterial evolution during chronic infection provides insights into how host selection pressures shape bacterial genomes. The human-restricted opportunistic pathogen nontypeable Haemophilus influenzae (NTHi) infects the lower airways of patients suffering chronic obstructive pulmonary disease (COPD) and contributes to disease progression. To identify bacterial genetic variation associated with bacterial adaptation to the COPD lung, we sequenced the genomes of 92 isolates collected from the sputum of 13 COPD patients over 1 to 9 years. Individuals were colonized by distinct clonal types (CTs) over time, but the same CT was often reisolated at a later time or found in different patients. Although genomes from the same CT were nearly identical, intra-CT variation due to mutation and recombination occurred. Recurrent mutations in several genes were likely involved in COPD lung adaptation. Notably, nearly a third of CTs were polymorphic for null alleles of ompP1 (also called fadL), which encodes a bifunctional membrane protein that both binds the human carcinoembryonic antigen-related cell adhesion molecule 1 (hCEACAM1) receptor and imports long-chain fatty acids (LCFAs). Our computational studies provide plausible three-dimensional models for FadL's interaction with hCEACAM1 and LCFA binding. We show that recurrent fadL mutations are likely a case of antagonistic pleiotropy, since loss of FadL reduces NTHi's ability to infect epithelia but also increases its resistance to bactericidal LCFAs enriched within the COPD lung. Supporting this interpretation, truncated fadL alleles are common in publicly available NTHi genomes isolated from the lower airway tract but rare in others. These results shed light on molecular mechanisms of bacterial pathoadaptation and guide future research toward developing novel COPD therapeutics.IMPORTANCE Nontypeable Haemophilus influenzae is an important pathogen in patients with chronic obstructive pulmonary disease (COPD). To elucidate the bacterial pathways undergoing in vivo evolutionary adaptation, we compared bacterial genomes collected over time from 13 COPD patients and identified recurrent genetic changes arising in independent bacterial lineages colonizing different patients. Besides finding changes in phase-variable genes, we found recurrent loss-of-function mutations in the ompP1 (fadL) gene. We show that loss of OmpP1/FadL function reduces this bacterium's ability to infect cells via the hCEACAM1 epithelial receptor but also increases its resistance to bactericidal fatty acids enriched within the COPD lung, suggesting a case of antagonistic pleiotropy that restricts ΔfadL strains' niche. These results show how H. influenzae adapts to host-generated inflammatory mediators in the COPD airways.

The involvement of the low-oxygen-activated locus of Burkholderia cenocepacia in adaptation during cystic fibrosis infection

Chronic infection with opportunistic pathogens including Burkholderia cepacia complex (Bcc) is a hallmark of cystic fibrosis (CF). We investigated the adaptive mechanisms facilitating chronic lung infection in sequential Bcc isolates from two siblings with CF (P1 and P2), one of whom also experienced intermittent blood-stream infections (P2). We previously showed increased lung cell attachment with colonisation time in both P1 and P2. WGS analysis confirmed that the isolates are closely related. Twelve genes showed three or more mutations, suggesting these were genes under selection. Single nucleotide polymorphisms (SNVs) in 45 regulatory genes were also observed. Proteomic analysis showed that the abundance of 149 proteins increased over 61-months in sputum isolates, and both time- and source-related alterations in protein abundance between the second patient’s isolates. A consistent time-dependent increase in abundance of 19 proteins encoded by a low-oxygen-activated (lxa) locus was observed in both sets of isolates. Attachment was dramatically reduced in a B. cenocepacia K56-2Δlxa-locus deletion mutant, further indicating that it encodes protein(s) involved in host-cell attachment. Time-related changes in virulence in Galleria mellonella or motility were not observed. We conclude that the lxa-locus, associated with anoxic persistence in vitro, plays a role in host-cell attachment and adaptation to chronic colonization in the hypoxic niche of the CF lung.

The OmpR Regulator of Burkholderia multivorans Controls Mucoid-to-Nonmucoid Transition and Other Cell Envelope Properties Associated with Persistence in the Cystic Fibrosis Lung

Bacteria from the Burkholderia cepacia complex grow in different natural and man-made environments and are feared opportunistic pathogens that cause chronic respiratory infections in cystic fibrosis patients. Previous studies showed that Burkholderia mucoid clinical isolates grown under stress conditions give rise to nonmucoid variants devoid of the exopolysaccharide cepacian. Here, we determined that a major cause of the nonmucoid morphotype involves nonsynonymous mutations and small indels in the ompR gene encoding a response regulator of a two-component regulatory system. In trans complementation of nonmucoid variants (NMVs) with the native gene restored exopolysaccharide production. The loss of functional Burkholderia multivorans OmpR had positive effects on growth, adhesion to lung epithelial cells, and biofilm formation in high-osmolarity medium, as well as an increase in swimming and swarming motilities. In contrast, phenotypes such as antibiotic resistance, biofilm formation at low osmolarity, and virulence in Galleria mellonella were compromised by the absence of functional OmpR. Transcriptomic studies indicated that loss of the ompR gene affects the expression of 701 genes, many associated with outer membrane composition, motility, stress response, iron acquisition, and the uptake of nutrients, consistent with starvation tolerance. Since the stresses here imposed on B. multivorans may strongly resemble the ones found in the cystic fibrosis (CF) airways and mutations in the ompR gene from longitudinally collected CF isolates have been found, this regulator might be important for the production of NMVs in the CF environment.IMPORTANCE Within the cystic fibrosis (CF) lung, bacteria experience high-osmolarity conditions due to an ion unbalance resulting from defects in CF transmembrane conductance regulator (CFTR) protein activity in epithelial cells. Understanding how bacterial CF pathogens thrive in this environment might help the development of new therapeutic interventions to prevent chronic respiratory infections. Here, we show that the OmpR response regulator of one of the species found in CF respiratory infections, Burkholderia multivorans, is involved in the emergence of nonmucoid colony variants and is important for osmoadaptation by regulating several cell envelope components. Specifically, genetic, phenotypic, genomic, and transcriptomic approaches uncover OmpR as a regulator of cell wall remodeling under stress conditions, with implications in several phenotypes such as exopolysaccharide production, motility, antibiotic resistance, adhesion, and virulence.

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.

A c-di-GMP-Modulating Protein Regulates Swimming Motility of Burkholderia cenocepacia in Response to Arginine and Glutamate

Burkholderia cenocepacia is an opportunistic bacterium that can thrive in different environments, including the amino acid-rich mucus of the cystic fibrosis (CF) lung. B. cenocepacia responds to the nutritional conditions that mimic the CF sputum by increasing flagellin expression and swimming motility. Individual amino acids also induce swimming but not flagellin expression. Here, we show that modulation of the second messenger cyclic dimeric guanosine monophosphate (c-di-GMP) levels by the PAS-containing c-di-GMP phosphodiesterase, BCAL1069 (CdpA), regulates the swimming motility of B. cenocepacia K56-2 in response to CF sputum nutritional conditions. Heterologous expression of WspR, a diguanylate cyclase, in B. cenocepacia K56-2 caused an increase in c-di-GMP levels and reduced swimming motility but did not affect flagellin expression or flagellar biosynthesis. After insertional mutagenesis of 12 putative genes encoding c-di-GMP metabolizing enzymes, one mutant of the locus BCAL1069 (cdpA), exhibited decreased swimming motility independent of flagellin expression in CF sputum nutritional conditions and an increase in intracellular c-di-GMP levels. The reduced swimming motility phenotype of the BCAL1069 mutant was observed in the presence of arginine and glutamate, but not of histidine, phenylalanine, or proline. The B. cenocepacia CdpA was also found to be involved in regulation of protease activity but not in biofilm formation. Altogether, these results highlight a role of B. cenocepacia BCAL1069 (CdpA) in sensing the nutritional conditions of the CF sputum and eliciting a pathogenic response that includes swimming motility toward amino acids and an increase in protease activity.

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.

An integrative strategy to identify the entire protein coding potential of prokaryotic genomes by proteogenomics

Accurate annotation of all protein-coding sequences (CDSs) is an essential prerequisite to fully exploit the rapidly growing repertoire of completely sequenced prokaryotic genomes. However, large discrepancies among the number of CDSs annotated by different resources, missed functional short open reading frames (sORFs), and overprediction of spurious ORFs represent serious limitations. Our strategy toward accurate and complete genome annotation consolidates CDSs from multiple reference annotation resources, ab initio gene prediction algorithms and in silico ORFs (a modified six-frame translation considering alternative start codons) in an integrated proteogenomics database (iPtgxDB) that covers the entire protein-coding potential of a prokaryotic genome. By extending the PeptideClassifier concept of unambiguous peptides for prokaryotes, close to 95% of the identifiable peptides imply one distinct protein, largely simplifying downstream analysis. Searching a comprehensive Bartonella henselae proteomics data set against such an iPtgxDB allowed us to unambiguously identify novel ORFs uniquely predicted by each resource, including lipoproteins, differentially expressed and membrane-localized proteins, novel start sites and wrongly annotated pseudogenes. Most novelties were confirmed by targeted, parallel reaction monitoring mass spectrometry, including unique ORFs and single amino acid variations (SAAVs) identified in a re-sequenced laboratory strain that are not present in its reference genome. We demonstrate the general applicability of our strategy for genomes with varying GC content and distinct taxonomic origin. We release iPtgxDBs for B. henselae, Bradyrhizobium diazoefficiens and Escherichia coli and the software to generate both proteogenomics search databases and integrated annotation files that can be viewed in a genome browser for any prokaryote.

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.

Structure of O-Antigen and Hybrid Biosynthetic Locus in Burkholderia cenocepacia Clonal Variants Recovered from a Cystic Fibrosis Patient

Burkholderia cenocepacia is an opportunistic pathogen associated with chronic lung infections and increased risk of death in patients with cystic fibrosis (CF). In this work, we investigated the lipopolysaccharide (LPS) of clinical variants of B. cenocepacia that were collected from a CF patient over a period of 3.5 years, from the onset of infection until death by necrotizing pneumonia (cepacia syndrome). We report the chemical structure of the LPS molecule of various sequential isolates and the identification of a novel hybrid O-antigen (OAg) biosynthetic cluster. The OAg repeating unit of the LPS from IST439, the initial isolate, is a [→2)-β-D-Ribf-(1→4)-α-D-GalpNAc-(1→] disaccharide, which was not previously described in B. cenocepacia. The IST439 OAg biosynthetic gene cluster contains 7 of 23 genes that are closely homologous to genes found in B. multivorans, another member of the Burkholderia cepacia complex. None of the subsequent isolates expressed OAg. Genomic sequencing of these isolates enabled the identification of mutations within the OAg cluster, but none of these mutations could be associated with the loss of OAg. This study provides support to the notion that OAg LPS modifications are an important factor in the adaptation of B. cenocepacia to chronic infection and that the heterogeneity of OAgs relates to variation within the OAg gene cluster, indicating that the gene cluster might have been assembled through multiple horizontal transmission events.