Amanda Bombassaro1
Gabriela X. Schneider1
Flávia F. Costa2
Aniele C. R. Leão2
Bruna S. Soley3
Fernanda Medeiros4
Nickolas M. da Silva2
Bruna J. F. S. Lima1
Raffael J. A. Castro5
Anamélia L. Bocca5
Valter A. Baura6Eduardo Balsanelli6Vania C. S. Pankievicz6Nyvia M. C. Hrysay7
Rosana H. Scola7
Leandro F. Moreno2
Conceição M. P. S. Azevedo8
Emanuel M. Souza6
Renata R. Gomes1Sybren de Hoog1,9*
Vânia A. Vicente1,2*- 1Microbiology, Parasitology and Pathology Post-graduation Program, Department of Basic Pathology, Federal University of Paraná, Curitiba, Brazil
- 2Engineering Bioprocess and Biotechnology Post-graduation Program, Department of Bioprocess Engineering and Biotechnology, Federal University of Paraná, Curitiba, Brazil
- 3Pharmacology Post-graduation Program, Department of Pharmacology, Federal University of Paraná, Curitiba, Brazil
- 4Graduation in Biology Sciences, Federal University of Paraná, Curitiba, Brazil
- 5Department of Cell Biology, University of Brasilia, Brasilia, Brazil
- 6Department of Biochemistry, Federal University of Paraná, Curitiba, Brazil
- 7Service of Neuromuscular and Demyelinating Diseases, Complex Histochemistry-Immunity Laboratory, Hospital of Clinics, Federal University of Paraná, Curitiba, Brazil
- 8Department of Medicine, Federal University of Maranhão, São Luís, Brazil
- 9Center of Expertise in Mycology of Radboud University Medical Center/Canisius Wilhelmina Hospital, Nijmegen, Netherlands
Among agents of chromoblastomycosis, Fonsecaea pugnacius presents a unique type of infection because of its secondary neurotropic dissemination from a chronic cutaneous case in an immunocompetent patient. Neurotropism occurs with remarkable frequency in the fungal family Herpotrichiellaceae, possibly associated with the ability of some species to metabolize aromatic hydrocarbons. In an attempt to understand this new disease pattern, were conducted genomic analysis of Fonsecaea pugnacius (CBS 139214) performed with de novo assembly, gene prediction, annotation and mitochondrial genome assembly, supplemented with animal infection models performed with Tenebrio molitor in Mus musculus lineages BALB/c and C57BL/6. The genome draft of 34.8 Mb was assembled with a total of 12,217 protein-coding genes. Several proteins, enzymes and metabolic pathways related to extremotolerance and virulence were recognized. The enzyme profiles of black fungi involved in chromoblastomycosis and brain infection were analyzed with the Carbohydrate-Active Enzymes (CAZY) and peptidases database (MEROPS). The capacity of the fungus to survive inside Tenebrio molitor animal model was confirmed by histopathological analysis and by presence of melanin and hyphae in host tissue. Although F. pugnacius was isolated from brain in a murine model following intraperitoneal infection, cytokine levels were not statistically significant, indicating a profile of an opportunistic agent. A dual ecological ability can be concluded from presence of metabolic pathways for nutrient scavenging and extremotolerance, combined with a capacity to infect human hosts.
Introduction
Melanized fungi that are known as ‘black yeasts and relatives’ and belonging to the family Herpotrichiellaceae (order Chaetothyriales) are associated with different clinical pictures such as mycetoma, phaeohyphomycosis, and chromoblastomycosis (Cañete-Gibas and Wiederhold, 2018). Chromoblastomycosis starts at the inoculation site of the etiological agent, leads to chronic acanthosis, and, triggered by the host’s immune response, develops fungal structures known as muriform cells (Vicente et al., 2017). Phaeohyphomycosis can be distinguished from other infectious syndromes by tissue invasion with pigmented hyphae (Thomas et al., 2018) and is often associated with necrosis.
Significant differences in pathogenicity and virulence between the main clinical species Cladophialophora bantiana, Exophiala dermatitidis, Fonsecaea pedrosoi, and Rhinocladiella mackenziei on the one hand, and closely related environmental species on the other, have been reported. These species can grow at human body temperature or higher and are able to cause systemic or disseminated disease, while many others, if causing infection, remain subcutaneous (Seyedmousavi et al., 2014; Teixeira et al., 2017). Cerebral infections by black fungi are characterized by abscesses with hyphae in tissue and are therefore classified as phaeohyphomycosis. Such infections often lead to death of the patient despite combined treatment with antifungal drugs and surgery (de Azevedo et al., 2015; Arcobello and Revankar, 2020).
Infection of humans by members of Herpotrichiellaceae is enabled by their stress tolerance and adaptability in their natural, environmental niche, the human host not being a preferential habitat (Gostinčar et al., 2018). The agents are saprobes in mostly still unclarified micro-habitats, and decompose organic matter for nutrition (Vicente et al., 2017). The high adaptability and invasive potential explain the relatively high frequency in animal hosts despite a low environmental occurrence. This potential for infection appears to be polyphyletic within the family Herpotrichiellaceae, as it differs between species (Vicente et al., 2017).
Fonsecaea sibling species differ significantly in their ecology and potential of infection (Vicente et al., 2014). Some, such as Fonsecaea pedrosoi and F. nubica, are recognized etiologic agents of chromoblastomycosis in human hosts, having an ability to form muriform cell in tissue, while F. monophora and F. pugnacius may also be involved in disseminated infection with hyphal growth in the brain (de Hoog et al., 2004; Tanabe et al., 2004; Najafzadeh et al., 2009, 2011; Vicente et al., 2012; de Azevedo et al., 2015). Nine cases of primary brain infection by these species have been confirmed (Lucasse et al., 1954; Nobrega et al., 2003; Surash et al., 2005; Takei et al., 2007; Koo et al., 2010; Raparia et al., 2010; Doymaz et al., 2015; Varghese et al., 2016; Helbig et al., 2018). Fonsecaea pugnacius is exceptional by combining features of chromoblastomycosis and secondary neurotropic dissemination. The single strain known to date of the species presented muriform cells in subcutaneous tissue, but hyphae in the cerebrum (de Azevedo et al., 2015). This duality of local and invasive morphologies has not been observed in any other species associated with chromoblastomycosis or brain infection, suggesting that Fonsecaea pugnacius presents a unique pathogenic profile different from that of Cladophialophora bantiana, the main agent of human brain infection which presumably follows a pulmonary route (Ozgun et al., 2019).
In the present study, we sequenced the genome of the type strain of F. pugnacius, CBS 139214, and performed genomic analysis in order to identify the relation between black fungi and neurotropism. In addition, we compared the enzymatous gene profile of F. pugnacius to other previously sequenced neurotropic species, including C. bantiana, E. dermatitidis and R. mackenziei. To obtain more insight into virulence, we evaluated animal infection models by F. pugnacius, using strain CBS 139214 isolated from a cutaneous lesion of a patient with disseminated neurotropic infection.
Materials and Methods
Genomic DNA Extraction, Sequencing and Assembly
Fonsecaea pugnacius strain CBS 139214 (type) originating from a skin lesion (de Azevedo et al., 2015) was obtained from the reference collection of Westerdijk Fungal Biodiversity Institute, Utrecht, Netherlands. The strain was grown in Sabouraud liquid medium during 7 days at 28°C for DNA extraction according to Vicente et al. (2008) using cetyltrimethylammonium bromide (CTAB) and phenol-chloroform/isoamyl alcohol and the Microbial DNA UltraCleanTM kit for purification. The Nextera kit (IlluminaTM) and Ion Plus Fragment Library Kit (Thermo Fisher Scientific) were used to prepare the DNA libraries for sequencing based on the producer’s guidelines. FastQC1 was used for quality control analyses of sequence reads generated. The SPAdes assembler v3.10.0 (Bankevich et al., 2012) and FGAP (Piro et al., 2014) were applied for de novo assembly and gap closure, respectively. The genome assembly was evaluated by BUSCO v4.0.2 using the ‘chaetothyriales_odb10’ dataset (Seppey et al., 2019) and the Bowtie2 program was used for assembly coverage measure (Langmead and Salzberg, 2012).
Gene Prediction, Annotation and Genomic Analysis
GeneMark-ES v4.39 (Besemer et al., 2001) was applied to predict the protein-coding genes using default parameter and RAFTS3 (Vialle et al., 2016) for automatic annotation with best hits comparison with self-score cutoff 0.5 using our internal database of sequences of Fonsecaea ssp. and related species from the Chaetothyriales (Vicente et al., 2017; Moreno et al., 2018). F. pugnacius functional characteristics were determined with GO enrichment analyses at a significance level of ≤0.05 according to Ashburner et al. (2000) using the InterProScan5 (Quevillon et al., 2005) to access the protein domain families. Phylogenomic trees based on orthologous clusters were obtained and generated using ORTHOFinder (Emms and Kelly, 2018). A phylogenomic tree was inferred for each and all-orthogroup trees were resolved by the OrthoFinder duplication-loss coalescent model (Emms and Kelly, 2018) and the final tree was constructed using the STAG method, present in the OrthoFinder pipeline (Emms and Kelly, 2018). The enzymatous gene profile of F. pugnacius and other melanized fungi causing brain infection (C. bantiana, E. dermatitidis, F. monophora, Verruconis gallopava, and R. mackenziei) were predicted with CAZY (Cantarel et al., 2009) and MEROPS databases (Rawlings et al., 2016). In addition, numbers, classes and similarities were analyzed using an all-vs.-all similarity 40% ≥ search and e-value of 10–4 (Vicente et al., 2017).
Mitochondrial Genome Assembly and Annotation
The SPAdes v3.6.2 program (Bankevich et al., 2012) was used for assembly and mappingd the mitochondrial genome from the F. pugnacius sequencing reads previously aligned against the complete mtDNA of Fonsecaea pedrosoi CBS 271.37. Mitochondrial genome annotations were done based on Vicente et al. (2017) and Moreno et al. (2018) using SILA (Vialle, 2013) and the final figure was produced by the software package Circos (Krzywinski et al., 2009).
Virulence Assays in Animal Models
Tenebrio molitor Infection
A Tenebrio molitor larval model was used to evaluate the virulence potential based on Fornari et al. (2018) using parameters of survival and melanization after infection. Larvae were inoculated with 1 × 106 cells/mL in PBS solution above the legs and the ventral portion sterile PBS solution and SHAM without physical damage (no treatment) as negative control, using 10 larvae per group of inoculation, in triplicate. The larvae were kept in darkness at 37°C, mortality was monitored daily for 10 days and dead larvae were collected at 4, 24, 72, 168, and 240 h post infection omitting the pupae in the calculation (Scorzoni et al., 2013). Survival curves were plotted and statistical analyses were performed using the Log-rank (Mantel–Cox) test with Graph Pad Prism software and statistical differences were set at p < 0.05 according to Maekawa et al. (2015) and Vicente et al. (2017). Melanization was determined measuring the OD at 405 nm (Scorzoni et al., 2013; Perdoni et al., 2014).
Fungal burden and histological analysis of infected caterpillars was performed according to Fornari et al. (2018). The samples were homogenized in PBS solution with a TissueLyser (Qiagen, Hilden, Germany), inoculated on Mycosel agar at 30°C for 14 days and the number of colony forming units (CFUs) of fungal per mL of solution estimated with some colonies re-isolated and sequenced to confirm the species ID. Moreover, the caterpillar samples were embedded in Adracanth gum solution (Fornari et al., 2018), immersed in liquid nitrogen and sectioned by steel blades in a cryostat (Leica CM 1850, Wetzlar, Germany), stained with hematoxylin and eosin (HE) and observed with Axio Imager Z2 (Carl Zeiss, Jena, Germany) equipped with Metafer 4/VSlide automated capture software (Metasystems, Altlussheim, Germany).
Murine Infection
Fungal burden and cytokine production evaluation was performed using immunocompetent mice as model, according Bocca et al. (2006), Badali et al. (2011), Rodrigues et al. (2015), and Schneider et al. (2019). The animals selected were male Balb/c (6–8 weeks) and C57/BL6 mice, maintained under standard laboratory conditions with controlled temperature (23–25°C) with water and food ad libitum, according to recommendations of the Federal University of Paraná Ethics Committee (current approval certificate 1002).
The experiments were performed in triplicate using groups of six animals infected with F. pugnacius CBS 139214 and one negative control inoculated with sterile phosphate-buffered saline (PBS according to described by Fornari et al. (2018) and Schneider et al. (2019). The animals were infected intraperitoneally or intradermally (per hind footpad) with 100 μL of 1 × 106 propagules or sterile PBS and were monitored weekly and sacrificed at 7, 14, and 21 days post-infection using CO2 anesthesia in an appropriate chamber (Fornari et al., 2018). Brain, lung, liver, kidney, spleen, footpad, and blood were aseptically collected for analysis.
For fungal burden determination, samples tissues were weighed, homogenized and diluted in PBS for culture as described by Vicente et al. (2017). Results were expressed as number of CFU ± standard error of mean (SEM) per gram of fresh tissue, counting colonies from the seventh day until the 15th day. Cytokine production was measured from homogenized tissue obtained from infected and non-infected (healthy) animals by ELISA (Vicente et al., 2017). The cytokines interleukin-1β (IL-1β), TNF-α, interleukin-6 (IL-6) and monocyte chemoattractant protein-1 (MCP-1/Ccl2) were measured with kits purchased from eBioscience and used according the manufacturer’s instructions. Results were expressed as pg of cytokine ± standard error of mean (SEM) per 100 milligrams of tissue. The infected tissue samples were fixed in 10% formalin, dehydrated in alcohol, and embedded in paraffin (Fornari et al., 2018). Serial 5-μm sections were stained with hematoxylin and eosin to visualize pathogen morphology.
Results and Discussion
De novo Assembly
The genome sequencing of F. pugnacius CBS 139214 was performed using Illumina MiSeq and Ion proton producing 5,424,908 paired-end reads and 1,853,059 mate-paired reads, respectively. The final high-quality assembly comprised 386 contigs with 34,872,293 bp and 52% of G + C content. The genome size estimated is 34.8 Mb with average coverage of 48.75X and using 97% of the reads for draft assembling. The genome completeness, checked using BUSCO, revealed that the assembly had 98.8% completeness. A total of 6188 complete BUSCO genes were found, including 6176 being single-copy BUSCOs, of the 6265 BUSCO groups searched, 12 (0.2%) were duplicated BUSCO genes, 39 (0.6%) were fragmented BUSCOs and 38 (0.6%) represented missing BUSCOs. Sequencing data were submitted to GenBank (accession number WJFF00000000). In addition, 12,217 protein-coding genes and 35 tRNAs were predicted for F. pugnacius (Table 1). Expected values for Fonsecaea siblings were between 11,681 in F. nubica and 12,527 in F. pedrosoi (Vicente et al., 2017).
Table 1. Fonsecaea pugnacius genome data assembly and related species of Herpotrichiellaceae.
Comparing F. pugnacius genome size of 34.8 Mb with related species of Herpotrichiellaceae (Table 1), some differences were observed. The Cladophialophora immunda genome is nearly 8.15 Mb larger, while the Exophiala dermatitidis genome is 8.43 Mb smaller, being the smallest genome in the family with 26.37 Mb (Teixeira et al., 2017). Genome sizes within Fonsecaea were similar, i.e., between 33.39 and 35.23 Mb (Vicente et al., 2017). Judging from the phylogenomic tree (Figure 1), F. pugnacius clustered with the clinical representatives of Fonsecaea, in accordance with previous phylogenetic analyses of Fonsecaea and Cladophialophora based on internal transcribed spacer (ITS), partial beta-tubulin protein-coding gene (BT2) and cell division cycle 42 (CDC42) sequences performed by de Azevedo et al. (2015).
Figure 1. Phylogenomic tree based on the concatenated alignment of genomes of Herpotrichiellaceae family. Species names are given between phylogenomic tree and genome size bars. The genomes size is indicated and the colors are representing the clades of Chaetothyriales summarized in colored boxes.
The mitochondrial genome was assembled and 3 contigs were obtained with a total of 25,098 bp and a GC% of 25.49; the largest contig was 14,011 bp. There were 44 proteins in the mtDNA, i.e., 26 hypothetical proteins and 18 with known function (Figure 2). The proteins encoded by mitochondrial genomes of herpotrichiellaceous species were very similar, most being proteins involved in ATP synthesis and respiratory metabolism (Vicente et al., 2017).
Figure 2. Fonsecaea pugnacius CBS 139214 mtDNA. Rectangles represent annotated genes: in red the rRNAs; in green the tRNAs; in black the other genes and orange inner circle shows reads coverage.
Protein-Coding Gene Annotation and General Characteristics
A total of 12,217 protein-coding genes were identified in F. pugnacius, of which 11,124 were annotated as hypothetical proteins and 1,093 proteins had inferred functions (Supplementary Table 1). In the GO annotation, the proteins were separated into three large groups: biological process, cellular components and molecular functions (Figure 3 and Supplementary Table 1). Among proteins to which functions were attributed, various proteins were shared among the Fonsecaea siblings, such as proteins related to virulence in transporter families, proteins from the glyoxylate cycle, genes encoding proteins related to oxidative stress and involved in the detoxification of reactive oxygen species (ROS), cytochrome P450 monooxygenases (CYPs/P450s), heat shock proteins, proteins of melanin pathways, enzymes able to degrade aromatic carbon compounds, and others.
Figure 3. The most numerous gene families of Fonsecaea pugnacius based on Gene Ontology annotation separating proteins annotated into three large groups: biological process, cellular components, and molecular functions.
Carrier families of zinc, iron, manganese, and sugar have also been identified, i.e., MFS transporters and ABC transporters. In addition to basal metabolism, these carriers play a role in survival strategies. MFS is the largest family of transporters, ubiquitous to all living organisms and involved in the active excretion of antifungal drugs (Vela-Corcía et al., 2019). Reportedly, they also enhance antifungal resistance in Candida albicans, Aspergillus fumigatus, and Cryptococcus neoformans (Costa et al., 2014).
Sequences of proteins involved in the glyoxylate cycle were annotated: isocitrate lyase and malate synthase (Supplementary Table 1). The glyoxylate pathway is a metabolic strategy for the synthesis of carbohydrates from carbon compounds, such as acetate and other degradation products from ethanol, fatty acids, and poly-b-hydroxybutyrate (White et al., 2017). The glyoxylate cycle consists of a modification of the citric acid cycle (TCA), as it shares the same initial reactions of citrate and isocitrate generation but continues with the formation of succinate and glyoxylate (Dunn et al., 2009). This pathway has been associated with fungal virulence, since it allows energy production in environments where complex carbon compounds are poorly found (Lorenz and Fink, 2001) and has been described in infectious fungi such as Rhinocladiella mackenziei (Moreno et al., 2018), Fonsecaea siblings related to chromoblastomycosis (Vicente et al., 2017), Beauveria bassiana (Yang et al., 2016), Talaromyces marneffei (Thirach et al., 2007), and Candida albicans (Lorenz and Fink, 2001).
In Paracoccidioides brasiliensis, agent of another implantation mycosis, paracoccidioidomycosis, an increase of transcriptional levels of isocitrate lyase and malate synthase genes was reported in an infection model (Derengowski et al., 2008). This suggests a possible mechanism of adaptation of the fungus in response to the internal environment of the phagosome, which is poor in complex sources of carbon. The same function could be assigned to Fonsecaea siblings to explain fungal persistence inside macrophages, which, according to Queiroz-Telles et al. (2017), seems to be fungistatic (rather than fungicidal) against agents of chromoblastomycosis. In addition, it may be related to the ability of these fungi to survive on low carbon sources during brain infection.
Characterization of the partial genome of F. pugnacius revealed genes encoding proteins related to oxidative stress and involved in the detoxification of ROS, such as alternative oxidase, manganese superoxide dismutase and cytoplasmic thioredoxins. Earlier reports observed that the alternative oxidase enzyme is present in the internal mitochondrial membrane of plants and some fungi and protozoa, in an alternative route of oxidation of the electron transport chain in cellular respiration (Duvenage et al., 2018). In the fungus Neurospora crassa, levels of the nuclear gene transcripts AOX1, which encode the alternative oxidase, increase when the cytochrome C oxidase pathway is inhibited. These results indicate activation of an alternative pathway, which the organism applies to correct conditions of oxidative stress and to decrease the production of ROS in respiration when competing with electrons of the classical oxidation pathway. Some studies have revealed that AOX1 gene expression can be stimulated under stress conditions, such as low temperature and ROS low level, AOX1 functioning as an antioxidant (Missall et al., 2004). Some pathogens resist to oxidative stress in the hostile environment of the phagosomes by the production of antioxidant enzymes that detoxify ROS, such as alternative oxidase, catalase and superoxide dismutase (Duvenage et al., 2018). ROS is an important cellular detrimental agent associated with the activation of immune response in human cells infected with fungi causing dermatomycoses or invasive mycoses (Castro et al., 2017).
According Vicente et al. (2017), many of the CYP/P450 enzymes revealed in herpotrichiellaceous fungi are abundantly present in Fonsecaea siblings. Likewise, they were observed in F. pugnacius. Cytochrome P450 monooxygenases are heme-thiolate proteins with roles in oxidative functions, e.g., degradation of xenobiotic compounds (Jawallapersand et al., 2014). Teixeira et al. (2017) noted that some black fungi are among the species of Ascomycota with the highest numbers of CYPs, with family expansion and diversification through gene duplication which might explain opportunism. The authors noted that these enzymes are involved in the metabolism of phenolic compounds and aromatic hydrocarbons, and suggested that similar compounds present in the human brain might explain their neurotropic predilection. Studies showed that fungi belonging to the genera Fusarium, Penicillium, Aspergillus, and the family Herpotrichellaceae are capable of degrading aromatic compounds (Conceição et al., 2005; Satow, 2005; Teixeira et al., 2017). Several species of black fungi have been isolated from hydrocarbon-contaminated environments (de Hoog et al., 2004). Fungi that are capable of assimilating monoaromatic hydrocarbons are enriched in the domestic environment (Sterflinger and Prillinger, 2001; Woertz et al., 2001; Prenafeta-Boldú et al., 2002). Prenafeta-Boldú et al. (2006) suggested physiological links between hydrocarbon assimilation by black fungi and certain patterns of brain infection. The brain contains small molecules that resemble alkylbenzene, phenylalanine metabolic products and lignin biodegradation intermediates, having structural similarity to neurotransmitters such as dopamine, which is catabolized in the brain (Fernstrom and Fernstrom, 2007). Tyrosine, used for the biosynthesis of the neurotransmitters dopamine, noradrenaline and adrenaline, has phenylalanine as precursor and it is involved in the formation of melanin and neuromelanin, dark pigments synthesized from L-dopamine for brain protection (Teixeira et al., 2017; Moreno et al., 2018).
Heat shock proteins are considered virulence factors because of their roles in thermotolerance and as molecular chaperones and are found in all prokaryotes and eukaryotes. They are classified based on approximate molecular weight (Tiwari et al., 2015). In F. pugnacius, the following families were identified: Hsp7, Hsp60, Hsp70, Hsp80, and Hsp90. Factors triggering the synthesis of heat shock proteins are oxidative or nutritional stress, UV radiation and exposure to chemical substances, indicating a protective role and aiding in cellular adaptation (Pockley, 2001).
Melanins confer resistance to heat, cold, enzymatic action and organic solvents, function as antioxidants and increase antifungal resistance (Nosanchuk and Casadevall, 2003). Main production route in Fonsecaea is the DHN pathway from acetate (acetyl-CoA) derived from glucose metabolism (Cunha et al., 2005; Casadevall and Eisenman, 2012). Melanin is a recognized virulence factor in several black and white pathogenic and opportunistic fungi, such as Candida albicans, Cryptococcus neoformans, Aspergillus fumigatus, Exophiala dermatitidis, Paracoccidioides brasiliensis, Histoplasma capsulatum, and Sporothrix schenckii (Jacobson, 2000; Langfelder et al., 2003; Morris-Jones et al., 2003; Nosanchuk and Casadevall, 2003). Homologous proteins were identified in the herpotrichiellaceous black fungi Exophiala dermatitidis (Youngchim et al., 2004). In F. pugnacius, melanin-associated proteins related to DHN and DOPA pathways were observed (Supplementary Table 1), as reported previously in F. monophora (Liu et al., 2019), such as tyrosinase, homogentisate dioxygenase and scytalone. Histopathological studies of organs such as brain, lung, liver and spleen in animal models did not clarify how the melanin production pathways are blocked in albino mutants of E. dermatitidis (Sudhadham et al., 2008), but the ability to block the oxidative burst increases significantly the pathogenic potential of E. dermatitidis, since the host is unable to eliminate this fungus (Kumar et al., 2019).
Virulence in Animal Models
The chromoblastomycosis agent F. pugnacius was described causing a secondary disseminated infection in an apparently immunocompetent patient. Among all agents of chromoblastomycosis, this is a unique type of infection, starting with a chronic skin disease and finally evolving to cerebritis. de Azevedo et al. (2015) reported that F. pugnacius was able to produce muriform cells in skin, but hyphae were present in the brain. Dissemination from skin to brain apparently led to conversion to another type of invasive morphology which has not been observed in agents of chromoblastomycosis. In F. monophora, which has also been reported from brain infection, the infection route was probably by inhalation, as no skin involvement was observed in any of the patients. Brain infection with a subcutaneous origin has thus far only been observed in F. pugnacius (de Azevedo et al., 2015).
Tenebrio molitor larvae were infected with inoculum concentrations of 5 × 106 cells/mL and observed for 10 days. The larvae infected with F. pugnacius exhibited higher mortality rates than control groups, PBS and SHAM (Figure 4A). Fonsecaea pugnacius presented a lower mortality rate than F. monophora, F. erecta, and F. pedrosoi, as reported by Fornari et al. (2018). This indicates that F. pugnacius infection presents a slower development compared to Fonsecaea siblings involved in chromoblastomycosis, as well as to environmental saprobes. The fungal burden inside the larvae was assessed and presented significant numbers of CFUs, despite the low mortality rates caused by F. pugnacius: CFU values were initially countless and decreased along 72, 164, and 240 h post-infection (Figure 4B).
Figure 4. Virulence test of Fonsecaea pugnacius using animal models: Tenebrio molitor and Balb/c (Mus musculus). (A) Mortality of the Tenebrio molitor larvaes at intervals of 4, 24, 72, 168, and 240 h post-infection. (B) Fungal burden of Tenebrio molitor larvae tissues at intervals of 4, 24, 72, 168, and 240 h post-infection. (C) Melanization of Tenebrio molitor hemolymph demonstrated by measuring the OD405 nm at 4, 24, 72, 168, and 240 h post-infection. (D) Presence of melanized hyphae in Tenebrio molitor larvae tissue 72 h post-infection. (E) Histopathology of footpad Balb/c tissue with presence of muriform cell after 21 days of intradermal infection. (F) Necrotic lesion in footpad of Balb/c mice after 7 days of intradermal inoculation.
Melanization of the larvae is an intracellular defense response and an effective barrier to infection. After 24 h, the larvae infected with F. pugnacius presented a dark pigmentation, caused by melanization in the hemolymph, which was observed during the entire 10-day period of analysis by visual observation and spectrophotometry (Figure 4C). Similar results were obtained by Fornari et al. (2018) in Fonsecaea siblings, showing maximum melanization with 24 h post infection. The capacity of the fungus to survive inside the larvae was confirmed by histopathology, revealing melanized hyphae in tissue that had developed within 4–72 h (Figure 4D).
Virulence tests using murine models Balb/c and C57/BL6 were conducted using two infection pathways: intradermal (per hind footpad) and intraperitoneal. In view of determination of fungal burden, F. pugnacius was recovered from kidney, lung and liver after 7 days of incubation, indicating a certain preference of the fungus for these organs. After 14 days of intraperitoneal inoculation, F. pugnacius was recovered from the brain. At 21 days after infection, a muriform cell was observed in histopathology of the footpad (Figure 4E). The animal host infected intraperitoneally presented 1 × 102 and 1 × 104 CFU/mL in blood and organs (lung, kidney and spleen) after 7 and 14 days, respectively. The animals infected intradermally presented 2 × 106 CFU/g in the plantar cushion. The clinical aspects of these animals were evaluated, but no lesions, tissue necrosis or morphological alterations of internal organs were observed, except for plantar cushion swellings with (sub)cutaneous lesions (Figure 4F). Vicente et al. (2017) obtained similar results with Balb/c mice infected with F. pedrosoi by intradermal inoculation.
Immunological assays revealed insignificant levels of IL-2, INF-γ, TNF-α, IL-6 and IL-10, compared to what was described for other chromoblastomycosis agents including F. pedrosoi (Dong et al., 2018) and F. monophora (Jiang et al., 2018). The low immune response could be associated with a low virulent ability of F. pugnacius. Although the fungus carries several genes with roles in pathogenicity and ability to survive in murine tissue, the immune system may not recognize their protein products, judging from absence of a cytokine response.
Potential Virulence Related to Enzymatous Genes
The observed brain infection might have been enhanced by the fungus ability to metabolize monoaromatic substrates as carbon source. Among these compounds, vanillic acid, phenyl acid, L-tyrosine, L-dopa, L-phenylanine, dopamine, sphingosine, and others are present in mammalian brain (Prenafeta-Boldú et al., 2006). Subsequently, we verified whether an enzymatic apparatus competent to degrade these compounds is present in the F. pugnacius genome, i.e., with carboxylases, reductases, aldolases, and kinases. The enzymatous gene profile was established with CAZymes and MEROPS databases. Carbohydrate-active enzyme analysis resulted in 476 genes encoding putative CAZymes, comprising 93 auxiliary activities (AA), 5 carbohydrate binding modules (CBM), 145 carbohydrate esterases (CE), 121 glycoside hydrolases (GH), 112 glycosyl transferases (GT), and zero polysaccharide lyases (PL). For MEROPS, 21 aspartic peptidases (A) were annotated, 59 cysteine peptidases (C), 14 trypsin peptidases (I), 133 metallo-peptidases (M), zero glutamic peptidases (G), zero asparagines peptidases (N), zero mixed peptidases (P), 189 serine peptidases (S) and 21 threonine peptidases (T) (Figure 5).
Figure 5. Species enzymatous gene profiles and epidemiology of human host pathogenicity are summarized in colored boxes. Cladophialophora bantiana: beige; Exophiala dermatitidis: light green; Fonsecaea pugnacius: medium green; Fonsecaea monophora: light blue; Rhinocladiella mackenziei: medium blue; Verruconis gallopava: dark blue. (A) Carbohydrate-active enzyme (CAZY) annotation separated by classes. (B) MEROPS annotation separated by classes. (C) The species of black fungi related to subcutaneous and cerebral human infection and others affected organs: lungs, lymphatic system, kidneys, spleen, liver, and thyroid gland.
We performed a comparative analysis of F. pugnacius with closely related neurotropic species: Cladophialophora bantiana, Exophiala dermatitidis, Fonsecaea monophora, Rhinocladiella mackenziei, and Verruconis gallopava (Figure 5 and Supplementary Table 2), which have been reported as agents of cutaneous and deep infection in skin, lungs, lymphatic system, kidneys, spleen, liver and thyroid glands (Cardeau-Desangles et al., 2013; de Azevedo et al., 2015; Geltner et al., 2015; Jennings et al., 2017; Mukhopadhyay et al., 2017; Stokes et al., 2017; Grewal et al., 2018; Sideris and Ge, 2018; Klasinc et al., 2019; Miossec et al., 2019; Wang et al., 2019; Alabdely and Alzahrani, 2020). The CAZy database is an exceptional reporter of fungal lifestyles degrading complex polysaccharides (Cantarel et al., 2009; Vicente et al., 2017). The GH class contains the majority of catalytic enzymes to degrade lignocelluloses (Cantarel et al., 2009) and seems to be expanded in species associated with disseminated infection (Figure 5). Through this analysis (Supplementary Table 3), we observed a high number of GH enzymes in all genomes analyzed, except for GH32 which was present in low numbers (C. bantiana: 3; F. pugnacius: 2; E. dermatitidis, F. monophora, V. gallopava: 1, R. mackenziei: 0). The CAZyme family GH32 are enzymes associated with energy storage and seem to provide energy for survival in extreme environments (Oren, 2011). Additionally, the GH18 family, observed in the species under study, are chitinases. β-1,3-Glucanases are able to degrade chitin that is normally present in animal exoskeletons. Likewise, 29 GH families involved in the degradation of plant biomass (Li et al., 2011) were observed in the strains evaluated, explaining the dual ability of these fungi to invade plant and animal tissues.
The high number glycosyltransferases (GTs) can also be linked to neurotropism, since these enzymes are responsible for the biosynthesis of glycoside (Cantarel et al., 2009). Some fungi are able to convert phenolic compounds into their corresponding glycosides (Tronina et al., 2013) and this ability may be derived from xenobiotic metabolism (Yang et al., 2014). The vertebrate brain contains a wide variety of gangliosides that are localized in specific cell types, such as on the surface of plasmatic membranes (Stanley, 2016). Gangliosides contain a fatty acid and a sphingosine base and are involved in neural functions such as memory formation, synaptic transmission, regeneration (Zitman et al., 2010). In cases of cerebral infection, the first symptoms described are locomotor difficulties and severe headache (Koo et al., 2010; de Azevedo et al., 2015). The GT41 family was the most numerous in our annotation, it is related to the metabolism of serine-threonine as indicated and involved in the process of subcutaneous infection (Vicente et al., 2017). Also, the GT2 family was abundant; proteins that act as chitin synthases (Breton et al., 2006) are numerous in many fungi (Lee et al., 2018; Stone et al., 2018).
The AA CAZyme family is composed of lignolytic enzymes and are commonly found in plant pathogens. Similar to the GH and CMB families, they have a role in breaking down plant cell-wall polysaccharides (Lowe et al., 2015). The CBM family is composed of lectins and sugar transporters, while GH are glycosyltransferases such as lignocelluloses or chitinases (Yang et al., 2014). However, the GH class was more numerous in the species related to disseminated and cerebral infection (Figure 5). The lectins from yeasts and fungi have been associated with early stages of human infection (Varrot et al., 2013), whereas in bacteria they seem to be involved in recognition of host glycans (Imberty and Varrot, 2008).
Polysaccharide lyases are enzymes able to cleave polysaccharide chains. This group presents many fold types (or classes), indicating that PLs are polyphyletic (Lombard et al., 2010). Many fungi that do not have enzymes from the PL CE and GH families are saprobes, as these classes have enzymes related to cell wall degradation in plants (Zhao et al., 2013). This enzymatous gene profile suggested a dual ecological ability of these agents, in line with their extremotolerance and adaptability to variable environmental niches, which is a prerequisite for their opportunistic profiles.
Peptidases play key roles in penetration of microorganisms into host tissue and are involved in pathogen-host interactions (Ohm et al., 2012). Herpotrichiellaceous agents produce a variety of extracellular peptidases for the degradation of environmental substrates, indicating a poorly specialized nutritional strategy (Sriranganadane et al., 2010; Vicente et al., 2017). The MEROPS S (serine) and M (metallo) peptidase families were the highest in number in the analyzed species causing brain infection. These two groups have been reported to be significantly enriched in transcriptome analyses of E. dermatitidis during infection (Poyntner et al., 2018). Among the Serine families, classes S33 and S9 were more numerous (Supplementary Table 3), both involved in prolyl metabolism. Class S33 is a prolyl aminopeptidase family, which is not essential for growth but may confer a selective advantage allowing the organism to use proline-rich substrates (Iqbal et al., 2018). Teixeira et al. (2017) observed an expansion of the protein-degrading peptidase enzyme family M38 (isoaspartyl dipeptidases) in the bantiana-clade, the most numerous of metallo-peptidases family in the analyzed strains. This family is unusual in that the majority of characterized proteins are not peptidases but are associated with β-aspartic dipeptidase acting in the release of iso-aspartate residues from peptides (Palmeira et al., 2018).
Differences between closely related taxa are expected in ecology-related genes. The enzymatic repertoire of these fungi shows their ability to degrade a wide variety of substrates (Figure 5). This may suggest generalist and opportunistic ecology comparable to Aspergillus spp. permitting transfer from the environment to the animal host (Vicente et al., 2017), rather than pathogenicity where focused adaptation (Moran et al., 2011). In addition, we analyzed genes codifying proteins related to degradation aromatic carbons pathway (Teixeira et al., 2017; Moreno et al., 2018) and the strains studied present a range of genes encoding homologous proteins (Table 2). The virulence of these strains is partially explained by general factors like the presence of melanin in the cell wall, thermotolerance and the ability to assimilate of monoaromatic hydrocarbons (Moreno et al., 2018). Furthermore, CYPs are involved in the degradation of aromatic hydrocarbons and xenobiotic metabolism, developing functions in the fungal pathogenicity and in the detoxification of exogenous compounds (Moreno et al., 2018).
Table 2. Proteins related to degradation aromatic carbons pathway in Fonsecaea pugnacius and homologs in black fungi related to subcutaneous and brain infection.
The family Herpotrichiellaceae contains numerous black fungi that present tolerance to various types of stress, showing great adaptability to extreme environmental conditions, presumably resulting from genomic information. Genomic studies of F. pugnacius showed the wide variety of genes involved in extreme tolerance and enzymes associated to occurrence of virulence factors. The survival capacity of fungi in animal models was confirmed by histopathological analysis and the presence of melanin in the host tissue. We have shown that F. pugnacius can colonize the brain and cause subcutaneous lesions with the formation of muriform cells in the murine model. An ecological capacity can be concluded from the presence of metabolic pathways for extremetolerance combined with the ability to infect human hosts. However, complementary molecular studies must be carried out in order to strengthen the connections between ecology and clinical profiles.
Gostinčar et al. (2018) found a link between (poly-)extremotolerance and opportunism, the ability to metabolize monoaromatic hydrocarbons improving human or animal disease. The order Chaetothyriales comprises opportunists, in which tolerance to various types of stress is associated with adaptability, presumably resulting in a large potential for habitat changes. The infection is as a side effect of the fungus adaptation to the human host, demonstrating that is not a favorable habitat, nor relevant to their evolutionary process. This defines opportunism against pathogenicity, where the infection is advantageous for the fitness of the species. Most organisms considered opportunistic are unable to transmit from host to host, so specific adaptations will be lost with the cure of the infection, explaining the lack of complex virulence characteristics. Therefore, opportunistic infections can be considered an evolutionary dead end, which is unlikely to lead to true pathogenicity (Gostinčar et al., 2018).
Data Availability Statement
The datasets generated for this study can be found in the GenBank/accession number: WJFF00000000/waiting to publish for to be public date.
Ethics Statement
All animal experiments in this study were approved by the Federal University of Paraná Ethics Committee (approval certificate 1002) and performed according to the Committee’s recommendations.
Author Contributions
AB, VV, RG, and SH contributed to the conception and design of the study. FC, AL, LM, and NS organized the database. AB, GS, BS, FM, BL, RC, VB, EB, VP, and NH performed the analysis. AB wrote the first draft of the manuscript. GS wrote sections of the manuscript. All the authors contributed to the manuscript revision, read and approved the submitted version.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer AM declared a past co-authorship with one of the authors SH to the handling editor.
Funding
This work was supported by the Brazilian Federal Agency for Support and Evaluation of Graduate: Education Coordination for the Improvement of Higher Education Personnel—CAPES (PVE project grant number 0592012), the PRINT (www.capes.gov.br), and the National Council for Scientific and Technological Development (http://cnpq.br/), Brazil; and by The National Institute of Science and Technology of Biological Nitrogen Fixation/CNPq/MCT (grant number 573828/2008-3). VV received fellowships from the National Council for Scientific and Technological Development (CNPq), Brasilia, Brazil.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.00822/full#supplementary-material
Footnotes
References
Alabdely, M. H., and Alzahrani, H. A. (2020). Inferior vena cava filter eroding the aorta. J. Appl. Hematol. 11:37. doi: 10.4103/joah.joah_71_19
Arcobello, J. T., and Revankar, S. G. (2020). Phaeohyphomycosis. Semin. Respir. Crit. Care Med. 41, 131–140. doi: 10.1055/s-0039-3400957
Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene Ontology: tool for the unification of biology. Nat. Genet. 25, 25–29. doi: 10.1038/75556
Badali, H., Prenafeta-Boldu, F. X., Guarro, J., Klaassen, C. H., Meis, J. F., and de Hoog, G. S. (2011). Cladophialophora psammophila, a novel species of Chaetothyriales with a potential use in the bioremediation of volatile aromatic hydrocarbons. Fungal Biol. 115, 1019–1029. doi: 10.1016/j.funbio.2011.04.005
Bankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., Dvorkin, M., Kulikov, A. S., et al. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477. doi: 10.1089/cmb.2012.0021
Besemer, J., Lomsadze, A., and Borodovsky, M. (2001). GeneMarkS: a self-training method for prediction of gene starts in microbial genomes. Implications for finding sequence motifs in regulatory regions. Nucleic Acids Res. 29, 2607–2618. doi: 10.1093/nar/29.12.2607
Bocca, A. L., Brito, P. P. M. S., Figueiredo, F., and Tosta, C. E. (2006). Inhibition of nitric oxide production by macrophages in chromoblastomycosis: a role for Fonsecaea pedrosoi melanin. Mycopathologia 161, 195–203. doi: 10.1007/s11046-005-0228-6
Breton, C., Šnajdrová, L., Jeanneau, C., Koča, J., and Imberty, A. (2006). Struc- tures and mechanisms of glycosyltransferases. Glycobiology 16, 29R–37R.
Cañete-Gibas, C. F., and Wiederhold, N. P. (2018). The black yeasts: an update on species identification and diagnosis. Curr. Fungal Infect. Rep. 12, 59–65. doi: 10.1007/s12281-018-0314-0
Cantarel, B. L., Coutinho, P. M., Rancurel, C., Bernard, T., Lombard, V., and Henrissat, B. (2009). The Carbohydrate-Active EnZymes database (CAZy): an expert resource for glycogenomics. Nucleic Acids Res. 37, 233–238. doi: 10.1093/nar/gkn663
Cardeau-Desangles, I., Fabre, A., Cointault, O., Guitard, J., Esposito, L., Iriart, X., et al. (2013). Disseminated Ochroconis gallopava infection in a heart transplant patient. Transplant Infect. Dis. 15, E115–E118. doi: 10.1111/tid.12084
Casadevall, A., and Eisenman, E. C. (2012). “Synthesis and assembly of fungal melanin,” in Applied Microbiology and Biotechnology, Vol. 93, ed. A. Steinbüchel (New York, NY: Springer Press), 931–940. doi: 10.1007/s00253-011-3777-2
Castro, R. J. A., Siqueira, M. I, Jerônimo, M. S., Basso, A. M. M., Veloso, P. H. Jr., Magalhães, K. G., et al. (2017). The major chromoblastomycosis etiologic agent Fonsecaea pedrosoi activates the NLRP3 inflammasome. Front. Immunol. 8:1572. doi: 10.3389/fimmu.2017.01572
Conceição, D. M., Angelis, D. A., Bidoia, E. D., and Angelis, D. F. (2005). Fungos filamentosos isolados do rio Atibaia, SP, e refinarias de petróleo biodegradadores de compostos fenólicos. Rev. Inst. Biol. 72, 99–106.
Costa, C., Dias, P. J., Sá-Correia, I., and Teixeira, M. C. (2014). MFS multidrug transporters in pathogenic fungi: do they have real clinical impact? Front. Physiol. 5:197. doi: 10.3389/fphys.2014.00197
Cunha, M. M., Franzen, A. J., Alviano, D. S., Zanardi, E., Alviano, C. S., Souza, W., et al. (2005). Inhibition of melanin synthesis pathway by tricyclazole increases susceptibility of Fonsecaea pedrosoi against mouse macrophages. Microsc. Res. Tech. 68, 377–384. doi: 10.002/jemt.20260
de Azevedo, C. M. P. S., Gomes, R. R., Vicente, V. A., Santos, D. W. C. L., Marques, S. G., do Nascimento, M. M. F., et al. (2015). Fonsecaea pugnacius, a novel agente of disseminated chromoblastomycosis. J. Clin. Microbiol. 53, 2674–2685. doi: 10.1128/JCM.00637-15
de Hoog, G. S., Attili-Angelis, D., Vicente, V. A., Van Den Ende, A. H., and Queiroz-Telles, F. (2004). Molecular ecology and pathogenic potential of Fonsecaea species. Med. Mycol. 42, 405–416. doi: 10.1080/13693780410001661464
Derengowski, L. S., Tavares, A. H., Silva, S., Procópio, L. S., Felipe, M. S., and Silva-Pereira, I. (2008). Upregulation of glyoxylate cycle genes upon Paracoccidioides brasiliensis internalization by murine macrophages and in vitro nutritional stress condition. Med. Mycol. 46, 125–134. doi: 10.1080/13693780701670509
Dong, B., Tong, Z., Li, R., Chen, S. C. A., Liu, W., Liu, W., et al. (2018). Transformation of Fonsecaea pedrosoi into sclerotic cells links to the refractoriness of experimental chromoblastomycosis in BALB/c mice via a mechanism involving a chitin-induced impairment of IFN-γ production. PLoS Negl. Trop. Dis. 12:e0006237. doi: 10.1371/journal.pntd.0006237
Doymaz, M. Z., Seyithanoglu, M. F., Hakyemez, I., Gultepe, B. S., Cevik, S., and Aslan, T. (2015). A case of cerebral phaeohyphomycosis caused by Fonsecaea monophora, a neurotropic dematiaceous fungus, and a review of the literature. Mycoses 58, 187–192. doi: 10.1111/myc.12290
Dunn, M. F., Ramirez-Trujillo, J. A., and Hernández-Lucas, I. (2009). Major roles of isocitrate lyase and malate synthase in bacterial and fungal pathogenesis. Microbiology 155, 3166–3175. doi: 10.1099/mic.0.030858-0
Duvenage, L., Walker, L. A., Bojarczuk, A., Johnston, A. S., MacCallum, D. M., Munro, C. A., et al. (2018). Alternative oxidase induction protects Candida albicans from respiratory stress and promotes hyphal growth. bioRxiv [Preprint]. doi: 10.1101/405670
Emms, D. M., and Kelly, S. (2018). OrthoFinder2: Fast and accurate phylogenomic orthology analysis from gene sequences. bioRxiv [Preprint]. doi: 10.1101/466201
Fernstrom, J., and Fernstrom, M. (2007). Tyrosine, phenylalanine, and catecholamine synthesis and function in the brain. J. Nutr. 137, 1539–1548. doi: 10.1093/jn/137.6.1539S
Fornari, G., Gomes, R. R., Degenhardt-Goldbach, J., Dos Santos, S. S., De Almeida, S. R., Dos Santos, G. D., et al. (2018). A model for trans-kingdom pathogenicity in Fonsecaea agents of human chromoblastomycosis. Front. Microbiol. 9:2211. doi: 10.3389/fmicb.2018.02211
Geltner, C., Sorschag, S., Willinger, B., Jaritz, T., Saric, Z., and Lass-Florl, C. (2015). Necrotizing mycosis due to Verruconis gallopava in an immunocompetent patient. Infection 43, 743–746. doi: 10.1007/s15010-015-0757-y
Gostinčar, C., Zajc, J., Lenassi, M., Plemenitas, A., de Hoog, G. S., Al-Hatmi, A. M. S., et al. (2018). Fungi between extremotolerance and opportunistic pathogenicity on humans. Fungal Divers. 93, 195–213. doi: 10.1007/s13225-018-0414-8
Grewal, H. K., Kumar, P. N., Shah, N., and Timpone, J. G. (2018). Verruconis Gallopava: recurrent infection in a renal transplant recipient. Infect. Dis. Clin. Pract. 26, e12–e15. doi: 10.1097/IPC.0000000000000569
Helbig, S., Thuermer, A., Dengl, M., Krukowski, P., and de With, K. (2018). Cerebral abscess by Fonsecaea monophora —The first case reported in Germany. Open Forum Infect. Dis. 5:ofy129. doi: 10.1093/ofid/ofy129.2018
Imberty, A., and Varrot, A. (2008). Microbial recognition of human cell surface glycoconjugates. Curr. Opin. Struct. Biol. 18, 567–576. doi: 10.1016/j.sbi.2008.08.001
Iqbal, M., Dubey, M., Gudmundssn, M., Viketoft, M., Jesen, D. F., Karlsson, M., et al. (2018). Comparative evolutionary histories of fungal proteases reveal gene gains in the mycoparasitic and nematode-parasitic fungus Clonostachys rósea. BMC Evol. Biol. 18:171. doi: 10.1186/s12862-018-1291-1
Jacobson, E. S. (2000). Pathogenic roles for fungal melanins. Clin. Microbiol. Rev. 13, 708–717. doi: 10.1128/CMR.13.4.708-717.2000
Jawallapersand, P., Mashele, S. S., Kovačič, L., Stojan, J., Komel, R., Pakala, S. B., et al. (2014). Cytochrome P450 Monooxygenase CYP53 Family in fungi: comparative structural and evolutionary analysis and its role as a common alternative anti-fungal drug target. PLoS One 9:107209. doi: 10.1371/journal.pone.0107209
Jennings, Z., Kable, K., Halliday, C. L., Nankivell, B. J., Kok, J., Wong, G., et al. (2017). Verruconis gallopava cardiac and endovascular infection with dissemination after renal transplantation: case report and lessons learned. Med. Mycol. Case Rep. 15, 5–8. doi: 10.1016/j.mmcr.2016.12.006
Jiang, M., Cai, W., Zhang, J., Xie, T., Xi, L., Li, X., et al. (2018). Melanization of a meristematic mutant of Fonsecaea monophora increase the pathogenesis in a BALB/c mice infection model. Med. Mycol. 56, 979–986. doi: 10.1093/mmy/myx148
Klasinc, R., Riesenhuber, M., Bacher, A., and Willinger, B. (2019). Invasive fungal infection caused by Exophiala dermatitidis in a patient after lung transplantation: case report and literature Review. Mycopathologia 184, 107–113. doi: 10.1007/s11046-018-0275-4
Koo, S., Klompas, M., and Marty, F. M. (2010). Fonsecaea monophora cerebral phaeohyphomycosis: case report of successful surgical excision and voriconazole treatment and review. Med. Mycol. 48, 769–774. doi: 10.3109/13693780903471081
Krzywinski, M., Schein, J., Birol, I., Connors, J., Gascoyne, R., Horsman, D., et al. (2009). Circos: an information aesthetic for comparative genomics. Genome Res. 19, 1639–1645. doi: 10.1101/gr.092759.109
Kumar, D., Anjum, N., Kumar, A., and Das, S. (2019). Role of Melanin production in fungal pathogenesis. Int. J. Sci. Res. 8, 20–22.
Langfelder, K., Streibel, M., Jahn, B., Haase, G. E., and Brakhage, A. A. (2003). Biosynthesis of fungal melanins and their importance for human pathogenic fungi. Fungal Genet. Biol. 38, 143–158. doi: 10.1016/S1087-1845(02)00526-1
Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with bowtie 2. Nat. Methods 9, 357–359. doi: 10.1038/nmeth.1923
Lee, C.-S., Kong, W.-S., and Park, Y.-J. (2018). Genome Sequencing and Genome-Wide Identification of Carbohydrate-Active Enzymes (CAZymes) in the White Rot Fungus Flammulina fennae. Biotechnol. Lett. 46, 300–312. doi: 10.4014/mbl.1808.08012
Li, D. C., Li, A. N., and Papageorgiou, A. C. (2011). Cellulases from thermophilic fungi: recent insights and biotechnological potential. Enzyme Res. 2011:308730. doi: 10.4061/2011/308730
Liu, Y., Huang, X., Liu, H., Xi, L., and Cooper, C. R. (2019). Increased virulence of albino mutant of Fonsecaea monophora in Galleria mellonella. Med. Mycol. 57, 1018–1023. doi: 10.1093/mmy/myz007
Lombard, V., Bernard, T., Rancurel, C., Brumer, H., Coutinho, P. M., and Henrissat, B. (2010). A hierarchical classification of polysaccharide lyases for glycogenomics. Biochem. J. 432, 437–444. doi: 10.1042/BJ20101185
Lorenz, M., and Fink, G. R. (2001). The glyoxylate cycle is required for fungal virulence. Nature 412, 83–86. doi: 10.1038/35083594
Lowe, R. G. T., McCorkelle, O., Bleackley, M., Collins, C., Faou, P., Mathivanan, S., et al. (2015). Extracellular peptidases of the cereal pathogen Fusarium graminearum. Front. Plant Sci. 6:962. doi: 10.3389/fpls.2015.00962
Lucasse, C., Chardome, J., and Magis, P. (1954). Cerebral mycosis from Cladosporium tricho.ides in a native of the Belgian Congo. Ann. Soc. Belge Med. Trop. 34, 475–478.
Maekawa, L. E., Rossoni, R. D., Barbosa, J. O., Jorge, A. O. C., Junqueira, J. C., and Valera, M. C. (2015). Different extracts of Zingiber officinale decrease Enterococcus faecalis infection in Galleria mellonella. Braz. Dent. J. 26, 105–109. doi: 10.1590/0103-6440201300199
Miossec, C., Jacob, S., Peipoch, L., Brard, M., Jolivet, E., Hochedez, P., et al. (2019). Cerebral phaeohyphomycosis due to Cladophialophora bantiana in a French Guianese child. J. Mycol. Méd. 30:100918. doi: 10.1016/j.mycmed.2019.100918
Missall, T. A., Lodge, J. K., and Mcewen, J. E. (2004). Mechanisms of resistance to oxidative and nitrosative stress: implications for fungal survival in mammalian hosts. Eukaryot. Cell 3, 835–846. doi: 10.1128/ec.3.4.835-846.2004
Moran, G. P., Coleman, D. C., and Sullivan, D. J. (2011). Comparative genomics and the evolution of pathogenicity in human pathogenic fungi. Eukaryot. Cell 10, 34–42. doi: 10.1128/EC.00242-10
Moreno, L. F., Ahmed, A. A. O., Brankovics, B., Cuomo, C. A., Menken, S. B. J., Taj-Aldeen, S. J., et al. (2018). Genomic understanding of an infectious brain disease from the desert. G3 8, 909–922. doi: 10.1534/g3.117.300421
Morris-Jones, R., Youngchim, S., Gomez, B. L., Aisen, P., Hay, R. J., Nosanchuk, J. D., et al. (2003). Synthesis of melanin-like pigments by Sporothrix schenckii in vitro and during mammalian infection. Infect. Immun. 71, 4026–4033. doi: 10.1128/iai.71.7.4026-4033.2003
Mukhopadhyay, S. L., Mahadevan, A., Bahubali, V. H., Dawn Bharath, R., Prabhuraj, A. R., Maji, S., et al. (2017). A rare case of multiple brain abscess and probably disseminated phaeohyphomycosis due to Cladophialophora bantiana in an immunosuppressed individual from India. J. Mycol. Méd. 27, 391–395. doi: 10.1016/j.mycmed.2017.04.002
Najafzadeh, M. J., Gueidan, C., Badali, H., van den Ende, A. H. G., Xi, L., and de Hoog, G. S. (2009). Genetic diversity and species delimitation in the opportunistic genus Fonsecaea. Med. Mycol. 47, 17–25. doi: 10.1080/13693780802527178
Najafzadeh, M. J., Vicente, V. A., Sun, J., Meis, J. F., and de Hoog, G. S. (2011). Fonsecaea multimorphosa sp. nov., a new species of Chaetothyriales isolated from a feline cerebral abscess. Fungal Biol. 115, 1066–1076. doi: 10.1016/j.funbio.2011.06.007
Nobrega, J. P. S., Rosemberg, S., Adami, A. M., Heins-Vaccari, E. M., Lacaz, C. S., and Brito, T. (2003). Feohifomicose cerebral (“cromoblastomicose”) por Fonsecaea pedrosoi: primeiro caso demonstrado por cultura do fungo no Brasil. Rev. Inst. Med. Trop. S. Paulo 45, 217–220. doi: 10.1590/S0036-46652003000400008
Nosanchuk, J. D., and Casadevall, A. (2003). The contribution of melanin to microbial pathogenesis. Cell. Microbiol. 5, 203–223. doi: 10.1046/j.1462-5814.2003.00268.x
Ohm, R. A., Feau, N., Henrissat, B., Schoch, C. L., Horwitz, B. A., Barry, K. W., et al. (2012). Diverse lifestyles and strategies of plant pathogenesis encoded in the genomes of eighteen Dothideomycetes fungi. PLoS Pathog. 8:e1003037. doi: 10.1371/journal.ppat.1003037
Oren, A. (2011). Thermodynamic limits to microbial life at high salt concentrations. Environ. Microbiol. 13, 1908–1923. doi: 10.1111/j.1462-2920.2010.02365.x
Ozgun, H. D., Jacobs, D. L., and Toms, S. A. (2019). “Cladophialophora bantiana,” in Fungal Infections of the Central Nervous System, eds M. Turgut, S. Challa, and A. Akhaddar (Cham: Springer).
Palmeira, V. F., Goulart, F. R. V., Granato, M. Q., Alviano, D. S., Alviano, C. S., Kneipp, L. F., et al. (2018). Fonsecaea pedrosoi sclerotic cells: secretion of aspartic-type peptidase and susceptibility to peptidase inhibitors. Front. Microbiol. 9:1383. doi: 10.3389/fmicb.2018.01383
Perdoni, F., Falleni, M., Tosi, D., Cirasola, D., Romagnoli, S., Braidotti, P., et al. (2014). A histological procedure to study fungal infection in the wax moth Galleria mellonella. Eur. J. Histochem. 58:2428. doi: 10.4081/ejh.2014.2428
Piro, V. C., Faoro, H., Weiss, V. A., Steffens, M. B. R., Pedrosa, F. O., Souza, E. M., et al. (2014). FGAP: an automated gap closing tool. BMC Res. Notes 7:371. doi: 10.1186/1756-0500-7-371
Pockley, A. G. (2001). Heat shock proteins in health and disease: therapeutic targets or therapeutic agents? Expert Rev. Mol. Med. 3, 1–21. doi: 10.1017/S1462399401003556
Poyntner, C., Mirastschijski, U., Sterflinger, K., and Tafer, H. (2018). Transcriptome study of an Exophiala dermatitidis PKS1 mutant on an ex vivo skin model: Is melanin important for infection? Front. Microbiol. 9:1457. doi: 10.3389/fmicb.2018.01457
Prenafeta-Boldú, F. X., Summerbell, R., and Sybren De Hoog, G. (2006). Fungi growing on aromatic hydrocarbons: biotechnology’s unexpected encounter with biohazard? FEMS Microbiol. Rev. 30, 109–130. doi: 10.1111/j.1574-6976.2005.00007.x
Prenafeta-Boldú, F. X., Vervoort, J. J. M., Grotenhuis, T., and van Froenestijn, J. W. (2002). Substrate interactions during the biodegradation of benzene, toluene, ethylbenzene, and xylene (BTEX) hydrocarbons by the fungus Cladophialophora sp. strain T1. Appl. Environ. Microbiol. 68, 2660–2665. doi: 10.1128/AEM.68.6.2660-2665.2002
Queiroz-Telles, F., De Hoog, S., Santos, D. W. C. L., Salgado, C. G., Vicente, V. A., Bonifaz, A., et al. (2017). Chromoblastomycosis. Clin. Microbiol. Rev. 30, 233–276. doi: 10.1128/CMR.00032-16
Quevillon, E., Silventoinen, V., Pillai, S., Harte, N., Mulder, N., Apweiler, R., et al. (2005). InterProScan: protein domains identifier. Nucleic Acids Res. 33, 116–120. doi: 10.1093/nar/gki442
Raparia, K., Powell, S. Z., Cernoch, P., and Takei, H. (2010). Cerebral mycosis: 7-year retrospective series in a tertiary center. Neuropathology 30, 218–223. doi: 10.1111/j.1440-1789.2009.01067.x
Rawlings, N. D., Barrett, A. J., and Finn, R. (2016). Twenty years of the MEROPS database of proteolytic enzymes, their substrates and inhibitors. Nucleic Acids Res. 44, 343–350. doi: 10.1093/nar/gkv1118
Rodrigues, A. M., de Hoog, G. S., and Camargo, Z. P. (2015). Molecular diagnosis of pathogenic Sporothrix Species. PLoS Negl. Trop. Dis. 9:e0004190. doi: 10.1371/journal.pntd.0004190
Satow, M. M. (2005). “Screening” de Fungos Degradadores de Hidrocarbonetos Complexos. Master’s dissertation, Universidade Estadual Paulista, Rio Claro.
Schneider, G. X., Gomes, R. R., Bombassaro, A., Zamarchi, K., Voidaleski, M. F., Costa, F. F., et al. (2019). New molecular markers distinguishing fonsecaea agents of chromoblastomycosis. Mycopathologia 184, 493–504. doi: 10.1007/s11046-019-00359-2
Scorzoni, L., de Lucas, M. P., Mesa-Arango, A. C., Fusco-Almeida, A. M., Lozano, E., Cuenca-Estrella, M., et al. (2013). Antifungal efficacy during Candida krusei infection in non-conventional models correlates with the yeast in vitro susceptibility profile. PLoS One 8:e60047. doi: 10.1371/journal.pone.0060047
Seppey, M., Manni, M., and Zdobnov, E. M. (2019). BUSCO: assessing genome assembly and annotation completeness. Methods Mol. Biol. 1962, 227–245. doi: 10.1007/978-1-4939-9173-0_14
Seyedmousavi, S., Netea, M. G., Mouton, J. W., Melchers, W. J., Verweij, P. E., and de Hoog, G. S. (2014). Black yeasts and their filamentous relatives: principles of pathogenesis and host defense. Clin. Microbiol. Rev. 27, 527–542. doi: 10.1128/CMR.00093-13
Sideris, E., and Ge, L. (2018). Chromoblastomycosis in immunosuppressed petients. Med. J. Aust. 209:295. doi: 10.5694/mja18.00350
Sriranganadane, D., Waridel, P., Salamin, K., Reichard, U., Grouzmann, E., Neuhaus, J., et al. (2010). Aspergillus protein degradation pathways with different secreted protease sets at neutral and acidic pH. J. Proteome Res. 9, 3511–3519. doi: 10.1021/pr901202z
Stanley, P. (2016). What have we learned from glycosyltransferase knockouts in mice? J. Mol. Biol. 428, 3166–3182. doi: 10.1016/j.jmb.2016.03.025
Sterflinger, K., and Prillinger, H. J. (2001). Molecular taxonomy and biodiversity of rock fungal communities in an urban environment (Vienna, Austria). Antonie Van Leeuwenhoek 80, 275–286. doi: 10.1023/A:1013060308809
Stokes, W., Fuller, J., Meier-Stephenson, V., Remington, L., and Meatherall, B. L. (2017). Case report of cerebral phaeohyphomycosis caused by Fonsecaea monophora. Off. J. Assoc. Med. Microbiol. Infect. Dis. Canada 2, 86–92. doi: 10.3138/jammi.2.1.013
Stone, B. A., Jacobs, A. K., Hrmova, M., Burton, R. A., and Fincher, G. B. (2018). Biosynthesis of plant cell wall and related polysaccharides by enzymes of the GT2 and GT48 families. Annu. Plant Rev. 41, 109–165. doi: 10.1002/9781119312994.apr0434
Sudhadham, M., Prakitsin, S., Sivichai, S., Chaiyarat, R., Dorrestein, G. M., Menken, S. B. J., et al. (2008). The neurotropic black yeast Exophiala dermatitidis has a possible origin in the tropical rain forest. Stud. Micol. 61, 145–155. doi: 10.3114/sim.2008.61.15
Surash, S., Tyagi, A., Hoog, G. S., Zeng, J. S., Barton, R. C., and Hobson, R. P. (2005). Cerebral phaeohyphomycosis caused by Fonsecaea monophora. Med. Mycol. 43, 465–472. doi: 10.1080/13693780500220373
Takei, H., Goodman, J. C., and Powell, S. Z. (2007). Cerebral phaeohyphomycosis caused by Cladophialophora bantiana and Fonsecaea monophora: report of three cases. Clin. Neuropathol. 26, 21–27. doi: 10.5414/NPP26021
Tanabe, H., Kawasaki, M., Mochizuki, T., and Ishizaki, H. (2004). Species identification and strain typing of Fonsecaea pedrosoi using ribosomal RNA gene internal transcribed spacer regions. Nippon Ishinkin Gakkai Zasshi 45, 105–112. doi: 10.3314/jjmm.45.105
Teixeira, M. M., Moreno, L. F., Stielow, B. J., Muszewska, A., Hainaut, M., Gonzaga, L., et al. (2017). Exploring the genomic diversity of black yeasts and relatives (Chaetothyriales, Ascomycota). Stud. Mycol. 86, 1–28. doi: 10.1016/j.simyco.2017.01.001
Thirach, S., Cooper, C. R. J., Vanittanakom, P., and Vanittanakom, N. (2007). The copper, zinc superoxide dismutase gene of Penicillium marneffei: cloning, characterization, and differential expression during phase transition and macrophage infection. Med. Mycol. 45, 409–417. doi: 10.1080/13693780701381271
Thomas, E., Bertolotti, A., Barreau, A., Klisnick, J., Tournebize, P., Borgherini, G., et al. (2018). From phaeohyphomycosis to disseminated chromoblastomycosis: a retrospective study of infections caused by dematiaceous fungi. Med. Malad. Infect. 48, 278–285. doi: 10.1016/j.medmal.2017.09.011
Tiwari, S., Thakur, R., and Shankar, J. (2015). Role of heat-shock proteins in cellular function and in the biology of fungi. Biotechnol. Res. Int. 2015:132635. doi: 10.1155/2015/132635
Tronina, T., Bartmanska, A., Milczarek, M., Wietrzyk, J., Poplonski, J., Roj, E., et al. (2013). Antioxidant and antiproliferative activity of glycosides obtained by biotransformation of xanthohumol. Bioorg. Med. Chem. Lett. 23, 1957–1960. doi: 10.1016/j.bmcl.2013.02.031
Varghese, P., Jalal, M. J. A., Ahmad, S., Khan, Z., Johny, M., Mahadevan, P., et al. (2016). Cerebral phaeohyphomycosis caused by Fonsecaea monophora: first report from India. Int. J. Surg. Med. 2, 44–49. doi: 10.5455/ijsm.neurosurgery01
Varrot, A., Basheer, A. M., and Imberty, A. (2013). Fungal lectins: structure, function and potential applications. Curr. Opin. Struct. Biol. 23, 678–685. doi: 10.1016/j.sbi.2013.07.007
Vela-Corcía, D., Aditya Srivastava, D., Dafa-Berger, A., Rotem, N., and Barda, O. (2019). MFS transporter from Botrytis cinerea provides tolerance to glucosinolate-breakdown products and is required for pathogenicity. Nat. Commun. 10:2886. doi: 10.1038/s41467-019-10860-3
Vialle, R. A. (2013). SILA - Um Sistema para Anotação Automática de Genomas Utilizando Técnicas Independentes de Alinhamento. Master’s dissertation, Universidade Federal do Paraná, Curitiba.
Vialle, R. A., Pedrosa, F. O., Weiss, V. A., Guizelini, D., Tibaes, J. H., Marchaukoski, J. N., et al. (2016). RAFTS3: rapid alignment-free tool for sequences similarity search. bioRxiv. [Preprint]. doi: 10.1101/055269
Vicente, V. A., Attili-Angelis, D., Pie, M. R., Queiroz-Telles, F., Cruz, L. M., Najafzadeh, M. J., et al. (2008). Environmental isolation of black yeast-like fungi involved in human infection. Stud. Mycol. 61, 137–144. doi: 10.3114/sim.2008.61.14
Vicente, V. A., Najafzadeh, M. J., Sun, J., Gomes, R. R., Robl, D., Marques, S. G., et al. (2014). Environmental siblings of black agents of human chromoblastomycosis. Fungal Divers. 62, 1–17. doi: 10.1007/s13225-013-0246-5
Vicente, V. A., Orélis-Ribeiro, R., Najafzadeh, M. J., Sun, J., Guerra, R. S., Miesch, S., et al. (2012). Black yeast-like fungi associated with lethargic Crab disease (LCD) in the mangrove-land crab, Ucides cordatus (Ocypodidae). Vet. Microbiol. 158, 109–122. doi: 10.1016/j.vetmic.2012.01.031
Vicente, V. A., Weiss, V. A., Bombassaro, A., Moreno, L. F., Costa, F. F., Raittz, R. T., et al. (2017). Comparative genomics of sibling species of Fonsecaea associated with human chromoblastomycosis. Front. Microbiol. 8:1924. doi: 10.3389/fmicb.2017.01924
Wang, C., Xing, H., Jiang, X., Zeng, J., Liu, Z., Chen, J., et al. (2019). Cerebral phaeohyphomycosis caused by Exophiala dermatitidis in a Chinese CARD9-Deficient Patient: a case report and literature review. Front. Neurol. 10:938. doi: 10.3389/fneur.2019.00938
White, C., Laird, D. W., and Hughes, L. H. (2017). From carbon waste to carbon product: converting oxalate to polyhydroxybutyrate using a mixed microbial culture. J. Environ. Chem. Eng. 5, 2362–2365. doi: 10.1016/j.jece.2017.04.040
Woertz, J. R., Kinney, K. A., McIntosh, N. D. P., and Szaniszlo, P. J. (2001). Removal of toluene in a vapor phase bioreactor containing a strain of the dimorphic black yeast Exophiala lecaniicorni. Biotechnol. Bioeng. 75, 550–558. doi: 10.1002/bit.10066
Yang, L., Wang, Z., Lei, H., Chen, R., Wang, X., Peng, Y., et al. (2014). Neuroprotective glucosides of magnolol and honokiol from microbial-specific glycosylation. Tetrahedron 70, 8244–8251. doi: 10.1016/j.tet.2014.09.033
Yang, Y. T., Lee, S. T., Nai, Y. S., Kim, S., and Kim, J. S. (2016). Up-regulation of carbon metabolism-related glyoxylate cycle and toxin production in Beauveria bassiana JEF-007 during infection of bean bug, Riptortus pedestris (Hemiptera: Alydidae). Fungal Biol. 120, 1236–1248. doi: 10.1016/j.funbio.2016.07.008
Youngchim, S., Morris-Jones, R., Hay, R. J., and Hamilton, A. J. (2004). Production of melanin by Aspergillus fumigates. J. Med. Microbiol. 53, 175–181. doi: 10.1099/jmm.0.05421-0
Zhao, Z., Liu, H., Wang, C., and Xu, J. (2013). Correction: comparative analysis of fungal genomes reveals different plant cell wall degrading capacity in fungi. BMC Genomics 14:274. doi: 10.1186/1471-2164-15-6
Keywords: black fungi, cerebral infection, genome assembly, virulence, dissemination, pathology, chromoblastomycosis, neurotropism
Citation: Bombassaro A, Schneider GX, Costa FF, Leão ACR, Soley BS, Medeiros F, da Silva NM, Lima BJFS, Castro RJA, Bocca AL, Baura VA, Balsanelli E, Pankievicz VCS, Hrysay NMC, Scola RH, Moreno LF, Azevedo CMPS, Souza EM, Gomes RR, de Hoog S and Vicente VA (2020) Genomics and Virulence of Fonsecaea pugnacius, Agent of Disseminated Chromoblastomycosis. Front. Genet. 11:822. doi: 10.3389/fgene.2020.00822
Received: 13 November 2019; Accepted: 08 July 2020;
Published: 04 August 2020.
Edited by:
Bridget Marie Barker, Northern Arizona University, United StatesReviewed by:
Anna Muszewska, Institute of Biochemistry and Biophysics (PAN), PolandAylin Döğen, Mersin University, Turkey
Macit Ilkit, Çukurova University, Turkey
Hamid Badali, The University of Texas Health Science Center at San Antonio, United States
Copyright © 2020 Bombassaro, Schneider, Costa, Leão, Soley, Medeiros, da Silva, Lima, Castro, Bocca, Baura, Balsanelli, Pankievicz, Hrysay, Scola, Moreno, Azevedo, Souza, Gomes, de Hoog and Vicente. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sybren de Hoog, Sybren.deHoog@radboudumc.nl; s.hoog@wi.knaw.nl; Vânia A. Vicente, vaniava63@gmail.com