Systematic Analysis of Copy Number Variations in the Pathogenic Yeast Candida parapsilosis Identifies a Gene Amplification in RTA3 That is Associated with Drug Resistance

We analyzed the genomes of 170 C. parapsilosis isolates and identified multiple copy number variations (CNVs). We identified two genes, RTA3 (CPAR2_104610) and ARR3 (CPAR2_601050), each of which was the target of multiple independent amplification events. Phylogenetic analysis shows that most of these amplifications originated only once. For ARR3, which encodes a putative arsenate transporter, 8 distinct CNVs were identified, ranging in size from 2.3 kb to 10.5 kb with 3 to 23 copies. For RTA3, 16 distinct CNVs were identified, ranging in size from 0.3 kb to 4.5 kb with 2 to ~50 copies. One unusual amplification resulted in a DUP-TRP/INV-DUP structure similar to some human CNVs. RTA3 encodes a putative phosphatidylcholine (PC) floppase which is known to regulate the inward translocation of PC in Candida albicans. We found that an increased copy number of RTA3 correlated with resistance to miltefosine, an alkylphosphocholine drug that affects PC metabolism. Additionally, we conducted an adaptive laboratory evolution experiment in which two C. parapsilosis isolates were cultured in increasing concentrations of miltefosine. Two genes, CPAR2_303950 and CPAR2_102700, coding for putative PC flippases homologous to S. cerevisiae DNF1 gained homozygous protein-disrupting mutations in the evolved strains. Overall, our results show that C. parapsilosis can gain resistance to miltefosine, a drug that has recently been granted orphan drug designation approval by the United States Food and Drug Administration for the treatment of invasive candidiasis, through both CNVs or loss-of-function alleles in one of the flippase genes.