Permalink
Cannot retrieve contributors at this time
Fetching contributors…
Cannot retrieve contributors at this time
| <?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.1d1 20130915//EN" "JATS-archivearticle1.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article" dtd-version="1.1d1"><front><journal-meta><journal-id journal-id-type="nlm-ta">elife</journal-id><journal-id journal-id-type="hwp">eLife</journal-id><journal-id journal-id-type="publisher-id">eLife</journal-id><journal-title-group><journal-title>eLife</journal-title></journal-title-group><issn publication-format="electronic">2050-084X</issn><publisher><publisher-name>eLife Sciences Publications, Ltd</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="publisher-id">00662</article-id><article-id pub-id-type="doi">10.7554/eLife.00662</article-id><article-categories><subj-group subj-group-type="display-channel"><subject>Research article</subject></subj-group><subj-group subj-group-type="heading"><subject>Genomics and evolutionary biology</subject></subj-group><subj-group subj-group-type="heading"><subject>Microbiology and infectious disease</subject></subj-group></article-categories><title-group><article-title>The evolution of drug resistance in clinical isolates of <italic>Candida albicans</italic></article-title></title-group><contrib-group><contrib contrib-type="author" id="author-20807" equal-contrib="yes"><name><surname>Ford</surname><given-names>Christopher B</given-names></name><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="aff" rid="aff2">2</xref><xref ref-type="fn" rid="equal-contrib1">†</xref><xref ref-type="other" rid="par-3"/><xref ref-type="fn" rid="con1"/><xref ref-type="fn" rid="conf2"/><xref ref-type="other" rid="dataro1"/></contrib><contrib contrib-type="author" id="author-4626" equal-contrib="yes"><name><surname>Funt</surname><given-names>Jason M</given-names></name><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="aff" rid="aff2">2</xref><xref ref-type="fn" rid="equal-contrib1">†</xref><xref ref-type="other" rid="par-1"/><xref ref-type="fn" rid="con2"/><xref ref-type="fn" rid="conf2"/></contrib><contrib contrib-type="author" id="author-4627" equal-contrib="yes"><name><surname>Abbey</surname><given-names>Darren</given-names></name><xref ref-type="aff" rid="aff4">4</xref><xref ref-type="fn" rid="equal-contrib2">‡</xref><xref ref-type="fn" rid="con4"/><xref ref-type="fn" rid="conf2"/></contrib><contrib contrib-type="author" id="author-4628" equal-contrib="yes"><name><surname>Issi</surname><given-names>Luca</given-names></name><xref ref-type="aff" rid="aff5">5</xref><xref ref-type="fn" rid="equal-contrib2">‡</xref><xref ref-type="fn" rid="con5"/><xref ref-type="fn" rid="conf2"/></contrib><contrib contrib-type="author" id="author-20802"><name><surname>Guiducci</surname><given-names>Candace</given-names></name><xref ref-type="aff" rid="aff2">2</xref><xref ref-type="fn" rid="con6"/><xref ref-type="fn" rid="conf2"/></contrib><contrib contrib-type="author" id="author-20803"><name><surname>Martinez</surname><given-names>Diego A</given-names></name><xref ref-type="aff" rid="aff2">2</xref><xref ref-type="fn" rid="con8"/><xref ref-type="fn" rid="conf2"/></contrib><contrib contrib-type="author" id="author-20805"><name><surname>Delorey</surname><given-names>Toni</given-names></name><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0001-6614-3803</contrib-id><xref ref-type="aff" rid="aff2">2</xref><xref ref-type="fn" rid="con9"/><xref ref-type="fn" rid="conf2"/></contrib><contrib contrib-type="author" id="author-20804"><name><surname>Li</surname><given-names>Bi yu</given-names></name><xref ref-type="aff" rid="aff2">2</xref><xref ref-type="fn" rid="con10"/><xref ref-type="fn" rid="conf2"/></contrib><contrib contrib-type="author" id="author-4630"><name><surname>White</surname><given-names>Theodore C</given-names></name><xref ref-type="aff" rid="aff6">6</xref><xref ref-type="fn" rid="con11"/><xref ref-type="fn" rid="conf2"/></contrib><contrib contrib-type="author" id="author-20806"><name><surname>Cuomo</surname><given-names>Christina</given-names></name><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0002-5778-960X</contrib-id><xref ref-type="aff" rid="aff2">2</xref><xref ref-type="fn" rid="con7"/><xref ref-type="fn" rid="conf2"/><xref ref-type="other" rid="dataro1"/></contrib><contrib contrib-type="author" id="author-4631"><name><surname>Rao</surname><given-names>Reeta P</given-names></name><xref ref-type="aff" rid="aff5">5</xref><xref ref-type="fn" rid="con12"/><xref ref-type="fn" rid="conf2"/></contrib><contrib contrib-type="author" id="author-8986"><name><surname>Berman</surname><given-names>Judith</given-names></name><xref ref-type="aff" rid="aff4">4</xref><xref ref-type="aff" rid="aff7">7</xref><xref ref-type="fn" rid="con13"/><xref ref-type="fn" rid="conf2"/></contrib><contrib contrib-type="author" corresp="yes" id="author-1388"><name><surname>Thompson</surname><given-names>Dawn A</given-names></name><xref ref-type="aff" rid="aff2">2</xref><xref ref-type="corresp" rid="cor1">*</xref><xref ref-type="fn" rid="con3"/><xref ref-type="fn" rid="conf2"/></contrib><contrib contrib-type="author" corresp="yes" id="author-1171"><name><surname>Regev</surname><given-names>Aviv</given-names></name><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="aff" rid="aff2">2</xref><xref ref-type="aff" rid="aff3">3</xref><xref ref-type="corresp" rid="cor2">*</xref><xref ref-type="other" rid="par-2"/><xref ref-type="other" rid="par-4"/><xref ref-type="other" rid="par-5"/><xref ref-type="other" rid="par-6"/><xref ref-type="fn" rid="con14"/><xref ref-type="fn" rid="conf1"/></contrib><aff id="aff1"><label>1</label><institution content-type="dept">Department of Biology</institution>, <institution>Broad Institute of MIT and Harvard</institution>, <addr-line><named-content content-type="city">Cambridge</named-content></addr-line>, <country>United States</country></aff><aff id="aff2"><label>2</label><institution>Broad Institute of MIT and Harvard</institution>, <addr-line><named-content content-type="city">Cambridge</named-content></addr-line>, <country>United States</country></aff><aff id="aff3"><label>3</label><institution content-type="dept">Department of Biology</institution>, <institution>Howard Hughes Medical Institute, Massachusetts Institute of Technology</institution>, <addr-line><named-content content-type="city">Cambridge</named-content></addr-line>, <country>United States</country></aff><aff id="aff4"><label>4</label><institution content-type="dept">Department of Genetics, Cell Biology and Development</institution>, <institution>University of Minnesota</institution>, <addr-line><named-content content-type="city">Minneapolis</named-content></addr-line>, <country>United States</country></aff><aff id="aff5"><label>5</label><institution content-type="dept">Department of Biology and Biotechnology</institution>, <institution>Worcester Polytechnic Institute</institution>, <addr-line><named-content content-type="city">Worcester</named-content></addr-line>, <country>United States</country></aff><aff id="aff6"><label>6</label><institution content-type="dept">School of Biological Sciences</institution>, <institution>University of Missouri at Kansas City</institution>, <addr-line><named-content content-type="city">Kansas City</named-content></addr-line>, <country>United States</country></aff><aff id="aff7"><label>7</label><institution content-type="dept">Department of Molecular Microbiology and Biotechnology</institution>, <institution>Tel Aviv University</institution>, <addr-line><named-content content-type="city">Tel Aviv</named-content></addr-line>, <country>Israel</country></aff></contrib-group><contrib-group content-type="section"><contrib contrib-type="editor"><name><surname>Dermitzakis</surname><given-names>Emmanouil T</given-names></name><role>Reviewing editor</role><aff><institution>University of Geneva Medical School</institution>, <country>Switzerland</country></aff></contrib></contrib-group><author-notes><corresp id="cor1"><label>*</label>For correspondence: <email>dawnt@broadinstitute.org</email> (DAT);</corresp><corresp id="cor2"><email>aregev@broadinstitute.org</email> (AR)</corresp><fn fn-type="con" id="equal-contrib1"><label>†</label><p>These authors contributed equally to this work</p></fn><fn fn-type="con" id="equal-contrib2"><label>‡</label><p>These authors contributed equally as second authors</p></fn></author-notes><pub-date publication-format="electronic" date-type="pub"><day>03</day><month>02</month><year>2015</year></pub-date><pub-date pub-type="collection"><year>2015</year></pub-date><volume>4</volume><elocation-id>e00662</elocation-id><history><date date-type="received"><day>22</day><month>02</month><year>2013</year></date><date date-type="accepted"><day>18</day><month>12</month><year>2014</year></date></history><permissions><copyright-statement>© 2014, Ford et al</copyright-statement><copyright-year>2014</copyright-year><copyright-holder>Ford et al</copyright-holder><license xlink:href="http://creativecommons.org/licenses/by/4.0/"><license-p>This article is distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License</ext-link>, which permits unrestricted use and redistribution provided that the original author and source are credited.</license-p></license></permissions><self-uri content-type="pdf" xlink:href="elife00662.pdf"/><abstract><object-id pub-id-type="doi">10.7554/eLife.00662.001</object-id><p><italic>Candida albicans</italic> is both a member of the healthy human microbiome and a major pathogen in immunocompromised individuals. Infections are typically treated with azole inhibitors of ergosterol biosynthesis often leading to drug resistance. Studies in clinical isolates have implicated multiple mechanisms in resistance, but have focused on large-scale aberrations or candidate genes, and do not comprehensively chart the genetic basis of adaptation. Here, we leveraged next-generation sequencing to analyze 43 isolates from 11 oral candidiasis patients. We detected newly selected mutations, including single-nucleotide polymorphisms (SNPs), copy-number variations and loss-of-heterozygosity (LOH) events. LOH events were commonly associated with acquired resistance, and SNPs in 240 genes may be related to host adaptation. Conversely, most aneuploidies were transient and did not correlate with drug resistance. Our analysis also shows that isolates also varied in adherence, filamentation, and virulence. Our work reveals new molecular mechanisms underlying the evolution of drug resistance and host adaptation.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.001">http://dx.doi.org/10.7554/eLife.00662.001</ext-link></p></abstract><abstract abstract-type="executive-summary"><object-id pub-id-type="doi">10.7554/eLife.00662.002</object-id><title>eLife digest</title><p>Nearly all humans are infected with the fungus <italic>Candida albicans</italic>. In most people, the infection does not produce any symptoms because their immune system is able to counteract the fungus' attempts to spread around the body. However, if the balance between fungal attack and body defence fails, the fungus is able to spread, which can lead to serious disease that is fatal in 42% of cases.</p><p>How does <italic>C. albicans</italic> outcompete the body's defences to cause disease? This is a pertinent question because the most effective antifungal medicines—including the drug fluconazole—do not kill the fungus; they only stop it from growing. This gives the fungus time to develop resistance to the drug by becoming able to quickly replace the fungal proteins the drug destroys, or to efficiently remove the drug from its cells.</p><p>In this study, Ford et al. studied the changes that occur in the DNA of <italic>C. albicans</italic> over time in patients who are being treated with fluconazole. Ford et al. took 43 samples of <italic>C. albicans</italic> from 11 patients with weakened immune systems. The experiments show that the fungus samples collected early on were more sensitive to the drug than the samples collected later.</p><p>In most cases, the genetic data suggest that the infections begin with a single fungal cell; the cells in the later samples are its offspring. Despite this, there is a lot of genetic variation between samples from the same patient, which indicates that the fungus is under pressure to become more resistant to the drug. There were 240 genes—including those that can alter the surface on the fungus cells to make it better at evading the host immune system—in which small changes occurred over time in three or more patients. Laboratory tests revealed that many of these genes are likely important for the fungus to survive in an animal host in the presence of the drug.</p><p><italic>C. albicans</italic> cells usually have two genetically distinct copies of every gene. Ford et al. found that for some genes—including some that make surface components or are involved in expelling drugs from cells—the loss of genetic information from one copy, so that both copies become identical, is linked to resistance to fluconazole. However, the gain of whole or partial chromosomes—which contain large numbers of genes—is not linked to resistance, but may provide additional genetic material for generating diversity in the yeast population that may help the cells to evolve resistance in the future.</p><p>These experiments have identified many new candidate genes that are important for drug resistance and evading the host immune system, and which could be used to guide the development of new therapeutics to treat these life-threatening infections.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.002">http://dx.doi.org/10.7554/eLife.00662.002</ext-link></p></abstract><kwd-group kwd-group-type="author-keywords"><title>Author keywords</title><kwd>Candida albicans</kwd><kwd>drug resistance</kwd><kwd>evolution</kwd><kwd>genomics</kwd><kwd>virulence</kwd></kwd-group><kwd-group kwd-group-type="research-organism"><title>Research organism</title><kwd>other</kwd></kwd-group><funding-group><award-group id="par-1"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000001</institution-id><institution>National Science Foundation</institution></institution-wrap></funding-source><award-id>Graduate Student Fellowship</award-id><principal-award-recipient><name><surname>Funt</surname><given-names>Jason M</given-names></name></principal-award-recipient></award-group><award-group id="par-2"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000011</institution-id><institution content-type="university">Howard Hughes Medical Institute</institution></institution-wrap></funding-source><award-id>Full investigator</award-id><principal-award-recipient><name><surname>Regev</surname><given-names>Aviv</given-names></name></principal-award-recipient></award-group><award-group id="par-3"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100005237</institution-id><institution>Helen Hay Whitney Foundation</institution></institution-wrap></funding-source><award-id>Postdoctoral Fellowship</award-id><principal-award-recipient><name><surname>Ford</surname><given-names>Christopher B</given-names></name></principal-award-recipient></award-group><award-group id="par-4"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000002</institution-id><institution content-type="university">National Institutes of Health</institution></institution-wrap></funding-source><award-id>8DP1CA174427</award-id><principal-award-recipient><name><surname>Regev</surname><given-names>Aviv</given-names></name></principal-award-recipient></award-group><award-group id="par-5"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000002</institution-id><institution content-type="university">National Institutes of Health</institution></institution-wrap></funding-source><award-id>2R01CA119176-01</award-id><principal-award-recipient><name><surname>Regev</surname><given-names>Aviv</given-names></name></principal-award-recipient></award-group><award-group id="par-6"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000879</institution-id><institution>Alfred P. Sloan Foundation</institution></institution-wrap></funding-source><principal-award-recipient><name><surname>Regev</surname><given-names>Aviv</given-names></name></principal-award-recipient></award-group><funding-statement>The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.</funding-statement></funding-group><custom-meta-group><custom-meta><meta-name>elife-xml-version</meta-name><meta-value>2.0</meta-value></custom-meta><custom-meta specific-use="meta-only"><meta-name>Author impact statement</meta-name><meta-value>Loss-of-heterozygosity mutations, but not aneuploidies, are linked to the evolution of drug resistance in <italic>Candida albicans</italic> isolated from immunocompromised patients.</meta-value></custom-meta></custom-meta-group></article-meta></front><body><sec sec-type="intro" id="s1"><title>Introduction</title><p>Virtually all humans are colonized with <italic>Candida albicans,</italic> but in some individuals this benign commensal organism becomes a serious, life-threatening pathogen. <italic>C. albicans</italic> possesses an arsenal of traits that promote its pathogenicity, including phenotypic switching (<xref ref-type="bibr" rid="bib1">Alby and Bennett, 2009</xref>), yeast–hyphae transition (<xref ref-type="bibr" rid="bib43">Kumamoto and Vinces, 2005</xref>) and the secretion of molecules that promote adhesion to abiotic surfaces (<xref ref-type="bibr" rid="bib9">Chandra et al., 2001</xref>). As a commensal, an intricate balance is maintained between the ability of <italic>C. albicans</italic> to invade host tissues and the host's defense mechanisms (<xref ref-type="bibr" rid="bib41">Kim and Sudbery, 2011</xref>; <xref ref-type="bibr" rid="bib42">Kumamoto and Pierce, 2011</xref>). Alteration of this delicate host–fungus balance can result in high levels of patient mortality (<xref ref-type="bibr" rid="bib60">Pittet et al., 1994</xref>; <xref ref-type="bibr" rid="bib11">Charles et al., 2003</xref>): systemic <italic>C. albicans</italic> infections are fatal in 42% of cases (<xref ref-type="bibr" rid="bib88">Wisplinghoff et al., 2003</xref>), despite the use of antifungal therapies, and <italic>C. albicans</italic> is the fourth most common infection in hospitals (<xref ref-type="bibr" rid="bib32">Gudlaugsson et al., 2003</xref>; <xref ref-type="bibr" rid="bib56">Pappas et al., 2003</xref>). While compromised immune function contributes to pathogenesis (<xref ref-type="bibr" rid="bib28">Gow and Hube, 2012</xref>), it is less clear how <italic>C. albicans</italic> evolves to better exploit the host environment during the course of infection.</p><p>Two classes of antifungals in clinical use target ergosterol, a major component of the fungal cell membrane: polyenes and azoles. Polyenes (e.g., Amphotericin B) are used sparingly due to toxicity (<xref ref-type="bibr" rid="bib64">Rex et al., 1994</xref>), whereas azoles (e.g., fluconazole) are used widely because they can be administered orally and have few side effects (<xref ref-type="bibr" rid="bib65">Rex et al., 2003</xref>). However, resistance to the azoles arises within the commensal population of the treated individual, primarily because azoles are fungistatic (inhibit growth but do not kill) (<xref ref-type="bibr" rid="bib14">Cowen et al., 2002</xref>). Epidemiological data suggest that the intensity of fluconazole use is driving the appearance of resistant isolates (<xref ref-type="bibr" rid="bib59">Pfaller et al., 1998</xref>). Studies of clinical isolates of <italic>C. albicans</italic> suggest that drug resistance can increase during an infection through the acquisition of aneuploidies (<xref ref-type="bibr" rid="bib75">Selmecki et al., 2009</xref>) due to genomic plasticity and rapid evolutionary selection during infection.</p><p>Previous studies have identified two molecular mechanisms of azole resistance in <italic>C. albicans</italic>. First, increased activity or level of the enzymes of the ergosterol pathway (e.g., <italic>ERG11</italic>) reduces direct impact of the drug on its target (<xref ref-type="bibr" rid="bib3">Asai et al., 1999</xref>; <xref ref-type="bibr" rid="bib55">Oliver et al., 2007</xref>). Second, increased efflux of the drug from cells by ABC transporters (encoded by <italic>CDR1</italic> and <italic>CDR2</italic>) (<xref ref-type="bibr" rid="bib13">Coste et al., 2006</xref>) or by the major facilitator superfamily efflux pump (encoded by <italic>MDR1</italic>) (<xref ref-type="bibr" rid="bib18">Dunkel et al., 2008</xref>) reduces the effective intracellular drug concentration. In both cases, such alterations can result from point mutations in genes encoding these proteins (<xref ref-type="bibr" rid="bib50">Marichal et al., 1999</xref>), in transcription factors that regulate mRNA expression levels (<xref ref-type="bibr" rid="bib48">MacPherson et al., 2005</xref>; <xref ref-type="bibr" rid="bib13">Coste et al., 2006</xref>; <xref ref-type="bibr" rid="bib18">Dunkel et al., 2008</xref>), or from increased copy number of the relevant genes, via genome rearrangements such as whole chromosome and segmental aneuploidies (<xref ref-type="bibr" rid="bib72">Selmecki et al., 2006</xref>; <xref ref-type="bibr" rid="bib74">2008</xref>; <xref ref-type="bibr" rid="bib75">2009</xref>). Indeed, the genomes of drug-resistant strains isolated following clinical treatment often exhibit large-scale changes, such as loss of heterozygosity (LOH) (<xref ref-type="bibr" rid="bib13">Coste et al., 2006</xref>; <xref ref-type="bibr" rid="bib19">Dunkel and Morschhauser, 2011</xref>), copy-number variation (CNV), including short segmental CNV, and whole chromosome aneuploidy (<xref ref-type="bibr" rid="bib73">Selmecki et al., 2010</xref>) accompanied by point mutations.</p><p>While we understand some aspects of the molecular basis of resistance, we understand less about the mechanisms that drive the evolution of drug resistance and overall pathogenicity in <italic>C. albicans</italic>. It is challenging to use forward genetic approaches in <italic>C. albicans</italic> due to its diploid genome and lack of a complete sexual cycle. Although <italic>C. albicans</italic> has conserved the genomic elements needed for mating, mating occurs instead through rare mating-competent haploids (<xref ref-type="bibr" rid="bib33">Hickman et al., 2013</xref>) or via a parasexual cycle consisting of mating of diploid strains to form tetraploids followed by chromosome loss to regenerate diploids (<xref ref-type="bibr" rid="bib4">Bennett and Johnson, 2005</xref>). An alternative approach is to use isolates sampled consecutively from the same patient to study the changes in the frequency of variants in natural populations under selection for drug resistance. Studies in evolved isolates have implicated multiple mechanisms in drug resistance, but have focused on large-scale aberrations such as aneuploidies and LOH (<xref ref-type="bibr" rid="bib74">Selmecki et al., 2008</xref>; <xref ref-type="bibr" rid="bib75">2009</xref>) or candidate genes (<xref ref-type="bibr" rid="bib58">Perea et al., 2001</xref>; <xref ref-type="bibr" rid="bib87">White et al., 2002</xref>), and do not comprehensively chart the genetic basis of adaptation.</p><p>Here, we used genome sequencing of isolates sampled consecutively from patients that were clinically treated with fluconazole to systematically analyze the genetic dynamics that accompany the appearance of drug resistance during oral candidiasis in human HIV patients. Most isolates from each individual patient were highly related, suggesting a clonal population structure and facilitating the identification of variation. Because each clinical sample was purified from a single colony, we cannot assess the population structure at any single time point. Instead, we have measured the occurrence of single-nucleotide polymorphisms (SNPs), CNV, and LOH events in each isolate and then compared them between isolates from the same patient and across patients' series. Consistent with previous studies, we found that LOH events were recurrent across patients' series and were associated with increased drug resistance. To identify SNPs with likely functional impact in the context of substantial genetic diversity, we focused on those events that were both persistent across isolates within a patient and were recurrent in the same gene across multiple patient series. We found 240 genes that recurrently contain persistent SNPs, many of which may be related not only to antifungal exposure but also to the complex process of adaptation to the host and antifungal exposure. In contrast, aneuploidies were prevalent in the isolates, yet they were more likely to be transient, and aneuploidy, per se, did not correlate with changes in drug resistance. Our work uses comparative analysis of a fungal pathogen to reveal new molecular mechanisms underlying drug resistance and host adaptation and provides a general model for such studies in other eukaryotic pathogens.</p></sec><sec sec-type="results" id="s2"><title>Results</title><sec id="s2-1"><title>Whole genome sequencing of 43 serial clinical isolates from 11 patients</title><p>To study the in vivo evolution of azole resistance in <italic>C. albicans</italic>, we analyzed 43 longitudinal isolates from 11 HIV-infected patients with oropharyngeal candidiasis (<xref ref-type="bibr" rid="bib85">White, 1997a</xref>; <xref ref-type="bibr" rid="bib58">Perea et al., 2001</xref>) (<xref ref-type="table" rid="tbl1">Table 1</xref>). The isolates were previously collected during incidences of infection and form a time series from each patient (2–16 isolates per series; <xref ref-type="fig" rid="fig1">Figure 1</xref>, <xref ref-type="fig" rid="fig2">Figure 2A</xref>). Each isolate was derived from a single colony, and thus, represents a single diploid genotype sampled from the within-host <italic>C. albicans</italic> population at the respective time point. In each series, the first isolate (‘progenitor’) was collected prior to any treatment with azole antifungals and the remaining isolates were collected at later, typically consecutive, time points, culminating in the final ‘endpoint’ isolate (<xref ref-type="table" rid="tbl1">Table 1</xref>).<table-wrap id="tbl1" position="float"><object-id pub-id-type="doi">10.7554/eLife.00662.003</object-id><label>Table 1.</label><caption><p>Isolate history and sequencing summary</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.003">http://dx.doi.org/10.7554/eLife.00662.003</ext-link></p></caption><table frame="hsides" rules="groups"><thead><tr><th>Publication name</th><th>PT</th><th>Strain</th><th>Entry date</th><th>Drug treatment</th><th>Dose (mg/day)</th><th>E-test MIC (ug/mL)</th><th>Depth of coverage</th><th>Reads</th><th>Percent aligned</th></tr></thead><tbody><tr><td>White, T.C.</td><td>1</td><td>1</td><td>9/10/90</td><td>Fluconazole</td><td>100</td><td>0.25</td><td>111.96</td><td>9,896,468</td><td>87.17%</td></tr><tr><td/><td/><td>2</td><td>12/14/90</td><td>Fluconazole</td><td>100</td><td>1</td><td>69.20</td><td>12,797,328</td><td>87.43%</td></tr><tr><td/><td/><td>3</td><td>12/21/90</td><td>Fluconazole</td><td>100</td><td>4</td><td>92.04</td><td>16,987,814</td><td>86.87%</td></tr><tr><td/><td/><td>4</td><td>12/31/90</td><td>Fluconazole</td><td>100</td><td>3</td><td>80.69</td><td>14,858,710</td><td>87.81%</td></tr><tr><td/><td/><td>5</td><td>2/8/91</td><td>Fluconazole</td><td>100</td><td>4</td><td>110.80</td><td>20,484,584</td><td>86.75%</td></tr><tr><td/><td/><td>6</td><td>2/22/91</td><td>Fluconazole</td><td>100</td><td>4</td><td>101.94</td><td>18,837,954</td><td>86.63%</td></tr><tr><td/><td/><td>7</td><td>3/25/91</td><td>Fluconazole</td><td>100</td><td>4</td><td>81.65</td><td>15,123,020</td><td>86.66%</td></tr><tr><td/><td/><td>8</td><td>4/8/91</td><td>Fluconazole</td><td>100</td><td>4</td><td>112.53</td><td>20,778,562</td><td>86.64%</td></tr><tr><td/><td/><td>9</td><td>6/4/91</td><td>Fluconazole</td><td>100</td><td>4</td><td>113.18</td><td>22,223,228</td><td>83.20%</td></tr><tr><td/><td/><td>11</td><td>7/15/91</td><td>Fluconazole</td><td>100</td><td>4</td><td>53.28</td><td>9,896,468</td><td>87.17%</td></tr><tr><td/><td/><td>12</td><td>11/26/91</td><td>Fluconazole</td><td>200</td><td>4</td><td>96.10</td><td>18,282,472</td><td>85.54%</td></tr><tr><td/><td/><td>13</td><td>12/13/91</td><td>Fluconazole</td><td>400</td><td>32</td><td>123.67</td><td>22,070,518</td><td>89.13%</td></tr><tr><td/><td/><td>14</td><td>1/28/92</td><td>Fluconazole</td><td>400</td><td>24</td><td>98.66</td><td>18,114,916</td><td>87.41%</td></tr><tr><td/><td/><td>15</td><td>2/21/92</td><td>Clotriminazole</td><td>50</td><td>24</td><td>120.90</td><td>22,401,374</td><td>86.57%</td></tr><tr><td/><td/><td>16</td><td>4/1/92</td><td>Fluconazole</td><td>400</td><td>96</td><td>87.44</td><td>16,061,560</td><td>87.17%</td></tr><tr><td/><td/><td>17</td><td>8/25/92</td><td>Fluconazole</td><td>800</td><td>96</td><td>97.83</td><td>18,317,118</td><td>85.91%</td></tr><tr><td>Perea, S. et al.</td><td>7</td><td>412</td><td>2/15/95</td><td>Fluconazole</td><td>0</td><td>0.25</td><td>93.15</td><td>17,417,588</td><td>86.69%</td></tr><tr><td/><td/><td>2307</td><td>11/22/95</td><td>Fluconazole</td><td>400</td><td>0.75</td><td>95.79</td><td>18,014,242</td><td>85.25%</td></tr><tr><td>Perea, S. et al.</td><td>9</td><td>1002</td><td>4/20/95</td><td>Fluconazole</td><td>100</td><td>0.125</td><td>188.49</td><td>34,834,970</td><td>86.74%</td></tr><tr><td/><td/><td><named-content content-type="author-callout-style1">2823</named-content></td><td>4/6/96</td><td>Fluconazole</td><td>800</td><td/><td>282.62</td><td>52,839,288</td><td>86.30%</td></tr><tr><td/><td/><td>3795</td><td>2/26/97</td><td>Fluconazole</td><td>800</td><td>128</td><td>77.63</td><td>13,901,062</td><td>88.78%</td></tr><tr><td>Perea, S. et al.</td><td>14</td><td>580</td><td>3/13/95</td><td>Fluconazole</td><td>0</td><td>1.5</td><td>77.08</td><td>14,711,804</td><td>85.00%</td></tr><tr><td/><td/><td>2440</td><td>1/3/96</td><td>Fluconazole</td><td>800</td><td>1.5</td><td>82.93</td><td>15,446,882</td><td>85.69%</td></tr><tr><td/><td/><td><italic>2501</italic><xref ref-type="table-fn" rid="tblfn1">*</xref></td><td>1/4/96</td><td>Fluconazole</td><td>800</td><td>96</td><td>88.59</td><td>17,480,274</td><td>81.98%</td></tr><tr><td>Perea, S. et al.</td><td>15</td><td>945</td><td>4/14/95</td><td>Fluconazole</td><td>300</td><td>4</td><td>108.59</td><td>20,591,044</td><td>85.19%</td></tr><tr><td/><td/><td>1619</td><td>7/11/95</td><td>Fluconazole</td><td>500</td><td>64</td><td>93.14</td><td>17,565,080</td><td>84.69%</td></tr><tr><td>Perea, S. et al.</td><td>16</td><td>3107</td><td>6/5/96</td><td>Fluconazole</td><td>800</td><td>4</td><td>97.01</td><td>18,361,266</td><td>84.84%</td></tr><tr><td/><td/><td>3119</td><td>6/5/96</td><td>Fluconazole</td><td>800</td><td>96</td><td>87.92</td><td>16,615,462</td><td>84.67%</td></tr><tr><td/><td/><td>3120</td><td>6/5/96</td><td>Fluconazole</td><td>800</td><td>96</td><td>105.95</td><td>19,442,016</td><td>86.79%</td></tr><tr><td/><td/><td><named-content content-type="author-callout-style1">3184</named-content></td><td>7/1/96</td><td>Fluconazole</td><td>800</td><td/><td>101.89</td><td>18,487,462</td><td>87.50%</td></tr><tr><td/><td/><td><named-content content-type="author-callout-style1">3281</named-content></td><td>7/16/96</td><td>Fluconazole</td><td>800</td><td/><td>76.44</td><td>14,327,376</td><td>85.69%</td></tr><tr><td>Perea, S. et al.</td><td>30</td><td>5106</td><td>1/7/98</td><td>Fluconazole</td><td>800</td><td>0.5</td><td>87.21</td><td>16,466,524</td><td>84.67%</td></tr><tr><td/><td/><td>5108</td><td>1/7/98</td><td>Fluconazole</td><td>800</td><td>0.75</td><td>82.32</td><td>17,480,274</td><td>81.98%</td></tr><tr><td>Perea, S. et al.</td><td>42</td><td><named-content content-type="author-callout-style1">1691</named-content></td><td>8/3/95</td><td>Fluconazole</td><td>100</td><td/><td>122.60</td><td>22,072,562</td><td>88.38%</td></tr><tr><td/><td/><td>3731</td><td>12/27/96</td><td>Fluconazole</td><td>400</td><td>256</td><td>119.90</td><td>21,436,034</td><td>88.72%</td></tr><tr><td/><td/><td>3733</td><td>12/27/96</td><td>Fluconazole</td><td>400</td><td>256</td><td>95.51</td><td>17,295,888</td><td>88.00%</td></tr><tr><td>Perea, S. et al.</td><td>43</td><td>1649</td><td>7/19/95</td><td>Fluconazole</td><td>0</td><td>0.125</td><td>102.10</td><td>19,545,530</td><td>84.08%</td></tr><tr><td/><td/><td>3034</td><td>5/15/96</td><td>Fluconazole</td><td>400</td><td>0.75</td><td>92.97</td><td>17,300,040</td><td>85.64%</td></tr><tr><td>Perea, S. et al.</td><td>59</td><td>3917</td><td>2/19/97</td><td>Fluconazole</td><td>800</td><td>2</td><td>113.27</td><td>21,549,704</td><td>83.86%</td></tr><tr><td/><td/><td>4617</td><td>8/28/97</td><td>Fluconazole</td><td>400</td><td>64</td><td>75.37</td><td>15,242,904</td><td>81.42%</td></tr><tr><td/><td/><td>4639</td><td>9/2/97</td><td>Fluconazole</td><td>400</td><td>128</td><td>115.32</td><td>25,468,190</td><td>75.69%</td></tr><tr><td>Perea, S. et al.</td><td>64</td><td>4018</td><td>4/2/97</td><td>Fluconazole</td><td>200</td><td/><td>110.16</td><td>20,118,736</td><td>86.78%</td></tr><tr><td/><td/><td><named-content content-type="author-callout-style1">4380</named-content></td><td>7/14/97</td><td>Fluconazole</td><td>200</td><td/><td>18.03</td><td>20,970,946</td><td>9.26%</td></tr></tbody></table><table-wrap-foot><fn><p>Strains and coverage.</p></fn><fn><p><named-content content-type="author-callout-style1">(Red) Not clonally derived from progenitor.</named-content></p></fn><fn id="tblfn1"><label>*</label><p>isolated on same day from same patient as previously published strain, 2500.</p></fn></table-wrap-foot></table-wrap><fig-group><fig id="fig1" position="float"><object-id pub-id-type="doi">10.7554/eLife.00662.004</object-id><label>Figure 1.</label><caption><title>Overview of study design.</title><p>(<bold>A</bold>) Background, persistent, transient, recurrent, and driver mutations in patient time courses. Shown is a schematic illustration of the genomes of isolates (gray bars) from two patient time courses (Patient A and B, left and right panels, respectively), ordered from the first isolate (progenitor, top) to the last (evolved, bottom). Background mutations (purple) exist in the all isolates; persistent mutations (yellow) are not in the progenitor, but found in all subsequent isolates after their first occurrence; transient mutations (pink) are not in the progenitor and only in some later isolates; recurrently polymorphic genes contain persistent mutations that occur in the same gene in more than one patient (black box). LOH events were also evaluated for persistence (light teal bar). Driver mutations, where a new persistent homozygous allele appears (e.g., G/T > A/A), are annotated in association with persistent LOH events (dark teal) and independent of these events (not shown). Each of these can be associated with a change in phenotype, such as drug resistance (boxes, right). (<bold>B</bold>) Sampling in the context of de novo mutation and selection bottlenecks. Each strain is a single clone (circle) isolated from an evolving population (represented by a phylogenetic tree). The population evolves and undergoes selective sweeps (dashed lines), with phenotypic changes occurring during the course of infection and treatment (i.e., drug resistance, black: high, white: low; gray scale at bottom). Persistent mutations (yellow lightning bolt) have likely swept through the population, whereas transient mutations (pink lightning bolt) have not. (<bold>C</bold>) Sampling in the context of selection on existing variation. Selection acts to vary the frequency of different pre-existing genotypes in the population. Persistent mutations (yellow lightning bolt) have risen in the population to a frequency that they are repeatedly sampled (large circles) whereas transient mutations (pink lightning bolt) have not (small circle).</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.004">http://dx.doi.org/10.7554/eLife.00662.004</ext-link></p></caption><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="elife00662f001"/></fig><fig id="fig1s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.00662.005</object-id><label>Figure 1—figure supplement 1.</label><caption><title>Analysis of discordant sites.</title><p>(<bold>A</bold>) Degree of concordance (Y axis) with Sequenom iPlex genotyping for 1973 SNP X strain combinations overall (leftmost red bar; 93.9%) and in each tested strain (X axis). (<bold>B</bold>) Shown are the classes of discordant sites by genotype as defined by Illumina (orange) or Sequenome (teal) (X axis) and the prevalence (Y axis) of that genotype call in Sequenom (blue) and Illumina (orange) based discordant calls. The most common discrepancies arose when Sequenom typing classified a site as homozygous, but Illumina sequencing identified it as heterozygous. (<bold>C</bold>–<bold>G</bold>) Comparison on distributions of quality features between concordant (blue bars) and discordant (green bars) sites: (<bold>C</bold>) depth of coverage, (<bold>D</bold>) RMS Mapping Quality (MQ) score, (<bold>E</bold>) PHRED scaled quality score for each base call, shown as log-normalized ‘QUAL’ scores, (<bold>F</bold>) quality by depth (QD) score for each variant site, and (<bold>G</bold>) the allele balance ratio (AB Score) for each variant site.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.005">http://dx.doi.org/10.7554/eLife.00662.005</ext-link></p></caption><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="elife00662fs001"/></fig></fig-group><fig-group><fig id="fig2" position="float"><object-id pub-id-type="doi">10.7554/eLife.00662.006</object-id><label>Figure 2.</label><caption><title>Most isolates from the same patient are clonal.</title><p>(<bold>A</bold>) Two possible models of infection may underlie serial isolates. In the ‘clonal model’ (top) each subsequent sample (circle) is related to the other isolates. In the non-clonal model (bottom) isolates in a series are un-related. (<bold>B</bold>) The phylogenetic relationship of the isolates (black) from 11 patients (blue) was inferred based on 201,793 informative SNP positions using maximum parsimony in PAUP*. Isolates from the same patient separated by a branch distance greater than 20,000 were considered non-clonal (3281, 2823, 3184, 1691, red). Most nodes were supported by 100% of 1000 bootstrap replicates (indicated by *), expect as indicated (in gray). Clade identifiers were included as appropriate.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.006">http://dx.doi.org/10.7554/eLife.00662.006</ext-link></p><p><supplementary-material id="SD1-data"><object-id pub-id-type="doi">10.7554/eLife.00662.008</object-id><label>Figure 2—source data 1.</label><caption><title>(A) SNP category summary and all patient-series SNPs SNP category summary.</title><p>Listed for each series <bold>(PT series SNP summary</bold>) are the number of filtered (‘Materials and methods’) coding and noncoding SNPs. Coding SNPs are further classified as synonymous or nonsynonymous. Noncoding SNPs are classified as intronic, promoter region (<800 bps from the start of an ORF), or general noncoding. <bold>Patient1–Patient 59:</bold> Listed is each base that is mutated in at least one isolate in the respective series. For this base, listed are the chromosomal position, the base in the SC5314 reference genome, the base in each isolate in the series (hyphen (‘-’): homozygous, same as reference; upper case: homozygous mutation; lower case: heterozygous mutation), whether the mutation is a background mutation, transient (trans) or persistent (pers), if it is upstream, downstream or within an ORF, and in the latter case, the effect on the amino acid sequence of the encoded protein. (<bold>B</bold>) <bold>Frequency of nonsynonymous SNP occurrence between serial isolates using different filters. All SNP arising aft prev:</bold> For each clinical series (PT1-PT59) listed are the number of ORFs in each chromosome (columns) containing for each isolate (rows) all the ‘newly arising’ SNPs, defined as those not present in the immediately preceding isolate (rows). <bold>All NS in ORF aft prev:</bold> the same as above, but only for NS SNPs. All SNPs are only outside of LOH regions. <bold>All instances of Pers NS SNPs</bold>: the same as above, but only for those NS SNPs that persist once they arose. <bold>All Rec SNP aft Prev:</bold> the same as above but restricted to those ORFs that contain persistent mutations in three or more clinical series.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.008">http://dx.doi.org/10.7554/eLife.00662.008</ext-link></p></caption><media xlink:href="elife00662s001.xlsx" mimetype="application" mime-subtype="xlsx"/></supplementary-material></p></caption><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="elife00662f002"/></fig><fig id="fig2s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.00662.007</object-id><label>Figure 2—figure supplement 1.</label><caption><title>SNP heterozygosity profiles for each strain.</title><p>The heterozygosity profiles shows, in chromosomal order (top), each variant locus that exists in at least one strain in the series (white is a heterozygous SNP, blue is homozygous for the SC5314 allele, red is a homozygous SNP relative to SC5314). (<bold>A</bold>) Patient 1; (<bold>B</bold>) Patient 7; (<bold>C</bold>) Patient 9; (<bold>D</bold>) Patient 14; (<bold>E</bold>) Patient 15; (<bold>F</bold>) Patient 16; (<bold>G</bold>) Patient 30; (<bold>H</bold>) Patient 42; (<bold>I</bold>) Patient 43; (<bold>J</bold>) Patient 59. Only Patient 9 (<bold>C</bold>), Patient 16 (<bold>F</bold>), Patient 42 (<bold>H</bold>), and Patient 64 (<italic>not shown**</italic>) contain un-related isolates. ** Patient 64 contained an isolate (4380) whose genome aligned poorly to the <italic>C. albicans</italic> reference, but aligned well to <italic>C. dubliniensis.</italic></p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.007">http://dx.doi.org/10.7554/eLife.00662.007</ext-link></p></caption><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="elife00662fs002"/></fig></fig-group></p><p>The progenitor isolates were more sensitive to fluconazole than subsequent isolates, as defined by the minimum inhibitory concentration (MIC) (<xref ref-type="table" rid="tbl1">Table 1</xref>, ‘Materials and methods’). Previous studies with some of these patient isolates identified several genomic alterations that may contribute to azole resistance, including segmental aneuploidy (<xref ref-type="bibr" rid="bib72">Selmecki et al., 2006</xref>), and LOH across large chromosomal segments (<xref ref-type="bibr" rid="bib13">Coste et al., 2006</xref>; <xref ref-type="bibr" rid="bib18">Dunkel et al., 2008</xref>), as well as targeted alterations including increased expression of drug efflux genes (<xref ref-type="bibr" rid="bib13">Coste et al., 2006</xref>), mutations in ergosterol biosynthetic genes (<xref ref-type="bibr" rid="bib3">Asai et al., 1999</xref>; <xref ref-type="bibr" rid="bib55">Oliver et al., 2007</xref>), and buffering by the chaperone heat shock protein 90 (Hsp90) (<xref ref-type="bibr" rid="bib15">Cowen and Lindquist, 2005</xref>).</p><p>We sequenced the genomic DNA of the isolates as well as the <italic>C. albicans</italic> lab strain, SC5314, using Illumina sequencing (53-283X coverage, 103X on average, ‘Materials and methods’, <xref ref-type="table" rid="tbl1">Table 1</xref>) and identified in each series point mutations, LOH events and aneuploidies that were not present in the first strain in that series. By convention, all mutations were defined relative to SC5314, the <italic>C. albicans</italic> genome reference strain. We validated our pipeline for detection of point mutations using Sequenom iPlex genotyping (<xref ref-type="bibr" rid="bib78">Storm et al., 2003</xref>) (‘Materials and methods’). We interrogated 1973 SNPs in 27 isolates from nine clinical series and found that the iPlex base calls matched 1853 (93.9%, <xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1A</xref>, <xref ref-type="table" rid="tbl2">Table 2</xref>) of the calls from our computational analysis of the sequencing data. Evaluation of the discordant sites showed somewhat lower quality scores by certain metrics but did not identify any metrics that could be used to systematically revise filtering in our computational pipeline without a radical reduction in sensitivity (<xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1B–G</xref>).<table-wrap id="tbl2" position="float"><object-id pub-id-type="doi">10.7554/eLife.00662.009</object-id><label>Table 2.</label><caption><p>Sequenom iPLEX genotyping assay validation</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.009">http://dx.doi.org/10.7554/eLife.00662.009</ext-link></p></caption><table frame="hsides" rules="groups"><thead><tr><th>Patient</th><th>Isolate</th><th>Total discordant</th><th>Total concordant</th><th>Total Assayed</th><th>% Concordant</th></tr></thead><tbody><tr><td>Patient_1</td><td>TWTC1</td><td align="right">2</td><td align="right">31</td><td align="right">33</td><td align="right">93.94%</td></tr><tr><td>Patient_1</td><td>TWTC2</td><td align="right">1</td><td align="right">32</td><td align="right">33</td><td align="right">96.97%</td></tr><tr><td>Patient_1</td><td>TWTC3</td><td align="right">1</td><td align="right">32</td><td align="right">33</td><td align="right">96.97%</td></tr><tr><td>Patient_1</td><td>TWTC12</td><td align="right">1</td><td align="right">32</td><td align="right">33</td><td align="right">96.97%</td></tr><tr><td>Patient_1</td><td>TWTC13</td><td align="right">1</td><td align="right">32</td><td align="right">33</td><td align="right">96.97%</td></tr><tr><td>Patient_1</td><td>TWTC15</td><td align="right">1</td><td align="right">32</td><td align="right">33</td><td align="right">96.97%</td></tr><tr><td>Patient_1</td><td>TWTC16</td><td align="right">1</td><td align="right">31</td><td align="right">32</td><td align="right">96.88%</td></tr><tr><td>Patient_1</td><td>TWTC17</td><td align="right">1</td><td align="right">32</td><td align="right">33</td><td align="right">96.97%</td></tr><tr><td>Patient_7</td><td>412</td><td align="right">3</td><td align="right">60</td><td align="right">63</td><td align="right">95.24%</td></tr><tr><td>Patient_7</td><td>2307</td><td align="right">4</td><td align="right">59</td><td align="right">63</td><td align="right">93.65%</td></tr><tr><td>Patient_9</td><td>1002</td><td align="right">16</td><td align="right">96</td><td align="right">112</td><td align="right">85.71%</td></tr><tr><td>Patient_9</td><td>3795</td><td align="right">9</td><td align="right">103</td><td align="right">112</td><td align="right">91.96%</td></tr><tr><td>Patient_14</td><td>580</td><td align="right">3</td><td align="right">49</td><td align="right">52</td><td align="right">94.23%</td></tr><tr><td>Patient_14</td><td>2440</td><td align="right">2</td><td align="right">27</td><td align="right">29</td><td align="right">93.10%</td></tr><tr><td>Patient_14</td><td>2501</td><td align="right">3</td><td align="right">33</td><td align="right">36</td><td align="right">91.67%</td></tr><tr><td>Patient_15</td><td>945</td><td align="right">8</td><td align="right">121</td><td align="right">129</td><td align="right">93.80%</td></tr><tr><td>Patient_15</td><td>1619</td><td align="right">10</td><td align="right">120</td><td align="right">130</td><td align="right">92.31%</td></tr><tr><td>Patient_16</td><td>3107</td><td align="right">2</td><td align="right">51</td><td align="right">53</td><td align="right">96.23%</td></tr><tr><td>Patient_16</td><td>3119</td><td align="right">3</td><td align="right">50</td><td align="right">53</td><td align="right">94.34%</td></tr><tr><td>Patient_16</td><td>3120</td><td align="right">2</td><td align="right">50</td><td align="right">52</td><td align="right">96.15%</td></tr><tr><td>Patient_30</td><td>5106</td><td align="right">3</td><td align="right">215</td><td align="right">218</td><td align="right">98.62%</td></tr><tr><td>Patient_30</td><td>5108</td><td align="right">19</td><td align="right">204</td><td align="right">223</td><td align="right">91.48%</td></tr><tr><td>Patient_43</td><td>1649</td><td align="right">7</td><td align="right">89</td><td align="right">96</td><td align="right">92.71%</td></tr><tr><td>Patient_43</td><td>3034</td><td align="right">8</td><td align="right">88</td><td align="right">96</td><td align="right">91.67%</td></tr><tr><td>Patient_59</td><td>3917</td><td align="right">3</td><td align="right">62</td><td align="right">65</td><td align="right">95.38%</td></tr><tr><td>Patient_59</td><td>4617</td><td align="right">2</td><td align="right">63</td><td align="right">65</td><td align="right">96.92%</td></tr><tr><td>Patient_59</td><td>4639</td><td align="right">4</td><td align="right">59</td><td align="right">63</td><td align="right">93.65%</td></tr><tr><td colspan="2">TOTAL</td><td align="right">120</td><td align="right">1853</td><td align="right">1973</td><td align="right">93.92%</td></tr></tbody></table></table-wrap></p></sec><sec id="s2-2"><title>Most series are clonal, but there is significant genetic diversity between isolates</title><p>We designated as <italic>background</italic> those polymorphisms those that are common to all isolates in a series, including the first (‘progenitor’) isolate (<xref ref-type="fig" rid="fig1">Figure 1A</xref>, purple) and use them to determine that isolates within most series were clonally related, suggesting a single (primary) infection source (<xref ref-type="fig" rid="fig1">Figure 1B,C</xref>, <xref ref-type="fig" rid="fig2">Figure 2</xref>, <xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1</xref>, ‘Materials and methods’). To distinguish between a single primary (clonal) infection (<xref ref-type="fig" rid="fig2">Figure 2A</xref>, top) and repeated, independent infections (<xref ref-type="fig" rid="fig2">Figure 2A</xref>, bottom), we determined the distance between every two isolates based on their SNP profile and used as a heuristic a neighbor-joining algorithm to construct a phylogenetic tree from this distance metric (‘Materials and methods’, <xref ref-type="fig" rid="fig2">Figure 2B</xref>). Patient 64 contained one <italic>C. albicans</italic> isolate (4018) and one C. <italic>dubliniensis</italic> isolate (4380); therefore, we have excluded this series from further analysis. Additionally, we detected at least one non-clonal <italic>C. albicans</italic> isolate in three of the remaining ten patient series (PT 9,16, 42; <xref ref-type="fig" rid="fig2">Figure 2B</xref>, red), indicating that at least ∼36% of the 11 patients sampled carried more than one unrelated <italic>Candida</italic> strain<italic>.</italic> We removed the four non-clonal samples (<xref ref-type="fig" rid="fig2">Figure 2B</xref>, red) from further consideration, and all subsequent analyses focused on samples from the 10 patients with at least two clonal isolates.</p><p>Despite these clonal relationships, the distance between isolates indicated significant genetic diversity <italic>within</italic> each patient series (<xref ref-type="fig" rid="fig2">Figure 2B</xref>), typically with each isolate differing by several thousand SNPs from its ‘progenitor’ isolate (<xref ref-type="supplementary-material" rid="SD1-data">Figure 2—source data 1</xref>). These data are consistent with two different evolutionary scenarios: accumulation of de novo mutations followed by selection (<xref ref-type="fig" rid="fig1">Figure 1B</xref>), or selection acting on pre-existing variation to vary the frequency of different genotypes in the population (<xref ref-type="fig" rid="fig1">Figure 1C</xref>). The large number of SNPs detected suggests that isolates from later time points in a series are not simply direct descendants of the earlier isolate; however, since mutation and mitotic recombination rates can be elevated under stressful conditions (e.g., drug treatment <xref ref-type="bibr" rid="bib26">Galhardo et al., 2007</xref>; <xref ref-type="bibr" rid="bib20">Forche et al., 2011</xref>), we cannot rule out the possibility that some of the variation may be due to de novo events occurring between time points. Formally distinguishing between these two models is not possible with the samples and data at hand. However, the role of pre-existing diversity is supported by the observation that different isolates collected on the same day from the same patient (patient 14 [2440 and 2501] and patient 16 [3107 and 3119]) differed by 9668 and 18,291 SNPs, respectively (<xref ref-type="supplementary-material" rid="SD1-data">Figure 2—source data 1</xref>) and had very different fluconazole MIC levels (<xref ref-type="table" rid="tbl1">Table 1</xref>) and different fitness phenotypes (see below), although in each case the strains were clearly genetically related (<xref ref-type="fig" rid="fig2">Figure 2B</xref>). Thus, we conclude that a population of related but divergent genotypes of the same lineage exists within a given patient. We next sought to identify potentially adaptive genetic changes by focusing on large-scale events (LOH and aneuploidies) as well as single-nucleotide polymorphisms.</p></sec><sec id="s2-3"><title>Genetic alterations absent from the progenitor isolate, persistent within a patient, and recurrent across patients are likely adaptive</title><p>Given the high number of SNPs, LOH events and aneuploidies, we next devised a strategy to identify those changes that are more likely to play an adaptive role in drug resistance and host adaptation. We previously filtered all <italic>background</italic> polymorphisms, defined as any SNP relative to the reference present in all isolates from a series. Next, we defined alterations as <italic>persistent</italic> if present within the same patient at all subsequent time points after the ‘non-progenitor’ isolate in which they are first identified. We reasoned that such persistent changes will include those variants that were driven to sufficiently high frequency by selection to ensure repeated sampling (<xref ref-type="fig" rid="fig1">Figure 1B,C</xref>, yellow lightning bolt), whereas non-persistent (transient) ones do not (<xref ref-type="fig" rid="fig1">Figure 1B,C</xref>, pink lightning bolt). We consider the special case of a genetic change detected only in the endpoint isolate as ‘persistent’ as well, since several of the time courses consist of only two or three isolates. We apply the persistence filter to better identify potentially adaptive aneuploidies, LOH events, and SNPs.</p><p>Next, we further focused on non-synonymous polymorphisms in coding regions and employed two different strategies to identify potentially adaptive changes. In the first strategy, to identify potential <italic>drivers</italic> of adaptation, we focused on non-synonymous SNPs that were homozygous for a genotype not found in the progenitor strain that persisted in the subsequent isolates (e.g., G/T > A/A) consistent with positive selection. In the second strategy, we analyzed genes that were <italic>recurrently</italic> polymorphic across patients, such that persistent, non-synonymous polymorphisms appeared within the same open reading frame (ORF) in different patient series (<xref ref-type="fig" rid="fig1">Figure 1A</xref> and <xref ref-type="supplementary-material" rid="SD4-data">Figure 5—source data 1A</xref>). For recurrence, we considered only those that were not included in LOH regions, as these regions artificially inflate the estimates of persistence and recurrence. Recurrence allows us to better handle polymorphisms from the endpoint isolate in a series for which ‘persistence’ does not provide a meaningful filter. Thus, we further considered polymorphisms occurring only in the terminal isolate in one patient if polymorphisms also recurred in the same ORF in a series from two other patients. For example, filtering for both persistence and recurrence across at least three series reduced the number of polymorphisms for patient 1 from 13,562 polymorphisms in 5022 genes to 23 recurrent genes (<xref ref-type="supplementary-material" rid="SD1-data">Figure 2—source data 1</xref>, <xref ref-type="supplementary-material" rid="SD4-data">Figure 5—source data 1A</xref>).</p></sec><sec id="s2-4"><title>LOH events are commonly associated with increased resistance</title><p>LOH events were detected in all of the series and were often persistent, recurrent, and associated with increased drug resistance (<xref ref-type="fig" rid="fig3 fig4">Figures 3 and 4</xref>, <xref ref-type="supplementary-material" rid="SD2-data">Figure 3—source data 1</xref>). For example, three of four LOH events in Patient 1 were persistent and associated with an increase in MIC and both of these events were recurrent, such that LOH events in these genomic regions coincided with increases in MIC in other patients. Highly recurrent LOH events occurred on the right arm of chromosome 3 (in Patients 1, 9, 14, 16, 42, and 59; <xref ref-type="fig" rid="fig3">Figure 3A</xref>, <xref ref-type="fig" rid="fig4">Figure 4A,B,D,F,H</xref>, <xref ref-type="supplementary-material" rid="SD2-data">Figure 3—source data 1</xref>) and on the left arm of chromosome 5 (in Patients 1, 14, 15, and 43; <xref ref-type="fig" rid="fig3">Figure 3A</xref>, <xref ref-type="fig" rid="fig4">Figure 4B,C,G</xref>, <xref ref-type="supplementary-material" rid="SD2-data">Figure 3—source data 1</xref>). These regions include key genes implicated in drug resistance: on Chromosome 3, genes encoding the Cdr1 and Cdr2 efflux pumps and the Mrr1 transcription factor that regulates the Mdr1 major facilitator superfamily efflux pump (<xref ref-type="bibr" rid="bib71">Schubert et al., 2011</xref>), and on Chromosome 5, genes encoding the drug target Erg11, and Tac1, a transcription factor that positively regulates expression of <italic>CDR1</italic> and <italic>CDR2</italic> (<xref ref-type="bibr" rid="bib13">Coste et al., 2006</xref>). The extent of persistence and recurrence of these two LOH events is statistically significant under a naïve binary model (p < 5 × 10<sup>−4</sup> for the Chr3R LOH; p < 0.01 for the Chr5L LOH). The recurrence of LOH events that coincide with changes in MIC suggests that they have been positively selected to rise in frequency relative to the progenitor strain. Notably, some of the recurrent LOH events may have been difficult to detect previously on SNP arrays (<xref ref-type="bibr" rid="bib21">Forche et al., 2008</xref>; <xref ref-type="bibr" rid="bib22">Forche et al., 2004</xref>; <xref ref-type="bibr" rid="bib24">Forche et al., 2005</xref>) due to the relative paucity of SNPs in those regions in the reference strain, SC5413, itself a clinical isolate.<fig id="fig3" position="float"><object-id pub-id-type="doi">10.7554/eLife.00662.010</object-id><label>Figure 3.</label><caption><title>LOH events were often persistent while aneuploidies were often transient.</title><p>For each time series shown are the genomes of all isolates (rows) from a patient, ordered from the first isolate (progenitor, top) to the last (evolved, bottom). Boxes on right indicate the MIC of the respective strain (black: high, white: low; gray scale at bottom). Persistent LOHs: blue, transient LOHs: pink; trisomies (all transient): green. The sequence coverage along each chromosome is indicated by black tickmarks. (<bold>A</bold>) Patient 1 has four LOH events, each coinciding with an increase in MIC (gray scale boxes, right). One LOH is transient (isolate 2, chromosome R, pink) and three are persistent (isolate 3, chromosome 3; isolate 13, chromosome 5; and isolate 16, chromosome 5, blue). The ploidy changes (isolates 6, 8, 13) are all transient. (<bold>B</bold>) Patient 7 has one LOH event (isolate 2307, chromosome 3, blue) which coincides with an increase in MIC.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.010">http://dx.doi.org/10.7554/eLife.00662.010</ext-link></p><p><supplementary-material id="SD2-data"><object-id pub-id-type="doi">10.7554/eLife.00662.011</object-id><label>Figure 3—source data 1.</label><caption><title>Persistent LOH regions LOH map.</title><p>For each isolate (strain column) in each series (patient column), listed are the coordinates of any persistent LOH in that isolate. Coordinates in blue are persistent LOH events, coordinates in red are transient LOH events.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.011">http://dx.doi.org/10.7554/eLife.00662.011</ext-link></p></caption><media xlink:href="elife00662s002.xlsx" mimetype="application" mime-subtype="xlsx"/></supplementary-material></p></caption><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="elife00662f003"/></fig><fig id="fig4" position="float"><object-id pub-id-type="doi">10.7554/eLife.00662.012</object-id><label>Figure 4.</label><caption><title>Persistent and transient LOH and aneuploidies.</title><p>For each time series shown are the genomes of all isolates (rows), ordered from the first isolate (progenitor, top) to the last (evolved, bottom). Boxes on right indicate the MIC of the respective strain (black: high, white: low, gray scale at bottom). Persistent LOHs: light blue, transient LOHs: pink; trisomies (all transient): green. The coverage along each chromosome is indicated by black tickmarks. (<bold>A</bold>) Patient 9; (<bold>B</bold>) Patient 14; (<bold>C</bold>) Patient 15; (<bold>D</bold>) Patient 16; (<bold>E</bold>) Patient 30; (<bold>F</bold>) Patient 42; (<bold>G</bold>) Patient 43; (<bold>H</bold>) Patient 59. Several LOHs are recurrent (right arm of chromosome 3, left arm of chromosome 5, and chromosome 1). Please note: data in <xref ref-type="supplementary-material" rid="SD2-data">Figure 3—source data 1</xref> also applies to this figure.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.012">http://dx.doi.org/10.7554/eLife.00662.012</ext-link></p><p><supplementary-material id="SD3-data"><object-id pub-id-type="doi">10.7554/eLife.00662.013</object-id><label>Figure 4—source data 1.</label><caption><title>Persistent LOH regions LOH map.</title><p>For each isolate (strain column) in each series (patient column), listed are the coordinates of any persistent LOH in that isolate. Coordinates in blue are persistent LOH events, coordinates in red are transient LOH events.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.013">http://dx.doi.org/10.7554/eLife.00662.013</ext-link></p></caption><media xlink:href="elife00662s003.xlsx" mimetype="application" mime-subtype="xlsx"/></supplementary-material></p></caption><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="elife00662f004"/></fig></p><p>Putative driver mutations (defined above as non-synonymous SNPs that were homozygous for a genotype not found in the progenitor strain that persisted in the subsequent isolates; e.g.<italic>,</italic> G/T > A/A) in these regions are suggestive of a point mutation followed by an LOH of the mutant allele that confers an advantage. There were 131 such mutations in 86 ORFs from 18 LOH regions from the 10 clonal patient series (<xref ref-type="table" rid="tbl1">Table 1</xref> and <xref ref-type="fig" rid="fig5">Figure 5</xref>; <xref ref-type="supplementary-material" rid="SD4-data">Figure 5—source data 1B,C</xref>). Some of the SNPs were in genes that encode proteins with key known roles in drug resistance and were associated with large LOH events. For example, a nonsynonymous homozygous change in the fluconazole drug target <italic>ERG11</italic> was associated with the formation of the persistent LOH on the left arm of chromosome 5 in Patient 1 (<xref ref-type="fig" rid="fig3">Figure 3A</xref>), consistent with previous reports (<xref ref-type="bibr" rid="bib86">White, 1997b</xref>), as was <italic>TAC1</italic> in Patient 42. In another example, the persistent and recurrent LOH on the right arm of chromosome 3 in Patients 9 and 16 (<xref ref-type="fig" rid="fig4">Figure 4A,D</xref>) was associated with the presence of a homozygous mutation in <italic>MRR1</italic> (<xref ref-type="bibr" rid="bib71">Schubert et al., 2011</xref>), a regulator of <italic>MDR1</italic> expression. Other mutations were in genes not previously related to fluconazole resistance, including cell adhesion (<italic>ALS3,5</italic> and <italic>7</italic> and <italic>HYR3;</italic> [<xref ref-type="bibr" rid="bib35">Hoyer et al., 1998</xref>; <xref ref-type="bibr" rid="bib76">Sheppard et al., 2004</xref>; <xref ref-type="bibr" rid="bib34">Hoyer et al., 2008</xref>]), filamentous growth (<italic>FGR14, FGR28,</italic> and <italic>EFH,</italic> [<xref ref-type="bibr" rid="bib83">Uhl et al., 2003</xref>; <xref ref-type="bibr" rid="bib12">Connolly et al., 2013</xref>]), and biofilm formation (<italic>BCR1</italic> and <italic>YAK1;</italic> [<xref ref-type="bibr" rid="bib53">Nobile and Mitchell, 2005</xref>; <xref ref-type="bibr" rid="bib29">Goyard et al., 2008</xref>; <xref ref-type="bibr" rid="bib54">Noble et al., 2010</xref>; <xref ref-type="bibr" rid="bib52">Nobile et al., 2012</xref>]). Thus, the detection of known genes involved in drug resistance confirms the approach works and that detection of genes involved in processes implicated in virulence, suggests that these process are co-evolving.<fig-group><fig id="fig5" position="float"><object-id pub-id-type="doi">10.7554/eLife.00662.014</object-id><label>Figure 5.</label><caption><title>Co-occurrence of nonsynonymous substitutions across isolates reveals functional clusters.</title><p>(<bold>A</bold>) For each of the recurrently mutated 240 genes (genes in which nonsynonymous persistent SNPs appear in more than three patients and are not within an LOH region), we constructed a patient-by-gene binary vector. We clustered the resulting patient-by-gene matrix using NMF clustering to reveal five coherent clusters (correlation matrix of the clusters left; red: positive correlation; blue: negative correlation; white: no correlation). (<bold>B</bold>) Co-occurrence clusters. For the genes in each cluster (rows), shown are their mutated occurrences in each patient (columns); green: gene is persistently mutated in patient, white: no persistent mutation, yellow circle: driver mutation. Functional enrichment of clusters was revealed using gene ontology, and genes matching the enriched cluster function are bolded. We have overlaid recurrent driver mutations (e.g., G/T > A/A) (n = 17) occurring outside of LOH regions (yellow circle, green box) and inside LOH regions (yellow circle, white box).</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.014">http://dx.doi.org/10.7554/eLife.00662.014</ext-link></p><p><supplementary-material id="SD4-data"><object-id pub-id-type="doi">10.7554/eLife.00662.016</object-id><label>Figure 5—source data 1.</label><caption><title>(A) Recurrence lists and clusters. 1 All Pers NS Genes.</title><p>Listed are all the ORFs with a persistent nonsynonymous SNPs, the series in which they occur as such (1 in relevant Patient 1-Patient 59 column), and the total number of series in which they recur (SUM column). <bold>2 All Pers NS in LOH:</bold> Listed are all the ORFs with a persistent nonsynonymous SNPs within an LOH region. <bold>3 All Pers NS not in LOH:</bold> Listed are all the ORFs with a persistent nonsynonymous SNPs NOT within an LOH region. <bold>4 Cluster Rec. genes not in LOH:</bold> NMF Clustering of the occurrence matrix from ‘All Pers NS not in LOH’. <bold>5 Cluster GO Enrichment:</bold> The GO enrichments for each of the clusters identified in ‘4 Cluster Rec. genes not in LOH’. (<bold>B</bold>) <bold>Driver mutations. Patient 1—59</bold>. Shown are all the positions where a nonsynonymous SNP changed from one homozygous genotype to another. Each column represents the base-call in that isolate of a given patient series. The formatting is consistent with <xref ref-type="supplementary-material" rid="SD1-data">Figure 2—source data 1</xref>. <bold>Drivers Recurrence in genome</bold>: for each of the driver candidates identified in the previous tabs, shown are the occurrence of a driver mutation in that ORF across each of the patient series. (<bold>C</bold>) <bold>Driver mutations in LOH regions.</bold> As above (<xref ref-type="supplementary-material" rid="SD4-data">Figure 5—source data 1B</xref>), but restricted to only driver mutations occurring within LOH regions. (<bold>D</bold>). <bold>Recurrence lists and clusters for MIC associated mutations</bold>. As above (<xref ref-type="supplementary-material" rid="SD4-data">Figure 5—source data 1A</xref>), but restricted to only recurrent mutations that occur in parallel with changes in MIC.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.016">http://dx.doi.org/10.7554/eLife.00662.016</ext-link></p></caption><media xlink:href="elife00662s004.xlsx" mimetype="application" mime-subtype="xlsx"/></supplementary-material></p></caption><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="elife00662f005"/></fig><fig id="fig5s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.00662.015</object-id><label>Figure 5—figure supplement 1.</label><caption><title>Co-occurrence of nonsynonymous SNPs occurring in conjunction with a shift in MIC.</title><p>(<bold>A</bold>) For each of the 166 recurrently mutated genes associated with a change in MIC, we constructed a patient-by-gene binary vector. We clustered the resulting patient by gene matrix using NMF clustering to reveal 5 coherent clusters (correlation matrix of the clusters left; red: positive correlation; blue: negative correlation; white: no correlation). (<bold>B</bold>) Co-occurrence clusters. For the genes in each cluster (rows), shown are their mutated occurrences in each patient (columns); green: gene is persistently mutated in patient, white: no persistent mutation. Functional enrichment of clusters was revealed using gene ontology, and genes matching the enriched cluster function are bolded.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.015">http://dx.doi.org/10.7554/eLife.00662.015</ext-link></p></caption><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="elife00662fs003"/></fig></fig-group></p></sec><sec id="s2-5"><title>Aneuploidies are not predictive of MIC, but may facilitate the appearance of drug resistance</title><p>Aneuploidies, either whole chromosomal or segmental, were evident in at least one isolate from 80% (8/10) of the clonal patient series, with the most prevalent aneuploidies involving Chromosome 5 (6 of 8 patients with at least one aneuploid isolate; <xref ref-type="fig" rid="fig3">Figure 3A</xref> and <xref ref-type="fig" rid="fig4">Figure 4B–G</xref>, green). In contrast to the persistence of most LOH events, persistent aneuploidies were rarer and were not consistently associated with adaptive increases in MIC levels (<xref ref-type="fig" rid="fig3 fig4">Figures 3 and 4</xref>). This is consistent with the irreversibility of LOH events in the absence of mating highlights the reversible nature (instability) of aneuploidy chromosomes.</p><p>While we cannot definitively infer an ordering of events from our singly sampled isolates, we hypothesize that aneuploidy could contribute to the evolution of LOH by increasing the likelihood of its occurrence. For example, in 4 of the 6 patients with a Chromosome 5 LOH, the isolate with an LOH event also harbors a Chromosome 5 trisomy or is preceded by an isolate with a Chromosome 5 trisomy. Thus, the additional copy may increase the likelihood of an LOH event on that chromosome. In three of these cases, <italic>ERG11,</italic> located on the region of Chromosome 5 with LOH, was mutated. Additionally, isolates in 2 of the 7 patients with a Chromosome 3 LOH were trisomic for this chromosome.</p></sec><sec id="s2-6"><title>Persistent SNPs in 240 recurrently polymorphic genes identify targets likely associated with drug resistance and host adaptation</title><p>We identified persistent nonsynonymous coding SNPs within 1470 genes outside LOH tracts, 167 of them harboring 336 driver-like polymorphisms (<xref ref-type="supplementary-material" rid="SD4-data">Figure 5—source data 1B</xref>). These again include <italic>ERG11</italic> in patients 9, 14, 30, and 59 and <italic>TAC1</italic> in patients 1, 7, 14, 15, 30 and 43 (<xref ref-type="fig" rid="fig5">Figure 5</xref>). Applying the recurrence filter (i.e., persistent nonsynonymous SNPs that appeared in the same ORFs in three or more patient series), we identified 240 polymorphic genes that are more likely to have contributed to adaptation (<xref ref-type="supplementary-material" rid="SD4-data">Figure 5—source data 1A</xref>). This number of genes is higher than expected by chance (empirical p < 10<sup>−4</sup> based on a Poisson model of background mutation, ‘Materials and methods’). Though the coding sequence for these 240 recurrent genes is longer than average (2.21 ± 1.53 kb vs 1.83 ± 1.29 kb for non-recurrent persistent genes, p < 3.68 × 10<sup>−5</sup>, <italic>t</italic>-test), and thus a larger target for mutation, our simulation accounts for gene length. Notably, 17 persistent recurrently polymorphic genes also had driver-like polymorphisms, eight of which were also homozygosed in an LOH tract in at least one patient series (<xref ref-type="fig" rid="fig5">Figure 5</xref>, <xref ref-type="supplementary-material" rid="SD4-data">Figure 5—source data 1A,B,C</xref>). Finally, polymorphisms in 166 of the 240 genes appeared together with an increase in MIC and are thus stronger candidates for making a significant functional contribution to resistance (<xref ref-type="fig" rid="fig5s1">Figure 5—figure supplement 1</xref>, <xref ref-type="supplementary-material" rid="SD4-data">Figure 5—source data 1D</xref>, empirical p < 10<sup>−5</sup> based a binomial model, ‘Materials and methods’).</p><p>The set of 240 recurrently mutated genes was enriched for fungal-type cell wall (18 genes, p < 0.0012) and cell surface genes (24 genes, p < 0.00012), including several members in each of three cell wall gene families important for biofilm formation and virulence (<xref ref-type="bibr" rid="bib34">Hoyer et al., 2008</xref>): the Hyr/Iff proteins (<italic>HYR1</italic> and <italic>3</italic>, <italic>IFF8</italic> and <italic>6</italic>), the <italic>ALS</italic> adhesins (<italic>ALS1</italic>-<italic>4</italic>,<italic>7</italic>,<italic>9</italic>), and the <italic>PGA-30</italic>-like proteins (seven genes) (<xref ref-type="supplementary-material" rid="SD4-data">Figure 5—source data 1A</xref>). All three families are specifically expanded in the genomes of pathogenic <italic>Candida</italic> species (<xref ref-type="bibr" rid="bib8">Butler et al., 2009</xref>). In addition, seven members of the <italic>FGR</italic> genes (<xref ref-type="bibr" rid="bib83">Uhl et al., 2003</xref>), involved in filamentous growth and specifically expanded in <italic>C. albicans</italic> (<xref ref-type="bibr" rid="bib8">Butler et al., 2009</xref>), are also among the 240 genes (<xref ref-type="supplementary-material" rid="SD4-data">Figure 5—source data 1A</xref>).</p><p>The most recurrently mutated gene outside of an LOH region was <italic>AXL1</italic> that encodes a putative endoprotease, whose transcript is upregulated in an RHE model of oral candidiasis and in clinical isolates from HIV+ patients with oral candidiasis (<xref ref-type="bibr" rid="bib90">Zakikhany et al., 2007</xref>). The gene is persistently mutated in eight series, (three of which were driver mutations), followed by ten genes mutated in seven series (<xref ref-type="supplementary-material" rid="SD4-data">Figure 5—source data 1A</xref>). <italic>ERG11</italic>, which encodes the drug target of fluconazole, was affected in 70% (7/10) of the patient series with persistent SNPs in four series (Patients 9, 14, 30, and 59) and mutations in the LOH events in three series (Patients 1, 15, and 43) (<xref ref-type="supplementary-material" rid="SD4-data">Figure 5—source data 1A</xref>). Likewise <italic>HYR3</italic>, a known virulence gene, was persistently mutated in nine of the patients, three of which occurred in LOH tracts, including one in which a new allele was homozygosed (<xref ref-type="fig" rid="fig5">Figure 5</xref>, <xref ref-type="supplementary-material" rid="SD4-data">Figure 5—source data 1A,B,C</xref>). More generally, 171 of the 240 genes were also mutated in an LOH tract in at least one additional patient (15/171 in three or more additional patients and 34/171 in two additional patients).</p><p>Next, we partitioned the 240 recurrently mutated genes into 5 ‘co-occurrence clusters’ based on the correlation in their mutation occurrence patterns (<xref ref-type="fig" rid="fig5">Figure 5</xref>, <xref ref-type="supplementary-material" rid="SD4-data">Figure 5—source data 1A</xref>). These correlations are significantly higher than expected in a null model (p < 5.2 × 10<sup>−182</sup>, permutation test, ‘Materials and methods’). The characterized genes in most of the clusters have coherent functions. Cluster 1 is enriched for cell wall and cell surface genes, Cluster 2 for cell cycle and stress genes, Cluster 3 for genes involved in drug response, and Cluster 5 for carbohydrate binding (<xref ref-type="fig" rid="fig5">Figure 5</xref>, <xref ref-type="supplementary-material" rid="SD4-data">Figure 5—source data 1A</xref>). Most of the genes in these clusters are not well characterized and represent new candidates involved drug resistance and adaptation to the host environment. The full list of genes and descriptions is given in <xref ref-type="supplementary-material" rid="SD4-data">Figure 5—source data 1A</xref>.</p></sec><sec id="s2-7"><title>Changes in virulence phenotypes in evolved drug-resistant isolates</title><p>To explore the possibility that some of the mutations reflect adaptation to other factors besides drug, we next measured phenotypes associated to virulence and interaction with the host (‘Materials and methods’). Adhesion, filamentation, and virulence in a <italic>C. elegans</italic> model of infection (<xref ref-type="bibr" rid="bib38">Jain et al., 2013</xref>) were measured for a large panel of isolates (<xref ref-type="fig" rid="fig6">Figure 6</xref>, <xref ref-type="fig" rid="fig6s1">Figure 6—figure supplement 1</xref>, <xref ref-type="supplementary-material" rid="SD5-data">Figure 6—source data 1</xref>). Additionally, we measured competitive fitness in standard tissue culture medium (RPMI) with and without drug in vitro (<xref ref-type="fig" rid="fig7">Figure 7</xref>).<fig-group><fig id="fig6" position="float"><object-id pub-id-type="doi">10.7554/eLife.00662.017</object-id><label>Figure 6.</label><caption><title>Filamentation, adhesion and virulence increase concurrently with fitness.</title><p>For each pair of consecutive isolates (green preceding blue), shown are the fitness, adhesion, filamentation, and virulence in a worm model of infection (each described in ‘Materials and methods’). A subset of fitness values are duplicated from <xref ref-type="fig" rid="fig7">Figure 7A</xref>, with selection coefficient (s) shown on the Y-axis. A subset of adhesion values are plotted from <xref ref-type="supplementary-material" rid="SD5-data">Figure 6—source data 1</xref>, with Abs590 nm on the Y-axis. A subset of images showing filamentation on spider media are shown, with the full set found in <xref ref-type="fig" rid="fig6s1">Figure 6—figure supplement 1</xref>. For virulence, shown are Kaplan–Meier plots of survival rates from <italic>C.elegans</italic> infection with the specified <italic>C. albicans</italic> isolates (‘Materials and methods’). For each isolate pair, significant changes in virulence were observed between the two isolates (in all cases, p < 0.001, log-rank test), with three of the four evolved isolates being more virulent than their corresponding progenitor. (<bold>A</bold>) Patient 30 isolates 5106 and 5108; (<bold>B</bold>) Patient 43 isolates 1649 and 3034; (<bold>C</bold>) Patient 1 isolates 12 and 13; (<bold>D</bold>) Patient 59 isolates 3917 and 4617.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.017">http://dx.doi.org/10.7554/eLife.00662.017</ext-link></p><p><supplementary-material id="SD5-data"><object-id pub-id-type="doi">10.7554/eLife.00662.019</object-id><label>Figure 6—source data 1.</label><caption><title>Adhesion values for the majority of isolates.</title><p>Adhesion was defined as described in ‘Materials and methods’ and measured eight times to determine the average adherence as measured by Abs590.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.019">http://dx.doi.org/10.7554/eLife.00662.019</ext-link></p></caption><media xlink:href="elife00662s005.xlsx" mimetype="application" mime-subtype="xlsx"/></supplementary-material></p></caption><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="elife00662f006"/></fig><fig id="fig6s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.00662.018</object-id><label>Figure 6—figure supplement 1.</label><caption><title>Filamentation increases in many patient series.</title><p>For several patient series, shown are the filamentation assay results after 7 days of growth on Spider Media (‘Materials and Methods’). These data, a subset of which is shown in <xref ref-type="fig" rid="fig6">Figure 6</xref>, demonstrate the heterogeneity seen between strains, as well as the general trend for filamentation to increase over time.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.018">http://dx.doi.org/10.7554/eLife.00662.018</ext-link></p></caption><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="elife00662fs004"/></fig></fig-group><fig id="fig7" position="float"><object-id pub-id-type="doi">10.7554/eLife.00662.020</object-id><label>Figure 7.</label><caption><title>Emergence of increased drug resistance often coincides with reduction in fitness in the absence of drug, but an increase in the presence of drug.</title><p>(<bold>A</bold>) For each patient (panel) shown is the fitness (‘Materials and methods’) of each strain (Y axis, mean ± STDV), ordered from the progenitor to evolved isolates (left to right, X axis). Fitness is calculated relative to an ENO1::YFP SC5314 reference isolate. The MIC of each strain is shown in the gray boxes on top (white: low; black: high, color bar at bottom). Green: isolates with aneuploidies; Blue: euploid isolates. (<bold>B</bold>) Shown is the mean difference between fitness in the absence and presence of drug (Y axis, error bars are ± STDV; n <underline>></underline> 3) for isolates (X axis) that showed a decrease in fitness (<xref ref-type="fig" rid="fig7">Figure 7A</xref>) in the absence of drug concomitant with an increase in MIC (asterisks), and flanking isolates in Patient 1 and 59 (ordered from the progenitor to evolved isolates, left to right, X axis). The difference in fitness is calculated as the difference in selection coefficient (<italic>s</italic>, Y axis) between matching competition experiments in RPMI and those in RPMI with one half the MIC for fluconazole (<xref ref-type="table" rid="tbl1">Table 1</xref>) for each isolate tested (X axis). Negative values indicate that the strain had higher fitness in the presence of fluconazole vs assays without fluconazole. For each assay, the fluconazole-resistant isolate 4639 <italic>ENO1::YFP</italic> was used as the reference strain.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.020">http://dx.doi.org/10.7554/eLife.00662.020</ext-link></p></caption><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="elife00662f007"/></fig></p><p>We found substantial variation in many of these phenotypes between isolates in the same series (<xref ref-type="fig" rid="fig6">Figure 6</xref> and <xref ref-type="fig" rid="fig7">Figure 7</xref>), supporting the notion that the isolates are samples from a broad range of genetic variants within clonal (single infection) populations. In general, increased fitness in vitro (in the absence of drug) correlated with an increase in traits associated with virulence (adhesion, filamentation, and virulence in nematode). For example, the later isolates in the series from patients 30 and 43 had increased fitness and higher virulence by all three measures (<xref ref-type="fig" rid="fig6">Figure 6A,B</xref>, <xref ref-type="fig" rid="fig7">Figure 7</xref>); whereas, a decrease in fitness in isolate 13 of patient 1 was accompanied by a decrease in virulence (<xref ref-type="fig" rid="fig6">Figures 6C</xref> and <xref ref-type="fig" rid="fig7">Figure 7</xref>). A notable exception was patient 59, where fitness in vitro decreased while virulence phenotypes increased in a later isolate (<xref ref-type="fig" rid="fig6">Figures 6D</xref> and <xref ref-type="fig" rid="fig7">Figure 7</xref>). This is consistent with the observations of Noble and co-workers that in vitro fitness is not always a reflection of virulence (<xref ref-type="bibr" rid="bib54">Noble et al., 2010</xref>).</p><p>Initially in a series, drug resistance (MIC) and in vitro fitness (in the absence of drug) were inversely related, suggesting that these are competing selective pressures. When MIC increases first appeared, they were usually accompanied by a <italic>decrease</italic> in fitness in the absence of drug (Patient 1, isolates 2, 13, and 16, Patients 9, 14, 15 16, and 59, <xref ref-type="fig" rid="fig7">Figure 7A</xref>). Consistent with the elevated MIC, these isolates exhibited <italic>increased</italic> relative fitness in the presence of the drug (<xref ref-type="fig" rid="fig7">Figure 7B</xref>). This is also consistent with a recent study (<xref ref-type="bibr" rid="bib70">Sasse et al., 2012</xref>) showing that resistance conferred in <italic>C. albicans</italic> by gain-of-function mutations in the transcription factors Mrr1, Tac1, and Upc2 is associated with reduced fitness under non-selective conditions in vitro as well as in vivo during colonization of a mammalian host. Consistent with subsequent selection of strains with compensatory variations, isolates from later time points were often more fit than those from earlier time points (measured in vitro<italic>,</italic> in the absence of drug) without further changes in MIC (e.g<italic>.</italic>, patient 1, isolates 5-7, isolate 14, <xref ref-type="fig" rid="fig7">Figure 7A</xref>), with notable exceptions (e.g., isolates 8 and 11). This general trend is consistent with previous studies in bacteria (<xref ref-type="bibr" rid="bib5">Bjorkman and Andersson, 2000</xref>); (<xref ref-type="bibr" rid="bib25">Gagneux et al., 2006</xref>) and in a single documented case in <italic>C. glabrata</italic> (<xref ref-type="bibr" rid="bib77">Singh-Babak et al., 2012</xref>), suggesting that compensatory mutations may subsequently arise to offset the major fitness cost of mutations conferring drug resistance. Nevertheless, substantial additional sampling will be required per time point to fully interpret such patterns.</p><p>In this context, it appears that aneuploidies (<xref ref-type="fig" rid="fig7">Figure 7A</xref>, green), while largely transient, may be an important intermediate giving rise to more stable adaptive genotypes in some cases, as was recently demonstrated in budding yeast adapting to a stressful environment in vitro (<xref ref-type="bibr" rid="bib89">Yona et al., 2012</xref>). For example, in Patient 1 isolate 13, an increase in MIC and a trisomy of 5 of 8 chromosomes accompanies a large decrease in fitness (in the absence of drug) relative to the preceding isolate 12 (<xref ref-type="fig" rid="fig7">Figure 7A</xref>) but has increased fitness in the presence of drug (<xref ref-type="fig" rid="fig7">Figure 7B</xref>). Isolate 14 has a similar MIC phenotype to isolate 13 but is euploid (<xref ref-type="fig" rid="fig3">Figure 3</xref>) and is much more fit (<xref ref-type="fig" rid="fig7">Figure 7A</xref>). Consistent with the general negative effect of aneuploidy on fitness (<xref ref-type="bibr" rid="bib80">Tang and Amon, 2013</xref>), the absence of the extra chromosomes resulted in improved overall fitness.</p></sec><sec id="s2-8"><title>Candidate mutated genes associated with drug resistance or virulence</title><p>The analysis of clinical isolates identified a range of new candidate genes that may affect drug resistance, fitness, and/or virulence. To test the contribution of some of the recurrently identified genes to specific <italic>C. albicans</italic> phenotypes, we profiled all 23 recurrently mutated loci for which a homozygous deletion mutant was available from a deletion strain collection (<xref ref-type="bibr" rid="bib54">Noble et al., 2010</xref>). We measured the MIC of fluconazole and the in vitro fitness in the absence of drug for each of these 23 mutants.</p><p>Deletion of one gene (<italic>orf19.4658</italic>) caused a twofold decrease in MIC, whereas the other 22 mutants tested did have no significant effect on MIC (<italic>data not shown</italic>). Deletion mutants are loss-of-function mutations, whereas the previously identified mechanisms of fluconazole resistance are ‘gain of function’, resulting in the increase in the amount or activity level of Erg11 (<xref ref-type="bibr" rid="bib3">Asai et al., 1999</xref>; <xref ref-type="bibr" rid="bib55">Oliver et al., 2007</xref>) or the efflux of drug transporters (<xref ref-type="bibr" rid="bib13">Coste et al., 2006</xref>; <xref ref-type="bibr" rid="bib18">Dunkel et al., 2008</xref>). Therefore, it is possible that the recurrent non-synonymous coding SNPs in the new loci, we identified in the clinical isolates confer resistance. Alternatively, these loci may not be involved in fluconazole resistance per se and may have a more general role in adaptation to the complex host environment.</p><p>Consistent with a role in host adaptation, 5 of the 22 deletion mutants reduced in vitro fitness in a culture medium thought to approximate in vivo conditions (<xref ref-type="fig" rid="fig8">Figure 8</xref>, ‘Materials and methods’). Three were significantly more fit than the WT parental strain (SN250, red, <xref ref-type="fig" rid="fig8">Figure 8</xref>), including <italic>CCN1</italic>, that encodes a G1 cyclin required for hyphal growth maintenance (<xref ref-type="bibr" rid="bib47">Loeb et al., 1999</xref>) and <italic>orf19.4471</italic>, an ortholog of <italic>Saccharomyces cerevisiae VPS64</italic>, which is required for cytoplasm-to-vacuole targeting of proteins (<xref ref-type="bibr" rid="bib6">Bonangelino et al., 2002</xref>), is involved in recycling pheromone receptors (<xref ref-type="bibr" rid="bib40">Kemp and Sprague, 2003</xref>), and is identified as an ‘aneuploidy-tolerating mutant’ (<xref ref-type="bibr" rid="bib82">Torres et al., 2010</xref>). Among the least fit were cell wall protein genes (<italic>HYR1</italic>, <italic>HYR3,</italic> and <italic>PIR1</italic>; <xref ref-type="bibr" rid="bib16">De Groot et al., 2003</xref>).<fig id="fig8" position="float"><object-id pub-id-type="doi">10.7554/eLife.00662.021</object-id><label>Figure 8.</label><caption><title>Deletion mutants of recurrently mutated genes reveal changes in relative fitness.</title><p>Shown is the fitness (‘Materials and methods’) for each deletion mutant strain and the corresponding wild-type strain (Y axis, mean ± STDV). The wild-type parental strain (SN250) is on the far left (red bar and dashed line). Fitness is calculated relative to an ENO1::YFP SC5314 reference isolate. Locus names are given for the mutant isolates (X axis). Asterisks denote statistical significance (* < 0.05, ** < 0.01, *** < 0.001, **** < 0.0001) by one-way ANOVA with Holm–Sidak correction for multiple comparisons.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00662.021">http://dx.doi.org/10.7554/eLife.00662.021</ext-link></p></caption><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="elife00662f008"/></fig></p></sec></sec><sec sec-type="discussion" id="s3"><title>Discussion</title><p>We sequenced the genomes of serial clinical isolates of <italic>C. albicans</italic> and analyzed them by comparing consecutive isolates from one patient to reach novel insights into drug resistance within the human host. This approach allowed us to distinguish (and remove from further analysis) isolates that were non-clonal and to estimate that at least ∼30% of the patients (3/10) carried at least one non-clonal strain of <italic>C. albicans.</italic> We used the clonal isolates to identify persistent SNPs, and the different series to identify those persistent SNPs that recurred within the same ORFs, thereby focusing the analysis on a small number of loci where the identified variants are more likely to be adaptive and excluding the substantial background of likely neutral variation that hitchhike along with selective beneficial mutations.</p><p>Our study identified substantial genetic diversity in each series, in contrast to the report of only 26 SNPs detected in a single clinical series of <italic>Candida glabrata</italic> isolates that spanned a 10-month period (<xref ref-type="bibr" rid="bib77">Singh-Babak et al., 2012</xref>). Several reasons may account for this difference. First, fluconazole, the antifungal drug used to treat the patients in our series, is fungistatic, such that many cells exposed to the drug arrest their growth but do not die. Thus, the range of diversity in the initial population is not entirely lost. In contrast, the <italic>C. glabrata</italic> study involved exposure to caspofungin, an echinocandin fungicidal drug. Therefore, most cells likely died upon drug exposure and only the rare survivors went on to seed the remaining population. Accordingly, the <italic>C. glabrata</italic> isolates may have been subjected to selection that would have removed much of the initial diversity in the population, whereas in the <italic>C. albicans</italic> series diversity persisted and selection acted mostly to change the relative proportions of different genotypes. Second, <italic>C. albicans</italic> is a highly heterozygous diploid whereas <italic>C. glabrata</italic> is haploid. Mutations can be more readily assimilated in a diploid than in a haploid organism, since deleterious mutations are potentially buffered by a functional version (<xref ref-type="bibr" rid="bib81">Thompson et al., 2006</xref>). Furthermore, because <italic>C. albicans</italic> genomes are initially much more diverse (with tens of thousands of heterozygous SNPs in a given isolate), LOH is a high frequency mechanism available to reveal mutations more readily. Third, <italic>C. albicans</italic> lab isolates likely undergo a stress-induced elevation of mutation and mitotic recombination rates (<xref ref-type="bibr" rid="bib62">Ponder et al., 2005</xref>; <xref ref-type="bibr" rid="bib20">Forche et al., 2011</xref>; <xref ref-type="bibr" rid="bib68">Rosenberg, 2011</xref>), and exposure to a mammalian host results in elevated frequencies of LOH and aneuploidy (<xref ref-type="bibr" rid="bib23">Forche et al. 2009</xref>). Thus, it is possible that <italic>C. albicans</italic> isolates within the human host also undergo elevated levels of LOH and of point mutations to generate a wider range of diversity. Thus, <italic>C. albicans</italic> like <italic>S. cerevisiae</italic> (<xref ref-type="bibr" rid="bib31">Gresham et al., 2008</xref>; <xref ref-type="bibr" rid="bib57">Pavelka et al., 2010</xref>; <xref ref-type="bibr" rid="bib89">Yona et al., 2012</xref>; <xref ref-type="bibr" rid="bib10">Chang et al., 2013</xref>) generates large scale genetic variation as a means of adaptation. This adds another level of variation to the genome and protein diversity (<xref ref-type="bibr" rid="bib69">Santos et al., 2004</xref>; <xref ref-type="bibr" rid="bib73">Selmecki et al., 2010</xref>) that <italic>C. albicans</italic> is able to tolerate.</p><p>LOH in several genes important for fluconazole resistance has been reported previously for <italic>ERG11</italic> (<xref ref-type="bibr" rid="bib55">Oliver et al., 2007</xref>), <italic>TAC1</italic> (<xref ref-type="bibr" rid="bib13">Coste et al., 2006</xref>), and <italic>MRR1</italic> (<xref ref-type="bibr" rid="bib71">Schubert et al., 2011</xref>), but the degree to which LOH is important in clinical infections was not known. In the ten patients studied here, LOH was commonly observed and was associated with changes in MIC. As we detected LOH of mutations in <italic>ERG11</italic> in three patients, it would be of interest to know if the LOH in these known genes was sufficient to increase MIC, or if other genes within the homozygous region make important contributions.</p><p>Aneuploidies appeared frequently within the drug-resistant isolates, consistent with previous reports (<xref ref-type="bibr" rid="bib72">Selmecki et al., 2006</xref>). Unlike LOH events, aneuploidies were often transient and not consistently correlated to increases in drug resistance. Perhaps these aneuploidies provide a mechanism akin to genetic assimilation (‘phenotype precedes genotype’), in which cells are provided with a phenotypic mechanism that facilitates survival until a more stable and/or less costly mechanism is attained. In this case, the ‘phenotypic’ mechanism would be genetic but unstable—the acquisition of one or more extra chromosomes. Nevertheless, aneuploidies may cause increased frequencies of LOH events through whole chromosome loss, as well as by increasing the likelihood of recombination events. A transient role for aneuploidy is consistent with recent findings from in vitro evolution studies in <italic>S. cerevisiae</italic> (<xref ref-type="bibr" rid="bib89">Yona et al., 2012</xref>) in which a transient aneuploidy was responsible for fitness at elevated temperature, but was eventually replaced by a more stable mutation.</p><p>In addition, a substantial number of persistent and recurrent SNPs, and clusters of co-occurring SNPs, implicate a broad range of pathways and functions that likely provide some growth advantage in the presence of the complex selective pressures found in the host. In particular, there was strong enrichment for cell wall gene families thought to be critical determinants of the transition from commensalism to pathogenesis (<xref ref-type="bibr" rid="bib28">Gow and Hube, 2012</xref>). The genes in several of these families (e.g., <italic>ALS1</italic>-<italic>4</italic>,<italic>7</italic>,<italic>9</italic> and <italic>HYR</italic>/<italic>IFF</italic> genes) frequently contain intragenic tandem repeats. Variation in intragenic repeat number modulates phenotypic diversity in adhesion and biofilm formation (<xref ref-type="bibr" rid="bib84">Verstrepen et al., 2005</xref>). This functional diversity of cell surface antigens has been proposed to allow rapid adaptation to the environment as well as evasion of the host immune system in fungi and other pathogens (<xref ref-type="bibr" rid="bib27">Gemayel et al., 2010</xref>). Notably, the cell wall deletion mutants (<italic>HYR1</italic>, <italic>HYR3,</italic> and <italic>PIR1</italic>) were among the least fit in vitro (<xref ref-type="fig" rid="fig8">Figure 8</xref>).</p><p>Indeed, many of the isolates evolved additional phenotypes, including changes in in vitro fitness, filamentation, adhesion, and in vivo virulence, and the data presented here points to candidate genes that underlie some of these evolved traits. For example, the evolved isolate 4617 in patient 59 had a dramatic increase in filamentation relative to the progenitor, which was concomitant with the appearance of persistent SNPs in genes associated with filamentous growth: <italic>CHO1</italic>, <italic>MNN2,</italic> and 7 different <italic>FGR</italic> (filamentous growth regulator) genes (<xref ref-type="bibr" rid="bib83">Uhl et al., 2003</xref>).</p><p>The evolution of drug resistance in <italic>C. albicans</italic> has many parallels with the somatic evolution of cancer cells undergoing chemotherapy or treated with specific inhibitors. These include variation on a background of clonal descent, lack of sexual recombination, acquisition of drug resistance, tolerance of aneuploidy and genome plasticity, and increased mutation and mitotic recombination rates under stress. Indeed, several recent studies have shown a similar spectrum of genetic alterations to those observed here during the somatic evolution of cancers in patients undergoing chemotherapy (<xref ref-type="bibr" rid="bib61">Podlaha et al., 2012</xref>; <xref ref-type="bibr" rid="bib44">Landau et al., 2013</xref>) or treated with specific inhibitors (<xref ref-type="bibr" rid="bib17">Ding et al., 2012</xref>) to those observed here.</p><p>Finally, our data and analyses provide a rich and novel resource for <italic>Candida</italic> researchers and a host of candidate genes for further functional studies. While our analysis focused on recurrent SNPs in ORFs, we nonetheless cataloged the many genetic alterations found in intergenic regions (<xref ref-type="supplementary-material" rid="SD1-data">Figure 2—source data 1</xref>), some of which could affect gene regulation. It will be especially interesting to analyze the similarities and differences in additional <italic>C. albicans</italic> genome sequences that are likely to become available in the near future. Our results suggest there may be complex population dynamics during the transition from commensal to pathogen and across the course of treatment. As sequencing capacity continues to grow, it will be especially interesting to more fully sample this population-level diversity during longitudinal collection to better understand these dynamics. In particular, it will be interesting to determine the degree to which specific mutations recur in different isolates prior to and after the acquisition of drug resistance.</p></sec><sec sec-type="materials|methods" id="s4"><title>Materials and methods</title><sec id="s4-1"><title>Isolates</title><p>Isolates were obtained from HIV-infected patients with oropharyngeal candidiasis, as previously described (<xref ref-type="bibr" rid="bib85">White, 1997a</xref>; <xref ref-type="bibr" rid="bib58">Perea et al., 2001</xref>). The patients were not on azole antifungal treatment at time of enrollment; subsequent samples were collected during recurrence of infection. Isolates were colony purified at collection and represent a single clone. The isolates are detailed in <xref ref-type="table" rid="tbl1">Table 1</xref>.</p></sec><sec id="s4-2"><title>Drug susceptibility</title><p>Minimal inhibitory concentrations (MIC) were determined for each strain (clinical and mutant) using fluconazole E-test strips (0.016–256 μg/ml, bioMérieux, Durham, NC) on RPMI 1640-agar plates (Remel, Lenexa, KS). Overnight YPD cultures were diluted in sterile 0.85% NaCl to an OD600 of 0.01 and 250 μl was plated using beads. After a 30-min pre-incubation, 2–3 E-test strips were applied and plates were incubated at 35°C for 48 hr. The susceptibility endpoint was read at the first growth-inhibition ellipse, and the median value is reported here.</p></sec><sec id="s4-3"><title>Illumina sequencing</title><p>Genomic DNA was prepared from different clinical time courses via a Qiagen Maxiprep kit (Qiagen, Valencia, CA) and sequenced using 101 base paired-end Illumina sequencing (<xref ref-type="bibr" rid="bib49">Mardis, 2008</xref>). Library preparation included an eight base barcode (<xref ref-type="bibr" rid="bib30">Grabherr et al., 2011</xref>); 43 samples from 11 patients were sequenced. All reads were mapped to the SC5314 reference genome (Candida Genome Database Assembly 21, gff downloaded on 4 January 2010) using the BWA alignment tool (version 0.5.9) (<xref ref-type="bibr" rid="bib45">Li and Durbin, 2009</xref>). To minimize false positive SNP calls near insertion/deletion (indel) events, poorly aligning regions were identified and realigned using the GATK RealignerTargetCreator and IndelRealigner (GATK version 1.4-14, [version 1.4-14]) (<xref ref-type="bibr" rid="bib51">McKenna et al., 2010</xref>). Coverage for each strain is reported in <xref ref-type="table" rid="tbl1">Table 1</xref>. Coverage was defined as the total number of bases with BWA mapping quality greater than 10 divided by the total number of sites in the nuclear genome. Isolate 4380 aligned poorly to the SC5314 genome; however, this sequence aligned at high identity to the <italic>C. dubliniensis</italic> genome and was therefore removed from further analysis. These data can be accessed from a genomics portal hosted by the Broad Institute at: <ext-link xmlns:xlink="http://www.w3.org/1999/xlink" ext-link-type="uri" xlink:href="http://www.broadinstitute.org/pubs/candidadrugresistance/">http://www.broadinstitute.org/pubs/candidadrugresistance/</ext-link>. Reads are deposited for access to the NCBI SRA under project accession number PRJNA257929.</p></sec><sec id="s4-4"><title>SNP calling</title><p>SNPs were identified using Unified Genotyper (GATK version 1.4.14) (<xref ref-type="bibr" rid="bib51">McKenna et al., 2010</xref>), using read alignments to the SC5314 reference sequence. Unreliable SNPs were identified using the GATK Variant Filtration module, with the version 3 best practice recommended annotation filters (QD < 2.0, MQ < 40.0, FS > 60.0, HaplotypeScore >13.0, MQRankSum < −12.5, ReadPosRankSum < −8.0) except that the HaplotypeScore was also filtered if greater than two standard deviations above the mean of all HaplotypeScore values. The combined list of SNP positions across all strains was used to evaluate those matching the reference allele; by emitting all sites using Unified Genotyper, high quality reference matches were identified as positions with quality of 30 or greater, with positions with extremes of read depth (top or bottom 0.5% quantile) eliminated. A matrix of all strains by all positions was created from the SNP calls, with reference calls added where identified. Non-clonal isolates (see below) were removed from further analysis.</p></sec><sec id="s4-5"><title>Sequenom iPLEX genotyping assay</title><p>We chose 1973 genetic locus X strain combination (523 unique sites across nine patients) for iPLEX genotyping as either (1) persistent within their time course (605 sites), (2) background mutations (1263 sites), or (3) transient mutations (105 sites). All selected loci were at least 150 bp away from any other SNP in either direction to avoid ambigious iPlex calls. Sites producing multiple iPlex results were eliminated from further consideration. 1,853 predictions were confirmed as correct by Sequenom genotyping and 120 were discordant (<xref ref-type="table" rid="tbl2">Table 2</xref>), to a calling accuracy of 93.9%.</p></sec><sec id="s4-6"><title>Determination of relatedness to determine clonality</title><p>We investigated the phylogenetic relationship of all strains using SNP calls to determine relatedness between strains; positions with missing data in 10% or more of strains were eliminated, resulting in a total of 201,793 parsimony informative positions. A distance based tree was estimated using maximum parsimony with PAUP* (4.0) (<xref ref-type="bibr" rid="bib79">Swofford, 2002</xref>); a step matrix was used to estimate the distance between homozygous and heterozygous positions, where each of the homozygotes is two steps apart from each other and one step from the heterozygote. SNP positions were resampled using 1000 bootstrap replicates, and the phylogeny re-estimated to test the branch support. We define isolates with a branch distance of greater than 20,000 as non-clonal.</p></sec><sec id="s4-7"><title>Copy-number determination</title><p>For each strain, we calculated a per-locus depth-of-coverage using GATK (<xref ref-type="bibr" rid="bib51">McKenna et al., 2010</xref>), with a minimal mapping quality of 10. The number of reads aligning to each 5 kb window across the nuclear genome was calculated and then normalized to the genome median. Each bin was then multiplied to the ploidy for the majority of the genome as determined by a FACS assay (below). We then applied a sliding window across each bin, defining a potential CNV if 70% of 10 consecutive bins had a normalized count >2.5×. Regional amplifications are identified if >15% of the chromosome is identified as having a CNV. Boundaries were confirmed by visual inspection in the Integrative Genome Viewer (<xref ref-type="bibr" rid="bib67">Robinson et al., 2011</xref>).</p></sec><sec id="s4-8"><title>High-resolution ploidy analysis by flow cytometry</title><p><italic>C. albicans</italic> cultures were grown to log phase. 200 μl of culture was centrifuged in a round bottom microtiter plate, and pellets were resuspended in 20 μl of 50 mM Tris pH8/50 mM EDTA (50/50 TE). 180 μl of 95% ethanol was added and suspensions were stored overnight at 4°C. Cells were centrifuged and pellets washed twice with 200 μl of 50/50 TE, then resuspended in 50 μl of RNAse A at 1 mg/ml in 50/50 TE and incubated 1 hr at 37°C. Cells were centrifuged and pellets resuspended in 50 μl of Proteinase K at 5 mg/ml in 50/50 TE for 30 min at 37°C. Cells were washed in 50/50 TE and pellets resuspended in 50 μl of a 1:85 dilution SYBR Green I (Invitrogen, Carlsbad, CA) in 50/50 TE and incubated overnight in the dark at 4°C. Cells were centrifuged and pellets were resuspended in 700 μl 50/50 TE and read on a FACS caliber flow cytometer (BD Biosciences, San Jose, CA). Flow data were fitted with a multi-Gaussian cell cycle model to produce estimates for whole genome ploidy.</p></sec><sec id="s4-9"><title>LOH determination</title><p>For each time course, we assembled the high quality SNPs (post-filtering, above) from multi-sample calling into the columns of a matrix, ordered by genome position, with the isolates in rows, ordered temporally. The genetic state of each locus in each sample was coded to distinguish loci homozygous for the haploid reference (−1), heterozygous SNPs (0), and homozygous SNPs for the non-reference state (1). We then applied a sliding window method across each chromosome, only looking at sites in which a SNP call was made in at least one isolate. An LOH event was defined as occurring if (1) at least one isolate had a heterozygosity content >40%, and (2) at least one other isolate had a heterozygosity content <5%. Window sizes were of length 500. Boundaries were trimmed such that if a window terminated in a heterozygous site in the isolate for which the LOH occurred, it was trimmed back until it was homozygous. If two 500+ windows were within 7 KB of each, the region was assessed to determine if the event was actually one event and merged if the heterozygous sites in the inter-window space had homozygosed. If two isolates had LOHs that overlapped but did not have precisely identical boundaries, the LOH regions were combined such that the LOH interval for both isolates was the same. All LOH regions were confirmed by visual inspection and are listed in <xref ref-type="fig" rid="fig3">Figure 3</xref> and <xref ref-type="supplementary-material" rid="SD3-data">Figure 4—source data 1</xref>.</p></sec><sec id="s4-10"><title>Classification of SNPs</title><p>For each time course, each SNP was classified for its position in the genome (<xref ref-type="supplementary-material" rid="SD1-data">Figure 2—source data 1</xref>). If the SNP fell within an ORF, the reference and altered amino acids were reported. If the SNP fell outside of an ORF, the distance to the closest flanking ORF(s) was reported, as well as the SNP's orientation with respect to these ORFs. SNP genotypes that are common to all isolates (including the ‘progenitor’) were classified as background mutations. Genotypes not present in the progenitor or evolved strain, but that occur in one or more intermediate strain, are classified as transient. Finally, genotypes that occur after the progenitor, and persist through the terminally evolved time point, are classified as persistent.</p><p>To determine if the number of persistent non-synonymous SNPs (nsSNPs) occurring in conjunction with changes in MIC was greater than expected, we developed a simple model to simulate the occurrence of nsSNPs outside of LOH regions at each time point. For each time point (i), a random variable <italic>X</italic><sub>i</sub> <italic>∼</italic> Pois(<italic>λ</italic><sub>i</sub>) was assigned, where <italic>λ</italic><sub>i</sub> represents the Poisson parameter for each time point:<list list-type="simple"><list-item><p>λ<sub>i</sub> = m/T × t<sub>i</sub>;</p></list-item><list-item><p>m = 1471, the number of ORFs with persistent nsSNP;</p></list-item><list-item><p>T = 23, the number of time points;</p></list-item><list-item><p>t<sub>i</sub> = the length of time (days) for time point (i) divided by the mean length of time.</p></list-item></list></p><p>The number of persistent nsSNP-containing ORFs for each of the 14 time points associated with a change in MIC was summed, and this was repeated 100,000 times to build a probability distribution, where <italic>p</italic> (observing <italic>x</italic> mutated ORFs) was determined by dividing the number of successes for each bin by the number of trials.</p><p>To determine the probability of observing <italic>x</italic> recurrent, persistently mutated ORFs outside of LOH regions, we developed an additional stochastic simulation model. For each patient series (i), at each non-LOH ORF (j), a random variable, X<sub>ij</sub> ∼ B(<italic>n, p</italic><sub>ij</sub>) was assigned, where <italic>n</italic> represents the number of trials, 1, and <italic>p</italic> represents the probability of a SNP occurring in that ORF:</p><p>p = m<sub>i</sub>/M<sub>i</sub> × h<sub>j</sub>;</p><p>m<sub>i</sub> = number of ORFs with persistent nsSNPs found outside of LOH events for patient series (i);</p><p>M<sub>i</sub> = number of ORFs outside of LOH events for patient series (i);</p><p>h<sub>j</sub> = log normalized ORF length divided by mean lognormalized ORF length for ORF (j).</p><p>This was repeated 10,000 times to build a distribution, and p (observing <italic>x</italic> recurrent nsSNPs) was determined by dividing the number of successes in each bin by the cumulative number of trials.</p></sec><sec id="s4-11"><title>Analysis of co-occurring mutations</title><p>For co-occurrence analysis we focused only SNPs that (1) had persistent nonsynonymous coding SNPs that did not occur in LOH regions and (2) recurred in three or more time courses. We generated for each such gene a binary patient vector, and we used NMF clustering (<xref ref-type="bibr" rid="bib7">Brunet et al., 2004</xref>) to identify the optimal number of clusters, based on local maximas. This was accomplished using the ‘NMFConsensus’ module (version 5) in GenePattern (<xref ref-type="bibr" rid="bib63">Reich et al., 2006</xref>). To determine the most appropriate number of clusters, k was selected such that it was the smallest value for which the cophenetic correlation begins decreasing. We then tested each of the co-occurrence gene clusters for functional enrichment (below). To determine if the degree of co-occurrence would have arisen by chance, we ran 1000 iterations of 1 million edge-pair swaps from the original binary matrix, calculating a Pearson correlation matrix for each of the 1000 iterations. We compared the distribution of Pearson correlations on the real and permuted vectors using a two-sample Kolmogorov–Smirnov (KS) test and Wilcoxon Rank Sum test.</p></sec><sec id="s4-12"><title>Functional enrichment</title><p>We calculated the overlap of each co-occurring cluster with Gene Ontology gene sets using the Gene Ontology toolset from the Candida Genome Database (<xref ref-type="bibr" rid="bib2">Arnaud et al., 2009</xref>; <xref ref-type="bibr" rid="bib36">Inglis et al., 2012</xref>, <xref ref-type="bibr" rid="bib37">2013</xref>). Bonferroni adjusted p-values as well as the False Discovery Rate are reported (<xref ref-type="supplementary-material" rid="SD4-data">Figure 5—source data 1A</xref>).</p></sec><sec id="s4-13"><title>Competition assay to assess fitness</title><p>We measured the relative fitness of the progenitor and evolved lines in RPMI Cell Culture medium (Gibco, Grand Island, NY), competing them against a reference strain (SC5314), expressing <italic>ENO1::YFP</italic>. Isolates stored at −80°C were revived on rich media petri plates and then grown overnight in 3-ml cultures of minimal media in a roller drum at 35°C. An aliquot of cells in each culture was removed, sonicated in a Branson 450 sonifier, and the concentration of cells was determined using a Cellometer M10 (Nexcelom, Lawrence, MA). The reference strain and experimental competitors were added to fresh RPMI medium in a 1:1 ratio and a final cell concentration of 1 × 10<sup>7</sup> cells/ml. The cultures were grown for 24 hr in a roller drum at 35°C. Cells were then counted as above, and 3 × 10<sup>6</sup> cells were transferred to fresh RPMI medium grown for 24 hr in a roller drum at 35°C (transfer cycle 1). This procedure was repeated (transfer cycle 2). This protocol represents 5–10 generations of growth, depending on the strain genotype. The ratio of the two competitors was quantified at the initial and final time points by flow cytometry (Accuri, San Jose, CA). 3 to 6 independent replicates for each fitness measurement were performed. The selective advantage, <italic>s</italic>, or disadvantage of the evolved population was calculated as previously described (<xref ref-type="bibr" rid="bib81">Thompson et al., 2006</xref>), where <italic>E</italic> and <italic>R</italic> are the numbers of evolved and reference cells in the final (<italic>f</italic>) and initial (<italic>i</italic>) populations, and <italic>T</italic> is the number of generations that reference cells have proliferated during the competition.</p></sec><sec id="s4-14"><title>Competition assay to assess fitness with and without fluconazole</title><p>Fitness assays were performed as described above except that the reference strain used was a derivative of the drug-resistant isolate 4639 from patient 59 (<xref ref-type="table" rid="tbl1">Table 1</xref>) expressing <italic>ENO1::YFP</italic> and competition experiments were performed in RPMI and in RPMI with one half the MIC for fluconazole (<xref ref-type="table" rid="tbl1">Table 1</xref>) initiated from replicate 1:1 mixtures of the same population of cells for each isolate tested.</p><p>In order to quantify fitness in the presence of fluconazole, we constructed a derivative of the fluconazole-resistant (128 μg/ml) isolate 3795 from patient 9 (<xref ref-type="table" rid="tbl1">Table 1</xref>) that expresses <italic>ENO1</italic>::<italic>YFP</italic> to use as the reference for competitive fitness assays in the presence and absence of fluconazole. This strain was chosen since it was a euploid strain with the highest MIC. There were two strains from patient 42 (3731 and 3733) with a higher MIC but these strains are aneuploid and thus the potential loss of additional chromosomes during the course of the competition could alter fitness and confound our results. The addition of the YFP marker reduced the fitness of the strain relative to the unaltered one and slightly reduced the fluconazole MIC as measured by the E-strip test. Fitness assays were performed as described above competition experiments except they were performed in RPMI and in RPMI with one half the MIC for fluconazole (<xref ref-type="table" rid="tbl1">Table 1</xref>) initiated from replicate 1:1 mixtures of the same population of cells for each isolate tested.</p></sec><sec id="s4-15"><title><italic>C. elegans</italic> survival assay</title><p>A <italic>C. elegans</italic> survival assay was performed as previously described (<xref ref-type="bibr" rid="bib39">Jain et al., 2009</xref>). Briefly, <italic>Escherichia coli</italic> OP50 and the different <italic>C. albicans</italic> clinical isolates were grown overnight respectively in LB at 37°C and YPD at 30°C. <italic>E. coli</italic> was then centrifuged and resuspended to a final concentration of 200 mg/ml, while <italic>C. albicans</italic> isolates were diluted with sterile water to OD<sub>600</sub> = 3. Small petri dishes (3.5 cm) containing NGM agar were spotted with a mixture of 10 μl streptomycin (stock solution 50 mg/ml), 2.5 μl of <italic>E. coli</italic>, 0.5 μl of <italic>C. albicans</italic>, and 7 μl of sterile water. The plates were incubated overnight at 25°C and 20 young synchronized N2 <italic>C. elegans</italic> adults were transferred on the spotted plates. Synchronous populations of adult worms were obtained by plating eggs on NGM plates seeded with <italic>E. coli</italic> OP50 at 20°C for 2–3 days. In this time frame, the eggs hatch and the larvae reach young adulthood. The survival assay was carried at 20°C, and worms were scored daily by gentle prodding with a platinum wire; dead worms were discarded while live ones were transferred to seeded plates grown overnight at 25°C. Worms accidentally killed while transferring or found dead on the edges of the plates were censored. Statistical analysis was performed using SPSS software; survival curves were obtained using the Kaplan–Meier method and p-values by using the log-rank test.</p></sec><sec id="s4-16"><title>Filamentation assay</title><p>Overnight cultures grown in YPD at 30°C were normalized to OD<sub>600</sub> = 1 with sterile water and spotted on Spider agar media (1% mannitol, 1% Difco nutrient broth, 0.2% K<sub>2</sub>HPO<sub>4</sub>). Plates were incubated at 37°C and colonies were photographed 3, 7, and 10 days post spotting. As a negative control for filamentation <italic>cph1/cph1 efg1/efg1</italic> (<xref ref-type="bibr" rid="bib46">Lo et al., 1997</xref>) double mutant strain was used.</p></sec><sec id="s4-17"><title>In vitro adhesion assay</title><p>The in vitro adhesion assay was performed as previously described for <italic>S. cerevisiae</italic> (<xref ref-type="bibr" rid="bib66">Reynolds and Fink, 2001</xref>). Briefly, cultures were grown in Synthetic Complete (SC) media + 0.15% glucose at 30°C overnight. Cells were then centrifuged at maximal speed and resuspended to OD<sub>600</sub> = 0.5 in fresh media. 200 ml of each culture were dispensed into 8 wells of a flat bottom 96-well plate and incubated at 37°C for 4 hr. The content of the plate was then decanted and 50 ml of crystal violet added to each well. After 45 min of incubation at room temperature, the content of the plate was decanted and the plate was rinsed ten times in DI water by alternate submerging and decanting. 200 ml of 75% methanol was added to each well and absorbance was measured after 30 min at OD<sub>590</sub>. An <italic>edt1/edt1</italic> knockout mutant (<xref ref-type="bibr" rid="bib90">Zakikhany et al., 2007</xref>) was used as a negative control for adhesion.</p></sec></sec></body><back><sec sec-type="additional-information"><title>Additional information</title><fn-group content-type="competing-interest"><title>Competing interests</title><fn fn-type="conflict" id="conf1"><p>AR: Senior editor, <italic>eLife</italic>.</p></fn><fn fn-type="conflict" id="conf2"><p>The other authors declare that no competing interests exist.</p></fn></fn-group><fn-group content-type="author-contribution"><title>Author contributions</title><fn fn-type="con" id="con1"><p>CBF, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con2"><p>JMF, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con3"><p>DAT, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con4"><p>DA, Acquisition of data, Analysis and interpretation of data</p></fn><fn fn-type="con" id="con5"><p>LI, Acquisition of data, Analysis and interpretation of data</p></fn><fn fn-type="con" id="con6"><p>CG, Conception and design, Acquisition of data, Analysis and interpretation of data</p></fn><fn fn-type="con" id="con7"><p>CC, Conception and design, Acquisition of data, Analysis and interpretation of data</p></fn><fn fn-type="con" id="con8"><p>DAM, Conception and design, Analysis and interpretation of data, Contributed unpublished essential data or reagents</p></fn><fn fn-type="con" id="con9"><p>TD, Conception and design, Acquisition of data, Contributed unpublished essential data or reagents</p></fn><fn fn-type="con" id="con10"><p>BL, Conception and design, Acquisition of data, Contributed unpublished essential data or reagents</p></fn><fn fn-type="con" id="con11"><p>TCW, Provided many of the fungal samples used in the study, Provided essential domain expertise necessary to interpret results</p></fn><fn fn-type="con" id="con12"><p>RPR, Conception and design, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con13"><p>JB, Conception and design, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con14"><p>AR, Conception and design, Analysis and interpretation of data, Drafting or revising the article</p></fn></fn-group></sec><sec sec-type="supplementary-material"><title>Additional files</title><sec sec-type="datasets"><title>Major dataset</title><p>The following dataset was generated:</p><p><related-object content-type="generated-dataset" id="dataro1"><name><surname>Cuomo</surname><given-names>C</given-names></name> and <name><surname>Ford</surname><given-names>C</given-names></name>, <year>2014</year><x>, </x><source>The evolution of gradual acquisition of drug resistance in clinical isolates of <italic>Candida albicans</italic></source><x>, </x><object-id pub-id-type="art-access-id">PRJNA257929</object-id><x>; </x><ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/sra">http://www.ncbi.nlm.nih.gov/sra</ext-link><x>, </x><comment>Publicly available at NCBI Sequence Read Archive.</comment></related-object></p></sec></sec><ref-list><title>References</title><ref id="bib1"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Alby</surname><given-names>K</given-names></name><name><surname>Bennett</surname><given-names>RJ</given-names></name></person-group><year>2009</year><article-title>Stress-induced phenotypic switching in <italic>Candida albicans</italic></article-title><source>Molecular Biology of the Cell</source><volume>20</volume><fpage>3178</fpage><lpage>3191</lpage><pub-id pub-id-type="doi">10.1091/mbc.E09-01-0040</pub-id></element-citation></ref><ref id="bib2"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Arnaud</surname><given-names>MB</given-names></name><name><surname>Costanzo</surname><given-names>MC</given-names></name><name><surname>Shah</surname><given-names>P</given-names></name><name><surname>Skrzypek</surname><given-names>MS</given-names></name><name><surname>Sherlock</surname><given-names>G</given-names></name></person-group><year>2009</year><article-title>Gene Ontology and the annotation of pathogen genomes: the case of <italic>Candida albicans</italic></article-title><source>Trends in Microbiology</source><volume>17</volume><fpage>295</fpage><lpage>303</lpage><pub-id pub-id-type="doi">10.1016/j.tim.2009.04.007</pub-id></element-citation></ref><ref id="bib3"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Asai</surname><given-names>K</given-names></name><name><surname>Tsuchimori</surname><given-names>N</given-names></name><name><surname>Okonogi</surname><given-names>K</given-names></name><name><surname>Perfect</surname><given-names>JR</given-names></name><name><surname>Gotoh</surname><given-names>O</given-names></name><name><surname>Yoshida</surname><given-names>Y</given-names></name></person-group><year>1999</year><article-title>Formation of azole-resistant <italic>Candida albicans</italic> by mutation of sterol 14-demethylase P450</article-title><source>Antimicrobial Agents and Chemotherapy</source><volume>43</volume><fpage>1163</fpage><lpage>1169</lpage></element-citation></ref><ref id="bib4"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bennett</surname><given-names>RJ</given-names></name><name><surname>Johnson</surname><given-names>AD</given-names></name></person-group><year>2005</year><article-title>Mating in Candida albicans and the search for a sexual cycle</article-title><source>Annual Review of Microbiology</source><volume>59</volume><fpage>233</fpage><lpage>255</lpage><pub-id pub-id-type="doi">10.1146/annurev.micro.59.030804.121310</pub-id></element-citation></ref><ref id="bib5"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Björkman</surname><given-names>J</given-names></name><name><surname>Andersson</surname><given-names>DI</given-names></name></person-group><year>2000</year><article-title>The cost of antibiotic resistance from a bacterial perspective</article-title><source>Drug Resistance Updates</source><volume>3</volume><fpage>237</fpage><lpage>245</lpage><pub-id pub-id-type="doi">10.1054/drup.2000.0147</pub-id></element-citation></ref><ref id="bib6"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bonangelino</surname><given-names>CJ</given-names></name><name><surname>Chavez</surname><given-names>EM</given-names></name><name><surname>Bonifacino</surname><given-names>JS</given-names></name></person-group><year>2002</year><article-title>Genomic screen for vacuolar protein sorting genes in <italic>Saccharomyces cerevisiae</italic></article-title><source>Molecular Biology of the Cell</source><volume>13</volume><fpage>2486</fpage><lpage>2501</lpage><pub-id pub-id-type="doi">10.1091/mbc.02-01-0005</pub-id></element-citation></ref><ref id="bib7"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Brunet</surname><given-names>JP</given-names></name><name><surname>Tamayo</surname><given-names>P</given-names></name><name><surname>Golub</surname><given-names>TR</given-names></name><name><surname>Mesirov</surname><given-names>JP</given-names></name></person-group><year>2004</year><article-title>Metagenes and molecular pattern discovery using matrix factorization</article-title><source>Proceedings of the National Academy of Sciences of USA</source><volume>101</volume><fpage>4164</fpage><lpage>4169</lpage><pub-id pub-id-type="doi">10.1073/pnas.0308531101</pub-id></element-citation></ref><ref id="bib8"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Butler</surname><given-names>G</given-names></name><name><surname>Rasmussen</surname><given-names>MD</given-names></name><name><surname>Lin</surname><given-names>MF</given-names></name><name><surname>Santos</surname><given-names>MA</given-names></name><name><surname>Sakthikumar</surname><given-names>S</given-names></name><name><surname>Munro</surname><given-names>CA</given-names></name><name><surname>Grabherr</surname><given-names>M</given-names></name><name><surname>Forche</surname><given-names>A</given-names></name><name><surname>Reedy</surname><given-names>JL</given-names></name><name><surname>Agrafioti</surname><given-names>I</given-names></name><name><surname>Arnaud</surname><given-names>MB</given-names></name><name><surname>Bates</surname><given-names>S</given-names></name><name><surname>Brown</surname><given-names>AJ</given-names></name><name><surname>Brunke</surname><given-names>S</given-names></name><name><surname>Costanzo</surname><given-names>MC</given-names></name><name><surname>Fitzpatrick</surname><given-names>DA</given-names></name><name><surname>de Groot</surname><given-names>PW</given-names></name><name><surname>Harris</surname><given-names>D</given-names></name><name><surname>Hoyer</surname><given-names>LL</given-names></name><name><surname>Hube</surname><given-names>B</given-names></name><name><surname>Klis</surname><given-names>FM</given-names></name><name><surname>Kodira</surname><given-names>C</given-names></name><name><surname>Lennard</surname><given-names>N</given-names></name><name><surname>Logue</surname><given-names>ME</given-names></name><name><surname>Martin</surname><given-names>R</given-names></name><name><surname>Neiman</surname><given-names>AM</given-names></name><name><surname>Nikolaou</surname><given-names>E</given-names></name><name><surname>Quail</surname><given-names>MA</given-names></name><name><surname>Quinn</surname><given-names>J</given-names></name><name><surname>Santos</surname><given-names>MC</given-names></name><name><surname>Schmitzberger</surname><given-names>FF</given-names></name><name><surname>Sherlock</surname><given-names>G</given-names></name><name><surname>Shah</surname><given-names>P</given-names></name><name><surname>Silverstein</surname><given-names>KA</given-names></name><name><surname>Skrzypek</surname><given-names>MS</given-names></name><name><surname>Soll</surname><given-names>D</given-names></name><name><surname>Staggs</surname><given-names>R</given-names></name><name><surname>Stansfield</surname><given-names>I</given-names></name><name><surname>Stumpf</surname><given-names>MP</given-names></name><name><surname>Sudbery</surname><given-names>PE</given-names></name><name><surname>Srikantha</surname><given-names>T</given-names></name><name><surname>Zeng</surname><given-names>Q</given-names></name><name><surname>Berman</surname><given-names>J</given-names></name><name><surname>Berriman</surname><given-names>M</given-names></name><name><surname>Heitman</surname><given-names>J</given-names></name><name><surname>Gow</surname><given-names>NA</given-names></name><name><surname>Lorenz</surname><given-names>MC</given-names></name><name><surname>Birren</surname><given-names>BW</given-names></name><name><surname>Kellis</surname><given-names>M</given-names></name><name><surname>Cuomo</surname><given-names>CA</given-names></name></person-group><year>2009</year><article-title>Evolution of pathogenicity and sexual reproduction in eight Candida genomes</article-title><source>Nature</source><volume>459</volume><fpage>657</fpage><lpage>662</lpage><pub-id pub-id-type="doi">10.1038/nature08064</pub-id></element-citation></ref><ref id="bib9"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chandra</surname><given-names>J</given-names></name><name><surname>Kuhn</surname><given-names>DM</given-names></name><name><surname>Mukherjee</surname><given-names>PK</given-names></name><name><surname>Hoyer</surname><given-names>LL</given-names></name><name><surname>McCormick</surname><given-names>T</given-names></name><name><surname>Ghannoum</surname><given-names>MA</given-names></name></person-group><year>2001</year><article-title>Biofilm formation by the fungal pathogen Candida albicans: development, architecture, and drug resistance</article-title><source>Journal of Bacteriology</source><volume>183</volume><fpage>5385</fpage><lpage>5394</lpage><pub-id pub-id-type="doi">10.1128/JB.183.18.5385-5394.2001</pub-id></element-citation></ref><ref id="bib10"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chang</surname><given-names>SL</given-names></name><name><surname>Lai</surname><given-names>HY</given-names></name><name><surname>Tung</surname><given-names>SY</given-names></name><name><surname>Leu</surname><given-names>JY</given-names></name></person-group><year>2013</year><article-title>Dynamic large-scale chromosomal rearrangements fuel rapid adaptation in yeast populations</article-title><source>PLOS Genetics</source><volume>9</volume><fpage>e1003232</fpage><pub-id pub-id-type="doi">10.1371/journal.pgen.1003232</pub-id></element-citation></ref><ref id="bib11"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Charles</surname><given-names>PE</given-names></name><name><surname>Doise</surname><given-names>JM</given-names></name><name><surname>Quenot</surname><given-names>JP</given-names></name><name><surname>Aube</surname><given-names>H</given-names></name><name><surname>Dalle</surname><given-names>F</given-names></name><name><surname>Chavanet</surname><given-names>P</given-names></name><name><surname>Milesi</surname><given-names>N</given-names></name><name><surname>Aho</surname><given-names>LS</given-names></name><name><surname>Portier</surname><given-names>H</given-names></name><name><surname>Blettery</surname><given-names>B</given-names></name></person-group><year>2003</year><article-title>Candidemia in critically ill patients: difference of outcome between medical and surgical patients</article-title><source>Intensive Care Medicine</source><volume>29</volume><fpage>2162</fpage><lpage>2169</lpage><pub-id pub-id-type="doi">10.1007/s00134-003-2002-x</pub-id></element-citation></ref><ref id="bib12"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Connolly</surname><given-names>LA</given-names></name><name><surname>Riccombeni</surname><given-names>A</given-names></name><name><surname>Grózer</surname><given-names>Z</given-names></name><name><surname>Holland</surname><given-names>LM</given-names></name><name><surname>Lynch</surname><given-names>DB</given-names></name><name><surname>Andes</surname><given-names>DR</given-names></name><name><surname>Butler</surname><given-names>G</given-names></name></person-group><year>2013</year><article-title>The APSES transcription factor Efg1 is a global regulator that controls morphogenesis and biofilm formation in Candida parapsilosis</article-title><source>Molecular Microbiology</source><volume>90</volume><fpage>36</fpage><lpage>53</lpage><pub-id pub-id-type="doi">10.1111/mmi.12345</pub-id></element-citation></ref><ref id="bib13"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Coste</surname><given-names>A</given-names></name><name><surname>Turner</surname><given-names>V</given-names></name><name><surname>Ischer</surname><given-names>F</given-names></name><name><surname>Morschhauser</surname><given-names>J</given-names></name><name><surname>Forche</surname><given-names>A</given-names></name><name><surname>Selmecki</surname><given-names>A</given-names></name><name><surname>Berman</surname><given-names>J</given-names></name><name><surname>Bille</surname><given-names>J</given-names></name><name><surname>Sanglard</surname><given-names>D</given-names></name></person-group><year>2006</year><article-title>A mutation in Tac1p, a transcription factor regulating CDR1 and CDR2, is coupled with loss of heterozygosity at chromosome 5 to mediate antifungal resistance in <italic>Candida albicans</italic></article-title><source>Genetics</source><volume>172</volume><fpage>2139</fpage><lpage>2156</lpage><pub-id pub-id-type="doi">10.1534/genetics.105.054767</pub-id></element-citation></ref><ref id="bib14"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cowen</surname><given-names>LE</given-names></name><name><surname>Anderson</surname><given-names>JB</given-names></name><name><surname>Kohn</surname><given-names>LM</given-names></name></person-group><year>2002</year><article-title>Evolution of drug resistance in Candida albicans</article-title><source>Annual Review of Microbiology</source><volume>56</volume><fpage>139</fpage><lpage>165</lpage><pub-id pub-id-type="doi">10.1146/annurev.micro.56.012302.160907</pub-id></element-citation></ref><ref id="bib15"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cowen</surname><given-names>LE</given-names></name><name><surname>Lindquist</surname><given-names>S</given-names></name></person-group><year>2005</year><article-title>Hsp90 potentiates the rapid evolution of new traits: drug resistance in diverse fungi</article-title><source>Science</source><volume>309</volume><fpage>2185</fpage><lpage>2189</lpage><pub-id pub-id-type="doi">10.1126/science.1118370</pub-id></element-citation></ref><ref id="bib16"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>De Groot</surname><given-names>PW</given-names></name><name><surname>Hellingwerf</surname><given-names>KJ</given-names></name><name><surname>Klis</surname><given-names>FM</given-names></name></person-group><year>2003</year><article-title>Genome-wide identification of fungal GPI proteins</article-title><source>Yeast</source><volume>20</volume><fpage>781</fpage><lpage>796</lpage><pub-id pub-id-type="doi">10.1002/yea.1007</pub-id></element-citation></ref><ref id="bib17"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ding</surname><given-names>L</given-names></name><name><surname>Ley</surname><given-names>TJ</given-names></name><name><surname>Larson</surname><given-names>DE</given-names></name><name><surname>Miller</surname><given-names>CA</given-names></name><name><surname>Koboldt</surname><given-names>DC</given-names></name><name><surname>Welch</surname><given-names>JS</given-names></name><name><surname>Ritchey</surname><given-names>JK</given-names></name><name><surname>Young</surname><given-names>MA</given-names></name><name><surname>Lamprecht</surname><given-names>T</given-names></name><name><surname>McLellan</surname><given-names>MD</given-names></name><name><surname>McMichael</surname><given-names>JF</given-names></name><name><surname>Wallis</surname><given-names>JW</given-names></name><name><surname>Lu</surname><given-names>C</given-names></name><name><surname>Shen</surname><given-names>D</given-names></name><name><surname>Harris</surname><given-names>CC</given-names></name><name><surname>Dooling</surname><given-names>DJ</given-names></name><name><surname>Fulton</surname><given-names>RS</given-names></name><name><surname>Fulton</surname><given-names>LL</given-names></name><name><surname>Chen</surname><given-names>K</given-names></name><name><surname>Schmidt</surname><given-names>H</given-names></name><name><surname>Kalicki-Veizer</surname><given-names>J</given-names></name><name><surname>Magrini</surname><given-names>VJ</given-names></name><name><surname>Cook</surname><given-names>L</given-names></name><name><surname>McGrath</surname><given-names>SD</given-names></name><name><surname>Vickery</surname><given-names>TL</given-names></name><name><surname>Wendl</surname><given-names>MC</given-names></name><name><surname>Heath</surname><given-names>S</given-names></name><name><surname>Watson</surname><given-names>MA</given-names></name><name><surname>Link</surname><given-names>DC</given-names></name><name><surname>Tomasson</surname><given-names>MH</given-names></name><name><surname>Shannon</surname><given-names>WD</given-names></name><name><surname>Payton</surname><given-names>JE</given-names></name><name><surname>Kulkarni</surname><given-names>S</given-names></name><name><surname>Westervelt</surname><given-names>P</given-names></name><name><surname>Walter</surname><given-names>MJ</given-names></name><name><surname>Graubert</surname><given-names>TA</given-names></name><name><surname>Mardis</surname><given-names>ER</given-names></name><name><surname>Wilson</surname><given-names>RK</given-names></name><name><surname>DiPersio</surname><given-names>JF</given-names></name></person-group><year>2012</year><article-title>Clonal evolution in relapsed acute myeloid leukaemia revealed by whole-genome sequencing</article-title><source>Nature</source><volume>481</volume><fpage>506</fpage><lpage>510</lpage><pub-id pub-id-type="doi">10.1038/nature10738</pub-id></element-citation></ref><ref id="bib18"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dunkel</surname><given-names>N</given-names></name><name><surname>Blass</surname><given-names>J</given-names></name><name><surname>Rogers</surname><given-names>PD</given-names></name><name><surname>Morschhauser</surname><given-names>J</given-names></name></person-group><year>2008</year><article-title>Mutations in the multi-drug resistance regulator MRR1, followed by loss of heterozygosity, are the main cause of MDR1 overexpression in fluconazole-resistant Candida albicans strains</article-title><source>Molecular Microbiology</source><volume>69</volume><fpage>827</fpage><lpage>840</lpage><pub-id pub-id-type="doi">10.1111/j.1365-2958.2008.06309.x</pub-id></element-citation></ref><ref id="bib19"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dunkel</surname><given-names>N</given-names></name><name><surname>Morschhauser</surname><given-names>J</given-names></name></person-group><year>2011</year><article-title>Loss of heterozygosity at an unlinked genomic locus is responsible for the phenotype of a Candida albicans sap4Δ sap5Δ sap6Δ mutant</article-title><source>Eukaryotic Cell</source><volume>10</volume><fpage>54</fpage><lpage>62</lpage><pub-id pub-id-type="doi">10.1128/EC.00281-10</pub-id></element-citation></ref><ref id="bib20"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Forche</surname><given-names>A</given-names></name><name><surname>Abbey</surname><given-names>D</given-names></name><name><surname>Pisithkul</surname><given-names>T</given-names></name><name><surname>Weinzierl</surname><given-names>MA</given-names></name><name><surname>Ringstrom</surname><given-names>T</given-names></name><name><surname>Bruck</surname><given-names>D</given-names></name><name><surname>Berman</surname><given-names>J</given-names></name></person-group><year>2011</year><article-title>Stress alters rates and types of loss of heterozygosity in Candida albicans</article-title><source>mBio</source><volume>2</volume><pub-id pub-id-type="doi">10.1128/mBio.00129-11</pub-id></element-citation></ref><ref id="bib21"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Forche</surname><given-names>A</given-names></name><name><surname>Alby</surname><given-names>K</given-names></name><name><surname>Schaefer</surname><given-names>D</given-names></name><name><surname>Johnson</surname><given-names>AD</given-names></name><name><surname>Berman</surname><given-names>J</given-names></name><name><surname>Bennett</surname><given-names>RJ</given-names></name></person-group><year>2008</year><article-title>The parasexual cycle in Candida albicans provides an alternative pathway to meiosis for the formation of recombinant strains</article-title><source>PLOS Biology</source><volume>6</volume><fpage>e110</fpage><pub-id pub-id-type="doi">10.1371/journal.pbio.0060110</pub-id></element-citation></ref><ref id="bib22"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Forche</surname><given-names>A</given-names></name><name><surname>Magee</surname><given-names>PT</given-names></name><name><surname>Magee</surname><given-names>BB</given-names></name><name><surname>May</surname><given-names>G</given-names></name></person-group><year>2004</year><article-title>Genome-wide single-nucleotide polymorphism map for Candida albicans</article-title><source>Eukaryotic Cell</source><volume>3</volume><fpage>705</fpage><lpage>714</lpage><pub-id pub-id-type="doi">10.1128/EC.3.3.705-714.2004</pub-id></element-citation></ref><ref id="bib23"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Forche</surname><given-names>A</given-names></name><name><surname>Magee</surname><given-names>PT</given-names></name><name><surname>Selmecki</surname><given-names>A</given-names></name><name><surname>Berman</surname><given-names>J</given-names></name><name><surname>May</surname><given-names>G</given-names></name></person-group><year>2009</year><article-title>Evolution in Candida albicans populations during a single passage through a mouse host</article-title><source>Genetics</source><volume>182</volume><fpage>799</fpage><lpage>811</lpage><pub-id pub-id-type="doi">10.1534/genetics.109.103325</pub-id></element-citation></ref><ref id="bib24"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Forche</surname><given-names>A</given-names></name><name><surname>May</surname><given-names>G</given-names></name><name><surname>Magee</surname><given-names>PT</given-names></name></person-group><year>2005</year><article-title>Demonstration of loss of heterozygosity by single-nucleotide polymorphism microarray analysis and alterations in strain morphology in Candida albicans strains during infection</article-title><source>Eukaryotic Cell</source><volume>4</volume><fpage>156</fpage><lpage>165</lpage><pub-id pub-id-type="doi">10.1128/EC.4.1.156-165.2005</pub-id></element-citation></ref><ref id="bib25"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gagneux</surname><given-names>S</given-names></name><name><surname>Long</surname><given-names>CD</given-names></name><name><surname>Small</surname><given-names>PM</given-names></name><name><surname>Van</surname><given-names>T</given-names></name><name><surname>Schoolnik</surname><given-names>GK</given-names></name><name><surname>Bohannan</surname><given-names>BJ</given-names></name></person-group><year>2006</year><article-title>The competitive cost of antibiotic resistance in Mycobacterium tuberculosis</article-title><source>Science</source><volume>312</volume><fpage>1944</fpage><lpage>1946</lpage><pub-id pub-id-type="doi">10.1126/science.1124410</pub-id></element-citation></ref><ref id="bib26"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Galhardo</surname><given-names>RS</given-names></name><name><surname>Hastings</surname><given-names>PJ</given-names></name><name><surname>Rosenberg</surname><given-names>SM</given-names></name></person-group><year>2007</year><article-title>Mutation as a stress response and the regulation of evolvability</article-title><source>Critical Reviews in Biochemistry and Molecular Biology</source><volume>42</volume><fpage>399</fpage><lpage>435</lpage><pub-id pub-id-type="doi">10.1080/10409230701648502</pub-id></element-citation></ref><ref id="bib27"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gemayel</surname><given-names>R</given-names></name><name><surname>Vinces</surname><given-names>MD</given-names></name><name><surname>Legendre</surname><given-names>M</given-names></name><name><surname>Verstrepen</surname><given-names>KJ</given-names></name></person-group><year>2010</year><article-title>Variable tandem repeats accelerate evolution of coding and regulatory sequences</article-title><source>Annual Review of Genetics</source><volume>44</volume><fpage>445</fpage><lpage>477</lpage><pub-id pub-id-type="doi">10.1146/annurev-genet-072610-155046</pub-id></element-citation></ref><ref id="bib28"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gow</surname><given-names>NA</given-names></name><name><surname>Hube</surname><given-names>B</given-names></name></person-group><year>2012</year><article-title>Importance of the Candida albicans cell wall during commensalism and infection</article-title><source>Current Opinion in Microbiology</source><volume>15</volume><fpage>406</fpage><lpage>412</lpage><pub-id pub-id-type="doi">10.1016/j.mib.2012.04.005</pub-id></element-citation></ref><ref id="bib29"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Goyard</surname><given-names>S</given-names></name><name><surname>Knechtle</surname><given-names>P</given-names></name><name><surname>Chauvel</surname><given-names>M</given-names></name><name><surname>Mallet</surname><given-names>A</given-names></name><name><surname>Prévost</surname><given-names>MC</given-names></name><name><surname>Proux</surname><given-names>C</given-names></name><name><surname>Coppée</surname><given-names>JY</given-names></name><name><surname>Schwarz</surname><given-names>P</given-names></name><name><surname>Dromer</surname><given-names>F</given-names></name><name><surname>Park</surname><given-names>H</given-names></name><name><surname>Filler</surname><given-names>SG</given-names></name><name><surname>Janbon</surname><given-names>G</given-names></name><name><surname>d'Enfert</surname><given-names>C</given-names></name></person-group><year>2008</year><article-title>The Yak1 kinase is involved in the initiation and maintenance of hyphal growth in Candida albicans</article-title><source>Molecular Biology of the Cell</source><volume>19</volume><fpage>2251</fpage><lpage>2266</lpage><pub-id pub-id-type="doi">10.1091/mbc.E07-09-0960</pub-id></element-citation></ref><ref id="bib30"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Grabherr</surname><given-names>MG</given-names></name><name><surname>Haas</surname><given-names>BJ</given-names></name><name><surname>Yassour</surname><given-names>M</given-names></name><name><surname>Levin</surname><given-names>JZ</given-names></name><name><surname>Thompson</surname><given-names>DA</given-names></name><name><surname>Amit</surname><given-names>I</given-names></name><name><surname>Adiconis</surname><given-names>X</given-names></name><name><surname>Fan</surname><given-names>L</given-names></name><name><surname>Raychowdhury</surname><given-names>R</given-names></name><name><surname>Zeng</surname><given-names>Q</given-names></name><name><surname>Chen</surname><given-names>Z</given-names></name><name><surname>Mauceli</surname><given-names>E</given-names></name><name><surname>Hacohen</surname><given-names>N</given-names></name><name><surname>Gnirke</surname><given-names>A</given-names></name><name><surname>Rhind</surname><given-names>N</given-names></name><name><surname>di Palma</surname><given-names>F</given-names></name><name><surname>Birren</surname><given-names>BW</given-names></name><name><surname>Nusbaum</surname><given-names>C</given-names></name><name><surname>Lindblad-Toh</surname><given-names>K</given-names></name><name><surname>Friedman</surname><given-names>N</given-names></name><name><surname>Regev</surname><given-names>A</given-names></name></person-group><year>2011</year><article-title>Full-length transcriptome assembly from RNA-Seq data without a reference genome</article-title><source>Nature Biotechnology</source><volume>29</volume><fpage>644</fpage><lpage>652</lpage><pub-id pub-id-type="doi">10.1038/nbt.1883</pub-id></element-citation></ref><ref id="bib31"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gresham</surname><given-names>D</given-names></name><name><surname>Desai</surname><given-names>MM</given-names></name><name><surname>Tucker</surname><given-names>CM</given-names></name><name><surname>Jenq</surname><given-names>HT</given-names></name><name><surname>Pai</surname><given-names>DA</given-names></name><name><surname>Ward</surname><given-names>A</given-names></name><name><surname>DeSevo</surname><given-names>CG</given-names></name><name><surname>Botstein</surname><given-names>D</given-names></name><name><surname>Dunham</surname><given-names>MJ</given-names></name></person-group><year>2008</year><article-title>The repertoire and dynamics of evolutionary adaptations to controlled nutrient-limited environments in yeast</article-title><source>PLOS Genetics</source><volume>4</volume><pub-id pub-id-type="doi">10.1371/Journal.Pgen.1000303</pub-id></element-citation></ref><ref id="bib32"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gudlaugsson</surname><given-names>O</given-names></name><name><surname>Gillespie</surname><given-names>S</given-names></name><name><surname>Lee</surname><given-names>K</given-names></name><name><surname>Vande Berg</surname><given-names>J</given-names></name><name><surname>Hu</surname><given-names>J</given-names></name><name><surname>Messer</surname><given-names>S</given-names></name><name><surname>Herwaldt</surname><given-names>L</given-names></name><name><surname>Pfaller</surname><given-names>M</given-names></name><name><surname>Diekema</surname><given-names>D</given-names></name></person-group><year>2003</year><article-title>Attributable mortality of nosocomial candidemia, revisited</article-title><source>Clinical Infectious Diseases</source><volume>37</volume><fpage>1172</fpage><lpage>1177</lpage><pub-id pub-id-type="doi">10.1086/378745</pub-id></element-citation></ref><ref id="bib33"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hickman</surname><given-names>MA</given-names></name><name><surname>Zeng</surname><given-names>G</given-names></name><name><surname>Forche</surname><given-names>A</given-names></name><name><surname>Hirakawa</surname><given-names>MP</given-names></name><name><surname>Abbey</surname><given-names>D</given-names></name><name><surname>Harrison</surname><given-names>BD</given-names></name><name><surname>Berman</surname><given-names>J</given-names></name></person-group><year>2013</year><article-title>The 'obligate diploid' Candida albicans forms mating-competent haploids</article-title><source>Nature</source><volume>494</volume><fpage>55</fpage><lpage>59</lpage><pub-id pub-id-type="doi">10.1038/nature11865</pub-id></element-citation></ref><ref id="bib34"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hoyer</surname><given-names>LL</given-names></name><name><surname>Green</surname><given-names>CB</given-names></name><name><surname>Oh</surname><given-names>SH</given-names></name><name><surname>Zhao</surname><given-names>X</given-names></name></person-group><year>2008</year><article-title>Discovering the secrets of the Candida albicans agglutinin-like sequence (ALS) gene family–a sticky pursuit</article-title><source>Medical Mycology</source><volume>46</volume><fpage>1</fpage><lpage>15</lpage><pub-id pub-id-type="doi">10.1080/13693780701435317</pub-id></element-citation></ref><ref id="bib35"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hoyer</surname><given-names>LL</given-names></name><name><surname>Payne</surname><given-names>TL</given-names></name><name><surname>Bell</surname><given-names>M</given-names></name><name><surname>Myers</surname><given-names>AM</given-names></name><name><surname>Scherer</surname><given-names>S</given-names></name></person-group><year>1998</year><article-title>Candida albicans ALS3 and insights into the nature of the ALS gene family</article-title><source>Current Genetics</source><volume>33</volume><fpage>451</fpage><lpage>459</lpage><pub-id pub-id-type="doi">10.1007/s002940050359</pub-id></element-citation></ref><ref id="bib36"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Inglis</surname><given-names>DO</given-names></name><name><surname>Arnaud</surname><given-names>MB</given-names></name><name><surname>Binkley</surname><given-names>J</given-names></name><name><surname>Shah</surname><given-names>P</given-names></name><name><surname>Skrzypek</surname><given-names>MS</given-names></name><name><surname>Wymore</surname><given-names>F</given-names></name><name><surname>Sherlock</surname><given-names>G</given-names></name></person-group><year>2012</year><article-title>The Candida genome database incorporates multiple Candida species: multispecies search and analysis tools with curated gene and protein information for Candida albicans and Candida glabrata</article-title><source>Nucleic Acids Research</source><volume>40</volume><fpage>D667</fpage><lpage>D674</lpage><pub-id pub-id-type="doi">10.1093/nar/gkr945</pub-id></element-citation></ref><ref id="bib37"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Inglis</surname><given-names>DO</given-names></name><name><surname>Skrzypek</surname><given-names>MS</given-names></name><name><surname>Arnaud</surname><given-names>MB</given-names></name><name><surname>Binkley</surname><given-names>J</given-names></name><name><surname>Shah</surname><given-names>P</given-names></name><name><surname>Wymore</surname><given-names>F</given-names></name><name><surname>Sherlock</surname><given-names>G</given-names></name></person-group><year>2013</year><article-title>Improved gene ontology annotation for biofilm formation, filamentous growth, and phenotypic switching in Candida albicans</article-title><source>Eukaryotic Cell</source><volume>12</volume><fpage>101</fpage><lpage>108</lpage><pub-id pub-id-type="doi">10.1128/EC.00238-12</pub-id></element-citation></ref><ref id="bib38"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jain</surname><given-names>C</given-names></name><name><surname>Pastor</surname><given-names>K</given-names></name><name><surname>Gonzalez</surname><given-names>AY</given-names></name><name><surname>Lorenz</surname><given-names>MC</given-names></name><name><surname>Rao</surname><given-names>RP</given-names></name></person-group><year>2013</year><article-title>The role of Candida albicans AP-1 protein against host derived ROS in in vivo models of infection</article-title><source>Virulence</source><volume>4</volume><fpage>67</fpage><lpage>76</lpage><pub-id pub-id-type="doi">10.4161/viru.22700</pub-id></element-citation></ref><ref id="bib39"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jain</surname><given-names>C</given-names></name><name><surname>Yun</surname><given-names>M</given-names></name><name><surname>Politz</surname><given-names>SM</given-names></name><name><surname>Rao</surname><given-names>RP</given-names></name></person-group><year>2009</year><article-title>A pathogenesis assay using Saccharomyces cerevisiae and Caenorhabditis elegans reveals novel roles for yeast AP-1, Yap1, and host dual oxidase BLI-3 in fungal pathogenesis</article-title><source>Eukaryotic Cell</source><volume>8</volume><fpage>1218</fpage><lpage>1227</lpage><pub-id pub-id-type="doi">10.1128/EC.00367-08</pub-id></element-citation></ref><ref id="bib40"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kemp</surname><given-names>HA</given-names></name><name><surname>Sprague</surname><given-names>GF</given-names><suffix>Jr</suffix></name></person-group><year>2003</year><article-title>Far3 and five interacting proteins prevent premature recovery from pheromone arrest in the budding yeast Saccharomyces cerevisiae</article-title><source>Molecular and Cellular Biology</source><volume>23</volume><fpage>1750</fpage><lpage>1763</lpage><pub-id pub-id-type="doi">10.1128/MCB.23.5.1750-1763.2003</pub-id></element-citation></ref><ref id="bib41"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kim</surname><given-names>J</given-names></name><name><surname>Sudbery</surname><given-names>P</given-names></name></person-group><year>2011</year><article-title>Candida albicans, a major human fungal pathogen</article-title><source>Journal of microbiology</source><volume>49</volume><fpage>171</fpage><lpage>177</lpage><pub-id pub-id-type="doi">10.1007/s12275-011-1064-7</pub-id></element-citation></ref><ref id="bib42"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kumamoto</surname><given-names>CA</given-names></name><name><surname>Pierce</surname><given-names>JV</given-names></name></person-group><year>2011</year><article-title>Immunosensing during colonization by Candida albicans: does it take a village to colonize the intestine?</article-title><source>Trends in Microbiology</source><volume>19</volume><fpage>263</fpage><lpage>267</lpage><pub-id pub-id-type="doi">10.1016/j.tim.2011.01.009</pub-id></element-citation></ref><ref id="bib43"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kumamoto</surname><given-names>CA</given-names></name><name><surname>Vinces</surname><given-names>MD</given-names></name></person-group><year>2005</year><article-title>Alternative Candida albicans lifestyles: growth on surfaces</article-title><source>Annual Review of Microbiology</source><volume>59</volume><fpage>113</fpage><lpage>133</lpage><pub-id pub-id-type="doi">10.1146/annurev.micro.59.030804.121034</pub-id></element-citation></ref><ref id="bib44"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Landau</surname><given-names>DA</given-names></name><name><surname>Carter</surname><given-names>SL</given-names></name><name><surname>Stojanov</surname><given-names>P</given-names></name><name><surname>McKenna</surname><given-names>A</given-names></name><name><surname>Stevenson</surname><given-names>K</given-names></name><name><surname>Lawrence</surname><given-names>MS</given-names></name><name><surname>Sougnez</surname><given-names>C</given-names></name><name><surname>Stewart</surname><given-names>C</given-names></name><name><surname>Sivachenko</surname><given-names>A</given-names></name><name><surname>Wang</surname><given-names>L</given-names></name><name><surname>Wan</surname><given-names>Y</given-names></name><name><surname>Zhang</surname><given-names>W</given-names></name><name><surname>Shukla</surname><given-names>SA</given-names></name><name><surname>Vartanov</surname><given-names>A</given-names></name><name><surname>Fernandes</surname><given-names>SM</given-names></name><name><surname>Saksena</surname><given-names>G</given-names></name><name><surname>Cibulskis</surname><given-names>K</given-names></name><name><surname>Tesar</surname><given-names>B</given-names></name><name><surname>Gabriel</surname><given-names>S</given-names></name><name><surname>Hacohen</surname><given-names>N</given-names></name><name><surname>Meyerson</surname><given-names>M</given-names></name><name><surname>Lander</surname><given-names>ES</given-names></name><name><surname>Neuberg</surname><given-names>D</given-names></name><name><surname>Brown</surname><given-names>JR</given-names></name><name><surname>Getz</surname><given-names>G</given-names></name><name><surname>Wu</surname><given-names>CJ</given-names></name></person-group><year>2013</year><article-title>Evolution and impact of subclonal mutations in chronic lymphocytic leukemia</article-title><source>Cell</source><volume>152</volume><fpage>714</fpage><lpage>726</lpage><pub-id pub-id-type="doi">10.1016/j.cell.2013.01.019</pub-id></element-citation></ref><ref id="bib45"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname><given-names>H</given-names></name><name><surname>Durbin</surname><given-names>R</given-names></name></person-group><year>2009</year><article-title>Fast and accurate short read alignment with Burrows-Wheeler transform</article-title><source>Bioinformatics</source><volume>25</volume><fpage>1754</fpage><lpage>1760</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/btp324</pub-id></element-citation></ref><ref id="bib46"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lo</surname><given-names>HJ</given-names></name><name><surname>Köhler</surname><given-names>JR</given-names></name><name><surname>DiDomenico</surname><given-names>B</given-names></name><name><surname>Loebenberg</surname><given-names>D</given-names></name><name><surname>Cacciapuoti</surname><given-names>A</given-names></name><name><surname>Fink</surname><given-names>GR</given-names></name></person-group><year>1997</year><article-title>Nonfilamentous C. albicans mutants are avirulent</article-title><source>Cell</source><volume>90</volume><fpage>939</fpage><lpage>949</lpage><pub-id pub-id-type="doi">10.1016/s0092-8674(00)80358-x</pub-id></element-citation></ref><ref id="bib47"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Loeb</surname><given-names>JD</given-names></name><name><surname>Sepulveda-Becerra</surname><given-names>M</given-names></name><name><surname>Hazan</surname><given-names>I</given-names></name><name><surname>Liu</surname><given-names>H</given-names></name></person-group><year>1999</year><article-title>A G1 cyclin is necessary for maintenance of filamentous growth in Candida albicans</article-title><source>Molecular and Cellular Biology</source><volume>19</volume><fpage>4019</fpage><lpage>4027</lpage></element-citation></ref><ref id="bib48"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>MacPherson</surname><given-names>S</given-names></name><name><surname>Akache</surname><given-names>B</given-names></name><name><surname>Weber</surname><given-names>S</given-names></name><name><surname>De Deken</surname><given-names>X</given-names></name><name><surname>Raymond</surname><given-names>M</given-names></name><name><surname>Turcotte</surname><given-names>B</given-names></name></person-group><year>2005</year><article-title>Candida albicans zinc cluster protein Upc2p confers resistance to antifungal drugs and is an activator of ergosterol biosynthetic genes</article-title><source>Antimicrobial Agents and Chemotherapy</source><volume>49</volume><fpage>1745</fpage><lpage>1752</lpage><pub-id pub-id-type="doi">10.1128/AAC.49.5.1745-1752.2005</pub-id></element-citation></ref><ref id="bib49"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mardis</surname><given-names>ER</given-names></name></person-group><year>2008</year><article-title>The impact of next-generation sequencing technology on genetics</article-title><source>Trends in Genetics</source><volume>24</volume><fpage>133</fpage><lpage>141</lpage><pub-id pub-id-type="doi">10.1016/j.tig.2007.12.007</pub-id></element-citation></ref><ref id="bib50"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Marichal</surname><given-names>P</given-names></name><name><surname>Koymans</surname><given-names>L</given-names></name><name><surname>Willemsens</surname><given-names>S</given-names></name><name><surname>Bellens</surname><given-names>D</given-names></name><name><surname>Verhasselt</surname><given-names>P</given-names></name><name><surname>Luyten</surname><given-names>W</given-names></name><name><surname>Borgers</surname><given-names>M</given-names></name><name><surname>Ramaekers</surname><given-names>FC</given-names></name><name><surname>Odds</surname><given-names>FC</given-names></name><name><surname>Bossche</surname><given-names>HV</given-names></name></person-group><year>1999</year><article-title>Contribution of mutations in the cytochrome P450 14alpha-demethylase (Erg11p, Cyp51p) to azole resistance in Candida albicans</article-title><source>Nature Reviews Microbiology</source><volume>145</volume><fpage>2701</fpage><lpage>2713</lpage></element-citation></ref><ref id="bib51"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>McKenna</surname><given-names>A</given-names></name><name><surname>Hanna</surname><given-names>M</given-names></name><name><surname>Banks</surname><given-names>E</given-names></name><name><surname>Sivachenko</surname><given-names>A</given-names></name><name><surname>Cibulskis</surname><given-names>K</given-names></name><name><surname>Kernytsky</surname><given-names>A</given-names></name><name><surname>Garimella</surname><given-names>K</given-names></name><name><surname>Altshuler</surname><given-names>D</given-names></name><name><surname>Gabriel</surname><given-names>S</given-names></name><name><surname>Daly</surname><given-names>M</given-names></name><name><surname>DePristo</surname><given-names>MA</given-names></name></person-group><year>2010</year><article-title>The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data</article-title><source>Genome Research</source><volume>20</volume><fpage>1297</fpage><lpage>1303</lpage><pub-id pub-id-type="doi">10.1101/gr.107524.110</pub-id></element-citation></ref><ref id="bib52"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nobile</surname><given-names>CJ</given-names></name><name><surname>Fox</surname><given-names>EP</given-names></name><name><surname>Nett</surname><given-names>JE</given-names></name><name><surname>Sorrells</surname><given-names>TR</given-names></name><name><surname>Mitrovich</surname><given-names>QM</given-names></name><name><surname>Hernday</surname><given-names>AD</given-names></name><name><surname>Johnson</surname><given-names>AD</given-names></name></person-group><year>2012</year><article-title>A recently evolved transcriptional network controls biofilm development in Candida albicans</article-title><source>Cell</source><volume>148</volume><fpage>126</fpage><lpage>138</lpage><pub-id pub-id-type="doi">10.1016/j.cell.2011.10.048</pub-id></element-citation></ref><ref id="bib53"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nobile</surname><given-names>CJ</given-names></name><name><surname>Mitchell</surname><given-names>AP</given-names></name></person-group><year>2005</year><article-title>Regulation of cell-surface genes and biofilm formation by the C. albicans transcription factor Bcr1p</article-title><source>Current Biology</source><volume>15</volume><fpage>1150</fpage><lpage>1155</lpage><pub-id pub-id-type="doi">10.1016/j.cub.2005.05.047</pub-id></element-citation></ref><ref id="bib54"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Noble</surname><given-names>SM</given-names></name><name><surname>French</surname><given-names>S</given-names></name><name><surname>Kohn</surname><given-names>LA</given-names></name><name><surname>Chen</surname><given-names>V</given-names></name><name><surname>Johnson</surname><given-names>AD</given-names></name></person-group><year>2010</year><article-title>Systematic screens of a Candida albicans homozygous deletion library decouple morphogenetic switching and pathogenicity</article-title><source>Nature Genetics</source><volume>42</volume><fpage>590</fpage><lpage>598</lpage><pub-id pub-id-type="doi">10.1038/ng.605</pub-id></element-citation></ref><ref id="bib55"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Oliver</surname><given-names>BG</given-names></name><name><surname>Song</surname><given-names>JL</given-names></name><name><surname>Choiniere</surname><given-names>JH</given-names></name><name><surname>White</surname><given-names>TC</given-names></name></person-group><year>2007</year><article-title>cis-Acting elements within the Candida albicans ERG11 promoter mediate the azole response through transcription factor Upc2p</article-title><source>Eukaryotic Cell</source><volume>6</volume><fpage>2231</fpage><lpage>2239</lpage><pub-id pub-id-type="doi">10.1128/EC.00331-06</pub-id></element-citation></ref><ref id="bib56"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pappas</surname><given-names>PG</given-names></name><name><surname>Rex</surname><given-names>JH</given-names></name><name><surname>Lee</surname><given-names>J</given-names></name><name><surname>Hamill</surname><given-names>RJ</given-names></name><name><surname>Larsen</surname><given-names>RA</given-names></name><name><surname>Powderly</surname><given-names>W</given-names></name><name><surname>Dismukes</surname><given-names>WE</given-names></name></person-group><year>2003</year><article-title>A prospective observational study of candidemia: epidemiology, therapy, and influences on mortality in hospitalized adult and pediatric patients</article-title><source>Clinical Infectious Diseases</source><volume>37</volume><fpage>634</fpage><lpage>643</lpage><pub-id pub-id-type="doi">10.1086/376906</pub-id></element-citation></ref><ref id="bib57"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pavelka</surname><given-names>N</given-names></name><name><surname>Rancati</surname><given-names>G</given-names></name><name><surname>Zhu</surname><given-names>J</given-names></name><name><surname>Bradford</surname><given-names>WD</given-names></name><name><surname>Saraf</surname><given-names>A</given-names></name><name><surname>Florens</surname><given-names>L</given-names></name><name><surname>Sanderson</surname><given-names>BW</given-names></name><name><surname>Hattem</surname><given-names>GL</given-names></name><name><surname>Li</surname><given-names>R</given-names></name></person-group><year>2010</year><article-title>Aneuploidy confers quantitative proteome changes and phenotypic variation in budding yeast</article-title><source>Nature</source><volume>468</volume><fpage>321</fpage><lpage>325</lpage><pub-id pub-id-type="doi">10.1038/nature09529</pub-id></element-citation></ref><ref id="bib58"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Perea</surname><given-names>S</given-names></name><name><surname>López-Ribot</surname><given-names>JL</given-names></name><name><surname>Kirkpatrick</surname><given-names>WR</given-names></name><name><surname>McAtee</surname><given-names>RK</given-names></name><name><surname>Santillan</surname><given-names>RA</given-names></name><name><surname>Martinez</surname><given-names>M</given-names></name><name><surname>Calabrese</surname><given-names>D</given-names></name><name><surname>Sanglard</surname><given-names>D</given-names></name><name><surname>Patterson</surname><given-names>TF</given-names></name></person-group><year>2001</year><article-title>Prevalence of molecular mechanisms of resistance to azole antifungal agents in Candida albicans strains displaying high-level fluconazole resistance isolated from human immunodeficiency virus-infected patients</article-title><source>Antimicrobial Agents and Chemotherapy</source><volume>45</volume><fpage>2676</fpage><lpage>2684</lpage><pub-id pub-id-type="doi">10.1128/AAC.45.10.2676-2684.2001</pub-id></element-citation></ref><ref id="bib59"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pfaller</surname><given-names>MA</given-names></name><name><surname>Lockhart</surname><given-names>SR</given-names></name><name><surname>Pujol</surname><given-names>C</given-names></name><name><surname>Swails-Wenger</surname><given-names>JA</given-names></name><name><surname>Messer</surname><given-names>SA</given-names></name><name><surname>Edmond</surname><given-names>MB</given-names></name><name><surname>Jones</surname><given-names>RN</given-names></name><name><surname>Wenzel</surname><given-names>RP</given-names></name><name><surname>Soll</surname><given-names>DR</given-names></name></person-group><year>1998</year><article-title>Hospital specificity, region specificity, and fluconazole resistance of Candida albicans bloodstream isolates</article-title><source>Journal of Clinical Microbiology</source><volume>36</volume><fpage>1518</fpage><lpage>1529</lpage></element-citation></ref><ref id="bib60"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pittet</surname><given-names>D</given-names></name><name><surname>Monod</surname><given-names>M</given-names></name><name><surname>Suter</surname><given-names>PM</given-names></name><name><surname>Frenk</surname><given-names>E</given-names></name><name><surname>Auckenthaler</surname><given-names>R</given-names></name></person-group><year>1994</year><article-title>Candida colonization and subsequent infections in critically ill surgical patients</article-title><source>Annals of Surgery</source><volume>220</volume><fpage>751</fpage><lpage>758</lpage><pub-id pub-id-type="doi">10.1097/00000658-199412000-00008</pub-id></element-citation></ref><ref id="bib61"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Podlaha</surname><given-names>O</given-names></name><name><surname>Riester</surname><given-names>M</given-names></name><name><surname>De</surname><given-names>S</given-names></name><name><surname>Michor</surname><given-names>F</given-names></name></person-group><year>2012</year><article-title>Evolution of the cancer genome</article-title><source>Trends in Genetics</source><volume>28</volume><fpage>155</fpage><lpage>163</lpage><pub-id pub-id-type="doi">10.1016/j.tig.2012.01.003</pub-id></element-citation></ref><ref id="bib62"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ponder</surname><given-names>RG</given-names></name><name><surname>Fonville</surname><given-names>NC</given-names></name><name><surname>Rosenberg</surname><given-names>SM</given-names></name></person-group><year>2005</year><article-title>A switch from high-fidelity to error-prone DNA double-strand break repair underlies stress-induced mutation</article-title><source>Molecular Cell</source><volume>19</volume><fpage>791</fpage><lpage>804</lpage><pub-id pub-id-type="doi">10.1016/j.molcel.2005.07.025</pub-id></element-citation></ref><ref id="bib63"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Reich</surname><given-names>M</given-names></name><name><surname>Liefeld</surname><given-names>T</given-names></name><name><surname>Gould</surname><given-names>J</given-names></name><name><surname>Lerner</surname><given-names>J</given-names></name><name><surname>Tamayo</surname><given-names>P</given-names></name><name><surname>Mesirov</surname><given-names>JP</given-names></name></person-group><year>2006</year><article-title>GenePattern 2.0</article-title><source>Nature Genetics</source><volume>38</volume><fpage>500</fpage><lpage>501</lpage><pub-id pub-id-type="doi">10.1038/ng0506-500</pub-id></element-citation></ref><ref id="bib64"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rex</surname><given-names>JH</given-names></name><name><surname>Bennett</surname><given-names>JE</given-names></name><name><surname>Sugar</surname><given-names>AM</given-names></name><name><surname>Pappas</surname><given-names>PG</given-names></name><name><surname>van der Horst</surname><given-names>CM</given-names></name><name><surname>Edwards</surname><given-names>JE</given-names></name><name><surname>Washburn</surname><given-names>RG</given-names></name><name><surname>Scheld</surname><given-names>WM</given-names></name><name><surname>Karchmer</surname><given-names>AW</given-names></name><name><surname>Dine</surname><given-names>AP</given-names></name></person-group><year>1994</year><article-title>A randomized trial comparing fluconazole with amphotericin B for the treatment of candidemia in patients without neutropenia. Candidemia Study Group and the National Institute</article-title><source>The New England Journal of Medicine</source><volume>331</volume><fpage>1325</fpage><lpage>1330</lpage><pub-id pub-id-type="doi">10.1056/NEJM199411173312001</pub-id></element-citation></ref><ref id="bib65"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rex</surname><given-names>JH</given-names></name><name><surname>Pappas</surname><given-names>PG</given-names></name><name><surname>Karchmer</surname><given-names>AW</given-names></name><name><surname>Sobel</surname><given-names>J</given-names></name><name><surname>Edwards</surname><given-names>JE</given-names></name><name><surname>Hadley</surname><given-names>S</given-names></name><name><surname>Brass</surname><given-names>C</given-names></name><name><surname>Vazquez</surname><given-names>JA</given-names></name><name><surname>Chapman</surname><given-names>SW</given-names></name><name><surname>Horowitz</surname><given-names>HW</given-names></name><name><surname>Zervos</surname><given-names>M</given-names></name><name><surname>McKinsey</surname><given-names>D</given-names></name><name><surname>Lee</surname><given-names>J</given-names></name><name><surname>Babinchak</surname><given-names>T</given-names></name><name><surname>Bradsher</surname><given-names>RW</given-names></name><name><surname>Cleary</surname><given-names>JD</given-names></name><name><surname>Cohen</surname><given-names>DM</given-names></name><name><surname>Danziger</surname><given-names>L</given-names></name><name><surname>Goldman</surname><given-names>M</given-names></name><name><surname>Goodman</surname><given-names>J</given-names></name><name><surname>Hilton</surname><given-names>E</given-names></name><name><surname>Hyslop</surname><given-names>NE</given-names></name><name><surname>Kett</surname><given-names>DH</given-names></name><name><surname>Lutz</surname><given-names>J</given-names></name><name><surname>Rubin</surname><given-names>RH</given-names></name><name><surname>Scheld</surname><given-names>WM</given-names></name><name><surname>Schuster</surname><given-names>M</given-names></name><name><surname>Simmons</surname><given-names>B</given-names></name><name><surname>Stein</surname><given-names>DK</given-names></name><name><surname>Washburn</surname><given-names>RG</given-names></name><name><surname>Mautner</surname><given-names>L</given-names></name><name><surname>Chu</surname><given-names>TC</given-names></name><name><surname>Panzer</surname><given-names>H</given-names></name><name><surname>Rosenstein</surname><given-names>RB</given-names></name><name><surname>Booth</surname><given-names>J</given-names></name><collab>National Institute of Allergy and Infectious Diseases Mycoses Study Group</collab></person-group><year>2003</year><article-title>A randomized and blinded multicenter trial of high-dose fluconazole plus placebo versus fluconazole plus amphotericin B as therapy for candidemia and its consequences in nonneutropenic subjects</article-title><source>Clinical Infectious Diseases</source><volume>36</volume><fpage>1221</fpage><lpage>1228</lpage><pub-id pub-id-type="doi">10.1086/374850</pub-id></element-citation></ref><ref id="bib66"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Reynolds</surname><given-names>TB</given-names></name><name><surname>Fink</surname><given-names>GR</given-names></name></person-group><year>2001</year><article-title>Bakers' yeast, a model for fungal biofilm formation</article-title><source>Science</source><volume>291</volume><fpage>878</fpage><lpage>881</lpage><pub-id pub-id-type="doi">10.1126/science.291.5505.878</pub-id></element-citation></ref><ref id="bib67"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Robinson</surname><given-names>JT</given-names></name><name><surname>Thorvaldsdóttir</surname><given-names>H</given-names></name><name><surname>Winckler</surname><given-names>W</given-names></name><name><surname>Guttman</surname><given-names>M</given-names></name><name><surname>Lander</surname><given-names>ES</given-names></name><name><surname>Getz</surname><given-names>G</given-names></name><name><surname>Mesirov</surname><given-names>JP</given-names></name></person-group><year>2011</year><article-title>Integrative genomics viewer</article-title><source>Nature Biotechnology</source><volume>29</volume><fpage>24</fpage><lpage>26</lpage><pub-id pub-id-type="doi">10.1038/nbt.1754</pub-id></element-citation></ref><ref id="bib68"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rosenberg</surname><given-names>SM</given-names></name></person-group><year>2011</year><article-title>Stress-induced loss of heterozygosity in Candida: a possible missing link in the ability to evolve</article-title><source>mBio</source><volume>2</volume><pub-id pub-id-type="doi">10.1128/mBio.00200-11</pub-id></element-citation></ref><ref id="bib69"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Santos</surname><given-names>MA</given-names></name><name><surname>Moura</surname><given-names>G</given-names></name><name><surname>Massey</surname><given-names>SE</given-names></name><name><surname>Tuite</surname><given-names>MF</given-names></name></person-group><year>2004</year><article-title>Driving change: the evolution of alternative genetic codes</article-title><source>Trends in Genetics</source><volume>20</volume><fpage>95</fpage><lpage>102</lpage><pub-id pub-id-type="doi">10.1016/j.tig.2003.12.009</pub-id></element-citation></ref><ref id="bib70"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sasse</surname><given-names>C</given-names></name><name><surname>Dunkel</surname><given-names>N</given-names></name><name><surname>Schafer</surname><given-names>T</given-names></name><name><surname>Schneider</surname><given-names>S</given-names></name><name><surname>Dierolf</surname><given-names>F</given-names></name><name><surname>Ohlsen</surname><given-names>K</given-names></name><name><surname>Morschhauser</surname><given-names>J</given-names></name></person-group><year>2012</year><article-title>The stepwise acquisition of fluconazole resistance mutations causes a gradual loss of fitness in Candida albicans</article-title><source>Molecular Microbiology</source><volume>86</volume><fpage>539</fpage><lpage>556</lpage><pub-id pub-id-type="doi">10.1111/j.1365-2958.2012.08210.x</pub-id></element-citation></ref><ref id="bib71"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schubert</surname><given-names>S</given-names></name><name><surname>Popp</surname><given-names>C</given-names></name><name><surname>Rogers</surname><given-names>PD</given-names></name><name><surname>Morschhäuser</surname><given-names>J</given-names></name></person-group><year>2011</year><article-title>Functional dissection of a Candida albicans zinc cluster transcription factor, the multidrug resistance regulator Mrr1</article-title><source>Eukaryotic Cell</source><volume>10</volume><fpage>1110</fpage><lpage>1121</lpage><pub-id pub-id-type="doi">10.1128/EC.05100-11</pub-id></element-citation></ref><ref id="bib72"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Selmecki</surname><given-names>A</given-names></name><name><surname>Forche</surname><given-names>A</given-names></name><name><surname>Berman</surname><given-names>J</given-names></name></person-group><year>2006</year><article-title>Aneuploidy and isochromosome formation in drug-resistant Candida albicans</article-title><source>Science</source><volume>313</volume><fpage>367</fpage><lpage>370</lpage><pub-id pub-id-type="doi">10.1126/science.1128242</pub-id></element-citation></ref><ref id="bib73"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Selmecki</surname><given-names>A</given-names></name><name><surname>Forche</surname><given-names>A</given-names></name><name><surname>Berman</surname><given-names>J</given-names></name></person-group><year>2010</year><article-title>Genomic plasticity of the human fungal pathogen Candida albicans</article-title><source>Eukaryotic Cell</source><volume>9</volume><fpage>991</fpage><lpage>1008</lpage><pub-id pub-id-type="doi">10.1128/EC.00060-10</pub-id></element-citation></ref><ref id="bib74"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Selmecki</surname><given-names>A</given-names></name><name><surname>Gerami-Nejad</surname><given-names>M</given-names></name><name><surname>Paulson</surname><given-names>C</given-names></name><name><surname>Forche</surname><given-names>A</given-names></name><name><surname>Berman</surname><given-names>J</given-names></name></person-group><year>2008</year><article-title>An isochromosome confers drug resistance in vivo by amplification of two genes, ERG11 and TAC1</article-title><source>Molecular Microbiology</source><volume>68</volume><fpage>624</fpage><lpage>641</lpage><pub-id pub-id-type="doi">10.1111/j.1365-2958.2008.06176.x</pub-id></element-citation></ref><ref id="bib75"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Selmecki</surname><given-names>AM</given-names></name><name><surname>Dulmage</surname><given-names>K</given-names></name><name><surname>Cowen</surname><given-names>LE</given-names></name><name><surname>Anderson</surname><given-names>JB</given-names></name><name><surname>Berman</surname><given-names>J</given-names></name></person-group><year>2009</year><article-title>Acquisition of aneuploidy provides increased fitness during the evolution of antifungal drug resistance</article-title><source>PLOS Genetics</source><volume>5</volume><fpage>e1000705</fpage><pub-id pub-id-type="doi">10.1371/journal.pgen.1000705</pub-id></element-citation></ref><ref id="bib76"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sheppard</surname><given-names>DC</given-names></name><name><surname>Yeaman</surname><given-names>MR</given-names></name><name><surname>Welch</surname><given-names>WH</given-names></name><name><surname>Phan</surname><given-names>QT</given-names></name><name><surname>Fu</surname><given-names>Y</given-names></name><name><surname>Ibrahim</surname><given-names>AS</given-names></name><name><surname>Edwards</surname><given-names>JE</given-names><suffix>Jr</suffix></name></person-group><year>2004</year><article-title>Functional and structural diversity in the Als protein family of Candida albicans</article-title><source>Journal of Biological Chemistry</source><volume>279</volume><fpage>30480</fpage><lpage>30489</lpage><pub-id pub-id-type="doi">10.1074/jbc.M401929200</pub-id></element-citation></ref><ref id="bib77"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Singh-Babak</surname><given-names>SD</given-names></name><name><surname>Babak</surname><given-names>T</given-names></name><name><surname>Diezmann</surname><given-names>S</given-names></name><name><surname>Hill</surname><given-names>JA</given-names></name><name><surname>Xie</surname><given-names>JL</given-names></name><name><surname>Chen</surname><given-names>YL</given-names></name><name><surname>Poutanen</surname><given-names>SM</given-names></name><name><surname>Rennie</surname><given-names>RP</given-names></name><name><surname>Heitman</surname><given-names>J</given-names></name><name><surname>Cowen</surname><given-names>LE</given-names></name></person-group><year>2012</year><article-title>Global analysis of the evolution and mechanism of echinocandin resistance in Candida glabrata</article-title><source>PLOS Pathogens</source><volume>8</volume><fpage>e1002718</fpage><pub-id pub-id-type="doi">10.1371/journal.ppat.1002718</pub-id></element-citation></ref><ref id="bib78"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Storm</surname><given-names>N</given-names></name><name><surname>Darnhofer-Patel</surname><given-names>B</given-names></name><name><surname>van den Boom</surname><given-names>D</given-names></name><name><surname>Rodi</surname><given-names>CP</given-names></name></person-group><year>2003</year><article-title>MALDI-TOF mass spectrometry-based SNP genotyping</article-title><source>Methods in Molecular Biology</source><volume>212</volume><fpage>241</fpage><lpage>262</lpage></element-citation></ref><ref id="bib79"><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Swofford</surname><given-names>DL</given-names></name></person-group><year>2002</year><article-title>PAUP*. Phylogenetic analysis using parsimony (*and other methods)</article-title><publisher-loc>Sunderland, MA</publisher-loc><publisher-name>Sinauer Associates</publisher-name></element-citation></ref><ref id="bib80"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tang</surname><given-names>YC</given-names></name><name><surname>Amon</surname><given-names>A</given-names></name></person-group><year>2013</year><article-title>Gene copy-number alterations: a cost-benefit analysis</article-title><source>Cell</source><volume>152</volume><fpage>394</fpage><lpage>405</lpage><pub-id pub-id-type="doi">10.1016/j.cell.2012.11.043</pub-id></element-citation></ref><ref id="bib81"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Thompson</surname><given-names>DA</given-names></name><name><surname>Desai</surname><given-names>MM</given-names></name><name><surname>Murray</surname><given-names>AW</given-names></name></person-group><year>2006</year><article-title>Ploidy controls the success of mutators and nature of mutations during budding yeast evolution</article-title><source>Current Biology</source><volume>16</volume><fpage>1581</fpage><lpage>1590</lpage><pub-id pub-id-type="doi">10.1016/j.cub.2006.06.070</pub-id></element-citation></ref><ref id="bib82"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Torres</surname><given-names>EM</given-names></name><name><surname>Dephoure</surname><given-names>N</given-names></name><name><surname>Panneerselvam</surname><given-names>A</given-names></name><name><surname>Tucker</surname><given-names>CM</given-names></name><name><surname>Whittaker</surname><given-names>CA</given-names></name><name><surname>Gygi</surname><given-names>SP</given-names></name><name><surname>Amon</surname><given-names>A</given-names></name></person-group><year>2010</year><article-title>Identification of aneuploidy-tolerating mutations</article-title><source>Cell</source><volume>143</volume><fpage>71</fpage><lpage>83</lpage><pub-id pub-id-type="doi">10.1016/j.cell.2010.08.038</pub-id></element-citation></ref><ref id="bib83"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Uhl</surname><given-names>MA</given-names></name><name><surname>Biery</surname><given-names>M</given-names></name><name><surname>Craig</surname><given-names>N</given-names></name><name><surname>Johnson</surname><given-names>AD</given-names></name></person-group><year>2003</year><article-title>Haploinsufficiency-based large-scale forward genetic analysis of filamentous growth in the diploid human fungal pathogen C.albicans</article-title><source>The EMBO Journal</source><volume>22</volume><fpage>2668</fpage><lpage>2678</lpage><pub-id pub-id-type="doi">10.1093/emboj/cdg256</pub-id></element-citation></ref><ref id="bib84"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Verstrepen</surname><given-names>KJ</given-names></name><name><surname>Jansen</surname><given-names>A</given-names></name><name><surname>Lewitter</surname><given-names>F</given-names></name><name><surname>Fink</surname><given-names>GR</given-names></name></person-group><year>2005</year><article-title>Intragenic tandem repeats generate functional variability</article-title><source>Nature Genetics</source><volume>37</volume><fpage>986</fpage><lpage>990</lpage><pub-id pub-id-type="doi">10.1038/ng1618</pub-id></element-citation></ref><ref id="bib85"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>White</surname><given-names>TC</given-names></name></person-group><year>1997a</year><article-title>Increased mRNA levels of ERG16, CDR, and MDR1 correlate with increases in azole resistance in Candida albicans isolates from a patient infected with human immunodeficiency virus</article-title><source>Antimicrobial Agents and Chemotherapy</source><volume>41</volume><fpage>1482</fpage><lpage>1487</lpage></element-citation></ref><ref id="bib86"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>White</surname><given-names>TC</given-names></name></person-group><year>1997b</year><article-title>The presence of an R467K amino acid substitution and loss of allelic variation correlate with an azole-resistant lanosterol 14alpha demethylase in Candida albicans</article-title><source>Antimicrobial Agents and Chemotherapy</source><volume>41</volume><fpage>1488</fpage><lpage>1494</lpage></element-citation></ref><ref id="bib87"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>White</surname><given-names>TC</given-names></name><name><surname>Holleman</surname><given-names>S</given-names></name><name><surname>Dy</surname><given-names>F</given-names></name><name><surname>Mirels</surname><given-names>LF</given-names></name><name><surname>Stevens</surname><given-names>DA</given-names></name></person-group><year>2002</year><article-title>Resistance mechanisms in clinical isolates of Candida albicans</article-title><source>Antimicrobial Agents and Chemotherapy</source><volume>46</volume><fpage>1704</fpage><lpage>1713</lpage><pub-id pub-id-type="doi">10.1128/AAC.46.6.1704-1713.2002</pub-id></element-citation></ref><ref id="bib88"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wisplinghoff</surname><given-names>H</given-names></name><name><surname>Seifert</surname><given-names>H</given-names></name><name><surname>Tallent</surname><given-names>SM</given-names></name><name><surname>Bischoff</surname><given-names>T</given-names></name><name><surname>Wenzel</surname><given-names>RP</given-names></name><name><surname>Edmond</surname><given-names>MB</given-names></name></person-group><year>2003</year><article-title>Nosocomial bloodstream infections in pediatric patients in United States hospitals: epidemiology, clinical features and susceptibilities</article-title><source>The Pediatric Infectious Disease Journal</source><volume>22</volume><fpage>686</fpage><lpage>691</lpage><pub-id pub-id-type="doi">10.1097/01.inf.0000078159.53132.40</pub-id></element-citation></ref><ref id="bib89"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yona</surname><given-names>AH</given-names></name><name><surname>Manor</surname><given-names>YS</given-names></name><name><surname>Herbst</surname><given-names>RH</given-names></name><name><surname>Romano</surname><given-names>GH</given-names></name><name><surname>Mitchell</surname><given-names>A</given-names></name><name><surname>Kupiec</surname><given-names>M</given-names></name><name><surname>Pilpel</surname><given-names>Y</given-names></name><name><surname>Dahan</surname><given-names>O</given-names></name></person-group><year>2012</year><article-title>Chromosomal duplication is a transient evolutionary solution to stress</article-title><source>Proceedings of the National Academy of Sciences of USA</source><volume>109</volume><fpage>21010</fpage><lpage>21015</lpage><pub-id pub-id-type="doi">10.1073/pnas.1211150109</pub-id></element-citation></ref><ref id="bib90"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zakikhany</surname><given-names>K</given-names></name><name><surname>Naglik</surname><given-names>JR</given-names></name><name><surname>Schmidt-Westhausen</surname><given-names>A</given-names></name><name><surname>Holland</surname><given-names>G</given-names></name><name><surname>Schaller</surname><given-names>M</given-names></name><name><surname>Hube</surname><given-names>B</given-names></name></person-group><year>2007</year><article-title>In vivo transcript profiling of Candida albicans identifies a gene essential for interepithelial dissemination</article-title><source>Cellular Microbiology</source><volume>9</volume><fpage>2938</fpage><lpage>2954</lpage><pub-id pub-id-type="doi">10.1111/j.1462-5822.2007.01009.x</pub-id></element-citation></ref></ref-list></back><sub-article article-type="article-commentary" id="SA1"><front-stub><article-id pub-id-type="doi">10.7554/eLife.00662.022</article-id><title-group><article-title>Decision letter</article-title></title-group><contrib-group content-type="section"><contrib contrib-type="editor"><name><surname>Dermitzakis</surname><given-names>Emmanouil T</given-names></name><role>Reviewing editor</role><aff><institution>University of Geneva Medical School</institution>, <country>Switzerland</country></aff></contrib></contrib-group></front-stub><body><boxed-text><p>eLife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see <ext-link ext-link-type="uri" xlink:href="http://elifesciences.org/review-process">review process</ext-link>). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.</p></boxed-text><p>Thank you for sending your work entitled “The mutational landscape of gradual acquisition of drug resistance in clinical isolates of <italic>Candida albicans</italic>” for consideration at <italic>eLife</italic>. Your article has been evaluated by a Senior editor and 3 reviewers, one of whom is a member of our Board of Reviewing Editors. The manuscript was considered appropriate for <italic>eLife</italic> but some serious concerns need to be addressed successfully before a final decision can be reached.</p><p>The following individuals responsible for the peer review of your submission want to reveal their identity: Manolis Dermitzakis (Reviewing editor), Jacques Fellay (peer reviewer), and Daniel Wilson (peer reviewer).</p><p>The Reviewing editor and the other reviewers discussed their comments before we reached this decision, and the Reviewing editor has assembled the following comments to help you prepare a revised submission.</p><p>The manuscript by Funt and colleagues provide an in-depth analysis of whole genome sequence of <italic>C. albicans</italic> progressively in HIV patients to infer selective events and describe patterns of mutations, and those that survive and likely have an impact on resistance. Overall, the study is large and the data collected is very important and relevant with a degree of phenotypic assessment paired with the sequence data. The experimental and sequence analysis methodologies seem to be state of the art.</p><p>1) One major concern arises from consideration of all sections of the manuscript. There is very little statistical modeling to assess how likely these combinations of certain mutational patterns are due to real selective events due to specific phenotypic consequences or just by chance. Most of the manuscript is presented in this style, and it is hard to know which signals are truly significant. We provide two of the many examples in the manuscript to illustrate this point:</p><p>In the Results section, the authors state: “We consider… two to three isolates”. While this is a reasonable argument, without proper statistics it is not possible to evaluate it. Perhaps some degree of linkage disequilibrium would be useful to consider.</p><p>In the Results section, the authors also mention: “The recurring appearance of LOH events that coincide with changes in MIC suggests that they have been positively selected…” This is a “coincidence”, as the authors say, but is it significant?</p><p>2) The last part of the Results describes some validation experiments where the specific mutations are brought to an experimental setting with deletion mutants. This is an experiment that was necessary but surprisingly it was not given the proper weight (it was only a small fraction of the total text), plus it was not as conclusive as one would like regarding the MIC phenotype.</p><p>3) Considering the central role in the paper of the various classes of mutations called from the sequencing reads, I think that one technical issue in particular deserves further exploration and/or discussion. The SNP calling accuracy was found to be 87.3% on average (using a Sequenome iPLEX genotyping assay; Results, “Sequenome iPLEX genotyping assay”, and <xref ref-type="table" rid="tbl2">Table 2</xref>). There was however an important inter-sample variability, with very poor results observed in half of the tested series (<65% calling accuracy). This is worrying. Could this be due to intra-patient diversity, with different <italic>Candida</italic> subpopulations being sampled for the two assays? Is there any other potential explanation? If not, these inconsistencies should be listed as an important limitation of the study.</p><p>4) The paper is presented as a study of the within-host mutations that occur and are selected for in the <italic>Candida albicans</italic> genome between longitudinally sampled genomes before and after the introduction of antibiotic treatment. However, the data suggest that the variants that distinguish <italic>C. albicans</italic> genomes at different time points did not arise during the sampling window but were largely, maybe entirely, pre-existing (standing variation). This view is based on: (i) the very considerable genetic diversity (thousands of SNP differences) between isolates of the same strain within the same host, sampled less than a few months apart, in one case on the same day; (ii) the implausible mutation rate that would be needed to explain this variation as newly arising by mutation; and, (iii) the lack of analysis demonstrating that the data shows any signal of “measurable evolution”.</p></body></sub-article><sub-article article-type="reply" id="SA2"><front-stub><article-id pub-id-type="doi">10.7554/eLife.00662.023</article-id><title-group><article-title>Author response</article-title></title-group></front-stub><body><p>We thank the reviewers’ for their thoughtful comments, which we addressed in full in the revised manuscript. We first highlight the key aspects of the revision and then address each point in the response below. We realize that this re-submission has been delayed: as we explain below this is due to our decision to re-sequence all strains to substantially greater depth to ensure the reproducibility and validity of our results. This was critical as this is a resource intended for the wider community.</p><p>1) Extensive deeper sequencing of each strain shows high concordance in genotype calls: the reviewers pointed out that there was a lower than desired agreement between our Sequenome iPLEX assay genotyping and the corresponding base calls from our sequencing data (87% in the original study on average per patient, 77% averaged over all samples with some strains near 100% and others much lower; new <xref ref-type="fig" rid="fig1s1">Figure 1–supplement 1A</xref> in red). The reviewer raised the possibility of a sample switching. We first ruled out the possibility of a sample switch by repeating our genotyping from newly isolated DNA, but many of the discrepancies remained. We then examined the discordant calls and found that those were typically related with lower sequence coverage. Although our initial analysis used best practices at the time of original data collection, in some regions of lower sequencing coverage our GATK pipeline defaulted to the genotype of the provided reference genome, resulting in erroneous SNP calls when comparing between strains. Because our strains are also a resource used widely by the community, we decided to re-sequence all the strains to high coverage (53-283X, an average of 103X); this was not feasible when our data was initially collected, but major increases in sequencing yields made it readily possible now. The new sequencing data resolved the discrepancies, such that concordance between Sequenom iPlex and Illumina data reached 94% across all patients (see also the response to the issue 3 below; <xref ref-type="fig" rid="fig1s1">Figure 1–figure supplement 1A</xref> in blue). All the key results and conclusions that we originally reported were all reproduced in this deeper analysis.</p><p>2) The new data has added many new interesting genes to our list of those most likely to have adaptive polymorphisms. Our repeated sequencing revealed 4,756 persistent SNPs across all isolates and 240 recurrent genes. While previous clustering analysis of persistent and recurrent SNPs failed to produce significant functional enrichments, clustering of the recurrence matrix of our new data reveals significant functional enrichment in four clusters (revised <xref ref-type="fig" rid="fig5">Figure 5</xref>). Additionally, of the 240 recurrent genes identified, 23 were available as deletion library mutants (up from five previously), and we evaluated them for phenotypic differences (revised <xref ref-type="fig" rid="fig8">Figure 8</xref>). These new data significantly extend our previous manuscript, and have allowed us to make additional conclusions regarding drug resistance and adaptation to the host environment.</p><p>3) We have clarified the two possible evolutionary scenarios—de novo mutation and selection or selection on largely pre-existing variation. We now clearly note that, although the data favors selection of pre-existing variation as the predominant mode, we cannot rule out that some of the observed variation was due to de novo events that occurred in the sampling window.</p><p>4) We have provided statistical modeling where possible, including stochastic simulations to model the occurrence of (1) persistent non synonymous SNP-containing ORFs associated with changes in MIC, and (2) the expected degree of recurrence amongst these ORFs.</p><p>5) We have heeded the reviewers’ suggestion to expand our description and interpretation of the validation experiments with the available deletion mutants for recurrent loci; as the number of recurrent loci grew, we now tested 23 deletion mutants, three of which were significantly more fit than the parental strain (SN250) as determined by growth in standard tissue culture media (RPMI).</p><p>6) We assessed the rate and time scale of evolution by applying Bayesian Evolutionary Analysis by Sampling Trees (BEAST) in the one patient series with a sufficient density of longitudinal samples and have included that analysis.</p><p>7) We addressed all other comments, including clarification of our use of the word “clonal” as detailed in the point-by-point response below.</p><p>Response to specific comments:</p><p><italic>The manuscript by Funt and colleagues provide an in-depth analysis of whole genome sequence of</italic> C. albicans <italic>progressively in HIV patients to infer selective events and describe patterns of mutations, and those that survive and likely have an impact on resistance. Overall, the study is large and the data collected is very important and relevant with a degree of phenotypic assessment paired with the sequence data. The experimental and sequence analysis methodologies seem to be state of the art</italic>.</p><p>We thank the reviewers for these comments. As we note below, realizing that the data set could serve as a gold standard (and indeed is being used by several labs already), we decided to leverage the advances in sequencing yields, to generate much deeper sequencing and better analysis. Below we address each of the concerns raised.</p><p><italic>1) One major concern arises from consideration of all sections of the manuscript. There is very little statistical modeling to assess how likely these combinations of certain mutational patterns are due to real selective events due to specific phenotypic consequences or just by chance. Most of the manuscript is presented in this style, and it is hard to know which signals are truly significant. We provide two of the many examples in the manuscript to illustrate this point</italic>:</p><p><italic>In the Results section, the authors state: “We consider… two to three isolates”. While this is a reasonable argument, without proper statistics it is not possible to evaluate it. Perhaps some degree of linkage disequilibrium would be useful to consider</italic>.</p><p><italic>In the Results section, the authors also mention: “The recurring appearance of LOH events that coincide with changes in MIC suggests that they have been positively selected…” This is a “coincidence”, as the authors say, but is it significant?</italic></p><p>We thank the reviewers for this comment. We appreciate the suggestion for linkage or LD analysis, however it is important to emphasize that linkage/LD, by their very nature, are a function of mating and meiosis; and while mating can be engineered in the laboratory, meiosis has never been observed in <italic>C. albicans.</italic> Instead, <italic>C. albicans</italic> is thought to proliferate primarily by asexual reproduction and thought to only rarely undergo a cryptic, parasexual lifestyle in which two diploid individuals fuse into a tetraploid, and followed by chromosome loss in subsequent mitotic divisions to a converge on diploidy. This feature of <italic>C. albicans</italic> biology is a major limitation to genetic dissection of <italic>C. albicans</italic> population genetics either in the lab or in natural settings, and is important to bear in mind when considering our work.</p><p>It is for this reason that we think that those mutations that are simultaneously concomitant with increases in drug resistance and that are also recurrent are likely a footprint of selection (regardless of whether all were selected from existing diversity within the population, all arose de novo, or perhaps, a combination thereof). Without a substantially increased number of samples and non-treated human-passaged <italic>C. albicans</italic> longitudinal isolates, it is infeasible to adequately model this system for particular point mutations and genes affected. Furthermore, our samples were archived over decades, by clinicians who have typically archived only a single clone (sampled from what is likely a complex heterogeneous population). Our assessment of the larger gene sizes for those genes with recurrence in the original manuscript was an open acknowledgement of the possible contribution of random mutations to our results.</p><p>In the revised manuscript, we have now more clearly described the inherent challenges of this system and the limitations of our analysis, and we have attempted to clarify its statistical basis wherever possible. First, we now highlight these challenges early in the manuscript, clearly acknowledging the inherent limitations at deriving direct statistical significance driven by a model of the underlying phenomena. Second, whenever possible, we articulate the statistical basis for our analyses. For example, consider the LOH presented in <xref ref-type="fig" rid="fig3 fig4">Figures 3 and 4</xref>. There are 30 LOH events, if we limit LOH events to chromosomes arms (of which they occur on 15 out of 16 possible chromosome arms). Of these 30 LOHs, 23 meet our standard for persistence. The right arm of chromosome 3 has a persistent LOH 7 times. If we assume a naïve binary distribution, then the probability of this occurring is less than<inline-formula><mml:math id="inf1"><mml:mrow><mml:mo> </mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:mn>23</mml:mn></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mn>7</mml:mn></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mn>15</mml:mn></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mn>7</mml:mn></mml:msup><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mn>14</mml:mn></mml:mrow><mml:mrow><mml:mn>15</mml:mn></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mn>16</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> = 4.76 × 10<sup>-4</sup>. The same calculation applies to the left of arm 5 (Chr5L), but with 5 LOH events for a probability of 0.0128. Unlike chromosome 3, Chr5L is known to be a functionally significant LOH (<xref ref-type="bibr" rid="bib13">Coste et al., 2006</xref>; <xref ref-type="bibr" rid="bib72">Selmecki et. al, 2006</xref> and <xref ref-type="bibr" rid="bib74">2008</xref>). We have added this analysis to the text.</p><p>Furthermore, we now use stochastic simulations to generate null distributions of: (1) the expected number of persistent non-synonymous SNP-containing ORFs associated with changes in MIC; and (2) the degree of recurrence amongst these ORFs. These models (described in the revised Methods) suggest that the number of observed cases of each of (1) and (2) is significantly larger than expected by chance (p<0.0001). The distributions are shown in Author response image 1 below, and the results are reported in the revised text.<fig id="fig9" position="float"><label>Author response image 1.</label><caption><p>Null distributions of (a) the expected number of persistent, recurrent non synonymous SNPs outside of LOH regions, and (b) persistent non-synonymous SNPs associated with MIC changes outside of LOH regions.</p></caption><graphic xlink:href="elife00662f009"/></fig></p><p><italic>2) The last part of the Results describes some validation experiments where the specific mutations are brought to an experimental setting with deletion mutants. This is an experiment that was necessary but surprisingly it was not given the proper weight (it was only a small fraction of the total text), plus it was not as conclusive as one would like regarding the MIC phenotype</italic>.</p><p>We thank the reviewers for the suggestion to further emphasize these results . The diploid genome and lack of a complete sexual cycle make genetic manipulation of <italic>C. albicans</italic> challenging. Therefore, we relied on an existing resource of available homozygous deletion mutants from a deletion strain collection of 674 loci to gain insight into the contribution of these loci to drug resistance and fitness. In the original manuscript, we tested all those loci recurrently mutated by our original analysis for which a deletion mutant was available.</p><p>Following our deeper sequencing, we could reliably detect a larger number of recurrently mutated loci, and hence could test 23 deletion strains (revised <xref ref-type="fig" rid="fig8">Figure 8</xref>). In this expanded panel of tests, consistent with a role in host adaptation, 8 out of 23 mutants affected in vitro fitness in a culture medium thought to approximate in vivo conditions. Strikingly, three mutants show an increased level of fitness relative to the WT parental isolate, SN250. We have clarified the text to convey this point. The results reported in the original submission were from fitness experiments conducted using the drug-resistant reference strain (also used in the fitness experiments illustrated in revised <xref ref-type="fig" rid="fig7">Figure 7B</xref>). In principle, any reference strain can be used in relative fitness experiments. However, in this revision we used the <italic>same</italic> reference strain to quantify the fitness of the clinical isolates and the deletion mutants, such that that results in <xref ref-type="fig" rid="fig7">Figure 7A</xref> and <xref ref-type="fig" rid="fig8">Figure 8</xref> can be readily comparable by the reader.</p><p><italic>3) Considering the central role in the paper of the various classes of mutations called from the sequencing reads, I think that one technical issue in particular deserves further exploration and/or discussion. The SNP calling accuracy was found to be 87.3% on average (using a Sequenome iPLEX genotyping assay; Results, “Sequenome iPLEX genotyping assay”, and</italic> <xref ref-type="table" rid="tbl2"><italic>Table 2</italic></xref><italic>). There was however an important inter-sample variability, with very poor results observed in half of the tested series (<65% calling accuracy). This is worrying. Could this be due to intra-patient diversity, with different</italic> Candida <italic>subpopulations being sampled for the two assays? Is there any other potential explanation? If not, these inconsistencies should be listed as an important limitation of the study.</italic></p><p>The same genomic DNA from a single archived clone was used for both the original sequencing and the Sequenom assay. We were originally concerned that a switching in samples was the cause of the discrepancy, but repeated genotyping and re-examination of the Sequenome iPLEX versus the sequencing data suggested that when sequence coverage was relatively low, our GATK pipeline defaulted to the genotype of the provided reference genome resulting in erroneous SNPs calls across isolates in the same time course. Given that the extensive sequence data set is a major strength of our study, and with the substantial increase in sequencing yields since the time when we originally sequenced these isolates, we chose to re-sequence all the strains to high coverage (53-283X, average depth of coverage: 103X).</p><p>The new sequencing data resolved the discordance between the Sequenom and sequencing genotype calls (<xref ref-type="table" rid="tbl2">Table 2</xref>) such that our average concordance rose to 94% (85%-98%, <xref ref-type="fig" rid="fig1s1">Figure 1–figure supplement 1A</xref>). To further understand the nature of our discordant sites, and importantly, determine if improved filtering of the SNP calls could improve concordance, we examined the discordant sites for any potential base bias (<xref ref-type="fig" rid="fig1s1">Figure 1–figure supplement 1</xref>). From this analysis, it is apparent that the majority of discordant sites are those that are identified as homozygous by Sequenom iPlex (<xref ref-type="fig" rid="fig1s1">Figure 1–figure supplement 1B</xref>, teal bars), but as heterozygous by Illumina sequencing (<xref ref-type="fig" rid="fig1s1">Figure 1–figure supplement 1B</xref>, orange bars). We compared several quality metrics for these sites in our Illumina data for discordant vs. concordant sites (<xref ref-type="fig" rid="fig1s1">Figure 1–figure supplement 1C-G</xref>: depth of coverage (DP), mapping quality (MQ), PHRED scaled quality scores (QUAL), the QD score (confidence of a variant call), and the ratio of alternative allele to reference allele (AB)). We found that while the distributions of some quality scores for discordant sites are somewhat lower, the degree of overlap is so extensive that any improvement in specificity will come at a substantial sensitivity cost. Given that the discordance rate is already low, and that we focus on SNPs that both persist within a series and recur across series (and are thus called many independent times), we opted not to use any stricter thresholds for SNP calling. We present this comparison in the revised text, in the Results section.</p><p>Importantly, all previous conclusions were reproduced in this enhanced analysis, and we were able to draw additional conclusions based on these higher quality data. For example, the new set of 240 mutated genes is now enriched for several functions known to influence pathogenicity, including fungal-type cell wall (18 genes, p<0.0012) and cell surface genes (24 genes, p<0.00012), with several members in each of three cell wall gene families important for biofilm formation and virulence (<xref ref-type="bibr" rid="bib34">Hoyer et al., 2008</xref>): the Hyr/Iff proteins (HYR1 and 3, IFF8 and 6), the ALS adhesins (ALS1-4,7,9), and the PGA-30-like proteins (7 genes) (<xref ref-type="supplementary-material" rid="SD4-data">Figure 5–source data 1</xref>). All three families are specifically expanded in the genomes of pathogenic <italic>Candida</italic> species (<xref ref-type="bibr" rid="bib8">Butler et al., 2009</xref>). In addition, seven members of the FGR gene family, involved in filamentous growth that is specifically expanded in <italic>C. albicans</italic> (<xref ref-type="bibr" rid="bib8">Butler et al., 2009</xref>), are also among the 240 genes (<xref ref-type="supplementary-material" rid="SD4-data">Figure 5–source data 1</xref>).</p><p><italic>4) The paper is presented as a study of the within-host mutations that occur and are selected for in the</italic> Candida albicans <italic>genome between longitudinally sampled genomes before and after the introduction of antibiotic treatment. However, the data suggest that the variants that distinguish</italic> C. albicans <italic>genomes at different time points did not arise during the sampling window but were largely, maybe entirely, pre-existing (standing variation). This view is based on: (i) the very considerable genetic diversity (thousands of SNP differences) between isolates of the same strain within the same host, sampled less than a few months apart, in one case on the same day; (ii) the implausible mutation rate that would be needed to explain this variation as newly arising by mutation; and, (iii) the lack of analysis demonstrating that the data shows any signal of “measurable evolution”.</italic></p><p>We completely agree with the reviewers that, given the evidence of considerable genetic diversity existing in the same host, selection may primarily be driving an increase in the frequency of variants that already existed in the population. We have revised the text throughout the manuscript to clarify this point (especially in the Introduction and Results sections). Nevertheless, there is also ample evidence that mutation rates are increased under stressful conditions (e.g., drug treatment) in many systems, including yeasts, akin to what occurs in cancer (<xref ref-type="bibr" rid="bib26">Galhardo et al., 2007</xref>). Mitotic recombination events that give rise to LOH occur at a much higher rate (3 x 10<sup>-5</sup>) than single-nucleotide mutation rates and are also increased by stress. Specifically relevant here, LOH rates are elevated in <italic>C. albicans</italic> treated with fluconazole (<xref ref-type="bibr" rid="bib20">Forche et al., 2011</xref>). Thus, we cannot rule out the possibility that some genetic changes, including LOH events, occurred de novo during the treatment and sampling window. We discuss this point in the revised manuscript, in the Discussion.</p><p>Importantly, we are unable to distinguish pre-existing from de novo events with the data in hand. The rate of collection and number of samples we have for the majority of patients likely does not constitute a dataset that can be analyzed by existing approaches to ‘measurably evolving population’ (MEP, as per Drummond et al., 2003), and we do not want to overstate our evolutionary scenario. However, even if, as the reviewers suggest, <italic>“all ‘persistent mutations’ are those pre-existing variants driven to sufficiently high frequency by selection to ensure repeated sampling, whereas transiently observed variants are not”,</italic> this scenario is still an evolutionary process, involving heritable variation and selection, and can still inform our functional understanding of host-pathogen interactions.</p><p>With these limitations in mind, we were intrigued by the reviewers’ suggestion to assess the rate and time scale of evolution in these strains. While the majority of our patient series represent only a few sequentially isolated clonal samples, Patient 1’s series consists of 16 clonal, sequential isolates sampled over a period of nearly two years, and thus these samples can be analyzed for evidence of measurable evolution. We have used a Bayesian Markov Chain Monte-Carlo approach (Bayesian Evolutionary Analysis by Sampling Trees, BEAST; Drummond et al., 2012) to assess the rate of mutation in these populations and the likely time to most recent common ancestor (TMRCA).</p><p>To test the hypothesis that the majority of mutations occurred before the sampling period, we sampled trees generated using a set of parameters based on available data. Because there are no estimates of the per-basepair, per-generation mutation rate in <italic>C. albicans</italic>, we began with the assumption that <italic>C. albicans</italic> mutates at a rate similar to the per-base pair, per-generation rate of <italic>S. cerevisiae</italic> (Lynch et al., 2008). We then adjusted this rate to a per-day rate based on the generation time observed for C<italic>. albicans</italic> in a mouse model (1.56 generations per day, <xref ref-type="bibr" rid="bib23">Forche et al, 2009</xref>), and determined the TMRCA for the isolates from Patient 1. Based on these assumptions, the TMRCA for these 16 isolates is ∼ 235 years (Author response image 2A). Indeed, using these parameters, many of the nodes separating two isolates are dated to be ∼150 years in the past. While little is known of the genetic diversity present <italic>in C. albicans</italic> populations at the time of colonization, it seems unlikely that the most recent common ancestor for the clonal population seen in Patient 1 evolved 235 years prior to isolation: if colonization occurs primarily during birth, it would mean that much of the diversity present in this patient actually evolved in the patient’s grandmother, or great grandmother.</p><p>To better understand the diversity present in the population, we repeated this analysis, but with a broader estimate of mutation rate (for the prior, the same mean as above, 5.5 x 10<sup>-10</sup> mutations/bp/day, but with a larger standard deviation), resulting in a 5 to 55 fold increase in the estimated substitution rate, and a corresponding reduction in the likely TMRCA (Author response image 2B,C). To determine which of these three models best fit the data, we performed a pairwise calculation of Bayes factors (Author response image 2D), and we find that model (C), with an estimated substitution rate of 8.48 x 10<sup>-8</sup> best approximates our data. These results are consistent with the reviewer’s comment that considerable variation developed before the onset of treatment and sampling; however, they also suggest that a substantial portion of the observed variation occurred during the treatment period in Patient 1. Importantly, they also suggest that an elevated mutation rate likely explains some of the diversity we observe, consistent with previous observations of in vivo hypermutability during infection and treatment (Oliver et al, Science 2000, <xref ref-type="bibr" rid="bib26">Galhardo et al., 2007</xref>, <xref ref-type="bibr" rid="bib20">Forche et al., 2011</xref>).<fig id="fig10" position="float"><label>Author response image 2.</label><caption><p>The results of three separate MCMC analyses based on differing priors for mutation rate are shown. This analysis is based only on SNPs from isolates of Patient 1.</p></caption><graphic xlink:href="elife00662f010"/></fig></p><p>Further work is necessary to determine: (1) the extent of population diversity at any given time point, (2) the extent of the bottlenecks during colonization and drug treatment, and (3) the impact of the host immune response and antifungal therapy on the mutation rate <italic>Candida albicans</italic>. In light of these issues—which are relevant to a variety of pathogens—this analysis is particularly provocative. However, as many of the parameters necessary to make these estimates are poorly informed and our sampling limits the analysis to only one patient series, we do not feel this analysis merits inclusion in the manuscript.</p></body></sub-article></article> |