Mixed strain pathogen populations accelerate the evolution of antibiotic resistance in patients.

Antibiotic resistance poses a global health threat, but the within-host drivers of resistance remain poorly understood. Pathogen populations are often assumed to be clonal within hosts, and resistance is thought to emerge due to selection for de novo variants. Here we show that mixed strain populations are common in the opportunistic pathogen P. aeruginosa. Crucially, resistance evolves rapidly in patients colonized by multiple strains through selection for pre-existing resistant strains. In contrast, resistance evolves sporadically in patients colonized by single strains due to selection for novel resistance mutations. However, strong trade-offs between resistance and growth rate occur in mixed strain populations, suggesting that within-host diversity can also drive the loss of resistance in the absence of antibiotic treatment. In summary, we show that the within-host diversity of pathogen populations plays a key role in shaping the emergence of resistance in response to treatment.

Rapid evolution and host immunity drive the rise and fall of carbapenem resistance during an acute Pseudomonas aeruginosa infection.

It is well established that antibiotic treatment selects for resistance, but the dynamics of this process during infections are poorly understood. Here we map the responses of Pseudomonas aeruginosa to treatment in high definition during a lung infection of a single ICU patient. Host immunity and antibiotic therapy with meropenem suppressed P. aeruginosa, but a second wave of infection emerged due to the growth of oprD and wbpM meropenem resistant mutants that evolved in situ. Selection then led to a loss of resistance by decreasing the prevalence of low fitness oprD mutants, increasing the frequency of high fitness mutants lacking the MexAB-OprM efflux pump, and decreasing the copy number of a multidrug resistance plasmid. Ultimately, host immunity suppressed wbpM mutants with high meropenem resistance and fitness. Our study highlights how natural selection and host immunity interact to drive both the rapid rise, and fall, of resistance during infection.

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.

Structural and biochemical approaches uncover multiple evolutionary trajectories of plant quinate dehydrogenases.

Quinate is produced and used by many plants in the biosynthesis of chlorogenic acids (CGAs). Chlorogenic acids are astringent and serve to deter herbivory. They also function as antifungal agents and have potent antioxidant properties. Quinate is produced at a branch point of shikimate biosynthesis by the enzyme quinate dehydrogenase (QDH). However, little information exists on the identity and biochemical properties of plant QDHs. In this study, we utilized structural and bioinformatics approaches to establish a QDH-specific primary sequence motif. Using this motif, we identified QDHs from diverse plants and confirmed their activity by recombinant protein production and kinetic assays. Through a detailed phylogenetic analysis, we show that plant QDHs arose directly from bifunctional dehydroquinate dehydratase-shikimate dehydrogenases (DHQD-SDHs) through different convergent evolutionary events, illustrated by our findings that eudicot and conifer QDHs arose early in vascular plant evolution whereas Brassicaceae QDHs emerged later. This process of recurrent evolution of QDH is further demonstrated by the fact that this family of proteins independently evolved NAD+ and NADP+ specificity in eudicots. The acquisition of QDH activity by these proteins was accompanied by the inactivation or functional evolution of the DHQD domain, as verified by enzyme activity assays and as reflected in the loss of key DHQD active site residues. The implications of QDH activity and evolution are discussed in terms of plant growth and development.

Identification and analysis of seven effector protein families with different adaptive and evolutionary histories in plant-associated members of the Xanthomonadaceae.

The Xanthomonadaceae family consists of species of non-pathogenic and pathogenic γ-proteobacteria that infect different hosts, including humans and plants. In this study, we performed a comparative analysis using 69 fully sequenced genomes belonging to this family, with a focus on identifying proteins enriched in phytopathogens that could explain the lifestyle and the ability to infect plants. Using a computational approach, we identified seven phytopathogen-enriched protein families putatively secreted by type II secretory system: PheA (CM-sec), LipA/LesA, VirK, and four families involved in N-glycan degradation, NixE, NixF, NixL, and FucA1. In silico and phylogenetic analyses of these protein families revealed they all have orthologs in other phytopathogenic or symbiotic bacteria, and are involved in the modulation and evasion of the immune system. As a proof of concept, we performed a biochemical characterization of LipA from Xac306 and verified that the mutant strain lost most of its lipase and esterase activities and displayed reduced virulence in citrus. Since this study includes closely related organisms with distinct lifestyles and highlights proteins directly related to adaptation inside plant tissues, novel approaches might use these proteins as biotechnological targets for disease control, and contribute to our understanding of the coevolution of plant-associated bacteria.

Selective Sweeps and Parallel Pathoadaptation Drive Pseudomonas aeruginosa Evolution in the Cystic Fibrosis Lung.

UNLABELLED: Pulmonary infections caused by Pseudomonas aeruginosa are a recalcitrant problem in cystic fibrosis (CF) patients. While the clinical implications and long-term evolutionary patterns of these infections are well studied, we know little about the short-term population dynamics that enable this pathogen to persist despite aggressive antimicrobial therapy. Here, we describe a short-term population genomic analysis of 233 P. aeruginosa isolates collected from 12 sputum specimens obtained over a 1-year period from a single patient. Whole-genome sequencing and antimicrobial susceptibility profiling identified the expansion of two clonal lineages. The first lineage originated from the coalescence of the entire sample less than 3 years before the end of the study and gave rise to a high-diversity ancestral population. The second expansion occurred 2 years later and gave rise to a derived population with a strong signal of positive selection. These events show characteristics consistent with recurrent selective sweeps. While we cannot identify the specific mutations responsible for the origins of the clonal lineages, we find that the majority of mutations occur in loci previously associated with virulence and resistance. Additionally, approximately one-third of all mutations occur in loci that are mutated multiple times, highlighting the importance of parallel pathoadaptation. One such locus is the gene encoding penicillin-binding protein 3, which received three independent mutations. Our functional analysis of these alleles shows that they provide differential fitness benefits dependent on the antibiotic under selection. These data reveal that bacterial populations can undergo extensive and dramatic changes that are not revealed by lower-resolution analyses. IMPORTANCE: Pseudomonas aeruginosa is a bacterial opportunistic pathogen responsible for significant morbidity and mortality in cystic fibrosis (CF) patients. Once it has colonized the lung in CF, it is highly resilient and rarely eradicated. This study presents a deep sampling examination of the fine-scale evolutionary dynamics of P. aeruginosa in the lungs of a chronically infected CF patient. We show that diversity of P. aeruginosa is driven by recurrent clonal emergence and expansion within this patient and identify potential adaptive variants associated with these events. This high-resolution sequencing strategy thus reveals important intraspecies dynamics that explain a clinically important phenomenon not evident at a lower-resolution analysis of community structure.

Phenotypic diversity within a Pseudomonas aeruginosa population infecting an adult with cystic fibrosis.

Chronic airway infections caused by Pseudomonas aeruginosa contribute to the progression of pulmonary disease in individuals with cystic fibrosis (CF). In the setting of CF, within-patient adaptation of a P. aeruginosa strain generates phenotypic diversity that can complicate microbiological analysis of patient samples. We investigated within- and between- sample diversity of 34 phenotypes among 235 P. aeruginosa isolates cultured from sputum samples collected from a single CF patient over the span of one year, and assessed colony morphology as a screening tool for predicting phenotypes, including antimicrobial susceptibilities. We identified 15 distinct colony morphotypes that varied significantly in abundance both within and between sputum samples. Substantial within sample phenotypic heterogeneity was also noted in other phenotypes, with morphotypes being unreliable predictors of antimicrobial susceptibility and other phenotypes. Emergence of isolates with reduced susceptibility to ß-lactams was observed during periods of clinical therapy with aztreonam. Our findings confirm that the P. aeruginosa population in chronic CF lung infections is highly dynamic, and that intra-sample phenotypic diversity is underestimated if only one or few colonies are analyzed per sample.

Lung microbiota across age and disease stage in cystic fibrosis.

Understanding the significance of bacterial species that colonize and persist in cystic fibrosis (CF) airways requires a detailed examination of bacterial community structure across a broad range of age and disease stage. We used 16S ribosomal RNA sequencing to characterize the lung microbiota in 269 CF patients spanning a 60 year age range, including 76 pediatric samples from patients of age 4-17, and a broad cross-section of disease status to identify features of bacterial community structure and their relationship to disease stage and age. The CF lung microbiota shows significant inter-individual variability in community structure, composition and diversity. The core microbiota consists of five genera - Streptococcus, Prevotella, Rothia, Veillonella and Actinomyces. CF-associated pathogens such as Pseudomonas, Burkholderia, Stenotrophomonas and Achromobacter are less prevalent than core genera, but have a strong tendency to dominate the bacterial community when present. Community diversity and lung function are greatest in patients less than 10 years of age and lower in older age groups, plateauing at approximately age 25. Lower community diversity correlates with worse lung function in a multivariate regression model. Infection by Pseudomonas correlates with age-associated trends in community diversity and lung function.

Analysis of the cystic fibrosis lung microbiota via serial Illumina sequencing of bacterial 16S rRNA hypervariable regions.

The characterization of bacterial communities using DNA sequencing has revolutionized our ability to study microbes in nature and discover the ways in which microbial communities affect ecosystem functioning and human health. Here we describe Serial Illumina Sequencing (SI-Seq): a method for deep sequencing of the bacterial 16S rRNA gene using next-generation sequencing technology. SI-Seq serially sequences portions of the V5, V6 and V7 hypervariable regions from barcoded 16S rRNA amplicons using an Illumina short-read genome analyzer. SI-Seq obtains taxonomic resolution similar to 454 pyrosequencing for a fraction of the cost, and can produce hundreds of thousands of reads per sample even with very high multiplexing. We validated SI-Seq using single species and mock community controls, and via a comparison to cystic fibrosis lung microbiota sequenced using 454 FLX Titanium. Our control runs show that SI-Seq has a dynamic range of at least five orders of magnitude, can classify >96% of sequences to the genus level, and performs just as well as 454 and paired-end Illumina methods in estimation of standard microbial ecology diversity measurements. We illustrate the utility of SI-Seq in a pilot sample of central airway secretion samples from cystic fibrosis patients.