Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mBio
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Editor's Pick Research Article | Ecological and Evolutionary Science

The Genome of the Human Pathogen Candida albicans Is Shaped by Mutation and Cryptic Sexual Recombination

Joshua M. Wang, Richard J. Bennett, Matthew Z. Anderson
Kirsten Nielsen, Editor
Joshua M. Wang
aDepartment of Molecular Microbiology and Immunology, Brown University, Providence, Rhode Island, USA
bDepartment of Microbiology, The Ohio State University, Columbus, Ohio, USA
cDepartment of Microbial Infection and Immunity, The Ohio State University, Columbus, Ohio, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Richard J. Bennett
aDepartment of Molecular Microbiology and Immunology, Brown University, Providence, Rhode Island, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Richard J. Bennett
Matthew Z. Anderson
aDepartment of Molecular Microbiology and Immunology, Brown University, Providence, Rhode Island, USA
bDepartment of Microbiology, The Ohio State University, Columbus, Ohio, USA
cDepartment of Microbial Infection and Immunity, The Ohio State University, Columbus, Ohio, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Matthew Z. Anderson
Kirsten Nielsen
University of Minnesota Medical School
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mBio.01205-18
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

The opportunistic fungal pathogen Candida albicans lacks a conventional sexual program and is thought to evolve, at least primarily, through the clonal acquisition of genetic changes. Here, we performed an analysis of heterozygous diploid genomes from 21 clinical isolates to determine the natural evolutionary processes acting on the C. albicans genome. Mutation and recombination shaped the genomic landscape among the C. albicans isolates. Strain-specific single nucleotide polymorphisms (SNPs) and insertions/deletions (indels) clustered across the genome. Additionally, loss-of-heterozygosity (LOH) events contributed substantially to genotypic variation, with most long-tract LOH events extending to the ends of the chromosomes suggestive of repair via break-induced replication. Consistent with a model of inheritance by descent, most polymorphisms were shared between closely related strains. However, some isolates contained highly mosaic genomes consistent with strains having experienced interclade recombination during their evolutionary history. A detailed examination of mitochondrial genomes also revealed clear examples of interclade recombination among sequenced strains. These analyses therefore establish that both (para)sexual recombination and mitotic mutational processes drive evolution of this important pathogen. To further facilitate the study of C. albicans genomes, we also introduce an online platform, SNPMap, to examine SNP patterns in sequenced isolates.

IMPORTANCE Mutations introduce variation into the genome upon which selection can act. Defining the nature of these changes is critical for determining species evolution, as well as for understanding the genetic changes driving important cellular processes. The heterozygous diploid fungus Candida albicans is both a frequent commensal organism and a prevalent opportunistic pathogen. A prevailing theory is that C. albicans evolves primarily through the gradual buildup of mitotic mutations, and a pressing issue is whether sexual or parasexual processes also operate within natural populations. Here, we establish that the C. albicans genome evolves by a combination of localized mutation and both short-tract and long-tract loss-of-heterozygosity (LOH) events within the sequenced isolates. Mutations are more prevalent within noncoding and heterozygous regions and LOH increases towards chromosome ends. Furthermore, we provide evidence for genetic exchange between isolates, establishing that sexual or parasexual processes have contributed to the diversity of both nuclear and mitochondrial genomes.

INTRODUCTION

A wide variety of genetic events contribute to the evolution of eukaryotic genomes. In asexual cells, haploid genomes evolve via the accumulation of point mutations as well as undergo recombination events that drive DNA expansions/contractions (indels). Heterozygous (het) diploid genomes also have the capacity to undergo loss-of-heterozygosity (LOH) events, in which genetic information is lost from one of the two chromosome (Chr) homologs. In addition, both haploid and diploid genomes may experience large-scale chromosomal changes such as gross rearrangements or acquisition of supernumerary chromosomes or other forms of aneuploidy (1, 2).

Many eukaryotic species also generate genetic diversity via sexual reproduction. Here, recombination between individuals provides an efficient mechanism for producing diverse progeny. Sexual reproduction can therefore promote adaptation to new environments more rapidly than asexual propagation (3, 4). However, associated fitness costs are incurred due to the energetic requirements of sex and to the fact that only 50% of parental alleles are passed on to sexual progeny (5–7). Sex can also be detrimental by breaking up highly beneficial allelic combinations (5, 8). Facultative sexuality, the ability to alternate between sexual and asexual forms of reproduction, promotes a more flexible lifestyle that can accelerate adaptation in response to environmental pressures (4, 9).

Candida albicans is a prevalent opportunistic fungal pathogen that was long thought to represent an obligate asexual species (10, 11). However, mating of diploid cells has been described in the laboratory and produces tetraploid cells that then return to the diploid state via a parasexual process of concerted chromosome loss (CCL) (12–15). Efficient mating requires that C. albicans cells undergo a phenotypic transition from the sterile “white” state to the mating-competent “opaque” state (16). Conjugation of opaque cells can occur via heterothallic or homothallic mating (17), and recombination during CCL involves Spo11, a conserved “meiosis-specific” factor involved in DNA double-strand-break formation across diverse eukaryotes (13, 18).

Clinical isolates of C. albicans exhibit a largely clonal population structure despite the potential for recombination via parasexual reproduction (19, 20). Multilocus sequence typing (MLST) separates C. albicans isolates into 17 clades, although discrepancies between MLST haplotypes and individual mutations suggested that recombination may act to generate new allelic variants (21). Analysis of a limited number of haploid mitochondrial loci also revealed allelic mixtures that suggest that recombination may have occurred within C. albicans populations (22, 23). However, despite these observations, C. albicans is still often described as an asexual species that does not undergo mating or recombination in nature (24, 25). Prior studies focused on a subset of genomic loci and presented conflicting evidence regarding the role of recombination in shaping C. albicans evolution (19, 21–23, 26), which can now be addressed by a detailed analysis of full-genome sequences.

In this work, we examined natural evolution of the diploid genome in the 21 sequenced C. albicans isolates that represent different clades, different sites of infection, and different countries of origin (27, 28). We reveal that emergent single nucleotide polymorphisms (SNPs) and indels often cluster together in the genome and occur more frequently in heterozygous (het) regions than in homozygous (hom) regions. Moreover, we highlight that both nuclear and mitochondrial genomes appear highly admixed in some isolates, which is indicative of genetic recombination having occurred between clades. Overall, our results establish that the C. albicans genome is a highly dynamic landscape shaped both by local mutations and by large-scale rearrangements and that sexual or parasexual mating has made a significant contribution to genotypic variation.

RESULTS

The availability of whole-genome sequencing data for 21 diverse C. albicans isolates (27) provided an opportunity to determine how genetic diversity is generated between strains in nature. The C. albicans diploid genome is ∼14 Mb and consists of eight chromosomes harboring ∼6,100 genes (27, 29, 30). SNPs occur at a frequency of ∼0.3% (i.e., an average of 1 SNP every 330 bp) between chromosome homologs in the standard laboratory strain SC5314 (27, 31). Among the 21 isolates, we found that SNP frequencies ranged from 0.5% between closely related strains within clade I to 1.1% between strains from different clades (see Table S1 in the supplemental material). A previous phylogenetic reconstruction using 112,223 SNP positions found that most strains matched their previously assigned fingerprinting clades and MLST subtypes, with the exception of P94015, which clustered separately from other clade I strains (Fig. 1A) (27). The finding of strong bootstrap values across the constructed phylogeny of these strains supports the idea of a primarily clonal lifestyle in which most polymorphisms are consistent with inheritance by descent (Fig. 1B). Accordingly, SNPs and indels fit a nonrandom distribution across the 21 sequenced isolates χ2 ((SNPs; 20, N = 302641) = 83118, P < 2.2E–16, (indels; 20, N = 19581) = 13825, P < 2.2E–16) (see Fig. S1 in the supplemental material).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Distribution of polymorphisms among 21 clinical isolates of C. albicans. (A) The numbers of heterozygous (purple) and homozygous (orange) strain-specific SNPs and insertion/deletions (indels) are plotted for each isolate. Clade designations for each isolate are color coded (right), and bootstrap values are placed within the phylogenetic tree (left). (B) A cartoon representation depicts the two types of sequence variants among the set of 21 isolates: single nucleotide polymorphisms (SNPs) shared by multiple isolates and SNPs that are specific to individual strains. Variants harbored by multiple strains suggest origins in a common ancestor, whereas strain-specific polymorphisms likely arose specifically in individual strain backgrounds. (C) Relative frequencies of strain-specific SNPs (blue) and strain-specific indels (orange) summed across all strains were plotted across the genome using 5-kb sliding windows. (D) Number of strains that show LOH for 50-kb windows across the C. albicans genome aligned to their chromosomal position. (E) The SNP frequency for each of the contiguous 445 heterozygous (het) and 468 homozygous (hom) regions across the 21 sequenced genomes was plotted. Normalized SNP density was calculated using the number of SNPs within a het or hom region divided by the region length and was significantly elevated for het regions (red) relative to hom regions (blue).

TABLE S1

Percentage of the genome harboring SNPs among the sequenced C. albicans isolates. Download Table S1, PDF file, 1 MB.
Copyright © 2018 Wang et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S1

The distribution of polymorphisms is nonrandom. The number of strains harboring each SNP (A) or indel (B) was determined and the frequency plotted. A best-fit line (red) was plotted for each distribution. Download FIG S1, PDF file, 0.8 MB.
Copyright © 2018 Wang et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Distribution of strain-specific polymorphisms across the C. albicans genome.LOH events can distort the patterns of SNPs inherited from ancestral strains (Fig. S2). To help limit these confounding effects, we restricted most analyses to “strain-specific” SNPs and indels that are unique to individual strains in our collection (Fig. 1B). Approximately 25% of all SNP positions and 10% of all indel positions were strain specific (66,086 and 6,474 events, respectively; Tables S2 and S3).

TABLE S2

Strain-specific SNPs among sequenced C. albicans isolates. Download Table S2, XLSX file, 1.9 MB.
Copyright © 2018 Wang et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S3

Strain-specific indels among sequenced C. albicans isolates. Download Table S3, CSV file, 0.2 MB.
Copyright © 2018 Wang et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S2

LOH can obfuscate patterns of inheritance by descent of shared polymorphisms. LOH of opposing alleles in a common ancestor can make it appear that mutations arose independently despite the mutations sharing a common origin. LOH of homolog A in clade III produces a SNP pattern different from that of clade II although t arose from the same ancestral strain. Download FIG S2, PDF file, 0.8 MB.
Copyright © 2018 Wang et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Analysis of the genome-wide distribution of strain-specific SNPs and indels across the 21 genomes revealed that these mutation types showed significant clustering with one another (Pearson, t = 11.64, df = 286, P  < 2.2E−16; Fig. 1C). Multiple SNPs often occurred within 100 bp of an indel (Fig. S3A), as was confirmed via Sanger sequencing of selected regions (see Fig. S4 at https://drive.google.com/drive/u/0/folders/1u0LiWcB-yufoL4iQOn_bQQ2ercxB-F6R). Enrichment of SNPs was observed immediately adjacent to indels (within 10 bp) but not within indels [Wilcoxon test, W(1.79E7), P  < 2.2E−16] (Fig. S3B). Three strains, P60002, P75010, and P94015, carried genes that encoded a large proportion of strain-specific mutations reflective of their longer branch lengths in the phylogenetic tree (Fig. 1A), which could potentially skew the analysis. However, even after removing these three strains from the analysis and reducing the four major clades to three representative strains each, we still observed a significant association between SNPs and indels [Wilcoxon test, W(3.04E7), P  < 2.2E-16] (see Fig. S5A at https://drive.google.com/drive/u/0/folders/1u0LiWcB-yufoL4iQOn_bQQ2ercxB-F6R). This association highlights that these mutagenic events often occur in close proximity to one another, as has been observed in other eukaryotic species (32), and that indel formation or the associated DNA repair processes may be mutagenic in C. albicans.

FIG S3

Unique SNPs and unique indels. Download FIG S3, PDF file, 1.2 MB.
Copyright © 2018 Wang et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Association between LOH recombination events and base substitution mutations.We next examined the global distribution of LOH tracts across the 21 genomes. LOH tracts were defined as contiguous 5-kb windows with a high frequency of homozygous SNPs (>0.4 events/kb; see Materials and Methods and reference 27). Plotting the incidence of LOH for all chromosomes revealed a clear pattern whereby the prevalence of LOH increased along each chromosome arm in progressing from centromere to telomere (Fig. 1D; see also Fig. S6 at https://drive.google.com/drive/u/0/folders/1u0LiWcB-yufoL4iQOn_bQQ2ercxB-F6R). In fact, the overwhelming majority of all long-tract LOH regions (155 of 170 regions larger than 50 kb) extended to the ends of the corresponding chromosomes (Table S4), indicating that most large LOH events initiate on one arm and then proceed to the end of the chromosome.

TABLE S4

Heterozygous and homozygous regions for the sequenced C. albicans isolates. Download Table S4, XLSX file, 0.3 MB.
Copyright © 2018 Wang et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Several studies have shown that the underlying genomic context can impact mutation rates, including the observation that mutation rates were higher in heterozygotes than in homozygotes during meiosis (33). We therefore examined mutational patterns in the 21 C. albicans genomes that consist of a mosaic of heterozygous and homozygous regions. We subdivided C. albicans genomes into heterozygous (het) and homozygous (hom) regions using defined criteria for all SNPs (see Materials and Methods), resulting in 468 het and 445 hom regions, respectively (Table S4). The het regions covered a total of 71.1% of the genome and hom regions 28.9%; het tracts were therefore considerably longer on average than hom tracts (∼480,000 bp versus 186,000 bp, respectively). Definition of het and hom regions using all SNPs allowed examination of the frequencies of strain-specific SNPs within these regions. Strain-specific SNPs comprised only 3% of all SNPs within these genomes and therefore did not contribute substantially to the designation of het and hom regions. Using these designations, het regions contained significantly higher densities of strain-specific SNPs than hom regions (1.4E−4 SNPs/bp versus 7.3E−5 SNPs/bp, respectively; Brunner-Munzel [BM] test value = −10.6, df = 786.6, P  < 2E−16; Fig. 1E). Again, to determine if the three outlier strains, P60002, P75010, and P94015, biased this analysis, we removed these three strains and reduced the four major clades to three representative strains each, and yet there were still significantly higher densities of strain-specific SNPs within het regions than within hom regions (BM test = −7.558, df = 377.0, P = 3.14E−13). Furthermore, all 21 isolates exhibited the same bias toward het regions containing more strain-specific SNPs than hom regions (two-tailed BM test = −1.11, df = 38.5, P = 0.28). This indicates that mutations preferentially accumulate in het regions rather than hom regions during natural evolution of C. albicans isolates, as was seen for het regions versus hom regions during meiotic propagation (33).

Identity by descent during C. albicans evolution.A hallmark of phylogenetic reconstructions in asexual species is the ability to track the relatedness of isolates based on inherited polymorphisms (21, 34, 35). Reconstruction often relies on a maximum parsimony model of “identity by descent,” in which more closely related strains share a greater percentage of shared polymorphisms (Fig. 2A). C. albicans SNP patterns generally follow this model as evidenced by strong bootstrap support at almost all nodes of the phylogenetic tree and visual examination of SNP patterns (Fig. 2B and C; see also Fig. S7 at https://drive.google.com/drive/u/0/folders/1u0LiWcB-yufoL4iQOn_bQQ2ercxB-F6R). Strikingly, however, certain regions of the genome exhibited clear violations of identity by descent (Fig. 2D and E; see also Fig. S8 at https://drive.google.com/drive/u/0/folders/1u0LiWcB-yufoL4iQOn_bQQ2ercxB-F6R). In heterozygous diploid genomes, these deviations could potentially arise through two mechanisms: (i) through sexual recombination between genetically distinct isolates or (ii) through multiple independent LOH events that obfuscate the actual pattern of descent (see Fig. S9 at https://drive.google.com/drive/u/0/folders/1u0LiWcB-yufoL4iQOn_bQQ2ercxB-F6R). In the latter case, multiple LOH events could cause loss or retention of SNPs through homozygosis of one chromosome homolog or the other, thereby generating a subset of isolates that appeared “recombinant,” i.e., appeared to have intermixed genetic content from two different relatives. Such a history can sometimes be inferred by a comparison of SNP patterns within the region of interest in multiple extant strains (Fig. 2A).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Mutational patterns following identity by descent and loss of heterozygosity are widespread. (A) Two patterns of polymorphisms exist within sequenced genomes. One pattern (depicted on the left) shows polymorphisms that are phylogenetically congruent (indicative of identity by descent), and the other (depicted on the right) shows polymorphisms that violate the phylogenetic relationship of the strains based on all variant positions (indicative of mechanisms other than direct inheritance). Shading denotes similarity in overall SNP patterns grouped by color. Heterozygous SNPs are purple, and homozygous SNPs are orange. (B and C) Polymorphisms for two loci (orf19.4822 [B] and orf19.4055 [C]) that display SNP patterns consistent with inheritance by descent across the sequenced C. albicans genomes are plotted. (D) The polymorphism pattern is plotted for a locus that does not follow inheritance by descent. Identical genotypes are color coded and connected to each other. (E) A cartoon depicting LOH of heterozygous positions in opposing directions provides the most parsimonious explanation for the observed SNP patterns across the C. albicans phylogeny.

We note that certain regions in the C. albicans genome are consistent with divergent short-tract LOH events having occurred during evolution (Fig. 2D and E; see also Fig. S8 at https://drive.google.com/drive/u/0/folders/1u0LiWcB-yufoL4iQOn_bQQ2ercxB-F6R and the supplemental material). In line with this, we identified 514 nonoverlapping 25-kb windows across the 21 sequenced isolates that do not harbor similar polymorphisms to closely related strains and that could represent regions that had experienced LOH. Interestingly, two strains, P60002 and P94015, contained 390 (75.9%) of these regions, although only 214 (54.9%) of the 390 incongruent regions in these two strains overlapped with LOH tracts (hom regions) in these isolates. In contrast, the majority (117 of 124) of the incongruent regions in all other strains overlapped with LOH tracts. This suggests that incongruence in polymorphisms in most strains likely results from divergent LOH events but that LOH cannot easily explain the majority of incongruent polymorphic patterns observed in P60002 and P94015.

Evidence for recombination in natural isolates of C. albicans.Previous studies provided conflicting messages regarding recombination in natural populations of C. albicans (19, 21–23, 26), and none examined whole-genome data for evidence of interclade mixing. We therefore examined the 21 sequenced genomes for mixed evolutionary histories. The similarity of genomic segments from each strain to the overall phylogenetic tree was determined by analysis of SNP patterns using 25-kb sliding windows. To aid visualization of SNP patterns, we developed a custom tool, SNPMap (http://snpmap.asc.ohio-state.edu/), which allows users to map the positions of individual mutations, mutation types, and het/hom tracts across user-defined regions of the 21 genomes.

Our analysis focused on P60002 and P94015, two isolates that had the weakest bootstrap support within the C. albicans phylogeny and that clustered with different strains by MLST or DNA fingerprinting analysis (27). The latter data could be consistent with their genomes being formed by recombination between isolates. In support of this, examination of Chr4 in P94015 identified one region with clear homology to clade I in close proximity to a region homologous to clade SA (Fig. 3A; see also Table S5). The region with homology to clade I (labeled P94015-A in Fig. 3A) shares a large number of SNPs that are present throughout clade I but absent in all other strains with the exception of P94015 (Table S5). Next to this region, a 1-kb segment (region P94015-B) was seen that lacks clear homology to any of the other sequenced isolates while, adjacent to this, the SNP pattern in P94015 was virtually identical to those of two clade SA isolates (region P94015-C). We similarly identified a region on ChrR in P60002 that switched homology between clades I and SA (Fig. 3B). A ∼4-kb stretch of the genome (region P60002-B; ChrR:1166000.1170000) harbored heterozygous SNPs that perfectly matched clade I SNPs and was flanked by segments that matched the SNP pattern among clade SA isolates (region P60002-A,C). These regions could therefore have been shaped by recombination between isolates.

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Evidence for recombination in C. albicans isolates. (A) All SNPs are shown for a 20-kb region of Chr4 for the 21 sequenced genomes. For each strain, dark gray bars indicate heterozygous genomic regions and light blue bars indicate regions that are mostly homozygous. SNPs that are heterozygous and homozygous compared to SC5314 are purple and orange, respectively. The SNP pattern in P94015 is highlighted to indicate one region with homology to clade I (P94015-A; outlined in red) next to a region without clear homology to any specific clade (P94015-B) followed by a region with homology to clade SA (P94015-C; outlined in gray). Zoomed segments from P94015-A and P94015-C showing identity to individual SNPs from clades I and SA, respectively, are displayed at the bottom of the figure panels. Clade designations for each isolate are color coded. All SNPs that are at the same position involved the same base substitutions. (B) The SNP patterns for an ∼20-kb region of ChrR highlights recombinant SNP patterns for P60002 and is labeled as described for panel A. The DNA segments corresponding to different clades are aligned next to the appropriate group and labeled according to homologous tracts. (C) SNPs for a 2.5-kb region of Chr3 were phased by comparison of clade SA and III strains that were heterozygous or homozygous for this region. Homozygous SNPs are connected across homologs by a black line. The SNP pattern of this region in P94015 matches the product of hybridization between the panel A homologs from clade SA and clade III (both the position of the SNP and the actual base substitution matched). Cartoons of phased SNP homologs are represented above the collective SNP pattern for each strain. Heterozygous and homozygous SNPs are purple and orange in the genome snapshots underneath the phased homolog cartoons, respectively. (D) Consensus SNP patterns for each clade were determined and used to assess similarities of 25-kb sliding windows in each isolate to the clade consensus patterns. The closest match for each window was color coded by clade. Two strains, P60002 and P94015, contained a large portion of regions assigned to multiple clades compared to all other strains. Dark gray, SA; blue, clade III; mustard, clade II; red, clade I; light gray, no clade consensus. (E) Each strain was compared to a consensus SNP pattern for the major represented clades using four K-clusters. The fraction of the genome corresponding to each of the four clusters (dark gray, SA; blue, clade III; mustard, clade II; red, clade I) is represented within the bar height for each strain. Light gray denotes an outlier clade.

TABLE S5

SNPs across the highlighted recombinant regions for P94015. Download Table S5, XLSX file, 0.01 MB.
Copyright © 2018 Wang et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S6

Genes with significant enrichment of strain-specific SNPs. Lists of genes at three different levels of SNPs per nucleotide cutoff are provided. Download Table S6, XLSX file, 0.02 MB.
Copyright © 2018 Wang et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Mating between isolates from different clades would be expected to generate hybrid DNA regions, with SNPs on one homolog of the recombinant strain matching those in one clade and SNPs on the other homolog matching those in a second clade. Identifying inherited SNPs following C. albicans mating in nature is complicated by the fact that, with the exception of SC5314 (31), haplotypes are not available for the 21 C. albicans genomes. Despite this, phasing of heterozygous SNPs for some isolates can be inferred using SNP patterns from related strains that have undergone LOH for that region (Fig. 3C). The region that experienced LOH would retain only the SNPs that reside on the same homolog (i.e., those that are phased). Using this approach, we phased SNPs for chromosomal regions that are heterozygous in some isolates but homozygous in closely related strains. The existence of multiple isolates that have undergone LOH for either of the two homologs strengthens the confidence of phasing assignments within a given clade.

We applied this approach to a region on Chr3 in P94105 that contains polymorphisms identical both to those on homolog A of a clade SA strain (12 of 12 SNPs are identical) and to those on homolog A from a clade III isolate (15 of 15 SNPs are identical) (Fig. 3C). Both the SNP positions and the nucleotide identities are conserved across this region in P94105 compared to the corresponding homologs from clade SA and clade III isolates. The identification of this hybrid region therefore provides a striking example of P94105 inheriting one homolog from a clade SA strain and one homolog from a clade III strain and establishes a nonclonal origin for this isolate.

To examine the global patterns of admixing among the set of 21 isolates, the distribution of all variant positions in each strain was compared to the consensus pattern for each clade using sliding 25-kb windows. The SNP patterns of most isolates resembled the consensus pattern for their assigned clade (98.5% of genomic windows matched their assigned clade), as expected for a population propagating clonally (Fig. 3D). In contrast, many regions within the P60002 and P94015 genomes showed homology to multiple different clades, producing highly mosaic genomes (Fig. 3D). Here, the genomes of P60002 and P94015 matched their assigned clades for only 58.3% and 76.7% of sliding windows, respectively (P = 1.14E−10). The majority of genomic regions in P94015 aligned with clade I (the genome is mostly color coded red in Fig. 3D), whereas numerous segments aligned to regions from three other major clades (SA, II, and III). In the case of P60002, clade SA regions made up the majority of the genome, with a smaller number of regions matching clade I or, to a lesser extent, clade II. In line with this, the branch point leading to P60002 is the least well-supported node in the phylogenetic reconstruction of all 21 isolates (27).

C. albicans strains containing mosaic genomes should display evidence of classical admixture and allow an estimation of clade contributions to the mosaicism. To assess admixture, we restricted analysis to the four major clades represented here and assigned the fraction of SNPs within each genome to each of the clades. Consistent with the previous sliding approach to assign sequence homology, P60002 and P94015 were identified as the only strains with significant levels of admixture (Fig. 3E). Both strains displayed genetic contributions from all four major clades, suggesting that relatively ancient recombination events generated these admixed genomes. The most parsimonious explanation for these highly mosaic genomes is that they are the products of mating and recombination between isolates from multiple C. albicans clades.

Analysis of mitochondrial genomes in C. albicans isolates.Haploid mitochondrial genomes provide a more simplified context for the search for evidence of recombination than heterozygous diploid genomes. In Saccharomyces cerevisiae, mitochondrial genomes are biparentally inherited following mating, and recombination between parental genomes can occur prior to zygote division (36). We therefore performed the first comparative analysis of global SNP patterns in C. albicans mitochondrial genomes using sequencing data from the set of 21 isolates. The mitochondrial genome in C. albicans is ∼41 kb in size, and a total of 1,847 SNPs (and 0 indels) were annotated within the 21 isolates, with an average SNP density of 1 polymorphism every 476 bp. The mitochondrial genome was highly heterogeneous and included areas of high SNP density (e.g., ChrM positions 15000 to 20000) and regions devoid of polymorphisms (e.g., ChrM: 6000–12000). Furthermore, of the 1,847 annotated mitochondrial SNPs, only 39 were strain specific in the set of 21 sequenced genomes (Table S2).

C. albicans mitochondrial genomes generally showed clade-specific SNP patterns that were again consistent with a clonal population structure, although the resolution of the SNP patterns was low due to relatively few clade-defining mitochondrial SNPs (Fig. 4). As with nuclear genomes, examination of mitochondrial genomes of P60002 and P94015 showed evidence of interclade mixing. For example, the mitochondrial genome of P94015 contained regions that aligned with mitochondrial segments from both clade SA and clade II (Fig. 4A). There, three polymorphisms were clade SA specific on ChrM: 1–6000 (region P94015-A; black dots), and all three are present in P94015 (Fig. 4A). An additional 2 of the 15 polymorphisms in region P94015-A are specific to this strain (Table S2). The remainder of the mitochondrial genome in P94015 (region P94015-B) harbors 144 polymorphisms matching the SNP pattern found in clade II (with the exception of two strain-specific SNPs). Recombination between clades I and II was even clearer in the mitochondrial genome of P60002, as the majority of this genome was identical to that of clade I (P60002 regions A and C), but a 4-kb segment (ChrM: 19000–22000) harbored 34 polymorphisms that were identical to clade II polymorphisms and were entirely absent in clade I (P60002 region B) (Fig. 4B; see also Fig. S10 at https://drive.google.com/drive/u/0/folders/1u0LiWcB-yufoL4iQOn_bQQ2ercxB-F6R). We also found that maximum parsimony approaches mostly separated the P90145 mitochondrial genome into clade II and SA regions and the P60002 mitochondrial genome into clade I and II regions (Fig. 4D), consistent with visual alignments (Fig. 4A and B). Results of direct Sanger sequencing of the mitochondrial genome (regions 5000 to 5700, 18500 to 19500, and 29000 to 30000) supported the SNP designations from whole-genome sequencing and therefore established that interclade recombination has occurred within the mitochondrial genomes of P60002 and P90145 (Fig. 4E; see also Fig. S11 at https://drive.google.com/drive/u/0/folders/1u0LiWcB-yufoL4iQOn_bQQ2ercxB-F6R).

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Mitochondrial genomes in C. albicans display recombinant genotypes. The mitochondrial (mt) genome sequences of 21 clinical isolates of C. albicans were compared. The positions of SNPs that differ from the SC5314 assembly are shown (excluding ChrM: 35000–41000 due to the absence of any SNPs in this region). (A to C) The mt genomes for P94015 (A), P60002 (B), and P75010 (C) are highlighted to show the mosaic configuration of SNPs relative to other clades for these strains. The mt genome in P94015 contains regions that share polymorphisms that are identical with those of clade SA (P94015-A) and clade II (P94015-B). A key region matching clade SA is marked with an asterisk, with the three clade SA-specific SNPs marked in black. In contrast, the P60002 mt genome aligns with sequences for clades I (P60002-A and P60002-C) and II (P60002-B and P60002-D). The mt genome of P75010 is homologous to clades II (P75010-A), I (P75010-B), and SA (P75010-C). Key clade II-specific SNPs in the P75010 alignment are marked in black. A 6-kb region devoid of SNPs is more lightly shaded. (D) The similarity of the mt genome from each isolate to consensus SNP patterns was determined for each clade in 2-kb sliding windows. The window was color coded to designate the clade with greatest similarity. Dark gray, SA; blue, clade III; mustard, clade II; red, clade I; light gray, no clade consensus. (E) Two strains from clade SA and clade I, along with P60002 and P94015, were subjected to Sanger sequencing across three separate 1-kb regions of their mitochondrial genomes. Chromatograms highlight variant positions in P60002 and P94015 consistent with recombination between clades as determined by genome analysis.

We further note that the single representative of clade E in our collection, strain P75010, also displays strong evidence of recombination in its mitochondrial genome (Fig. 4C). The first ∼12 kb of P75010 (region P75010-A) aligns closely with clade II and harbors all three clade II-specific SNPs in this region (region P75010-A; black dots). The ∼7-kb region of P75010 (region P75010-B) that follows matches that of clade I strains and is followed by a region with clear homology to clade SA, harboring 22 of 25 clade SA-specific SNPs. Identity to several clades implies that multiple recombination events gave rise to this complex SNP pattern. Taken together, these results reveal that recombination events have occurred within C. albicans mitochondrial genomes and provide clear evidence that sexual/parasexual processes have occurred during C. albicans evolution.

DISCUSSION

Our analysis of the C. albicans genome structure reveals a number of important aspects concerning mutational patterns occurring during the natural evolution of the species. We highlight that (i) there is a significant association between the positions of emergent SNPs and indels; (ii) diverse LOH events contribute to genetic inheritance, including long-tract LOH events that extend to the ends of the chromosomes; (iii) heterozygous regions of the genome accumulate more mutations than coding and homozygous regions; (iv) a subset of strains exhibit mosaic nuclear and mitochondrial genomes; and (v) analysis of specific chromosomal regions reveals clear evidence for interclade recombination.

Mutations driving natural evolution of the C. albicans genome.Mutation rates differ across eukaryotic genomes in a context-dependent manner (37, 38). We found that mutation patterns arising in natural C. albicans populations similarly exhibit a nonrandom distribution across the genome. Our analysis focused on the distribution of strain-specific SNPs, as these SNPs are likely to have emerged since these strains diverged from one another. We found that the location of strain-specific SNPs was biased toward heterozygous rather than homozygous regions of the genome. This is consistent with recent studies that showed that mutations arising during meiosis occurred more frequently within heterozygous regions of eukaryotic genomes than within homozygous regions, although the mechanism driving this bias is unknown (39, 40). Our results extend these findings by indicating preferential accumulation of mutations within heterozygous regions during mitotic growth in C. albicans, suggesting that elevated mutation rates may be a common feature associated with heterozygosity.

Our analysis also reveals that emergent SNPs and indels cluster together within the C. albicans genome, with a significant enrichment of SNPs within 10 bp of emergent indels (Fig S3). This is similar to what has been observed in other eukaryotic species where indels were found to promote elevated substitution rates close to the founding indel (32, 41–44).

Evolutionary impact of loss of heterozygosity.Loss-of-heterozygosity events are a frequent occurrence in C. albicans genomes (1, 45). Previous analysis of the 21 C. albicans genomes noted that LOH tracts can be highly variable in size, with several isolates having experienced whole-chromosome LOH (27). Here, we made the observation that the vast majority of very large LOH tracts (defined as segmental LOH tracts of >50 kb) initiated between the centromere and the telomere and extended to the chromosome ends. Such LOH events are likely due to break-induced replication, in which one chromosome homolog is used as a template to repair a double-strand DNA break in the other homolog, although reciprocal crossover events can also generate long LOH tracts (25, 46, 47). The breakpoints for most of these long-tract LOH events differed between individual strains, suggesting that independent LOH events had occurred. Short-tract LOH (<50 kb) was also common, with approximately half of these events being shared between related strains (see Table S4 in the supplemental material; see also Fig. S6 at https://drive.google.com/drive/u/0/folders/1u0LiWcB-yufoL4iQOn_bQQ2ercxB-F6R). Homozygosis of certain regions may confer a selective advantage given that LOH has been shown to alter C. albicans traits, from growth rates to virulence to drug resistance (28, 48–50).

The cumulative effects of all of these mutational forces on C. albicans genomes result in accelerated evolution of heterozygous regions relative to homozygous regions. Eventually, the rapid accumulation of deleterious mutations during clonal growth would be expected to result in a fitness decline due to Muller’s ratchet (51, 52). The occurrence of LOH may counterbalance these forces both by culling mutations from the genome and by reducing heterozygosity to low evolutionary rates. LOH events appear at different points during the evolution of individual strains, with a number of LOH tracts differing even between closely related strains. Indeed, the accumulation of strain-specific SNPs within more-ancestral LOH tracts demonstrates that these LOH events are not recent occurrences and additional mutations have emerged since their origin.

Evidence for genetic exchange in clinical C. albicans isolates.Studies of C. albicans population structure point to a largely clonal mode of reproduction, and yet there is also evidence of mixed evolutionary histories indicative of a species reproducing sexually/parasexually (20, 21, 23, 26, 53–55). Furthermore, genes harbored at the mating locus show evidence for ongoing selection consistent with a conserved role in regulating sexual/parasexual reproduction (56). In this study, we interrogated whole-genome data for evidence of genetic admixture, while noting that LOH events can complicate analysis of recombination (Fig. 2; see also Fig. S8 at https://drive.google.com/drive/u/0/folders/1u0LiWcB-yufoL4iQOn_bQQ2ercxB-F6R). We reveal that a subset of isolates contained mosaic genomes, consistent with those genomes being the products of mating between different C. albicans clades. Both nuclear and mitochondrial genomes of P60002 and P94015 show recombinant genotypes, supporting the idea of a (para)sexual origin for these strains. Genetic information from multiple clades contributed to these genomes, and recombination tracts ranged in length from a few kilobases to hundreds of kilobases (Fig. 3 and 4). The existence of a subset of C. albicans strains with mosaic genomes is similar to what has been observed in wild and domesticated strains of S. cerevisiae, where both nonmosaic and recombinant mosaic genomes have been identified (57). Analysis of admixed genomes in P60002 and P94015 suggests that recombination events may be relatively ancient, as recombination involves multiple clades and more-recent mutational events have obscured the precise evolutionary histories of these strains.

Disagreement between strain phylogeny data determined by MLST and Ca3 fingerprinting may also be indicative of recombination in the population. On the basis of Ca3 fingerprinting, P94015 should cluster with other clade I strains and is supported by MLST analysis, which groups P94015 and other MLST 6 strains closer to MLST 1/clade I than other groups (27). And yet, whole-genome sequencing clusters P94015 with clade SA rather than with other clade I strains. This could reflect how recombination has distorted the phylogenetic relationship between strains on the basis of analysis of a small subset of loci. Analysis of additional C. albicans isolates will help define the prevalence of recombination in the species and whether these recombination events are ancient or are more recent occurrences.

Critically, we identified regions in C. albicans genomes that exactly matched the pattern of recombinant SNPs expected from mating events occurring between two extant clades. This was exemplified by one region in P94015, which consisted of multiple SNPs that exactly matched those present in clade I and that were followed, after a short gap, by a run of SNPs that precisely matched those in clade SA (Fig. 3A). Moreover, recombination events were clearly evident in the mitochondrial genomes of at least 3 isolates (P60002, P75010, and P94015) among the 21 examined. Taken together, these studies provide the clearest evidence to date that C. albicans populations have been shaped by (para)sexual exchange.

In summary, C. albicans genomes reveal multiple signatures of the forces that have shaped genetic diversity within the species. Both short and long LOH events have played a major role in increasing population diversity, with large tracts extending to chromosome termini that impact hundreds to thousands of polymorphisms. Base substitution mutations and indels cluster within heterozygous regions of the genome, suggestive of faster evolution of these regions. Finally, recombination between isolates has generated mosaic nuclear and mitochondrial genomes, which could have enabled adaptation in the species. The diploid heterozygous genome of C. albicans is therefore a highly dynamic platform on which selection can act.

MATERIALS AND METHODS

Variant calling, processing, and display.Whole-genome sequencing, variant identification, and loss-of-heterozygosity (LOH) windows were identified in a previous study (27). Briefly, BWA 0.5.9 (58) read alignments were filtered with a minimum mapping quality of 30 using SAMtools (59). To reduce the incidence of false-positive SNPs called near indels, poorly aligned regions were realigned using GATK RealignerTargetCreator and IndelRealigner (GATK version 1.4-14) (60). Both prior SNP variants and variants present in the mitochondrial genomes and indels for both the mitochondrial and nuclear genomes described here used GATK UnifiedGenotyper, and results were filtered using GATK VariantFiltration and hard filters (QD, <2.0; MQ, <40.0; FS,  >60.0; MQRankSum, less than 12.5; ReadPosRankSum, less than −8.0). SNPs were called homozygous if greater than 90% of the reads contained the nonreference nucleotide. This high threshold also reduced miscalling due to trisomic chromosomes/chromosomal regions.

To assess heterozygous and homozygous variant calls, the number of reads at each variant position was divided by the total number of reads at that position in that strain. On average, each variant position had 52.17 reads, with an interquartile range of 33 to 71 reads. The distribution of the allelic ratio at each variant is plotted in Fig. S12 at https://drive.google.com/drive/u/0/folders/1u0LiWcB-yufoL4iQOn_bQQ2ercxB-F6R. Dividing the mean number of reads of the allele by the total number of reads resulted in a value of 0.4499 for heterozygous positions and a value of 0.9889 for homozygous regions. Thus, a 90% threshold was used for homozygous regions whereas heterozygous regions span from 10 to 50%.

From this data set, strain-specific SNPs and indels were parsed into a separate set for additional analysis. Strain-specific variant features were required to be uniquely called in only 1 of the 21 strains at its genomic position. The data sets are available online in a searchable interface, using R shiny for the backend (https://snpmap.asc.ohio-state.edu/). Manual interrogation of 100 variants using the Integrated Genome Viewer (61) confirmed the variant presence and call quality metrics in all 100. Manual interrogation of SNP-indel pairs in the pileup confirmed 98% (49 of 50) are valid using the same criteria.

To adjust for mutations in homozygous regions that might have occurred but were then lost due to LOH, the following calculation was used for each homozygous region: {[(homozygous SNPs/all SNPs) + 1] * homozygous SNPs}/length. This effectively doubles the number of homozygous SNPs in hom regions to account for heterozygous SNPs that are lost due to LOH. Even with this doubling, significantly more emergent SNPs occurred in het regions than in hom regions.

Phylogenetic construction of strain relatedness.The methods used to construct the phylogeny of the sequenced strains were previously described by Hirakawa et al. (27). Briefly, a representation of the phylogenetic relatedness of the strains was constructed using all whole-genome SNP calls, which totaled 113,339. A distance-based tree was estimated relying on maximum parsimony and a stepwise matrix where homozygous positions are two steps away compared to one step for heterozygous positions. SNP positions were resampled in 1,000 bootstrapped samples, and each node indicates the bootstrap support.

Determination of LOH.Previously defined heterozygous and homozygous genomic regions were used as detailed previously (27). Briefly, the SNP density of homozygous and heterozygous positions was calculated across the genome in nonoverlapping 5-kb windows for each isolate. The resulting 5-kb windows were managed as follows: (i) homozygous regions shared between individual strains and SC5314 were identified and marked as homozygous; (ii) a single 5-kb window adjacent to these regions lacking any polymorphisms was merged into homozygous tracts (if present); (iii) contiguous, adjacent windows with a significantly higher frequency of homozygous SNPs than SC5314 homozygous regions (>0.4 SNPs/kb) were merged, allowing one intervening window lacking sufficient polymorphism into homozygous tracts. These regions were defined as homozygous whereas the remaining regions covered by 5-kb windows of the genome were designated heterozygous and contained significantly more heterozygous SNPs. The borders between homozygous and heterozygous regions were manually inspected for accuracy.

Introgression/tree violations.Two independent methods were used to assign clade designations for nuclear and mitochondrial genomes in nonoverlapping 25,000-bp or 750-bp windows, respectively. The first process assigned the clade that most closely corresponded to the SNP pattern of the query strain. However, clade I contains fewer SNP calls due to alignment to the SC5314 reference genome (also a clade I isolate), which can introduce an artificial bias. Therefore, a data frame was constructed where each row represented any SNP contained within any strain in the query window. A SNP could be counted only as a single row, so the identity of that SNP position was recorded as “0” for the absence of the SNP, “1” for a heterozygous position, and “2” for a homozygous position. The correlations between the target strain’s resulting numeric SNP profile and that of each of the other strains were individually calculated. Scores from strains within the same clade were averaged as follows and taken as the absolute value: clade score = abs(mean(cor(query strain SNP, strain X SNP))). The clade with the highest similarity score (and thus the greatest proportion of shared SNPs) was selected as the most similar clade for that window. This process was repeated across the full genome. As a follow-up, this process was repeated by removing all strains within the query strain’s clade and the scores were recalculated.

The second phylogenetic approach used an expression matrix listing the 21 strains against all possible SNP positions present for each nonoverlapping window. For each cell, we assigned a value of 0 if the respective SNP was not present in the respective strain, a value of 1 if one copy of the SNP was present, and a value of 2 if both copies were present. A distance matrix was constructed from these data using a binary method (R dist), finally resulting in a phylogenetic tree (R hclust). The appropriate number of K clusters for the phylogenetic output of each window was estimated by traversing all 21 possible values (1 to 21) incrementally using R cutree. The first k-value that was chosen was the value that allowed the target strain to cluster with at least two other strains that were members of the same clade as defined by the current phylogenetic relationship among the sequenced strains. These criteria effectively eliminated clade E because it contains only one strain. Additionally, this approach was assessed manually for congruence across candidate regions.

Testing for admixtures.The NGSAdmix tool was used to perform measurements of genetic admixtures in the sequenced clinical isolates (62). The SNP table for each strain mapped to SC5314 was the input, and the data were compared to clade profiles produced using the three representative strains from each clade highlighted in the text for other comparisons as follows: for clade I, strains 12 C, L26, and P78048; for clade II, strains P57072, P76055, and P76067; for clade III, P34048, P78042, and P57055; for clade SA, strains P87, GC75, and P75063. The data were framed for four K-clusters and visualized.

Sanger sequencing.The association between the SNPs and indels identified in the next-generation-sequencing (NGS) data was tested by amplification of specific regions that were either shared among a number of strains (P75063) or strain specific (P60002). Two regions were PCR amplified using primers 5′ AGTCGGTGATGTCTATAGTG 3′ and 5′ GCTGTCCTTGGATCATTGAT 3′ to amplify Chr7:48017.48664 in P75063 and 5′ TTCTGCTGTTGCTGCTGCTA 3′/5′ CTGTCAACTGTCAACCAAAG 3′ to amplify ChrR: 19979596.1998145 in P60002. Amplicons were purified and subjected to Sanger sequencing.

Mitochondrial SNPs were verified using 3 sets of primers to amplify different regions of ChrM across the 21 natural isolates. Two isolates were analyzed from each clade, including strains P60002, P75010, and P94015. Primers 5′ TTAGTAGTGTCGGTGTCTTC 3′ and 5′ AGAGAGGGTTTTGGTTAGGG 3′ were used to PCR amplify ChrM: 4899.6076, 5′ GAATCTCAGAGACTACACGT 3′/5′ GTGGTATACGACGAGGCATT 3′ were used for ChrM: 18265.20660, and 5′ TGGGAAGTAGAGGCTGAAGA 3′/5′ AGGGGCATTATAAGGAGGAG 3′ were used for ChrM: 28094.29548. PCR amplicons were purified and subjected to Sanger sequencing.

Statistical testing.Statistical analyses were performed using Student’s t test unless otherwise indicated. All statistical tests were performed in the R 3.2.5 programming environment (63).

Accession number(s).The genome sequences used in this study are available under BioProject identifier (ID) PRJNA193498 (https://www.ncbi.nlm.nih.gov/bioproject). SNP data are available from dbSNP (https://www.ncbi.nlm.nih.gov/projects/SNP/) under noninclusive accession no. 1456786277 to 1457237021.

TEXT S1

Supplemental results. Download Text S1, PDF file, 0.3 MB.
Copyright © 2018 Wang et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

ACKNOWLEDGMENTS

We thank Iuliana Ene and Christophe D’Enfert for comments on the manuscript, Christina Cuomo for technical advice, and members of the R.J.B. and M.Z.A. laboratories for helpful discussions.

This work was supported by National Institutes of Health grants AI081704 and AI122011 (to R.J.B.), a PATH award from the Burroughs Wellcome Fund (to R.J.B.), and by a Karen T. Romer Undergraduate Teaching and Research Award (to J.M.W.).

FOOTNOTES

    • Received 5 June 2018
    • Accepted 10 August 2018
    • Published 18 September 2018
  • Copyright © 2018 Wang et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Bennett RJ,
    2. Forche A,
    3. Berman J
    . 2014. Rapid mechanisms for generating genome diversity: whole ploidy shifts, aneuploidy, and loss of heterozygosity. Cold Spring Harb Perspect Med 4:a019604. doi:10.1101/cshperspect.a019604.
    OpenUrlAbstract/FREE Full Text
  2. 2.↵
    1. Dunham MJ,
    2. Badrane H,
    3. Ferea T,
    4. Adams J,
    5. Brown PO,
    6. Rosenzweig F,
    7. Botstein D
    . 2002. Characteristic genome rearrangements in experimental evolution of Saccharomyces cerevisiae. Proc Natl Acad Sci U S A 99:16144–16149. doi:10.1073/pnas.242624799.
    OpenUrlAbstract/FREE Full Text
  3. 3.↵
    1. Goddard MR
    . 2016. Molecular evolution: sex accelerates adaptation. Nature 531:176–177. doi:10.1038/nature17304.
    OpenUrlCrossRefPubMed
  4. 4.↵
    1. McDonald MJ,
    2. Rice DP,
    3. Desai MM
    . 2016. Sex speeds adaptation by altering the dynamics of molecular evolution. Nature 531:233–236. doi:10.1038/nature17143.
    OpenUrlCrossRefPubMed
  5. 5.↵
    1. Lehtonen J,
    2. Jennions MD,
    3. Kokko H
    . 2012. The many costs of sex. Trends Ecol Evol 27:172–178. doi:10.1016/j.tree.2011.09.016.
    OpenUrlCrossRefPubMedWeb of Science
  6. 6.↵
    1. Roze D,
    2. Otto SP
    . 2012. Differential selection between the sexes and selection for sex. Evolution 66:558–574. doi:10.1111/j.1558-5646.2011.01459.x.
    OpenUrlCrossRefPubMedWeb of Science
  7. 7.↵
    1. Grimberg B,
    2. Zeyl C
    . 2005. The effects of sex and mutation rate on adaptation in test tubes and to mouse hosts by Saccharomyces cerevisiae. Evolution 59:431–438. doi:10.1111/j.0014-3820.2005.tb01001.x.
    OpenUrlCrossRefPubMedWeb of Science
  8. 8.↵
    1. Nielsen K,
    2. Heitman J
    . 2007. Sex and virulence of human pathogenic fungi. Adv Genet 57:143–173. doi:10.1016/S0065-2660(06)57004-X.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Becks L,
    2. Agrawal AF
    . 2010. Higher rates of sex evolve in spatially heterogeneous environments. Nature 468:89–92. doi:10.1038/nature09449.
    OpenUrlCrossRefPubMedWeb of Science
  10. 10.↵
    1. Horn DL,
    2. Neofytos D,
    3. Anaissie EJ,
    4. Fishman JA,
    5. Steinbach WJ,
    6. Olyaei AJ,
    7. Marr KA,
    8. Pfaller MA,
    9. Chang CH,
    10. Webster KM
    . 2009. Epidemiology and outcomes of candidemia in 2019 patients: data from the prospective antifungal therapy alliance registry. Clin Infect Dis 48:1695–1703. doi:10.1086/599039.
    OpenUrlCrossRefPubMedWeb of Science
  11. 11.↵
    1. Sardi JC,
    2. Scorzoni L,
    3. Bernardi T,
    4. Fusco-Almeida AM,
    5. Mendes Giannini MJ
    . 2013. Candida species: current epidemiology, pathogenicity, biofilm formation, natural antifungal products and new therapeutic options. J Med Microbiol 62:10–24. doi:10.1099/jmm.0.045054-0.
    OpenUrlCrossRefPubMedWeb of Science
  12. 12.↵
    1. Bennett RJ,
    2. Johnson AD
    . 2003. Completion of a parasexual cycle in Candida albicans by induced chromosome loss in tetraploid strains. EMBO J 22:2505–2515. doi:10.1093/emboj/cdg235.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    1. Forche A,
    2. Alby K,
    3. Schaefer D,
    4. Johnson AD,
    5. Berman J,
    6. Bennett RJ
    . 2008. The parasexual cycle in Candida albicans provides an alternative pathway to meiosis for the formation of recombinant strains. PLoS Biol 6:e110. doi:10.1371/journal.pbio.0060110.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Magee BB,
    2. Magee PT
    . 2000. Induction of mating in Candida albicans by construction of MTLa and MTLalpha strains. Science 289:310–313. doi:10.1126/science.289.5477.310.
    OpenUrlAbstract/FREE Full Text
  15. 15.↵
    1. Hull CM,
    2. Raisner RM,
    3. Johnson AD
    . 2000. Evidence for mating of the "asexual" yeast Candida albicans in a mammalian host. Science 289:307–310. doi:10.1126/science.289.5477.307.
    OpenUrlAbstract/FREE Full Text
  16. 16.↵
    1. Miller MG,
    2. Johnson AD
    . 2002. White-opaque switching in Candida albicans is controlled by mating-type locus homeodomain proteins and allows efficient mating. Cell 110:293–302. doi:10.1016/S0092-8674(02)00837-1.
    OpenUrlCrossRefPubMedWeb of Science
  17. 17.↵
    1. Alby K,
    2. Schaefer D,
    3. Bennett RJ
    . 2009. Homothallic and heterothallic mating in the opportunistic pathogen Candida albicans. Nature 460:890–893. doi:10.1038/nature08252.
    OpenUrlCrossRefPubMedWeb of Science
  18. 18.↵
    1. Keeney S
    . 2008. Spo11 and the formation of DNA double-strand breaks in meiosis. Genome Dyn Stab 2:81–123. doi:10.1007/7050_2007_026.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Tavanti A,
    2. Davidson AD,
    3. Fordyce MJ,
    4. Gow NA,
    5. Maiden MC,
    6. Odds FC
    . 2005. Population structure and properties of Candida albicans, as determined by multilocus sequence typing. J Clin Microbiol 43:5601–5613. doi:10.1128/JCM.43.11.5601-5613.2005.
    OpenUrlAbstract/FREE Full Text
  20. 20.↵
    1. Nebavi F,
    2. Ayala FJ,
    3. Renaud F,
    4. Bertout S,
    5. Eholie S,
    6. Moussa K,
    7. Mallie M,
    8. de Meeus T
    . 2006. Clonal population structure and genetic diversity of Candida albicans in AIDS patients from Abidjan (Cote d'Ivoire). Proc Natl Acad Sci U S A 103:5632–3668. doi:10.1073/pnas.0601704103.
    OpenUrlFREE Full Text
  21. 21.↵
    1. Odds FC,
    2. Bougnoux ME,
    3. Shaw DJ,
    4. Bain JM,
    5. Davidson AD,
    6. Diogo D,
    7. Jacobsen MD,
    8. Lecomte M,
    9. Li SY,
    10. Tavanti A,
    11. Maiden MC,
    12. Gow NA,
    13. d'Enfert C
    . 2007. Molecular phylogenetics of Candida albicans. Eukaryot Cell 6:1041–1052. doi:10.1128/EC.00041-07.
    OpenUrlAbstract/FREE Full Text
  22. 22.↵
    1. Anderson JB,
    2. Wickens C,
    3. Khan M,
    4. Cowen LE,
    5. Federspiel N,
    6. Jones T,
    7. Kohn LM
    . 2001. Infrequent genetic exchange and recombination in the mitochondrial genome of Candida albicans. J Bacteriol 183:865–872. doi:10.1128/JB.183.3.865-872.2001.
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    1. Jacobsen MD,
    2. Rattray AM,
    3. Gow NA,
    4. Odds FC,
    5. Shaw DJ
    . 2008. Mitochondrial haplotypes and recombination in Candida albicans. Med Mycol 46:647–654. doi:10.1080/13693780801986631.
    OpenUrlCrossRefPubMed
  24. 24.↵
    1. Zhang L,
    2. Yan L,
    3. Jiang J,
    4. Wang Y,
    5. Jiang Y,
    6. Yan T,
    7. Cao Y
    . 2014. The structure and retrotransposition mechanism of LTR-retrotransposons in the asexual yeast Candida albicans. Virulence 5:655–664. doi:10.4161/viru.32180.
    OpenUrlCrossRef
  25. 25.↵
    1. Diogo D,
    2. Bouchier C,
    3. d'Enfert C,
    4. Bougnoux ME
    . 2009. Loss of heterozygosity in commensal isolates of the asexual diploid yeast Candida albicans. Fungal Genet Biol 46:159–168. doi:10.1016/j.fgb.2008.11.005.
    OpenUrlCrossRefPubMedWeb of Science
  26. 26.↵
    1. Bougnoux ME,
    2. Pujol C,
    3. Diogo D,
    4. Bouchier C,
    5. Soll DR,
    6. d'Enfert C
    . 2008. Mating is rare within as well as between clades of the human pathogen Candida albicans. Fungal Genet Biol 45:221–231. doi:10.1016/j.fgb.2007.10.008.
    OpenUrlCrossRefPubMed
  27. 27.↵
    1. Hirakawa MP,
    2. Martinez DA,
    3. Sakthikumar S,
    4. Anderson MZ,
    5. Berlin A,
    6. Gujja S,
    7. Zeng Q,
    8. Zisson E,
    9. Wang JM,
    10. Greenberg JM,
    11. Berman J,
    12. Bennett RJ,
    13. Cuomo CA
    . 2015. Genetic and phenotypic intra-species variation in Candida albicans. Genome Res 25:413–425. doi:10.1101/gr.174623.114.
    OpenUrlAbstract/FREE Full Text
  28. 28.↵
    1. Wu W,
    2. Lockhart SR,
    3. Pujol C,
    4. Srikantha T,
    5. Soll DR
    . 2007. Heterozygosity of genes on the sex chromosome regulates Candida albicans virulence. Mol Microbiol 64:1587–1604. doi:10.1111/j.1365-2958.2007.05759.x.
    OpenUrlCrossRefPubMedWeb of Science
  29. 29.↵
    1. van Het Hoog M,
    2. Rast TJ,
    3. Martchenko M,
    4. Grindle S,
    5. Dignard D,
    6. Hogues H,
    7. Cuomo C,
    8. Berriman M,
    9. Scherer S,
    10. Magee BB,
    11. Whiteway M,
    12. Chibana H,
    13. Nantel A,
    14. Magee PT
    . 2007. Assembly of the Candida albicans genome into sixteen supercontigs aligned on the eight chromosomes. Genome Biol 8:R52. doi:10.1186/gb-2007-8-4-r52.
    OpenUrlCrossRefPubMed
  30. 30.↵
    1. Butler G,
    2. Rasmussen MD,
    3. Lin MF,
    4. Santos MAS,
    5. Sakthikumar S,
    6. Munro CA,
    7. Rheinbay E,
    8. Grabherr M,
    9. Forche A,
    10. Reedy JL,
    11. Agrafioti I,
    12. Arnaud MB,
    13. Bates S,
    14. Brown AJP,
    15. Brunke S,
    16. Costanzo MC,
    17. Fitzpatrick DA,
    18. de Groot PWJ,
    19. Harris D,
    20. Hoyer LL,
    21. Hube B,
    22. Klis FM,
    23. Kodira C,
    24. Lennard N,
    25. Logue ME,
    26. Martin R,
    27. Neiman AM,
    28. Nikolaou E,
    29. Quail MA,
    30. Quinn J,
    31. Santos MC,
    32. Schmitzberger FF,
    33. Sherlock G,
    34. Shah P,
    35. Silverstein KAT,
    36. Skrzypek MS,
    37. Soll D,
    38. Staggs R,
    39. Stansfield I,
    40. Stumpf MPH,
    41. Sudbery PE,
    42. Srikantha T,
    43. Zeng Q,
    44. Berman J,
    45. Berriman M,
    46. Heitman J,
    47. Gow NAR,
    48. Lorenz MC,
    49. Birren BW,
    50. Kellis M,
    51. Cuomo CA
    . 2009. Evolution of pathogenicity and sexual reproduction in eight Candida genomes. Nature 459:657–662. doi:10.1038/nature08064.
    OpenUrlCrossRefPubMedWeb of Science
  31. 31.↵
    1. Muzzey D,
    2. Schwartz K,
    3. Weissman JS,
    4. Sherlock G
    . 2013. Assembly of a phased diploid Candida albicans genome facilitates allele-specific measurements and provides a simple model for repeat and indel structure. Genome Biol 14:R97. doi:10.1186/gb-2013-14-9-r97.
    OpenUrlCrossRefPubMed
  32. 32.↵
    1. Tian D,
    2. Wang Q,
    3. Zhang P,
    4. Araki H,
    5. Yang S,
    6. Kreitman M,
    7. Nagylaki T,
    8. Hudson R,
    9. Bergelson J,
    10. Chen JQ
    . 2008. Single-nucleotide mutation rate increases close to insertions/deletions in eukaryotes. Nature 455:105–108. doi:10.1038/nature07175.
    OpenUrlCrossRefPubMedWeb of Science
  33. 33.↵
    1. Nikpay M,
    2. Goel A,
    3. Won H-H,
    4. Hall LM,
    5. Willenborg C,
    6. Kanoni S,
    7. Saleheen D,
    8. Kyriakou T,
    9. Nelson CP,
    10. Hopewell JC,
    11. Webb TR,
    12. Zeng L,
    13. Dehghan A,
    14. Alver M,
    15. Armasu SM,
    16. Auro K,
    17. Bjonnes A,
    18. Chasman DI,
    19. Chen S,
    20. Ford I,
    21. Franceschini N,
    22. Gieger C,
    23. Grace C,
    24. Gustafsson S,
    25. Huang J,
    26. Hwang S-J,
    27. Kim YK,
    28. Kleber ME,
    29. Lau KW,
    30. Lu X,
    31. Lu Y,
    32. Lyytikäinen L-P,
    33. Mihailov E,
    34. Morrison AC,
    35. Pervjakova N,
    36. Qu L,
    37. Rose LM,
    38. Salfati E,
    39. Saxena R,
    40. Scholz M,
    41. Smith AV,
    42. Tikkanen E,
    43. Uitterlinden A,
    44. Yang X,
    45. Zhang W,
    46. Zhao W,
    47. de Andrade M,
    48. de Vries PS,
    49. van Zuydam NR, et al
    . 2015. A comprehensive 1,000 genomes-based genome-wide association meta-analysis of coronary artery disease. Nat Genet 47:1121–1130. doi:10.1038/ng.3396.
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. Carbone I,
    2. Kohn L
    . 2004. Inferring process from pattern in fungal population genetics. Appl Mycol Biotechnol 4:30. doi:10.1016/S1874-5334(04)80005-4.
    OpenUrlCrossRef
  35. 35.↵
    1. Harvey PH,
    2. Pagel MD
    . 1991. The comparative method in evolutionary biology, Vol 239. Oxford University Press, Oxford, United Kingdom.
  36. 36.↵
    1. Chen XJ,
    2. Butow RA
    . 2005. The organization and inheritance of the mitochondrial genome. Nat Rev Genet 6:815–825. doi:10.1038/nrg1708.
    OpenUrlCrossRefPubMedWeb of Science
  37. 37.↵
    1. Hodgkinson A,
    2. Eyre-Walker A
    . 2011. Variation in the mutation rate across mammalian genomes. Nat Rev Genet 12:756–766. doi:10.1038/nrg3098.
    OpenUrlCrossRefPubMed
  38. 38.↵
    1. Yue JX,
    2. Li J,
    3. Aigrain L,
    4. Hallin J,
    5. Persson K,
    6. Oliver K,
    7. Bergstrom A,
    8. Coupland P,
    9. Warringer J,
    10. Lagomarsino MC,
    11. Fischer G,
    12. Durbin R,
    13. Liti G
    . 2017. Contrasting evolutionary genome dynamics between domesticated and wild yeasts. Nat Genet 49:913–924. doi:10.1038/ng.3847.
    OpenUrlCrossRef
  39. 39.↵
    1. Xie Z,
    2. Wang L,
    3. Wang L,
    4. Wang Z,
    5. Lu Z,
    6. Tian D,
    7. Yang S,
    8. Hurst LD
    . 2016. Mutation rate analysis via parent-progeny sequencing of the perennial peach. I. A low rate in woody perennials and a higher mutagenicity in hybrids. Proc Biol Sci 283:20161016. doi:10.1098/rspb.2016.1016.
    OpenUrlAbstract/FREE Full Text
  40. 40.↵
    1. Yang S,
    2. Wang L,
    3. Huang J,
    4. Zhang X,
    5. Yuan Y,
    6. Chen JQ,
    7. Hurst LD,
    8. Tian D
    . 2015. Parent-progeny sequencing indicates higher mutation rates in heterozygotes. Nature 523:463–467. doi:10.1038/nature14649.
    OpenUrlCrossRefPubMed
  41. 41.↵
    1. Jovelin R,
    2. Cutter AD
    . 2013. Fine-scale signatures of molecular evolution reconcile models of indel-associated mutation. Genome Biol Evol 5:978–986. doi:10.1093/gbe/evt051.
    OpenUrlCrossRefPubMed
  42. 42.↵
    1. Abyzov A,
    2. Li S,
    3. Kim DR,
    4. Mohiyuddin M,
    5. Stutz AM,
    6. Parrish NF,
    7. Mu XJ,
    8. Clark W,
    9. Chen K,
    10. Hurles M,
    11. Korbel JO,
    12. Lam HY,
    13. Lee C,
    14. Gerstein MB
    . 2015. Analysis of deletion breakpoints from 1,092 humans reveals details of mutation mechanisms. Nat Commun 6:7256. doi:10.1038/ncomms8256.
    OpenUrlCrossRefPubMed
  43. 43.↵
    1. Longman-Jacobsen N,
    2. Williamson JF,
    3. Dawkins RL,
    4. Gaudieri S
    . 2003. In polymorphic genomic regions indels cluster with nucleotide polymorphism: quantum genomics. Gene 312:257–261. doi:10.1016/S0378-1119(03)00621-8.
    OpenUrlCrossRefPubMedWeb of Science
  44. 44.↵
    1. McDonald MJ,
    2. Wang WC,
    3. Huang HD,
    4. Leu JY
    . 2011. Clusters of nucleotide substitutions and insertion/deletion mutations are associated with repeat sequences. PLoS Biol 9:e1000622. doi:10.1371/journal.pbio.1000622.
    OpenUrlCrossRefPubMed
  45. 45.↵
    1. Selmecki A,
    2. Forche A,
    3. Berman J
    . 2010. Genomic plasticity of the human fungal pathogen Candida albicans. Eukaryot Cell 9:991–1008. doi:10.1128/EC.00060-10.
    OpenUrlAbstract/FREE Full Text
  46. 46.↵
    1. Forche A,
    2. Abbey D,
    3. Pisithkul T,
    4. Weinzierl MA,
    5. Ringstrom T,
    6. Bruck D,
    7. Petersen K,
    8. Berman J
    . 2011. Stress alters rates and types of loss of heterozygosity in Candida albicans. mBio 2:e00129-11. doi:10.1128/mBio.00129-11.
    OpenUrlAbstract/FREE Full Text
  47. 47.↵
    1. Feri A,
    2. Loll-Krippleber R,
    3. Commere PH,
    4. Maufrais C,
    5. Sertour N,
    6. Schwartz K,
    7. Sherlock G,
    8. Bougnoux ME,
    9. d'Enfert C,
    10. Legrand M
    . 2016. Analysis of repair mechanisms following an induced double-strand break uncovers recessive deleterious alleles in the Candida albicans diploid genome. mBio 7:e01109-16. doi:10.1128/mBio.01109-16.
    OpenUrlAbstract/FREE Full Text
  48. 48.↵
    1. Ford CB,
    2. Funt JM,
    3. Abbey D,
    4. Issi L,
    5. Guiducci C,
    6. Martinez DA,
    7. Delorey T,
    8. Li BY,
    9. White TC,
    10. Cuomo C,
    11. Rao RP,
    12. Berman J,
    13. Thompson DA,
    14. Regev A
    . 2015. The evolution of drug resistance in clinical isolates of Candida albicans. Elife 4:e00662. doi:10.7554/eLife.00662.
    OpenUrlCrossRefPubMed
  49. 49.↵
    1. Hickman MA,
    2. Zeng G,
    3. Forche A,
    4. Hirakawa MP,
    5. Abbey D,
    6. Harrison BD,
    7. Wang YM,
    8. Su CH,
    9. Bennett RJ,
    10. Wang Y,
    11. Berman J
    . 2013. The ‘obligate diploid’ Candida albicans forms mating-competent haploids. Nature 494:55–59. doi:10.1038/nature11865.
    OpenUrlCrossRefPubMedWeb of Science
  50. 50.↵
    1. Selmecki A,
    2. Forche A,
    3. Berman J
    . 2006. Aneuploidy and isochromosome formation in drug-resistant Candida albicans. Science 313:367–370. doi:10.1126/science.1128242.
    OpenUrlAbstract/FREE Full Text
  51. 51.↵
    1. Muller HJ
    . 1964. The relation of recombination to mutational advance. Mutat Res 106:2–9.
    OpenUrlCrossRefPubMed
  52. 52.↵
    1. Tucker AE,
    2. Ackerman MS,
    3. Eads BD,
    4. Xu S,
    5. Lynch M
    . 2013. Population-genomic insights into the evolutionary origin and fate of obligately asexual Daphnia pulex. Proc Natl Acad Sci U S A 110:15740–15745. doi:10.1073/pnas.1313388110.
    OpenUrlAbstract/FREE Full Text
  53. 53.↵
    1. Pujol C,
    2. Reynes J,
    3. Renaud F,
    4. Raymond M,
    5. Tibayrenc M,
    6. Ayala FJ,
    7. Janbon F,
    8. Mallie M,
    9. Bastide JM
    . 1993. The yeast Candida albicans has a clonal mode of reproduction in a population of infected human immunodeficiency virus-positive patients. Proc Natl Acad Sci U S A 90:9456–9459. doi:10.1073/pnas.90.20.9456.
    OpenUrlAbstract/FREE Full Text
  54. 54.↵
    1. Schmid J,
    2. Cannon RD,
    3. Holland B
    . 2004. A futile act? Thoughts on the reproductive biology of Candida albicans. Mycologist 18:158–163. doi:10.1017/S0269915X0400415X.
    OpenUrlCrossRef
  55. 55.↵
    1. Odds FC
    . 2010. Molecular phylogenetics and epidemiology of Candida albicans. Future Microbiol 5:67–79. doi:10.2217/fmb.09.113.
    OpenUrlCrossRefPubMedWeb of Science
  56. 56.↵
    1. Zhang N,
    2. Magee BB,
    3. Magee PT,
    4. Holland BR,
    5. Rodrigues E,
    6. Holmes AR,
    7. Cannon RD,
    8. Schmid J
    . 2015. Selective advantages of a parasexual cycle for the yeast Candida albicans. Genetics 200:1117–1132. doi:10.1534/genetics.115.177170.
    OpenUrlAbstract/FREE Full Text
  57. 57.↵
    1. Liti G,
    2. Carter DM,
    3. Moses AM,
    4. Warringer J,
    5. Parts L,
    6. James SA,
    7. Davey RP,
    8. Roberts IN,
    9. Burt A,
    10. Koufopanou V,
    11. Tsai IJ,
    12. Bergman CM,
    13. Bensasson D,
    14. O’Kelly MJT,
    15. van Oudenaarden A,
    16. Barton DBH,
    17. Bailes E,
    18. Nguyen AN,
    19. Jones M,
    20. Quail MA,
    21. Goodhead I,
    22. Sims S,
    23. Smith F,
    24. Blomberg A,
    25. Durbin R,
    26. Louis EJ
    . 2009. Population genomics of domestic and wild yeasts. Nature 458:337–341. doi:10.1038/nature07743.
    OpenUrlCrossRefPubMedWeb of Science
  58. 58.↵
    1. Li H,
    2. Durbin R
    . 2010. Fast and accurate long-read alignment with burrows-Wheeler transform. Bioinformatics 26:589–595. doi:10.1093/bioinformatics/btp698.
    OpenUrlCrossRefPubMedWeb of Science
  59. 59.↵
    1. Li H,
    2. Handsaker B,
    3. Wysoker A,
    4. Fennell T,
    5. Ruan J,
    6. Homer N,
    7. Marth G,
    8. Abecasis G,
    9. Durbin R, 1000 Genome Project Data Processing Subgroup
    . 2009. The sequence alignment/map format and SAMtools. Bioinformatics 25:2078–2079. doi:10.1093/bioinformatics/btp352.
    OpenUrlCrossRefPubMedWeb of Science
  60. 60.↵
    1. McKenna A,
    2. Hanna M,
    3. Banks E,
    4. Sivachenko A,
    5. Cibulskis K,
    6. Kernytsky A,
    7. Garimella K,
    8. Altshuler D,
    9. Gabriel S,
    10. Daly M,
    11. DePristo MA
    . 2010. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res 20:1297–1303. doi:10.1101/gr.107524.110.
    OpenUrlAbstract/FREE Full Text
  61. 61.↵
    1. Robinson JT,
    2. Thorvaldsdottir H,
    3. Winckler W,
    4. Guttman M,
    5. Lander ES,
    6. Getz G,
    7. Mesirov JP
    . 2011. Integrative genomics viewer. Nat Biotechnol 29:24–26. doi:10.1038/nbt.1754.
    OpenUrlCrossRefPubMedWeb of Science
  62. 62.↵
    1. Skotte L,
    2. Korneliussen TS,
    3. Albrechtsen A
    . 2013. Estimating individual admixture proportions from next generation sequencing data. Genetics 195:693–702. doi:10.1534/genetics.113.154138.
    OpenUrlAbstract/FREE Full Text
  63. 63.↵
    Team RC. 2016. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.
View Abstract
PreviousNext
Back to top
Download PDF
Citation Tools
The Genome of the Human Pathogen Candida albicans Is Shaped by Mutation and Cryptic Sexual Recombination
Joshua M. Wang, Richard J. Bennett, Matthew Z. Anderson
mBio Sep 2018, 9 (5) e01205-18; DOI: 10.1128/mBio.01205-18

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mBio article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
The Genome of the Human Pathogen Candida albicans Is Shaped by Mutation and Cryptic Sexual Recombination
(Your Name) has forwarded a page to you from mBio
(Your Name) thought you would be interested in this article in mBio.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
The Genome of the Human Pathogen Candida albicans Is Shaped by Mutation and Cryptic Sexual Recombination
Joshua M. Wang, Richard J. Bennett, Matthew Z. Anderson
mBio Sep 2018, 9 (5) e01205-18; DOI: 10.1128/mBio.01205-18
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

Candida
evolution
loss of heterozygosity
parasex
recombination

Related Articles

Cited By...

About

  • About mBio
  • Editor in Chief
  • Board of Editors
  • AAM Fellows
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Article Types
  • Ethics
  • Contact Us

Follow #mBio

@ASMicrobiology

       

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2150-7511