Genome Organization and Copy-Number Variation Reveal Clues to Virulence Evolution in Coccidioides posadasii.

The human fungal pathogen Coccidioides spp. causes valley fever, a treatment-refractory and sometimes deadly disease prevalent in arid regions of the western hemisphere. Fungal virulence in the mammalian host hinges on a switch between growth as hyphae and as large spherules containing infectious spores. How these virulence programs are encoded in the genome remains poorly understood. Drawing on Coccidioides genomic resources, we first discovered a new facet of genome organization in this system: spherule-gene islands, clusters of genes physically linked in the genome that exhibited specific mRNA induction in the spherule phase. Next, we surveyed copy-number variation genome-wide among strains of C. posadasii. Emerging from this catalog were spherule-gene islands with striking presence-absence differentiation between C. posadasii populations, a pattern expected from virulence factors subjected to different selective pressures across habitats. Finally, analyzing single-nucleotide differences across C. posadasii strains, we identified signatures of natural selection in spherule-expressed genes. Together, our data establish spherule-gene islands as candidate determinants of virulence and targets of selection in Coccidioides.

Single nucleotide polymorphisms and chromosomal copy number variation may impact the Sporothrix brasiliensis antifungal susceptibility and sporotrichosis clinical outcomes.

Feline-transmitted sporotrichosis has garnered attention due to the recent high incidence and the lack of efficient control in the epicenter of the epidemic, Rio de Janeiro, Brazil. Sporothrix brasiliensis is the major pathogen involved in feline-to-human sporotrichosis in Brazil and displays more virulent genotypes than the closely related species S. schenckii. Over the last two decades, several reports of antifungal-resistant strains have emerged. Sequencing and comparison analysis of the outbreak strains allowed us to observe that the azole non-wild-type S. brasiliensis strain CFP 1054 had significant chromosomal variations compared to wild-type strains. One of these variants includes a region of 231 Kb containing 75 duplicated genes, which were overrepresented for lipid and isoprenoid metabolism. We also identified an additional strain (CFP 1055) that was resistant to itraconazole and amphotericin B, which had a single nucleotide polymorphism in the tac1 gene. The patients infected with these two strains showed protracted clinical course and sequelae. Even though our sample size is modest, these results suggest the possibility of identifying specific point mutations and large chromosomal duplications potentially associated with antifungal resistance and clinical outcomes of sporotrichosis.

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.

Paracoccidioides Genomes Reflect High Levels of Species Divergence and Little Interspecific Gene Flow.

The fungus Paracoccidioides is a prevalent human pathogen endemic to South America. The genus is composed of five species. In this report, we use 37 whole-genome sequences to study the allocation of genetic variation in Paracoccidioides We tested three genome-wide predictions of advanced speciation, namely, that all species should be reciprocally monophyletic, that species pairs should be highly differentiated along the whole genome, and that there should be low rates of interspecific gene exchange. We find support for these three hypotheses. Species pairs with older divergences show no evidence of gene exchange, while more recently diverged species pairs show evidence of modest rates of introgression. Our results indicate that as divergence progresses, species boundaries become less porous among Paracoccidioides species. Our results suggest that species in Paracoccidioides are at different stages along the divergence continuum.IMPORTANCEParacoccidioides is the causal agent of a systemic mycosis in Latin America. Most of the inference of the evolutionary history of Paracoccidioides has used only a few molecular markers. In this report, we evaluate the extent of genome divergence among Paracoccidioides species and study the possibility of interspecific gene exchange. We find that all species are highly differentiated. We also find that the amount of gene flow between species is low and in some cases is even completely absent in spite of geographic overlap. Our study constitutes a systematic effort to identify species boundaries in fungal pathogens and to determine the extent of gene exchange among fungal species.

Investigating the Role of Animal Burrows on the Ecology and Distribution of Coccidioides spp. in Arizona Soils.

The lack of knowledge regarding the ecology of Coccidioides spp. makes both modeling the potential for disease outbreaks and predicting the distribution of the organism in the environment challenging. No single ecological parameter explains the biogeography of the pathogen. Previous investigations suggest an association with desert mammals, but these results should be confirmed with modern molecular techniques. Therefore, we used molecular tools to analyze soils associated with animal activity (i.e., burrows) to better define the ecology and biogeography of Coccidioides spp. in Arizona. Soils were collected from locations predicted to have favorable habitat outside of the established endemic regions to better understand the ecological niche of the organism in this state. Our central hypothesis is that soils taken from within animal burrows will have a higher abundance of Coccidioides spp. when compared to soils not directly associated with animal burrows. Our results show that there is a positive relationship with Coccidioides spp. and animal burrows. The organism was detected in two locations in northern Arizona at sites not known previously to harbor the fungus. Moreover, this fungus is able to grow on keratinized tissues (i.e., horse hair). These results provide additional evidence that there is a relationship between Coccidioides spp. and desert animals, which sheds new light on Coccidioides' ecological niche. These results also provide evidence that the geographic range of the organism may be larger than previously thought, and the concept of endemicity should be reevaluated for Coccidioides.

Population Structure and Genetic Diversity among Isolates of Coccidioides posadasii in Venezuela and Surrounding Regions.

Coccidioides posadasii is a pathogenic fungus that causes coccidioidomycosis in many arid regions of the Americas. One of these regions is bordered by the Caribbean Sea, and the surrounding landscape may play an important role in the dispersion of C. posadasii across South America through southeastern Mexico, Honduras, Guatemala, and Venezuela. Comparative phylogenomic analyses of C. posadasii reveal that clinical strains from Venezuela are genetically distinct from the North American populations found in (i) Arizona and (ii) Texas, Mexico, and the rest of South America (TX/MX/SA). We find evidence for admixture between the Venezuela and the North American populations of C. posadasii in Central America. Additionally, the proportion of Venezuelan alleles in the admixed population decreases as latitude (and distance from Venezuela) increases. Our results indicate that the population in Venezuela may have been subjected to a recent bottleneck and shows a strong population structure. This analysis provides insight into potential for Coccidioides spp. to invade new regions.IMPORTANCE Valley Fever is a fungal disease caused by two species of fungi: Coccidioides immitis and C. posadasii These fungi are found throughout the arid regions of North and South America; however, our understanding of genetic diversity and disease in South America is limited. In this report, we analyze 10 new genomes of Coccidioides posadasii from regions bordering the Caribbean Sea. We show that these populations are distinct and that isolates from Venezuela are likely a result of a recent bottleneck. These data point to patterns that might be observed when investigating recently established populations.

Fantastic yeasts and where to find them: the hidden diversity of dimorphic fungal pathogens.

Dimorphic fungal pathogens are a significant cause of human disease worldwide. Notably, the dimorphic fungal pathogens within the order Onygenales are considered primary pathogens, causing disease in healthy hosts. Current changes in taxonomy are underway due to advances in molecular phylogenetics, population genetics, and new emerging dimorphic fungal pathogens causing human disease. In this review, we highlight evolutionary relationships of dimorphic fungal pathogens that cause human disease within the order Onygenales and provide rationale to support increased investment in studies understanding the evolutionary relationships of these pathogens to improve rapid diagnostics, help identify mechanisms of antifungal resistance, understand adaptation to human host, and factors associated with virulence.

Gene exchange between two divergent species of the fungal human pathogen, Coccidioides.

The fungal genus Coccidioides is composed of two species, Coccidioides immitis and Coccidioides posadasii. These two species are the causal agents of coccidioidomycosis, a pulmonary disease also known as valley fever. The two species are thought to have shared genetic material due to gene exchange in spite of their long divergence. To quantify the magnitude of shared ancestry between them, we analyzed the genomes of a population sample from each species. Next, we inferred what is the expected size of shared haplotypes that might be inherited from the last common ancestor of the two species and find a cutoff to find what haplotypes have conclusively been exchanged between species. Finally, we precisely identified the breakpoints of the haplotypes that have crossed the species boundary and measure the allele frequency of each introgression in this sample. We find that introgressions are not uniformly distributed across the genome. Most, but not all, of the introgressions segregate at low frequency. Our results show that divergent species can share alleles, that species boundaries can be porous, and highlight the need for a systematic exploration of gene exchange in fungal species.

Population structure and genetic diversity among isolates of coccidioides posadasii in Venezuela and surrounding regions

Coccidioides posadasii is a pathogenic fungus that causes coccidioidomycosis in many arid regions of the Americas. One of these regions is bordered by the Caribbean Sea, and the surrounding landscape may play an important role in the dispersion of C. posadasii across South America through southeastern Mexico, Honduras, Guatemala, and Venezuela. Comparative phylogenomic analyses of C. posadasii reveal that clinical strains from Venezuela are genetically distinct from the North American populations found in (i) Arizona and (ii) Texas, Mexico, and the rest of South America (TX/MX/SA). We find evidence for admixture between the Venezuela and the North American populations of C. posadasii in Central America. Additionally, the proportion of Venezuelan alleles in the admixed population decreases as latitude (and distance from Venezuela) increases. Our results indicate that the population in Venezuela may have been subjected to a recent bottleneck and shows a strong population structure. This analysis provides insight into potential for Coccidioides spp. to invade new regions. IMPORTANCE Valley Fever is a fungal disease caused by two species of fungi: Coccidioides immitis and C. posadasii. These fungi are found throughout the arid regions of North and South America; however, our understanding of genetic diversity and disease in South America is limited. In this report, we analyze 10 new genomes of Coccidioides posadasii from regions bordering the Caribbean Sea. We show that these populations are distinct and that isolates from Venezuela are likely a result of a recent bottleneck. These data point to patterns that might be observed when investigating recently established populations.

Population structure and genetic diversity among isolates of coccidioides posadasii in Venezuela and surrounding regions

Coccidioides posadasii is a pathogenic fungus that causes coccidioidomycosis in many arid regions of the Americas. One of these regions is bordered by the Caribbean Sea, and the surrounding landscape may play an important role in the dispersion of C. posadasii across South America through southeastern Mexico, Honduras, Guatemala, and Venezuela. Comparative phylogenomic analyses of C. posadasii reveal that clinical strains from Venezuela are genetically distinct from the North American populations found in (i) Arizona and (ii) Texas, Mexico, and the rest of South America (TX/MX/SA). We find evidence for admixture between the Venezuela and the North American populations of C. posadasii in Central America. Additionally, the proportion of Venezuelan alleles in the admixed population decreases as latitude (and distance from Venezuela) increases. Our results indicate that the population in Venezuela may have been subjected to a recent bottleneck and shows a strong population structure. This analysis provides insight into potential for Coccidioides spp. to invade new regions. IMPORTANCE Valley Fever is a fungal disease caused by two species of fungi: Coccidioides immitis and C. posadasii. These fungi are found throughout the arid regions of North and South America; however, our understanding of genetic diversity and disease in South America is limited. In this report, we analyze 10 new genomes of Coccidioides posadasii from regions bordering the Caribbean Sea. We show that these populations are distinct and that isolates from Venezuela are likely a result of a recent bottleneck. These data point to patterns that might be observed when investigating recently established populations.

Species boundaries in the human pathogen Paracoccidioides.

The use of molecular taxonomy for identifying recently diverged species has transformed the study of speciation in fungi. The pathogenic fungus Paracoccidioides spp has been hypothesized to be composed of five phylogenetic species, four of which compose the brasiliensis species complex. Nuclear gene genealogies support this divergence scenario, but mitochondrial loci do not; while all species from the brasiliensis complex are differentiated at nuclear coding loci, they are not at mitochondrial loci. We addressed the source of this incongruity using 11 previously published gene fragments, 10 newly-sequenced nuclear non-coding loci, and 10 microsatellites. We hypothesized and further demonstrated that the mito-nuclear incongruence in the brasiliensis species complex results from interspecific hybridization and mitochondrial introgression, a common phenomenon in eukaryotes. Additional population genetic analyses revealed possible nuclear introgression but much less than that seen in the mitochondrion. Our results are consistent with a divergence scenario of secondary contact and subsequent mitochondrial introgression despite the continued persistence of species boundaries. We also suggest that yeast morphology slightly-but significantly-differs across all five Paracoccidioides species and propose to elevate four of these phylogenetic species to formally described taxonomic species.

Use of Population Genetics to Assess the Ecology, Evolution, and Population Structure of Coccidioides.

During the past 20 years, a general picture of the genetic diversity and population structure of Coccidioides, the causal agent of coccidioidomycosis (Valley fever), has emerged. The genus consists of 2 genetically diverse species, C. immitis and C. posadasii, each of which contains 1 or more distinct populations with limited gene flow. Genotypic data indicate that C. immitis is divided into 2 subpopulations (central and southern California populations) and C. posadasii is divided into 3 subpopulations (Arizona, Mexico, and Texas/South America populations). However, admixture within and among these populations and the current paucity of environmental isolates limit our understanding of the population genetics of Coccidioides. We assessed population structure of Coccidioides in Arizona by analyzing 495 clinical and environmental isolates. Our findings confirm the population structure as previously described and indicate a finer scale population structure in Arizona. Environmental isolates appear to have higher genetic diversity than isolates from human patients.

Functional characterization of Sporothrix schenckii glycosidases involved in the N-linked glycosylation pathway.

Protein glycosylation pathways are conserved metabolic processes in eukaryotic organisms and are required for cell fitness. In fungal pathogens, the N-linked glycosylation pathway is indispensable for proper cell wall composition and virulence. In Sporothrix schenckii sensu stricto, the causative agent of sporotrichosis, little is known about this glycosylation pathway. Here, using a genome-wide screening for putative members of the glycosyl hydrolase (CAZy - GH) families 47 and 63, which group enzymes involved in the processing step during N-linked glycan maturation, we found seven homologue genes belonging to family 47 and one to family 63. The eight genes were individually expressed in C. albicans null mutants lacking either MNS1 (for members of family 47) or CWH41 (for the member of family 63). Our results indicate that SsCWH41 is the functional ortholog of CaCWH41, whereas SsMNS1 is the functional ortholog of CaMNS1. The remaining genes of family 47 encode Golgi mannosidases and endoplasmic reticulum degradation-enhancing alpha-mannosidase-like proteins (EDEMs). Since these GH families gather proteins used as target for drugs to control cell growth, identification of these genes could help in the design of antifungals that could be used to treat sporotrichosis and other fungal diseases. In addition, to our knowledge, we are the first to report that Golgi mannosidases and EDEMs are expressed and characterized in yeast cells.

Comparative genomics of the major fungal agents of human and animal Sporotrichosis: Sporothrix schenckii and Sporothrix brasiliensis.

BACKGROUND: The fungal genus Sporothrix includes at least four human pathogenic species. One of these species, S. brasiliensis, is the causal agent of a major ongoing zoonotic outbreak of sporotrichosis in Brazil. Elsewhere, sapronoses are caused by S. schenckii and S. globosa. The major aims on this comparative genomic study are: 1) to explore the presence of virulence factors in S. schenckii and S. brasiliensis; 2) to compare S. brasiliensis, which is cat-transmitted and infects both humans and cats with S. schenckii, mainly a human pathogen; 3) to compare these two species to other human pathogens (Onygenales) with similar thermo-dimorphic behavior and to other plant-associated Sordariomycetes. RESULTS: The genomes of S. schenckii and S. brasiliensis were pyrosequenced to 17x and 20x coverage comprising a total of 32.3 Mb and 33.2 Mb, respectively. Pair-wise genome alignments revealed that the two species are highly syntenic showing 97.5% average sequence identity. Phylogenomic analysis reveals that both species diverged about 3.8-4.9 MYA suggesting a recent event of speciation. Transposable elements comprise respectively 0.34% and 0.62% of the S. schenckii and S. brasiliensis genomes and expansions of Gypsy-like elements was observed reflecting the accumulation of repetitive elements in the S. brasiliensis genome. Mitochondrial genomic comparisons showed the presence of group-I intron encoding homing endonucleases (HE's) exclusively in S. brasiliensis. Analysis of protein family expansions and contractions in the Sporothrix lineage revealed expansion of LysM domain-containing proteins, small GTPases, PKS type1 and leucin-rich proteins. In contrast, a lack of polysaccharide lyase genes that are associated with decay of plants was observed when compared to other Sordariomycetes and dimorphic fungal pathogens, suggesting evolutionary adaptations from a plant pathogenic or saprobic to an animal pathogenic life style. CONCLUSIONS: Comparative genomic data suggest a unique ecological shift in the Sporothrix lineage from plant-association to mammalian parasitism, which contributes to the understanding of how environmental interactions may shape fungal virulence. . Moreover, the striking differences found in comparison with other dimorphic fungi revealed that dimorphism in these close relatives of plant-associated Sordariomycetes is a case of convergent evolution, stressing the importance of this morphogenetic change in fungal pathogenesis.

Differences in cell morphometry, cell wall topography and gp70 expression correlate with the virulence of Sporothrix brasiliensis clinical isolates.

Sporotrichosis is a chronic infectious disease affecting both humans and animals. For many years, this subcutaneous mycosis had been attributed to a single etiological agent; however, it is now known that this taxon consists of a complex of at least four pathogenic species, including Sporothrix schenckii and Sporothrix brasiliensis. Gp70 was previously shown to be an important antigen and adhesin expressed on the fungal cell surface and may have a key role in immunomodulation and host response. The aim of this work was to study the virulence, morphometry, cell surface topology and gp70 expression of clinical isolates of S. brasiliensis compared with two reference strains of S. schenckii. Several clinical isolates related to severe human cases or associated with the Brazilian zoonotic outbreak of sporotrichosis were genotyped and clustered as S. brasiliensis. Interestingly, in a murine subcutaneous model of sporotrichosis, these isolates showed a higher virulence profile compared with S. schenckii. A single S. brasiliensis isolate from an HIV-positive patient not only showed lower virulence but also presented differences in cell morphometry, cell wall topography and abundant gp70 expression compared with the virulent isolates. In contrast, the highly virulent S. brasiliensis isolates showed reduced levels of cell wall gp70. These observations were confirmed by the topographical location of the gp70 antigen using immunoelectromicroscopy in both species. In addition, the gp70 molecule was sequenced and identified using mass spectrometry, and the sequenced peptides were aligned into predicted proteins using Blastp with the S. schenckii and S. brasiliensis genomes.

Polymorphism analysis of the CTLA-4 gene in paracoccidioidomycosis patients.

The CTLA-4 protein is expressed in activated T cells and plays an essential role in the immune response through its regulatory effect on T cell activation. Polymorphisms of the CTLA-4 gene have been correlated with autoimmune, neoplastic and infectious illnesses. This work aimed to verify possible associations between single nucleotide polymorphisms (SNPs) in CTLA-4, -318C/T in the promoter and +49A/G in exon 1 and paracoccidioidomycosis (PCM) caused by Paracoccidioides brasiliensis. For this purpose, 66 chronic form PCM patients and 76 healthy controls had their allele, genotype and haplotype frequencies determined. The genetic admixture structure of the patients and controls was evaluated to eliminate ancestral bias. The comparison of frequencies indicated no significant differences between patients and controls that could link the SNPs to PCM. Groups were admixture matched with no difference observed in population ancestry inference, indicating that the absence of association between CTLA-4 polymorphisms and PCM could not be attributed to ancestral bias. This study showed that there was no association between the CTLA-4 SNPs -318 and +49 and the resistance or susceptibility to PCM.

Polymorphism analysis of the CTLA-4 gene in paracoccidioidomycosis patients

The CTLA-4 protein is expressed in activated T cells and plays an essential role in the immune response through its regulatory effect on T cell activation. Polymorphisms of the CTLA-4 gene have been correlated with autoimmune, neoplastic and infectious illnesses. This work aimed to verify possible associations between single nucleotide polymorphisms (SNPs) in CTLA-4, -318C/T in the promoter and +49A/G in exon 1 and paracoccidioidomycosis (PCM) caused by Paracoccidioides brasiliensis. For this purpose, 66 chronic form PCM patients and 76 healthy controls had their allele, genotype and haplotype frequencies determined. The genetic admixture structure of the patients and controls was evaluated to eliminate ancestral bias. The comparison of frequencies indicated no significant differences between patients and controls that could link the SNPs to PCM. Groups were admixture matched with no difference observed in population ancestry inference, indicating that the absence of association between CTLA-4 polymorphisms and PCM could not be attributed to ancestral bias. This study showed that there was no association between the CTLA-4 SNPs -318 and +49 and the resistance or susceptibility to PCM.

Comparative genomics allowed the identification of drug targets against human fungal pathogens.

BACKGROUND: The prevalence of invasive fungal infections (IFIs) has increased steadily worldwide in the last few decades. Particularly, there has been a global rise in the number of infections among immunosuppressed people. These patients present severe clinical forms of the infections, which are commonly fatal, and they are more susceptible to opportunistic fungal infections than non-immunocompromised people. IFIs have historically been associated with high morbidity and mortality, partly because of the limitations of available antifungal therapies, including side effects, toxicities, drug interactions and antifungal resistance. Thus, the search for alternative therapies and/or the development of more specific drugs is a challenge that needs to be met. Genomics has created new ways of examining genes, which open new strategies for drug development and control of human diseases. RESULTS: In silico analyses and manual mining selected initially 57 potential drug targets, based on 55 genes experimentally confirmed as essential for Candida albicans or Aspergillus fumigatus and other 2 genes (kre2 and erg6) relevant for fungal survival within the host. Orthologs for those 57 potential targets were also identified in eight human fungal pathogens (C. albicans, A. fumigatus, Blastomyces dermatitidis, Paracoccidioides brasiliensis, Paracoccidioides lutzii, Coccidioides immitis, Cryptococcus neoformans and Histoplasma capsulatum). Of those, 10 genes were present in all pathogenic fungi analyzed and absent in the human genome. We focused on four candidates: trr1 that encodes for thioredoxin reductase, rim8 that encodes for a protein involved in the proteolytic activation of a transcriptional factor in response to alkaline pH, kre2 that encodes for α-1,2-mannosyltransferase and erg6 that encodes for Δ(24)-sterol C-methyltransferase. CONCLUSIONS: Our data show that the comparative genomics analysis of eight fungal pathogens enabled the identification of four new potential drug targets. The preferred profile for fungal targets includes proteins conserved among fungi, but absent in the human genome. These characteristics potentially minimize toxic side effects exerted by pharmacological inhibition of the cellular targets. From this first step of post-genomic analysis, we obtained information relevant to future new drug development.

Phylogenetic analysis reveals a high level of speciation in the Paracoccidioides genus.

Paracoccidioidomycosis (PCM) is a systemic disease endemic to most of Latin America, with greatest impact in rural areas. The taxonomic status of one of the best studied Paracoccidioides isolates (Pb01) as P. brasiliensis remains unresolved due to its genomic differences from the other three previously described phylogenetic species (S1, PS2 and PS3; Carrero et al., 2008. Fungal Genet. Biol. 45, 605). Using the genealogic concordance method of phylogenetic species recognition (GCPSR) via maximum parsimony and Bayesian analysis, we identified a clade of 17 genotypically similar isolates, including Pb01, which are distinct from the S1/PS2/P3 clade. Consistent with GCPSR, this "Pb01-like" group can be considered a new phylogenetic species, since it is strongly supported by all independent and concatenated genealogies. "Pb01-like" species exhibit great sequence and morphological divergence from the S1/PS2/PS3 species clade, and we estimate that these groups last shared a common ancestor approximately 32 million years ago. In addition, recombination analysis revealed independent events inside both main groups suggesting reproductive isolation. Consequently, we recommend the formal description of the "Pb01-like" cluster as the new species Paracoccidioides lutzii, a tribute to Adolpho Lutz, discoverer of P. brasiliensis in 1908.