Whole genome sequencing of bacteria in cystic fibrosis as a model for bacterial genome adaptation and evolution

Staphylococcus aureus and Cystic Fibrosis-A Close Relationship. What Can We Learn from Sequencing Studies?

Staphylococcus aureus is next to Pseudomonas aeruginosa the most isolated pathogen from the airways of cystic fibrosis (CF) patients, who are often infected by a dominant S. aureus clone for extended periods. To be able to persist, the pathogen has to adapt to the hostile niche of the airways to counteract host defence, antibiotic therapy and the competition with coinfecting pathogens. S. aureus is equipped with many virulence factors including adhesins, toxins that are localized on the chromosome, on plasmids or are phage-related. S. aureus is especially versatile and adaptation and evolution of the pathogen occurs by the acquisition of new genes by horizontal gene transfer (HGT), changes in nucleotides (single nucleotide variations, SNVs) that can cause a selective advantage for the bacteria and become fixed in subpopulations. Methicillin-resistant S. aureus are a special threat to CF patients due to the more severe lung disease occurring in infected patients. Today, with decreasing costs for sequencing, more and more studies using S. aureus isolates cultured from CF patients are being published, which use whole genome sequencing (WGS), multilocus sequence typing (MLST) or spa-sequence typing (spa-typing) to follow the population dynamics of S. aureus, elucidate the underlying mechanisms of phenotypic variants, newly acquired resistance or adaptation to the host response in this particular niche. In the first part of this review, an introduction to the genetic make-up and the pathogenesis of S. aureus with respect to CF is provided. The second part presents an overview of recent studies and their findings using genotypic methods such as single or multilocus sequencing and whole genome sequencing, which identify factors contributing to the adaptation of S. aureus and its evolution in the airways of individuals with CF.

Aspergillus fumigatus In-Host HOG Pathway Mutation for Cystic Fibrosis Lung Microenvironment Persistence

The prevalence of Aspergillus fumigatus colonization in individuals with cystic fibrosis (CF) and subsequent fungal persistence in the lung is increasingly recognized. However, there is no consensus for clinical management of A. fumigatus in CF individuals, due largely to uncertainty surrounding A. fumigatus CF pathogenesis and virulence mechanisms. To address this gap in knowledge, a longitudinal series of A. fumigatus isolates from an individual with CF were collected over 4.5 years. Isolate genotypes were defined with whole-genome sequencing that revealed both transitory and persistent A. fumigatus in the lung. Persistent lineage isolates grew most readily in a low-oxygen culture environment, and conidia were more sensitive to oxidative stress-inducing conditions than those from nonpersistent isolates. Closely related persistent isolates harbored a unique allele of the high-osmolarity glycerol (HOG) pathway mitogen-activated protein kinase kinase, Pbs2 (pbs2C2). Data suggest this novel pbs2C2 allele arose in vivo and is necessary for the fungal response to osmotic stress in a low-oxygen environment through hyperactivation of the HOG (SakA) signaling pathway. Hyperactivation of the HOG pathway through pbs2C2 comes at the cost of decreased conidial stress resistance in the presence of atmospheric oxygen levels. These novel findings shed light on pathoadaptive mechanisms of A. fumigatus in CF, lay the foundation for identifying persistent A. fumigatus isolates that may require antifungal therapy, and highlight considerations for successful culture of persistent Aspergillus CF isolates. IMPORTANCE Aspergillus fumigatus infection causes a spectrum of clinical manifestations. For individuals with cystic fibrosis (CF), allergic bronchopulmonary aspergillosis (ABPA) is an established complication, but there is a growing appreciation for A. fumigatus airway persistence in CF disease progression. There currently is little consensus for clinical management of A. fumigatus long-term culture positivity in CF. A better understanding of A. fumigatus pathogenesis mechanisms in CF is expected to yield insights into when antifungal therapies are warranted. Here, a 4.5-year longitudinal collection of A. fumigatus isolates from a patient with CF identified a persistent lineage that harbors a unique allele of the Pbs2 mitogen-activated protein kinase kinase (MAPKK) necessary for unique CF-relevant stress phenotypes. Importantly for A. fumigatus CF patient diagnostics, this allele provides increased fitness under CF lung-like conditions at a cost of reduced in vitro growth under standard laboratory conditions. These data illustrate a molecular mechanism for A. fumigatus CF lung persistence with implications for diagnostics and antifungal therapy.

From genotype to phenotype: adaptations of Pseudomonas aeruginosa to the cystic fibrosis environment

Pseudomonas aeruginosa is one of the main microbial species colonizing the lungs of cystic fibrosis patients and is responsible for the decline in respiratory function. Despite the hostile pulmonary environment, P. aeruginosa is able to establish chronic infections thanks to its strong adaptive capacity. Various longitudinal studies have attempted to compare the strains of early infection with the adapted strains of chronic infection. Thanks to new '-omics' techniques, convergent genetic mutations, as well as transcriptomic and proteomic dysregulations have been identified. As a consequence of this evolution, the adapted strains of P. aeruginosa have particular phenotypes that promote persistent infection.

Genome-based classification of Burkholderia cepacia complex provides new insight into its taxonomic status

Background

Accurate classification of different Burkholderia cepacia complex (BCC) species is essential for therapy, prognosis assessment and research. The taxonomic status of BCC remains problematic and an improved knowledge about the classification of BCC is in particular needed.

Methods

We compared phylogenetic trees of BCC based on 16S rRNA, recA, hisA and MLSA (multilocus sequence analysis). Using the available whole genome sequences of BCC, we inferred a species tree based on estimated single-copy orthologous genes and demarcated species of BCC using dDDH/ANI clustering.

Results

We showed that 16S rRNA, recA, hisA and MLSA have limited resolutions in the taxonomic study of closely related bacteria such as BCC. Our estimated species tree and dDDH/ANI clustering clearly separated 116 BCC strains into 36 clusters. With the appropriate reclassification of misidentified strains, these clusters corresponded to 22 known species as well as 14 putative novel species.

Conclusions

This is the first large-scale and systematic study of the taxonomic status of the BCC and could contribute to further insights into BCC taxonomy. Our study suggested that conjunctive use of core phylogeny based on single-copy orthologous genes, as well as pangenome-based dDDH/ANI clustering would provide a preferable framework for demarcating closely related species.

Reviewer

This article was reviewed by Dr. Xianwen Ren.

Complete nontuberculous mycobacteria whole genomes using an optimized DNA extraction protocol for long-read sequencing

Background

Nontuberculous mycobacteria (NTM) are a major cause of pulmonary and systemic disease in at-risk populations. Gaps in knowledge about transmission patterns, evolution, and pathogenicity during infection have prompted a recent surge in genomic NTM research. Increased availability and affordability of whole genome sequencing (WGS) techniques provide new opportunities to sequence and construct complete bacterial genomes faster and at a lower cost. However, extracting large quantities of pure genomic DNA is particularly challenging with NTM due to its slow growth and recalcitrant cell wall. Here we report a DNA extraction protocol that is optimized for long-read WGS of NTM, yielding large quantities of highly pure DNA with no additional clean-up steps.

Results

Our DNA extraction method was compared to 6 other methods with variations in timing of mechanical disruption and enzymatic digestion of the cell wall, quantity of matrix material, and reagents used in extraction and precipitation. We tested our optimized method on 38 clinical isolates from the M. avium and M. abscessus complexes, which yielded optimal quality and quantity measurements for Oxford Nanopore Technologies sequencing. We also present the efficient completion of circularized M. avium subspecies hominissuis genomes using our extraction technique and the long-read sequencing MinION platform, including the identification of a novel plasmid.

Conclusions

Our optimized extraction protocol and assembly pipeline was both sufficient and efficient for genome closure. We expect that our finely-tuned extraction method will prove to be a valuable tool in long-read sequencing and completion of mycobacterial genomes going forward. Utilization of comprehensive, long-read based approaches will advance the understanding evolution and pathogenicity of NTM infections.

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

Chronic bacterial infections of the lung are the leading cause of morbidity and mortality in cystic fibrosis patients. Tracking bacterial evolution during chronic infections can provide insights into how host selection pressures—including immune responses and therapeutic interventions—shape bacterial genomes. We carried out genomic and phenotypic analyses of 215 serially collected Burkholderia cenocepacia isolates from 16 cystic fibrosis patients, spanning a period of 2–20 yr and a broad range of epidemic lineages. Systematic phenotypic tests identified longitudinal bacterial series that manifested progressive changes in liquid media growth, motility, biofilm formation, and acute insect virulence, but not in mucoidy. The results suggest that distinct lineages follow distinct evolutionary trajectories during lung infection. Pan-genome analysis identified 10,110 homologous gene clusters present only in a subset of strains, including genes restricted to different molecular types. Our phylogenetic analysis based on 2148 orthologous gene clusters from all isolates is consistent with patient-specific clades. This suggests that initial colonization of patients was likely by individual strains, followed by subsequent diversification. Evidence of clonal lineages shared by some patients was observed, suggesting inter-patient transmission. We observed recurrent gene losses in multiple independent longitudinal series, including complete loss of Chromosome III and deletions on other chromosomes. Recurrently observed loss-of-function mutations were associated with decreases in motility and biofilm formation. Together, our study provides the first comprehensive genome-phenome analyses of B. cenocepacia infection in cystic fibrosis lungs and serves as a valuable resource for understanding the genomic and phenotypic underpinnings of bacterial evolution.

Variation of Burkholderia cenocepacia virulence potential during cystic fibrosis chronic lung infection

During long-term lung infection in cystic fibrosis (CF) patients, Burkholderia cenocepacia faces multiple selective pressures in this highly stressful and fluctuating environment. As a consequence, the initial infecting strain undergoes genetic changes that result in the diversification of genotypes and phenotypes. Whether this clonal expansion influences the pathogenic potential is unclear. The virulence potential of 39 sequential B. cenocepacia (recA lineage IIIA) isolates, corresponding to 3 different clones retrieved from 3 chronically infected CF patients was compared in this study using the non-mammalian infection hosts Galleria mellonella and Caenorhabditis elegans. The isolates used in this retrospective study were picked randomly from selective agar plates as part of a CF Center routine, from the onset of infection until patients' death after 3.5 and 7.5 y or the more recent isolation date after 12.5 y of chronic infection. The infection models proved useful to assess virulence potential diversification, but for some isolates the relative values diverged in C. elegans and G. mellonella. Results also reinforce the concept of the occurrence of clonal diversification and co-existence of multiple phenotypes within the CF lungs, also with respect to pathogenicity. No clear trend of decrease (or increase) of the virulence potential throughout long-term infection was found but there is an apparent tendency for a clone/patient-dependent decrease of virulence when the G. mellonella model was used. The sole avirulent variant in both infection hosts was found to lack the small third replicon previously associated to virulence. Although possible, the in vivo loss of this nonessential megaplasmid was found to be a rare event (1 among a total of 64 isolates examined).

Specific Gene Loci of Clinical Pseudomonas putida Isolates

Pseudomonas putida are ubiquitous inhabitants of soils and clinical isolates of this species have been seldom described. Clinical isolates show significant variability in their ability to cause damage to hosts because some of them are able to modulate the host’s immune response. In the current study, comparisons between the genomes of different clinical and environmental strains of P. putida were done to identify genetic clusters shared by clinical isolates that are not present in environmental isolates. We show that in clinical strains specific genes are mostly present on transposons, and that this set of genes exhibit high identity with genes found in pathogens and opportunistic pathogens. The set of genes prevalent in P. putida clinical isolates, and absent in environmental isolates, are related with survival under oxidative stress conditions, resistance against biocides, amino acid metabolism and toxin/antitoxin (TA) systems. This set of functions have influence in colonization and survival within human tissues, since they avoid host immune response or enhance stress resistance. An in depth bioinformatic analysis was also carried out to identify genetic clusters that are exclusive to each of the clinical isolates and that correlate with phenotypical differences between them, a secretion system type III-like was found in one of these clinical strains, a determinant of pathogenicity in Gram-negative bacteria.