Noroviruses are a major cause of acute gastroenteritis worldwide and can establish chronic infection in immunocompromised individuals. To investigate the mechanisms of norovirus evolution during chronic infection, we selected seven representative patients from a National Institutes of Health study cohort who sustained norovirus infection for periods ranging from 73 to 1,492 days. Six patients shed viruses belonging to a single genotype (GII.2[PNA], GII.4 New Orleans[P4], GII.4 Den Haag[P4], GII.3[P21], GII.6[P7], or GII.14[P7]) over the period examined, while one patient sequentially shed two genotypes (GII.6[P7] followed by GII.4 Sydney[P31]). Norovirus genomes from consecutive stool samples were sequenced at high resolution (>3,300 reads/nucleotide position) using the Illumina platform and subjected to bioinformatics analysis. Norovirus sequences could be resolved into one or more discrete clonal RNA genomes that persisted within these patients over time. Phylogenetic analyses inferred that clonal populations originated from a single founder virus and not by reinfection with community strains. Estimated evolutionary rates of clonal populations during persistent infection were similar to those of noroviruses from acute infection in the global database, suggesting that inherently higher RNA-dependent polymerase error rates were not associated with the ability to persist. The high-resolution analysis of norovirus diversity and evolution at the population level described here should allow a better understanding of adaptive mutations sustained during chronic infection. IMPORTANCE Noroviruses are an important cause of chronic diarrhea in patients with compromised immune systems. Presently, there are no effective therapies to clear the virus, which can persist for years in the intestinal tract. The goal of our study was to develop a better understanding of the norovirus strains that are associated with these long-term infections. With the remarkable diversity of norovirus strains detected in the immunocompromised patient cohort we studied, it appears that most, if not all, noroviruses circulating in nature may have the capacity to establish a chronic infection when a person is unable to mount an effective immune response. Our work is the most comprehensive genetic data set generated to date in which near full-length genomes from noroviruses associated with chronic infection were analyzed by high-resolution next-generation sequencing. Analysis of this data set led to our discovery that certain patients in our cohort were shedding noroviruses that could be subdivided into distinct haplotypes or populations of viruses that were co-evolving independently. The ability to track haplotypes of noroviruses during chronic infection will allow us to fine-tune our understanding of how the virus adapts and maintains itself in the human host, and how selective pressures such as antiviral drugs can affect these distinct populations.

Antibiotic use has been linked to changes in the population structure of human pathogens and the clonal expansion of multidrug-resistant (MDR) strains among healthcare- and community-acquired infections. Here we present a compelling example in a veterinary pathogen, Rhodococcus equi, the causative agent of a severe pulmonary infection affecting foals worldwide. We show that the erm(46) gene responsible for emerging macrolide resistance among equine R. equi isolates in the United States is part of a 6.9-kb transposable element, TnRErm46, actively mobilized by an IS481 family transposase. TnRErm46 is carried on an 87-kb conjugative plasmid, pRErm46, transferable between R. equi strains at frequencies up to 10-3 The erm(46) gene becomes stabilized in R. equi by pRErm46's apparent fitness neutrality and wholesale TnRErm46 transposition onto the host genome. This includes the conjugally exchangeable pVAPA virulence plasmid, enabling the possibility of cotransfer of two essential traits for survival in macrolide-treated foals in a single mating event. Despite its high horizontal transfer potential, phylogenomic analyses show that erm(46) is paradoxically confined to a specific R. equi clone, 2287. R. equi 2287 also carries a unique rpoBS531F mutation conferring high-level resistance to rifampin, systematically administered together with macrolides against rhodococcal pneumonia on equine farms. Our data illustrate that under sustained combination therapy, several independent "founder" genetic events are concurrently required for resistance, limiting not only its emergence but also, crucially, horizontal spread, ultimately determining multiresistance clonality.IMPORTANCE MDR clades arise upon acquisition of resistance traits, but the determinants of their clonal expansion remain largely undefined. Taking advantage of the unique features of Rhodococcus equi infection control in equine farms, involving the same dual antibiotic treatment since the 1980s (a macrolide and rifampin), this study sheds light into the determinants of multiresistance clonality and the importance of combination therapy in limiting the dissemination of mobile resistance elements. Clinically effective therapeutic alternatives against R. equi foal pneumonia are currently lacking, and the identified macrolide-rifampin MDR clone 2287 has serious implications. Still at early stages of evolution and local spread, R. equi 2287 may disseminate globally, posing a significant threat to the equine industry and, also, public health due to the risk of zoonotic transmission. The characterization of the 2287 clone and its resistance determinants will enable targeted surveillance and control interventions to tackle the emergence of MDR R. equi.

An outbreak of the fungal pathogen Cryptococcus gattii began in the Pacific Northwest (PNW) in the late 1990s. This outbreak consists of three clonal subpopulations: VGIIa/major, VGIIb/minor, and VGIIc/novel. Both VGIIa and VGIIc are unique to the PNW and exhibit increased virulence. In this study, we sequenced the genomes of isolates from these three groups, as well as global isolates, and analyzed a total of 53 isolates. We found that VGIIa/b/c populations show evidence of clonal expansion in the PNW. Whole-genome sequencing provided evidence that VGIIb originated in Australia, while VGIIa may have originated in South America, and these were likely independently introduced. Additionally, the VGIIa outbreak lineage may have arisen from a less virulent clade that contained a mutation in the MSH2 ortholog, but this appears to have reverted in the VGIIa outbreak strains, suggesting that a transient mutator phenotype may have contributed to adaptation and evolution of virulence in the PNW outbreak. PNW outbreak isolates share genomic islands, both between the clonal lineages and with global isolates, indicative of sexual recombination. This suggests that VGII C. gattii has undergone sexual reproduction, either bisexual or unisexual, in multiple locales contributing to the production of novel, virulent subtypes. We also found that the genomes of two basal VGII isolates from HIV(+) patients contain an introgression tract spanning three genes. Introgression substantially contributed to intra-VGII polymorphism and likely occurred through sexual reproduction with VGI. More broadly, these findings illustrate how both microevolution and sexual reproduction play central roles in the development of infectious outbreaks from avirulent or less virulent progenitors. Importance: Cryptococcus gattii is the causative agent responsible for ongoing infections in the Pacific Northwest of the United States and western Canada. The incidence of these infections increased dramatically in the 1990s and remains elevated. These infections are attributable to three clonal lineages of C. gattii, VGIIa, VGIIb, and VGIIc, with only VGIIa identified once previously in the Pacific Northwest prior to the start of the outbreak, albeit in a less virulent form. This study addresses the origin and emergence of this outbreak, using whole-genome sequencing and comparison of both outbreak and global isolates. We show that VGIIa arose mitotically from a less virulent clonal group, possibly via the action of a mutator phenotype, while VGIIb was likely introduced from Australia, and VGIIc appears to have emerged in the United States or in an undersampled locale via sexual reproduction. This work shows that multiple processes can contribute to the emergence of an outbreak.

Carbapenem-resistant Klebsiella pneumoniae (CRKp) is an urgent public health threat. Worldwide dissemination of CRKp has been largely attributed to clonal group (CG) 258. However, recent evidence indicates the global emergence of a CRKp CG307 lineage. Houston, TX, is the first large city in the United States with detected cocirculation of both CRKp CG307 and CG258. We sought to characterize the genomic and clinical factors contributing to the parallel endemic spread of CG258 and CG307. CRKp isolates were collected as part of the prospective, Consortium on Resistance against Carbapenems in Klebsiella and other Enterobacterales 2 (CRACKLE-2) study. Hybrid short-read and long-read genome assemblies were generated from 119 CRKp isolates (95 originated from Houston hospitals). A comprehensive characterization of phylogenies, gene transfer, and plasmid content with pan-genome analysis was performed on all CRKp isolates. Plasmid mating experiments were performed with CG307 and CG258 isolates of interest. Dissection of the accessory genomes suggested independent evolution and limited horizontal gene transfer between CG307 and CG258 lineages. CG307 contained a diverse repertoire of mobile genetic elements, which were shared with other non-CG258 K. pneumoniae isolates. Three unique clades of Houston CG307 isolates clustered distinctly from other global CG307 isolates, indicating potential selective adaptation of particular CG307 lineages to their respective geographical niches. CG307 strains were often isolated from the urine of hospitalized patients, likely serving as important reservoirs for genes encoding carbapenemases and extended-spectrum β-lactamases. Our findings suggest parallel cocirculation of high-risk lineages with potentially divergent evolution. IMPORTANCE The prevalence of carbapenem-resistant Klebsiella pneumoniae (CRKp) infections in nosocomial settings remains a public health challenge. High-risk clones such as clonal group 258 (CG258) are particularly concerning due to their association with blaKPC carriage, which can severely complicate antimicrobial treatments. There is a recent emergence of clonal group 307 (CG307) worldwide with little understanding of how this successful clone has been able to adapt while cocirculating with CG258. We provide the first evidence of potentially divergent evolution between CG258 and CG307 with limited sharing of adaptive genes. Houston, TX, is home to the largest medical center in the world, with a large influx of domestic and international patients. Thus, our unique geographical setting, where two pandemic strains of CRKp are circulating, provides an indication of how differential accessory genome content can drive stable, endemic populations of CRKp. Pan-genomic analyses such as these can reveal unique signatures of successful CRKp dissemination, such as the CG307-associated plasmid (pCG307_HTX), and provide invaluable insights into the surveillance of local carbapenem-resistant Enterobacterales (CRE) epidemiology.

Highly invasive, community-acquired Klebsiella pneumoniae infections have recently emerged, resulting in pyogenic liver abscesses. These infections are caused by hypervirulent K. pneumoniae (hvKP) isolates primarily of capsule serotype K1 or K2. Hypervirulent K1 isolates belong to clonal complex 23 (CC23), indicating that this clonal lineage has a specific genetic background conferring hypervirulence. Here, we apply whole-genome sequencing to a collection of K. pneumoniae isolates to characterize the phylogenetic background of hvKP isolates with an emphasis on CC23. Most of the hvKP isolates belonged to CC23 and grouped into a distinct monophyletic clade, revealing that CC23 is a unique clonal lineage, clearly distinct from nonhypervirulent strains. Separate phylogenetic analyses of the CC23 isolates indicated that the CC23 lineage evolved recently by clonal expansion from a single common ancestor. Limited grouping according to geographical origin was observed, suggesting that CC23 has spread globally through multiple international transmissions. Conversely, hypervirulent K2 strains clustered in genetically unrelated groups. Strikingly, homologues of a large virulence plasmid were detected in all hvKP clonal lineages, indicating a key role in K. pneumoniae hypervirulence. The plasmid encodes two siderophores, aerobactin and salmochelin, and RmpA (regulator of the mucoid phenotype); all these factors were found to be restricted to hvKP isolates. Genomic comparisons revealed additional factors specifically associated with CC23. These included a distinct variant of a genomic island encoding yersiniabactin, colibactin, and microcin E492. Furthermore, additional novel genomic regions unique to CC23 were revealed which may also be involved in the increased virulence of this important clonal lineage.

Homologous recombination (HR) enables the exchange of genetic material between and within species. Recent studies suggest that this process plays a major role in the microevolution of microbial genomes, contributing to core genome homogenization and to the maintenance of cohesive population structures. However, we still have a very poor understanding of the possible adaptive roles of intraspecific HR and of the factors that determine its differential impact across clades and lifestyles. Here we used a unified methodological framework to assess HR in 338 complete genomes from 54 phylogenetically diverse and representative prokaryotic species, encompassing different lifestyles and a broad phylogenetic distribution. Our results indicate that lifestyle and presence of restriction-modification (RM) machineries are among the main factors shaping HR patterns, with symbionts and intracellular pathogens having the lowest HR levels. Similarly, the size of exchanged genomic fragments correlated with the presence of RM and competence machineries. Finally, genes exchanged by HR showed functional enrichments which could be related to adaptations to different environments and ecological strategies. Taken together, our results clarify the factors underlying HR impact and suggest important adaptive roles of genes exchanged through this mechanism. Our results also revealed that the extent of genetic exchange correlated with lifestyle and some genomic features. Moreover, the genes in exchanged regions were enriched for functions that reflected specific adaptations, supporting identification of HR as one of the main evolutionary mechanisms shaping prokaryotic core genomes.IMPORTANCE Microbial populations exchange genetic material through a process called homologous recombination. Although this process has been studied in particular organisms, we lack an understanding of its differential impact over the genome and across microbes with different life-styles. We used a common analytical framework to assess this process in a representative set of microorganisms. Our results uncovered important trends. First, microbes with different lifestyles are differentially impacted, with endosymbionts and obligate pathogens being those less prone to undergo this process. Second, certain genetic elements such as restriction-modification systems seem to be associated with higher rates of recombination. Most importantly, recombined genomes show the footprints of natural selection in which recombined regions preferentially contain genes that can be related to specific ecological adaptations. Taken together, our results clarify the relative contributions of factors modulating homologous recombination and show evidence for a clear a role of this process in shaping microbial genomes and driving ecological adaptations.

The Gram-negative bacterium Pseudomonas aeruginosa is a common cause of chronic airway infections in individuals with the heritable disease cystic fibrosis (CF). After prolonged colonization of the CF lung, P. aeruginosa becomes highly resistant to host clearance and antibiotic treatment; therefore, understanding how this bacterium evolves during chronic infection is important for identifying beneficial adaptations that could be targeted therapeutically. To identify potential adaptive traits of P. aeruginosa during chronic infection, we carried out global transcriptomic profiling of chronological clonal isolates obtained from 3 individuals with CF. Isolates were collected sequentially over periods ranging from 3 months to 8 years, representing up to 39,000 in vivo generations. We identified 24 genes that were commonly regulated by all 3 P. aeruginosa lineages, including several genes encoding traits previously shown to be important for in vivo growth. Our results reveal that parallel evolution occurs in the CF lung and that at least a proportion of the traits identified are beneficial for P. aeruginosa chronic colonization of the CF lung.

Beneath the seafloor, microbial life subsists in isolation from the surface world under persistent energy limitation. The nature and extent of genomic evolution in subseafloor microbes have been unknown. Here, we show that the genomes of Thalassospira bacterial populations cultured from million-year-old subseafloor sediments evolve in clonal populations by point mutation, with a relatively low rate of homologous recombination and elevated numbers of pseudogenes. Ratios of nonsynonymous to synonymous substitutions correlate with the accumulation of pseudogenes, consistent with a role for genetic drift in the subseafloor strains but not in type strains of Thalassospira isolated from the surface world. Consistent with this, pangenome analysis reveals that the subseafloor bacterial genomes have a significantly lower number of singleton genes than the type strains, indicating a reduction in recent gene acquisitions. Numerous insertion-deletion events and pseudogenes were present in a flagellar operon of the subseafloor bacteria, indicating that motility is nonessential in these million-year-old subseafloor sediments. This genomic evolution in subseafloor clonal populations coincided with a phenotypic difference: all subseafloor isolates have a lower rate of growth under laboratory conditions than the Thalassospira xiamenensis type strain. Our findings demonstrate that the long-term physical isolation of Thalassospira, in the absence of recombination, has resulted in clonal populations whereby reduced access to novel genetic material from neighbors has resulted in the fixation of new mutations that accumulate in genomes over millions of years. IMPORTANCE The nature and extent of genomic evolution in subseafloor microbial populations subsisting for millions of years below the seafloor are unknown. Subseafloor populations have ultralow metabolic rates that are hypothesized to restrict reproduction and, consequently, the spread of new traits. Our findings demonstrate that genomes of cultivated bacterial strains from the genus Thalassospira isolated from million-year-old abyssal sediment exhibit greatly reduced levels of homologous recombination, elevated numbers of pseudogenes, and genome-wide evidence of relaxed purifying selection. These substitutions and pseudogenes are fixed into the population, suggesting that the genome evolution of these bacteria has been dominated by genetic drift. Thus, reduced recombination, stemming from long-term physical isolation, resulted in small clonal populations of Thalassospira that have accumulated mutations in their genomes over millions of years.

Maturation of dengue viruses (DENVs) alters the structure, immunity, and infectivity of the virion and highly mature particles represent the dominant form in vivo. The production of highly mature virions principally relies on the structure and function of the viral premature membrane protein (prM) and its cleavage by the host protease furin. We redeveloped a reliable clonal cell line (VF1) which produces single-round mature DENVs without the need for DENV reverse genetics. More importantly, using protein engineering and directed evolution of the prM cleavage site, we engineered genetically stable mature DENVs in all serotypes independent of cell or host, usually with minimal impact on viral yield. Using these complementary strategies to regulate maturation, we demonstrate that the resulting mature DENVs are antigenically distinct from their isogenic partially mature forms. Given the clinical importance of mature DENVs in immunity, our study provides reliable strategies and reagents for the production of stable, high-titer mature DENVs for DENV antibody neutralization and vaccination immunity studies. Biologically, our data from directed evolution across host species reveals distinct maturation-dependent selective pressures between mammalian and insect cells, verifying the substrate preference between mammalian and insect furin, while hinting at an evolutionary equilibrium of DENV prM cleavage site between its host and vector in nature. IMPORTANCE Mature DENVs represent the dominant form in vivo and are the target for vaccine development. Here, we used multiple strategies, including protein engineering and natural and directed evolution to generate DENV1, -2, -3, and -4 variants that are highly mature without compromising replication efficiency compared to the parental strains. Given the clinical importance of mature DENVs in immunity, this work provides a roadmap for engineering highly mature DENV that could apply to future vaccine development. Our directed-evolution data also shed light on the divergent evolutionary relationship of DENVs between its host and vector.

Microbial population polymorphisms are commonly observed in natural environments, including long-term infected hosts. However, the underlying processes promoting and stabilizing diversity are difficult to unravel and are not well understood. Here, we use chronic infection of cystic fibrosis airways by the opportunistic pathogen Pseudomonas aeruginosa as a system for investigating bacterial diversification processes during the course of infection. We analyze clonal bacterial isolates sampled during a 32-year period and map temporal and spatial variations in population diversity to different infection sites within the infected host. We show that the ancestral infecting strain diverged into distinct sublineages, each with their own functional and genomic signatures and rates of adaptation, immediately after initial colonization. The sublineages coexisted in the host for decades, suggesting rapid evolution of stable population polymorphisms. Critically, the observed generation and maintenance of population diversity was the result of partitioning of the sublineages into physically separated niches in the CF airway. The results reveal a complex within-host population structure not previously realized and provide evidence that the heterogeneity of the highly structured and complex host environment promotes the evolution and long-term stability of pathogen population diversity during infection.

Community-acquired methicillin-resistant Staphylococcus aureus (CA-MRSA) was recognized in Europe and worldwide in the late 1990s. Within a decade, several genetically and geographically distinct CA-MRSA lineages carrying the small SCCmec type IV and V genetic elements and the Panton-Valentine leukocidin (PVL) emerged around the world. In Europe, the predominant CA-MRSA strain belongs to clonal complex 80 (CC80) and is resistant to kanamycin/amikacin and fusidic acid. CC80 was first reported in 1993 but was relatively rare until the late 1990s. It has since been identified throughout North Africa, the Middle East, and Europe, with recent sporadic reports in sub-Saharan Africa. While strongly associated with skin and soft tissue infections, it is rarely found among asymptomatic carriers. Methicillin-sensitive S. aureus (MSSA) CC80 strains are extremely rare except in sub-Saharan Africa. In the current study, we applied whole-genome sequencing to a global collection of both MSSA and MRSA CC80 isolates. Phylogenetic analyses strongly suggest that the European epidemic CA-MRSA lineage is derived from a PVL-positive MSSA ancestor from sub-Saharan Africa. Moreover, the tree topology suggests a single acquisition of both the SCCmec element and a plasmid encoding the fusidic acid resistance determinant. Four canonical SNPs distinguish the derived CA-MRSA lineage and include a nonsynonymous mutation in accessory gene regulator C (agrC). These changes were associated with a star-like expansion into Europe, the Middle East, and North Africa in the early 1990s, including multiple cases of cross-continent imports likely driven by human migrations.

Cryptococcus gattii is a fungal pathogen of humans, causing pulmonary infections in otherwise healthy hosts. To characterize genomic variation among the four major lineages of C. gattii (VGI, -II, -III, and -IV), we generated, annotated, and compared 16 de novo genome assemblies, including the first for the rarely isolated lineages VGIII and VGIV. By identifying syntenic regions across assemblies, we found 15 structural rearrangements, which were almost exclusive to the VGI-III-IV lineages. Using synteny to inform orthology prediction, we identified a core set of 87% of C. gattii genes present as single copies in all four lineages. Remarkably, 737 genes are variably inherited across lineages and are overrepresented for response to oxidative stress, mitochondrial import, and metal binding and transport. Specifically, VGI has an expanded set of iron-binding genes thought to be important to the virulence of Cryptococcus, while VGII has expansions in the stress-related heat shock proteins relative to the other lineages. We also characterized genes uniquely absent in each lineage, including a copper transporter absent from VGIV, which influences Cryptococcus survival during pulmonary infection and the onset of meningoencephalitis. Through inclusion of population-level data for an additional 37 isolates, we identified a new transcontinental clonal group that we name VGIIx, mitochondrial recombination between VGII and VGIII, and positive selection of multidrug transporters and the iron-sulfur protein aconitase along multiple branches of the phylogenetic tree. Our results suggest that gene expansion or contraction and positive selection have introduced substantial variation with links to mechanisms of pathogenicity across this species complex.

Invasive alien species often have reduced genetic diversity and must adapt to new environments. Given the success of many invasions, this is sometimes called the genetic paradox of invasion. Phytophthora ramorum is invasive, limited to asexual reproduction within four lineages, and presumed clonal. It is responsible for sudden oak death in the United States, sudden larch death in Europe, and ramorum blight in North America and Europe. We sequenced the genomes of 107 isolates to determine how this pathogen can overcome the invasion paradox. Mitotic recombination (MR) associated with transposons and low gene density has generated runs of homozygosity (ROH) affecting 2,698 genes, resulting in novel genotypic diversity within the lineages. One ROH enriched in effectors was fixed in the NA1 lineage. An independent ROH affected the same scaffold in the EU1 lineage, suggesting an MR hot spot and a selection target. Differences in host infection between EU1 isolates with and without the ROH suggest that they may differ in aggressiveness. Non-core regions (not shared by all lineages) had signatures of accelerated evolution and were enriched in putative pathogenicity genes and transposons. There was a striking pattern of gene loss, including all effectors, in the non-core EU2 genome. Positive selection was observed in 8.0% of RxLR and 18.8% of Crinkler effector genes compared with 0.9% of the core eukaryotic gene set. We conclude that the P. ramorum lineages are diverging via a rapidly evolving non-core genome and that the invasive asexual lineages are not clonal, but display genotypic diversity caused by MR.IMPORTANCE Alien species are often successful invaders in new environments, despite the introduction of a few isolates with a reduced genetic pool. This is called the genetic paradox of invasion. We found two mechanisms by which the invasive forest pathogen causing sudden oak and sudden larch death can evolve. Extensive mitotic recombination producing runs of homozygosity generates genotypic diversity even in the absence of sexual reproduction, and rapid turnover of genes in the non-core, or nonessential portion of genome not shared by all isolates, allows pathogenicity genes to evolve rapidly or be eliminated while retaining essential genes. Mitotic recombination events occur in genomic hot spots, resulting in similar ROH patterns in different isolates or groups; one ROH, independently generated in two different groups, was enriched in pathogenicity genes and may be a target for selection. This provides important insights into the evolution of invasive alien pathogens and their potential for adaptation and future persistence.

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.

The clonal population structure of Candida albicans suggests that (para)sexual recombination does not play an important role in the lifestyle of this opportunistic fungal pathogen, an assumption that is strengthened by the fact that most C. albicans strains are heterozygous at the mating type locus (MTL) and therefore mating-incompetent. On the other hand, mating might occur within clonal populations and allow the combination of advantageous traits that were acquired by individual cells to adapt to adverse conditions. We have investigated if parasexual recombination may be involved in the evolution of highly drug-resistant strains exhibiting multiple resistance mechanisms against fluconazole, an antifungal drug that is commonly used to treat infections by C. albicans Growth of strains that were heterozygous for MTL and different fluconazole resistance mutations in the presence of the drug resulted in the emergence of derivatives that had become homozygous for the mutated allele and the mating type locus and exhibited increased drug resistance. When MTLa/a and MTLα/α cells of these strains were mixed in all possible combinations, we could isolate mating products containing the genetic material from both parents. The initial mating products did not exhibit higher drug resistance than their parental strains, but further propagation under selective pressure resulted in the loss of the wild-type alleles and increased fluconazole resistance. Therefore, fluconazole treatment not only selects for resistance mutations but also promotes genomic alterations that confer mating competence, which allows cells in an originally clonal population to exchange individually acquired resistance mechanisms and generate highly drug-resistant progeny.IMPORTANCE Sexual reproduction is an important mechanism in the evolution of species, since it allows the combination of advantageous traits of individual members in a population. The pathogenic yeast Candida albicans is a diploid organism that normally propagates in a clonal fashion, because heterozygosity at the mating type locus (MTL) inhibits mating between cells. Here we show that C. albicans cells that have acquired drug resistance mutations during treatment with the commonly used antifungal agent fluconazole rapidly develop further increased resistance by genome rearrangements that result in simultaneous loss of heterozygosity for the mutated allele and the mating type locus. This enables the drug-resistant cells of a population to switch to the mating-competent opaque morphology and mate with each other to combine different individually acquired resistance mechanisms. The tetraploid mating products reassort their merged genomes and, under selective pressure by the drug, generate highly resistant progeny that have retained the advantageous mutated alleles. Parasexual propagation, promoted by stress-induced genome rearrangements that result in the acquisition of mating competence in cells with adaptive mutations, may therefore be an important mechanism in the evolution of C. albicans populations.

The increasing clinical importance of human infections (frequently severe) caused by Clostridium difficile PCR ribotype 078 (RT078) was first reported in 2008. The severity of symptoms (mortality of ≤30%) and the higher proportion of infections among community and younger patients raised concerns. Farm animals, especially pigs, have been identified as RT078 reservoirs. We aimed to understand the recent changes in RT078 epidemiology by investigating a possible role for antimicrobial selection in its recent evolutionary history. Phylogenetic analysis of international RT078 genomes (isolates from 2006 to 2014, n = 400), using time-scaled, recombination-corrected, maximum likelihood phylogenies, revealed several recent clonal expansions. A common ancestor of each expansion had independently acquired a different allele of the tetracycline resistance gene tetM Consequently, an unusually high proportion (76.5%) of RT078 genomes were tetM positive. Multiple additional tetracycline resistance determinants were also identified (including efflux pump tet40), frequently sharing a high level of nucleotide sequence identity (up to 100%) with sequences found in the pig pathogen Streptococcus suis and in other zoonotic pathogens such as Campylobacter jejuni and Campylobacter coli Each RT078 tetM clonal expansion lacked geographic structure, indicating rapid, recent international spread. Resistance determinants for C. difficile infection-triggering antimicrobials, including fluoroquinolones and clindamycin, were comparatively rare in RT078. Tetracyclines are used intensively in agriculture; this selective pressure, plus rapid, international spread via the food chain, may explain the increased RT078 prevalence in humans. Our work indicates that the use of antimicrobials outside the health care environment has selected for resistant organisms, and in the case of RT078, has contributed to the emergence of a human pathogen.IMPORTANCEClostridium difficile PCR ribotype 078 (RT078) has multiple reservoirs; many are agricultural. Since 2005, this genotype has been increasingly associated with human infections in both clinical settings and the community. Investigations of RT078 whole-genome sequences revealed that tetracycline resistance had been acquired on multiple independent occasions. Phylogenetic analysis revealed a rapid, recent increase in numbers of closely related tetracycline-resistant RT078 (clonal expansions), suggesting that tetracycline selection has strongly influenced its recent evolutionary history. We demonstrate recent international spread of emergent, tetracycline-resistant RT078. A similar tetracycline-positive clonal expansion was also identified in unrelated nontoxigenic C. difficile, suggesting that this process may be widespread and may be independent of disease-causing ability. Resistance to typical C. difficile infection-associated antimicrobials (e.g., fluoroquinolones, clindamycin) occurred only sporadically within RT078. Selective pressure from tetracycline appears to be a key factor in the emergence of this human pathogen and the rapid international dissemination that followed, plausibly via the food chain.

Vibrio parahaemolyticus is the leading cause of seafood-related infections with illnesses undergoing a geographic expansion. In this process of expansion, the most fundamental change has been the transition from infections caused by local strains to the surge of pandemic clonal types. Pandemic clone sequence type 3 (ST3) was the only example of transcontinental spreading until 2012, when ST36 was detected outside the region where it is endemic in the U.S. Pacific Northwest causing infections along the U.S. northeast coast and Spain. Here, we used genome-wide analyses to reconstruct the evolutionary history of the V. parahaemolyticus ST36 clone over the course of its geographic expansion during the previous 25 years. The origin of this lineage was estimated to be in ~1985. By 1995, a new variant emerged in the region and quickly replaced the old clone, which has not been detected since 2000. The new Pacific Northwest (PNW) lineage was responsible for the first cases associated with this clone outside the Pacific Northwest region. After several introductions into the northeast coast, the new PNW clone differentiated into a highly dynamic group that continues to cause illness on the northeast coast of the United States. Surprisingly, the strains detected in Europe in 2012 diverged from this ancestral group around 2000 and have conserved genetic features present only in the old PNW lineage. Recombination was identified as the major driver of diversification, with some preliminary observations suggesting a trend toward a more specialized lifestyle, which may represent a critical element in the expansion of epidemics under scenarios of coastal warming.IMPORTANCEVibrio parahaemolyticus and Vibrio cholerae represent the only two instances of pandemic expansions of human pathogens originating in the marine environment. However, while the current pandemic of V. cholerae emerged more than 50 years ago, the global expansion of V. parahaemolyticus is a recent phenomenon. These modern expansions provide an exceptional opportunity to study the evolutionary process of these pathogens at first hand and gain an understanding of the mechanisms shaping the epidemic dynamics of these diseases, in particular, the emergence, dispersal, and successful introduction in new regions facilitating global spreading of infections. In this study, we used genomic analysis to examine the evolutionary divergence that has occurred over the course of the most recent transcontinental expansion of a pathogenic Vibrio, the spreading of the V. parahaemolyticus sequence type 36 clone from the region where it is endemic on the Pacific coast of North America to the east coast of the United States and finally to the west coast of Europe.

Carbapenem-resistant Klebsiella pneumoniae (CRKP) strains belonging to sequence type 258 (ST258) are frequent causes of hospital-associated outbreaks and are a major contributor to the spread of carbapenemases. This genetic lineage emerged several decades ago and remains a major global health care challenge. In this study, genomic epidemiology was used to investigate the emergence, evolution, and persistence of ST258 carbapenem-resistant K. pneumoniae outbreak-causing lineages at a large tertiary care hospital over 8 years. A time-based phylogenetic analysis of 136 ST258 isolates demonstrated the succession of multiple genetically distinct ST258 sublineages over the 8-year period. Ongoing genomic surveillance identified the emergence and persistence of several distinct clonal ST258 populations. Patterns of multidrug resistance determinants and plasmid replicons were consistent with continued evolution and persistence of these populations. Five ST258 outbreaks were documented, including three that were caused by the same clonal lineage. Mutations in genes encoding effectors of biofilm production and iron acquisition were identified among persistent clones. Two emergent lineages bearing K. pneumoniae integrative conjugative element 10 (ICEKp10) and harboring yersiniabactin and colibactin virulence factors were identified. The results show how distinct ST258 subpopulations have evolved and persisted within the same hospital over nearly a decade.IMPORTANCE The carbapenem class of antibiotics is invaluable for the treatment of selected multidrug-resistant Gram-negative pathogens. The continued transmission of carbapenem-resistant bacteria such as ST258 K. pneumoniae is of serious global public health concern, as treatment options for these infections are limited. This genomic epidemiologic investigation traced the natural history of ST258 K. pneumoniae in a single health care setting over nearly a decade. We found that distinct ST258 subpopulations have caused both device-associated and ward-associated outbreaks, and some of these populations remain endemic within our hospital to the present day. The finding of virulence determinants among emergent ST258 clones supports the idea of convergent evolution of drug-resistant and virulent CRKP strains and highlights the need for continued surveillance, prevention, and control efforts to address emergent and evolving ST258 populations in the health care setting.

Adaptation via natural selection is an important driver of evolution, and repeatable adaptations of replicate populations, under conditions of a constant environment, have been extensively reported. However, isolated groups of populations in nature tend to harbor both genetic and physiological divergence due to multiple selective pressures that they have encountered. How this divergence affects adaptation of these populations to a new common environment remains unclear. To determine the impact of prior genetic and physiological divergence in shaping adaptive evolution to accommodate a new common environment, an experimental evolution study with the sulfate-reducing bacterium Desulfovibrio vulgaris Hildenborough (DvH) was conducted. Two groups of replicate populations with genetic and physiological divergence, derived from a previous evolution study, were propagated in an elevated-temperature environment for 1,000 generations. Ancestor populations without prior experimental evolution were also propagated in the same environment as a control. After 1,000 generations, all the populations had increased growth rates and all but one had greater fitness in the new environment than the ancestor population. Moreover, improvements in growth rate were moderately affected by the divergence in the starting populations, while changes in fitness were not significantly affected. The mutations acquired at the gene level in each group of populations were quite different, indicating that the observed phenotypic changes were achieved by evolutionary responses that differed between the groups. Overall, our work demonstrated that the initial differences in fitness between the starting populations were eliminated by adaptation and that phenotypic convergence was achieved by acquisition of mutations in different genes.IMPORTANCE Improving our understanding of how previous adaptation influences evolution has been a long-standing goal in evolutionary biology. Natural selection tends to drive populations to find similar adaptive solutions for the same selective conditions. However, variations in historical environments can lead to both physiological and genetic divergence that can make evolution unpredictable. Here, we assessed the influence of divergence on the evolution of a model sulfate-reducing bacterium, Desulfovibrio vulgaris Hildenborough, in response to elevated temperature and found a significant effect at the genetic but not the phenotypic level. Understanding how these influences drive evolution will allow us to better predict how bacteria will adapt to various ecological constraints.

Strains of Pseudomonas aeruginosa with deficiencies in DNA mismatch repair have been studied in the context of chronic infection, where elevated mutational rates ("hypermutation") may facilitate the acquisition of antimicrobial resistance. Whether P. aeruginosa hypermutation can also play an adaptive role in the more dynamic context of acute infection remains unclear. In this work, we demonstrate that evolved mismatch repair deficiencies may be exploited by P. aeruginosa to facilitate rapid acquisition of antimicrobial resistance in acute infection, and we directly document rapid clonal succession by such a hypermutating lineage in a patient. Whole-genome sequencing (WGS) was performed on nine serially cultured blood and respiratory isolates from a patient in whom ceftazidime-avibactam (CZA) resistance emerged in vivo over the course of days. The CZA-resistant clone was differentiated by 14 mutations, including a gain-of-function G183D substitution in the PDC-5 chromosomal AmpC cephalosporinase conferring CZA resistance. This lineage also contained a substitution (R656H) at a conserved position in the ATPase domain of the MutS mismatch repair (MMR) protein, and elevated mutational rates were confirmed by mutational accumulation experiments with WGS of evolved lineages in conjunction with rifampin resistance assays. To test whether MMR-deficient hypermutation could facilitate rapid acquisition of CZA resistance, in vitro adaptive evolution experiments were performed with a mutS-deficient strain. These experiments demonstrated rapid hypermutation-facilitated acquisition of CZA resistance compared with the isogenic wild-type strain. Our results suggest a possibly underappreciated role for evolved MMR deficiency in facilitating rapid adaptive evolution of P. aeruginosa in the context of acute infection.IMPORTANCE Antimicrobial resistance in bacteria represents one of the most consequential problems in modern medicine, and its emergence and spread threaten to compromise central advances in the treatment of infectious diseases. Ceftazidime-avibactam (CZA) belongs to a new class of broad-spectrum beta-lactam/beta-lactamase inhibitor combinations designed to treat infections caused by multidrug-resistant bacteria. Understanding the emergence of resistance to this important new drug class is of critical importance. In this work, we demonstrate that evolved mismatch repair deficiency in P. aeruginosa, an important pathogen responsible for significant morbidity and mortality among hospitalized patients, may facilitate rapid acquisition of resistance to CZA in the context of acute infection. These findings are relevant for both diagnosis and treatment of antimicrobial resistance emerging in acute infection in the hypermutator background and additionally have implications for the emergence of more virulent phenotypes.