Next Article in Journal
Characterization of the Barley Net Blotch Pathosystem at the Center of Origin of Host and Pathogen
Previous Article in Journal
Immunosuppressive Compounds Affect the Fungal Growth and Viability of Defined Aspergillus Species
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Pathogens
Volume 8
Issue 4
10.3390/pathogens8040274
pathogens-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Transcriptional Profile of Trichophyton rubrum Co-Cultured with Human Keratinocytes Shows New Insights about Gene Modulation by Terbinafine

by
Monise Fazolin Petrucelli
1,
Josie Budag Matsuda
1,
Kamila Peroni
2,
Pablo Rodrigo Sanches
3,
Wilson Araújo Silva, Jr.
2,3,4,5,
Rene Oliveira Beleboni
1,
Nilce Maria Martinez-Rossi
3,
Mozart Marins
1 and
Ana Lúcia Fachin
1,*
1
Biotechnology Unit, University of Ribeirão Preto-UNAERP, Av. Costábile Romano 2201, Ribeirão Preto 14960-900, SP, Brazil
2
National Institute of Science and Technology in Stem Cell and Cell Therapy, Center for Cell-Based Therapy, Ribeirão Preto 14051-140, SP, Brazil
3
Department of Genetics, Ribeirão Preto Medical School, University of São Paulo, Ribeirão Preto 14049-900, SP, Brazil
4
Center for Integrative System Biology-CISBi-NAP/USP, University of São Paulo, Ribeirão Preto 14049-900, SP, Brazil
5
Center for Medical Genomics, University Hospital of the Ribeirão Preto Medical School, University of São Paulo, Ribeirão Preto 14015-010, SP, Brazil
*
Author to whom correspondence should be addressed.
Pathogens 2019, 8(4), 274; https://doi.org/10.3390/pathogens8040274
Submission received: 6 November 2019 / Revised: 20 November 2019 / Accepted: 22 November 2019 / Published: 29 November 2019
Browse Figures
Review Reports Versions Notes

Abstract

:
The dermatophyte Trichophyton rubrum is the main causative agent of dermatophytoses worldwide. Although a superficial mycosis, its incidence has been increasing especially among diabetic and immunocompromised patients. Terbinafine is commonly used for the treatment of infections caused by dermatophytes. However, cases of resistance of T. rubrum to this allylamine were reported even with the efficacy of this drug. The present study is the first to evaluate the effect of terbinafine using a co-culture model of T. rubrum and human keratinocytes, mimicking a fungus-host interaction, in conjunction with RNA-seq technique. Our data showed the repression of several genes involved in the ergosterol biosynthesis cascade and the induction of genes encoding major facilitator superfamily (MFS)- and ATP-binding cassette superfamily (ABC)-type membrane transporter which may be involved in T. rubrum mechanisms of resistance to this drug. We observed that some genes reported in the scientific literature as candidates of new antifungal targets were also modulated. In addition, we found the modulation of several genes that are hypothetical in T. rubrum but that possess known orthologs in other dermatophytes. Taken together, the results indicate that terbinafine can act on various targets related to the physiology of T. rubrum other than its main target of ergosterol biosynthetic pathway.

1. Introduction

Dermatophytoses are superficial infections of keratinized tissues that are caused by a group of filamentous fungi, called dermatophytes [1]. Among these dermatophytes, Trichophyton rubrum is the main causative agent of dermatophytoses in the world [2,3]. Despite their restriction to the superficial layers of the epidermis, dermatophytoses have gained particular importance because of the increasing number of cases among diabetic patients [4] and patients using immunosuppressive drugs [5]. Furthermore, knowledge of the molecular mechanisms involved in the fungus-host interaction is still limited because of technical difficulties of the available models that mimic this interaction. However, new models such as co-culture of T. rubrum with keratinocyte cell lines and culture media containing keratin and elastin were introduced recently as strategies to evaluate compounds with antifungal activities [6,7,8].
Among commercial compounds, terbinafine, an antifungal agent of the allylamine class, is commonly used for the treatment of infections caused by dermatophytes, with the highest activity being observed against species of the genus Trichophyton, Microsporum, and Epidermophyton [9]. Its main target is the inhibition of the ergosterol biosynthetic pathway through a specific mechanism of action that consists of inhibiting the enzyme squalene epoxidase [10,11]. The inhibition of the early steps of this pathway triggers fungicidal activity against susceptible species, particularly filamentous fungi [11]. However, despite its efficacy, cases of resistance of T. rubrum to this allylamine have been frequently reported [12,13,14,15,16]. Studies suggest that the resistance mechanisms of species of the genus Trichophyton to terbinafine are associated with specific mutations in the squalene epoxidase gene [14,15] and with membrane transporters that mediate the efflux of the drug [17].
Previous studies showed the effect of terbinafine on the expression of T. rubrum genes using microarray [18] and subtractive suppressive cDNA library in culture medium [19]. However, DNA microarrays lack sensitivity for genes that are expressed at lower or very high levels and therefore have a smaller dynamic rang [20].
In order to present a more in-depth view of the effect of terbinafine on T. rubrum genes modulation, the present study is the first to evaluate the effect of this drug using large-scale analysis of the expression profile of T. rubrum co-cultured with human keratinocytes in an attempt to mimic fungus-host interactions. In addition to this model, we applied the dual RNA-seq technique that is widely used for the study of these interactions [21] in bacteria [22], viruses [23], and fungi [24,25,26], and enable the entire transcriptome to be surveyed in a high-throughput and quantitative manner [20].
Our data demonstrated the modulation of several T. rubrum genes involved in the biosynthesis of ergosterol in addition to the target gene of terbinafine that is already known, genes encoding major facilitator superfamily (MFS)- and ABC-type membrane transporters that are associated with the phenomenon of multidrug resistance, genes that have been discussed in the scientific literature as possible novel targets for antifungal compounds and genes which do not have a known function in T. rubrum yet, but have been modulated in response to the drug.

2. Results

2.1. Viability of HaCaT Keratinocytes after Co-Culture in the Presence of Terbinafine

The lactate dehydrogenase (LDH) assay was carried out to evaluate the viability of keratinocytes after 24 h of co-culture. The percentage of LDH release was 30% for co-culture treated or not with terbinafine, showing that terbinafine does not interfere in cell viability during co-cultivation. Triton X-100 (1%) was used as positive control and resulted in 100% of LDH release.

2.2. RNA-seq Analysis of T. rubrum Genes in Response to Co-Culture in the Presence of Terbinafine

Sequencing generated an average of 34 and 39 million raw reads for the co-culture and co-culture with terbinafine libraries, respectively. Low-quality reads were removed and the remaining reads were aligned to the T. rubrum reference genome (CBS 118892) of the Broad Institute’s Dermatophyte Comparative Database and Homo sapiens HG19 genome (UCSC Genome Bioinformatics Site, Santa Cruz, CA, USA). An average of 5% and 6.7% quality reads were obtained for the co-culture and co-culture with terbinafine libraries, respectively, when the sequences were compared with the T. rubrum genome. The reads generated for all conditions tested, as well as the percentages of alignment of each library, are shown in Supplementary Table S1.

2.3. Effect of Terbinafine on the Differential Expression of T. rubrum Genes Co-Cultured with Hacat Keratinocytes

Table 1; Table 2 show the 50 most induced and the 50 most repressed T. rubrum genes, respectively, after exposure to terbinafine for 24 h of co-culture. There were 277 T. rubrum genes that were differentially expressed after 24 h of co-culture in the presence of terbinafine. Of these genes, 163 were hypothetical or not annotated genes in T. rubrum, corresponding to approximately 59% of all modulated genes. Although considered hypothetical in T. rubrum, some orthologs were identified in other dermatophytes. It is important to note that some genes discussed in this paper are not listed in Table 1; Table 2 because the fold change values does not reach the 50 most induced or repressed genes. However, the complete list of all induced and repressed genes and their respective orthologs are shown in Supplementary Table S2; Table S3, respectively.

2.4. Functional Categorization of Differentially Expressed Genes

In an attempt to determine which molecular mechanisms were involved during co-culture of T. rubrum with HaCaT keratinocytes treated with terbinafine, the differentially expressed genes were categorized according to their molecular functions. The most representative categories are those with p < 0.05 and are shown in Figure 1.
The majority of the repressed genes were grouped into functional classes related to ribosomal functions, translation and oxidation-reduction processes. In contrast, most of the induced genes were grouped into functional categories related to ATP-binding, integral to membrane, and transmembrane transport.

2.5. Validation by qPCR

Fifteen genes were selected for validation: genes involved in ergosterol biosynthesis (TERG_01703; TERG_05717), MFS and ABC transporters (TERG_02822; TERG_08613), genes involved in T. rubrum metabolism and nutrient transport (TERG_08278; TERG_02609; TERG_06106; TERG_02542), and some potential antifungal targets that were also modulated in the presence of the drug (TERG_03223; TERG_08480; TERG_03861; TERG_08058; TERG_05698; TERG_02694).
The gene expression results obtained by RNA-seq showed a strong correlation (r = 0.93, p < 0.001) with the gene modulation values obtained by qPCR, demonstrating the reproducibility and accuracy of the sequencing. The comparison of the gene modulation values obtained by qPCR and RNA-seq is shown in Figure 2.

3. Discussion

Analysis of the transcriptional profile of T. rubrum co-cultured with HaCaT keratinocytes in the presence of terbinafine revealed a total of 277 T. rubrum genes that were differentially expressed after 24 h of culture. Within this set of genes, approximately 59% are hypothetical or not yet categorized in T. rubrum. However, approximately 48% of these hypothetical genes in T. rubrum possess orthologs with known functions in other dermatophytes. These findings contribute to further studies in order to identify the function of these hypothetical genes in the terbinafine response, thus enabling the discovery of targets which have not been described for this drug yet. Furthermore, the genomes of seven dermatophytes were sequenced and deposited, including T. rubrum, which will permit different types of analysis [27]. However, the assembly of the T. rubrum genome is not yet complete and may still be altered [28]. We therefore suggest that the expressive number of hypothetical genes that were found to be modulated in this study might be related to limitations in the annotation of the data currently available for the T. rubrum genome. As was observed in this work, a microarray analysis of gene expression in T. rubrum cultured in medium with keratin also demonstrated strong induction of genes encoding hypothetical proteins [29].
Based on the sequencing data generated, about 6.7% of the quality reads of the co-culture with terbinafine library could be aligned to the T. rubrum reference genome, indicating a predominance of human reads in this library. However, these reads reached coverage of 92.2% of the 8616 annotated T. rubrum genes, considering genes with at least one count read. It was expected that a restricted number of reads aligned with the T. rubrum genome in the co-culture library, since the difficulty in obtaining equivalent amounts of fungal and human RNA during sequencing of mixed transcriptomes has been discussed in the literature as one of the main challenges since the amount of RNA in the cell differs between the two cell types. Whereas a human cell contains about 20–25 pg RNA, a fungal cell contains only about 0.5–1 pg [21,30].
Regarding the molecular mechanisms involved in the co-culture of T. rubrum with keratinocytes in the presence of terbinafine, we observed that the functional categories that grouped the most of the modulated genes were those related to ribosomal functions, translation, oxidation-reduction process, integral component of membrane, membrane transport and ATP-binding. Within the set of genes which were grouped into these functional categories, we chose some of them for validation: serine/threonine (TERG_08278), sucrase ferredoxin (TERG_02609), sulfate permease 2 (TERG_06106), and an integral membrane protein (TERG_02542). In addition to these genes, we validated genes involved in the phenomenon of multidrug resistance, genes that participate in ergosterol biosynthesis, and genes suggested as target candidates for new antifungal drugs, which will be discussed in more detail below.

3.1. Transmembrane Transporters and Antifungal Resistance

A major problem in the treatment of fungal infections is the phenomenon of multidrug resistance, which is characterized by the acquisition of chemical and structural resistance to different compounds. Several biochemical mechanisms are involved in the process of multidrug resistance in fungi. The most common mechanisms are reduced absorption of the drug, structural alterations at the target site, and increased drug efflux [31]. The last mechanism is mediated mainly by the action of transmembrane transporters that belong to two superfamilies: the ATP-binding cassette superfamily (ABC) and the major facilitator superfamily (MFS) [32].
In the present study, we observed the induction of seven genes that encode transmembrane transporters belonging to the ABC and MFS superfamily: TERG_00574 (Log2 Fold Change: 1.83), TERG_02822 (Log2 Fold Change: 1.54), TERG_02912 (Log2 Fold Change: 1.89), TERG_05575 (Log2 Fold Change: 1.54), TERG_06679 (Log2 Fold Change: 1.95), TERG_08613 (Log2 Fold Change: 1.63), and TERG_08751 (Log2 Fold Change: 2.07). The modulation of these genes might be associated with the efflux of this allylamine and with resistance of T. rubrum to terbinafine.
The efflux of drugs seems to be the main resistance mechanism in dermatophytes [33]. Indeed, Fachin et al. [34] observed an increase in the expression of the TruMDR2 ABC-transporter gene in T. rubrum in the presence of terbinafine. Furthermore, Maranhão et al. [35] found that deletion of the TrMDR2 gene in the mutant strain ∆TruMDR2 reduced the infective capacity of T. rubrum, characterized by low growth of the fungus on human nails when compared to the growth of the wild-type strain. Taking together to our results, these datas suggests that transmembrane transporters play an important role in the pathogenicity of T. rubrum.

3.2. The Role of Terbinafine in Inhibition of Ergosterol Biosynthesis

The main target of terbinafine is the enzyme squalene epoxidase, which is encoded by the ERG1 gene, an early step in the late pathway of ergosterol biosynthesis. Thus, antifungal agents of the allylamine class target the ERG1 gene and inhibit the ergosterol biosynthesis, resulting in the accumulation of squalene, which is toxic, and consequently causing the death of the fungal cell [10,33,36].
As expected, we observed in the present study the repression of the ERG1 gene (TERG_05717/ Log2 Fold Change: −1.65) (Figure 2). However, our data suggest an important role of terbinafine in the repression of various other genes involved in the ergosterol biosynthesis cascade. In addition to the main target of terbinafine, other genes involved in the biosynthesis of ergosterol were also repressed: C-8 sterol isomerase (ERG2 gene) (TERG_06755/ Log2 Fold Change: −1.91), sterol C-24 reductase (ERG4 gene) (TERG_02979/ Log2 Fold Change: −2.22), and C-4 sterol methyl oxidase (ERG 25 gene) (TERG_08545/Log2 Fold Change: −2.13). Similar results werereported by [18] who evaluated the transcriptional profile of response to terbinafine in T. rubrum using a cDNA microarray technique. However, in that case the authors not observed modulation of the ERG1 gene.
We also highlight the repression of TERG_01703 (CYP51) (Log2 Fold Change: −2.13) and TERG_04906 (CYP61) (Log2 Fold Change: −1.57), which encode lanosterol 14α-demethylase (ERG11) and C-22 sterol desaturase (ERG5), respectively, considered as targets for azolic antifungal agents. The enzyme lanosterol 14α-demethylase is mainly studied in fungi because of its fundamental importance for membrane permeability [37]. In fungi, CYP51 catalyzes the removal of the C-14α-methyl group from lanosterol. This demethylation is considered a checkpoint in the transformation of lanosterol to other sterols, including ergosterol [38]. Since ergosterol is a fundamental component of fungal membranes, many drugs were designed to inhibit the enzyme 14α-demethylase and to prevent the enzyme from functioning [39], such as azole antifungal agents. Significant increases in the transcription of the ERG2, ERG3, ERG5, ERG6, ERG25 and ERG11 genes were described in many fungal species in response to treatment with azolic antifungal agents, including T. rubrum [40,41].
Considering these results, our data suggest an important role of terbinafine in the repression of various other genes involved in the ergosterol biosynthesis cascade (ERG2, ERG4, ERG5, ERG25, ERG11) besides its main target that is already known.

3.3. Genes Identified after Co-Culture of T. rubrum in the Presence of Terbinafine

In addition to the inhibition of genes involved in ergosterol biosynthesis and the induction of transmembrane transporter genes, we found other genes that were modulated during co-culture of T. rubrum in the presence of terbinafine that are important for the pathogenicity of this dermatophyte. Some of these genes were identified as possible targets to be explored by new antifungal drugs and will be discussed in more detail below.

3.3.1. Sulfite Efflux Pump SSU1 (TERG_02694)

In dermatophytes, keratinolytic proteases are activated by the cleavage of disulfide bridges through sulfite secretion (sulfitolysis), which is mediated by a sulfite efflux pump. It is important to facilitate the degradation of keratin into oligopeptides and amino acids that will be assimilated by the fungus. In addition, the sulfite efflux pump may also play a role in sulfite detoxification in this dermatophyte [17,42]. Considering its importance during the infection process, we highlight the repression of TERG_02694 (Log2 Fold Change: −1.74), which encodes the sulfite efflux pump SSU1 (Figure 2).
The sulfite efflux pump SSU1 is encoded by the TruSsu1 gene and was shown by Lechenne et al. [43] that the reduction of disulfide bridges in T. rubrum depends on the activity of a sulfite efflux pump encoded by the TruSsu1 gene, which belongs to the tellurite-resistance/dicarboxylate transporter (TDT) family. Thus, sulfite efflux pumps in keratinolytic fungi may represent new drug targets in dermatology, since they are involved in the activation of keratinolytic proteases that are important for fungus nutrition during infection.

3.3.2. Thioredoxin (TERG_08480)

In addition to the genes encoding sulfite efflux pumps, the thioredoxin gene (TERG_08480; Log2 Fold Change: 1.82) has also been proposed as a possible antifungal target [44] and was up-regulated by terbinafine in this work (Figure 2).
In response to fungal attacks, the host initiates a front line defense system, that it includes the production of reactive oxygen species (ROS). It has been observed that some pathogens neutralize the effects of oxidative stress through the synthesis of enzymes and antioxidant molecules such as superoxide dismutases, catalases, peroxidades, thioredoxin, among others [45].
Kaufman et al. [46] identified a fungal thioredoxin in Trichophyton mentagrophytes using differential cDNA screening technique to identify genes encoding putative virulence factors. The fungal thioredoxin was induced by superficial (keratin) and deep (elastin) skin elements suggesting that the product of this gene could be important in superficial and deep dermatophytosis. In Cryptococcus neoformans, the thioredoxin gene is believed to be involved in the mechanism of protection against oxidative stress produced by the host and is considered a virulence factor [44].
There is no evidence yet showing some relation between the thioredoxin gene and terbinafine or allylamine class antifungals. However, we highlight in this paper that thioredoxin gene was up-regulated in the presence of terbinafine when compared to co-culture in the absence of the drug, suggesting that T. rubrum probably induces this virulence factor necessary for growth in the host in response to the stress generated by the drug action.

3.3.3. Glycosyl Hydrolase (TERG_02742)

We also highlight the induction of TERG_02742 (Log2 Fold Change: 1.72), which encodes a glycosyl hydrolase (Figure 2). Glycosyl hydrolases are classified into more than 45 families [47]. Among these families, chitinases belong to family 18 of glycosyl hydrolases and were chemically and genetically validated in pathogenic fungi as a potential drug target. Chitinases hydrolyze chitin by cleavage of the β-1,4 glycosidic bond and play an important role in the biosynthesis of the hyphal cell wall of filamentous fungi. However, chitinases are also present in the human genome [48].
Further studies still need to be performed to verify which family belongs to glycosyl hydrolase modulated in this work in order to identify its possible function in T. rubrum as a new therapeutic target in response to drugs.

3.3.4. Transcription Factor C2H2 (TERG_03861)

The gene that encode the transcription factor C2H2 (TERG_03861; Log2 Fold Change: 2.39), also considered as a possible new virulence factor, was induced during co-culture in the presence of terbinafine (Figure 2). The role of this transcription factor in dermatophytes has not yet been elucidated. However, a query of the NCBI (National Center of Biotechnology) database showed that this protein exhibits 100% of similarity with the C2H2 transcription factor described in Aspergillus fumigatus. In A. fumigatus, the C2H2 transcription factor seems to function as a negative regulator of host cell damage and stimulation, as well as of virulence during invasive pulmonary disease [49]. Moreover, A. fumigatus C2H2 showed limited homology with the Candida albicans C2H2 Bcr1 transcription factor. In the pathogenic yeast C. albicans, the C2H2 Bcr1 transcription factor governs the expression of multiple cell surface proteins that mediate biofilm formation, as well as adherence to host cells [50,51]. The similarity of the C2H2 transcription factors of A. fumigatus with the C2H2 transcription factor of T. rubrum may suggest another virulence factor that is still unknown in this dermatophyte. However, studies still need to be done in order to verify if the transcription factor C2H2 modulated in this work is a virulence factor for T. rubrum and what is its function during the infectious process.

3.3.5. Βeta-lactamase and Metallo-β-lactamase (TERG_ 05698 and TERG_08360)

Other possible virulence factors that could be better studied in fungi are fungal lactamases. In this work, the genes encoding β-lactamase (Log2 Fold Change: −2.33) (Figure 2) and metallo-β-lactamase (Log2 Fold Change: −1.83) in T. rubrum were repressed by exposure to terbinafine
To combat antibiosis, bacteria developed resistance mechanisms such as efflux pumps and hydrolytic enzymes, including β-lactamases [52]. Parallel to the presence in bacteria, β-lactamase-encoding genes are also abundant in different fungal families. In contrast to bacteria, the function of these genes in fungi is still poorly understood. There are only two studies on the hydrolytic function of lactamase (metallo-β-lactamase, MBL)-encoding genes in Fusarium verticillioides and Fusarium pseudograminearum [53,54]. These studies raised the hypothesis that, asin bacteria, hydrolytic enzymes in fungi might be related to the degradation of and resistance to xenobiotic compounds [55].
Further studies should be carried out to evaluate the role of lactamases in dermatophytes, since the participation of fungal lactamases in dermatophytoses and in response to antifungal agents is still unknown.

3.3.6. N-acetylglucosamine-6-phosphate Deacetylase (TERG_03223)

N-acetylglucosamine (GlcNAc) is the main component of the structural polymers of bacteria, plants, and animals. In eukaryotes, GlcNAc is used as a substrate by chitin synthase, an enzyme that produces chitin, a GlcNAc homopolymer [56,57].
In the present study, we observed the induction of N-acetylglucosamine-6-phosphate deacetylase (TERG_03223; Log2 Fold Change: 1.96) (Figure 2). This enzyme catalyzes the deacetylation of GlcNAc. In bacteria, this is an important process for lipopolysaccharide synthesis and cell wall recycling [56]. The pathogenic fungus Candida albicans utilizes GlcNAc as a carbon source for growth. Deletion of the CaNAG2 gene, which encodes N-acetylglucosamine-phosphate deacetylase, attenuated C. albicans virulence in a mouse model of systemic infection [58].
We therefore suggest that the induction of N-acetylglucosamine-6-phosphate deacetylase in the co-culture with terbinafine may indicate a potential virulence factor of T. rubrum to maintain the integrity of its cell wall during the infectious process due to effects of this allylamine on the plasma membrane of the fungus.

4. Materials and Methods

4.1. Trichophyton rubrum, Media and Growth Conditions

The T. rubrum strain CBS 118892 (CBS-KNAW Fungal Biodiversity Center, Utrecht, The Netherlands) sequenced by the Broad Institute (Cambridge, MA, USA) was used. The strain was cultured on Sabouraud dextrose agar (Oxoid, Hampshire, UK) for 15 days at 28 ℃. After this period of incubation, the mycelial surface of the fungus was scraped with a sterile spatula to prepare conidia solution in saline (0.9% of NaCl). The conidia solution was filtered on glass wool and the solution concentration was adjusted in a hemocytometer to 1 × 107 conidia/mL.

4.2. Keratinocytes, Media and Growth Conditions

The immortal human keratinocyte cell line HaCaT purchased from Cell Lines Service GmbH (Eppelheim, Germany) was used. The cells were cultured on RPMI-1640 medium (Sigma Aldrich, St. Louis, MO, USA) supplemented with 10% fetal bovine serum at 37 ℃ in an atmosphere containing 5% CO2. Penicillin (100 U/mL) and streptomycin (100 µg/mL) were added to the culture medium to prevent bacterial contamination.

4.3. Chemicals

Terbinafine was purchased from Sigma Aldrich (St. Louis, MO, USA) and a stock solution (1 mg/mL) was prepared by dilution of the compound in ultrapure water.

4.4. Co-Culture Assay and Conditions

HaCaT cells (2.5 × 105 cells/mL) and T. rubrum conidia (1 × 107 conidia/mL) were used for the co-culture assay as described by Komoto et al. [6]. To evaluate the effect of terbinafine, 0.0162 µM (¼ of Minimal Inhibitory Concentration value) of this antifungal agent was added to the co-culture as determined by Komoto et al. [6]. As control, keratinocytes were cultured in the presence and absence of terbinafine to evaluate the effect of the antifungal agent on keratinocytes. The assay of each condition was carried out in triplicate in three independent experiments. The viability of HaCaT cells was evaluated after 24 h of co-culture treated or not with terbinafine by quantifying the release of lactate dehydrogenase (LDH) (TOX7 kit, Sigma Aldrich) according to manufacturer instructions and described by Santiago et al. [59]. Triton X-100 (Sigma Aldrich) was used as control. Absorbance was read at 490 nm in a microplate reader (Multiscan FC, Thermo Fisher Scientific, West Columbia, SC, USA).

4.5. RNA Isolation and Integrity Analysis

RNA was isolated and its integrity was analyzed as described by Petrucelli et al. [26]. Only RNA with an RNA integrity number (RIN) > 7 was accepted for construction of the cDNA libraries. The selected RNAs were quantified in a Quantus™ Fluorometer (Promega Corporation, Madison, WI, USA) and then used for the library construction.

4.6. Library Construction and Sequencing

The cDNA libraries were constructed in triplicate for each condition: keratinocyte culture, keratinocyte culture with terbinafine, co-culture and co-culture with terbinafine. We compared the library of co-culture with terbinafine vs. the library of co-culture in order to assess the differentially expressed genes between this conditions. We also compared the library of keratinocyte culture vs. the library of keratinocyte culture with terbinafine to verify if the drug had some effect in the modulation of human keratinocytes genes. However, there was no observed gene expression modulation (data not shown). The libraries were constructed using the TrueSeq® RNA Sample Preparation Kit v2 (Illumina, San Diego, CA, USA) according to manufacturer recommendations. The libraries were validated following the Library Quantitative PCR (qPCR) Quantification Guide (Illumina). A pool of 11 pM of each library was applied to each lane of the flowcell (Illumina) and clustering was performed using cBot (Illumina) according to manufacturer instructions.
Single-end and paired-end sequencing was performed on the Genome Analyzer IIx and Hiseq 2000 (Illumina), respectively, according to manufacturer recommendations. The RNA-seq data are deposited in the GEO database [60] under the accession number GSE135095.

4.7. Sequence Data Analysis

The sequence data were analyzed as described by Petrucelli et al. [26]. According to the distribution of the genes, those exhibiting statistical significance <0.05 and a log2 fold change ≥1.5 or ≤−1.5 were classified as differentially expressed genes (Supplementary Figure S1). Functional categorization of the T. rubrum genes was based on the Gene Ontology (GO) categories [61] using the Blast2GO algorithm. The BayGO algorithm was used for functional enrichment [62] and categories with a p value < 0.05 were considered the most representative.

4.8. Validation of the Data

Fifteen genes were validated by qPCR. For reverse transcription, 1 μg of total RNA from the same conditions assayed and used for sequencing was treated with DNAse 1 Amplification Grade® (Sigma Aldrich, St. Louis, MO, USA) to remove any contamination with genomic DNA. Messenger RNA was converted to cDNA using the High-Capacity cDNA Reverse Transcription® Kit (Applied Biosystems, Foster City, CA, USA) according to manufacturer recommendations. The real-time quantitative (RT)-PCR assays were carried out in triplicate using the SYBR Taq Ready Mix Kit (Sigma Aldrich) in an Mx3300 qPCR System (Stratagene, San Diego, CA, USA). The cycling conditions were: initial denaturation at 94 °C for 10 min, followed by 40 cycles at 94 °C for 2 min, at 60 °C for 60 s and at 72 °C for 1 min. A dissociation curve was constructed at the end of each PCR cycle to check the amplification products. The gene expression levels were calculated using the comparative 2−ΔΔCT method. The 18S [63] and β-tubulin [64] genes were used as normalizer genes. The results are reported as the mean ± standard deviation of three experiments. Pearson’s correlation test was used to evaluate the correlation between the qPCR and RNA-seq techniques. The primers used in the qPCR assays are shown in Supplementary Table S4.

5. Conclusions

The RNA sequencing of the T. rubrum co-culture in the presence of terbinafine showed that several other genes of the ergosterol biosynthesis cascade were repressed. We highlight the induction of transmembrane transporters, which may be associated with the efflux of this allylamine. In addition, other genes that could be involved in the pathogenicity of T. rubrum were also modulated in the presence of the drug. These findings offer new perspectives for the discovery of novel antifungal targets or even for structural modifications in the terbinafine molecule that would allow to increase its spectrum of antifungal activity. In addition, we observed the modulation of hypothetical genes with still unknown functions in T. rubrum but that possess orthologs in other dermatophytes.
However, further studies are necessary to confirm the relationship of these modulated genes with the pathogenicity and resistance of T. rubrum to terbinafine.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-0817/8/4/274/s1, Figure S1: Distribution of differentially expressed genes after 24 h of co-culture. The red points indicate differentially expressed genes, Table S1: General features of dual RNA-seq sequences against reference genomes. CO I, CO II, CO III: co-culture libraries; H I, H II, H III: human keratinocyte libraries. COT I, COT II, COT III: co-culture treated with terbinafine libraries. HT I, HT II, HT III: human keratinocyte treated with terbinafine. The libraries were constructed in triplicate, with I, II and III corresponding to the sample number of each condition. PE: paired-end sequence; SR: single read sequence, Table S2: Complete list of up-regulated genes in T. rubrum after 24 h of co-culture treated with terbinafine, Table S3: Complete list of down-regulated genes in T. rubrum after 24 h of co-culture treated with terbinafine, Table S4: Primers used for qPCR analysis.

Author Contributions

M.F.P. performed the laboratory experiments and wrote the manuscript. J.B.M. performed the bioinformatics analysis and qPCR experiments. K.P. and W.A.S.J. constructed the libraries and performed the sequencing. P.R.S. performed the bioinformatics analysis. R.O.B. discussed the manuscript and contributed with reagents and materials. N.M.M.-R. supervised the bioinformatics analysis. M.M. supervised the research and contributed reagents and materials. A.L.F. designed the project, supervised the research, contributed reagents/materials/analysis tools, and wrote the manuscript. All authors read and approved the final version of the manuscript.

Funding

This study was supported by grants from Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP, 2016/22701-9) and CAPES (PhD fellowships granted to M.F.P.).

Conflicts of Interest

The authors declare that they have no competing interests.

References

  1. Bouchara, J.P.; Mignon, B.; Chaturvedi, V. Dermatophytes and Dermatophytoses: A Thematic Overview of State of the Art, and the Directions for Future Research and Developments. Mycopathologia 2017, 182, 1–4. [Google Scholar] [CrossRef] [PubMed]
  2. Aly, R. Ecology and epidemiology of dermatophyte infections. J. Am. Acad. Dermatol. 1994, 31, S21–S25. [Google Scholar] [CrossRef]
  3. Havlickova, B.; Czaika, V.A.; Fredrich, M. Epidemiological trends in skin mycosis worldwide. Mycosis 2008, 51, 2–15. [Google Scholar] [CrossRef] [PubMed]
  4. Gupta, A.K.; Foley, K.A.; Versteeg, S.G. New Antifungal Agents and New Formulations Against Dermatophytes. Mycopathologia 2017, 182. [Google Scholar] [CrossRef]
  5. Kaur, R.; Panda, P.S.; Sardana, K.; Khan, S. Mycological Pattern of Dermatomycoses in a Tertiary Care Hospital. J. Trop. Med. 2015, 2015. [Google Scholar] [CrossRef]
  6. Komoto, T.T.; Bitencourt, T.A.; Silva, G.; Beleboni, R.O.; Marins, M.; Fachin, A.L. Gene Expression Response of Trichophyton rubrum during Coculture on Keratinocytes Exposed to Antifungal Agents. Evid. Based Complement. Altern. Med. 2015, 2015. [Google Scholar] [CrossRef]
  7. Cantelli, B.A.M.; Bitencourt, T.A.; Komoto, T.T.; Beleboni, R.O.; Marins, M.; Fachin, A.L. Caffeic acid and licochalcone A interfere with the glyoxylate cycle of Trichophyton rubrum. Biomed. Pharmacother. 2017, 96, 1389–1394. [Google Scholar] [CrossRef]
  8. Bitencourt, T.A.; Macedo, C.; Franco, M.E.; Rocha, M.C.; Moreli, I.S.; Cantelli, B.A.M.; Sanches, P.R.; Beleboni, R.O.; Malavazi, I.; Passos, G.A.; et al. Trans-chalcone activity against Trichophyton rubrum relies on an interplay between signaling pathways related to cell wall integrity and fatty acid metabolism. BMC Genom. 2019, 20, 411. [Google Scholar] [CrossRef]
  9. Abdel-Rahman, S. Newland Update on terbinafine with a focus on dermatophytoses. Clin. Cosmet. Investig. Dermatol. 2009, 49. [Google Scholar] [CrossRef]
  10. Ryder, N.S. Terbinafine: Mode of action and properties of the squalene epoxidase inhibition. Br. J. Dermatol. 1992, 126, 2–7. [Google Scholar] [CrossRef]
  11. Odds, F.C.; Brown, A.J.P.; Gow, N.A.R. Antifungal agents: Mechanisms of action. Trends Microbiol. 2003, 11, 272–279. [Google Scholar] [CrossRef]
  12. Mukherjee, P.K.; Leidich, S.D.; Isham, N.; Leitner, I.; Ryder, N.S.; Ghannoum, M.A. Clinical Trichophyton rubrum strain exhibiting primary resistance to terbinafine. Antimicrob. Agents Chemother. 2003, 47, 82–86. [Google Scholar] [CrossRef] [PubMed]
  13. Osborne, C.S.; Leitner, I.; Favre, B.; Neil, S.; Osborne, C.S.; Leitner, I.; Favre, B.; Ryder, N.S. Amino Acid Substitution in Trichophyton rubrum Squalene Epoxidase Associated with Resistance to Terbinafine Amino Acid Substitution in Trichophyton rubrum Squalene Epoxidase Associated with Resistance to Terbinafine. Antimicrob. Agents Chemother. 2005, 49, 2840–2844. [Google Scholar] [CrossRef] [PubMed]
  14. Yamada, T.; Maeda, M.; Alshahni, M.M.; Tanaka, R.; Yaguchi, T.; Bontems, O.; Salamin, K.; Fratti, M.; Monod, M. Terbinafine Resistance of Trichophyton Clinical Isolates Caused by Specific Point Mutations in the Squalene Epoxidase Gene. Antimicrob. Agents Chemother. 2017, 61. [Google Scholar] [CrossRef] [PubMed]
  15. Rudramurthy, S.M.; Shankarnarayan, S.A.; Dogra, S.; Shaw, D.; Mushtaq, K.; Paul, R.A.; Narang, T.; Chakrabarti, A. Mutation in the Squalene Epoxidase Gene of Trichophyton interdigitale and Trichophyton rubrum Associated with Allylamine Resistance. Antimicrob. Agents Chemother. 2018, 62. [Google Scholar] [CrossRef] [PubMed]
  16. Schøsler, L.; Andersen, L.K.; Arendrup, M.C.; Sommerlund, M. Recurrent terbinafine resistant Trichophyton rubrum infection in a child with congenital ichthyosis. Pediatr. Dermatol. 2018, 35, 259–260. [Google Scholar] [CrossRef]
  17. Peres, N.T.D.A.; Maranhão, F.C.A.; Rossi, A.; Martinez-Rossi, N.M. Dermatophytes: Host-pathogen interaction and antifungal resistance. An. Bras. Dermatol. 2010, 85, 657–667. [Google Scholar] [CrossRef]
  18. Zhang, W.; Yu, L.; Yang, J.; Wang, L.; Peng, J.; Jin, Q. Transcriptional profiles of response to terbinafine in Trichophyton rubrum. Appl. Microbiol. Biotechnol. 2009, 82, 1123–1130. [Google Scholar] [CrossRef]
  19. Maranhão, F.C.A.; Paião, F.G.; Martinez-Rossi, N.M. Isolation of transcripts over-expressed in human pathogen Trichophyton rubrum during growth in keratin. Microb. Pathog. 2007, 43, 166–172. [Google Scholar] [CrossRef]
  20. Wang, Z.; Gerstein, M.; Snyder, M. RNA-Seq: A revolutionary tool for transcriptomics. Nat. Rev. Genet. 2009, 10, 57–63. [Google Scholar] [CrossRef]
  21. Wolf, T.; Kämmer, P.; Brunke, S.; Linde, J. Two’s company: Studying interspecies relationships with dual RNA-seq. Curr. Opin. Microbiol. 2018, 42, 7–12. [Google Scholar] [CrossRef] [PubMed]
  22. Aprianto, R.; Slager, J.; Holsappel, S.; Veening, J.W. Time-resolved dual RNA-seq reveals extensive rewiring of lung epithelial and pneumococcal transcriptomes during early infection. Genome Biol. 2016, 17, 198. [Google Scholar] [CrossRef] [PubMed]
  23. Wesolowska-Andersen, A.; Everman, J.L.; Davidson, R.; Rios, C.; Herrin, R.; Eng, C.; Janssen, W.J.; Liu, A.H.; Oh, S.S.; Kumar, R.; et al. Dual RNA-seq reveals viral infections in asthmatic children without respiratory illness which are associated with changes in the airway transcriptome. Genome Biol. 2017, 18, 12. [Google Scholar] [CrossRef]
  24. Tierney, L.; Linde, J.; Müller, S.; Brunke, S.; Molina, J.C.; Hube, B.; Schöck, U.; Guthke, R.; Kuchler, K. An interspecies regulatory network inferred from simultaneous RNA-seq of Candida albicans invading innate immune cells. Front. Microbiol. 2012, 3, 85. [Google Scholar] [CrossRef] [PubMed]
  25. Meyer, F.E.; Shuey, L.S.; Naidoo, S.; Mamni, T.; Berger, D.K.; Myburg, A.A.; van den Berg, N.; Naidoo, S. Dual RNA-Sequencing of Eucalyptus nitens during Phytophthora cinnamomi Challenge Reveals Pathogen and Host Factors Influencing Compatibility. Front. Plant Sci. 2016, 7, 191. [Google Scholar] [CrossRef]
  26. Petrucelli, M.F.; Peronni, K.; Sanches, P.R.; Komoto, T.T.; Matsuda, J.B.; da Silva Junior, W.A.; Beleboni, R.O.; Martinez-Rossi, N.M.; Marins, M.; Fachin, A.L. Dual RNA-Seq analysis of trichophyton rubrum and HaCat keratinocyte co-culture highlights important genes for fungal-host interaction. Genes 2018, 9, 362. [Google Scholar] [CrossRef]
  27. Martinez, D.A.; Oliver, B.G.; Gräser, Y.; Goldberg, J.M.; Li, W.; Martinez-Rossi, N.M.; Monod, M.; Shelest, E.; Barton, R.C.; Birch, E.; et al. Comparative genome analysis of Trichophyton rubrum and related dermatophytes reveals candidate genes involved in infection. MBio 2012, 3, e00259-12. [Google Scholar] [CrossRef]
  28. Franco, M.E.; Bitencourt, T.A.; Marins, M.; Fachin, A.L. In silico characterization of tandem repeats in Trichophyton rubrum and related dermatophytes provides new insights into their role in pathogenesis. Database 2017, 2017, 1–10. [Google Scholar] [CrossRef]
  29. Bitencourt, T.A.; Macedo, C.; Franco, M.E.; Assis, A.F.; Komoto, T.T.; Stehling, E.G.; Beleboni, R.O.; Malavazi, I.; Marins, M.; Fachin, A.L. Transcription profile of Trichophyton rubrum conidia grown on keratin reveals the induction of an adhesin-like protein gene with a tandem repeat pattern. BMC Genom. 2016, 17, 249. [Google Scholar] [CrossRef]
  30. Westermann, A.J.; Barquist, L.; Vogel, J. Resolving host–pathogen interactions by dual RNA-seq. PLoS Pathog. 2017, 13, e1006033. [Google Scholar] [CrossRef]
  31. Martinez-Rossi, N.M.; Peres, N.T.A.; Rossi, A. Antifungal resistance mechanisms in dermatophytes. Mycopathologia 2008, 166, 369–383. [Google Scholar] [CrossRef] [PubMed]
  32. Costa, C.; Dias, P.J.; Sá-correia, I.; Teixeira, M.C. MFS multidrug transporters in pathogenic fungi: Do they have real clinical impact? Front. Physiol. 2014, 5, 197. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Scorzoni, L.; de Paula e Silva, A.C.A.; Marcos, C.M.; Assato, P.A.; de Melo, W.C.M.A.; de Oliveira, H.C.; Costa-Orlandi, C.B.; Mendes-Giannini, M.J.S.; Fusco-Almeida, A.M. Antifungal therapy: New advances in the understanding and treatment of mycosis. Front. Microbiol. 2017, 8, 36. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Fachin, A.L.; Ferreira-Nozawa, M.S.; Maccheroni, W.; Martinez-Rossi, N.M. Role of the ABC transporter TruMDR2 in terbinafine, 4-nitroquinoline N-oxide and ethidium bromide susceptibility in Trichophyton rubrum. J. Med. Microbiol. 2006, 55, 1093–1099. [Google Scholar] [CrossRef] [Green Version]
  35. Maranhao, F.C.A.; Paiao, F.G.; Fachin, A.L.; Martinez-Rossi, N.M. Membrane transporter proteins are involved in Trichophyton rubrum pathogenesis. J. Med. Microbiol. 2009, 58, 163–168. [Google Scholar] [CrossRef] [Green Version]
  36. Kumar Nigam, P. Antifungal drugs and resistance: Current concepts. Dermatol. Online 2015, 6, 212–221. [Google Scholar] [CrossRef]
  37. Daum, G.; Lees, N.D.; Bard, M.; Dickson, R. Biochemistry, cell biology and molecular biology of lipids ofSaccharomyces cerevisiae. Yeast 1998, 14, 1471–1510. [Google Scholar] [CrossRef]
  38. Lepesheva, G.I.; Waterman, M.R. Sterol 14α-demethylase cytochrome P450 (CYP51), a P450 in all biological kingdoms. Biochim. Biophys. Acta Gen. Subj. 2007, 1770, 467–477. [Google Scholar] [CrossRef] [Green Version]
  39. Becher, R.; Wirsel, S.G.R. Fungal cytochrome P450 sterol 14α-demethylase (CYP51) and azole resistance in plant and human pathogens. Appl. Microbiol. Biotechnol. 2012, 95, 825–840. [Google Scholar] [CrossRef]
  40. Yu, L.; Zhang, W.; Wang, L.; Yang, J.; Liu, T.; Peng, J.; Leng, W.; Chen, L.; Li, R.; Jin, Q. Transcriptional profiles of the response to ketoconazole and amphotericin B in Trichophyton rubrum. Antimicrob. Agents Chemother. 2007, 51, 144–153. [Google Scholar] [CrossRef] [Green Version]
  41. Sun, X.; Wang, W.; Wang, K.; Yu, X.; Liu, J.; Zhou, F.; Xie, B.; Li, S. Sterol C-22 Desaturase ERG5 Mediates the Sensitivity to Antifungal Azoles in Neurospora crassa and Fusarium verticillioides. Front. Microbiol. 2013, 4, 127. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Chinnapun, D. Virulence factors involved in pathogenicity of dermatophytes. Walailak J. Sci. Technol. 2015, 12, 573–580. [Google Scholar]
  43. Lechenne, B.; Reichard, U.; Zaugg, C.; Fratti, M.; Kunert, J.; Boulat, O.; Monod, M. Sulphite efflux pumps in Aspergillus fumigatus and dermatophytes. Microbiology 2007, 153, 905–913. [Google Scholar] [CrossRef] [Green Version]
  44. Missall, T.A.; Lodge, J.K. Function of the thioredoxin proteins in Cryptococcus neoformans during stress or virulence and regulation by putative transcriptional modulators. Mol. Microbiol. 2005, 57, 847–858. [Google Scholar] [CrossRef] [PubMed]
  45. Ahmad, A.; Molepo, J.; Patel, M. Challenges in the Development of Antifungal Agents Against Candida: Scope of Phytochemical Research. Curr. Pharm. Des. 2016, 22, 4135–4150. [Google Scholar] [CrossRef]
  46. Kaufman, G.; Berdicevsky, I.; Woodfolk, J.A.; Horwitz, B.A. Markers for Host-Induced Gene Expression in Trichophyton Dermatophytosis. Infect. Immun. 2005, 73, 6584–6590. [Google Scholar] [CrossRef] [Green Version]
  47. Davies, G.; Henrissat, B. Structures and mechanisms of glycosyl hydrolases. Structure 1995, 3, 853–859. [Google Scholar] [CrossRef] [Green Version]
  48. Castelli, M.V.; Butassi, E.; Monteiro, M.C.; Svetaz, L.A.; Vicente, F.; Zacchino, S.A. Novel antifungal agents: A patent review (2011—Present). Expert Opin. Ther. Pat. 2014, 24, 323–338. [Google Scholar] [CrossRef]
  49. Ejzykowicz, D.E.; Solis, N.V.; Gravelat, F.N.; Chabot, J.; Li, X.; Sheppard, D.C.; Filler, S.G. Role of Aspergillus fumigatus DvrA in Host Cell Interactions and Virulence. Eukaryot. Cell 2010, 9, 1432–1440. [Google Scholar] [CrossRef] [Green Version]
  50. Nobile, C.J.; Mitchell, A.P. Regulation of Cell-Surface Genes and Biofilm Formation by the C. albicans Transcription Factor Bcr1p. Curr. Biol. 2005, 15, 1150–1155. [Google Scholar] [CrossRef] [Green Version]
  51. Nobile, C.J.; Andes, D.R.; Nett, J.E.; Smith, F.J.; Yue, F.; Phan, Q.T.; Edwards, J.E.; Filler, S.G.; Mitchell, A.P. Critical Role of Bcr1-Dependent Adhesins in C. albicans Biofilm Formation In Vitro and In Vivo. PLoS Pathog. 2006, 2, e63. [Google Scholar] [CrossRef] [PubMed]
  52. Dallenne, C.; Da Costa, A.; Decré, D.; Favier, C.; Arlet, G. Development of a set of multiplex PCR assays for the detection of genes encoding important β-lactamases in Enterobacteriaceae. J. Antimicrob. Chemother. 2010, 65, 490–495. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Kettle, A.J.; Carere, J.; Batley, J.; Benfield, A.H.; Manners, J.M.; Kazan, K.; Gardiner, D.M. A γ-lactamase from cereal infecting Fusarium spp. catalyses the first step in the degradation of the benzoxazolinone class of phytoalexins. Fungal Genet. Biol. 2015, 83, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Glenn, A.E.; Davis, C.B.; Gao, M.; Gold, S.E.; Mitchell, T.R.; Proctor, R.H.; Stewart, J.E.; Snook, M.E. Two Horizontally Transferred Xenobiotic Resistance Gene Clusters Associated with Detoxification of Benzoxazolinones by Fusarium Species. PLoS ONE 2016, 11, e0147486. [Google Scholar] [CrossRef] [PubMed]
  55. Gao, M.; Glenn, A.E.; Blacutt, A.A.; Gold, S.E. Fungal Lactamases: Their Occurrence and Function. Front. Microbiol. 2017, 8, 1775. [Google Scholar] [CrossRef] [PubMed]
  56. Yadav, V.; Panilaitis, B.; Shi, H.; Numuta, K.; Lee, K.; Kaplan, D.L. N-acetylglucosamine 6-Phosphate Deacetylase (nagA) Is Required for N-acetyl Glucosamine Assimilation in Gluconacetobacter xylinus. PLoS ONE 2011, 6, e18099. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Liu, N.; Tu, J.; Dong, G.; Wang, Y.; Sheng, C. Emerging New Targets for the Treatment of Resistant Fungal Infections. J. Med. Chem. 2018. [Google Scholar] [CrossRef]
  58. Yamada-Okabe, T.; Sakamori, Y.; Mio, T.; Yamada-Okabe, H. Identification and characterization of the genes for N-acetylglucosamine kinase and N-acetylglucosamine-phosphate deacetylase in the pathogenic fungus Candida albicans. Eur. J. Biochem. 2001, 268, 2498–2505. [Google Scholar] [CrossRef]
  59. Santiago, K.; Bomfim, G.F.; Criado, P.R.; Almeida, S.R. Monocyte-derived dendritic cells from patients with dermatophytosis restrict the growth of Trichophyton rubrum and induce CD4-T cell activation. PLoS ONE 2014, 9, e110879. [Google Scholar] [CrossRef]
  60. Edgar, R. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30, 207–210. [Google Scholar] [CrossRef] [Green Version]
  61. Blake, J.A.; Harris, M.A. The Gene Ontology (GO) Project: Structured Vocabularies for Molecular Biology and Their Application to Genome and Expression Analysis. In Current Protocols in Bioinformatics; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2008; pp. 1–9. ISBN 0471250953. [Google Scholar]
  62. Vêncio, R.Z.; Koide, T.; Gomes, S.L.; de B Pereira, C.A. BayGO: Bayesian analysis of ontology term enrichment in microarray data. BMC Bioinform. 2006, 7, 86. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Bitencourt, T.A.; Komoto, T.T.; Massaroto, B.G.; Miranda, C.E.S.; Beleboni, R.O.; Marins, M.; Fachin, A.L. Trans-chalcone and quercetin down-regulate fatty acid synthase gene expression and reduce ergosterol content in the human pathogenic dermatophyte Trichophyton rubrum. BMC Complement. Altern. Med. 2013, 13, 229. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Jacob, T.R.; Peres, N.T.A.; Persinoti, G.F.; Silva, L.G.; Mazucato, M.; Rossi, A.; Martinez-Rossi, N.M. rpb2 is a reliable reference gene for quantitative gene expression analysis in the dermatophyte Trichophyton rubrum. Med. Mycol. 2012, 50, 368–377. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Pathogens 08 00274 g001 550
Figure 1. Gene Ontology-based functional categorization of differentially expressed genes. Main representative functional categories (p < 0.05) of differentially expressed genes in T. rubrum co-cultured with keratinocytes treated with terbinafine.
Figure 1. Gene Ontology-based functional categorization of differentially expressed genes. Main representative functional categories (p < 0.05) of differentially expressed genes in T. rubrum co-cultured with keratinocytes treated with terbinafine.
Pathogens 08 00274 g001
Pathogens 08 00274 g002 550
Figure 2. Comparison of gene modulation obtained by RNA-seq and qPCR. The error bars represent the standard error of three independent replicates. Pearson’s test indicated very strong correlation between the two techniques (r = 0.93; p < 0.001).
Figure 2. Comparison of gene modulation obtained by RNA-seq and qPCR. The error bars represent the standard error of three independent replicates. Pearson’s test indicated very strong correlation between the two techniques (r = 0.93; p < 0.001).
Pathogens 08 00274 g002
Table
Table 1. The 50 most up-regulated genes in 24 h of co-culture treated with terbinafine.
Table 1. The 50 most up-regulated genes in 24 h of co-culture treated with terbinafine.
IDLog2 Fold ChangeGene Product NameOrthologous
TERG_005235.50Hypothetical protein-
TERG_081824.79Hypothetical protein-
TERG_016364.72Hypothetical proteinMicrosporum gypseum CBS 118893 ADP-ribosylglycohydrolase (1715 nt)
TERG_029024.64Hypothetical protein-
TERG_005404.44Hypothetical proteinTrichophyton verrucosum HKI 0517 oxidoreductase, zinc-binding, putative (1134 nt)
TERG_120294.19Hypothetical protein-
TERG_020674.03Hypothetical protein-
TERG_026093.93Sucrase/ferredoxin Trichophyton equinum CBS 127.97 actin patches distal protein 1 (950 nt)
TERG_063473.93Hypothetical protein-
TERG_117473.79Hypothetical protein-
TERG_085033.45Hypothetical proteinMicrosporum canis CBS 113480 secalin (2547 nt)
TERG_081213.38Protein kinaseMicrosporum gypseum CBS 118893 protein kinase subdomain-containing protein (891 nt)
TERG_082783.33Serine/threonine Trichophyton tonsurans CBS 112818 serine/threonine protein kinase (2270 nt)
TERG_058433.32Hypothetical proteinTrichophyton equinum CBS 127.97 F-box domain-containing protein (2405 nt)
TERG_081943.24Hypothetical proteinMicrosporum canis CBS 113480 serine/threonine protein kinase (2427 nt)
TERG_001973.23Aldose 1-epimeraseTrichophyton verrucosum HKI 0517 aldose 1-epimerase family protein, putative (1631 nt)
TERG_028993.17Hypothetical protein-
TERG_051113.12Hypothetical protein-
TERG_032933.03Hypothetical protein-
TERG_031322.91Hypothetical protein-
TERG_009592.79Hypothetical proteinArthroderma benhamiae CBS 112371 RNA binding protein, putative (3143 nt)
TERG_032522.78Hypothetical protein-
TERG_033042.77Hypothetical proteinTrichophyton verrucosum HKI 0517 AAA family ATPase, putative (2368 nt)
TERG_064022.75Hypothetical proteinTrichophyton verrucosum HKI 0517 Ser/Thr protein phosphatase family protein (921 nt)
TERG_024482.74Hypothetical protein-
TERG_119322.74Hypothetical protein-
TERG_023032.74 Ankyrin repeat proteinArthroderma benhamiae CBS 112371 ankyrin repeat protein (5579 nt)
TERG_049512.73Hypothetical proteinTrichophyton equinum CBS 127.97 U-box domain-containing protein (2461 nt)
TERG_064452.72Hypothetical protein-
TERG_069922.69Pyridine nucleotide-disulfide oxidoreductaseTrichophyton tonsurans CBS 112818 pyridine nucleotide-disulphide oxidoreductase (1785 nt)
TERG_029002.69Hypothetical protein-
TERG_042342.66Hypothetical proteinTrichophyton verrucosum HKI 0517 hydrophobin, putative (562 nt)
TERG_005832.63Hypothetical protein-
TERG_080462.56Hypothetical proteinMicrosporum gypseum CBS 118893 beta-lactamase (852 nt)
TERG_047212.54Glutamate carboxypeptidaseTrichophyton equinum CBS 127.97 glutamate carboxypeptidase (2412 nt)
TERG_059092.51Hypothetical protein-
TERG_052392.46DNA polymerase lambdaTrichophyton verrucosum HKI 0517 DNA polymerase POL4, putative (2133 nt)
TERG_019562.46Hypothetical proteinArthroderma benhamiae CBS 112371 C2H2 finger domain protein, putative (3359 nt)
TERG_060652.43Hypothetical proteinTrichophyton verrucosum HKI 0517 conserved glycine-rich protein (903 nt)
TERG_062072.39Hypothetical proteinTrichophyton verrucosum HKI 0517 proline oxidase PrnD (1941 nt)
TERG_070342.39Hypothetical protein-
TERG_038612.39C2H2 transcription factorTrichophyton tonsurans CBS 112818 C2H2 transcription factor (895 nt)
TERG_034432.34Hypothetical proteinTrichophyton equinum CBS 127.97 ankyrin repeat protein (1561 nt)
TERG_036282.33Serine/threonine protein kinaseTrichophyton tonsurans CBS 112818 serine/threonine protein kinase (2187 nt)
TERG_006422.32Hypothetical proteinTrichophyton equinum CBS 127.97 HHE domain-containing protein (552 nt)
TERG_079822.32Hypothetical protein-
TERG_054692.32Hypothetical protein-
TERG_123292.32Hypothetical proteinTrichophyton tonsurans CBS 112818 Ku70/Ku80 beta-barrel domain-containing protein (2242 nt)
TERG_069902.30Hypothetical protein-
TERG_021312.29Hypothetical protein-
Table
Table 2. The 50 most down-regulated genes in 24 h of co-culture treated with terbinafine.
Table 2. The 50 most down-regulated genes in 24 h of co-culture treated with terbinafine.
IDLog2 Fold ChangeGene Product NameOrthologous
TERG_01731−3.68 Hypothetical protein -
TERG_11886−3.49 Hypothetical protein Trichophyton tonsurans CBS 112818 copper radical oxidase (2955 nt)
TERG_06315−3.32 Hypothetical protein Arthroderma benhamiae CBS 112371 integral membrane protein (1704 nt)
TERG_00499−3.32 Hypothetical protein -
TERG_02811−3.26 Hypothetical protein Arthroderma benhamiae CBS 112371 acetyl xylan esterase (Axe1), putative (955 nt)
TERG_01900−3.25Aquaglyceroporin Arthroderma benhamiae CBS 112371 aquaglyceroporin, putative (1058 nt)
TERG_00520−3.19 Hypothetical protein -
TERG_03105−3.16 Hypothetical protein -
TERG_07234−2.76 Hypothetical protein -
TERG_04164−2.60 Hypothetical protein -
TERG_07011−2.58 Hypothetical protein Arthroderma benhamiae CBS 112371 conserved fungal protein (800 nt)
TERG_12474−2.48 Hypothetical protein Trichophyton tonsurans CBS 112818 ABC-transporter (2471 nt)
TERG_06276−2.47Chromate ion transporterTrichophyton tonsurans CBS 112818 chromate ion transporter (1906 nt)
TERG_01901−2.44Glycerol kinaseTrichophyton equinum CBS 127.97 glycerol kinase (1680 nt)
TERG_00765−2.38 Hypothetical protein -
TERG_01406−2.37 Hypothetical protein Trichophyton equinum CBS 127.97 phospholipase D (939 nt)
TERG_12475−2.37 Hypothetical protein Arthroderma benhamiae CBS 112371 ABC-transporter, putative (2336 nt)
TERG_07199−2.35 Hypothetical protein -
TERG_05698−2.33Beta-lactamaseTrichophyton verrucosum HKI 0517 transesterase (LovD), putative (1428 nt)
TERG_07539−2.30 Hypothetical protein Trichophyton tonsurans CBS 112818 multidrug resistance protein (2253 nt)
TERG_11621−2.27 Hypothetical protein -
TERG_02161−2.25DOC FamilyMicrosporum canis CBS 113480 DOC family protein (456 nt)
TERG_02979−2.22Delta(24(24(1)))-sterol reductaseTrichophyton tonsurans CBS 112818 Delta(24(24(1)))-sterol reductase (1551 nt)
TERG_07798−2.22 Hypothetical protein -
TERG_02722−2.21 Hypothetical protein Trichophyton equinum CBS 127.97 WSC domain containing protein (2371 nt)
TERG_00754−2.20 Hypothetical protein -
TERG_04382−2.15C-14 sterol reductaseTrichophyton tonsurans CBS 112818 c-14 sterol reductase (1597 nt)
TERG_05808−2.14 Hypothetical protein -
TERG_03083−2.143-dehydroquinate synthaseTrichophyton equinum CBS 127.97 pentafunctional AROM polypeptide (4908 nt)
TERG_01703−2.13Cytochrome P450 51Trichophyton equinum CBS 127.97 cytochrome P450 51 (1790 nt)
TERG_08666−2.13 Hypothetical protein -
TERG_08545−2.13C-4 methylsterol oxidase Trichophyton equinum CBS 127.97 C-4 methyl sterol oxidase Erg25 (940 nt)
TERG_01676−2.106,7-dimethyl-8-ribityllumazine synthaseTrichophyton tonsurans CBS 112818 6,7-dimethyl-8-ribityllumazine synthase (761 nt)
TERG_04041−2.10Sad1/UNC domain-containing proteinTrichophyton verrucosum HKI 0517 Sad1/UNC domain protein (2625 nt)
TERG_07810−2.09 Hypothetical protein Trichophyton tonsurans CBS 112818 phospholipase (3533 nt)
TERG_00613−2.09 Hypothetical protein -
TERG_04793−2.04CyclinTrichophyton tonsurans CBS 112818 cyclin (1194 nt)
TERG_06265−2.03 Hypothetical protein Trichophyton equinum CBS 127.97 LPS glycosyltransferase (1368 nt)
TERG_08359−2.01FAD-dependent monooxygenaseTrichophyton tonsurans CBS 112818 FAD-dependent monooxygenase (1321 nt)
TERG_02842−2.01 Hypothetical protein Trichophyton equinum CBS 127.97 6-hydroxy-D-nicotine oxidase (1598 nt)
TERG_03204−2.0060S ribosomal protein L7Trichophyton tonsurans CBS 112818 60S ribosomal protein L7 (1079 nt)
TERG_07797−1.98Isoflavone reductase Trichophyton equinum CBS 127.97 amino acid permease (2474 nt)
TERG_01994−1.98 Hypothetical protein Trichophyton equinum CBS 127.97 OPT oligopeptide transporter protein (3035 nt)
TERG_00032−1.98Mitochondrial dicarboxylate carrierTrichophyton tonsurans CBS 112818 mitochondrial dicarboxylate transporter (1162 nt)
TERG_01604−1.9760S ribosomal protein L36Trichophyton tonsurans CBS 112818 60S ribosomal protein L36 (527 nt)
TERG_03148−1.95 Hypothetical protein Trichophyton equinum CBS 127.97 molybdenum cofactor sulfurase (1551 nt)
TERG_05236−1.9460S ribosomal protein L35Trichophyton tonsurans CBS 112818 60S ribosomal protein L35 (1143 nt)
TERG_01399−1.93 Hypothetical protein -
TERG_06755−1.91C-8 sterol isomeraseTrichophyton verrucosum HKI 0517 C-8 sterol isomerase (Erg-1), putative (628 nt)
TERG_02542−1.91 Integral membrane proteinTrichophyton verrucosum HKI 0517 integral membrane protein Pth11-like, putative (1773 nt)

Share and Cite

MDPI and ACS Style

Petrucelli, M.F.; Matsuda, J.B.; Peroni, K.; Sanches, P.R.; Silva, W.A., Jr.; Beleboni, R.O.; Martinez-Rossi, N.M.; Marins, M.; Fachin, A.L. The Transcriptional Profile of Trichophyton rubrum Co-Cultured with Human Keratinocytes Shows New Insights about Gene Modulation by Terbinafine. Pathogens 2019, 8, 274. https://doi.org/10.3390/pathogens8040274

AMA Style

Petrucelli MF, Matsuda JB, Peroni K, Sanches PR, Silva WA Jr., Beleboni RO, Martinez-Rossi NM, Marins M, Fachin AL. The Transcriptional Profile of Trichophyton rubrum Co-Cultured with Human Keratinocytes Shows New Insights about Gene Modulation by Terbinafine. Pathogens. 2019; 8(4):274. https://doi.org/10.3390/pathogens8040274

Chicago/Turabian Style

Petrucelli, Monise Fazolin, Josie Budag Matsuda, Kamila Peroni, Pablo Rodrigo Sanches, Wilson Araújo Silva, Jr., Rene Oliveira Beleboni, Nilce Maria Martinez-Rossi, Mozart Marins, and Ana Lúcia Fachin. 2019. "The Transcriptional Profile of Trichophyton rubrum Co-Cultured with Human Keratinocytes Shows New Insights about Gene Modulation by Terbinafine" Pathogens 8, no. 4: 274. https://doi.org/10.3390/pathogens8040274

APA Style

Petrucelli, M. F., Matsuda, J. B., Peroni, K., Sanches, P. R., Silva, W. A., Jr., Beleboni, R. O., Martinez-Rossi, N. M., Marins, M., & Fachin, A. L. (2019). The Transcriptional Profile of Trichophyton rubrum Co-Cultured with Human Keratinocytes Shows New Insights about Gene Modulation by Terbinafine. Pathogens, 8(4), 274. https://doi.org/10.3390/pathogens8040274

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop