Next Article in Journal
Phylogenetic and Taxonomic Analyses of Three New Wood-Inhabiting Fungi of Xylodon (Basidiomycota) in a Forest Ecological System
Next Article in Special Issue
Metabolomics Analysis and Antioxidant Potential of Endophytic Diaporthe fraxini ED2 Grown in Different Culture Media
Previous Article in Journal
Evaluation of the Antifungal Activities of Photorhabdus akhurstii and Its Secondary Metabolites against Phytopathogenic Colletotrichum gloeosporioides
Previous Article in Special Issue
A Comprehensive Assessment of the Secretome Responsible for Host Adaptation of the Legume Root Pathogen Aphanomyces euteiches
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JoF
Volume 8
Issue 4
10.3390/jof8040404
jof-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
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Differential Gene Expression of Mucor lusitanicus under Aerobic and Anaerobic Conditions

by
Mónika Homa
1,
Sandugash Ibragimova
1,
Csilla Szebenyi
1,
Gábor Nagy
2,
Nóra Zsindely
3,
László Bodai
2,
Csaba Vágvölgyi
1,
Gábor Nagy
1 and
Tamás Papp
1,*
1
Department of Microbiology, Faculty of Science and Informatics, University of Szeged, H-6726 Szeged, Hungary
2
Department of Biochemistry and Molecular Biology, Faculty of Science and Informatics, University of Szeged, H-6726 Szeged, Hungary
3
Department of Genetics, Faculty of Science and Informatics, University of Szeged, H-6726 Szeged, Hungary
*
Author to whom correspondence should be addressed.
J. Fungi 2022, 8(4), 404; https://doi.org/10.3390/jof8040404
Submission received: 21 February 2022 / Revised: 12 April 2022 / Accepted: 14 April 2022 / Published: 15 April 2022
(This article belongs to the Special Issue Fungi: What Have We Learned from Omics?)
Browse Figures
Review Reports Versions Notes

Abstract

:
Mucor lusitanicus and some other members of the fungal order Mucorales display the phenomenon of morphological dimorphism. This means that these fungi aerobically produce filamentous hyphae, developing a coenocytic mycelium, but they grow in a multipolar yeast-like form under anaerobiosis. Revealing the molecular mechanism of the reversible yeast-hyphal transition can be interesting for both the biotechnological application and in the understanding of the pathomechanism of mucormycosis. In the present study, transcriptomic analyses were carried out after cultivating the fungus either aerobically or anaerobically revealing significant changes in gene expression under the two conditions. In total, 539 differentially expressed genes (FDR < 0.05, |log2FC| ≥ 3) were identified, including 190 upregulated and 349 downregulated transcripts. Within the metabolism-related genes, carbohydrate metabolism was proven to be especially affected. Anaerobiosis also affected the transcription of transporters: among the 14 up- and 42 downregulated transporters, several putative sugar transporters were detected. Moreover, a considerable number of transcripts related to amino acid transport and metabolism, lipid transport and metabolism, and energy production and conversion were proven to be downregulated when the culture had been transferred into an anaerobic atmosphere.

1. Introduction

Mucor lusitanicus (also known as M. circinelloides f. lusitanicus) is a filamentous fungus belonging to the order Mucorales [1]. Several Mucor species have biotechnological significance as producers of extracellular enzymes [2], carotenoids [3], unsaturated fatty acids, and organic acids [4,5,6]. Certain species, such as members of the genera Rhizopus and Mucor, are also known as agents of frequently fatal opportunistic fungal infections, the so-called mucormycoses [7,8]. Such infections associated with rapid progression and high mortality rates most frequently occur in immunocompromised patients and those with diabetic ketoacidosis and are notorious for the difficulties in their diagnosis and treatment [7,9]. Furthermore, M. lusitanicus and the closely related M. circinelloides (i.e., M. circinelloides f. circinelloides) are frequently used as a model organism during studies among others on carotenoid biosynthesis [10,11], biofuel production [12,13], morphogenesis [14,15,16], RNA interference [17], and pathogenicity mechanisms [18,19,20], as well as the molecular mechanisms of all these biological processes [21,22,23].
Mucor lusitanicus, similarly to several other mucoralean fungi, has the interesting ability of morphological dimorphism [14,24]. This fungus aerobically produces filamentous hyphae developing a coenocytic mycelium, but it grows in a multipolar yeast-like form under anaerobiosis [24,25]. Blockage of the mitochondrial energy-producing functions, such as electron transfer or oxidative phosphorylation, can also induce yeast-like growth, even under aerobic conditions [26]. The morphological transitions reversibly follow the changes of the environmental conditions [27].
cAMP-dependent protein kinase A (PKA) and calcineurin have been found as the main factors in the regulation of the dimorphic transition [15,16,26,28,29]. The addition of cAMP to the cultures can induce yeast-like growth even under aerobic conditions and the amount of PKA increases in the yeast-like cells, indicating that PKA may have a special role in the induction and/or maintenance of the yeast morphology [15]. At the same time, calcineurin plays a crucial role in the yeast-hyphal transition, and in the case of its deletion or inhibition, yeast cells fail to switch to hyphal growth [16].
Control and optimization of fungal morphology and growth are essential for various fermentation-based applications [24,25]. Morphological dimorphism may also have consequences on virulence, and it is regarded as a virulence factor for several pathogenic fungi [30]. In the case of M. lusitanicus, Lee et al. [16] found that the yeast-like form is less virulent than the spores or the hyphal form and raised the possibility of targeting the morphological transition to address the infection (e.g., by using calcineurin inhibitors) [16,20]. Thus, revealing the molecular mechanism of the yeast-hyphal transition can be interesting for both the potential biotechnological applications and developments aimed at controlling mucormycosis. In the present study, transcriptomic analyses were carried out after cultivating the fungus either aerobically or anaerobically. Transcriptomes were compared and genes differentially expressed were analyzed.

2. Materials and Methods

2.1. Strain and Culture Conditions

The M. lusitanicus double auxotrophic strain MS12 (leuA− and pyrG−), a derivative of the strain CBS 277.49 [31], was used in the present study. For spore collection, the strain was grown for 7 days on malt extract agar (MEA; 1% glucose, 0.25% yeast extract, 5% malt extract, 2% agar; pH 5.4) plates at 25 °C; then, the spores were washed with PBS by gentle scraping, and the spore suspensions obtained were maintained at 4 °C until use. For aerobic growth, 30 mL liquid minimal medium (YNB; 1% glucose, 0.15% ammonium sulphate, 0.15% sodium glutamate, 0.05% yeast nitrogen base without amino acids (BD Difco, Becton Dickinson, Franklin Lakes, NJ, USA), supplemented with 0.05% leucine and 0.05% uracil; pH 6.8) was inoculated to obtain a final concentration of 104 sporangiospores/mL. The cultures were kept at 25 °C for 24 h with constant shaking (i.e., 190 rpm). For the mycelium–yeast transition, the 24 h old culture grown aerobically was transferred into a 50 mL flask into 50 mL fresh YNB medium and placed for 4 h into a BBL GasPak Anaerobic System (Becton Dickinson, Franklin Lakes, NJ, USA) at 25 °C.

2.2. RNA Extraction and Sequencing

For RNA purification, the 30 mL cultures detailed above were filtered through 0.45 µm MCE membrane filters (Millipore, Burlington, MA, USA) and the extraction was performed by using the Direct-zol™ RNA MiniPrep Kit (Zymo Research, Irvine, CA, USA), following the recommendations of the manufacturer. The samples were kept at −80 °C until their use. The quality of total RNA samples was analyzed by capillary gel electrophoresis in an Agilent 2100 Bioanalyzer using Agilent RNA 6000 Nano Kit and only samples with RIN > 9 were processed further. RNA-seq libraries were prepared from three biological and two technical replicates per treatment condition (i.e., six samples per treatment were used in total). Thus, six libraries were generated and all of them were sequenced. Libraries were pooled before sequencing, but since these were indexed, libraries sequence reads belonging to different samples could be demultiplexed (separated to different files) after sequencing. From 800 ng total RNA samples, polyA-RNA was isolated using NEBNext Poly(A) mRNA Magnetic Isolation Module (New England Biolabs, Ipswich, MA, USA), then indexed; strand-specific sequencing libraries were prepared using NEBNext Ultra II Directional RNA Library Prep Kit (New England Biolabs, Ipswich, MA, USA) following the protocol of the manufacturer. Sequencing libraries were validated and quantitated with an Agilent 2100 Bioanalyzer using Agilent DNA 1000 Kit; the average fragment length of the libraries was in the range of 271–289 bp. Libraries were pooled in equimolar concentrations, and after denaturing, the library pool was loaded at 15 pM concentration with 1% PhiX Control V3 (Illumina, San Diego, CA, USA) in a MiSeq Reagent Kit v3 (150-cycle) kit and sequenced in an Illumina MiSeq instrument producing 2 × 75 bp paired-end reads.

2.3. Analysis of the RNA-Seq Data

FASTQ files were generated with the GenerateFASTQ 1.1.0.64 application on Illumina BaseSpace. Sequence quality checks and adapter trimming were performed using TrimGalore/Cutadapt. Sequence reads were aligned to the M. lusitanicus CBS 277.49 v2.0 reference genome assembly available from the Joint Genome Institute (JGI) website (https://genome.jgi.doe.gov/Mucci2/Mucci2.home.html (accessed on 5 March 2021) [32]) using HISAT2 (parameters: –max-intronlen 45000, –rna-strandness RF). From the generated binary alignment files, gene specific read counts were calculated in R using the summarizeOverlaps (parameters: ignore.strand = FALSE, singleEnd = FALSE, fragments = TRUE, preprocess.reads = invertStrand) Bioconductor package with a TxDb transcript metadata object generated from the transcript annotation available from JGI using the makeTxDbFromGFF function of the GenomicFeatures Bioconductor package. Differential gene expression analysis was performed with DESeq2 parameter: p < 0.05, and a cpm > 1 read count filter. Gene Ontology analysis was performed with the topGO Bioconductor package (parameters: algorithm = “classic”, statistic = “fisher”). A gene was considered differentially expressed if the absolute log2 fold-change was greater than or equal to 3. Differentially expressed genes (DEGs) were functionally classified into KOG categories and subcategories on the basis of the manual curation of the results from the InterProScan analysis using InterProScan version 5.47–82.0 (Table S1). Additionally, this table has been supplemented with the predicted homologs of each DEG identified in the following whole proteomes of Schizosaccharomyces pombe 972 h- [33], Saccharomyces cerevisiae S288C [34], Candida albicans SC5314 [35], Cryptococcus neoformans H99 [36], M. circinelloides f. circinelloides 1006PhL [37], Neurospora crassa OR74A [38], Aspergillus fumigatus Af293 [39], Aspergillus nidulans FGSC A4 [40], Scedosporium apiospermum IHEM 14462 [41], and Anaeromyces robustus S4 [42]. To this end, we grouped proteins of the above species into protein clusters by performing an all-vs-all search using MMseqs2 [43], setting sensitivity to 5.7 and coverage to 0.8. then, the Markov cluster algorithm [44,45] was used with an inflation parameter 2.0. Orthologs were found on the basis of reciprocal best hit search using MMSeqs2 with easy-rbh function.

3. Results

Our data revealed prominent changes in the gene expression when the culture has been transferred from aerobic to anaerobic conditions. The average number of raw sequence read pairs per sample was 3.95 million (3,603,894, 3,992,427, and 3,384,845, and 4,722,290, 3,349,089, and 4,617,836 reads in the aerobic and anaerobic sample triplicates, respectively) and 3.94 million (3,596,442, 3,984,122, and 3,376,382, and 4,715,087, 3,343,617, and 4,608,645 reads in the aerobic and anaerobic sample triplicates, respectively) after trimming. From the 11,791 annotated features in the applied transcriptome annotation, 9476 genes were expressed at the level of FPKM > 1 in all aerobic samples (9580, 9587, and 9578 expressed genes in the three replicates, respectively) and 9049 genes were expressed in all anaerobic samples (9259, 9249, and 9261 expressed genes in the three replicates, respectively); 8933 genes were expressed in all replicates of both conditions. In total, 539 DEGs (FDR < 0.05, |log2FC| ≥ 3) were identified, including 190 upregulated and 349 downregulated transcripts (Figure 1 and Table S1). During functional characterization, DEGs were classified into four main functional categories. Among both down- and upregulated transcripts, the group of poorly characterized transcripts proved to be the most prominent, followed by transcripts related to “Metabolism”, “Cellular process and signaling”, and “Information storage and processing”.
Among metabolism-related genes, carbohydrate metabolism is affected the most during mycelium to yeast transition: 23 transcripts predicted to be involved in carbohydrate transport and metabolism were found to be upregulated, while 32 were downregulated. Among these transcripts, several encoded proteins were related to cell wall metabolism. For instance, transcripts encoding putative chitin deacetylases (ID 155630, 186155, 155048, 155566, and 156744) and predicted chitinases (108620 and 107467) were upregulated, while genes for other putative glycoside hydrolase/deacetylase enzymes (i.e., IDs 133575, 155735, 156524, 136560, 131995, and 155780), predicted chitin synthases (i.e., 151786 and 85917), and a putative chitinase (ID 105626) were found to be downregulated. Additionally, anaerobiosis also affected the transcription of transporters: among the 14 up- and 42 downregulated transporters, several putative sugar transporters were detected (Figure 2). Namely, we observed the upregulation of 155841, 151596, 145859, 105934, and 134811, which are putative sugar transporters belonging to the major facilitator superfamily (MFS) without any orthologs in the proteome of S. pombe, S. cerevisiae, C. albicans, or A. nidulans, and the downregulation of 153150 (an ortholog of a predicted transmembrane transporter in A. nidulans FGSC A4: AN2466 and a putative MFS glucose transporter in C. albicans SC5314 C1_02110C_ A), 155500 (an ortholog of A. nidulans FGSC A4: AN5104 with predicted transmembrane transporter activity), and 108626, among other provisionally classified proteins as carbohydrate transporters.
Additionally, a prominent number of transcripts related to amino acid transport and metabolism, lipid transport and metabolism, and energy production and conversion proved to be downregulated when the culture was transferred into an anaerobic atmosphere. In this latter group, transcripts related to oxidative metabolism were identified, such as predicted aconitases (i.e., 153298 and 186385) that catalyze citrate–isocitrate transition in TCA cycle, aldehyde dehydrogenases (i.e., 179978 and 157230), a predicted ATP-citrate lyase (i.e., 110808) that is potentially involved in the synthesis of Acetyl-CoA from citrate, predicted mitochondrial carrier proteins (e.g., 135605, 91563, 157132, 152871, 154970, and 155050), putative mitochondrial aldehyde dehydrogenases (e.g., 157230), putative glycerol-3-phosphate dehydrogenases (e.g., 153523, an ortholog of A. nidulans FGSCA4: GfdA), elements of the branched-chain α-ketoacid dehydrogenase complex (i.e., 139791 and 148011 encoding the E1 alpha and beta subunits, respectively, and 105619 encoding the E2 subunit), a putative mitochondrial external NADH dehydrogenase (i.e., 73711) that catalyzes the oxidation of cytosolic NADH, other oxidoreductases (e.g., 153534, 136103, and 122685), and a putative NADH-cytochrome b-5 reductase (i.e., 156238). Simultaneously, a small number of transcripts related to energy production and conversion were found to be upregulated, for instance, a predicted malate dehydrogenase (i.e., 186772), a putative acetyl-CoA hydrolase (i.e., 155896), a predicted fumarate reductase (104895), and a putative mitochondrial tricarboxylate/dicarboxylate carrier protein (i.e., 108484). In addition to this, we identified five predicted alcohol dehydrogenases (ADHs) among secondary metabolism-related DEGs: expression of 90838, 157439, and 160397 were downregulated, while that of 120424 and 155149 was upregulated under anaerobic conditions.
A relatively small number of transcripts with putative regulatory roles in dimorphism were also identified among the DEGs in our datasets. For example, a predicted PKAC subunit (47751) was found to be upregulated, while transcripts encoding a putative serine/threonine protein phosphatase 2A catalytic subunit (162179), a predicted cyclophilin-type peptidyl-prolyl cis-trans isomerase (154978), and a putative guanine nucleotide exchange factor (GEF; ProtID: 86922) were downregulated.
We also found 13 differentially upregulated and 18 downregulated transcription factors, and 12 up- and 20 downregulated carbohydrate-active enzymes (CAZymes) (Table S1).

4. Discussion

Anaerobiosis, together with the presence by a fermentable hexose, induces the mycelium-to-yeast transition and a switch from respiration to fermentation in Mucor species [46]. Therefore, we expected that this condition would primarily affect the properties of the cell wall and the basic metabolic and energy conservation processes of the cell.
The fungal cell wall is a strictly regulated, dynamic structure, and its major physical and chemical properties, such as thickness, structure, and composition, are mainly regulated by the cell cycle and affected by various environmental conditions, such as temperature, pH, carbon source, and oxygen level [47]. According to the recent review of Lecointe et al. [48], the main components of the cell wall of Mucorales are chitin and its deacetylated form, chitosan, β-glucans, mucoran, and mucoric acid. The chitin content of the cell walls varies according to the morphological phase of the fungus: it was proven to be lower in the yeasts while higher in the filamentous fungi [49]. Additionally, compared to other fungal species, many Basidiomycota and Mucoromycota possess a much higher chitosan level relative to chitin [50]. For instance, in C. neoformans, the cell wall has minor quantities of chitin and higher quantities of chitosan, which, among others, helps to direct cell wall integrity and bud separation [51]. In contrast to this, the chitin and chitosan content of the cell wall of Mucor rouxii proved to be very similar in the yeast and the hyphal cultures [52]. The molecular mass of mucoran was reported to be almost two times higher in the mycelial form compared to the yeast-like cells [53,54]. Conversely, mucoric acid has a similar molecular mass in sporangiophores, mycelium, or yeast-like cells [55].
According to our DEG dataset, anaerobiosis did indeed affect the expression of several cell wall-related genes. Our analysis revealed the upregulation of transcripts encoding putative chitin deacetylases and predicted chitinases. This result could suggest an altered chitosan content in the cell wall of yeast cells under anaerobic growth, which may contribute to the production of multiple budding sites, as was previously demonstrated in C. neoformans [51]. According to López-Fernández et al. [56], the chitosan content affects the virulence properties of M. lusitanicus. Thus, the fact that yeast cells of M. lusitanicus are less virulent than its hyphal form [16], among others, might be explained by the different chitosan content of the cell wall.
Some of the various glycoside hydrolases/deacetylases, the predicted chitin synthases, and a putative chitinase were also downregulated in our dataset. It was previously demonstrated by Lopez-Matas et al. [57] that the transcript of Mcchs1 (ID 114551) encoding a class II chitin synthase accumulates only during the exponentially growing hyphal stage, while it was not expressed in the yeast form. Although the transcript level of Mcchs1 proved to be unchanged during the mycelium-yeast transition of M. lusitanicus in our study, differentially downregulated transcripts encoding other chitin synthases (i.e., 151786 and 85917) were detected. These results correspond with the assumption of Lopez-Matas et al. [57] that various chitin synthase activities may have different roles in the dimorphic growth of Mucor spp. Furthermore, the transcript level of a yet uncharacterized protein coding gene (ID 118900), which has a predicted role in “cell wall, membrane and envelope biogenesis” and shows sequence similarity with the A. nidulans FGSC A4: AN1554, a protein with a predicted role in the regulation of chitin synthase activity [58], was found to be decreased under anaerobic conditions. All these results mentioned above suggest that the dimorphic switch in M. lusitanicus is accompanied by changes in the chitin and chitosan content of the cell wall.
Anaerobic conditions induced marked changes in the expression of energy production- and conversion-related genes. As was expected, numerous transcripts related to oxidative metabolism (e.g., TCA cycle and oxidative phosphorylation) proved to be downregulated when the culture has been transferred into an anaerobic atmosphere. These findings agree with previous studies describing that M. lusitanicus switches from oxidative to fermentative metabolism under anaerobic conditions in the presence of hexoses, independently of its morphology [29]. In line with this, we identified five (i.e., three down- and two upregulated) predicted class V alcohol dehydrogenases (ADHs) among secondary metabolism-related DEGs, an enzyme that is responsible for the conversion of acetaldehyde into ethanol in M. lusitanicus [59]. These results suggest that alcoholic fermentation is not restricted to anaerobic cultures, and that ethanol production is differentially regulated in M. lusitanicus depending on the presence or absence of oxygen. Our observation corresponds with the previous description of M. lusitanicus as a Crabtree-positive microorganism. This means that it can produce ethanol under anaerobic and aerobic conditions as well. Aerobic ethanol production occurs when glucose concentration is sufficiently high in the environment, but once glucose is depleted, ethanol is consumed as a carbon source [25]. We identified a total of eight putative ADHs in the genome of M. lusitanicus CBS277.49 v2.0 (Protein IDs: 90838, 120424, 157439, 160397,155149, 139628, 140177, 152844, and 34200). Although their functions are unknown, four were differentially expressed in our dataset: 90838 and 160397 were found to be downregulated, while 120424 and 155149 were upregulated. Out of these transcripts, 155149 was previously characterized by Rangel-Porras et al. [59] as the major (if not the only) ADH encoded by adh1 in M. lusitanicus, that was expressed in the cytosol of both mycelia and yeast cells. Adh1 has also been mentioned as a marker linked to fermentative metabolism under anaerobic growth in M. lusitanicus by Patiño-Medina et al. [60]. In Mucor rouxii, ADH activity was proven to be essential for anaerobic growth, most probably because of its necessity for the regeneration of NAD+ for glycolysis [61]. Furthermore, the higher transcript level of adh1 in the yeast form of M. lusitanicus was also verified by the study of Valle-Maldonado et al. [14]. Moreover, as pyruvate is not just the key intermediate in sugar metabolism and a precursor for the synthesis of several amino acids, but also an intermediate for ethanol production, we assume that the upregulation of the predicted malate dehydrogenase 186772, an ortholog of A. nidulans FGSC A4: AN6168 (MaeA) with a predicted role in the conversion of malate to pyruvate, also correlates with the high levels of anaerobic ethanol production of M. lusitanicus. As a potential proof for the elevated ethanol production in the anaerobic culture, we also observed the upregulation of a putative acetyl-CoA hydrolase, which is found to be the ortholog of A. nidulans FGSC A4: AN1547 CoA-transferase (CoaT) with an ethanol- and acetate-inducible promoter [60].
Finally, other metabolism-related transcripts, such as 104895 and 187130 encoding a putative fumarate reductase and a multicopper oxidase Fet3A, respectively, were also upregulated in the absence of oxygen, suggesting that they might have a potential role in anaerobic growth. The fumarate reductase is the ortholog of the soluble fumarate reductase of S. cerevisiae S288C: YEL047C (FRD1), which has been proven to be required with isoenzyme Osm1p for anaerobic growth on glucose as a carbon source [61]. The fet3a has been recently characterized by Navarro-Mendoza et al. [62] as a member of a gene family of three ferroxidases. These authors observed that Fet3A was specifically expressed during yeast growth under anaerobic conditions, while Fet3B (ID 50174) and Fet3C (ID 91148) were expressed in the mycelium during aerobic growth [62]. In contrast to this, we did not observe the downregulation of either Fet3B or Fet3C in the anaerobic culture.
The role of cAMP-PKA signal transduction pathway in the dimorphism of Mucor species has been reported several times previously [26,28,47,63]. The recent work of Moriwaki-Takano et al. [64] provides the first detailed hypothetical scheme detailing the elements of this pathway and their possible interactions. Briefly, this work suggests that the intracellular bicarbonic ion (HCO3−) level is increased, suppressing the transcription of Nce103 (carbonic anhydrase), Ras3GTPase, and Pde (phosphodiesterase), while inducing the Cyr1 adenylate cyclase in M. circinelloides grown in CO2 atmosphere. These latter two components are assumed to be responsible for the high cAMP level in yeast cells. cAMP induces PKA activity, which is suggested to inhibit hyphal growth through the inhibition of the elongation factor Efg1 [64]. On the basis of this study, we expected the elements of the cAMP-PKA pathway to be differentially regulated under anaerobic conditions, but only a single PKAC subunit (47751)-encoding gene was proven to be upregulated under anaerobic conditions in our dataset. Interestingly, the expression of this subunit, PKAC3 has been investigated before and found to be upregulated after a 5 h shift from anaerobic yeast-like growth to aerobic hyphal growth [65]. It is worth mentioning that orthologues of efg1 (122144) and cyr1 (126447) were found to be down- and upregulated, respectively, with log2FC values ≤ 3, i.e., −0.998 and 0.572, respectively. Similarly, small changes in the expression levels of these genes were previously detected by Moriwaki-Takano et al. [64] after 6 h of cultivation under aerobic and anaerobic conditions.
Besides the cAMP-PKA pathway, calcineurin was proven to be a key regulator of dimorphism in Mucor species [15,16]. Investigating the expression of the calcineurin pathway components during the hyphal–yeast transition in M. lusitanicus, neither the catalytic and regulatory subunits of calcineurin, nor calmodulin or FKBP12-encoding transcripts were identified as differentially regulated.
In the present study, we investigated the transcriptional response of M. lusitanicus when transferred from aerobic to anaerobic atmosphere, and thus indirectly assessed the genes involved in the mycelium-to-yeast transition as well. While some of the transcriptional changes identified correspond well with previous published studies, our results also provide new valuable data to this field of research.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/jof8040404/s1, Table S1: List and description of differentially regulated genes. Worksheet 1 summarizes differentially upregulated genes; differentially downregulated genes are listed in worksheet 2; worksheet 3 contains differentially regulated transcription factors; worksheet 4 summarizes differentially regulated transporters; and finally, differentially regulated CAZymes are listed in worksheet 5.

Author Contributions

Conceptualization, T.P., M.H., and G.N. (Gábor Nagy, [email protected]); methodology, T.P., M.H., N.Z., G.N. (Gábor Nagy, [email protected]), G.N. (Gábor Nagy, [email protected]), and L.B.; validation, T.P., M.H., G.N. (Gábor Nagy, [email protected]), G.N. (Gábor Nagy, [email protected]), and L.B.; formal analysis, M.H. and G.N. (Gábor Nagy, [email protected]); investigation, M.H., G.N. (Gábor Nagy, [email protected]), C.S., N.Z., and S.I.; data curation, M.H., N.Z., G.N. (Gábor Nagy, [email protected]), and L.B.; writing—original draft preparation, M.H., T.P., G.N. (Gábor Nagy, [email protected]), C.V., and L.B.; writing—review and editing, M.H., T.P., and L.B.; visualization, M.H.; supervision, T.P. and C.V.; funding acquisition, T.P. and G.N. (Gábor Nagy, [email protected]). All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the NKFI project K131796 and the grant ITM NKFIA TKP-2021-EGA-28. G.N. is grateful for support of the Premium Postdoctoral Fellowship Program of the Hungarian Academy of Sciences (460050). L.B. was supported by the János Bolyai Research Scholarship (BO/00522/19/8) of the Hungarian Academy of Sciences and by the ÚNKP-21-5-SZTE-574—New National Excellence Program of the Ministry for Innovation and Technology from the source of the National Research, Development, and Innovation Fund. C.S. is supported by the ÚNKP-20-4-I New National Excellence Program of the Ministry for Innovation and Technology from the source of the National Research, Development and Innovation Fund.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Trancriptomic datasets obtained during the study were deposited in the NCBI Sequence Read Archive (SRA) database (https://www.ncbi.nlm.nih.gov/sra) under the accession number PRJNA729254. Genes found to be differentially expressed (DEGs) are listed in Table S1.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Wagner, L.; Stielow, J.B.; de Hoog, G.S.; Bensch, K.; Schwartze, V.U.; Voigt, K.; Alastruey-Izquierdo, A.; Kurzai, O.; Walther, G. A new species concept for the clinically relevant Mucor circinelloides complex. Pers. Mol. Phylogeny Evol. Fungi 2020, 44, 67–97. [Google Scholar] [CrossRef] [PubMed]
  2. Voigt, K.; Wolf, T.; Ochsenreiter, K.; Nagy, G.; Kaerger, K.; Shelest, E.; Papp, T. 15 Genetic and metabolic aspects of primary and secondary metabolism of the Zygomycetes. Biochem. Mol. Biol. 2016, 3, 361–385. [Google Scholar] [CrossRef]
  3. Avalos, J.; Nordzieke, S.; Parra, O.; Pardo-Medina, J.; Carmen Limón, M. Carotenoid Production by Filamentous Fungi and Yeasts BT—Biotechnology of Yeasts and Filamentous Fungi; Sibirny, A.A., Ed.; Springer International Publishing: Cham, Switzerland, 2017; pp. 225–279. ISBN 978-3-319-58829-2. [Google Scholar]
  4. Tang, X.; Zan, X.; Zhao, L.; Chen, H.; Chen, Y.Q.; Chen, W.; Song, Y.; Ratledge, C. Proteomics analysis of high lipid-producing strain Mucor circinelloides WJ11: An explanation for the mechanism of lipid accumulation at the proteomic level. Microb. Cell Fact. 2016, 15, 35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Londoño-Hernández, L.; Ramírez-Toro, C.; Ruiz, H.A.; Ascacio-Valdés, J.A.; Aguilar-Gonzalez, M.A.; Rodríguez-Herrera, R.; Aguilar, C.N. Rhizopus oryzae—Ancient microbial resource with importance in modern food industry. Int. J. Food Microbiol. 2017, 257, 110–127. [Google Scholar] [CrossRef]
  6. Papp, T.; Nyilasi, I.; Csernetics, Á.; Nagy, G.; Takó, M.; Vágvölgyi, C. Improvement of Industrially Relevant Biological Activities in Mucoromycotina Fungi BT—Gene Expression Systems in Fungi: Advancements and Applications; Schmoll, M., Dattenböck, C., Eds.; Springer International Publishing: Cham, Switzerland, 2016; pp. 97–118. ISBN 978-3-319-27951-0. [Google Scholar]
  7. Riley, T.T.; Muzny, C.A.; Swiatlo, E.; Legendre, D.P. Breaking the mold: A review of mucormycosis and current pharmacological treatment options. Ann. Pharmacother. 2016, 50, 747–757. [Google Scholar] [CrossRef]
  8. Jeong, W.; Keighley, C.; Wolfe, R.; Lee, W.L.; Slavin, M.A.; Kong, D.C.M.; Chen, S.C.-A. The epidemiology and clinical manifestations of mucormycosis: A systematic review and meta-analysis of case reports. Clin. Microbiol. Infect. Off. Publ. Eur. Soc. Clin. Microbiol. Infect. Dis. 2019, 25, 26–34. [Google Scholar] [CrossRef] [Green Version]
  9. Skiada, A.; Lanternier, F.; Groll, A.H.; Pagano, L.; Zimmerli, S.; Herbrecht, R.; Lortholary, O.; Petrikkos, G.L. Diagnosis and treatment of mucormycosis in patients with hematological malignancies: Guidelines from the 3rd European Conference on Infections in Leukemia (ECIL 3). Haematologica 2013, 98, 492–504. [Google Scholar] [CrossRef] [Green Version]
  10. Iturriaga, E.A.; Alvarez, M.I.; Eslava, A.P.; Papp, T. Expression Vectors and Gene Fusions for the Directed Modification of the Carotenoid Biosynthesis Pathway in Mucor circinelloides BT—Microbial Carotenoids: Methods and Protocols; Barreiro, C., Barredo, J.-L., Eds.; Springer New York: New York, NY, USA, 2018; pp. 239–256. ISBN 978-1-4939-8742-9. [Google Scholar]
  11. Naz, T.; Nosheen, S.; Li, S.; Nazir, Y.; Mustafa, K.; Liu, Q.; Garre, V.; Song, Y. Comparative analysis of β-carotene production by Mucor circinelloides strains CBS 277.49 and WJ11 under light and dark conditions. Metabolites 2020, 10, 38. [Google Scholar] [CrossRef] [Green Version]
  12. Carvalho, A.K.F.; Bento, H.B.S.; Izário Filho, H.J.; de Castro, H.F. Approaches to convert Mucor circinelloides lipid into biodiesel by enzymatic synthesis assisted by microwave irradiations. Renew. Energy 2018, 125, 747–754. [Google Scholar] [CrossRef]
  13. Carvalho, A.K.F.; da Conceição, L.R.V.; Silva, J.P.V.; Perez, V.H.; de Castro, H.F. Biodiesel production from Mucor circinelloides using ethanol and heteropolyacid in one and two-step transesterification. Fuel 2017, 202, 503–511. [Google Scholar] [CrossRef]
  14. Valle-Maldonado, M.I.; Jácome-Galarza, I.E.; Gutiérrez-Corona, F.; Ramírez-Díaz, M.I.; Campos-García, J.; Meza-Carmen, V. Selection of reference genes for quantitative real time RT-PCR during dimorphism in the zygomycete Mucor circinelloides. Mol. Biol. Rep. 2015, 42, 705–711. [Google Scholar] [CrossRef]
  15. Lee, S.C.; Li, A.; Calo, S.; Heitman, J. Calcineurin plays key roles in the dimorphic transition and virulence of the human pathogenic zygomycete Mucor circinelloides. PLoS Pathog. 2013, 9, e1003625. [Google Scholar] [CrossRef] [Green Version]
  16. Lee, S.C.; Li, A.; Calo, S.; Inoue, M.; Tonthat, N.K.; Bain, J.M.; Louw, J.; Shinohara, M.L.; Erwig, L.P.; Schumacher, M.A.; et al. Calcineurin orchestrates dimorphic transitions, antifungal drug responses and host-pathogen interactions of the pathogenic mucoralean fungus Mucor circinelloides. Mol. Microbiol. 2015, 97, 844–865. [Google Scholar] [CrossRef] [Green Version]
  17. Torres-Martínez, S.; Ruiz-Vazquez, R.M. RNAi pathways in Mucor: A tale of proteins, small RNAs and functional diversity. Fungal Genet. Biol. 2016, 90, 44–52. [Google Scholar] [CrossRef] [Green Version]
  18. Li, C.H.; Cervantes, M.; Springer, D.J.; Boekhout, T.; Ruiz-Vazquez, R.M.; Torres-Martinez, S.R.; Heitman, J.; Lee, S.C. Sporangiospore size dimorphism is linked to virulence of Mucor circinelloides. PLoS Pathog. 2011, 7, e1002086. [Google Scholar] [CrossRef]
  19. Chang, Z.; Heitman, J. Drug-resistant epimutants exhibit organ-specific stability and induction during murine infections caused by the human fungal pathogen Mucor circinelloides. mBio 2019, 10, e02579-19. [Google Scholar] [CrossRef] [Green Version]
  20. Vellanki, S.; Billmyre, R.B.; Lorenzen, A.; Campbell, M.; Turner, B.; Huh, E.Y.; Heitman, J.; Lee, S.C. A novel resistance pathway for calcineurin inhibitors in the human-pathogenic mucorales mucor circinelloides. mBio 2020, 11, e02949-19. [Google Scholar] [CrossRef] [Green Version]
  21. Nagy, G.; Vaz, A.G.; Szebenyi, C.; Takó, M.; Tóth, E.J.; Csernetics, Á.; Bencsik, O.; Szekeres, A.; Homa, M.; Ayaydin, F.; et al. CRISPR-Cas9-mediated disruption of the HMG-CoA reductase genes of Mucor circinelloides and subcellular localization of the encoded enzymes. Fungal Genet. Biol. 2019, 129, 30–39. [Google Scholar] [CrossRef]
  22. Nagy, G.; Szebenyi, C.; Csernetics, Á.; Vaz, A.G.; Tóth, E.J.; Vágvölgyi, C.; Papp, T. Development of a plasmid free CRISPR-Cas9 system for the genetic modification of Mucor circinelloides. Sci. Rep. 2017, 7, 16800. [Google Scholar] [CrossRef] [Green Version]
  23. Garcia, A.; Adedoyin, G.; Heitman, J.; Lee, S.C. Construction of a recyclable genetic marker and serial gene deletions in the human pathogenic Mucorales Mucor circinelloides. G3 2017, 7, 2047–2054. [Google Scholar] [CrossRef] [Green Version]
  24. Iturriaga, A.E.; Papp, T.; Breum, J.; Arnau, J. Strain and Culture Conditions Improvement for β-Carotene Production with Mucor BT—Microbial Processes and Products; Barredo, J.-L., Ed.; Humana Press: Totowa, NJ, USA, 2005; pp. 239–256. ISBN 978-1-59259-847-2. [Google Scholar]
  25. McIntyre, M.; Breum, J.; Arnau, J.; Nielsen, J. Growth physiology and dimorphism of Mucor circinelloides (syn. racemosus) during submerged batch cultivation. Appl. Microbiol. Biotechnol. 2002, 58, 495–502. [Google Scholar] [CrossRef] [PubMed]
  26. Orlowski, M. Mucor dimorphism. Microbiol. Rev. 1991, 55, 234–258. [Google Scholar] [CrossRef] [PubMed]
  27. de la Cruz, J.O.; García-Soto, J.; Uriostegui, C.; Carranza, L.; Novoa, G.; Reyna, G.; Martínez-Cadena, G. Differential expression of Rho1GTPase and Rho3GTPase during isotropic and polarized growth of Mucor circinelloides. Can. J. Microbiol. 2007, 53, 168–176. [Google Scholar] [CrossRef] [PubMed]
  28. Wolff, A.M.; Appel, K.F.; Petersen, J.B.; Poulsen, U.; Arnau, J. Identification and analysis of genes involved in the control of dimorphism in Mucor circinelloides (syn. racemosus). FEMS Yeast Res. 2002, 2, 203–213. [Google Scholar] [CrossRef]
  29. Lübbehüsen, T.L.; Nielsen, J.; McIntyre, M. Morphology and physiology of the dimorphic fungus Mucor circinelloides (syn. M. racemosus) during anaerobic growth. Mycol. Res. 2003, 107, 223–230. [Google Scholar] [CrossRef] [Green Version]
  30. Brunke, S.; Mogavero, S.; Kasper, L.; Hube, B. Virulence factors in fungal pathogens of man. Curr. Opin. Microbiol. 2016, 32, 89–95. [Google Scholar] [CrossRef]
  31. Benito, E.P.; Díaz-Mínguez, J.M.; Iturriaga, E.A.; Campuzano, V.; Eslava, A.P. Cloning and sequence analysis of the Mucor circinelloides pyrG gene encoding orotidine-5’-monophosphate decarboxylase: Use of pyrG for homologous transformation. Gene 1992, 116, 59–67. [Google Scholar] [CrossRef]
  32. Corrochano, L.M.; Kuo, A.; Marcet-Houben, M.; Polaino, S.; Salamov, A.; Villalobos-Escobedo, J.M.; Grimwood, J.; Álvarez, M.I.; Avalos, J.; Bauer, D.; et al. Expansion of signal transduction pathways in fungi by extensive genome duplication. Curr. Biol. 2016, 26, 1577–1584. [Google Scholar] [CrossRef] [Green Version]
  33. Wood, V.; Gwilliam, R.; Rajandream, M.-A.; Lyne, M.; Lyne, R.; Stewart, A.; Sgouros, J.; Peat, N.; Hayles, J.; Baker, S.; et al. The genome sequence of Schizosaccharomyces pombe. Nature 2002, 415, 871–880. [Google Scholar] [CrossRef] [Green Version]
  34. Goffeau, A.; Barrell, B.G.; Bussey, H.; Davis, R.W.; Dujon, B.; Feldmann, H.; Galibert, F.; Hoheisel, J.D.; Jacq, C.; Johnston, M.; et al. Life with 6000 genes. Science 1996, 274, 546–567. [Google Scholar] [CrossRef] [Green Version]
  35. Jones, T.; Federspiel, N.A.; Chibana, H.; Dungan, J.; Kalman, S.; Magee, B.B.; Newport, G.; Thorstenson, Y.R.; Agabian, N.; Magee, P.T.; et al. The diploid genome sequence of Candida albicans. Proc. Natl. Acad. Sci. USA 2004, 101, 7329–7334. [Google Scholar] [CrossRef] [Green Version]
  36. Janbon, G.; Ormerod, K.L.; Paulet, D.; Byrnes III, E.J.; Yadav, V.; Chatterjee, G.; Mullapudi, N.; Hon, C.-C.; Billmyre, R.B.; Brunel, F.; et al. Analysis of the genome and transcriptome of Cryptococcus neoformans var. grubii reveals complex RNA expression and microevolution leading to virulence attenuation. PLoS Genet. 2014, 10, e1004261. [Google Scholar] [CrossRef]
  37. Lee, S.C.; Billmyre, R.B.; Li, A.; Carson, S.; Sykes, S.M.; Huh, E.Y.; Mieczkowski, P.; Ko, D.C.; Cuomo, C.A.; Heitman, J. Analysis of a food-borne fungal pathogen outbreak: Virulence and genome of a Mucor circinelloides isolate from yogurt. mBio 2014, 5, e01390-14. [Google Scholar] [CrossRef] [Green Version]
  38. Galagan, J.E.; Calvo, S.E.; Borkovich, K.A.; Selker, E.U.; Read, N.D.; Jaffe, D.; FitzHugh, W.; Ma, L.-J.; Smirnov, S.; Purcell, S.; et al. The genome sequence of the filamentous fungus Neurospora crassa. Nature 2003, 422, 859–868. [Google Scholar] [CrossRef]
  39. Nierman, W.C.; Pain, A.; Anderson, M.J.; Wortman, J.R.; Kim, H.S.; Arroyo, J.; Berriman, M.; Abe, K.; Archer, D.B.; Bermejo, C.; et al. Genomic sequence of the pathogenic and allergenic filamentous fungus Aspergillus fumigatus. Nature 2005, 438, 1151–1156. [Google Scholar] [CrossRef]
  40. Galagan, J.E.; Calvo, S.E.; Cuomo, C.; Ma, L.-J.; Wortman, J.R.; Batzoglou, S.; Lee, S.-I.; Baştürkmen, M.; Spevak, C.C.; Clutterbuck, J.; et al. Sequencing of Aspergillus nidulans and comparative analysis with A. fumigatus and A. oryzae. Nature 2005, 438, 1105–1115. [Google Scholar] [CrossRef]
  41. Vandeputte, P.; Ghamrawi, S.; Rechenmann, M.; Iltis, A.; Giraud, S.; Fleury, M.; Thornton, C.; Delhaès, L.; Meyer, W.; Papon, N.; et al. Draft genome sequence of the pathogenic fungus Scedosporium apiospermum. Genome Announc. 2014, 2, e00988-14. [Google Scholar] [CrossRef] [Green Version]
  42. Haitjema, C.H.; Gilmore, S.P.; Henske, J.K.; Solomon, K.V.; de Groot, R.; Kuo, A.; Mondo, S.J.; Salamov, A.A.; LaButti, K.; Zhao, Z.; et al. A parts list for fungal cellulosomes revealed by comparative genomics. Nat. Microbiol. 2017, 2, 17087. [Google Scholar] [CrossRef] [Green Version]
  43. Hauser, M.; Steinegger, M.; Söding, J. MMseqs software suite for fast and deep clustering and searching of large protein sequence sets. Bioinformatics 2016, 32, 1323–1330. [Google Scholar] [CrossRef] [Green Version]
  44. van Dongen, S.M. Graph Clustering by Flow Simulation. Ph.D. Thesis, University of Utrecht, Utrecht, The Netherlands, 2000. [Google Scholar]
  45. Enright, A.J.; Van Dongen, S.; Ouzounis, C.A. An efficient algorithm for large-scale detection of protein families. Nucleic Acids Res. 2002, 30, 1575–1584. [Google Scholar] [CrossRef]
  46. Lax, C.; Pérez-Arques, C.; Navarro-Mendoza, M.I.; Cánovas-Márquez, J.T.; Tahiri, G.; Pérez-Ruiz, J.A.; Osorio-Concepción, M.; Murcia-Flores, L.; Navarro, E.; Garre, V.; et al. Genes, pathways, and mechanisms involved in the virulence of Mucorales. Genes 2020, 11, 317. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Kubo, H.; Mihara, H. cAMP promotes hyphal branching in Mucor globosus. Mycoscience 2007, 48, 187–189. [Google Scholar] [CrossRef]
  48. Lecointe, K.; Cornu, M.; Leroy, J.; Coulon, P.; Sendid, B. Polysaccharides cell wall architecture of Mucorales. Front. Microbiol. 2019, 10, 469. [Google Scholar] [CrossRef] [PubMed]
  49. Garcia-Rubio, R.; de Oliveira, H.C.; Rivera, J.; Trevijano-Contador, N. The fungal cell wall: Candida, Cryptococcus, and Aspergillus species. Front. Microbiol. 2020, 10, 2993. [Google Scholar] [CrossRef]
  50. Brown, H.E.; Esher, S.K.; Alspaugh, J.A. Chitin: A “hidden figure” in the fungal cell wall. Curr. Top. Microbiol. Immunol. 2020, 425, 83–111. [Google Scholar] [CrossRef]
  51. Doering, T.L. How sweet it is! Cell wall biogenesis and polysaccharide capsule formation in Cryptococcus neoformans. Annu. Rev. Microbiol. 2009, 63, 223–247. [Google Scholar] [CrossRef] [Green Version]
  52. Bartnicki-Garcia, S.; Nickerson, W.J. Isolation, composition, and structure of cell walls of filamentous and yeast-like forms of Mucor rouxii. Biochim. Biophys. Acta 1962, 58, 102–119. [Google Scholar] [CrossRef]
  53. Dow, J.M.; Olona, P.M.; Villa, V.D. Glucuronosyl transferase from the dimorphic fungus Mucor rouxii. Exp. Mycol. 1982, 6, 329–334. [Google Scholar] [CrossRef]
  54. Dow, J.M.; Darnall, D.W.; Villa, V.D. Two distinct classes of polyuronide from the cell walls of a dimorphic fungus, Mucor rouxii. J. Bacteriol. 1983, 155, 1088–1093. [Google Scholar] [CrossRef] [Green Version]
  55. Bartnicki-Garcia, S.; Lindberg, B. Partial characterization of mucoran: The glucuronomannan component. Carbohydr. Res. 1972, 23, 75–85. [Google Scholar] [CrossRef]
  56. López-Fernández, L.; Sanchis, M.; Navarro-Rodríguez, P.; Nicolás, F.E.; Silva-Franco, F.; Guarro, J.; Garre, V.; Navarro-Mendoza, M.I.; Pérez-Arques, C.; Capilla, J. Understanding Mucor circinelloides pathogenesis by comparative genomics and phenotypical studies. Virulence 2018, 9, 707–720. [Google Scholar] [CrossRef] [Green Version]
  57. López-Matas, M.A.; Eslava, A.P.; Díaz-Mínguez, J.M. Mcchs1, a member of a chitin synthase gene family in Mucor circinelloides, is differentially expressed during dimorphism. Curr. Microbiol. 2000, 40, 169–175. [Google Scholar] [CrossRef]
  58. de Groot, P.W.J.; Brandt, B.W.; Horiuchi, H.; Ram, A.F.J.; de Koster, C.G.; Klis, F.M. Comprehensive genomic analysis of cell wall genes in Aspergillus nidulans. Fungal Genet. Biol. 2009, 46 (Suppl. S1), S72–S81. [Google Scholar] [CrossRef] [Green Version]
  59. Rangel-Porras, R.A.; Meza-Carmen, V.; Martinez-Cadena, G.; Torres-Guzmán, J.C.; González-Hernández, G.A.; Arnau, J.; Gutiérrez-Corona, J.F. Molecular analysis of an NAD-dependent alcohol dehydrogenase from the zygomycete Mucor circinelloides. Mol. Genet. Genom. 2005, 274, 354–363. [Google Scholar] [CrossRef]
  60. Fleck, C.B.; Brock, M. Characterization of an acyl-CoA: Carboxylate CoA-transferase from Aspergillus nidulans involved in propionyl-CoA detoxification. Mol. Microbiol. 2008, 68, 642–656. [Google Scholar] [CrossRef]
  61. Arikawa, Y.; Enomoto, K.; Muratsubaki, H.; Okazaki, M. Soluble fumarate reductase isoenzymes from Saccharomyces cerevisiae are required for anaerobic growth. FEMS Microbiol. Lett. 1998, 165, 111–116. [Google Scholar] [CrossRef] [Green Version]
  62. Navarro-Mendoza, M.I.; Pérez-Arques, C.; Murcia, L.; Martínez-García, P.; Lax, C.; Sanchis, M.; Capilla, J.; Nicolás, F.E.; Garre, V. Components of a new gene family of ferroxidases involved in virulence are functionally specialized in fungal dimorphism. Sci. Rep. 2018, 8, 7660. [Google Scholar] [CrossRef] [Green Version]
  63. Ocampo, J.; McCormack, B.; Navarro, E.; Moreno, S.; Garre, V.; Rossi, S. Protein kinase A regulatory subunit isoforms regulate growth and differentiation in Mucor circinelloides: Essential role of PKAR4. Eukaryot. Cell 2012, 11, 989–1002. [Google Scholar] [CrossRef] [Green Version]
  64. Moriwaki-Takano, M.; Iwakura, R.; Hoshino, K. Dimorphic mechanism on cAMP mediated signal pathway in Mucor circinelloides. Appl. Biochem. Biotechnol. 2021, 193, 1252–1265. [Google Scholar] [CrossRef]
  65. Fernández Núñez, L.; Ocampo, J.; Gottlieb, A.M.; Rossi, S.; Moreno, S. Multiple isoforms for the catalytic subunit of PKA in the basal fungal lineage Mucor circinelloides. Fungal Biol. 2016, 120, 1493–1508. [Google Scholar] [CrossRef]
Jof 08 00404 g001 550
Figure 1. Differentially up- (A) and downregulated (B) genes of M. lusitanicus grown anaerobically compared to those after aerobic culturing of the fungus and their distribution among the main functional categories of the genes.
Figure 1. Differentially up- (A) and downregulated (B) genes of M. lusitanicus grown anaerobically compared to those after aerobic culturing of the fungus and their distribution among the main functional categories of the genes.
Jof 08 00404 g001
Jof 08 00404 g002 550
Figure 2. Proportion of the up- and downregulated transporter genes of M. lusitanicus grown anaerobically compared to those after aerobic culturing of the fungus.
Figure 2. Proportion of the up- and downregulated transporter genes of M. lusitanicus grown anaerobically compared to those after aerobic culturing of the fungus.
Jof 08 00404 g002
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Homa, M.; Ibragimova, S.; Szebenyi, C.; Nagy, G.; Zsindely, N.; Bodai, L.; Vágvölgyi, C.; Nagy, G.; Papp, T. Differential Gene Expression of Mucor lusitanicus under Aerobic and Anaerobic Conditions. J. Fungi 2022, 8, 404. https://doi.org/10.3390/jof8040404

AMA Style

Homa M, Ibragimova S, Szebenyi C, Nagy G, Zsindely N, Bodai L, Vágvölgyi C, Nagy G, Papp T. Differential Gene Expression of Mucor lusitanicus under Aerobic and Anaerobic Conditions. Journal of Fungi. 2022; 8(4):404. https://doi.org/10.3390/jof8040404

Chicago/Turabian Style

Homa, Mónika, Sandugash Ibragimova, Csilla Szebenyi, Gábor Nagy, Nóra Zsindely, László Bodai, Csaba Vágvölgyi, Gábor Nagy, and Tamás Papp. 2022. "Differential Gene Expression of Mucor lusitanicus under Aerobic and Anaerobic Conditions" Journal of Fungi 8, no. 4: 404. https://doi.org/10.3390/jof8040404

APA Style

Homa, M., Ibragimova, S., Szebenyi, C., Nagy, G., Zsindely, N., Bodai, L., Vágvölgyi, C., Nagy, G., & Papp, T. (2022). Differential Gene Expression of Mucor lusitanicus under Aerobic and Anaerobic Conditions. Journal of Fungi, 8(4), 404. https://doi.org/10.3390/jof8040404

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