Comparing genomic variant identification protocols for Candida auris.

Genomic analyses are widely applied to epidemiological, population genetic and experimental studies of pathogenic fungi. A wide range of methods are employed to carry out these analyses, typically without including controls that gauge the accuracy of variant prediction. The importance of tracking outbreaks at a global scale has raised the urgency of establishing high-accuracy pipelines that generate consistent results between research groups. To evaluate currently employed methods for whole-genome variant detection and elaborate best practices for fungal pathogens, we compared how 14 independent variant calling pipelines performed across 35 Candida auris isolates from 4 distinct clades and evaluated the performance of variant calling, single-nucleotide polymorphism (SNP) counts and phylogenetic inference results. Although these pipelines used different variant callers and filtering criteria, we found high overall agreement of SNPs from each pipeline. This concordance correlated with site quality, as SNPs discovered by a few pipelines tended to show lower mapping quality scores and depth of coverage than those recovered by all pipelines. We observed that the major differences between pipelines were due to variation in read trimming strategies, SNP calling methods and parameters, and downstream filtration criteria. We calculated specificity and sensitivity for each pipeline by aligning three isolates with chromosomal level assemblies and found that the GATK-based pipelines were well balanced between these metrics. Selection of trimming methods had a greater impact on SAMtools-based pipelines than those using GATK. Phylogenetic trees inferred by each pipeline showed high consistency at the clade level, but there was more variability between isolates from a single outbreak, with pipelines that used more stringent cutoffs having lower resolution. This project generated two truth datasets useful for routine benchmarking of C. auris variant calling, a consensus VCF of genotypes discovered by 10 or more pipelines across these 35 diverse isolates and variants for 2 samples identified from whole-genome alignments. This study provides a foundation for evaluating SNP calling pipelines and developing best practices for future fungal genomic studies.

Prognostic Value of IGF2 mRNA-Binding Protein 3 (IGF2BP3) Intratumoral Expression in Melanoma Patients at the Time of Diagnosis: Comparative Analysis of RT-qPCR Versus Immunohistochemistry.

Screening for prognostic biomarkers is crucial for clinical melanoma management. Insulin-like growth factor-II mRNA-binding protein 3 (IGF2BP3) has emerged as a potential melanoma diagnostic and prognostic biomarker. It is commonly tested by immunohistochemistry (IHC). Our study retrospectively examines IGF2BP3 mRNA and protein expression in primary melanomas, their correlation with clinicopathologic factors, clinical outcome, and selected miRNAs expression, and their efficiency in predicting melanoma progression and survival. RT-qPCR and IHC on IGF2BP3 expression were performed in 61 cryopreserved and 63 formalin-fixed paraffin-embedded primary melanomas, respectively, and correlated to clinicopathologic factors, distant metastasis-free survival (DMFS), and melanoma -specific survival (MSS). The correlation between RT-qPCR and IHC was significant but moderate. IGF2BP3 mRNA showed a stronger association with clinicopathologic factors (Breslow thickness, ulceration, mitosis rate, growth phase, development of metastasis, and melanoma-specific survival) than its protein counterpart. Interestingly, higher IGF2BP3 mRNA expression was detected in primary melanomas that further metastasized to distant sites and was an independent prognostic factor for the risk of unfavorable DMFS and MSS. RT-qPCR outperformed IHC in sensitivity and in predicting worse clinical outcomes. Therefore, RT-qPCR may successfully be implemented for routine IGF2BP3 assessing for the selection of melanoma patients with a higher risk of developing distant metastasis and dying of melanoma.

The Prognostic Value of miR-125b, miR-200c and miR-205 in Primary Cutaneous Malignant Melanoma Is Independent of BRAF Mutational Status.

BRAF mutations are present in around 50% of cutaneous malignant melanomas and are related to a poor outcome in advanced-stage melanoma patients. miRNAs are epigenetic regulators that modulate different cellular processes in cancer, including melanoma development and progression. However, there are no studies on the potential associations of the genetic alterations of the BRAF gene with miRNA expression in primary cutaneous melanomas. Here, in order to analyze the influence of BRAF mutations in the ability of selected miRNAs to predict clinical outcome and patient survival at the time of diagnosis, we studied the prognostic value of miR-125b, miR-200c and miR-205 expression depending on the BRAF mutational status in fresh, frozen primary tumor specimens. For this purpose, RNA was extracted for studying both BRAF mutations by Sanger sequencing and miRNA expression. Our results indicate that, although there seems to be a slight preference for their predictive ability in the BRAF mutated group, the expression of these three miRNAs serves effectively to predict the clinical outcome of melanoma patients independently of BRAF mutational status at the time of primary tumor diagnosis.

Chromosome-Level Genome Assembly of a Human Fungal Pathogen Reveals Synteny among Geographically Distinct Species.

Histoplasma capsulatum, a dimorphic fungal pathogen, is the most common cause of fungal respiratory infections in immunocompetent hosts. Histoplasma is endemic in the Ohio and Mississippi River Valleys in the United States and is also distributed worldwide. Previous studies have revealed at least eight clades, each specific to a geographic location: North American classes 1 and 2 (NAm 1 and NAm 2), Latin American groups A and B (LAm A and LAm B), Eurasian, Netherlands, Australian and African, and an additional distinct lineage (H81) comprised of Panamanian isolates. Previously assembled Histoplasma genomes are highly fragmented, with the highly repetitive G217B (NAm 2) strain, which has been used for most whole-genome-scale transcriptome studies, assembled into over 250 contigs. In this study, we set out to fully assemble the repeat regions and characterize the large-scale genome architecture of Histoplasma species. We resequenced five Histoplasma strains (WU24 [NAm 1], G217B [NAm 2], H88 [African], G186AR [Panama], and G184AR [Panama]) using Oxford Nanopore Technologies long-read sequencing technology. Here, we report chromosomal-level assemblies for all five strains, which exhibit extensive synteny among the geographically distant Histoplasma isolates. The new assemblies revealed that RYP2, a major regulator of morphology and virulence, is duplicated in G186AR. In addition, we mapped previously generated transcriptome data sets onto the newly assembled chromosomes. Our analyses revealed that the expression of transposons and transposon-embedded genes are upregulated in yeast phase compared to mycelial phase in the G217B and H88 strains. This study provides an important resource for fungal researchers and further highlights the importance of chromosomal-level assemblies in analyzing high-throughput data sets. IMPORTANCE Histoplasma species are dimorphic fungi causing significant morbidity and mortality worldwide. These fungi grow as mold in the soil and as budding yeast within the human host. Histoplasma can be isolated from soil in diverse regions, including North America, South America, Africa, and Europe. Phylogenetically distinct species of Histoplasma have been isolated and sequenced. However, for the commonly used strains, genome assemblies have been fragmented, leading to underutilization of genome-scale data. This study provides chromosome-level assemblies of the commonly used Histoplasma strains using long-read sequencing technology. Comparative analysis of these genomes shows largely conserved gene order within the chromosomes. Mapping existing transcriptome data on these new assemblies reveals clustering of transcriptionally coregulated genes. The results of this study highlight the importance of obtaining chromosome-level assemblies in understanding the biology of human fungal pathogens.

In vivo emergence of high-level resistance during treatment reveals the first identified mechanism of amphotericin B resistance in Candida auris.

OBJECTIVE: Candida auris has emerged as a health-care-associated and multidrug-resistant fungal pathogen of great clinical concern. As many as 50% of C. auris clinical isolates are reported to be resistant to amphotericin B, but no mechanisms contributing to this resistance have been identified. Here we describe a clinical case in which high-level amphotericin B resistance was acquired in vivo during therapy and undertake molecular and genetic studies to identify and characterize the genetic determinant of resistance. METHODS: Whole-genome sequencing was performed on four C. auris isolates obtained from a single patient case. Cas9-mediated genetic manipulations were then used to generate mutant strains harbouring mutations of interest, and these strains were subsequently subjected to amphotericin B susceptibility testing and comprehensive sterol profiling. RESULTS: A novel mutation in the C. auris sterol-methyltransferase gene ERG6 was found to be associated with amphotericin B resistance, and this mutation alone conferred a >32-fold increase in amphotericin B resistance. Comprehensive sterol profiling revealed an abrogation of ergosterol biosynthesis and a corresponding accumulation of cholesta-type sterols in isolates and strains harbouring the clinically derived ERG6 mutation. CONCLUSIONS: Together these findings definitively demonstrate mutations in C. auris ERG6 as the first identified mechanism of clinical amphotericin B resistance in C. auris and represent a significant step forward in the understanding of antifungal resistance in this emerging public health threat.

Diversity, multifaceted evolution, and facultative saprotrophism in the European Batrachochytrium salamandrivorans epidemic.

While emerging fungi threaten global biodiversity, the paucity of fungal genome assemblies impedes thoroughly characterizing epidemics and developing effective mitigation strategies. Here, we generate de novo genomic assemblies for six outbreaks of the emerging pathogen Batrachochytrium salamandrivorans (Bsal). We reveal the European epidemic currently damaging amphibian populations to comprise multiple, highly divergent lineages demonstrating isolate-specific adaptations and metabolic capacities. In particular, we show extensive gene family expansions and acquisitions, through a variety of evolutionary mechanisms, and an isolate-specific saprotrophic lifecycle. This finding both explains the chytrid's ability to divorce transmission from host density, producing Bsal's enigmatic host population declines, and is a key consideration in developing successful mitigation measures.

In Vitro Antifungal Resistance of Candida auris Isolates from Bloodstream Infections, South Africa.

Candida auris is a multidrug-resistant fungal pathogen that is endemic in South African hospitals. We tested bloodstream C. auris isolates that were submitted to a reference laboratory for national laboratory-based surveillance for candidemia in 2016 and 2017. We confirmed the species identification by phenotypic/molecular methods. We tested susceptibility to amphotericin B, anidulafungin, caspofungin, micafungin, itraconazole, posaconazole, voriconazole, fluconazole, and flucytosine using broth microdilution and Etest methods. We interpreted MICs using tentative breakpoints. We sequenced the genomes of a subset of isolates and compared them to the C. auris B8441 reference strain. Of 400 C. auris isolates, 361 (90%) were resistant to at least one antifungal agent, 339 (94%) to fluconazole alone (MICs of ≥32 µg/ml), 19 (6%) to fluconazole and amphotericin B (MICs of ≥2 µg/ml), and 1 (0.3%) to amphotericin B alone. Two (0.5%) isolates from a single patient were pan-resistant (resistant to fluconazole, amphotericin B, and echinocandins). Of 92 isolates selected for whole-genome sequencing, 77 clustered in clade III, including the pan-resistant isolates, 13 in clade I, and 2 in clade IV. Eighty-four of the isolates (91%) were resistant to at least one antifungal agent; both resistant and susceptible isolates had mutations. The common substitutions identified across the different clades were VF125AL, Y132F, K177R, N335S, and E343D in ERG11; N647T in MRR1; A651P, A657V, and S195G in TAC1b; S639P in FKS1HP1; and S58T in ERG3. Most South African C. auris isolates were resistant to azoles, although resistance to polyenes and echinocandins was less common. We observed mutations in resistance genes even in phenotypically susceptible isolates.

Cutaneous Lymphadenoma Is a Distinct Trichoblastoma-like Lymphoepithelial Tumor With Diffuse Androgen Receptor Immunoreactivity, Notch1 Ligand in Reed-Sternberg-like Cells, and Common EGFR Somatic Mutations.

The term "cutaneous lymphadenoma" was coined in this journal for an unusual lymphoepithelial cutaneous adnexal neoplasm, possibly with immature pilosebaceous differentiation. Some authors further proposed that cutaneous lymphadenoma was an adamantinoid trichoblastoma. However, although a hair follicle differentiation is widely accepted, the fact that this is a lymphoepithelial tumor is not appropriately explained by the trichoblastoma hypothesis. Our goal was to further clarify the phenotypic and genotypic features of cutaneous lymphadenoma in a series of 11 cases. Histologically, a lobular architecture surrounded by a dense fibrous stroma was present in all cases. The lobules were composed of epithelial cells admixtured with small lymphocytes and isolated or clustered large Reed-Sternberg-like (RS-L) cells. The epithelial cells were diffusely positive for the hair follicle stem cell markers CK15, PHLDA1, and for androgen receptor. No immunostaining for markers of sebaceous differentiation was found. Intraepithelial lymphocytes were predominantly CD3+, CD4+, FoxP3+ T cells. RS-L cells showed both strong Jagged-1 and Notch1 cytoplasmic immunostaining. Androgen-regulated NKX3.1 nuclear immunostaining was present in a subset of large intralobular cells in all cases. Double immunostaining showed coexpression of NKX3.1 and CD30 in a subset of RS-L cells. No immunostaining for lymphocytic or epithelial markers was present in RS-L cells. EGFR, PIK3CA, and FGFR3 somatic mutations were found by next-generation sequencing in 56% of the cases. We consider that cutaneous lymphadenoma is a distinct benign lymphoepithelial tumor with androgen receptor and hair follicle bulge stem cell marker expression, RS-L cell-derived Notch1 ligand, and common EGFR gene mutations.

Best practices on the differential expression analysis of multi-species RNA-seq.

Advances in transcriptome sequencing allow for simultaneous interrogation of differentially expressed genes from multiple species originating from a single RNA sample, termed dual or multi-species transcriptomics. Compared to single-species differential expression analysis, the design of multi-species differential expression experiments must account for the relative abundances of each organism of interest within the sample, often requiring enrichment methods and yielding differences in total read counts across samples. The analysis of multi-species transcriptomics datasets requires modifications to the alignment, quantification, and downstream analysis steps compared to the single-species analysis pipelines. We describe best practices for multi-species transcriptomics and differential gene expression.

Genome-Wide Analysis of Experimentally Evolved Candida auris Reveals Multiple Novel Mechanisms of Multidrug Resistance.

Candida auris is globally recognized as an opportunistic fungal pathogen of high concern, due to its extensive multidrug resistance (MDR). Still, molecular mechanisms of MDR are largely unexplored. This is the first account of genome-wide evolution of MDR in C. auris obtained through serial in vitro exposure to azoles, polyenes, and echinocandins. We show the stepwise accumulation of copy number variations and novel mutations in genes both known and unknown in antifungal drug resistance. Echinocandin resistance was accompanied by a codon deletion in FKS1 hot spot 1 and a substitution in FKS1 "novel" hot spot 3. Mutations in ERG3 and CIS2 further increased the echinocandin MIC. Decreased azole susceptibility was linked to a mutation in transcription factor TAC1b and overexpression of the drug efflux pump Cdr1, a segmental duplication of chromosome 1 containing ERG11, and a whole chromosome 5 duplication, which contains TAC1b The latter was associated with increased expression of ERG11, TAC1b, and CDR2 but not CDR1 The simultaneous emergence of nonsense mutations in ERG3 and ERG11 was shown to decrease amphotericin B susceptibility, accompanied with fluconazole cross-resistance. A mutation in MEC3, a gene mainly known for its role in DNA damage homeostasis, further increased the polyene MIC. Overall, this study shows the alarming potential for and diversity of MDR development in C. auris, even in a clade until now not associated with MDR (clade II), stressing its clinical importance and the urge for future research.IMPORTANCECandida auris is a recently discovered human fungal pathogen and has shown an alarming potential for developing multi- and pan-resistance toward all classes of antifungals most commonly used in the clinic. Currently, C. auris has been globally recognized as a nosocomial pathogen of high concern due to this evolutionary potential. So far, this is the first study in which the stepwise progression of multidrug resistance (MDR) in C. auris is monitored in vitro Multiple novel mutations in known resistance genes and genes previously not or vaguely associated with drug resistance reveal rapid MDR evolution in a C. auris clade II isolate. Additionally, this study shows that in vitro experimental evolution can be a powerful tool to discover new drug resistance mechanisms, although it has its limitations.

Clade-specific chromosomal rearrangements and loss of subtelomeric adhesins in Candida auris.

Candida auris is an emerging fungal pathogen of rising concern due to global spread, the ability to cause healthcare-associated outbreaks, and antifungal resistance. Genomic analyses revealed that early contemporaneously detected cases of C. auris were geographically stratified into four major clades. While Clades I, III, and IV are responsible for ongoing outbreaks of invasive and multidrug-resistant infections, Clade II, also termed the East Asian clade, consists primarily of cases of ear infection, is often susceptible to all antifungal drugs, and has not been associated with outbreaks. Here, we generate chromosome-level assemblies of twelve isolates representing the phylogenetic breadth of these four clades and the only isolate described to date from Clade V. This Clade V genome is highly syntenic with those of Clades I, III, and IV, although the sequence is highly divergent from the other clades. Clade II genomes appear highly rearranged, with translocations occurring near GC-poor regions, and large subtelomeric deletions in most chromosomes, resulting in a substantially different karyotype. Rearrangements and deletion lengths vary across Clade II isolates, including two from a single patient, supporting ongoing genome instability. Deleted subtelomeric regions are enriched in Hyr/Iff-like cell-surface proteins, novel candidate cell wall proteins, and an ALS-like adhesin. Cell wall proteins from these families and other drug-related genes show clade-specific signatures of selection in Clades I, III, and IV. Subtelomeric dynamics and the conservation of cell surface proteins in the clades responsible for global outbreaks causing invasive infections suggest an explanation for the different phenotypes observed between clades.

Blastomycosis in Africa and the Middle East: A Comprehensive Review of Reported Cases and Reanalysis of Historical Isolates Based on Molecular Data.

BACKGROUND: Blastomycosis has been reported from countries in Africa and the Middle East, but a decades-long debate has persisted regarding whether this is the same disease known in North America and caused by Blastomyces dermatitidis and Blastomyces gilchristii. METHODS: We reviewed published cases of human and veterinary blastomycosis from Africa and the Middle East. We abstracted epidemiological and clinical features of cases, including sites of disease, diagnosis, management, outcomes, and, where available, genetic and antigenic typing of case isolates. In addition, we sequenced nucleic acids from 9 clinical isolates from Africa deposited in global collections as B. dermatitidis; for 5, we sequenced the internal transcribed spacer regions, and for the other 4 we sequenced the whole genomes. RESULTS: We identified 172 unique human patients with blastomycosis, including 159 patients from 25 African countries and 12 patients from 5 Middle Eastern countries, and also identified 7 reports of veterinary blastomycosis. In humans, cutaneous disease predominated (n = 100/137, 73%), followed by pulmonary (n = 73/129, 57%) and osteoarticular involvement (n = 61/128, 48%). Unusual direct microscopy/histopathological presentations included short hyphal fragments in tissues (n = 23/129, 18%). There were 34 genotyped case isolates that comprised 4 species: Blastomyces percursus (n = 22, 65%), from 8 countries throughout all regions; Blastomyces emzantsi (n = 9, 26%), from South Africa; B. dermatitidis (n = 1, 3%), from the Democratic Republic of Congo; and B. gilchristii (n = 2, 6%), from South Africa and Zimbabwe. CONCLUSIONS: Blastomycosis occurs throughout Africa and the Middle East and is caused predominantly by B. percursus and, at least in South Africa, B. emzantsi, resulting in distinct clinical and pathological patterns of disease.

An oxindole efflux inhibitor potentiates azoles and impairs virulence in the fungal pathogen Candida auris.

Candida auris is an emerging fungal pathogen that exhibits resistance to multiple drugs, including the most commonly prescribed antifungal, fluconazole. Here, we use a combinatorial screening approach to identify a bis-benzodioxolylindolinone (azoffluxin) that synergizes with fluconazole against C. auris. Azoffluxin enhances fluconazole activity through the inhibition of efflux pump Cdr1, thus increasing intracellular fluconazole levels. This activity is conserved across most C. auris clades, with the exception of clade III. Azoffluxin also inhibits efflux in highly azole-resistant strains of Candida albicans, another human fungal pathogen, increasing their susceptibility to fluconazole. Furthermore, azoffluxin enhances fluconazole activity in mice infected with C. auris, reducing fungal burden. Our findings suggest that pharmacologically targeting Cdr1 in combination with azoles may be an effective strategy to control infection caused by azole-resistant isolates of C. auris.

Mitochondrial Genome Sequences of the Emerging Fungal Pathogen Candida auris.

Candida auris is an emerging fungal pathogen capable of causing invasive infections in humans. Since its first appearance around 1996, it has been isolated in countries spanning five continents. C. auris is a yeast that has the potential to cause outbreaks in hospitals, can survive in adverse conditions, including dry surfaces and high temperatures, and has been frequently misidentified by traditional methods. Furthermore, strains have been identified that are resistant to two and even all three of the main classes of antifungals currently in use. Several nuclear genome assemblies of C. auris have been published representing different clades and continents, yet until recently, the mitochondrial genomes (mtDNA chromosomes) of this species and the closely related species of C. haemulonii, C. duobushaemulonii, and C. pseudohaemulonii had not been analyzed in depth. We used reads from PacBio and Illumina sequencing to obtain a de novo reference assembly of the mitochondrial genome of the C. auris clade I isolate B8441 from Pakistan. This assembly has a total size of 28.2 kb and contains 13 core protein-coding genes, 25 tRNAs and the 12S and 16S ribosomal subunits. We then performed a comparative analysis by aligning Illumina reads of 129 other isolates from South Asia, Japan, South Africa, and South America with the B8441 reference. The clades of the phylogenetic tree we obtained from the aligned mtDNA sequences were consistent with those derived from the nuclear genome. The mitochondrial genome revealed a generally low genetic variation within clades, although the South Asian clade displayed two sub-branches including strains from both Pakistan and India. In particular, the 86 isolates from Colombia and Venezuela had mtDNA sequences that were all identical at the base level, i.e., a single conserved haplotype or mitochondrial background that exhibited characteristic differences from the Pakistan reference isolate B8441, such as a unique 25-nt insert that may affect function.

Experimental Evolution Identifies Adaptive Aneuploidy as a Mechanism of Fluconazole Resistance in Candida auris.

Candida auris is a newly emerging fungal pathogen of humans and has attracted considerable attention from both the clinical and basic research communities. Clinical isolates of C. auris are often resistant to one or more antifungal agents. To explore how antifungal resistance develops, we performed experimental evolution assays using a fluconazole-susceptible isolate of C. auris (BJCA001). After a series of passages through medium containing increasing concentrations of fluconazole, fungal cells acquired resistance. By sequencing and comparing the genomes of the parental fluconazole-susceptible strain and 26 experimentally evolved strains of C. auris, we found that a portion of fluconazole-resistant strains carried one extra copy of chromosome V. In the absence of fluconazole, C. auris cells rapidly became susceptible and lost the extra copy of chromosome V. Genomic and transcriptome sequencing (RNA-Seq) analyses indicate that this chromosome carries a number of drug resistance-related genes, which were transcriptionally upregulated in the resistant, aneuploid strains. Moreover, missense mutations were identified in the genes TAC1B, RRP6, and SFT2 in all experimentally evolved strains. Our findings suggest that the gain of an extra copy of chromosome V is associated with the rapid acquisition of fluconazole resistance and may represent an important evolutionary mechanism of antifungal resistance in C. auris.

Super resolution for root imaging.

PREMISE: High-resolution cameras are very helpful for plant phenotyping as their images enable tasks such as target vs. background discrimination and the measurement and analysis of fine above-ground plant attributes. However, the acquisition of high-resolution images of plant roots is more challenging than above-ground data collection. An effective super-resolution (SR) algorithm is therefore needed for overcoming the resolution limitations of sensors, reducing storage space requirements, and boosting the performance of subsequent analyses. METHODS: We propose an SR framework for enhancing images of plant roots using convolutional neural networks. We compare three alternatives for training the SR model: (i) training with non-plant-root images, (ii) training with plant-root images, and (iii) pretraining the model with non-plant-root images and fine-tuning with plant-root images. The architectures of the SR models were based on two state-of-the-art deep learning approaches: a fast SR convolutional neural network and an SR generative adversarial network. RESULTS: In our experiments, we observed that the SR models improved the quality of low-resolution images of plant roots in an unseen data set in terms of the signal-to-noise ratio. We used a collection of publicly available data sets to demonstrate that the SR models outperform the basic bicubic interpolation, even when trained with non-root data sets. DISCUSSION: The incorporation of a deep learning-based SR model in the imaging process enhances the quality of low-resolution images of plant roots. We demonstrate that SR preprocessing boosts the performance of a machine learning system trained to separate plant roots from their background. Our segmentation experiments also show that high performance on this task can be achieved independently of the signal-to-noise ratio. We therefore conclude that the quality of the image enhancement depends on the desired application.

Updates and Comparative Analysis of the Mitochondrial Genomes of Paracoccidioides spp. Using Oxford Nanopore MinION Sequencing.

The mitochondrial genome of the Paracoccidioides brasiliensis reference isolate Pb18 was first sequenced and described by Cardoso et al. (2007), as a circular genome with a size of 71.3 kb and containing 14 protein coding genes, 25 tRNAs, and the large and small subunits of ribosomal RNA. Later in 2011, Desjardins et al. (2011) obtained partial assemblies of mitochondrial genomes of P. lutzii (Pb01), P. americana (Pb03), and P. brasiliensis sensu stricto (Pb18), although with a size of only 43.1 kb for Pb18. Sequencing errors or other limitations resulting from earlier technologies, and the advantages of NGS (short and long reads), prompted us to improve and update the mtDNA sequences and annotations of two Paracoccidioides species. Using Oxford Nanopore and Illumina read sequencing, we generated high-quality complete de novo mitochondrial genome assemblies and annotations for P. brasiliensis (Pb18) and P. americana (Pb03). Both assemblies were characterized by an unusually long spacer or intron region (>50 kb) between exons 2 and 3 of the nad5 gene, which was moderately conserved between Pb03 and Pb18 but not similar to other reported sequences, except for an unassigned contig in the 2011 assembly of Pb03. The reliability of the insert missing from previous mtDNA genome assemblies was confirmed by inspection of the individual Nanopore read sequences containing nad5 coding DNA, and experimentally by PCR for Pb18. We propose that the insert may aid replication initiation and may be excised to produce a smaller structural variant. The updated mtDNA genomes should enable more accurate SNP and other comparative or evolutionary analyses and primer/probe designs. A comparative analysis of the mtDNA from 32 isolates of Paracoccidioides spp., using the SNPs of the aligned mitochondrial genomes, showed groupings within the brasiliensis species complex that were largely consistent with previous findings from only five mitochondrial loci.

Understanding the Emergence of Multidrug-Resistant Candida: Using Whole-Genome Sequencing to Describe the Population Structure of Candida haemulonii Species Complex.

The recent emergence of a multidrug-resistant yeast, Candida auris, has drawn attention to the closely related species from the Candida haemulonii complex that include C. haemulonii, Candida duobushaemulonii, Candida pseudohaemulonii, and the recently identified Candida vulturna. Here, we used antifungal susceptibility testing and whole-genome sequencing (WGS) to investigate drug resistance and genetic diversity among isolates of C. haemulonii complex from different geographic areas in order to assess population structure and the extent of clonality among strains. Although most isolates of all four species were genetically distinct, we detected evidence of the in-hospital transmission of C. haemulonii and C. duobushaemulonii in one hospital in Panama, indicating that these species are also capable of causing outbreaks in healthcare settings. We also detected evidence of the rising azole resistance among isolates of C. haemulonii and C. duobushaemulonii in Colombia, Panama, and Venezuela linked to substitutions in ERG11 gene as well as amplification of this gene in C. haemulonii in isolates in Colombia suggesting the presence of evolutionary pressure for developing azole resistance in this region. Our results demonstrate that these species need to be monitored as possible causes of outbreaks of invasive infection.

Mutations in TAC1B: a Novel Genetic Determinant of Clinical Fluconazole Resistance in Candida auris.

Candida auris has emerged as a multidrug-resistant pathogen of great clinical concern. Approximately 90% of clinical C. auris isolates are resistant to fluconazole, the most commonly prescribed antifungal agent, and yet it remains unknown what mechanisms underpin this fluconazole resistance. To identify novel mechanisms contributing to fluconazole resistance in C. auris, fluconazole-susceptible C. auris clinical isolate AR0387 was passaged in media supplemented with fluconazole to generate derivative strains which had acquired increased fluconazole resistance in vitro Comparative analyses of comprehensive sterol profiles, [3H]fluconazole uptake, sequencing of C. auris genes homologous to genes known to contribute to fluconazole resistance in other species of Candida, and relative expression levels of C. aurisERG11, CDR1, and MDR1 were performed. All fluconazole-evolved derivative strains were found to have acquired mutations in the zinc-cluster transcription factor-encoding gene TAC1B and to show a corresponding increase in CDR1 expression relative to the parental clinical isolate, AR0387. Mutations in TAC1B were also identified in a set of 304 globally distributed C. auris clinical isolates representing each of the four major clades. Introduction of the most common mutation found among fluconazole-resistant clinical isolates of C. auris into fluconazole-susceptible isolate AR0387 was confirmed to increase fluconazole resistance by 8-fold, and the correction of the same mutation in a fluconazole-resistant isolate, AR0390, decreased fluconazole MIC by 16-fold. Taken together, these data demonstrate that C. auris can rapidly acquire resistance to fluconazole in vitro and that mutations in TAC1B significantly contribute to clinical fluconazole resistance.IMPORTANCECandida auris is an emerging multidrug-resistant pathogen of global concern, known to be responsible for outbreaks on six continents and to be commonly resistant to antifungals. While the vast majority of clinical C. auris isolates are highly resistant to fluconazole, an essential part of the available antifungal arsenal, very little is known about the mechanisms contributing to resistance. In this work, we show that mutations in the transcription factor TAC1B significantly contribute to clinical fluconazole resistance. These studies demonstrated that mutations in TAC1B can arise rapidly in vitro upon exposure to fluconazole and that a multitude of resistance-associated TAC1B mutations are present among the majority of fluconazole-resistant C. auris isolates from a global collection and appear specific to a subset of lineages or clades. Thus, identification of this novel genetic determinant of resistance significantly adds to the understanding of clinical antifungal resistance in C. auris.

Tracing the Evolutionary History and Global Expansion of Candida auris Using Population Genomic Analyses.

Candida auris has emerged globally as a multidrug-resistant yeast that can spread via nosocomial transmission. An initial phylogenetic study of isolates from Japan, India, Pakistan, South Africa, and Venezuela revealed four populations (clades I, II, III, and IV) corresponding to these geographic regions. Since this description, C. auris has been reported in more than 30 additional countries. To trace this global emergence, we compared the genomes of 304 C. auris isolates from 19 countries on six continents. We found that four predominant clades persist across wide geographic locations. We observed phylogeographic mixing in most clades; clade IV, with isolates mainly from South America, demonstrated the strongest phylogeographic substructure. C. auris isolates from two clades with opposite mating types were detected contemporaneously in a single health care facility in Kenya. We estimated a Bayesian molecular clock phylogeny and dated the origin of each clade within the last 360 years; outbreak-causing clusters from clades I, III, and IV originated 36 to 38 years ago. We observed high rates of antifungal resistance in clade I, including four isolates resistant to all three major classes of antifungals. Mutations that contribute to resistance varied between the clades, with Y132F in ERG11 as the most widespread mutation associated with azole resistance and S639P in FKS1 for echinocandin resistance. Copy number variants in ERG11 predominantly appeared in clade III and were associated with fluconazole resistance. These results provide a global context for the phylogeography, population structure, and mechanisms associated with antifungal resistance in C. aurisIMPORTANCE In less than a decade, C. auris has emerged in health care settings worldwide; this species is capable of colonizing skin and causing outbreaks of invasive candidiasis. In contrast to other Candida species, C. auris is unique in its ability to spread via nosocomial transmission and its high rates of drug resistance. As part of the public health response, whole-genome sequencing has played a major role in characterizing transmission dynamics and detecting new C. auris introductions. Through a global collaboration, we assessed genome evolution of isolates of C. auris from 19 countries. Here, we described estimated timing of the expansion of each C. auris clade and of fluconazole resistance, characterized discrete phylogeographic population structure of each clade, and compared genome data to sensitivity measurements to describe how antifungal resistance mechanisms vary across the population. These efforts are critical for a sustained, robust public health response that effectively utilizes molecular epidemiology.