Comparative Genomics of Serial Candida glabrata Isolates and the Rapid Acquisition of Echinocandin Resistance during Therapy

Candida glabrata hospital isolate from Lebanon reveals micafungin resistance associated with increased chitin and resistance to a cell-surface-disrupting agent

OBJECTIVES

This study aimed to identify the resistance mechanisms to micafungin and fluconazole in a clinical isolate of Candida glabrata.

METHODS

The isolate was whole-genome sequenced to identify amino acid changes in key proteins involved in antifungal resistance, and the isolate was further characterised by pathogenicity-related phenotypic assays that supported the sequencing results.

RESULTS

Amino acid substitutions were detected in 8 of 17 protein candidates. Many of these substitutions were novel, including in CHS3, CHS3B, and KRE5, which are involved in the development of micafungin resistance. Regarding fluconazole resistance, overexpression of efflux pumps was observed. Our isolate did not exhibit an increased virulence potential compared with the control strain; however, a significant increase in chitin content and potential to resist the cell surface disruptant sodium dodecyl sulphate was observed.

CONCLUSIONS

This clinical Candida glabrata isolate experienced a change in cell wall architecture, which correlates with the development of micafungin resistance.

Genomic description of acquired fluconazole- and echinocandin-resistance in patients with serial Candida glabrata isolates

Candida glabrata is one of the most common causes of systemic candidiasis, often resistant to antifungal medications. To describe the genomic context of emerging resistance, we conducted a retrospective analysis of 82 serially collected isolates from 33 patients from population-based candidemia surveillance in the United States. We used whole-genome sequencing to determine the genetic relationships between isolates obtained from the same patient. Phylogenetic analysis demonstrated that isolates from 29 patients were clustered by patient. The median SNPs between isolates from the same patient was 30 (range: 7-96 SNPs), while unrelated strains infected four patients. Twenty-one isolates were resistant to echinocandins, and 24 were resistant to fluconazole. All echinocandin-resistant isolates carried a mutation either in the FKS1 or FKS2 HS1 region. Of the 24 fluconazole-resistant isolates, 17 (71%) had non-synonymous polymorphisms in the PDR1 gene, which were absent in susceptible isolates. In 11 patients, a genetically related resistant isolate was collected after recovering susceptible isolates, indicating in vivo acquisition of resistance. These findings allowed us to estimate the intra-host diversity of C. glabrata and propose an upper boundary of 96 SNPs for defining genetically related isolates, which can be used to assess donor-to-host transmission, nosocomial transmission, or acquired resistance. IMPORTANCE In our study, mutations associated to azole resistance and echinocandin resistance were detected in Candida glabrata isolates using a whole-genome sequence. C. glabrata is the second most common cause of candidemia in the United States, which rapidly acquires resistance to antifungals, in vitro and in vivo.

In vitro activity of ibrexafungerp against clinically relevant echinocandin-resistant Candida strains

Invasive candidiasis is a major hospital-acquired infection. Usually, echinocandins are considered first-line treatment. However, resistant phenotypes have emerged. Ibrexafungerp (IBX) is a new antifungal substance with potent anti-Candida activity. We challenged IBX with a library of 192 pheno-/genotypically echinocandin-resistant Candida isolates, focusing on the substance susceptibility, its activity on certain FKS hotspot (HS) mutated strains, and applying WTULs (wild-type upper limits). Therefore, a 9-year-old strain and patient data collection provided by the German National Reference Center for Invasive Fungal Infections were analyzed. Species identification was confirmed through ITS-sequencing. Molecular susceptibility testing was performed by sequencing HS of the FKS gene. Anidulafungin (AND) and IBX EUCAST-broth-microdilution was conducted. The four most common echinocandin-resistance mediating mutations were found in Candida glabrata [112/192 isolates; F659-(43×) and S663-(48×)] and Candida albicans [63/192 isolates; F641-(15×) and S645-(39×)]. Mutations at the HS-start sequence were associated with higher IBX MIC-values (F659 and F641 (MIC 50/90 mg/L: >4/>4 and 2/4 mg/L) in comparison to AND (F659 and F641 (MIC 50/90: 1/4 and 0.25/1 mg/L). MIC-values in HS-center mutations were almost equal [MIC50/90 in S663: 2/4 (AND and IBX); in S645: 0.5/1 (AND) and 0.25/1 (IBX) mg/L]. In total, 61 vs 78 of 192 echinocandin-resistant isolates may be classified as IBX wild type by applying WTULs, whereas the most prominent effect was seen in C. albicans [48% (30/63) vs 70% (44/63)]. IBX shows in vitro activity against echinocandin-resistant Candida and thus is an addition to the antifungal armory. However, our data suggest that this effect is more pronounced in C. albicans and strains harboring mutations, affecting the HS-center.

Recent gene selection and drug resistance underscore clinical adaptation across Candida species

Understanding how microbial pathogens adapt to treatments, humans and clinical environments is key to infer mechanisms of virulence, transmission and drug resistance. This may help improve therapies and diagnostics for infections with a poor prognosis, such as those caused by fungal pathogens, including Candida. Here we analysed genomic variants across approximately 2,000 isolates from six Candida species (C. glabrata, C. auris, C. albicans, C. tropicalis, C. parapsilosis and C. orthopsilosis) and identified genes under recent selection, suggesting a highly complex clinical adaptation. These involve species-specific and convergently affected adaptive mechanisms, such as adhesion. Using convergence-based genome-wide association studies we identified known drivers of drug resistance alongside potentially novel players. Finally, our analyses reveal an important role of structural variants and suggest an unexpected involvement of (para)sexual recombination in the spread of resistance. Our results provide insights on how opportunistic pathogens adapt to human-related environments and unearth candidate genes that deserve future attention.

Multiple genome analysis of Candida glabrata clinical isolates renders new insights into genetic diversity and drug resistance determinants

The emergence of drug resistance significantly hampers the treatment of human infections, including those caused by fungal pathogens such as Candida species. Candida glabrata ranks as the second most common cause of candidiasis worldwide, supported by rapid acquisition of resistance to azole and echinocandin antifungals frequently prompted by single nucleotide polymorphisms (SNPs) in resistance associated genes, such as PDR1 (azole resistance) or FKS1/2 (echinocandin resistance). To determine the frequency of polymorphisms and genome rearrangements as the possible genetic basis of C. glabrata drug resistance, we assessed genomic variation across 94 globally distributed isolates with distinct resistance phenotypes, whose sequence is deposited in GenBank. The genomes of three additional clinical isolates were sequenced, in this study, including two azole resistant strains that did not display Gain-Of-Function (GOF) mutations in the transcription factor encoding gene PDR1. Genomic variations in susceptible isolates were used to screen out variants arising from genome diversity and to identify variants exclusive to resistant isolates. More than half of the azole or echinocandin resistant isolates do not possess exclusive polymorphisms in PDR1 or FKS1/2, respectively, providing evidence of alternative genetic basis of antifungal resistance. We also identified copy number variations consistently affecting a subset of chromosomes. Overall, our analysis of the genomic and phenotypic variation across isolates allowed to pinpoint, in a genome-wide scale, genetic changes enriched specifically in antifungal resistant strains, which provides a first step to identify additional determinants of antifungal resistance. Specifically, regarding the newly sequenced strains, a set of mutations/genes are proposed to underlie the observed unconventional azole resistance phenotype.

Using genomics to understand the mechanisms of virulence and drug resistance in fungal pathogens

Fungal pathogens pose an increasingly worrying threat to human health, food security and ecosystem diversity. To tackle fungal infections and improve current diagnostic and therapeutic tools it is necessary to understand virulence and antifungal drug resistance mechanisms in diverse species. Recent advances in genomics approaches have provided a suitable framework to understand these phenotypes, which ultimately depend on genetically encoded determinants. In this work, we review how the study of genome sequences has been key to ascertain the bases of virulence and drug resistance traits. We focus on the contribution of comparative genomics, population genomics and directed evolution studies. In addition, we discuss how different types of genomic mutations (small or structural variants) contribute to intraspecific differences in virulence or drug resistance. Finally, we review current challenges in the field and anticipate future directions to solve them. In summary, this work provides a short overview of how genomics can be used to understand virulence and drug resistance in fungal pathogens.

Population genetics and microevolution of clinical Candida glabrata reveals recombinant sequence types and hyper-variation within mitochondrial genomes, virulence genes, and drug targets

Candida glabrata is the second most common etiological cause of worldwide systemic candidiasis in adult patients. Genome analysis of 68 isolates from 8 hospitals across Scotland, together with 83 global isolates, revealed insights into the population genetics and evolution of C. glabrata. Clinical isolates of C. glabrata from across Scotland are highly genetically diverse, including at least 19 separate sequence types that have been recovered previously in globally diverse locations, and 1 newly discovered sequence type. Several sequence types had evidence for ancestral recombination, suggesting transmission between distinct geographical regions has coincided with genetic exchange arising in new clades. Three isolates were missing MATα1, potentially representing a second mating type. Signatures of positive selection were identified in every sequence type including enrichment for epithelial adhesins thought to facilitate fungal adhesin to human epithelial cells. In patent microevolution was identified from 7 sets of recurrent cases of candidiasis, revealing an enrichment for nonsynonymous and frameshift indels in cell surface proteins. Microevolution within patients also affected epithelial adhesins genes, and several genes involved in drug resistance including the ergosterol synthesis gene ERG4 and the echinocandin target FKS1/2, the latter coinciding with a marked drop in fluconazole minimum inhibitory concentration. In addition to nuclear genome diversity, the C. glabrata mitochondrial genome was particularly diverse, with reduced conserved sequence and conserved protein-encoding genes in all nonreference ST15 isolates. Together, this study highlights the genetic diversity within the C. glabrata population that may impact virulence and drug resistance, and 2 major mechanisms generating this diversity: microevolution and genetic exchange/recombination.

Experimental Evolution of Multidrug Resistance in Neurospora crassa under Antifungal Azole Stress

Multidrug resistance, defined as the resistance to multiple drugs in different categories, has been an increasing serious problem. Limited antifungal drugs and the rapid emergence of antifungal resistance prompt a thorough understanding of how the occurrence of multidrug resistance develops and which mechanisms are involved. In this study, experimental evolution was performed under single-azole-drug stress with the model filamentous fungus Neurospora crassa. By about 30 weeks of continuous growth on agar plates containing ketoconazole or voriconazole with weekly transfer, four evolved multidrug-resistant strains 30thK1, 30thK2, 26thV1, and 24thV2 were obtained. Compared to the ancestral strain, all four strains increased resistance not only to commonly used azoles, including ketoconazole, voriconazole, itraconazole, fluconazole, and triadimefon, but also to antifungal drugs in other categories, including terbinafine (allylamine), amorolfine (morpholine), amphotericin B (polyene), polyoxin B (chitin synthesis inhibitor), and carbendazim (β-tubulin inhibitor). After 8 weeks of growth on agar plates without antifungal drugs with weekly transfer, these evolved strains still displayed multidrug-resistant phenotype, suggesting the multidrug resistance could be stably inherited. Transcriptional measurement of drug target genes and drug transporter genes and deletion analysis of the efflux pump gene cdr4 in the evolved strains suggest that overexpression of cdr4 played a major role in the resistance mechanisms for azoles and terbinafine in the evolved strains, particularly for 30thK2 and 26thV1, and evolved drug-resistant strains had less intracellular ketoconazole accumulation and less disruption of ergosterol accumulations under ketoconazole stress compared to wild type. Mutations specifically present in evolved drug-resistant strains were identified by genome re-sequencing, and drug susceptibility test of knockout mutants for most of mutated genes suggests that mutations in 16 genes, functionally novel in drug resistance, potentially contribute to multidrug resistance in evolved strains.

Opportunistic yeast pathogen Candida spp.: Secreted and membrane-bound virulence factors

Candidiasis is a fungal infection caused by Candida spp. especially Candida albicans, C. glabrata, C. parapsilosis and C. tropicalis. Although the medicinal therapeutic strategies have rapidly improved, the mortality rate due to candidiasis has continuously increased. The secreted and membrane-bound virulence factors (VFs) are responsible for fungal invasion, damage and translocation through the host enterocytes besides the evasion from host immune system. VFs such as agglutinin-like sequences (Als), heat shock protein 70, phospholipases, secreted aspartyl proteinases (Sap), lipases, enolases and phytases are mostly hydrolases which degrade the enterocyte membrane components except for candidalysin, the VF acts as a peptide toxin to induce necrotic cell lysis. To date, structural studies of the VFs remain underexplored, hindering their functional analyses. Among the VFs, only secreted aspartyl proteinases and agglutinin-like sequences have their structures deposited in Protein Data Bank (PDB). Therefore, this review scrutinizes the mechanisms of these VFs by discussing the VF-deficient studies of several Candida spp. and their abilities to produce these VFs. Nonetheless, their latest reported sequential and structural analyses are discussed to impart a wider perception of the host-pathogen interactions and potential vaccine or antifungal drug targets. This review signifies that more VFs structural investigations and mining in the emerging Candida spp. are required to decipher their pathogenicity and virulence mechanisms compared to the prominent C. albicans.

LAY ABSTRACT

Candida virulence factors (VFs) including mainly enzymes and proteins play vital roles in breaching the human intestinal barrier and causing deadly candidiasis. Limited VFs' structural studies hinder deeper comprehension of their mechanisms and thus the design of vaccines and antifungal drugs against fungal infections.

The Emergence of Echinocandin-Resistant Candida glabrata Exhibiting High MICs and Related FKS Mutations in Turkey

The frequency of invasive fungal infections shows a rising trend as well as a high morbidity and mortality. Among the causative agents, a shift toward the non-albicans Candida species including Candida glabrata species complex is being observed in several centers. Echinocandin resistance is increasingly published; however, isolates presenting with an in vitro resistance have not yet been reported from Turkey. We, herein, report the first FKS mutant and phenotypically echinocandin-resistant C. glabrata clinical strains from a single center in Turkey. In a 43-year-old female patient, several enterocutaneous fistulae developed after a long term hospitalization period and several complicated surgeries. She eventually required parenteral nutrition via a tunneled central venous catheter (CVC). Following a number of bacteremic and fungemic episodes as well as intensive antimicrobial interventions (including fluconazole, caspofungin and anidulafungin), a CVC-related candidemia caused by C. glabrata was detected. The isolated strain yielded high minimum inhibitory concentration (MIC) values for echinocandins and was categorized as resistant. A resistance-related mutation was detected in FKS2 HS1 (D666V). Blood cultures remained negative after the removal of the CVC and treatment with caspofungin and high-dose fluconazole. Following this first case, two additional C. glabrata strains with high echinocandin MICs were isolated from the urine cultures of two unrelated patients from different wards with different mutations in FKS2 HS1 (S663P and delF659). Our findings indicate that routine antifungal susceptibility testing is crucial and underlines the need for attention for the increasing trend of acquired echinocandin resistance in C. glabrata.

Emergence of resistant Candida glabrata in Germany

Background

Candida glabrata is the second leading fungal pathogen causing candidaemia and invasive candidiasis in Europe. This yeast is recognized for its rapid ability to acquire antifungal drug resistance.

Objectives

We systematically evaluated 176 C. glabrata isolates submitted to the German National Reference Center for Invasive Fungal Infections (NRZMyk) between 2015 and 2019 with regard to echinocandin and fluconazole susceptibility.

Methods

Susceptibility testing was performed using a reference protocol (EUCAST) and a range of commercial assays. Hot spot regions of the echinocandin target FKS genes were sequenced using Sanger sequencing.

Results

In total, 84 of 176 isolates were initially classified as anidulafungin-resistant based on EUCAST testing. Of those, 71 harboured mutations in the glucan synthase encoding FKS genes (13% in FKS1, 87% in FKS2). Significant differences in anidulafungin MICs were found between distinct mutation sites. 11 FKS wild-type (WT) isolates initially classified as resistant exhibited anidulafungin MICs fluctuating around the interpretation breakpoint upon re-testing with multiple assays. Two FKS WT isolates consistently showed high anidulafungin MICs and thus must be considered resistant despite the absence of target gene mutations. Over one-third of echinocandin-resistant strains displayed concomitant fluconazole resistance. Of those, isolates linked to bloodstream infection carrying a change at Ser-663 were associated with adverse clinical outcome.

Conclusions

Resistant C. glabrata strains are emerging in Germany. Phenotypic echinocandin testing can result in misclassification of susceptible strains. FKS genotyping aids in detecting these strains, however, echinocandin resistance may occur despite a wild-type FKS genotype.

Cell wall protein variation, break-induced replication, and subtelomere dynamics in Candida glabrata

Candida glabrata is an opportunistic pathogen of humans, responsible for up to 30% of disseminated candidiasis. Adherence of C. glabrata to host cells is mediated by adhesin-like proteins (ALPs), about half of which are encoded in the subtelomeres. We performed a de novo assembly of two C. glabrata strains, BG2 and BG3993, using long single-molecule real-time (SMRT) reads, and constructed high-quality telomere-to-telomere assemblies of all 13 chromosomes to assess differences between C. glabrata strains. We documented variation between strains, and in agreement with earlier studies, found high (~0.5%-1%) frequencies of SNVs across the genome, including within subtelomeric regions. We documented changes in ALP gene structure and complement: there are large length differences in ALP genes in different strains, resulting from copy number variation in tandem repeats. We compared strains to characterize chromosome rearrangement events including within the poorly characterized subtelomeric regions. We show that rearrangements within the subtelomere regions all affect ALP-encoding genes, and 14/16 involve just the most terminal ALP gene. We present evidence that these rearrangements are mediated by break-induced replication. This study highlights the constrained nature of subtelomeric changes impacting ALP gene complement and subtelomere structure.

DNA damage response of major fungal pathogen Candida glabrata offers clues to explain its genetic diversity

How cells respond to DNA damage is key to maintaining genome integrity or facilitating genetic change. In fungi, DNA damage responses have been extensively characterized in the model budding yeast Saccharomyces cerevisiae, which is generally not pathogenic. However, it is not clear how closely these responses resemble those in fungal pathogens, in which genetic change plays an important role in the evolutionary arms race between pathogen and host and the evolution of antifungal drug resistance. A close relative of S. cerevisiae, Candida glabrata, is an opportunistic pathogen that displays high variability in chromosome structure among clinical isolates and rapidly evolves antifungal drug resistance. The mechanisms facilitating such genomic flexibility and evolvability in this organism are unknown. Recently we characterized the DNA damage response of C. glabrata and identified several features that distinguish it from the well characterized DNA damage response of S. cerevisiae. First, we discovered that, in contrast to the established paradigm, C. glabrata effector kinase Rad53 is not hyperphosphorylated upon DNA damage. We also uncovered evidence of an attenuated DNA damage checkpoint response, wherein in the presence of DNA damage C. glabrata cells did not accumulate in S-phase and proceeded with cell division, leading to aberrant mitoses and cell death. Finally, we identified evidence of transcriptional rewiring of the DNA damage response of C. glabrata relative to S. cerevisiae, including an upregulation of genes involved in mating and meiosis-processes that have not been reported in C. glabrata. Together, these results open new possibilities and raise tantalizing questions of how this major fungal pathogen facilitates genetic change.

A Noncanonical DNA Damage Checkpoint Response in a Major Fungal Pathogen

In order to preserve genome integrity, all cells must mount appropriate responses to DNA damage, including slowing down or arresting the cell cycle to give the cells time to repair the damage and changing gene expression, for example to induce genes involved in DNA repair. The Rad53 protein kinase is a conserved central mediator of these responses in eukaryotic cells, and its extensive phosphorylation upon DNA damage is necessary for its activation and subsequent activity.

DNA damage checkpoints are key guardians of genome integrity. Eukaryotic cells respond to DNA damage by triggering extensive phosphorylation of Rad53/CHK2 effector kinase, whereupon activated Rad53/CHK2 mediates further aspects of checkpoint activation, including cell cycle arrest and transcriptional changes. Budding yeast Candida glabrata, closely related to model eukaryote Saccharomyces cerevisiae, is an opportunistic pathogen characterized by high genetic diversity and rapid emergence of drug-resistant mutants. However, the mechanisms underlying this genetic variability are unclear. We used Western blotting and mass spectrometry to show that, unlike S. cerevisiae, C. glabrata cells exposed to DNA damage did not induce C. glabrata Rad53 (CgRad53) phosphorylation. Furthermore, flow cytometry analysis showed that, unlike S. cerevisiae, C. glabrata cells did not accumulate in S phase upon DNA damage. Consistent with these observations, time-lapse microscopy showed C. glabrata cells continuing to divide in the presence of DNA damage, resulting in mitotic errors and cell death. Finally, transcriptome sequencing (RNAseq) analysis revealed transcriptional rewiring of the DNA damage response in C. glabrata and identified several key protectors of genome stability upregulated by DNA damage in S. cerevisiae but downregulated in C. glabrata, including proliferating cell nuclear antigen (PCNA). Together, our results reveal a noncanonical fungal DNA damage response in C. glabrata, which may contribute to rapidly generating genetic change and drug resistance.

Antifungal Susceptibility of Clinical Yeast Isolates from a Large Canadian Reference Laboratory and Application of Whole-Genome Sequence Analysis To Elucidate Mechanisms of Acquired Resistance

To understand the epidemiology and susceptibility patterns of yeast infections in Ontario, Canada, we examined 4,715 clinical yeast isolates submitted to our laboratory for antifungal susceptibility testing from 2014 to 2018. Candida albicans was the most frequently submitted species (43.0%), followed by C. glabrata (21.1%), C. parapsilosis (15.0%), and C. tropicalis (6.2%). Twenty-three other Candida spp. (11.6%) and 4 non-Candida species (3.1%) were also identified. Few changes in species distribution were observed from 2014 to 2018, but the total numbers of yeast isolates sent for testing increased, with an annual 7.4% change. According to CLSI clinical breakpoints, resistance rates remained low overall. Moderate fluconazole resistance was noted among C. glabrata (9%), C. parapsilosis (9%), and C. tropicalis (12%) isolates. Only 1% of C. glabrata isolates were resistant to caspofungin, micafungin, and anidulafungin. Whole-genome sequence analysis confirmed 11 cases of acquired resistance to azoles or echinocandins via in-host evolution. There were mutations in the gene for the catalytic subunit of 1,3-beta-glucan synthase-mediated echinocandin resistance in 3 of 3 C. albicans strains, 3 of 4 C. glabrata strains, and 1 strain of C. tropicalis Azole resistance was likely caused by a homozygous ERG3 mutation in 1 C. albicans strain and a previously undescribed chromosomal-duplication event involving ERG11 and TAC1 orthologs in 1 C. tropicalis strain. While antifungal resistance rates remain low among yeast isolates in Ontario, ongoing surveillance is necessary to inform empirical therapy for optimal patient management and to guide antifungal stewardship.

The Quiet and Underappreciated Rise of Drug-Resistant Invasive Fungal Pathogens

Human fungal pathogens are attributable to a significant economic burden and mortality worldwide. Antifungal treatments, although limited in number, play a pivotal role in decreasing mortality and morbidities posed by invasive fungal infections (IFIs). However, the recent emergence of multidrug-resistant Candida auris and Candida glabrata and acquiring invasive infections due to azole-resistant C. parapsilosis, C. tropicalis, and Aspergillus spp. in azole-naïve patients pose a serious health threat considering the limited number of systemic antifungals available to treat IFIs. Although advancing for major fungal pathogens, the understanding of fungal attributes contributing to antifungal resistance is just emerging for several clinically important MDR fungal pathogens. Further complicating the matter are the distinct differences in antifungal resistance mechanisms among various fungal species in which one or more mechanisms may contribute to the resistance phenotype. In this review, we attempt to summarize the burden of antifungal resistance for selected non-albicansCandida and clinically important Aspergillus species together with their phylogenetic placement on the tree of life. Moreover, we highlight the different molecular mechanisms between antifungal tolerance and resistance, and comprehensively discuss the molecular mechanisms of antifungal resistance in a species level.

Genome Comparisons of Candida glabrata Serial Clinical Isolates Reveal Patterns of Genetic Variation in Infecting Clonal Populations

Candida glabrata is an opportunistic fungal pathogen that currently ranks as the second most common cause of candidiasis. Although the mechanisms underlying virulence and drug resistance in C. glabrata are now starting to be elucidated, we still lack a good understanding of how this yeast adapts during the course of an infection. Outstanding questions are whether the observed genomic plasticity of C. glabrata plays a role during infection, or what levels of genetic variation exist within an infecting clonal population. To shed light onto the genomic variation within infecting C. glabrata populations, we compared the genomes of 11 pairs and one trio of serial clinical isolates, each obtained from a single patient. Our results provide a catalog of genetic variations existing within clonal infecting isolates, and reveal an enrichment of non-synonymous changes in genes encoding cell-wall proteins. Genetic variation and the presence of non-synonymous mutations and copy number variations accumulated within the host, suggest that clonal populations entail a non-negligible level of genetic variation that may reflect selection processes that occur within the human body. As we show here, these genomic changes can underlie phenotypic differences in traits that are relevant for infection.