Clinical Spectrum, Molecular Characterization, Antifungal Susceptibility Testing of Exophiala spp. From India and Description of a Novel Exophiala Species, E. arunalokei sp. nov

Introduction

Exophiala species are dematiaceous (black-pigmented) saprophytic molds that are frequently isolated from hot and humid environments that are nutrient deficient but rich in hydrocarbons (Lian and De Hoog, 2010; Zhao et al., 2010; Woo et al., 2013). This genus is a member of the family Herpotrichiellaceae and encompasses >30 species. They are often referred to as black yeasts due to the early phenotypic character of their colonies (Matsumoto et al., 1987; Woo et al., 2013). Apart from isolation from environmental samples, they also colonize the human intestine and respiratory tract particularly in the airway of patients with cystic fibrosis (De Hoog et al., 2005; Kondori et al., 2011). However, they must not be disregarded as mere commensals or contaminants as they are also important opportunist pathogens causing systemic infections in immunosuppressed patients and subcutaneous or even fatal disseminated infections in otherwise healthy individuals (Fothergill, 1996; Zeng et al., 2007). Although Exophiala species are relatively uncommon, various reports of this genus as an etiologic agent of phaeohyphomycosis are being described from India in recent years (Sood et al., 2014; Venkateshwar et al., 2014; Bhardwaj et al., 2016; Haridasan et al., 2017). Despite the apparent virulence of these species and worrisome clinical scenario, there are no systematic studies on the Indian isolates of Exophiala species. Additionally, there is sparse information on the susceptibility profile of Indian isolates of Exophiala species against the commonly available antifunga agents. In this study, we aimed to describe the clinical spectrum of infections due to Exophiala from India along with phenotypic and molecular characterization and antifungal susceptibility of Exophiala species. We also describe a potentially novel Exophiala species from a case of subcutaneous phaeohyphomycosis. The novel species exhibited distinct phenotypic characteristics and three independent DNA regions were sequenced to show that it is distinct compared to all other Exophiala species. Based on these findings, we propose a new species, Exophiala arunalokei sp. nov.

Material and Methods

Ethics

The study was approved by the Institute ethical committee of Post Graduate Institute of Medical Education and Research (PGIMER), Chandigarh. Since all clinical data were retrospectively collected from patient records informed consent was waived.

Patients and Strains

The isolates of Exophiala species included in this study were retrieved from the National Culture Collection of Pathogenic Fungi (NCCPF), Chandigarh, India. These isolates were received from various parts of the country during a period of 12 years (2008 through 2020) and clinical data were collected from patient records and NCCPF requisition forms. The identification of all strains was confirmed by molecular methods as described below.

Conventional Identification

Phenotypic characterization was performed by subculture on malt extract agar (MEA, HiMedia) at 30°C for 7 days followed by subculture on Sabouraud’s dextrose agar (SDA, HiMedia) incubated at 25°C. Colony characters such as texture, color, and rate of growth were recorded. Microscopic features were studied by preparation of the lactophenol cotton blue (LCB) from slide culture (De Hoog et al., 2000). Also, the culture of Exophiala arunalokei sp. nov. was inoculated on potato dextrose agar and Oatmeal agar and incubated at 28, 37, and 45°C for up to 2 weeks to study temperature tolerance and to observe a difference in growth characteristics on different media.

Scanning Electron Microscopy (SEM)

The SEM was performed for one isolate each of E. jeanselmei, E. oligosperma, and E. arunalokei sp. nov. The coverslips from slide culture were used for SEM as previously described (Alves et al., 2013). Briefly, the culture growth over the coverslip was fixed using 2.5% glutaraldehyde followed by gentle washing with phosphate-buffered saline. Then gradient dehydration was performed using different concentration of ethanol (60, 70, 80% and absolute ethanol) followed by air-drying of the sample. The coverslip was loaded onto SEM stub with carbon tape and sputter-coated with platinum at 20 mA current for 30 s and placed onto the SEM stage for reading using the scanning electron microscope JSM I T300 InTouchScope (Jeol, Tokyo, Japan) with low vacuum mode (pressure range from 10 to 650 Pa). All samples were screened at ×250 magnification to identify the area of growth followed by subsequent visualization at ×1,000, ×2,500, and ×5,000 magnification.

Molecular Identification

For DNA extraction, freshly grown culture growth was transferred to mortar and pestle followed by treatment with liquid nitrogen. Crushed mycelia were suspended in 500 µl of lysis buffer followed by incubation at 56°C for 5 min. DNA was then extracted by the phenol-chloroform-isoamyl alcohol (25:24:1) (PCI) method as per the methods described previously (Shivaprakash et al., 2011). Briefly, the PCI mixture (Sigma Aldrich Inc., MO, USA) was added and mixed well with the lysed specimen followed by centrifugation at 13,000 rpm for 15 min. Then we transferred the aqueous phase to a sterile microcentrifuge tube and DNA was precipitated with absolute ethanol at −20°C overnight. DNA was then centrifuged followed by 70% ethanol wash and then air-dried. It was suspended by thoroughly mixing in 50 µl of sterile, deionized water (Milli Q, Millipore systems) and the yield and quality of the DNA were assessed by spectrophotometric measurement at a wavelength of 260 nm (NanoDrop 2000, Thermo Scientific, Wilmington, DE, USA). Molecular identification was done by sequencing the internal transcribed spacer (ITS) region of ribosomal DNA as described previously (White et al., 1990). The potentially novel Exophiala spp was confirmed by sequencing of the β-tubulin gene and β-actin gene using previously published primers (Woo et al., 2013).

Phylogenetic Characterization

For identification, the ITS sequences of all isolates were aligned with reference sequences at both the online sequence alignment portals, GenBank and International Society of Human and Animal Mycology (ISHAM)-ITS reference DNA barcoding database. The sequences and their relative strains representing each Exophiala species obtained from GenBank were used for alignment by the Molecular Evolutionarily Genetics Analysis (MEGA) Version 7.0 software (Kumar et al., 2016). A phylogenetic tree was constructed by neighbor joining (NJ) analysis and the evolutionary relationship among individual isolates was described as horizontal branches. The GenBank accession numbers of all the ITS gene sequences of included Exophiala strains are listed in Supplementary Table S1.

Matrix Assisted Laser Desorption–Time of Flight (MALDI-TOF)

The fresh growth of representative isolates belonging to each sequence confirmed species of Exophiala were used for off plate protein extraction using formic acid as per published protocol (Paul et al., 2018). For the MALDI-TOF analysis, six isolates were included (five clinical isolates and one reference isolate of Exophiala spinifera, ID B2762 CDC, Atlanta). E. mesophiala isolate could not be revived and was therefore not included in the MALDI-TOF analysis. MALDI-TOF identification was carried out on MALDI Microflex LT mass spectrometer (Bruker Daltonik GmbH, Bremen, Germany). The spectra of the isolates were added to the Bruker MALDI Biotyper database as described previously (Paul et al., 2018). Briefly, ionized spectral mass ranging from 2 to 20 kDa were acquired at a laser frequency of 60 Hz and MALDI BioTyper v3.1 was used for automated smoothing, normalization, baseline subtraction, and peak selection of different spectra. MALDI-TOF MS-based typing was performed using principal component analysis (PCA) for common species and a minimum of three MSPs (main spectrum profile) were used to construct a dendrogram using MALDI Biotyper software (MALDI Biotyper 3.1, Bruker Daltonics, MC, Italy) after smoothening, baseline subtraction, and normalization of spectra.

Antifungal Susceptibility Testing

Antifungal susceptibility testing was performed for amphotericin B, fluconazole, itraconazole, voriconazole, posaconazole, isavuconazole, caspofungin, and anidulafungin as per M38A3 protocol of Clinical Laboratory Standards Institute (CLSI, 2017). Briefly, the antifungals were diluted in the standard RPMI 1640 medium and dispensed into 96-well microdilution plates. The final concentration was 0.016–16 μg/ml for amphotericin-B, voriconazole, itraconazole, posaconazole, isavuconazole, and caspofungin; 0.08–8 μg/ml for anidulafungin and 0.063–64 μg/ml for fluconazole. Inoculum suspension in sterile saline with tween 40 was prepared from all isolates grown on SDA plates for 7 days. They were transferred to sterile tubes and optical density (OD) at 530 nm was adjusted to range from 0.09 to 0.13. The final suspension of the isolates ranged from 0.4 × 10⁴ to 3.1 × 10⁴ as confirmed hemocytometer counting and was diluted (1:50) in RPMI 1640 for final use. Microtiter plates were incubated at 35°C and the minimum inhibitory concentration (MIC) endpoints were noted (CLSI, 2017).

Literature Review

We conducted an electronic search on Medline and Cochrane Library for studies in the English language published up to 1^st January 2021. The search was based on the Mesh terms for the keywords (Exophiala) AND (India). Only articles reporting human infections were included. We included case reports, brief communications, letter to the editor, and published abstracts. The list of references in the original articles and reviews were searched manually searched to include any additional studies missed by the electronic search. The data extraction was performed using predesigned excel data extraction sheet and the author’s name, year of publication, state/region of study, patient’s clinical and demographic details, mycological findings, treatment, and clinical outcome was documented.

Results

A total of 24 isolates of Exophiala species from 23 patients could be revived from our culture repository. This included E. dermatitidis (n = 8 isolates), E. jeanselmei (n = 7 isolates), E. spinifera (n = 5 isolates), E. mesophila (n = 1 isolate), E. oligosperma (n = 1 isolate), E. xenobiotica (n = 1 isolate), and E. arunalokei sp. nov. (n = 1 isolate). The clinical profile and microbiological findings of cases are shown in Supplementary Table S1 and Figure 1.

FIGURE 1

Figure 1 Representative clinical images from subcutaneous infection due to Exophiala species. (A). Case no. NCCPF 106025 multiple subcutaneous nodular lesions over dorsum of left-hand (B). Case no. NCCPF 106030 nasal polyps (arrow) in patient with allergic fungal rhino sinusitis (C). Case no. NCCPF 106028 single subcutaneous cystic swelling skin over chest (left side).

Phenotypic Identification

The gross and microscopic morphological characteristics of the representative isolates from various Exophiala species are shown in Figure 2. Salient identification features of Exophiala species prevalent in India is provided in Table 1.

FIGURE 2

Figure 2 Microscopic findings on lactophenol cotton blue mount from slide culture and gross morphology on obverse surface of Sabouraud’s dextrose agar after 3–7 days at 25°C (inset). (A) E. dermatitidis, (B) E. mesophile, (C) E. jeanselmei, (D) E. spinifera, (E) E. xenobiotica, (F) E. mesophila.

TABLE 1

Table 1 Macroscopic and microscopic details of various Exophiala species characterized in the present study.

E. arunalokei sp. nov.

Growth on SDA was black and yeast-like initially, after 7 days of incubation at 25°C filamentous growth appearing greenish-gray and velvety of ~2 mm in diameter was noted with black pigmentation on reverse surface. Microscopic examination of the slide culture revealed hyaline hyphae and few hyaline yeast-like cells: ovoidal, 4–6 × 2–3.6 µm in size. Vegetative hyphae were subhyaline, olive to black-brown, branched, and septate (1–2.5 µm in diameter). Cylindrical annellophores (6–11.2 × 1–3.2 µm) with many predominantly polyblastic oval secondary conidia were seen terminally with a slight prominence at the tip (less evident than seen in E. spinifera) or laterally from hyphae. Annelloconidia were hyaline single-celled, oval, 2–4 × 1–2 µm in size. Even on SEM, the size of the conidiophores was 7.1–11.36 × 2.0–2.7 µm and the conidia ranged from 2.2–3.93 × 1.5–2.09 µm. Figure 3 shows the phenotypic characteristics of E. arunalokei sp. nov. as seen on culture, light microscopy, and scanning electron microscopy. E. arunalokei sp. nov. grew rapidly on PDA and OMA compared to SDA agar. Optimal growth of E. arunalokei sp. nov. was seen at 28°C, with less growth at 37°C and no growth at 45°C.

FIGURE 3

Figure 3 Phenotypic characteristics of E arunalokei sp. nov. Gross morphology on obverse surface culture growth after 14 days at 28°C on (A) Potato dextrose agar, (B) Oatmeal agar, and (C) Sabouraud’s dextrose agar. (D) Microscopic findings from slide culture (day 10) on lactophenol cotton blue mount showing branched annellophore, (E) terminal prominence (arrow), (F) LCB from slide culture on day 21. Scanning electron microscopy (SEM) of (G) E oligosperma, (H) E. jeanselmei, and (I) E. arunalokei sp. nov. day 10, ×2,500; (J) ×5,000 (day 10); (K) day 21, ×2,500, (L) day 21, ×2,500 with size measurements.

The amplicon size of the β-tubulin and β-actin genes of E. arunalokei sp. nov. and the reference strains were ~300 and ~500 bp, respectively. E. arunalokei sp. nov. was found to be most closely related to E. xenobiotica [90% similarity with β-actin gene (JN625237.1) and 80% similarity with β-tubulin gene (JN625236.1)] on phylogenetic analysis (Figure 4) concurrent with the findings of ITS gene sequencing (91% similarity with E. xenobiotica CBS204.50, KP132163).

FIGURE 4

Figure 4 Phylogenetic tree showing the relationship of E. arunalokei to closely related species using (A) β-actin gene and (B) β-tubulin gene sequencing data by using the Neighbor-Joining method clustered together in the bootstrap test (1,000 replicates); conducted in MEGA 7.

Taxonomy of Exophiala arunalokei sp. nov.: Exophiala arunalokei Singh and Rudramurthy sp. nov. MycoBank accession no. is awaited. Teleomorph: unknown. Known distribution: India. Etymology: a.ru.na.lo’ke.i. N.L. masc. gen. n. arunalokei, of Arunalokei in honor of Professor Arunaloke Chakrabarti, for his many academic and scientific contributions in the field of medical mycology. Specimen examined: in India; from a subcutaneous lesion on the face of a human presenting with subcutaneous phaeohyphomycosis in 2020. Based on the microscopic morphological features and sequence analysis of the ITS region, beta-actin and beta-tubulin region presented above, the isolate is proposed as the strain of a novel species in the genus Exophiala. Holotype: dried culture at the NCCPF, Chandigarh, India (E. arunalokei NCCPF 106033).

PCR and Sequence Analysis

The ITS sequences of the 24 Indian isolates of Exophiala species recovered from our collection measures ~600 bp. Figure 5 depicts the phylogram prepared from the ITS sequences of all 24 clinical isolates of the present study and 34 strains from various species of Exophiala reported from clinical and environmental samples. The aligned sequences of various previously published type strains E. jeanselmei (Types A, B, and C) and E. dermatitidis (Types A0, A1, A2, and B) were also used to study the evolutionary relationship of E. jeanselmei and E. dermatitidis isolates in our study with previously described type isolates is shown in Supplementary Figures 1 and 2 (Suh et al., 2012). Two isolates of E. jeanselmei (NCCPF 106028 and 106021) from cases of subcutaneous infection showed sequence similarity with E. jeanselmei Type B while the remaining five isolates clustered separately but close to E. jeanselmei Type A, suggesting the possibility of a distinct type or subtype. The phylogram of ITS sequences of E. dermatitidis isolates showed that majority of the isolates clustered with E. dermatitidis Type B. However, the isolates number NCCPF 106012 and NCCPF 106013 were both recovered from the same patient with endocarditis, and they clustered with E. dermatitidis type A isolates.

FIGURE 5

Figure 5 Evolutionary relationships of 24 isolates of Exophialia species with Neighbor-joining algorithm with 1,000 bootstrap replicates. Phaeococcomyces mexicanus CBS 127164 is taken as out-group. ITS, internal transcribed spacer.

MALDI-TOF Analysis

None of the six isolates used in initial evaluation could be correctly identified using the existing Bruker database. A total of 24 spectra for each isolate were subjected to smoothing, normalization, baseline subtraction and added to our in-house database to expand the spectral library. The novel species E. arunalokei sp. nov., had a distinct spectral profile which was most similar to E. jeanselmei (Supplementary Figure 3A). Subsequently, E. dermatitidis (n = 7), E. jeanselmei (n = 5), E. spinifera (n = 6), and one isolate each of E. oligosperma, E. xenobiotica, and E. arunalokei sp. nov., were identified using the expanded in-house database and all could be correctly identified at a log score value of ≥1.7. Since only one isolate each of E. oligosperma, E. xenobiotica, and E. arunalokei sp. nov. were available the same isolate was used during in-house spectral database expansion and validation. The PCA dendrogram obtained by MALDI-TOF MS of the common Exophiala species revealed clusters which concurred with the sequencing results and dendrogram with representative species is shown in Supplementary Figure 3B.

Antifungal Susceptibility Testing (AFST)

A total of 24 isolates were subjected to AFST and the MIC (mg/L) ranged from 0.03 to 4 for amphotericin, 0.03 to 2 for all tested triazoles (itraconazole, voriconazole, and posaconazole), 0.125 to 16 for micafungin and caspofungin, and 0.03 to 16 for anidulafungin. The MIC of echinocandins was >4 mg/L in 62.5% (n = 5) isolates of E. dermatitidis and in one (16.7%) isolate of E. spinifera. The MIC of voriconazole was high (>2 mg/L) in an isolate each of E. oligosperma and E. mesophila (Table 2). The isolate of E. arunalokei sp. nov., had low MICs for all antifungals tested with a high MIC of 4 mg/L for amphotericin B.

TABLE 2

Table 2 Antifungal susceptibility testing of clinical isolates of Exophiala species.

Review of Literature and Data Analysis

A total of 27 articles were retrieved after an online and manual literature search. Of these, three were excluded based on title and abstract review. Eight isolates of seven published studies had been submitted at NCCPF and were included in this study as shown previously in Table 1. Duplication of publication was seen in two studies and it was counted as a single case (Capoor et al., 2007; Khan et al., 2007). The analysis of case details of 39 cases of Exophiala infections [23 included in the present study (Supplementary Table S1) and 16 from the literature review (Supplementary Table S2)] was performed.

A total of 39 cases due to Exophiala species have been reported from India with a median age of 35 years (range 8 to 73 years). The most common species reported was E. jeanselmei (n = 16, 41.02%), followed by E. dermatitidis (n = 10, 25.6%), E. spinifera (n = 9, 23.1%), E. oligosperma (n = 1, 2.6%), E. mesophila (n = 1, 2.6%), E. xenobiotica (n = 1, 2.6%), and E. arunalokei (n = 1, 2.6%). The most common clinical presentation was subcutaneous phaeohyphomycosis (n = 16, 41.02%), followed by invasive/disseminated infection (n = 9, 23.1%), eumycetoma (n = 5, 12.8%), skin infection (n = 3, 7.69%), allergic fungal rhinosinusitis (n = 1, 2.6%), and endophthalmitis (n = 1, 2.6%), while the details of four cases were not available.

The most common risk factor associated with infection (details available in 33) was trauma (n = 11, 33.3%), post-renal transplant immunosuppression (n = 4, 12.1%), history of surgical procedures (n = 3, 9.1%), diabetes (n = 3, 9.1%), malnourishment (n = 2, 6.1%), tuberculosis (n = 1, 3.03%), and chronic hepatitis B (n = 1, 3.03%) while 8 (24.2%) were immunocompetent. In case of infection due to E. dermatitidis, a significant association with a history of surgical procedures was seen compared to other species (30 vs. 0%; p-value: 0.013). E. mesophila infection was significantly associated with a history of tuberculosis compared to other species (100 vs. 0%; p-value: 0.026). Analysis of species wise clinical presentation revealed a significant association of E. dermatitidis infection with the invasive disease compared to other species (50 vs. 10.7%, p-value: 0.032).

The details of patient management were available for 24 cases. Antifungals were used in 19 cases of which the majority received itraconazole (n = 7), followed by voriconazole (n = 4) and ketoconazole (n = 3) while a combination of different antifungals such as amphotericin, voriconazole, itraconazole, fluconazole, and terbinafine was used in six cases. Surgical intervention involving the drainage or debridement of the infected site was performed in 11 and cryotherapy in one case. Outcome details were available for 26 cases of which 24 (92%) recovered while 2 (7.69%) died.

Discussion

Exophiala species are difficult to characterize using phenotypic tests and different molecular methods have been recommended for accurate species identification (Zeng et al., 2007; Woo et al., 2013; Silva et al., 2017). Disagreements while comparing morphological and molecular identification methods have been previously reported for this genus (Zeng et al., 2007). In the present study, we describe the largest Indian casuistic of Exophiala spp. isolated from clinical samples and characterized phenotypically and by ITS sequencing for accurate species-level identification. Since molecular characterization of isolates is not always possible at hospital based laboratories we also performed the MALDI-TOF MS base identification of our isolates and none of them could be correctly identified up to species level using the existing Bruker database. Hence, we expanded our in-house library of protein spectra.

The clinical details of cases included in the present study and the data recovered by a systematic review of previous Indian studies were analyzed to elaborate on the spectrum of clinical Exophiala species from India. A large diversity of species was demonstrated, with E. jeanselmei being the most prevalent followed by E. dermatitidis. This is also the first Indian study, to report a case of human infection due to E. xenobiotica. Phylogenetic analysis of sequencing data from three independent regions i.e., ITS and two housekeeping genes (β-tubulin, and β-actin gene) unambiguously showed that the isolate NCCPF 106033 formed branches distinct from other species and we propose the novel species E. arunalokei. This species was related most closely to E. xenobiotica evidenced in all phylogenetic trees with high bootstrap support. Compared to E. jeanselmei, E. arunalokei grew slowly and had smaller conidia with branched conidiophores with a slight prominence at the end. E. arunalokei culture growth was velvety compared to E. xenobiotica and conidia were also smaller in size.

The sequence typing of all isolates of E. jeanselmei and E. dermatitidis were performed as reported previously by Suh et al. (Suh et al., 2012) and most Indian isolates of E. jeanselmei were found to be related to E. jeanselmei type A, while E. dermatitidis were Type B. Previous studies analyzing ITS rDNA of E. dermatitidis have revealed that type B is commonly constituted by environmental strains while type A isolate are prevalent in the clinical setting (Matos et al., 2003; Suh et al., 2012). Interestingly, while all E. dermatitidis isolates in our study caused invasive disease the isolates infecting immunocompetent hosts clustered in type B. We also identified two E. dermatitidis Type A2 strains in patients with post-renal transplant immunosuppression. Previously Type A2 has been reported from the sputum of a cystic fibrosis patient from Germany (E. dermatitidis CBS 149.90 AF050268). E. jeanselmei Type B has been previously reported from subcutaneous infections in humans from Korea and England (Suh et al., 2012). We too found that the two Type B E. jeanselmei isolates in our study, caused subcutaneous infection.

Although clinical breakpoints are not available for Exophiala spp. are not available the AFST results show low MICs to all antifungals amongst all the described species, except one isolate each of E. dermatitidis and E. oligosperma which show high voriconazole MIC. Thus, monitoring the trends of antifungal susceptibility is essential to anticipate the development of antifungal resistance in the future and modify the patient management protocols accordingly. Azoles, especially itraconazole and posaconazole and amphotericin B are the most active antifungals. However, there are no established breakpoints for Exophiala species and reports of AFST profile remain sparse (Zeng et al., 2007; Revankar and Sutton, 2010; Silva et al., 2017).

Association of E. dermatitidis infection with the invasive disease is concordant with previous reports where this species has been reported as the most virulent in this genus (Revankar and Sutton, 2010; Silva et al., 2017). In case of infection due to E. dermatitidis, a significant association with a history of surgical procedures was seen compared to other species (p-value: 0.013) which could indicate a possibility of nosocomial acquisition. Although there was a statistical association between tuberculosis and E. mesophila infection (p value: 0.026), since it was only one case, an association at this juncture is debatable. Subclinical airway colonization may be involved in such cases as observed in cases with cystic fibrosis (Kondori et al., 2011). However, no such cases were seen in the present study suggesting the need for prospective studies to assess the prevalence of Exophiala spp. in respiratory secretions from such patients in our region. The wide spectrum of infections due to Exophiala species in both immunocompetent and immunocompromised patients is evident from this study. In addition to antifungal therapy, concomitant surgical management should be performed whenever possible.

In a study from Hong Kong, characterization from 12 clinical isolates of Exophiala spp. revealed that superficial and subcutaneous infections were more common and the most common species recovered was E. oligosperma (n = 3; two from cases of onychomycosis and one from case of pneumonia) followed by E. jeanselmei (n = 2; both subcutaneous infection), E. lecanii-corni (n = 2; both subcutaneous infection), E. bergeri (n = 1; onychomycosis), E. xenobiotica (n = 1, peritonitis), and a novel species E. hongkongensis (case of onychomycosis) (Woo et al., 2013). Exophiala strains (n = 185) were characterized in another study from the United States of America, where 50.8% were seen in cutaneous/subcutaneous infection while 39.9% were found to be isolated from deep infections. E. dermatitidis, E. xenobiotica, and E. oligosperma comprised more than two-thirds of the isolates and while systemic infections were mainly caused by E. dermatitidis, E. oligosperma, and E. xenobiotica, E. spinifera, E. jeanselmei, and E. mesophila caused cutaneous/subcutaneous infections (Zeng et al., 2007).

This, to the best of our knowledge, is the first study from India elucidating the molecular and clinical characteristics of Exophiala species. We find that the combination of ITS sequence phylogenetic analysis can supplement traditional morphological identification of Exophiala species as well as the evaluation of its distribution, determination of types and subtypes. It also enabled us to describe a novel species E. arunalokei. Apart from sequence based identification, MALDI-TOF can also serve to be quite useful provided the spectral library is enhanced. We have expanded our in-house database using the spectra acquired from representative isolates in this study however, further expansion and validation of this library is essential for identification with better log scores. The antifungal susceptibility data presented here can also provide an insight into the complete picture of Exophiala species from India. This can help in establishment of epidemiological cutoffs for MIC interpretation and should be further investigated in multicenter studies.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author Contributions

SS, SR, and AP designed and supervised the study. BH, RI, VH, and AS provided the isolates and performed data collection. SS, SA, and SG carried out the laboratory experiments. SS carried out the data analysis and wrote the manuscript. SR, AG, and HK revised the manuscript. All authors contributed to the article and approved the submitted version.

Conflict of Interest

Author BH was employed by the company Hitech Laboratories.

The remaining 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.

Acknowledgments

We would like to acknowledge the Department of Science & Technology (DST)- Science and Engineering Research Board (SERB), New Delhi, India and Council of Scientific and Industrial Research, New Delhi, India. We would also like to acknowledge Hansraj Chaudhary, MSc for technical assistance.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2021.686120/full#supplementary-material

References

Alves, E., Lucas, G. C., Pozza, E. A., de Carvalho Alves, M. (2013). “Scanning Electron Microscopy for Fungal Sample Examination,” in Laboratory Protocols in Fungal Biology. Eds. Gupta, V., Tuohy, M., Ayyachamy, M., Turner, K., O’ Donovan, A. (New York, NY: Springer New York), 133–150. doi: 10.1007/978-1-4614-2356-0_8

CrossRef Full Text | Google Scholar

Bhardwaj, S., Capoor, M. R., Kolte, S., Purohit, G., Dawson, L., Gupta, K., et al. (2016). Phaeohyphomycosis Due to Exophiala Jeanselmei: An Emerging Pathogen in India—Case Report and Review. Mycopathologia 181, 279–284. doi: 10.1007/s11046-015-9955-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Capoor, M., Khanna, G., Nair, D., Hasan, A., Rajni, Deb, M., et al. (2007). Eumycetoma Pedis Due to Exophiala Jeanselmei. Indian J. Med. Microbiol. 25, 155. doi: 10.4103/0255-0857.32726

PubMed Abstract | CrossRef Full Text | Google Scholar

CLSI (2017). Clinical and Laboratory Standards Institute (Clsi) M38-Ed3 Filamentous Fungi Antifungal Susceptibility Test. Available at: https://clsi.org/standards/products/microbiology/documents/m38/.

Google Scholar

De Hoog, G. S., Guarro, J., Gene, J., Figueras, M. (2000). “Atlas of Clinical Fungi,” in Centraalbureau Voor Schimmelcultures/Universitat Rovira I Virgili, 2nd ed, vol. 00. (Utrecht: The Netherlands.Publisher), Vii + 1126pp.

Google Scholar

De Hoog, G. S., Matos, T., Sudhadham, M., Luijsterburg, K. F., Haase, G. (2005). Intestinal Prevalence of the Neurotropic Black Yeast Exophiala (Wangiella) Dermatitidis in Healthy and Impaired Individuals. Mycoses 48, 142–145. doi: 10.1111/j.1439-0507.2004.01083.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Fothergill, A. W. (1996). Identification of Dematiaceous Fungi and Their Role in Human Disease. Clin. Infect. Dis. 22 Suppl 2, S179–S184. doi: 10.1093/clinids/22.Supplement_2.S179

CrossRef Full Text | Google Scholar

Haridasan, S., Parameswaran, S., Bheemanathi, S. H., Chandrasekhar, L., Suseela, B. B., Singh, R., et al. (2017). Subcutaneous Phaeohyphomycosis in Kidney Transplant Recipients: A Series of Seven Cases. Transpl. Infect. Dis. 19 (6), e12788. doi: 10.1111/tid.12788

CrossRef Full Text | Google Scholar

Khan, S. A., Hasan, A. S., Capoor, M. R., Varshney, M. K., Trikha, V. (2007). Calcaneal Osteomyelitis Caused by Exophiala Jeanselmei in an Immunocompetent Child. J. Bone Jt. Surg. 89, 859–862. doi: 10.2106/00004623-200704000-00024

CrossRef Full Text | Google Scholar

Kondori, N., Gilljam, M., Lindblad, A., Jönsson, B., Moore, E. R. B., Wennerås, C. (2011). High Rate of Exophiala Dermatitidis Recovery in the Airways of Patients With Cystic Fibrosis Is Associated With Pancreatic Insufficiency. J. Clin. Microbiol. 49, 1004–1009. doi: 10.1128/JCM.01899-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, S., Stecher, G., Tamura, K. (2016). Mega7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 33, 1870–1874. doi: 10.1093/molbev/msw054

PubMed Abstract | CrossRef Full Text | Google Scholar

Lian, X., De Hoog, G. S. (2010). Indoor Wet Cells Harbour Melanized Agents of Cutaneous Infection. Med. Mycol. 48, 622–628. doi: 10.3109/13693780903405774

PubMed Abstract | CrossRef Full Text | Google Scholar

Matos, T., Haase, G., Gerrits van den Ende, A. H. G., de Hoog, G. S. (2003). Molecular Diversity of Oligotrophic and Neurotropic Members of the Black Yeast Genus Exophiala, With Accent on E. Dermatitidis. Antonie. Van. Leeuwenhoek. 83, 293–303. doi: 10.1023/A:1023373329502

PubMed Abstract | CrossRef Full Text | Google Scholar

Matsumoto, T., Padhye, A. A., Ajello, L. (1987). Medical Significance of the So-Called Black Yeasts. Eur. J. Epidemiol. 3, 87–95. doi: 10.1007/BF00239744

PubMed Abstract | CrossRef Full Text | Google Scholar

Paul, S., Singh, P., Sharma, S., Prasad, G. S., Rudramurthy, S., Chakrabarti, A., et al. (2018). Maldi-Tof MS-Based Identification of Melanized Fungi is Faster and Reliable After the Expansion of In-House Database. Proteomics Clin. Appl. 13, 1800070. doi: 10.1002/prca.201800070

CrossRef Full Text | Google Scholar

Revankar, S. G., Sutton, D. A. (2010). Melanized Fungi in Human Disease. Clin. Microbiol. Rev. 23, 884–928. doi: 10.1128/CMR.00019-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Shivaprakash, M. R., Appannanavar, S. B., Dhaliwal, M., Gupta, A., Gupta, S., Gupta, A., et al. (2011). Colletotrichum Truncatum: An Unusual Pathogen Causing Mycotic Keratitis and Endophthalmitis. J. Clin. Microbiol. 49, 2894–2898. doi: 10.1128/JCM.00151-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Silva, W. C., Gonçalves, S. S., Santos, D. W. C. L., Padovan, A. C. B., Bizerra, F. C., Melo, A. S. A. (2017). Species Diversity, Antifungal Susceptibility and Phenotypic and Genotypic Characterisation of Exophiala Spp. Infecting Patients in Different Medical Centres in Brazil. Mycoses 60, 328–337. doi: 10.1111/myc.12597

PubMed Abstract | CrossRef Full Text | Google Scholar

Sood, S., Vaid, V. K., Sharma, M., Bhartiya, H. (2014). Cerebral Phaeohyphomycosis by Exophiala Dermatitidis. Indian J. Med. Microbiol. 32, 188–190. doi: 10.4103/0255-0857.129830

PubMed Abstract | CrossRef Full Text | Google Scholar

Suh, M. K., Lee, H. C., Kim, D. M., Ha, G. Y., Choi, J. S. (2012). Molecular Phylogenetics of Exophiala Species Isolated From Korea. Ann. Dermatol. 24, 287–294. doi: 10.5021/ad.2012.24.3.287

PubMed Abstract | CrossRef Full Text | Google Scholar

Venkateshwar, S., Ambroise, M. M., Asir, G. J., Mudhigeti, N., Ramdas, A., Authy, K., et al. (2014). A Rare Case Report of Subcutaneous Phaeohyphomycotic Cyst Caused by Exophiala Oligosperma in an Immunocompetent Host With Literature Review. Mycopathologia 178, 117–121. doi: 10.1007/s11046-014-9762-4

PubMed Abstract | CrossRef Full Text | Google Scholar

White, T. J., Bruns, S., Lee, S., Taylor, J. (1990). “Amplification and Direct Sequencing of Fungal Ribosomal RNA Genes for Phylogenetics,” in Pcr Protocols: A Guide to Methods and Applications Academic Press (Academic Press; Spi edition). 315–322.

Google Scholar

Woo, P. C. Y., Ngan, A. H. Y., Tsang, C. C. C., Ling, I. W. H., Chan, J. F. W., Leung, S. Y., et al. (2013). Clinical Spectrum of Exophiala Infections and a Novel Exophiala Species, Exophiala Hongkongensis. J. Clin. Microbiol. 51, 260–267. doi: 10.1128/JCM.02336-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, J. S., Sutton, D. A., Fothergill, A. W., Rinaldi, M. G., Harrak, M. J., De Hoog, G. S. (2007). Spectrum of Clinically Relevant Exophiala Species in the United States. J. Clin. Microbiol. 45, 3713–3720. doi: 10.1128/JCM.02012-06

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, J., Zeng, J., de Hoog, G. S., Attili-Angelis, D., Prenafeta-Boldú, F. X. (2010). Isolation and Identification of Black Yeasts by Enrichment on Atmospheres of Monoaromatic Hydrocarbons. Microb. Ecol. 60, 149–156. doi: 10.1007/s00248-010-9651-4

PubMed Abstract | CrossRef Full Text | Google Scholar