Bimodal distribution of azole susceptibility in Sporothrix brasiliensis isolates in Brazil

Sporothrix brasiliensis is an emerging zoonotic fungal pathogen that can be difficult to treat. Antifungal susceptibility testing was performed on the mold phase of a convenience sample of 61 Sporothrix spp. isolates from human and cat sporotrichosis cases in Brazil using the Clinical and Laboratory Standards Institute standard M38. A bimodal distribution of azole susceptibility was observed with 50% (28/56) of S. brasiliensis isolates showing elevated itraconazole minimum inhibitory concentrations ≥16 µg/mL. Phylogenetic analysis found the in vitro resistant isolates were not clonal and were distributed across three different S. brasiliensis clades. Single nucleotide polymorphism (SNP) analysis was performed to identify potential mechanisms of in vitro resistance. Two of the 28 resistant isolates (MIC ≥16 mg/L) had a polymorphism in the cytochrome P450 gene, cyp51, corresponding to the well-known G448S substitution inducing azole resistance in Aspergillus fumigatus. SNPs corresponding to other known mechanisms of azole resistance were not identified in the remaining 26 in vitro resistant isolates.

Emergence of zoonotic sporotrichosis in Brazil: a genomic epidemiology study

BACKGROUND

Zoonotic sporotrichosis is a neglected fungal disease, whereby outbreaks are primarily driven by Sporothrix brasiliensis and linked to cat-to-human transmission. To understand the emergence and spread of sporotrichosis in Brazil, the epicentre of the current epidemic in South America, we aimed to conduct whole-genome sequencing (WGS) to describe the genomic epidemiology.

METHODS

In this genomic epidemiology study, we included Sporothrix spp isolates from sporotrichosis cases from Brazil, Colombia, and the USA. We conducted WGS using Illumina NovaSeq on isolates collected by three laboratories in Brazil from humans and cats with sporotrichosis between 2013 and 2022. All isolates that were confirmed to be Sporothrix genus by internal transcribed spacer or beta-tubulin PCR sequencing were included in this study. We downloaded eight Sporothrix genome sequences from the National Center for Biotechnology Information (six from Brazil, two from Colombia). Three Sporothrix spp genome sequences from the USA were generated by the US Centers for Disease Control and Prevention as part of this study. We did phylogenetic analyses and correlated geographical and temporal case distribution with genotypic features of Sporothrix spp isolates.

FINDINGS

72 Sporothrix spp isolates from 55 human and 17 animal sporotrichosis cases were included: 67 (93%) were from Brazil, two (3%) from Colombia, and three (4%) from the USA. Cases spanned from 1999 to 2022. Most (61 [85%]) isolates were S brasiliensis, and all were reported from Brazil. Ten (14%) were Sporothrix schenckii and were reported from Brazil, USA, and Colombia. For S schenckii isolates, two distinct clades were observed wherein isolates clustered by geography. For S brasiliensis isolates, five clades separated by more than 100 000 single-nucleotide polymorphisms were observed. Among the five S brasiliensis clades, clades A and C contained isolates from both human and cat cases, and clade A contained isolates from six different states in Brazil. Compared with S brasiliensis isolates, larger genetic diversity was observed among S schenckii isolates from animal and human cases within a clade.

INTERPRETATION

Our results suggest that the ongoing epidemic driven by S brasiliensis in Brazil represents several, independent emergence events followed by animal-to-animal and animal-to human transmission within and between Brazilian states. These results describe how S brasiliensis can emerge and spread within a country.

FUNDING

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, Brazil; the São Paulo Research Foundation; Productivity in Research fellowships by the National Council for Scientific and Technological Development, and Ministry of Science and Technology of Brazil.

Genomic epidemiology and antifungal-resistant characterization of Candida auris, Colombia, 2016-2021

Since 2016, in Colombia, ongoing transmission of Candida auris has been reported in multiple cities. Here, we provide an updated description of C. auris genomic epidemiology and the dynamics of antifungal resistance in Colombia. We sequenced 99 isolates from C. auris cases with collection dates ranging from June 2016 to January 2021; the resulting sequences coupled with 103 previously generated sequences from C. auris cases were described in a phylogenetic analysis. All C. auris cases were clade IV. Of the 182 isolates with antifungal susceptibility data, 67 (37%) were resistant to fluconazole, and 39 (21%) were resistant to amphotericin B. Isolates predominately clustered by country except for 16 isolates from Bogotá, Colombia, which grouped with isolates from Venezuela. The largest cluster (N = 166 isolates) contained two subgroups. The first subgroup contained 26 isolates, mainly from César; of these, 85% (N = 22) were resistant to fluconazole. The second subgroup consisted of 47 isolates from the north coast; of these, 81% (N = 38) were resistant to amphotericin B. Mutations in the ERG11 and TAC1B genes were identified in fluconazole-resistant isolates. This work describes molecular mechanisms associated with C. auris antifungal resistance in Colombia. Overall, C. auris cases from different geographic locations in Colombia exhibited high genetic relatedness, suggesting continued transmission between cities since 2016. These findings also suggest a lack of or minimal introductions of different clades of C. auris into Colombia.

IMPORTANCE

Candida auris is an emerging fungus that presents a serious global health threat and has caused multiple outbreaks in Colombia. This work discusses the likelihood of introductions and local transmission of C. auris and provides an updated description of the molecular mechanisms associated with antifungal resistance in Colombia. Efforts like this provide information about the evolving C. auris burden that could help guide public health strategies to control C. auris spread.

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

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

Genomic description of human clinical Aspergillus fumigatus isolates, California, 2020

Aspergillus fumigatus, an environmental mold, causes life-threatening infections. Studies on the phylogenetic structure of human clinical A. fumigatus isolates are limited. Here, we performed whole genome sequencing of 24 A. fumigatus isolates collected from 18 patients in U.S. healthcare facilities in California. Single-nucleotide polymorphism (SNP) differences between patient isolates ranged from 187 to 70 829 SNPs. For five patients with multiple isolates, we calculated the within-host diversities. Three patients had a within-host diversity that ranged from 4 to 10 SNPs and two patients ranged from 2 to 16 977 SNPs. Findings revealed highly diverse A. fumigatus strains among patients and two patterns of diversity for isolates that come from the same patient, low and extremely high diversity.

Investigation of a Candida auris outbreak in a Skilled Nursing Facility — Virginia, United States, October 2020–June 2021

Candida auris, an emerging multi-drug resistant organism, is an urgent public health threat. We report on a C. auris outbreak investigation at a Virginia ventilator skilled nursing facility. During October 2020-June 2021, we identified 28 cases among residents in the ventilator unit. Genomic evidence suggested ≥2 distinct C. auris introductions to the facility. We identified multiple infection and prevention control challenges, highlighting the importance of strengthening multi-drug resistant organism prevention efforts at ventilator skilled nursing facilities.

P473 Genomic epidemiology of antifungal-resistant Candida auris in Colombia

Poster session 1, September 21, 2022, 12:30 PM - 1:30 PM

 

Introduction

Candida auris is a public health threat. Five major clades of C. auris have been identified (Clades I–V). In Colombia, C. auris infections were first reported in 2016 with ongoing transmission reported from multiple cities. Here, we describe C. auris genomic epidemiology in Colombia detailing cases from 2016–2021.

Methods

A total of 99 isolates from C. auris cases were collected between June 2016 to January 2021 in Colombia, representing 11 geographic locations. Species confirmation, antifungal susceptibility testing, and whole-genome sequencing (WGS) were performed. In all, 37 genomic sequences generated previously from isolates from C. auris cases in Colombia, Venezuela, Panama, Israel, and United States were also analyzed MycoSNP workflow was used to assess sequence quality, map reads to the reference, and identify single-nucleotide polymorphisms (SNPs). Pairwise distances and a neighbor-joining tree were generated. IQtree was used to generate a maximum-likelihood tree with bootstrap values.

Results

Phylogenetic analysis identified 1 493 SNP positions. Isolates from Colombia clustered to Clade IV and predominately grouped by country except for 16 fluconazole-resistant isolates from Bogota, Colombia that grouped with five isolates from Venezuela. In this cluster, 20 (95%) were resistant to fluconazole and 5 (24%) were resistant to fluconazole and the echinocandin micafungin. Remaining isolates from Bogota did not group in this cluster and were susceptible to fluconazole and micafungin.

A total of 98 isolates from Colombia clustered together. Within this Colombian cluster, there were two subgroups that had bootstrap support of 100% and were separated by 13 SNPs. The first subgroup was a cluster that contained 18 isolates from the north coast; 17 (94%) isolates were resistant to amphotericin B. A second subgroup consisted of 26 isolates from Cesar and Norte de Santander, and 22 (84%) isolates were resistant to fluconazole.

Conclusions

Based on the phylogenetic reconstruction, C. auris in Colombia continues to be of Clade IV. Amphotericin B-resistant isolates were predominately from the north coast, fluconazole-resistant isolates were from a wider geographic area in Colombia, and echinocandin-resistant isolates were from Bogota. Within the Colombian cluster comprising two subgroups, we observed high genetic relatedness between isolates from different geographic locations suggesting transmission among cities.

MycoSNP: A Portable Workflow for Performing Whole-Genome Sequencing Analysis of Candida auris

Candida auris is an urgent public health threat characterized by high drug-resistant rates and rapid spread in healthcare settings worldwide. As part of the C. auris response, molecular surveillance has helped public health officials track the global spread and investigate local outbreaks. Here, we describe whole-genome sequencing analysis methods used for routine C. auris molecular surveillance in the United States; methods include reference selection, reference preparation, quality assessment and control of sequencing reads, read alignment, and single-nucleotide polymorphism calling and filtration. We also describe the newly developed pipeline MycoSNP, a portable workflow for performing whole-genome sequencing analysis of fungal organisms including C. auris.

Candida auris Whole-Genome Sequence Benchmark Dataset for Phylogenomic Pipelines

Candida auris is a multidrug-resistant pathogen that represents a serious public health threat due to its rapid global emergence, increasing incidence of healthcare-associated outbreaks, and high rates of antifungal resistance. Whole-genome sequencing and genomic surveillance have the potential to bolster C. auris surveillance networks moving forward. Laboratories conducting genomic surveillance need to be able to compare analyses from various national and international surveillance partners to ensure that results are mutually trusted and understood. Therefore, we established an empirical outbreak benchmark dataset consisting of 23 C. auris genomes to help validate comparisons of genomic analyses and facilitate communication among surveillance networks. Our outbreak benchmark dataset represents a polyclonal phylogeny with three subclades. The genomes in this dataset are from well-vetted studies that are supported by multiple lines of evidence, which demonstrate that the whole-genome sequencing data, phylogenetic tree, and epidemiological data are all in agreement. This C. auris benchmark set allows for standardized comparisons of phylogenomic pipelines, ultimately promoting effective C. auris collaborations.

Hypertension and the roles of the 9p21.3 risk locus: Classic findings and new association data

BACKGROUND

The band 9p21.3 contains an established genomic risk zone for cardiovascular disease (CVD). Since the initial 2007 Wellcome Trust Case Control Consortium study (WTCCC), the increased CVD risk associated with 9p21.3 has been confirmed by multiple studies in different continents. However, many years later there was still no confirmed report of a corresponding association of 9p21.3 with hypertension, a major CV risk factor, nor with blood pressure (BP).

THEORY

In this contribution, we review the bipartite haplotype structure of the 9p21.3 risk locus: one block is devoid of protein-coding genes but contains the lead CVD risk SNPs, while the other block contains the first exon and regulatory DNA of the gene for the cell cycle inhibitor p15. We consider how findings from molecular biology offer possibilities of an involvement of p15 in hypertension etiology, with expression of the p15 gene modulated by genetic variation from within the 9p21.3 risk locus.

RESULTS

We present original results from a Colombian study revealing moderate but persistent association signals for BP and hypertension within the classic 9p21.3 CVD risk locus. These SNPs are mostly confined to a 'hypertension island' that spans less than 60 kb and coincides with the p15 haplotype block. We find confirmation in data originating from much larger, recent European BP studies, albeit with opposite effect directions.

CONCLUSION

Although more work will be needed to elucidate possible mechanisms, previous findings and new data prompt reconsidering the question of how variation in 9p21.3 might influence hypertension components of cardiovascular risk.

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.

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.

The LUFS domain, its transcriptional regulator proteins, and drug resistance in the fungal pathogen Candida auris

The LUFS domain (LUG/LUH, Flo8, single-strand DNA-binding protein [SSBP]) is a well-conserved and apparently ancient region found in diverse proteins and taxa. This domain, which has as its most obvious structural feature a series of three helices, has been identified in transcriptional regulator proteins of animals, plants, and fungi. Recently, in these pages (Wang et al., Protein Sci., 2019, 28:788-793), the first crystal structure of a LUFS domain was reported, for the human SSBP2, a transcriptional repressor. We briefly address how the new insights into LUFS structures might contribute to a better understanding of an important transcriptional activator of yeasts that contains the LUFS domain, Flo8, and consider how a focus on the LUFS domain and its variation could help us to understand etiologies of drug resistance in a recently emerged pathogenic fungus, Candida auris.

Molecular Epidemiology of Candida auris in Colombia Reveals a Highly Related, Countrywide Colonization With Regional Patterns in Amphotericin B Resistance

Background

Candida auris is a multidrug-resistant yeast associated with hospital outbreaks worldwide. During 2015-2016, multiple outbreaks were reported in Colombia. We aimed to understand the extent of contamination in healthcare settings and to characterize the molecular epidemiology of C. auris in Colombia.

Methods

We sampled patients, patient contacts, healthcare workers, and the environment in 4 hospitals with recent C. auris outbreaks. Using standardized protocols, people were swabbed at different body sites. Patient and procedure rooms were sectioned into 4 zones and surfaces were swabbed. We performed whole-genome sequencing (WGS) and antifungal susceptibility testing (AFST) on all isolates.

Results

Seven of the 17 (41%) people swabbed were found to be colonized. Candida auris was isolated from 37 of 322 (11%) environmental samples. These were collected from a variety of items in all 4 zones. WGS and AFST revealed that although isolates were similar throughout the country, isolates from the northern region were genetically distinct and more resistant to amphotericin B (AmB) than the isolates from central Colombia. Four novel nonsynonymous mutations were found to be significantly associated with AmB resistance.

Conclusions

Our results show that extensive C. auris contamination can occur and highlight the importance of adherence to appropriate infection control practices and disinfection strategies. Observed genetic diversity supports healthcare transmission and a recent expansion of C. auris within Colombia with divergent AmB susceptibility.

Genome Diversity, Recombination, and Virulence across the Major Lineages of Paracoccidioides

The Paracoccidioides genus includes two species of thermally dimorphic fungi that cause paracoccidioidomycosis, a neglected health-threatening human systemic mycosis endemic to Latin America. To examine the genome evolution and the diversity of Paracoccidioides spp., we conducted whole-genome sequencing of 31 isolates representing the phylogenetic, geographic, and ecological breadth of the genus. These samples included clinical, environmental and laboratory reference strains of the S1, PS2, PS3, and PS4 lineages of P. brasiliensis and also isolates of Paracoccidioides lutzii species. We completed the first annotated genome assemblies for the PS3 and PS4 lineages and found that gene order was highly conserved across the major lineages, with only a few chromosomal rearrangements. Comparing whole-genome assemblies of the major lineages with single-nucleotide polymorphisms (SNPs) predicted from the remaining 26 isolates, we identified a deep split of the S1 lineage into two clades we named S1a and S1b. We found evidence for greater genetic exchange between the S1b lineage and all other lineages; this may reflect the broad geographic range of S1b, which is often sympatric with the remaining, largely geographically isolated lineages. In addition, we found evidence of positive selection for the GP43 and PGA1 antigen genes and genes coding for other secreted proteins and proteases and lineage-specific loss-of-function mutations in cell wall and protease genes; these together may contribute to virulence and host immune response variation among natural isolates of Paracoccidioides spp. These insights into the recent evolutionary events highlight important differences between the lineages that could impact the distribution, pathogenicity, and ecology of Paracoccidioides. IMPORTANCE Characterization of genetic differences between lineages of the dimorphic human-pathogenic fungus Paracoccidioides can identify changes linked to important phenotypes and guide the development of new diagnostics and treatments. In this article, we compared genomes of 31 diverse isolates representing the major lineages of Paracoccidioides spp. and completed the first annotated genome sequences for the PS3 and PS4 lineages. We analyzed the population structure and characterized the genetic diversity among the lineages of Paracoccidioides, including a deep split of S1 into two lineages (S1a and S1b), and differentiated S1b, associated with most clinical cases, as the more highly recombining and diverse lineage. In addition, we found patterns of positive selection in surface proteins and secreted enzymes among the lineages, suggesting diversifying mechanisms of pathogenicity and adaptation across this species complex. These genetic differences suggest associations with the geographic range, pathogenicity, and ecological niches of Paracoccidioides lineages.

From NGS assembly challenges to instability of fungal mitochondrial genomes: A case study in genome complexity

The presence of repetitive or non-unique DNA persisting over sizable regions of a eukaryotic genome can hinder the genome's successful de novo assembly from short reads: ambiguities in assigning genome locations to the non-unique subsequences can result in premature termination of contigs and thus overfragmented assemblies. Fungal mitochondrial (mtDNA) genomes are compact (typically less than 100 kb), yet often contain short non-unique sequences that can be shown to impede their successful de novo assembly in silico. Such repeats can also confuse processes in the cell in vivo. A well-studied example is ectopic (out-of-register, illegitimate) recombination associated with repeat pairs, which can lead to deletion of functionally important genes that are located between the repeats. Repeats that remain conserved over micro- or macroevolutionary timescales despite such risks may indicate functionally or structurally (e.g., for replication) important regions. This principle could form the basis of a mining strategy for accelerating discovery of function in genome sequences. We present here our screening of a sample of 11 fully sequenced fungal mitochondrial genomes by observing where exact k-mer repeats occurred several times; initial analyses motivated us to focus on 17-mers occurring more than three times. Based on the diverse repeats we observe, we propose that such screening may serve as an efficient expedient for gaining a rapid but representative first insight into the repeat landscapes of sparsely characterized mitochondrial chromosomes. Our matching of the flagged repeats to previously reported regions of interest supports the idea that systems of persisting, non-trivial repeats in genomes can often highlight features meriting further attention.

The Dynamic Genome and Transcriptome of the Human Fungal Pathogen Blastomyces and Close Relative Emmonsia

Three closely related thermally dimorphic pathogens are causal agents of major fungal diseases affecting humans in the Americas: blastomycosis, histoplasmosis and paracoccidioidomycosis. Here we report the genome sequence and analysis of four strains of the etiological agent of blastomycosis, Blastomyces, and two species of the related genus Emmonsia, typically pathogens of small mammals. Compared to related species, Blastomyces genomes are highly expanded, with long, often sharply demarcated tracts of low GC-content sequence. These GC-poor isochore-like regions are enriched for gypsy elements, are variable in total size between isolates, and are least expanded in the avirulent B. dermatitidis strain ER-3 as compared with the virulent B. gilchristii strain SLH14081. The lack of similar regions in related species suggests these isochore-like regions originated recently in the ancestor of the Blastomyces lineage. While gene content is highly conserved between Blastomyces and related fungi, we identified changes in copy number of genes potentially involved in host interaction, including proteases and characterized antigens. In addition, we studied gene expression changes of B. dermatitidis during the interaction of the infectious yeast form with macrophages and in a mouse model. Both experiments highlight a strong antioxidant defense response in Blastomyces, and upregulation of dioxygenases in vivo suggests that dioxide produced by antioxidants may be further utilized for amino acid metabolism. We identify a number of functional categories upregulated exclusively in vivo, such as secreted proteins, zinc acquisition proteins, and cysteine and tryptophan metabolism, which may include critical virulence factors missed before in in vitro studies. Across the dimorphic fungi, loss of certain zinc acquisition genes and differences in amino acid metabolism suggest unique adaptations of Blastomyces to its host environment. These results reveal the dynamics of genome evolution and of factors contributing to virulence in Blastomyces.

Dimorphic fungal pathogens including Blastomyces are the cause of major fungal diseases in North and South America. The genus Emmonsia includes species infecting small mammals as well as a newly emerging pathogenic species recently reported in HIV-positive patients in South Africa. Here, we synthesize both genome sequencing of four isolates of Blastomyces and two species of Emmonsia as well as deep sequencing of Blastomyces RNA to draw major new insights into the evolution of this group and the pathogen response to infection. We investigate the trajectory of genome evolution of this group, characterizing the phylogenetic relationships of these species, a remarkable genome expansion that formed large isochore-like regions of low GC content in Blastomyces, and variation of gene content, related to host interaction, among the dimorphic fungal pathogens. Using RNA-Seq, we profile the response of Blastomyces to macrophage and mouse pulmonary infection, identifying key pathways and novel virulence factors. The identification of key fungal genes involved in adaptation to the host suggests targets for further study and therapeutic intervention in Blastomyces and related dimorphic fungal pathogens.

Genome update of the dimorphic human pathogenic fungi causing paracoccidioidomycosis

Paracoccidiodomycosis (PCM) is a clinically important fungal disease that can acquire serious systemic forms and is caused by the thermodimorphic fungal Paracoccidioides spp. PCM is a tropical disease that is endemic in Latin America, where up to ten million people are infected; 80% of reported cases occur in Brazil, followed by Colombia and Venezuela. To enable genomic studies and to better characterize the pathogenesis of this dimorphic fungus, two reference strains of P. brasiliensis (Pb03, Pb18) and one strain of P. lutzii (Pb01) were sequenced [1]. While the initial draft assemblies were accurate in large scale structure and had high overall base quality, the sequences had frequent small scale defects such as poor quality stretches, unknown bases (N's), and artifactual deletions or nucleotide duplications, all of which caused larger scale errors in predicted gene structures. Since assembly consensus errors can now be addressed using next generation sequencing (NGS) in combination with recent methods allowing systematic assembly improvement, we re-sequenced the three reference strains of Paracoccidioides spp. using Illumina technology. We utilized the high sequencing depth to re-evaluate and improve the original assemblies generated from Sanger sequence reads, and obtained more complete and accurate reference assemblies. The new assemblies led to improved transcript predictions for the vast majority of genes of these reference strains, and often substantially corrected gene structures. These include several genes that are central to virulence or expressed during the pathogenic yeast stage in Paracoccidioides and other fungi, such as HSP90, RYP1-3, BAD1, catalase B, alpha-1,3-glucan synthase and the beta glucan synthase target gene FKS1. The improvement and validation of these reference sequences will now allow more accurate genome-based analyses. To our knowledge, this is one of the first reports of a fully automated and quality-assessed upgrade of a genome assembly and annotation for a non-model fungus.

The fungal genus Paracoccidioides is the causal agent of paracoccidioidomycosis (PCM), a neglected tropical disease that is endemic in several countries of South America. Paracoccidioides is a pathogenic dimorphic fungus that is capable of converting to a virulent yeast form after inhalation by the host. Therefore the molecular biology of the switch to the yeast phase is of particular interest for understanding the virulence of this and other human pathogenic fungi, and ultimately for reducing the morbidity and mortality caused by such fungal infections. We here present the strategy and methods we used to update and improve accuracy of three reference genome sequences of Paracoccidioides spp. utilizing state-of-the-art Illumina re-sequencing, assembly improvement, re-annotation, and quality assessment. The resulting improved genome resource should be of wide use not solely for advancing research on the genetics and molecular biology of Paracoccidioides and the closely related pathogenic species Histoplasma and Blastomyces, but also for fungal diagnostics based on sequencing or molecular assays, characterizing rapidly changing proteins that may be involved in virulence, SNP-based population analyses and other tasks that require high sequence accuracy. The genome update and underlying strategy and methods also serve as a proof of principle that could encourage similar improvements of other draft genomes.