StuA-Regulated Processes in the Dermatophyte Trichophyton rubrum: Transcription Profile, Cell-Cell Adhesion, and Immunomodulation

Fungal infections represent a significant concern worldwide, contributing to human morbidity and mortality. Dermatophyte infections are among the most significant mycoses, and Trichophyton rubrum appears to be the principal causative agent. Thus, an understanding of its pathophysiology is urgently required. Several lines of evidence have demonstrated that the APSES family of transcription factors (Asm1p, Phd1p, Sok2p, Efg1p, and StuA) is an important point of vulnerability in fungal pathogens and a potential therapeutic target. These transcription factors are unique to fungi, contributing to cell differentiation and adaptation to environmental cues and virulence. It has recently been demonstrated that StuA plays a pleiotropic role in dermatophyte pathophysiology. It was suggested that it functions as a mediator of crosstalk between different pathways that ultimately contribute to adaptive responses and fungal-host interactions. The complex regulation of StuA and its interaction pathways are yet to be unveiled. Thus, this study aimed to gain a deeper understanding of StuA-regulated processes in T. rubrum by assessing global gene expression following growth on keratin or glucose sources. The data showed the involvement of StuA in biological processes related to central carbon metabolism and glycerol catabolism, reactive oxygen species metabolism, and cell wall construction. Changes in carbohydrate metabolism may be responsible for the significant alteration in cell wall pattern and consequently in cell-cell interaction and adhesion. Loss of StuA led to impaired biofilm production and promoted proinflammatory cytokine secretion in a human keratinocyte cell line. We also observed the StuA-dependent regulation of catalase genes. Altogether, these data demonstrate the multitude of regulatory targets of StuA with a critical role in central metabolism that may ultimately trigger a cascade of secondary effects with substantial impact on fungal physiology and virulence traits.

Introduction

The interaction between fungi and their hosts involves a plethora of molecular pathways that synergize to sense and respond to changes in the surrounding milieu (Burmester et al., 2011; Peres et al., 2016; Martinez-Rossi et al., 2017). Hence, signal transduction plays an important role to enable cell responsiveness in a stimulus-specific manner. Fungal cells have the refined ability to discriminate environmental signals and tune their metabolism to deal with changing conditions. Transduction pathways are highly conserved between divergent fungal species due to their adaptive value from an evolutionary perspective (Lengeler et al., 2000). Moreover, a dynamic and synchronized fungal-host transcription program is activated to cope with different stimuli during fungal infection. Extracellular molecules can trigger alterations in the host immune response, impacting on the fungus metabolism and promoting cell structural modifications, thereby altering the infection outcome (Munoz et al., 2019). Altogether, invasion success strictly depends on a balance between virulence, responsiveness, and the ability to inhibit the host defense mechanisms.

Dermatophyte fungal infection consists of a coordinated process that involves adhesion, penetration, and colonization of keratinized tissues. Moreover, dermatophytoses are chronic, opportunistic diseases responsible for high morbidity, with both infection and treatment dramatically decreasing the quality of life of affected individuals (Sharma et al., 2017). Species from the genus Trichophyton are mainly responsible for causing cutaneous infections in humans, and T. rubrum is the most common clinical isolate worldwide (Seebacher et al., 2008; Brown et al., 2012; Almeida et al., 2019). Anthropophilic species modulate a mild inflammatory response, causing the infection to slowly and persistently progress (Peres et al., 2010). For this reason, a broader understanding of the various mechanisms involved in the interaction between dermatophytes and their hosts is crucial for the development of efficient therapeutic strategies. In particular, transcription factors are promising candidates for new antifungal treatments, given their complex role and involvement in many aspects of virulence and fungal adaptive strategies (Shelest, 2008).

The APSES family of transcription factors (Asm1p, Phd1p, Sok2p, Efg1p, and StuA) is a conserved class of proteins that occur exclusively in fungi, and are involved in the control of cell differentiation and development in ascomycetes. An important correlation between these factors and virulence has been observed in several phytopathogenic species (Nishimura et al., 2009; Lysoe et al., 2011; Pasquali et al., 2013), as well as in Candida species in which Efg1affects filamentation and subsequent biofilm formation (Ramage et al., 2002; Connolly et al., 2013). In Aspergillus fumigatus, StuA is involved in the regulation of morphogenesis and biosynthesis of secondary metabolites (Sheppard et al., 2005). Furthermore, null mutants of stuA in the dermatophyte Anthroderma benhamiae highlighted its role in keratin degradation and sexual development (Krober et al., 2017).

To further investigate the role of StuA in dermatophytes, a stuA knockout (ΔstuA) T. rubrum strain was previously constructed by our research group. This strain exhibited several macroscopic and microscopic morphological changes during growth in different conditions including keratin, and human skin and nail fragments (Lang et al., 2020). The study revealed that this transcription factor is directly involved in fungal morphogenesis, germination, conidiation, pigmentation, mycelial hydrophobicity, and susceptibility to different stressors. The strain also exhibited decreased activity of keratinolytic proteases, insufficient capacity to form aerial hyphae, and an overall loss of fitness.

Given what is already know about StuA, it is now crucial to have a thorough understanding of its function. In this study, RNA sequencing (RNA-seq) analysis of the ΔstuA T. rubrum strain was performed to assess changes in global gene expression promoted by the presence of keratin or glucose medium compared with wild-type fungus. This approach is expected to reveal significant changes in metabolism triggered by the lack of StuA. The potential for biofilm formation and interleukin secretion of the ΔstuA strain were also evaluated.

Material and Methods

Strains and Culture Conditions

T. rubrum wild-strain CBS 118892 was provided by the Westerdijk Fungal Biodiversity Institute (Utrecht, Netherlands). Its complete genomic sequence is publicly available at the Dermatophyte Comparative Database of the Broad Institute (ftp://ftp.broadinstitute.org/pub/annotation/fungi). The null mutant strain (ΔstuA) was constructed by replacing the wild stuA gene with an inactivation cassette containing hph, which confers resistance to hygromycin (Lang et al., 2020). The wild-type and mutant strains were initially grown on malt extract agar solid medium (2% glucose, 2% malt extract, 0.1% peptone, 2% agar, pH 5.7) at 28°C for 20 days. Approximately 1 × 10⁶ conidia of each strain were inoculated into 100 mL of Sabouraud Dextrose Broth. The cultures were maintained at 28°C for 96 h under constant agitation (120 rpm).

The resulting mycelia were washed with sterile water and transferred to 100 mL of minimal medium (Cove, 1966) containing 70 mM sodium nitrate (Sigma-Aldrich, St. Louis, MO, USA) and either 50 mM glucose (Sigma-Aldrich) or 0.5% bovine keratin (m/v). Cultures were incubated for 24, 48, or 96 h at 28°C under agitation. Biological material from three independent glucose or keratin replicates was filtered and stored at −80°C.

RNA Extraction and Sequencing

Total RNA was extracted from mycelia using the Illustra RNAspin mini isolation kit (GE Healthcare, Chicago, IL, USA). The RNA concentration and quality were estimated using a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) and the Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA). cDNA synthesis was carried out with the TruSeq RNA library Kit (Illumina, San Diego, CA, USA) according to the manufacturer’s protocol. cDNA libraries were sequenced on a NextSeq 500 system (Illumina), generating 150-bp paired-end reads.

Data Analysis

Analysis of RNA sequencing (RNA-seq) data was performed as previously described (Martins et al., 2020). Briefly, a filter for quality control was applied through the FastQC tool, and Trimmomatic (Bolger et al., 2014) was used to remove adapters and Illumina-specific sequences. A cutoff threshold based on an average base quality score of 15 over a window of 4 bases was determined, and reads shorter than 36 bases in post-trimming were excluded. Then, the trimmed paired-end reads were aligned to the reference genome retrieved from Ensembl Fungi database (https://fungi.ensembl.org/Trichophyton_rubrum_cbs_118892_gca_000151425/Info/Index) using the STAR software (Dobin et al., 2013). Reads mapping to multiple locations were excluded using the STAR´s – out Filter Multimap N max 1 parameter, and gene-level read-counts were quantified using the STAR’s -quant Mode Gene Counts parameter. A visual inspection was carried out using the Integrative Genomics Viewer software (Thorvaldsdottir et al., 2013). Biological replicates were inspected through principal component analysis plots.

Differentially expressed genes (DEGs) were determined by comparing the RNA-seq data from wild-type strain, with the data obtained for the ΔstuA strain. The cultivation, RNA preparation, and sequencing of the wild-type and stuA strains were performed together. The RNA-seq data from the wild-type strain were used as the reference of gene modulation in these comparative analyses. The DESeq2 Bioconductor package (Love et al., 2014) was used to identify the DEGs. After that, the Benjamini-Hochberg correction (Benjamini and Hochberg, 1995) was applied (P < 0.05), and a cutoff threshold of ± 1.5 log2-fold change was set to reveal statistically significant expression differences. We conducted a functional categorization of DEGs according to gene ontology (GO) terms assigned by the Blast2GO algorithm. Lastly, the most representative categories for each experimental condition were identified by enrichment analysis using the BayGO algorithm (Vencio et al., 2006).

Validation by Real-Time Quantitative Reverse Transcription Polymerase Chain Reaction (RT-qPCR)

We selected a set of genes according to their prominence in fold-change and their potential implications in processes regulated by the StuA factor. Quantification was performed using the StepOnePlus Real-Time PCR system (Applied Biosystems, Waltham, MA, USA), and reactions were prepared with SYBR Green and PCR Master Mix (Life Technologies, Waltham, MA, USA). We used ROX dye as a normalizer for the fluorescent signal, and the reference genes rpb2 and gapdh served as internal controls (Jacob et al., 2012). The 2^−ΔΔCt relative quantification method was applied to calculate relative fold-changes. The oligonucleotide sequences, concentrations, and efficiency of each reaction are shown in Supplementary Table S1. All reactions were performed in triplicate.

Cytokine Secretion in Co-Culture With Human Keratinocytes (HaCaT)

The co-culture assay was performed using 1 × 10⁶ conidia/mL and 2.5 × 10⁵ HaCaT cells/mL, which were incubated for 24 h at 37°C in 5% CO₂, as previously described (Lang et al., 2020). To gain a better understanding of fungus–host interaction and the contribution of StuA to the chronicity of dermatophyte infections, the levels of interleukin (IL)-1β, IL-4, and interferon (INF)-γ in cell supernatants were assessed by enzyme-linked immunosorbent assay (Peprotech, Cranbury, NJ, USA) in accordance with the manufacturer’s recommendations.

In Vitro Characterization of Biofilm Formation

The characterization of biofilm formation by T. rubrum strains was carried out as previously described (Costa-Orlandi et al., 2014), with slight modifications. Briefly, approximately 1 × 10⁶ conidia/mL obtained from 10-day-old plates of T. rubrum were seeded in 96-well plates, and incubated for 4 h for pre-adhesion. Then, the supernatant was removed and replaced by 200 µL of RPMI 1640 medium supplemented with L-glutamine and without sodium bicarbonate (Gibco), buffered with Mops (Sigma-Aldrich), and added with 2% of glucose (w/v).

The plates were incubated at 37°C without agitation for 24, 48, 72, 96, and 120 h. The metabolic activity of biofilms was determined using the XTT reduction assay (2.3-bis(2-methoxy-4-nitro-5-sulfophenyl)-5-[carbonyl (phenylamino)]-2H-tetrazoliumhydroxide). At the end of each incubation period, 50 μL of XTT solution (1 mg∙mL⁻¹ in phosphate-buffered saline solution) plus 4 μL of menadione solution (1 mM in acetone) were added to each well followed by incubation at 37°C for 3 h. The colorimetric change in fungal biofilms, which correlates to cell viability, was determined using a microplate reader (Multiskan FC, Thermo Fisher Scientific) at 450 nm. The biomass of T. rubrum-derived biofilms was measured by crystal violet staining after 96 h. The culture medium was removed from each well and adhered cells were washed thrice with phosphate-buffered saline solution, then the wells were allowed to dry at room temperature for about 1 min. Subsequently, 100 μL of 0.5% crystal violet solution was added to each well for 5 min. The wells were then washed four times with sterile water, followed by decolorization by the addition of 100 μL of 95% ethanol solution. The plates were analyzed using a microplate reader (Multiskan FC, Thermo Fisher Scientific) at a wavelength of 550 nm.

Statistical Analyses

Statistical significance was evaluated by the unpaired Student’s t-test or one-way analysis of variance (ANOVA) followed by Tukey’s post-hoc test. The Pearson correlation coefficient was determined for the RT-qPCR and RNA-seq fold-change values. Prism v. 5.1 (GraphPad Software, San Diego, CA, USA) was used for the statistical analyses and to generate the graphs.

Results

Transcriptional Profile of the ΔstuA Mutant in the Presence of Glucose or Keratin

RNA-seq data revealed the global transcriptional changes in T. rubrum after deletion of the stuA and growth in glucose or keratin medium for three different time points. Next-generation sequencing generated 7,889,656 and 13,990,682 high-quality reads for each condition, which resulted in approximately 7,482,678 and 12,889,091 mapped paired-end sequences (Supplementary Table S2). Upregulated and downregulated transcripts were defined by a 1.5-fold change cutoff and a stringent statistical significance threshold of P < 0.05.

Comparison of the ΔstuA T. rubrum mutant and the wild-type strain (control) revealed 481, 749, and 503 genes modulated upon 24, 48, and 96 h, respectively, in the presence of glucose. As for keratin, 460, 689, 798 genes were modulated at 24, 48, and 96 h, respectively. A total of 1,170 DEGs were identified in samples grown in glucose and 1,230 DEGs in samples grown with keratin. Among these genes, 134 and 209 were modulated in the ΔstuA strain in a time-independent manner in response to glucose and keratin, respectively; considering only the three times analyzed.

As expected, most of the observed changes were time-dependent. Notably, the greatest number of DEGs was observed at 48 h of culture in glucose and at 96 h of culture with keratin. A previous study reported that dermatophyte growth on keratin triggers an alkalization of the culture medium at 72−96 h, which corresponds to high proteolysis activity and probably correlates with the initiation of the dermatophyte infection process (Maranhão et al., 2007). Figure 1 schematically shows the distribution of these data.

FIGURE 1

Figure 1 Time-dependent distribution of differentially expressed genes (DEGs) in the ΔstuA strain compared with that in the wild-type during growth on glucose or keratin sources. (A) Venn diagram illustrates the number of transcripts modulated in ΔstuA T. rubrum grown in glucose medium for 24, 48, or 96 h compared to wild-type. (B) Number of DEGs at each time point during ΔstuA T. rubrum growth in glucose. (C) Venn diagram illustrates the number of transcripts modulated at each time point in ΔstuA T. rubrum grown in keratin-supplemented medium. (D) Number of DEGs in ΔstuA T. rubrum grown on keratin for 24, 48, or 96 h.

A survey of the top 10 most significantly modulated genes (Supplementary Tables S3, S4) revealed several transcripts with undefined molecular functions, most of which were exclusive to one single experimental condition. Thus, it seems that future improvements in the annotation of T. rubrum and of other related species it is still needed for a better understanding of root processes, such as the ones studied here. Nonetheless, critical metabolic pathways, including the reactive oxygen species (ROS) metabolism, cell wall organization, and genes involved in protein folding and proteostasis, seemed highly relevant in this analysis. The complete DEGs list for each culture condition is shown in Supplementary Table S5.

Functional Categorization of DEGs

Following transcriptome analyses, the identified DEGs were further characterized to determine the main metabolic pathways and processes to which they are associated. This approach was expected to precisely define the StuA regulator activity implicated in these particular molecular contexts. Thus, the Blast2Go tool (Conesa et al., 2005) was used to assess the distribution of DEGs as an effect of the presence of the ΔstuA mutation.

A summary of main GO terms is hierarchically presented in Figure 2, and the full list of the enriched GO terms is shown in Supplementary Table S6. Carbohydrate metabolism, phosphorylation, oxidation-reduction, and regulation of transcription, microtubule nucleation, and chitin processes, were among the processes associated with the identified DEGs. Moreover, these main classes comprised genes involved in biological processes related to cell cycle, carbon metabolism, ROS metabolism, phosphorylation, and cell wall construction (Table 1).

FIGURE 2

Figure 2 A summary of the Gene Ontology (GO) terms related to biological processes (BP) classification. Red and green bars represent up- and down-regulated genes, respectively. Categories determined for DEGs in ΔstuA T. rubrum strain during its growth in (A) glucose- and (B) keratin-supplemented medium at each time point.

TABLE 1

Table 1 Differentially expressed genes belonging to some gene ontology terms – metabolism, oxidative stress response, phosphorylation,and reproduction – in ΔstuA mutant vs. wild‐type Trichophyton rubrum.

Validation of Transcriptome Data by RT-qPCR

Next, potential StuA-regulated genes were strategically chosen to validate the RNA-seq data. Therefore, transcripts modulated under different conditions were selected and their individual expression was assessed. RT-qPCR analysis was carried out for 16 genes related to processes such as lipid, carbohydrate, and amino acid metabolism; cell wall organization; membrane transport; signal transduction; cell cycle regulation; fungal resistance; and pathogenesis. Expression profiles are shown in Figure 3, and an overview of these results is displayed in Supplementary Table S7. A positive correlation between the RNA-seq and RT-qPCR values was observed (Pearson’s correlation, r ≥ 0.70, P < 0.001).

FIGURE 3

Figure 3 Panel of genes validated by RT-qPCR. Transcript levels are represented as fold-change values at each time point in the ΔstuA strain during grown in (A) glucose- and (B) keratin-supplemented medium relative to control (wild-type strain). Statistical significance was determined using an unpaired t-test (*P < 0.05; **P < 0.01; ***P < 0.001).

StuA Affects Cell-Cell Adhesion

Given the involvement of StuA on cell wall organization and cell wall in interactions with host tissues, the biofilm formation by the ΔstuA strain were further characterized. Biofilms are heterogeneous structures surrounded by an extracellular matrix. The formation of biofilms is associated with a set of advantages such as protection against hostile environments, microbial communication, gene expression regulation, and improvement in nutrition, growth, and excretion (Percival et al., 2012; Costa-Orlandi et al., 2017; de Barros et al., 2020). The data revealed that ΔstuA mutant had a different kinetic of mature biofilm formation and marked differences in hyphal development compared to wild-type T. rubrum (Figures 4A, B). In the wild-type strain, a significant increase in metabolic activity was observed at 72 h, with a slight increase at 96 h, followed by an activity decay at 120 h (Figure 4C), consistent with previous data (Garcia et al., 2020). In contrast, in the ΔstuA strain, an increase in metabolic activity was only detected at 120 h (Figure 4D). Moreover, the biomass production at 96 h was significantly reduced in the mutant strain (Figure 4E). This finding agreed with the metabolic activity data determined by XTT assay (Figures 4C, D).

FIGURE 4

Figure 4 StuA supports cell-cell adhesion. (A) Morphology and (B) comparative kinetics of biofilm formation in ΔstuA and wild-type strains of T. rubrum. Metabolic activity of the biofilm produced by (C) wild-type and (D) ΔstuA strains of T. rubrum at each time point. Statistical significance was determined using ANOVA followed by Tukey’s post-hoc test (*P < 0.05; **P < 0.01; ***P < 0.001). (E) Biomass production in ΔstuA and wild-type strains of T. rubrum determined at 96 h by crystal violet staining. Statistical significance was determined by unpaired T-test (*P < 0.05).

stuA Deletion Induces a Proinflammatory Response of Keratinocytes

Lastly, we evaluated the impact of stuA deletion on the cell wall composition regarding the host immune response. Conceivably, the components of T. rubrum cell wall interact with keratinocytes and, due to this interplay, a set of cytokines is released as an inflammatory response by host cells. HaCaT cells were incubated with conidia cells from wild-type and ΔstuA strains for 24 h. The culture supernatants were collected, and cytokine production was assessed. The levels of IL-1β and IFN-γ were evaluated, and a pronounced increase in proinflammatory cytokine production was observed in the ΔstuA strain (Figure 5). In turn, similarly low levels of the anti-inflammatory cytokine IL-4 was observed in both strains (data not shown).

FIGURE 5

Figure 5 stuA deletion induces a proinflammatory profile in human keratinocytes. Levels of pro-inflammatory cytokines produced by HaCaT cells after coculture with conidia from wild-type and ΔstuA strains of T. rubrum for 24 h. Comparison of (A) interleukin (IL)-1β and (B) interferon (INF)-γ secretion induced by each strain. Statistical significance was determined using ANOVA followed by Tukey’s post-hoc test (*P < 0.05; **P < 0.01).

Discussion

Under the lens of high-throughput approaches, it is possible to uncover some of the mechanisms governed by transcription factors. In this work, high-throughput RNA-seq was performed to unveil the regulatory role of StuA in T. rubrum. Our data revealed that deletion of stuA triggers a series of stress response events involving the central carbon metabolism, glycerol catabolism, ROS metabolism, and cell wall remodeling.

Previous work has demonstrated the involvement of StuA in central carbon metabolism related to glycolysis, tricarboxylic acid cycle, and amino acid synthesis, and reported slower growth of the mutant strain on glucose, sucrose, and glycerol. In contrast, supplementation of media with glutamate restored growth rates (IpCho et al., 2010). Thus, it is reasonable to hypothesize that glutamate metabolism might work as a bypass that leads to the production of osmolytes such as glycerol in ΔstuA mutant strains. Our data confirmed the role of StuA in central carbon metabolism and showed a trend towards upregulation of glutamate metabolism (TERG_00694 and TERG_06509) when the strain was grown on keratin sources. Increased glycerol catabolism was observed in the ΔstuA T. rubrum strain growth on keratin, whereas reduced glycerol catabolism was observed in glucose-supplemented culture conditions (Table 1). Previously, it has been shown that the growth of wild-type and ΔstuA strains from T. rubrum in liquid keratin media generated equivalent biomass. Simultaneously, the mutant strain showed a considerable reduction in keratinolytic activity (Lang et al., 2020). Therefore, it is conceivable that a shift in an alternative metabolism route with the activation of glycerol catabolism occurred as an adaptation to the change in preferred sources of carbon and energy, thereby buffering the stress caused by gene deletion. Pathogenic microorganisms have demonstrated the ability to adjust their metabolism according to their surrounding milieu by activating glycerol metabolism, a specific adaptation that was deemed as pathometabolism (Eisenreich et al., 2015; Blotz and Stulke, 2017). Here, we showed that dermatophytes can exhibit this metabolic adaptation mechanism, which could be another survival strategy in the host.

The present data also suggests that the impact of the ΔstuA mutation on central carbon metabolism and osmolyte synthesis may be responsible for secondary effects in fungal cells. Indeed, a recent study suggested the involvement of StuA in osmoregulation, with impact on cell wall construction and conidia germination (Lang et al., 2020). The consumption of osmolytes is closely related to conidial germination and viability (Hallsworth and Magan, 1995), as they are responsible for generating a driving force for conidial development that ultimately controls the turgor pressure to counteract the rigid cell wall barrier (Foster et al., 2017). Crosstalk between some pathways is evidenced during these processes, mainly represented by the mitogen-activated protein kinase (MAPK) cell wall integrity and high-osmolarity glycerol (Hog1) pathways (Hohmann, 2009). Indeed, the present data suggested that some genes belonging to these pathways appeared to be regulated by StuA (GO:0016310 – phosphorylation), such as genes encoding serine/threonine-protein kinases -STE family (TERG_08066; TERG_04058; TERG_04731), and calmodulin-dependent kinases (CAMK; TERG_06366; TERG_06761), as shown in Table 1. Another line of evidence suggests that MAPK and CAMK are associated with conidial germination, fungal development, and ROS metabolism (Riquelme et al., 2018).

It has also been demonstrated that ΔstuA mutant strains in species such as A. fumigatus and T. rubrum have a rapid germination rate and low conidiation production (Sheppard et al., 2005; Lang et al., 2020). In microbes, growth rate and conidia germination are used to determine stress status. It is assumed that the optimal growth rate might trigger lower cellular resilience, with little competitive ability. In contrast, a rapid germination rate stretches the metabolism too far and can impact on cell function and structure, mainly due to increased ROS production (Hallsworth, 2018). Previously, a reduction in conidiation was demonstrated in an ΔstuA Aspergillus nidulans mutant (Miller et al., 1992; Melin et al., 2003), suggesting a cooperative relationship between stuA, bristle (brlA), and abacus (abaA) for asexual development (Miller et al., 1992). It has been demonstrated that StuA is responsible for regulating the abaA-brlA pathway in A. nidulans, indicated by the presence of genes that display StuAp response elements (StREs) upstream these regulatory developmental pathways (Dutton et al., 1997). Thus, we verified that genes encoding regulators of cell division, such as ascospore- and cyclin-related genes, were mainly upregulated in the ΔstuA T. rubrum strain. In contrast, genes involved in the asexual reproduction FluG-pathway (TERG_02133 and TERG_02368) were decreased in the ΔstuA T. rubrum strain (Table 1).

In addition, RNA-seq analyses revealed a reduction in the transcript levels of catalase genes in ΔstuA T. rubrum strain (TERG_06053, TERG_02005), suggesting a StuA-dependent regulation of catalase genes (Table 1). This finding agrees with other reports that showed the involvement of StuA regulation in catalase expression and an increase in susceptibility to oxidative stress (Scherer et al., 2002; Sheppard et al., 2005).

The ΔstuA T. rubrum strain presented hyphae thickness, increased conidia size, hypertrophy, and structures similar to that of chlamydoconidia, unlike the wild-type strain (Lang et al., 2020). In addition, the conidial wall thickness of the mutant strain was about 1.6 times greater than that of the wild-type strain (results not shown). Here, we showed that this mutant presented differently expressed genes encoding chitinases, glucanases, and cell wall integrity-related genes (Table 1). In addition, a LysM-domain protein was downregulated in our stuA mutant strain. Similarly, a previous report also showed downregulation of a LysM-encoding protein in a strain harboring another mutated APSES transcription factor (Sarmiento-Villamil et al., 2018). Fungal LysM-containing proteins contribute to host immune evasion by damping the free fungal chitin molecules (de Jonge et al., 2010). In addition, a previous study also suggested the involvement of LysM-domain genes in dermatophytes degradation of keratin (Lopes et al., 2019).

Given the relationship between cell wall composition and immune response, we sought to investigate the inflammatory response triggered by the ΔstuA mutant strain. Our data showed a shift in HaCaT cells towards a proinflammatory profile upon coculture with the mutant strain compared with the wild-type strain. These findings suggest a role for StuA in immune evasion and/or chronicity of skin fungal infections. Noteworthy, the faster germination rate of the ΔstuA mutant strain, as previously shown (Lang et al., 2020), might also be related to a more robust proinflammatory response and consequent fungal clearance (Rosowski et al., 2018). Spore germination is followed by cell wall remodeling. The altered cell wall composition could also impact on the adhesion of conidia/hyphae to other surfaces or fungal cells, as previously demonstrated for a null sitA strain of A. fumigatus, which identified considerable differences in cell wall surface proteins (Bom et al., 2015). Moreover, a previous report showed the involvement of Efg1 and C2H2 transcription factors, amino acid metabolism, and cell wall genes in the biofilm environment (Fanning and Mitchell, 2012). This study demonstrated that more than 50% of cell wall genes from A. fumigatus are changed in the biofilm habitat, which are expected to be related to adherence properties and sensory for adherence-induced responses. Among the surface genes potentially involved in the biofilm matrix are included those encoding the beta-glucan and polysaccharide, hydrophobins, and adhesins. (Fanning and Mitchell, 2012). Herein, we showed that deletion of stuA downregulated the genes coding the 1,3-beta-glucan synthase component FKS1, cell wall serine-threonine-rich galactomannnoprotein Mp1, and a putative hydrophobin (Table S5). Moreover, genes encoding C2H2 zinc finger domain-containing proteins were differentially modulated in the ΔstuA strain (Table S5).

Hence, we decided to evaluate the kinetics of biofilm formation and biomass production in the ΔstuA strain of T. rubrum. The involvement of StuA in spore germination and development has been widely discussed. Therefore, to rule out any of these features as interference factors in biofilm experiments, we previously evaluated the growth of wild-type and ΔstuA mutants in RPMI medium, and we obtained similar dried weights (Lang et al., 2020). Hence, based on the biofilm kinetics results, we suggested that the ΔstuA strain needed more time to reach the maturation stage than the wild-type strain. This result suggests that the mutant strain takes more time to disperse cells, colonize other surfaces, and restart the cycle of biofilm formation. Previous reports have shed light on the cross-interaction between biofilm formation and concealing of cell wall ligands, which ultimately impacts on host recognition (Lin et al., 2015; Kernien et al., 2017). In addition, biofilm formation relies on MAPK cascade signaling (Chen and Thorner, 2007) and on hydrophobin spore coating, as previously shown for RodA and DewC in A. fumigatus (Carrion Sde et al., 2013; Brown et al., 2016). In this sense, the causative effect of StuA deletion on biofilm formation and immunomodulation may occur due to deep changes in overall metabolism in the ΔstuA strain, which trigger important differences in cell wall pattern construction. We proposed a model that illustrates the current knowledge on the main processes potentially regulated by StuA in T. rubrum (Figure 6).

FIGURE 6

Figure 6 Schematic overview of the biological processes potentially regulated by StuA in the dermatophyte Trichophyton rubrum. The regulation of central carbon metabolism is cross-related to other StuA targeted-processes. Dashed arrows indicate a correlation between some downstream effects. Interwoven pathways yet to be elucidated might involve turning in and off some routes, such as the control of keratinolytic activity and the germination rates (green arrows). A conceivable hypothesis of the regulation of germination rates and cell cycle is associated with a stress response mechanism (yellow dashed arrow).

In summary, this study demonstrated that StuA plays a crucial role in overall fungal metabolism, cell wall construction, and host recognition. How these paths cross with StuA in this complex regulatory circuit remains to be investigated. Nevertheless, based on the present and previous data, the tight role of StuA in central metabolism is fundamental for cellular physiology and virulence traits.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The generated RNA-seq dataset is available at the Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo), under the accession numbers GSE163357 and GSE134406.

Author Contributions

NM-R, TB, JN-d-R, and AR designed the research. TB and JN-d-R performed the experimental design and laboratory experiments. PS performed the bioinformatics analyses. MM and EL assisted in RNA-seq construction. JB assisted in biofilm assays. TB, JN-d-R, NM-R, and AR wrote the manuscript. All authors reviewed the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by grants from the Brazilian Agencies: São Paulo Research Foundation - FAPESP [proc. No. 2019/22596-9, and Fellowships No. 2015/23435-8 to TB, No. 2018/15458-6 to JN-d-R, and No. 2018/11319-1 to MM]; National Council for Scientific and Technological Development - CNPq [Grants No. 305797/2017-4 and 304989/2017-7]; Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) - Finance Code 001, and Fundação de Apoio ao Ensino, Pesquisa e Assistência - FAEPA.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank V. M. Oliveira for assistance in fungus cultivation and RNA preparations, M. Mazucato, and M. D. Martins for technical support, and C. B. Costa-Orlandi for discussions about biofilm assays.

Supplementary Material

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

References

Almeida, F., Rodrigues, M. L., Coelho, C. (2019). The Still Underestimated Problem of Fungal Diseases Worldwide. Front. Microbiol. 10, 214. doi: 10.3389/fmicb.2019.00214

PubMed Abstract | CrossRef Full Text | Google Scholar

Benjamini, Y., Hochberg, Y. (1995). Controlling the False Discovery Rate - a Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B-Stat. Methodol. 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x

CrossRef Full Text | Google Scholar

Blotz, C., Stulke, J. (2017). Glycerol Metabolism and Its Implication in Virulence in Mycoplasma. FEMS Microbiol. Rev. 41, 640–652. doi: 10.1093/femsre/fux033

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolger, A. M., Lohse, M., Usadel, B. (2014). Trimmomatic: A Flexible Trimmer for Illumina Sequence Data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | CrossRef Full Text | Google Scholar

Bom, V. L., De Castro, P. A., Winkelstroter, L. K., Marine, M., Hori, J. I., Ramalho, L. N., et al. (2015). The Aspergillus fumigatus Sita Phosphatase Homologue is Important for Adhesion, Cell Wall Integrity, Biofilm Formation, and Virulence. Eukaryot. Cell 14, 728–744. doi: 10.1128/EC.00008-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Brown, G. D., Denning, D. W., Gow, N. A., Levitz, S. M., Netea, M. G., White, T. C. (2012). Hidden Killers: Human Fungal Infections. Sci. Transl. Med. 4, 165rv113. doi: 10.1126/scitranslmed.3004404

CrossRef Full Text | Google Scholar

Brown, N. A., Ries, L. N. A., Reis, T. F., Rajendran, R., Santos, D., Correa dos Santos, R. A., et al. (2016). Rnaseq Reveals Hydrophobins That are Involved in the Adaptation of Aspergillus nidulans to Lignocellulose. Biotech. Biofuels 9, 1–17. doi: 10.1186/s13068-016-0558-2

CrossRef Full Text | Google Scholar

Burmester, A., Shelest, E., Glockner, G., Heddergott, C., Schindler, S., Staib, P., et al. (2011). Comparative and Functional Genomics Provide Insights Into the Pathogenicity of Dermatophytic Fungi. Genome Biol. 12, R7. doi: 10.1186/gb-2011-12-1-r7

PubMed Abstract | CrossRef Full Text | Google Scholar

Carrion Sde, J., Leal, S. M., Jr., Ghannoum, M. A., Aimanianda, V., Latge, J. P., Pearlman, E. (2013). The Roda Hydrophobin on Aspergillus fumigatus Spores Masks Dectin-1- and Dectin-2-Dependent Responses and Enhances Fungal Survival in Vivo. J. Immunol. 191, 2581–2588. doi: 10.4049/jimmunol.1300748

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, R. E., Thorner, J. (2007). Function and Regulation in MAPK Signaling Pathways: Lessons Learned From the Yeast Saccharomyces cerevisiae. Biochim. Biophys. Acta 1773, 1311–1340. doi: 10.1016/j.bbamcr.2007.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Conesa, A., Gotz, S., Garcia-Gomez, J. M., Terol, J., Talon, M., Robles, M. (2005). Blast2GO: A Universal Tool for Annotation, Visualization and Analysis in Functional Genomics Research. Bioinformatics 21, 3674–3676. doi: 10.1093/bioinformatics/bti610

PubMed Abstract | CrossRef Full Text | Google Scholar

Connolly, L. A., Riccombeni, A., Grozer, Z., Holland, L. M., Lynch, D. B., Andes, D. R., et al. (2013). The APSES Transcription Factor Efg1 is a Global Regulator That Controls Morphogenesis and Biofilm Formation in Candida parapsilosis. Mol. Microbiol. 90, 36–53. doi: 10.1111/mmi.12345

PubMed Abstract | CrossRef Full Text | Google Scholar

Costa-Orlandi, C. B., Sardi, J. C. O., Pitangui, N. S., De Oliveira, H. C., Scorzoni, L., Galeane, M. C., et al. (2017). Fungal Biofilms and Polymicrobial Diseases. J. Fungi (Basel) 3, 22. doi: 10.3390/jof3020022

CrossRef Full Text | Google Scholar

Costa-Orlandi, C. B., Sardi, J. C., Santos, C. T., Fusco-Almeida, A. M., Mendes-Giannini, M. J. (2014). In Vitro Characterization of Trichophyton rubrum and T. mentagrophytes Biofilms. Biofouling 30, 719–727. doi: 10.1080/08927014.2014.919282

PubMed Abstract | CrossRef Full Text | Google Scholar

Cove, D. J. (1966). The Induction and Repression of Nitrate Reductase in the Fungus Aspergillus nidulans. Biochim. Biophys. Acta 113, 51–56. doi: 10.1016/S0926-6593(66)80120-0

PubMed Abstract | CrossRef Full Text | Google Scholar

de Barros, P. P., Rossoni, R. D., De Souza, C. M., Scorzoni, L., Fenley, J. C., Junqueira, J. C. (2020). Candida Biofilms: An Update on Developmental Mechanisms and Therapeutic Challenges. Mycopathologia 185, 415–424. doi: 10.1007/s11046-020-00445

PubMed Abstract | CrossRef Full Text | Google Scholar

de Jonge, R., Van Esse, H. P., Kombrink, A., Shinya, T., Desaki, Y., Bours, R., et al. (2010). Conserved Fungal LysM Effector Ecp6 Prevents Chitin-Triggered Immunity in Plants. Science 329, 953–955. doi: 10.1126/science.1190859

PubMed Abstract | CrossRef Full Text | Google Scholar

Dobin, A., Davis, C. A., Schlesinger, F., Drenkow, J., Zaleski, C., Jha, S., et al. (2013). STAR: Ultrafast Universal RNA-Seq Aligner. Bioinformatics 29, 15–21. doi: 10.1093/bioinformatics/bts635

PubMed Abstract | CrossRef Full Text | Google Scholar

Dutton, J. R., Johns, S., Miller, B. L. (1997). StuAp is a Sequence-Specific Transcription Factor That Regulates Developmental Complexity in Aspergillus nidulans. EMBO J. 16, 5710–5721. doi: 10.1093/emboj/16.18.5710

PubMed Abstract | CrossRef Full Text | Google Scholar

Eisenreich, W., Heesemann, J., Rudel, T., Goebel, W. (2015). Metabolic Adaptations of Intracellullar Bacterial Pathogens and Their Mammalian Host Cells During Infection (“Pathometabolism”). Microbiol. Spectr. 3, 27–58. doi: 10.1128/microbiolspec.MBP-0002-2014

CrossRef Full Text | Google Scholar

Fanning, S., Mitchell, A. P. (2012). Fungal Biofilms. PLoS Pathog. 8, e1002585. doi: 10.1371/journal.ppat.1002585

PubMed Abstract | CrossRef Full Text | Google Scholar

Foster, A. J., Ryder, L. S., Kershaw, M. J., Talbot, N. J. (2017). The Role of Glycerol in the Pathogenic Lifestyle of the Rice Blast Fungus Magnaporthe oryzae. Environ. Microbiol. 19, 1008–1016. doi: 10.1111/1462-2920.13688

PubMed Abstract | CrossRef Full Text | Google Scholar

Garcia, L. M., Costa-Orlandi, C. B., Bila, N. M., Vaso, C. O., Goncalves, L. N. C., Fusco-Almeida, A. M., et al. (2020). A Two-Way Road: Antagonistic Interaction Between Dual-Species Biofilms Formed by Candida albicans/Candida parapsilosis and Trichophyton rubrum. Front. Microbiol. 11, 1980. doi: 10.3389/fmicb.2020.01980

PubMed Abstract | CrossRef Full Text | Google Scholar

Hallsworth, J. E. (2018). Stress-Free Microbes Lack Vitality. Fungal Biol. 122, 379–385. doi: 10.1016/j.funbio.2018.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Hallsworth, J. E., Magan, N. (1995). Manipulation of Intracellular Glycerol and Erythritol Enhances Germination of Conidia At Low Water Availability. Microbiology 141 ( Pt 5), 1109–1115. doi: 10.1099/13500872-141-5-1109

PubMed Abstract | CrossRef Full Text | Google Scholar

Hohmann, S. (2009). Control of High Osmolarity Signalling in the Yeast Saccharomyces cerevisiae. FEBS Lett. 583, 4025–4029. doi: 10.1016/j.febslet.2009.10.069

PubMed Abstract | CrossRef Full Text | Google Scholar

IpCho, S. V., Tan, K. C., Koh, G., Gummer, J., Oliver, R. P., Trengove, R. D., et al. (2010). The Transcription Factor StuA Regulates Central Carbon Metabolism, Mycotoxin Production, and Effector Gene Expression in the Wheat Pathogen Stagonospora nodorum. Eukaryot. Cell 9, 1100–1108. doi: 10.1128/EC.00064-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Jacob, T. R., Peres, N. T., Persinoti, G. F., Silva, L. G., Mazucato, M., Rossi, A., et al. (2012). rpb2 is a Reliable Reference Gene for Quantitative Gene Expression Analysis in the Dermatophyte Trichophyton rubrum. Med. Mycol. 50, 368–377. doi: 10.3109/13693786.2011.616230

PubMed Abstract | CrossRef Full Text | Google Scholar

Kernien, J. F., Snarr, B. D., Sheppard, D. C., Nett, J. E. (2017). The Interface Between Fungal Biofilms and Innate Immunity. Front. Immunol. 8 1968, 1–10. doi: 10.3389/fimmu.2017.01968

CrossRef Full Text | Google Scholar

Krober, A., Etzrodt, S., Bach, M., Monod, M., Kniemeyer, O., Staib, P., et al. (2017). The Transcriptional Regulators SteA and StuA Contribute to Keratin Degradation and Sexual Reproduction of the Dermatophyte Arthroderma benhamiae. Curr. Genet. 63, 103–116. doi: 10.1007/s00294-016-0608-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Lang, E. A. S., Bitencourt, T. A., Peres, N. T. A., Lopes, L., Silva, L. G., Cazzaniga, R. A., et al. (2020). The stuA Gene Controls Development, Adaptation, Stress Tolerance, and Virulence of the Dermatophyte Trichophyton rubrum. Microbiol. Res. 241, 126592. doi: 10.1016/j.micres.2020.126592

PubMed Abstract | CrossRef Full Text | Google Scholar

Lengeler, K. B., Davidson, R. C., D’souza, C., Harashima, T., Shen, W. C., Wang, P., et al. (2000). Signal Transduction Cascades Regulating Fungal Development and Virulence. Microbiol. Mol. Biol. Rev. 64, 746–785. doi: 10.1128/MMBR.64.4.746-785.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, C. J., Sasse, C., Gerke, J., Valerius, O., Irmer, H., Frauendorf, H., et al. (2015). Transcription Factor Soma is Required for Adhesion, Development and Virulence of the Human Pathogen Aspergillus fumigatus. PLoS Pathog. 11, e1005205. doi: 10.1371/journal.ppat.1005205

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopes, L., Bitencourt, T. A., Lang, E. A. S., Sanches, P. R., Peres, N. T. A., Rossi, A., et al. (2019). Genes Coding for Lysm Domains in the Dermatophyte Trichophyton rubrum: A Transcription Analysis. Med. Mycol. 58, 372–379. doi: 10.1093/mmy/myz068

CrossRef Full Text | Google Scholar

Love, M. I., Huber, W., Anders, S. (2014). Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data With Deseq2. Genome Biol. 15, 550. doi: 10.1186/s13059-014-0550-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Lysoe, E., Pasquali, M., Breakspear, A., Kistler, H. C. (2011). The Transcription Factor FgStuAp Influences Spore Development, Pathogenicity, and Secondary Metabolism in Fusarium graminearum. Mol. Plant Microbe Interact. 24, 54–67. doi: 10.1094/MPMI-03-10-0075

PubMed Abstract | CrossRef Full Text | Google Scholar

Maranhão, F. C. A., Paião, F. G., Martinez-Rossi, N. M. (2007). Isolation of Transcripts Over-Expressed in Human Pathogen Trichophyton rubrum During Growth in Keratin. Microb. Pathog. 43, 166–172. doi: 10.1016/j.micpath.2007.05.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Martinez-Rossi, N. M., Peres, N. T., Rossi, A. (2017). Pathogenesis of Dermatophytosis: Sensing the Host Tissue. Mycopathologia 182, 215–227. doi: 10.1007/s11046-016-0057-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Martins, M. P., Rossi, A., Sanches, P. R., Bortolossi, J. C., Martinez-Rossi, N. M. (2020). Comprehensive Analysis of the Dermatophyte Trichophyton rubrum Transcriptional Profile Reveals Dynamic Metabolic Modulation. Biochem. J. 477, 873–885. doi: 10.1042/BCJ20190868

PubMed Abstract | CrossRef Full Text | Google Scholar

Melin, P., Schnurer, J., Wagner, E. G. (2003). Characterization of Phia, a Gene Essential for Phialide Development in Aspergillus nidulans. Fungal Genet. Biol. 40, 234–241. doi: 10.1016/S1087-1845(03)00108-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, K. Y., Wu, J., Miller, B. L. (1992). StuA is Required for Cell Pattern Formation in Aspergillus. Genes Dev. 6, 1770–1782. doi: 10.1101/gad.6.9.1770

PubMed Abstract | CrossRef Full Text | Google Scholar

Munoz, J. F., Delorey, T., Ford, C. B., Li, B. Y., Thompson, D. A., Rao, R. P., et al. (2019). Coordinated Host-Pathogen Transcriptional Dynamics Revealed Using Sorted Subpopulations and Single Macrophages Infected With Candida albicans. Nat. Commun. 10, 1607. doi: 10.1038/s41467-019-09599-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Nishimura, M., Fukada, J., Moriwaki, A., Fujikawa, T., Ohashi, M., Hibi, T., et al. (2009). Mstu1, an APSES Transcription Factor, is Required for Appressorium-Mediated Infection in Magnaporthe grisea. Biosci. Biotechnol. Biochem. 73, 1779–1786. doi: 10.1271/bbb.90146

PubMed Abstract | CrossRef Full Text | Google Scholar

Pasquali, M., Spanu, F., Scherm, B., Balmas, V., Hoffmann, L., Hammond-Kosack, K. E., et al. (2013). FcStuA From Fusarium culmorum Controls Wheat Foot and Root Rot in a Toxin Dispensable Manner. PLoS One 8, e57429. doi: 10.1371/journal.pone.0057429

PubMed Abstract | CrossRef Full Text | Google Scholar

Percival, S. L., Emanuel, C., Cutting, K. F., Williams, D. W. (2012). Microbiology of the Skin and the Role of Biofilms in Infection. Int. Wound J. 9, 14–32. doi: 10.1111/j.1742-481X.2011.00836.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Peres, N. T., Maranhao, F. C., Rossi, A., Martinez-Rossi, N. M. (2010). Dermatophytes: Host-Pathogen Interaction and Antifungal Resistance. Bras. Dermatol. 85, 657–667. doi: 10.1590/S0365-05962010000500009

CrossRef Full Text | Google Scholar

Peres, N. T., Silva, L. G., Santos Rda, S., Jacob, T. R., Persinoti, G. F., Rocha, L. B., et al. (2016). In Vitro and Ex Vivo Infection Models Help Assess the Molecular Aspects of the Interaction of Trichophyton rubrum With the Host Milieu. Med. Mycol. 54, 420–427. doi: 10.1093/mmy/myv113

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramage, G., Vandewalle, K., Lopez-Ribot, J. L., Wickes, B. L. (2002). The Filamentation Pathway Controlled by the Efg1 Regulator Protein is Required for Normal Biofilm Formation and Development in Candida albicans. FEMS Microbiol. Lett. 214, 95–100. doi: 10.1111/j.1574-6968.2002.tb11330.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Riquelme, M., Aguirre, J., Bartnicki-Garcia, S., Braus, G. H., Feldbrugge, M., Fleig, U., et al. (2018). Fungal Morphogenesis, From the Polarized Growth of Hyphae to Complex Reproduction and Infection Structures. Microbiol. Mol. Biol. Rev. 82, 1–47. doi: 10.1128/MMBR.00068-17

CrossRef Full Text | Google Scholar

Rosowski, E. E., Raffa, N., Knox, B. P., Golenberg, N., Keller, N. P., Huttenlocher, A. (2018). Macrophages Inhibit Aspergillus fumigatus Germination and Neutrophil-Mediated Fungal Killing. PLoS Pathog. 14, e1007229. doi: 10.1371/journal.ppat.1007229

PubMed Abstract | CrossRef Full Text | Google Scholar

Sarmiento-Villamil, J. L., Garcia-Pedrajas, N. E., Baeza-Montanez, L., Garcia-Pedrajas, M. D. (2018). The APSES Transcription Factor Vst1 is a Key Regulator of Development in Microsclerotium- and Resting Mycelium-Producing Verticillium Species. Mol. Plant Pathol. 19, 59–76. doi: 10.1111/mpp.12496

PubMed Abstract | CrossRef Full Text | Google Scholar

Scherer, M., Wei, H., Liese, R., Fischer, R. (2002). Aspergillus nidulans Catalase-Peroxidase Gene (Cpea) is Transcriptionally Induced During Sexual Development Through the Transcription Factor StuA. Eukaryot. Cell 1, 725–735. doi: 10.1128/EC.1.5.725-735.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

Seebacher, C., Bouchara, J. P., Mignon, B. (2008). Updates on the Epidemiology of Dermatophyte Infections. Mycopathologia 166, 335–352. doi: 10.1007/s11046-008-9100-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, R., Adhikari, L., Sharma, R. L. (2017). Recurrent Dermatophytosis: A Rising Problem in Sikkim, a Himalayan State of India. Indian J. Pathol. Microbiol. 60, 541–545. doi: 10.4103/IJPM.IJPM_831_16

PubMed Abstract | CrossRef Full Text | Google Scholar

Shelest, E. (2008). Transcription Factors in Fungi. FEMS Microbiol. Lett. 286, 145–151. doi: 10.1111/j.1574-6968.2008.01293.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheppard, D. C., Doedt, T., Chiang, L. Y., Kim, H. S., Chen, D., Nierman, W. C., et al. (2005). The Aspergillus fumigatus StuA Protein Governs the Up-Regulation of a Discrete Transcriptional Program During the Acquisition of Developmental Competence. Mol. Biol. Cell 16, 5866–5879. doi: 10.1091/mbc.e05-07-0617

PubMed Abstract | CrossRef Full Text | Google Scholar

Thorvaldsdottir, H., Robinson, J. T., Mesirov, J. P. (2013). Integrative Genomics Viewer (IGV): High-Performance Genomics Data Visualization and Exploration. Brief Bioinform. 14, 178–192. doi: 10.1093/bib/bbs017

PubMed Abstract | CrossRef Full Text | Google Scholar

Vencio, R. Z., Koide, T., Gomes, S. L., Pereira, C. A. (2006). Baygo: Bayesian Analysis of Ontology Term Enrichment in Microarray Data. BMC Bioinf. 7, 86. doi: 10.1186/1471-2105-7-86

CrossRef Full Text | Google Scholar