Jose F. Muñoz

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.

Novel taxa of thermally dimorphic systemic pathogens in the Ajellomycetaceae (Onygenales).

Recent discoveries of novel systemic fungal pathogens with thermally dimorphic yeast‐like phases have challenged the current taxonomy of the Ajellomycetaceae, a family currently comprising the genera Blastomyces, Emmonsia, Emmonsiellopsis, Helicocarpus, Histoplasma, Lacazia and Paracoccidioides. Our morphological, phylogenetic and phylogenomic analyses demonstrated species relationships and their specific phenotypes, clarified generic boundaries and provided the first annotated genome assemblies to support the description of two new species. A new genus, Emergomyces, accommodates Emmonsia pasteuriana as type species, and the new species Emergomyces africanus, the aetiological agent of case series of disseminated infections in South Africa. Both species produce small yeast cells that bud at a narrow base at 37°C and lack adiaspores, classically associated with the genus Emmonsia. Another novel dimorphic pathogen, producing broad‐based budding cells at 37°C and occurring outside North America, proved to belong to the genus Blastomyces, and is described as Blastomyces percursus.

Inhibition of PbGP43 expression may suggest that gp43 is a virulence factor in Paracoccidioides brasiliensis

Glycoprotein gp43 is an immunodominant diagnostic antigen for paracoccidioidomycosis caused by Paracoccidioides brasiliensis. It is abundantly secreted in isolates such as Pb339. It is structurally related to beta-1,3-exoglucanases, however inactive. Its function in fungal biology is unknown, but it elicits humoral, innate and protective cellular immune responses; it binds to extracellular matrix-associated proteins. In this study we applied an antisense RNA (aRNA) technology and Agrobacterium tumefaciens-mediated transformation to generate mitotically stable PbGP43 mutants (PbGP43 aRNA) derived from wild type Pb339 to study its role in P. brasiliensis biology and during infection. Control PbEV was transformed with empty vector. Growth curve, cell vitality and morphology of PbGP43 aRNA mutants were indistinguishable from those of controls. PbGP43 expression was reduced 80–85% in mutants 1 and 2, as determined by real time PCR, correlating with a massive decrease in gp43 expression. This was shown by immunoblotting of culture supernatants revealed with anti-gp43 mouse monoclonal and rabbit polyclonal antibodies, and also by affinity-ligand assays of extracellular molecules with laminin and fibronectin. In vitro, there was significantly increased TNF-α production and reduced yeast recovery when PbGP43 aRNA1 was exposed to IFN-γ-stimulated macrophages, suggesting reduced binding/uptake and/or increased killing. In vivo, fungal burden in lungs of BALB/c mice infected with silenced mutant was negligible and associated with decreased lung ΙΛ−10 and IL-6. Therefore, our results correlated low gp43 expression with lower pathogenicity in mice, but that will be definitely proven when PbGP43 knockouts become available. This is the first study of gp43 using genetically modified P. brasiliensis.

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

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. ABSTRACT 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.

Involvement of the 90 kDa heat shock protein during adaptation of Paracoccidioides brasiliensis to different environmental conditions.

HSP90 is a molecular chaperone that participates in folding, stabilization, activation, and assembly of several proteins, all of which are key regulators in cell signaling. In dimorphic pathogenic fungi such as Paracoccidioides brasiliensis, the adaptation to a higher temperature, acid pH and oxidative stress, is an essential event for fungal survival and also for the establishing of the infectious process. To further understand the role of this protein, we used antisense RNA technology to generate a P. brasiliensis isolate with reduced PbHSP90 gene expression (PbHSP90-aRNA). Reduced expression of HSP90 decreased yeast cell viability during batch culture growth and increased susceptibility to acid pH environments and imposed oxidative stress. Also, PbHSP90-aRNA yeast cells presented reduced viability upon interaction with macrophages. The findings presented here suggest a protective role for HSP90 during adaptation to hostile environments, one that promotes survival of the fungus during host-pathogen interactions.

50 Years of Emmonsia Disease in Humans: The Dramatic Emergence of a Cluster of Novel Fungal Pathogens

CITATION: Schwartz, I. S. et al. 2015. 50 years of Emmonsia disease in humans : the dramatic emergence of a cluster of novel fungal pathogens. PLoS Pathogens, 11(11): e1005198, doi:10.1371/journal.ppat.1005198.

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 (Ns), 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.

Identification and Analysis of the Role of Superoxide Dismutases Isoforms in the Pathogenesis of Paracoccidioides spp.

The ability of Paracoccidioides to defend itself against reactive oxygen species (ROS) produced by host effector cells is a prerequisite to survive. To counteract these radicals, Paracoccidioides expresses, among different antioxidant enzymes, superoxide dismutases (SODs). In this study, we identified six SODs isoforms encoded by the Paracoccidioides genome. We determined gene expression levels of representative isolates of the phylogenetic lineages of Paracoccidioides spp. (S1, PS2, PS3 and Pb01-like) using quantitative RT-PCR. Assays were carried out to analyze SOD gene expression of yeast cells, mycelia cells, the mycelia-to-yeast transition and the yeast-to-mycelia germination, as well as under treatment with oxidative agents and during interaction with phagocytic cells. We observed an increased expression of PbSOD1 and PbSOD3 during the transition process, exposure to oxidative agents and interaction with phagocytic cells, suggesting that these proteins could assist in combating the superoxide radicals generated during the host-pathogen interaction. Using PbSOD1 and PbSOD3 knockdown strains we showed these genes are involved in the response of the fungus against host effector cells, particularly the oxidative stress response, and in a mouse model of infection. Protein sequence analysis together with functional analysis of knockdown strains seem to suggest that PbSOD3 expression is linked with a pronounced extracellular activity while PbSOD1 seems more related to intracellular requirements of the fungus. Altogether, our data suggests that P. brasiliensis actively responds to the radicals generated endogenously during metabolism and counteracts the oxidative burst of immune cells by inducing the expression of SOD isoforms.

Paracoccidioides brasiliensis PbP27 gene: knockdown procedures and functional characterization

Paracoccidioides brasiliensis PbP27 gene encodes a protein localized in both the fungal cytoplasm and cell wall. The parasitic infectious form produces this protein preferentially with the genes expression varying between the fungus phylogenetic species. The biological function of the native p27 has yet to be determined during either growth of the yeast or host infection. Therefore, in this study, through the use of antisense RNA technology and Agrobacterium tumefaciens-mediated transformation, we generated mitotically stable PbP27 mutants (PbP27 aRNA) with the goal to evaluate the role of p27 in the biology and virulence of this fungus. PbP27 expression was reduced 60-75% in mutants, as determined by real-time PCR in correlation with a decrease in p27 expression. No alterations in the growth curve or in the ability to shift from mycelia to yeast or from yeast to mycelia were observed in PbP27 aRNA strains; however, we did observe a reduction in cell vitality. Moreover, a decrease in cell viability of PbP27 aRNA yeast cells after interaction with IFN-γ-stimulated macrophages was detected. Based on these results, we propose that p27 plays a role in yeast cell architecture and represents one of the mechanisms employed by this fungus for its interaction with the monocyte/macrophage system.

Paracoccidioides spp. catalases and their role in antioxidant defense against host defense responses

Dimorphic human pathogenic fungi interact with host effector cells resisting their microbicidal mechanisms. Yeast cells are able of surviving within the tough environment of the phagolysosome by expressing an antioxidant defense system that provides protection against host-derived reactive oxygen species (ROS). This includes the production of catalases (CATs). Here we identified and analyzed the role of CAT isoforms in Paracoccidioides, the etiological agent of paracoccidioidomycosis. Firstly, we found that one of these isoforms was absent in the closely related dimorphic pathogen Coccidioides and dermatophytes, but all of them were conserved in Paracoccidioides, Histoplasma and Blastomyces species. We probed the contribution of CATs in Paracoccidioides by determining the gene expression levels of each isoform through quantitative RT-qPCR, in both the yeast and mycelia phases, and during the morphological switch (transition and germination), as well as in response to oxidative agents and during interaction with neutrophils. PbCATP was preferentially expressed in the pathogenic yeast phase, and was associated to the response against exogenous H2O2. Therefore, we created and analyzed the virulence defects of a knockdown strain for this isoform, and found that CATP protects yeast cells from H2O2 generated in vitro and is relevant during lung infection. On the other hand, CATA and CATB seem to contribute to ROS homeostasis in Paracoccidioides cells, during endogenous oxidative stress. CAT isoforms in Paracoccidioides might be coordinately regulated during development and dimorphism, and differentially expressed in response to different stresses to control ROS homeostasis during the infectious process, contributing to the virulence of Paracoccidioides.