Strain and Culture Conditions
A total of 47 strains belonging to 15 taxa of dermatophytes were obtained from both human and animal dermatophytosis cases as well as from culture collections (Table 1). For comparison, three non-dermatophyte fungal species that are frequently found on human skin as contaminants or pathogens were selected: Aspergillus fumigatus CCF 3522, Paecilomyces variotii CCF 3230 and Scopulariopsis brevicaulis CCF 6436. All strains were maintained on Malt Extract Agar (MEA; Sigma-Aldrich, St. Louis, USA) at 6°C.
Table 1
List of strains examined in this study
Species | Strain designations1 | Year of isolation / Origin | Analysis type | GenBank/ENA/DDBJ accession numbers |
| ITS | tubb | tef1a |
Trichophyton benhamiae var. benhamiae | IHEM 4710 | unknown/ USA | GC-MS | LR794129 | LR794285 | LR794260 |
| CCF 6486 | 1995/ USA | GC-MS | LR794130 | LR794286 | LR794261 |
| IHEM 3287 | 1970/ Belgium | GC-MS | LR794129 | - | LR794260 |
Trichophyton benhamiae var. luteum | CLIS 6433/20 | 2020/ Czechia | GC-MS; GCxGC-TOFMS | OR619675 | OR792381 | OR825132 |
| CLIS 9200/20 | 2020/ Czechia | GC-MS; GCxGC-TOFMS | OR619678 | OR792383 | OR825134 |
| CLIS 9106/20 | 2020/ Czechia | GC-MS; GCxGC-TOFMS | OR619677 | OR792382 | OR825133 |
| SK 2597/20 | 2020/ Czechia | GCxGC-TOFMS | LR794131 | - | - |
| CLIS 8587/19 | 2019/ Czechia | GCxGC-TOFMS | LR794131 | - | - |
Trichophyton europaeum | CLIS 5085/20 | 2020/ Czechia | GC-MS | OR619670 | OR792376 | OR825119 |
| CLIS 4898/20 | 2020/ Czechia | GC-MS | OR619668 | OR792374 | OR825117 |
| SK 586/20 | 2020/ Czechia | GC-MS | OR619693 | - | - |
Trichophyton japonicum | CLIS 9483/18 | 2018/ Czechia | GC-MS | OR619679 | - | - |
| L 1430/20 | 2020/ Czechia | GC-MS | OR619684 | OR792386 | OR825128 |
| SK 3908/19 | 2019/ Czechia | GC-MS | OR619692 | OR792388 | OR825130 |
Trichophyton persicum | CCF 6543 | 2017/ Iran | GC-MS | MG356864 | MG356864 | MW959139 |
Trichophyton erinacei | CCF 6504 | 2019/ Czechia | GC-MS; GCxGC-TOFMS | OR619666 | OR792373 | OR825116 |
| CCF 6399 | 2019/ Czechia | GC-MS; GCxGC-TOFMS | OR619665 | OR792372 | OR825115 |
| CLIS 3198/20 | 2020/ Czechia | GC-MS; GCxGC-TOFMS | MZ314454 | MZ320337 | MZ320327 |
| CCF 4472 | 2012/ Czechia | GCxGC-TOFMS | LR794136 | - | - |
| CCF 6563 | 2017/ Japan | GCxGC-TOFMS | OR619664 | - | OR825114 |
Trichophyton verrucosum | L 28/19 | 2019/ Czechia | GC-MS | OR619685 | OR792387 | OR825129 |
| CLIS 7396/20 | 2020/ Czechia | GC-MS | OR619676 | - | - |
Trichophyton quinckeanum | CLIS 9771/20 | 2020/ Czechia | GC-MS | OR619680 | OR792384 | OR825123 |
| CLIS 6248/20 | 2020/ Czechia | GC-MS | OR619674 | OR792380 | OR825122 |
Microsporum canis | D 106/20 | 2020/ Czechia | GC-MS | OR619681 | OR792390 | OR825135 |
| CLIS 2072/20 | 2020/ Czechia | GC-MS | OR619667 | - | - |
Trichophyton mentagrophytes | CLIS 6191/20 | 2020/ Czechia | GC-MS | OR619673 | OR792379 | OM568761 |
| PL 700/20 | 2020/ Czechia | GC-MS | OR619689 | OM314973 | OM568761 |
| CCF 6572 | 2017/ Czechia | GC-MS; GCxGC-TOFMS | OM283516 | OM568760 | OM314972 |
| CCF 6579 | 2017/ Czechia | GC-MS; GCxGC-TOFMS | OM283522 | OM568766 | OM314978 |
| CCF 6573 | 2016/ Czechia | GCxGC-TOFMS | OM283517 | OM314973 | OM568761 |
| CCF 6574 | 2016/ Czechia | GCxGC-TOFMS | OM283517 | OM314973 | OM568761 |
| CCF 6584 | 2015/ Czechia | GCxGC-TOFMS | OM283517 | OM314973 | OM568761 |
Trichophyton mentagrophytes var. indotineae | CCF 6599 | 2016/ Czechia | GC-MS | OM283512 | OM568780 | OM314968 |
| CCF 6597 | 2019/ Japan | GC-MS | OM283543 | OM568778 | OM314999 |
| CCF 6598 | 2019/ Japan | GC-MS | OM283544 | OM568779 | OM315000 |
Trichophyton rubrum | LY 8 | before 2016/ Denmark | GC-MS | OR619687 | - | - |
| LY 19 | before 2016/ Denmark | GC-MS | OR619686 | - | - |
| SK 701/20 | 2020/ Czechia | GC-MS | OR619694 | OR792389 | OR825131 |
| IDE 241/20 | 2020/ Czechia | GC-MS | OR619683 | OR792385 | OR825127 |
Trichophyton tonsurans | CLIS 5203/20 | 2020/ Czechia | GC-MS | OR619671 | OR792377 | OR825120 |
| CLIS 5011/20 | 2020/ Czechia | GC-MS | OR619669 | OR792375 | OR825118 |
| CLIS 5356/20 | 2020/ Czechia | GC-MS | OR619672 | OR792378 | OR825121 |
Epidermophyton floccosum | PL 782/20 | 2020/ Czechia | GC-MS | OR619688 | OR792392 | - |
Nannizzia gypsea | SK 2458/20 | 2020/ Czechia | GC-MS | OR619691 | OR792394 | OR825126 |
| SK 2316/20 | 2020/ Czechia | GC-MS | OR619690 | OR792393 | OR825125 |
| D 653/20 | 2020/ Czechia | GC-MS | OR619682 | OR792391 | OR825124 |
1Acronyms used for culture collection/clinician providing the strains: IHEM: BCCM/IHEM Fungi Collection: Human & Animal Health, Ixelles, Belgium (Dr. Dirk Stubbe); CCF: Culture Collection of Fungi, Department of Botany, Charles University, Prague, Czech Republic; CLIS: Public Health Institute in Ostrava, Ostrava, Czech Republic (RNDr. S. Dobiášová), SK: Department of Dermatology and Venereology, First Faculty of Medicine, General University Hospital in Prague, Charles University and General University Hospital in Prague, Prague, Czech Republic (MUDr. M. Skořepová, CSc.); LY/PL: Public Health Institute in Ústí Nad Labem, Prague, Czech Republic (P. Lysková, PhD.), L: Regional Hospital Liberec (MUDr. J. Doležalová); D: Laboratory of Medical Parasitology and Mycology, Hospital České Budějovice, České Budějovice, Czech Republic (MUDr. N. Mallátová); RK: Teikyo University Institute of Medical Mycology (TIMM), Tokyo, Japan (Prof. R. Kano); IR-ARM: Cellular and Molecular Research Center, Medical Basic Sciences Research Institute, Ahvaz Jundishapur University of Medical Sciences, Ahvaz, Iran (Assoc. Prof. A. Rezaei-Matehkolaei); IDE: Department of Microbiology, Palacký University in Olomouc, Faculty of Medicine (doc.MUDr. P. Hamal, PhD.); USA: University of Illinois at Urbana-Champaign, USA (Dr. Lois L. Hoyer) |
Cultivation and pre-treatment: GC-MS analysis
The strains were pre-cultivated in 500 ml Florence flasks with 100 ml of Czapek-Dox broth (Sigma-Aldrich, St. Louis, USA), 0.1 g of unprocessed autoclaved sheep wool, 0.088 g/l of thiamine hydrochloride (Sigma-Aldrich, St. Louis, USA) and 0.096 g/l of glycine (Sigma-Aldrich, St. Louis, USA). The flasks were shaken at 200 rpm and incubated at 37°C for 10 days. Control samples were prepared identically, but they were not inoculated. After vigorous shaking, 0.5 ml of the pre-cultivation media was inoculated into 20 ml glass vials. These vials contained a PTFE/silicone septum and an autoclaved mixture of 0.1 g of unprocessed sheep wool with 2 ml of 2% agar (Sigma-Aldrich, St. Louis, USA) for firming properties. The vials were then cultivated at 30°C for 30 days. This length of cultivation time was chosen because it provided better detectability of VOCs in the headspace compared to the two-week and three-week cultivations performed in the initial pilot study (data not shown). Three biological replicates per each strain were studied.
Cultivation and pre-treatment: GC×GC-TOF MS analysis
Compared to GC-MS analysis, the dataset was reduced to focus on a more detailed analysis of three specific taxa due to computational requirements. The zoophilic species T. erinacei, T. mentagrophytes var. mentagrophytes, and T. benhamiae var. luteum were chosen based on their significance in veterinary medicine. These species were cultivated on 60 mm plates using a 1.5% MEA medium for one week. Afterwards, 2 ml of 0.1% Tween 80 (Lach-Ner, Neratovice, Czech Republic) PBS (Phosphate Buffer Saline; Sigma-Aldrich, St. Louis, USA) was added to the plates. The plates were then swabbed with a sterile cotton swab and resuspended in 6 ml of sterile 0.1% Tween 80 PBS. Each suspension was vortexed and diluted to approximately the same density. The suspension was then left at room temperature for 5 minutes to allow the heavier hyphal elements to settle. From the upper part of the suspension, 1 ml was taken and added to 100 ml of Czapek-Dox broth medium with 0.1 g of unprocessed autoclaved sheep wool, 0.088 g/l of thiamine hydrochloride (Sigma-Aldrich, St. Louis, USA), and 0.096 g/l of glycine (Sigma-Aldrich, St. Louis, USA). The cultivation was performed in triplicate. Control samples were prepared using the same media without the inoculum and processed in the same manner. The strains were cultivated at 37°C with shaking at 200 rpm. The presence of glucose in the medium was measured daily using a glucose strip test (DIAPHAN®, Erba Lachema, Czech R.). The cultures were filtered two days after glucose depletion. The length of the cultivation period varied from 10 to 13 days for T. benhamiae var. luteum, 7 to 15 days for T. erinacei, and 7 to 9 days for T. mentagrophytes. Once the fungi had depleted the glucose, they were expected to primarily digest keratin and produce the highest number of secondary metabolites. Filtration was done using a double layer of sterile Miracloth filtration material (Merck, Rahway, USA) to remove large mycelium fragments. The remaining conidia and small hyphal elements were then removed using the Merck-Millex 0.22 µm syringe filter (Merck, Rahway, USA). The sterile filtrate was frozen at -35°C until analysis. The defrosted filtrate from the cultivations was pipetted into 4 ml portions in 20 ml headspace vials. A quality control sample was prepared by combining 1 ml from each sample, and the mixture was aliquoted into headspace vials with 5 ml in each and frozen. The quality control sample was used in each measurement to monitor any deviations in the analysis process.
GC-MS analysis
Sample analysis was performed using a gas chromatograph 450-GC and ion trap mass spectrometer 240-MS (GC-MS; both Varian, Walnut Creek, CA, USA). Prior to the analysis, the HS-SPME (Headspace-Solid Phase Microextraction) technique was performed using a CombiPal autosampler (CTC Analytics AG, MN, USA) equipped with a heating station, SPME fiber holder, and a 65 µm PDMS/DVB SPME fiber (Supelco, Bellefonte, PA, USA). The measuring conditions for the GC-MS analysis, extraction, and desorption steps were adopted from Semerád et al. [36] with modifications. A vial containing a strain cultivated for 30 days on sheep wool was transferred to the heating station and conditioned for 30 minutes at 40°C. Then, the SPME fiber was inserted into the headspace of the vial for 30 minutes at 40°C to extract the analytes. The injector and SPME liner (Topaz, Restek, Bellefonte, PA, USA) temperature was set to 250°C to desorb the analytes from the SPME fiber. After each injection, the SPME fiber was left in the injector for 10 minutes at 250°C to prevent sample cross-contamination. The separation of selected analytes was performed on a DB-5MS column (30 m × 0.25 mm id, 0.25 µm df; Agilent Technologies, Santa Clara, CA, USA). Helium was used as the carrier gas for GC in a constant flow of 1.0 ml/min. The temperature gradient was as follows: 40°C (2 minutes isothermal), 200°C (10°C/min), 260°C (25°C/min, 5 minutes isothermal). The total time for one GC-MS analysis was 25 minutes and 24 seconds. The mass spectrometer temperatures for the ion trap, manifold, and transfer line were set to 240, 50, and 280°C, respectively. The collected data ranged from 50–500 m/z. The software used for data acquisition was MS Workstation 6.9.1 (Agilent Technologies, Santa Clara, USA).
Data analysis (GC-MS analysis)
Data were processed using software AMDIS v.2.66 (NIST, Gaithersburg, USA) and Mass Profiler Professional v. 15.1 (Agilent Technologies, Santa Clara, CA, USA). The Mass Profiler Professional software generated a table containing 387 variables. Then, in short, redundant peaks, which were contamination from the column, SPME fiber, and medium, were removed from the data by comparing the sample data with controls. Peak areas were log2 transformed and normalized by probabilistic quotient normalization (PQN) [37]. Compound dentification was conducted by comparing the obtained spectra against the NIST20 library (RRID:SCR_006452), similarity values greater than 750 were tentatively identified, and values below 750 were considered unknowns (Unk). Data from the GC-MS analysis are provided in Supplementary 2. From the results of the GC-MS analysis, categorical data were created by classifying compounds in the headspace of dermatophytes as either present or absent in order to qualitatively investigate the VOCs spectra. The information for each species was averaged. The data were divided into three categories based on abundance level (present in all, present in some replicates, not present). Principal component analysis (PCA) and hierarchical clustering analysis (HCA) using Ward's method with Euclidean distances were performed normalized data. PCA analysis was conducted on the data from the GC-MS analysis and the data from the GC×GC-TOF MS analysis to compare the results of these methods.
Statistical analysis and data visualization were performed using R v. 4.2.1 [38]. For the statistical analysis, the packages stats, reshape2 [39], pvclust [40], FactoMineR [41], and factoextra [42] were utilized. For data visualization, the packages from Tidyverse [43], especially ggplot2 [44], viridis [45], gridExtra [46], and cowplot [47] were employed.
HS-SPME-GC×GC-TOF MS analysis
All extractions were performed using the triphasic DVB/CAR/PDMS SPME fiber (Supelco, Bellefonte, PA, USA). The samples underwent an initial incubation time of 10 minutes at a temperature of 40°C with an agitator speed of 250 rpm. After that, the compounds from the headspace were extracted for 20 minutes under the same temperature conditions. Subsequently, the VOCs were desorbed for 3 minutes at the injector port, which was set at 250°C prior to the GC-GC analyses. The GC-GC data acquisition was conducted on a Pegasus BT 4D GC-GC-TOF MS instrument equipped with a cryogenic modulator (LECO Corp., St Joseph, MI, USA). A mid-polar first-dimension column (1D) (Rxi-624Sil MS, 30 m x 0.25 mm id x 1.4 µm df; Restek Corp., Bellefonte, PA, USA) was combined with a highly polar second dimension (2D) column (Stabilwax, 2 m x 0.25 mm id x 0.5 µm df; Restek Corp., Bellefonte, PA, USA) for all chromatographic separations. The primary GC oven ramped from 40°C to 220°C with a 5°C/min gradient, resulting in a total run time of 36 minutes. The secondary oven was regulated with a temperature offset of + 20°C relative to the primary oven. The transfer line was maintained at an isothermal temperature of 250°C. A modulation period of 2.5 seconds was used for all separations, consisting of a 0.75-second hot jet and a 0.5-second cold jet. Mass spectra were acquired at a speed of 200 spectra/second between mass channels (m/z) 33 and 450. All measurements were performed in triplicate.
Data analysis (GC×GC-TOF MS)
The acquired data was processed using ChromaTOF Tile software (ver. 1.2.6.0, LECO Corp, RRID:SCR_023077), with a 1D × 2D tile size of 5 modulations (12.5 s) × 15 spectra (75 ms). The software generated a peak table containing 398 variables for the 45 chromatographic measurements. Analyte identification was conducted by comparing the obtained spectra against the NIST17 library (RRID:SCR_006452). Analytes with similarity values greater than 800 were tentatively identified, while analytes with similarity values below 800 were considered unknowns (Unk). Data from the analysis are provided in Supplementary 3. Prior to chemometric analysis, missing values, which accounted for 2.4% of the data, were imputed using the column median. The data was then scaled using the z-score method, which involved centring each column to a mean of 0 and scaling it by the standard deviation. Finally, the data was normalized by probabilistic quotient normalization (PQN) [37]. Statistical analysis and data visualization were performed using in-house MATLAB v. R2019b scripts (MathWorks, Natick, MA, USA, RRID: SCR_001622). Feature selection was performed using the machine-learning algorithm Random Forest [48].
Phylogenetic Analysis
Phylogenetic analysis was conducted using three loci (ITS, tubb, and tef1-α) and the strains that were used for GC-MS analysis. The strains can be found in Table 1, along with the accession numbers for DNA sequences. DNA was extracted from strains cultivated on Sabouraud broth (Sigma-Aldrich, St. Louis, USA) using the Quick-DNA Fungal/Bacterial miniprep Kit (Zymo Research, Irvine, USA). MyTaq polymerase (Meridian Bioscience, Wilford, USA) was used for the PCR. The primers for the three genetic loci were chosen according to Čmoková et al. [49]. The PCR products were sequenced using the Sanger method on ABI Prism 3130XL (Applied Biosystems, Waltham, USA). The sequences were aligned using the FFT-NS-i option implemented in the MAFFT online service [50]. The alignments were trimmed, concatenated, and analysed using Maximum likelihood (ML) and Bayesian inference (BI) methods. The best partition scheme and substitution models were determined through computations. Partitioning schemes and substitution models were selected using PartitionFinder 2 based on the Bayesian information criterion [51]. Introns, exons and segments of the ITS region were treated as independent datasets. The optimal partitioning schemes of the analysed sequences can be found in Supplementary 1. The ML analysis was performed using IQ-TREE version 2.1. 2 [52]. Nodal support was determined using nonparametric bootstrapping (BS) with 1000 replicates. The phylogenetic tree’s graphical outputs were generated using iTOL (Interactive Tree Of Life) [53].