NapA Mediates a Redox Regulation of the Antioxidant Response, Carbon Utilization and Development in Aspergillus nidulans

The redox-regulated transcription factors (TFs) of the bZIP AP1 family, such as yeast Yap1 and fission yeast Pap1, are activated by peroxiredoxin proteins (Prxs) to regulate the antioxidant response. Previously, Aspergillus nidulans mutants lacking the Yap1 ortholog NapA have been characterized as sensitive to H2O2 and menadione. Here we study NapA roles in relation to TFs SrrA and AtfA, also involved in oxidant detoxification, showing that these TFs play different roles in oxidative stress resistance, catalase gene regulation and development, during A. nidulans life cycle. We also uncover novel NapA roles in repression of sexual development, normal conidiation, conidial mRNA accumulation, and carbon utilization. The phenotypic characterization of ΔgpxA, ΔtpxA, and ΔtpxB single, double and triple peroxiredoxin mutants in wild type or ΔnapA backgrounds shows that none of these Prxs is required for NapA function in H2O2 and menadione resistance. However, these Prxs participate in a minor NapA-independent H2O2 resistance pathway and NapA and TpxA appear to regulate conidiation along the same route. Using transcriptomic analysis we show that during conidial development NapA-dependent gene expression pattern is different from canonical oxidative stress patterns. In the course of conidiation, NapA is required for regulation of at least 214 genes, including ethanol utilization genes alcR, alcA and aldA, and large sets of genes encoding proteins involved in transcriptional regulation, drug detoxification, carbohydrate utilization and secondary metabolism, comprising multiple oxidoreductases, membrane transporters and hydrolases. In agreement with this, ΔnapA mutants fail to grow or grow very poorly in ethanol, arabinose or fructose as sole carbon sources. Moreover, we show that NapA nuclear localization is induced not only by oxidative stress but also by growth in ethanol and by carbon starvation. Together with our previous work, these results show that SakA-AtfA, SrrA and NapA oxidative stress-sensing pathways regulate essential aspects of spore physiology (i.e., cell cycle arrest, dormancy, drug production and detoxification, and carbohydrate utilization).


INTRODUCTION
The proposed role of reactive oxygen species (ROS) as essential cell differentiation signals (Hansberg and Aguirre, 1990;Aguirre et al., 2005) led us to study the mechanisms by which eukaryotic cells detoxify ROS, using the model fungus Aspergillus nidulans. To perceive and transmit oxidative stress signals, fungi utilize phosphorelay systems connected to MAP kinases specialized in transducing stress signals or SAPKs. Schizosaccharomyces pombe paradigmatic SAPK Sty1/Spc1 has been characterized as a MAPK involved in cell-cycle control (Shiozaki and Russell, 1995) that is activated by osmotic (Millar et al., 1995;Degols et al., 1996), oxidative (Degols et al., 1996), heat shock (Nguyen and Shiozaki, 1999), nitrogen limitation , and UV light (Degols and Russell, 1997) stress. As indicated in Figure S1 and Table S1, the phosphorelay system linked to Sty1/Spc1 is composed by histidine kinases (HK) Mak1, Mak2, and Mak3 (Buck et al., 2001), the phosphotransfer protein (HPt) Mpr1 and the response regulator (RR) Mcs4. Sty1/Spc1 in turn regulates transcription factor Atf1. Despite the architecture similarity to Saccharomyces cerevisiae Sln1-Ypd1-Ssk1-Hog1 system (de Nadal et al., 2011), S. pombe phosphorelay transmits oxidative, not osmotic stress signals (Nguyen et al., 2000). A second phosphorelay component, the transcription factor Prr1 is also required for oxidative stress responses, independently of Sty1/Spc1 (Quinn et al., 2011). In addition to Atf1 and Prr1, transcription factor Pap1, a homolog of S. cerevisiae Yap1 (Moye- Rowley et al., 1989), is critical for the antioxidant response in this fungus. The oxidation signal is perceived by different peroxiredoxins or Prxs and then transmitted to Pap1 or Yap1, which once oxidized accumulate in the nucleus to regulate the expression of multiple genes involved in the antioxidant response.
All peroxiredoxins belong to a conserved family of peroxidases that reduce peroxide and contain a conserved "peroxidatic" cysteine. Peroxides oxidize this Cys to sulphenic acid, which then reacts with another "resolving" Cys to form a disulfide bond, subsequently reduced by a suitable electron donor to complete a catalytic cycle. Prxs are classified into 2-Cys, atypical 2-Cys and 1-Cys families. 2-Cys are homodimeric and contain peroxidatic and resolving Cys residues in the same subunit. However, the disulfide bond is formed between two different subunits. In atypical 2-Cys an intermolecular disulfide is formed within the same subunit. 1-Cys Prxs form a disulfide with a resolving Cys present in other proteins or small thiol molecules (Rhee, 2016). Until now, typical 2-Cys Prxs have not been found in filamentous fungi. The role of S. cerevisiae peroxiredoxin Gpx3 in Yap1 activation, which also requires Yap1-binding protein Ybp1, was the first description of Prx function in H 2 O 2 sensing (Delaunay et al., 2002). However, under certain conditions peroxiredoxin Tsa1 can also mediate Yap1 activation by H 2 O 2 (Tachibana et al., 2009). In S. pombe, the 2-Cys peroxiredoxin Tpx1 transmits the redox signal to Pap1 (Vivancos et al., 2004(Vivancos et al., , 2005. A. nidulans contains 15 HKs and the function of most of them is unknown. Genetic evidence indicates that HK NikA transmits osmostress and fungicide signals to (HPt) YpdA and to SrkA RR, which is coupled to the SAPK SakA/HogA (Han and Prade, 2002;Kawasaki et al., 2002), as well as to the SAPKindependent RR SrrA Vargas-Perez et al., 2007). Upstream MAPKK PbsB and MAPKKK SskB regulate SakA (Furukawa et al., 2005), which is able to replace Sty1/Spc1 functions in S. pombe, and in A. nidulans is phosphorylated in response to multiple types of stress, including osmotic, oxidative (Kawasaki et al., 2002), nutrient starvation (Lara-Rojas et al., 2011) and hypoxia (Sánchez and Aguirre, unpublished). Stressactivated SakA translocates to the nucleus, where it interacts with transcription factor AtfA, required for induction of multiple genes and both, sakA and atfA mutants are sensitive to oxidative stress (Lara-Rojas et al., 2011). Additionally, SakA and AtfA are required for osmotic-induced gene expression (Hagiwara et al., 2009).
TF SrrA is also needed for oxidative stress resistance (Vargas-Perez et al., 2007) and both, SakA and SrrA play important roles during development. SakA represses sexual development and is activated during asexual development (Kawasaki et al., 2002).
sakA intact conidia progressively lose their viability and this is consistent with the fact that phosphorylated SakA accumulates in asexual spores (conidia) in an AtfA-dependent manner, and its dephosphorylation is necessary for germination to take place (Lara-Rojas et al., 2011). Likewise, srrA mutants show severely decreased asexual sporulation and produce conidia that very rapidly lose their viability (Vargas-Perez et al., 2007). In addition to transcription factors (TFs) AtfA and SrrA, the Yap1/Pap1 functional homolog NapA has been shown to be required for resistance to H 2 O 2 in A. nidulans (Asano et al., 2007). Unrelated protein AN8863, putatively involved in nucleosome assembly, was later also referred to as NapA (Araújo-Bazan et al., 2008). Here we keep using NapA to name the A. nidulans Yap1/Pap1 homolog (AN7513) because it has been used this way in other publications (Lessing et al., 2007;Thön et al., 2010), and because the name "ap" preceded by the first letter of the species name (i.e., nidulans apA) has been widely used in many other filamentous fungi, where the role of Yap1/Pap1 homologs in oxidative stress resistance has been demonstrated (Lessing et al., 2007;Qiao et al., 2008;Temme and Tudzynski, 2009;Tian et al., 2011;Cartwright and Scott, 2013). Notably, in several plant pathogens Yap1/Pap1 homologs are involved not only in regulation of the antioxidant response but also in plant virulence (Molina and Kahmann, 2007;Guo et al., 2011;Huang et al., 2011).
Here we compared the relative contribution of (TFs) AtfA, SrrA and NapA to the antioxidant response and development in A. nidulans and uncovered novel NapA roles in regulation of sexual and asexual development, carbon utilization and gene regulation during asexual sporulation. strains in a veA + background were obtained from sexual crosses with strain FGSCA4. The presence of wild type veA allele was confirmed by PCR using genomic DNA from selected progeny and the primers veAforward and veAreverse, as reported (Han et al., 2010). Menadione was filter sterilized and like H 2 O 2 , added to agar medium at 50 • C before solidification. H 2 O 2containing plates were used the day they were prepared or stored at 4 • C for no more than 24 h. Since H 2 O 2 can react with medium components, the actual concentration in plates cannot be estimated. To ensure experimental reproducibility, the same batch of H 2 O 2 containing medium was used when comparing different strains. Spore suspensions containing 1 × 10 3 or 1 × 10 4 conidia were used to inoculate plates by dropping the suspension on the center of plates containing different stressors or media. Higher H 2 O 2 resistance is observed at higher spore densities, presumably due to the high catalase A activity levels found in conidia (Navarro et al., 1996;Navarro and Aguirre, 1998). For catalase activity, 30 µg of total protein extracts prepared from conidia or mycelia were separated on native polyacrylamide gels to determine catalase activity as reported (Navarro et al., 1996;Kawasaki et al., 1997). Briefly gels are incubated in 5% methanol with shaking for 5 min and then rinsed with tap water 3 times. After this, the gel is incubated in a 0.03% hydrogen peroxide solution (100 µl of commercial 30% solution in 100 ml of deionized water) for 5 min and rinsed with water. Finally, the gel is incubated in the staining solution until the bands of activity are visible. The staining solution is made by mixing equal volumes of a 2% (w/v) FeCl 3 solution and a 2% (w/v) K 3 Fe(CN) 6 . The different catalases were mapped before using single, double and triple mutants affected in catA, catB, and catC genes (Kawasaki and Aguirre, 2001).
Deletion of napA, gpxA, tpxA, tpxB and alcA Genes, and Tagging of NapA Genomic DNA was used as template to produce the gene-deletion constructs by double joint PCR (Yu et al., 2004). For napA gene (AN7513) replacement construct, the napA ORF was amplified with primers 5 ′ For-napA and 5 ′ Rev-napA (see Table S3). The 3 ′ napA fragment was amplified with primers 3 ′ For-napA and 3 ′ Rev-napA. Aspergillus fumigatus pyrG marker was amplified with primers pyrGforward and pyrGreverse, using plasmid PFNO3 as template (Nayak et al., 2006). The three fragments were purified, mixed and used in a fusion PCR with primers 5 ′ nest-napA and 3 ′ nest-napA. The final 4900 bp napA-AfpyrG-napA cassette was purified and used to transform A. nidulans strain 11035 by electroporation (Sanchez and Aguirre, 1996;Sánchez et al., 1998). Five PyrG + transformants were obtained and analyzed by Southern blot to confirm the elimination of napA. After confirming the proper deletion event ( Figure S2), strain TFL9 was chosen and crossed with strain CLK43 to get rid of the kuA deletion, and progeny strain CFL7 was confirmed by PCR and used in further experiments.
The alcA deletion construct containing the riboB gene, as a selective marker, was generated using primers 5 ′ For-alcA and 5 ′ Rev-alcA for 5 ′ region and 3 ′ For-alcA and 3 ′ Rev-alcA for the 3 ′ region. A. fumigatus riboB marker was amplified with primers 5Ribo and 6Ribo as before. The 4167 bp band obtained with primers NestForalcA and NestRevaclA was used to transform strain 11035 by electroporation. Six transformants contained the expected event ( Figure S4). Transformant TAM20 was crossed with strain CLK43 and progeny strain CAM17 was confirmed by PCR and chosen for further experiments. Sexual crosses generated double, triple and quadruple mutants, as indicated in Table S2.
To generate NapA::GFP C-terminal construct, three PCR products were used. First, 5 ′ and entire napA ORF were amplified with primers GSP1napA and GSP2napA. Second, a 3 ′ napA fragment was amplified with primers GSP3napA and GSP4napA. Third, GFP and A. fumigatus pyrG marker were amplified with primers GFP1napA and GFP2napA as before. Purified fragments were used in a fusion PCR with primers GSP1napA and GSP4napA. The 6625 bp napA-GFP-AfpyrG cassette was used to transform A. nidulans strain A1155 by electroporation. Transformant TFL14 strain was confirmed by PCR and used for further experiments ( Figure S5).

Microscopy
Fluorescence microscopy images were captured in vivo. Spores from NapA::GFP strain were grown for 18 h on coverslips containing liquid glucose-MM at 37 • C with no shaking. After this, samples were shifted to the same medium with or without 2 mM H 2 O 2 , incubated for 0-120 min and observed using a NIKON Eclipse E600 microscope to detect DAPI and GFP fluorescence. Images were captured with a cooled camera Neo Andor sCMOS. For DAPI staining, samples of conidia or mycelia were fixed in methanol/acetone for 10 min, washed in water and stained for 5 min with 0.1 mg/ml of 4 ′ ,6-diamidino-2phenylindole (DAPI).

Transcriptomic Analysis
Conidia from 6-day old cultures from strains CLK43 (WT) and CFL7 ( napA) were collected and washed 3 times with 10 ml of cold water. Excess liquid was removed by centrifugation and conidia were immediately frozen with liquid nitrogen. Total RNA was extracted by cryogenic grinding using the Tissue Lyser (Qiagen) and purified using the RNAeasy Mini kit (Qiagen), following the manufacturer's protocol. RNA integrity number (RIN) for RNA quality was generated using the Agilent 2100 Bioanalyzer System (Agilent technologies). Two independent samples from each strain (Biological replicates) were processed for cDNA synthesis using Illumina's kit TrueSeqV2 and sequenced at the next-generation sequencing core facility at IBT-UNAM. using Illumina's platform. An average of 10168 reads of 72 bp per sample were obtained, representing nearly 96.3% A. nidulans genome lengths. Biological replicates showed a good level of correlation (r_0.966 for WT and r_0.968 for napA). Differential gene expression was inferred based on total mapping counts using the EdgeR package. Genes showing a value of log fold change (LFC) ≥2 and a false discovery rate (FDR) ≤0.05 were considered as differentially expressed (DE). The differential expression analysis was deposited under the GEO identifier GSE94747, as part of the BioProject PRJNA373914. The sequencing raw reads for all experiments were deposited in the SRA database under the SRP099165 identifier.

Real-Time PCR
Total RNA was isolated using the Plant RNA purification kit (Quiagen) according to the manufacturer's instructions and treated with DNase Turbo DNA-free kit (Ambion). The RNA integrity (RIN) and concentration were determinate using the 2200 TapeStation (Agilent Genomics). 1 µg of RNA was used to synthetize cDNA templates for PCR amplification, using SuperScript III reverse transcriptase (Invitrogen) according to the manufacturer's instructions. The expression of each gene was measured using 1:10 cDNA sample dilutions. Specific primers used for RT-PCR (Table S3) were devised to produce cDNA amplicons around 130 pb. To estimate genomic DNA contamination, primers were designed to amplify a larger product, including an intron of each gene, from a genomic template. RT-PCR products were observed in 2% agarose gels.
Triplicates from each sample were performed. RT-PCR was done using Platinum SYBR Green qPCR SuperMix-UDG with ROX kit (Invitrogen). RT-PCR was performed in a StepOne Real-Time PCR System (Applied Biosystems). The program used included an initial UDG incubation for 2 min at 50 • C, followed by a 2 min denaturation step at 95 • C, and 40 amplification cycles at 95 • C for 15 s, followed by 30 s at 60 • C. For relative quantification, we used the Ct comparative method and the data was analyzed with StepOne Software V2.3 (Applied Biosystems), using histone 2B (H2B) gene (AN3469) as reference.

RESULTS
TFs NapA, SrrA, and AtfA Play Differential Roles in A. nidulans Antioxidant Response Previous work has shown that (TFs) SrrA (Vargas-Perez et al., 2007), AtfA (Lara-Rojas et al., 2011) and NapA (Asano et al., 2007) are involved in A. nidulans antioxidant response. To analyse the relative contribution of each TF in this process, we compared the sensitivity of conidia and mycelia from napA, srrA, and atfA mutants to H 2 O 2 and menadione. Results show that conidia from napA mutant were the most sensitive to H 2 O 2 followed by srrA and atfA mutants, while only napA mutants were sensitive to menadione ( Figure 1A). In contrast, the same assay carried out with mycelia showed a slightly higher H 2 O 2 sensitivity for srrA mutants followed by napA, while the atfA mutant was as resistant to H 2 O 2 as the WT strain. Again only the napA mutant was sensitive to menadione ( Figure 1B).
To further dissect NapA, SrrA, and AtfA contribution to the antioxidant response, we used napA, srrA, and atfA single mutants to perform sexual crosses and obtain double and triple mutants, which were confirmed by PCR analysis (not shown). Results in Figure 1A show that conidia from double and triple mutants did not show additive phenotypes and those carrying the napA deletion were as sensitive to H 2 O 2 and menadione as the single napA mutant, except in the case of the srrA napA strain whose H 2 O 2 sensitivity was similar to the one displayed by the srrA mutant ( Figure 1A). In the test performed with mycelia, double and triple mutants carrying the srrA deletion behaved as the single srrA mutant ( Figure 1B). Notably, srrA mutant growth defects (Vargas-Perez et al., 2007) are enhanced by the presence of the atfA deletion, as colonies from srrA atfA mutants show a higher decrease in growth and highly irregular colony borders. These growth and conidiation defects were even more drastic in the mutant lacking the 3 TFs (Figures 1A,B).
These results show that TFs NapA, SrrA, and AtfA play different roles at different stages of A. nidulans life cycle. For H 2 O 2 spore resistance NapA plays a more prominent role than SrrA, which in turn is more important than AtfA. For mycelial H 2 O 2 resistance SrrA is somewhat more critical than NapA, while AtfA plays no role in this process, and only NapA is necessary for menadione resistance in both spores and mycelia. During vegetative growth, neither NapA nor AtfA are individually required for normal radial growth. However, srrA mutant growth defects are enhanced by the deletion of AtfA, indicating that AtfA contributes to normal radial growth. SrrA is almost essential for normal conidiation, NapA is needed for full conidiation (see further) and AtfA is dispensable for the process. Regarding conidial function, SrrA and AtfA are critical for conidial viability, and NapA is necessary for the accumulation of multiple mRNAs in conidia (see further).

NapA, SrrA, and AtfA Regulate Different Catalase Genes
To further understand NapA roles in the antioxidant response, we compared the H 2 O 2 sensitivity of conidia from mutants lacking NapA, the spore-specific catalase CatA (Navarro et al., 1996;Navarro and Aguirre, 1998) or the mycelium inducible catalase CatB (Kawasaki et al., 1997). Notably, napA conidia were much more sensitive to H 2 O 2 than catA conidia, while catB conidia showed only a minor sensitivity at 4 mM H 2 O 2 ( Figure S6). As shown in Figure 2A, the presence of CatA activity was completely dependent on AtfA, while being independent of NapA. In contrast, H 2 O 2 induction of mycelial catalase CatB was largely dependent on both NapA ( Figure 2B) and SrrA (Vargas-Perez et al., 2007). Since CatB is also highly induced during the stationary phase of growth (Kawasaki and Aguirre, 2001), we asked if NapA, SrrA or AtfA were required for this induction. Unexpectedly, none of these TFs was needed for this process. However, the induction of catalase-peroxidase CatD/CpeA activity (Kawasaki and Aguirre, 2001;Scherer et al., 2002) required of SrrA ( Figure 2C). In summary, these results show that NapA and SrrA are both required for CatB induction by H 2 O 2 , AtfA is required for CatA expression in conidia and none of them is required for CatB induction during the late stationary phase of growth. This and SrrA regulation of CatD/CpeA activity (Kawasaki and Aguirre, 2001) confirm the differential roles that these TFs perform during A. nidulans antioxidant response.

NapA Localizes in Nuclei in Response to Oxidative Stress
The high sensitivity of napA conidia to H 2 O 2 , not explained by a simple lack of spore catalase activity (Figure 2A), suggests that NapA is required for other conidial functions, including those involved in resistance to menadione. It is well known that S. cerevisiae, S. pombe and other fungal NapA homologs show nuclear accumulation when oxidized in response to H 2 O 2 . Peroxiredoxins Gpx3 and Tpx1 perceive the oxidation signal and relay it to Yap1 and Pap1, respectively. To determine NapA localization we introduced a GFP tag at its C-terminus and showed that the H 2 O 2 resistance of the corresponding strain FIGURE 2 | Transcriptional factors NapA, SrrA and AtfA play differential roles in catalase regulation. (A) CatA activity levels are not affected in napA conidia. Conidial protein extracts from strains CLK43 (WT), TFL atfA-04 ( atfA) and CFL7 ( napA) were analyzed for catalase in-gel activity. (B) CatB activity induction by H 2 O 2 requires NapA. Mycelia from the strains CLK43 (WT) and CFL7 ( napA) were grown for 12 h in liquid MM and then H 2 O 2 was added at the indicated times and concentrations and protein extracts were used for catalase activity determination. (C) NapA, AtfA, and SrrA are not required for CatB induction during the stationary phase of growth. Mycelia from the strains CLK43 (WT), CFL7 ( napA), TFL atfA-02 ( atfA), COSsrrA3 ( srrA) and TLK12 ( catB) was grown for 12, 24, and 48 h in liquid MM medium. The experiment was repeated at least 3 times; a representative experiment is shown.
was not affected ( Figure S5), indicating the functionality of this NapA::GFP fusion. Then we analyzed NapA::GFP localization using different H 2 O 2 concentrations. Initially, we found a very low basal signal of NapA::GFP expressed from the napA promoter, which increased in nuclei after 30 min of treatment with different H 2 O 2 concentrations (not shown). Under these conditions 2 mM was the minimum concentration needed to induce NapA::GFP nuclear accumulation and therefore we used this concentration for a time-course analysis. As seen in Figure 3A, NapA::GFP starts to show nuclear localization after 20 min of treatment, which gradually increases up to 120 min. When after this time the H 2 O 2 treatment was stop by transferring mycelia to a medium lacking H 2 O 2 , it took 6 h to observe that NapA::GFP was no longer observed in nuclei, showing instead a cytoplasm localization. This indicates that under these conditions NapA-mediated adaptation to oxidative stress is a relatively slow process. Since napA mutants are sensitive to menadione, we also determined if menadione was able to induce NapA::GFP nuclear localization. Indeed, 10 µM menadione induced NapA nuclear accumulation after a 90 min treatment ( Figure 3B).

NapA Function in the Antioxidant
Response Is Independent of Peroxiredoxins GpxA, TpxA, and TpxB As indicated before, S. cerevisiae peroxiredoxin Gpx3 and S. pombe Tpx1 relay the oxidation signal to Yap1 and Pap1, respectively. To determine if homologous peroxiredoxins were involved in NapA function, we searched the A. nidulans genome and found genes AN2846, AN10223, and AN3973 (Oh et al., 2010), which according to their closest homolog in S. cerevisiae are named as gpxA (Thön et al., 2010), tpxA and tpxB, respectively. We generated strains in which one of these peroxiredoxin genes was deleted and by sexual crosses generated triple as well as quadruple mutants containing the napA deletion. In sharp contrast to napA mutants, mutants lacking either a single or all 3 peroxiredoxins were not sensitive to H 2 O 2 or menadione. However, the simultaneous inactivation of the 3 peroxiredoxins resulted in an enhancement of the sensitivity to H 2 O 2 caused by the inactivation of NapA (Figure 4). This suggests that these 3 proteins might play a partially redundant minor function in H 2 O 2 resistance, different from NapA.

NapA Is Involved in Developmental Regulation
napA mutants did not show any obvious defects except that they seemed to develop more and paler cleistothecia than the WT strain. To examine this in more detail, we induced sexual development in confluent cultures and determined the number of cleistothecia as reported before (Kawasaki et al., 2002). Indeed, a napA mutant produced about 3 times more cleistothecia than the WT strain ( Figure 5A). Notably, napA young cleistothecia were white, while older cleistothecia developed FIGURE 3 | NapA is induced by oxidative stress and gradually accumulates in nuclei in response to oxidative stress. (A) NapA nuclear localization increases during incubation with H 2 O 2 . Conidia from strain CAM20 (NapA::GFP) were grown for 18 h in MM and then exposed to 2 mM H 2 O 2 for the indicated times (0-120 min), observed in vivo and photographed every 10 min using Epifluorescence microscopy. Lower panel shows NapA::GFP and nuclei (DAPI) signal in mycelia treated with H 2 O 2 for 120 min, fixed and photographed. (B) Menadione also induces nuclear localization of NapA. Conidia from strain CAM20 (NapA::GFP) were grown for 18 h in MM and then exposed to 10 µM menadione for the indicated times (0-90 min). Larger square areas in each picture show enlargements of the areas indicated by smaller squares. Bars = 10 µm. some pigmentation ( Figure 5B). When observed under the microscope, it became clear that napA cleistothecia are not pigmented ( Figure 5C, right panel) but produce pigmented ascospores, which are viable (not shown). A. nidulans veA gene encodes a member of a fungal protein family required for cleistothecium formation (Kim et al., 2002) and regulation of secondary metabolism (Bayram et al., 2008). Since our laboratory strains contain a veA1 allele that causes higher production of conidia and lower numbers of cleistothecia (Käfer, 1965), we also evaluated napA sexual development phenotypes in the presence of a wild-type veA allele. As expected, results show higher numbers of cleistothecia in veA + strain FGSC4 and even higher in napA veA + strains. Again, young cleistothecia were not pigmented, while older cleistothecia were non-pigmented but produced pigmented ascospores (not shown). Likewise, the presence of wild type veA gene did not modify napA oxidative stress sensitivity. These results indicate that NapA represses sexual development and is needed for the expression of genes, possibly a polyketide synthetase gene, associated with the synthesis of cleistothecial melanin. When conidiation levels were determined, we found that a napA mutant produced half of the conidia of the WT strain (Figure 6), indicating that NapA plays a role in the conidiation process. We also examined the conidiation levels of single, double and triple peroxiredoxin mutants. As shown in Figure 6, the inactivation of gpxA and tpxB had a very minor impact in conidiation. In contrast, the inactivation of tpxA reduced conidiation to levels similar to those observed in gpxA tpxA tpxB and napA mutants. Moreover, a quadruple gpxA tpxA tpxB napA mutant showed conidiation numbers slightly higher than those seen in single tpxA and napA mutants, indicating that TpxA and NapA functions in conidiation are not additive and therefore suggesting that these two proteins work in the same conidiation pathway. Notably, tpxA orthologs are induced by H 2 O 2 and regulated by NapA homologs AfYap1 and Bap1 in A. fumigatus (Lessing et al., 2007) and Botrytis cinerea (Temme and Tudzynski, 2009), respectively. These results suggest that TpxA plays NapA independent and dependent roles in oxidative stress resistance and conidiation, respectively.

NapA Is Required for Gene Regulation in Conidia
To address the role of NapA in the production of fully functional conidia, we decided to carry out a transcriptomic analysis. For this, we isolated total RNA from napA and WT intact conidia from 6-day old colonies. Two independent samples from each strain (Biological replicates) were processed for cDNA synthesis and DNA sequencing using Illumina's platform. An average of 10168 reads of 72 bp per sample representing nearly 96.3% A. nidulans genome lengths per sample were obtained. The two biological replicates showed a good level of correlation (r_0.966 for WT and r_0.968 for napA) and principal component analysis showed a clear separation between napA and WT samples. Results show changes in the expression of 284 genes with a logFC differences higher than 1.5 but lower than 2 (Table S4) and 214 genes with logFC changes higher than 2 (Tables 1, 2). The following analysis will be based mostly on these 214 genes, from which 13 genes were up regulated and 201 were down regulated in napA conidia, as compared to WT conidia. Using gene ontology (GO) enrichment analysis, Aspergillus genome database AspGD (Cerqueira et al., 2014) and manual annotations, we grouped the genes regulated by NapA as shown in Tables 1, 2.
Our results show that the number of genes negatively regulated by NapA is small (13), compared with the number of genes positively regulated by NapA (201). Among genes most highly repressed by NapA, there is the putative HET (heterokaryon incompatibility) containing domain protein AN7834, the cysteine-rich secreted protein AN7836 and the putative ABC transporter AN7839, which seems to represent the only NapA-regulated ABC-type transporter. In addition, this group includes genes for 3 additional putative transporters (MSF type), 2 oxidoreductases, 1 hydrolase, 1 acetyltransferase, 1 protein with a nucleoside phosphorylase domain and 2 proteins with no identifiable domains ( Figure 7A, Table 1).
Twenty eight genes encode proteins that have recognizable domains but did not belong to the previous categories. Among these, AN2587 shows the highest LogFC value (8.7), and encodes a putative membrane protein with a cysteine-rich CFEM domain, recently shown to be involved in Fe 3+ heme acquisition (Nasser et al., 2016). This and the fact NapA also regulates ferric reductase AN7981, putative heme binding protein AN2778 and NRPS SidC (AN0607, siderophore biosynthesis; Table S4) genes supports a NapA role in iron acquisition. Also notable is NapA regulation of several genes involved in the biosynthesis of secondary metabolites (Tables 2 and Table S4), such as sterigmatocystin (stcU/AN7806), penicillin (ipnA/AN2622), monodictyphenone and prenyl xanthones (Bok et al., 2009;Sanchez et al., 2011;Andersen et al., 2013) (xptB/AN12402, mdpB/AN10049, mdpC/AN0146 and mdpD/AN0147), all with clear antimicrobial activity. Other genes are involved in the biosynthesis of asperfuranone (afoE/AN1034) or unknown metabolites (PKS AN11191). In contrast, NapA has been reported to repress  the production of secondary metabolites in A. parasiticus (Reverberi et al., 2008) and A. nidulans (Yin et al., 2013) during growing conditions, suggesting that NapA regulates secondary metabolism in opposite ways during growth and conidiation. Finally, NapA regulates at least 54 genes encoding hypothetical proteins with no recognizable domains ( Figure 7B, Table 2).
In summary, our transcriptomic results support a model in which during conidiation NapA regulates functions associated with the successful germination of conidia in natural environments. In the case of A. nidulans saprophytic life style this would include a large group of enzymes to degrade plant complex carbohydrates, proteins to transport the corresponding derived sugars and enzymes needed to transform them into acetyl-CoA. NapA also regulates the production of secondary metabolites, such as penicillin, which would prevent growth from competing organisms, as well as enzymes and transporters to detoxify drugs produced by competitors. In the same line, NapA regulates enzymes and proteins involved in iron scavenging, such as those involved in siderophore biosynthesis, iron reduction and transport, and heme acquisition. As such functions are also important during pathogenic interactions, we propose that these and not NapA ROS-detoxification roles, might be relevant for NapA critical virulence roles (Molina and Kahmann, 2007;Guo et al., 2011;Huang et al., 2011), specially when conidia mediate the infection process. It called our attention that the mRNA levels of genes alcR (AN8978 ; Table S4), alcA (AN8979), alcU (AN8982 ; Table S4) and aldA (AN0554), all members of the well-studied ethanol system (Fillinger and Felenbok, 1996), were clearly reduced in napA conidia. The ethanol system is composed by genes alcP, alcR, alcO, alcA, alcM, alcS, and alcU, all clustered in chromosome VII and aldA located in chromosome VIII. alc genes share the same regulation, being strongly induced by the physiological inducer acetaldehyde and the non-physiological inducer 2-butanone (Flipphi et al., 2002). They are also subject to strict control by the transcriptional activator AlcR and are repressed by glucose via the CreA repressor (Fillinger and Felenbok, 1996). However, only AlcR, alcohol dehydrogenase I AlcA, and aldehyde dehydrogenase AldA are required for ethanol utilization, and the specific function of the other alc genes is unknown. Yet, AlcS is a membrane protein with homology to acetate transporters that is nevertheless dispensable for growth on ethanol, acetaldehyde or acetate (Flipphi et al., 2006). Therefore, we decided to confirm alcR, alcA and aldA expression results by using qPCR. Results in Figure 7C show that indeed, the accumulation of the corresponding mRNAs in conidia depends on NapA.
To test the physiological significance of NapA-mediated alc gene regulation, we examined the ability of napA mutants to grow on ethanol as sole carbon source. As control we included a creA d 204 mutant, which shows derepressed expression of alcA in the presence of glucose (Shroff et al., 1996). Since AlcA transforms allyl alcohol (AA) into the toxic compound acrolein, this creA mutant is unable to grow on glucose plus AA (Figures 8A,B). Notably, napA conidia and mycelia were unable to grow on ethanol, supporting NapA function in proper alcA expression. However, the napA mutant was also unable to grow on glucose plus AA (Figures 8A,B). AA/acrolein causes oxidative stress and induces Yap1 activation in S. cerevisiae (Kwolek-Mirek et al., 2009;Golla et al., 2015) and therefore it might result particularly toxic for the napA mutant. To test this, we generated alcA and napA alcA mutants and examined them for ethanol utilization and AA resistance. As seen in Figure 8C, the alcA mutant is not sensitive to H 2 O 2 , is unable to grow on ethanol and is fully resistant to AA. In contrast, a napA alcA mutant is sensitive to H 2 O 2 , unable to grow on ethanol and more sensitive to AA than the alcA mutant, supporting the idea that napA mutants are per se more sensitive to AA. The failure to grow in 1% ethanol was not remediated by the presence of the antioxidant N-acetylcysteine (5 mM), suggesting that such failure is related to the reduced expression of the alc genes. However, it cannot be excluded that this might be also related to an increased ethanol sensitivity of the napA mutants. In summary, these results show that NapA is required for expression of the alc regulon during asexual development and for growth on ethanol, adding a new level of regulation to this well know pathway.
These results prompted us to determine if NapA and TFs SrrA and AtfA were also required for conidia germination and growth on ethanol or other alternative carbon sources. As seen in Figure 9, NapA was also required for arabinose and fructose utilization, while SrrA was partially required for glycerol utilization and AtfA was dispensable for the utilization of all carbon sources tested. These results indicate that NapA plays an

NapA Is Localized in Nuclei during Growth in Ethanol or Glucose Starvation
Having found that NapA is necessary for proper growth in arabinose and fructose, we explored NapA localization during growth in poor carbon sources, as well as during glucose starvation. Results in Figure 10 show that as it occurs during growth in glucose, NapA expression is very low during growth in arabinose and fructose as sole carbon sources. In contrast, NapA::GFP expression was induced during growth in ethanol, as well as during glucose starvation, and NapA::GFP nuclear localization was detected between 60 and 120 min after shifting to ethanol or glucose lacking media. In addition, we found that a osmotic stress treatment with 1 M NaCl is able to induce NapA nuclear localization under normal glucose growth conditions ( Figure S7), despite the fact that napA mutants are not sensitive to osmotic stress.

DISCUSSION
AtfA, SrrA, and NapA Play Non-redundant and Differential Roles in the Antioxidant Response At least three pathways, mediated by TFs, AtfA, SrrA and NapA, are involved in fungal responses to oxidative and other types of stress and we have studied them here in the same organism. We have reported that four catalases present in A. nidulans display differential regulation during growth, stress and development (Navarro et al., 1996;Navarro and Aguirre, 1998;Kawasaki and Aguirre, 2001). Notably, catalase CatA is found only in asexual and sexual spores, despite the fact that catA mRNA accumulates under many stress conditions Lara-Rojas et al., 2011). Our results here show that TFs, AtfA, SrrA and NapA are essential components in such catalase gene regulation. NapA and SrrA are both required for CatB induction by H 2 O 2 , only AtfA is required for CatA expression in conidia and none of them is required for CatB induction during late stationary phase of growth, while CatD/CpeA activity levels are influenced by SrrA. The phenotypic analysis of oxidative stress resistance using napA, srrA, and atfA single, double and triple mutants also showed that TFs NapA, SrrA, and AtfA play differential roles in oxidative stress resistance.
The functions of Ap1-like TFs in different fungi show both similarities and differences. For example and in contrast to our results, yapA deletion in symbiotic fungus Epichloë festucae causes sensitivity to H 2 O 2 and t-BOOH in conidia but not in mycelia and YapA, instead of AtfA, is required for expression of the spore-specific catalase CatA (Cartwright and Scott, 2013). Neurospora crassa mutants lacking the NapA homolog NcAp-1 were reported as showing no sensitivity to osmotic stress and only a slight sensitivity to H 2 O 2 (Takahashi et al., 2010), although recent work showed that mutants were sensitive to osmotic stress, cadmium and H 2 O 2 (Tian et al., 2011). In Magnaporthe oryzae inactivation of NapA homolog MoAP1 causes only a mild sensitivity to H 2 O 2 (Guo et al., 2011), while Ustilago maydis yap1 null mutants are sensitive to H 2 O 2 and their virulence is significantly reduced (Molina and Kahmann, 2007). A. fumigatus, yap1 null mutants are also sensitive to H 2 O 2 and menadione but not affected in pathogenicity (Lessing et al., 2007).
Like in the case of Yap1 and Pap1, the activation of NapA homologs by H 2 O 2 results in nuclear accumulation in all filamentous fungi where this has been studied, and in U. maydis two Yap1 conserved cysteines  were shown to be crucial for both, nuclear accumulation and functionality (Molina and Kahmann, 2007). An additional level of redox regulation has been proposed in A. nidulans, where the CCAAT-binding complex (CBC) represses napA expression under low ROS levels, while high ROS levels result in oxidation of two cysteines in HapC CBC subunit, and the transcriptional activation of napA (Thön et al., 2010). We found that in addition to H 2 O 2 and menadione, osmotic stress ( Figure S7), glucose starvation stress and growth on ethanol also induced NapA nuclear localization, supporting the idea that different types of stress can all lead to oxidative stress (Hansberg and Aguirre, 1990).

GpxA, TpxA, and TpxB Peroxiredoxin Function in the Antioxidant Response and Conidiation
We analyzed the role of Gpx3 peroxiredoxin homolog GpxA and Tpx1 homologs TpxA and TpxB in NapA antioxidant and developmental functions. We found that none of these peroxiredoxins was required for H 2 O 2 or menadione resistance and therefore are unnecessary for NapA activation. This is consistent with results in E. festucae, where peroxiredoxins Gpx3 (GpxA) and Tpx1 (TpxA) were not needed for YapA H 2 O 2 -induced nuclear accumulation (Cartwright and Scott, 2013). On the contrary, M. oryzae mutants lacking Gpx3 functional homolog MoHYR1 are sensitive to H 2 O 2 , fail to express several genes related to the antioxidant response and show reduced virulence. However, initial data showing that MoYap1 mutants are no affected in pathogenicity indicates that MoHYR1 functions in virulence are not mediated by MoYap1 (Huang et al., 2011). In A. fumigatus AfYap1 also accumulates in the nuclei in response to H 2 O 2 and notably, GpxA homolog AspF3 and TpxA require AfYap1 for its H 2 O 2 -mediated induction (Lessing et al., 2007). These and our results suggest that NapA might be required for H 2 O 2 -induced expression of gpxA, tpxA and tpxB. A minor direct or indirect role for these peroxiredoxins in H 2 O 2 detoxification is indicated by the fact that in a napA background, the simultaneous inactivation of Gpx3, TpxA, and TpxB resulted in increased sensitivity to H 2 O 2 . Our results suggest that TpxA and NapA regulate conidiation through the same pathway. One possibility is that NapA mediates TpxA induction during conidiation, although we did not detect tpxA mRNA in our transcriptomic experiments, and TpxA could in turn regulate other activities needed for full sporulation.

NapA Regulates Development
We reported that the regulated production of ROS is essential for fungal sexual development (Lara-Ortíz et al., 2003;Cano-Domínguez et al., 2008), a finding demonstrated in several fungi (Malagnac et al., 2004;Siegmund et al., 2013;Dirschnabel et al., 2014). In addition to their roles in the antioxidant response, AtfA and SrrA play different roles in gene regulation during development. AtfA mediates SakA roles in sexual development and conidial viability (Kawasaki et al., 2002;Lara-Rojas et al., 2011), while SrrA is required for normal conidiation and conidial viability (Vargas-Perez et al., 2007). The results reported here uncovered a novel developmental role for NapA, by showing that NapA represses sexual development and is needed for full conidiation, as well as for the accumulation of multiple mRNAs in conidia. The requirement for increased ROS levels during sexual development is consistent with the fact that napA mutants show an increased production of fruiting bodies, as NapA is required to maintain low ROS levels and to express genes involved in the biosynthesis of cleistothecial melanin, a well-known antioxidant.
The lower conidiation observed in napA mutants might be related to their lower expression of gmcA gene (AN8547), encoding a putative glucose-methanol-choline oxidoreductase required for early stages of conidiophore development (Etxebeste et al., 2012). The fact that GmcA ortholog AFUA_3G01580 is also induced by H 2 O 2 in an AfYap1-dependent manner in A. fumigatus (Lessing et al., 2007) suggests that gmcA regulation by NapA is conserved, at least in the Aspergilli. Moreover, NapA roles in development might be conserved in fungi. In M. oryzae MoAP1 deletion causes only mild sensitivity to H 2 O 2 but mutants show a drastic reduction in formation of aerial mycelium and conidiation. Notably, the same phenotypes are observed in mutants affected in MoAP1 regulated genes MGG_01230 and MGG_15157, encoding succinic semialdehyde dehydrogenase MoSsadh and acetyltransferase MoAct, respectively (Guo et al., 2011), which however we did not detect as NapA-dependent. In the dimorphic fungus Talaromyces marneffei, yapA mutants are sensitive to H 2 O 2 and menadione, show decreased radial growth, produce conidiophores with fewer phialides and conidia, and conidia show decreased germination rates, while yeast cells fail to undergo binary fission (Dankai et al., 2016). Overall, the developmental roles of TFs long associated only with ROS detoxification support the role of ROS as developmental signals (Hansberg and Aguirre, 1990;Aguirre et al., 2005).
Overall, such relatively conserved H 2 O 2 -induced Ap1mediated gene expression patterns are different from the one we observe during conidiation. Indeed, we do not find genes involved in major NADPH or GSH generation pathways. Instead, we find a large number of genes involved in drug efflux and detoxification, including several genes for enzymes with putative cytochrome P450 activity, as well as other oxidases. Notably, Yap1 (Lee et al., 1999) and Pap1 (Chen et al., 2008) regulate genes coding for efflux pumps and dehydrogenases that seem necessary for defense against multiple drugs. In fact, pap1 was first identified as a gene whose overexpression conferred resistance to drugs like brefeldin A, staurosporine or caffeine (Toda et al., 1991). Furthermore, Calvo et al. (2012) have shown that in response to H 2 O 2 Pap1 requires the transcription factor Prr1 to activate the antioxidant but not the drug tolerance genes, providing a possible mechanism to explain the major role that NapA shows in the regulation of drug tolerance genes during conidiation. Indeed, the possibility that NapA plays overlapping but different functions during oxidative stress, conidiation and sexual development deserves further research.
Among the NapA-dependent genes, we also identified two members of the cupin superfamily, which includes metaldependent and independent enzymes, as well as catalytically inactive proteins associated with abiotic stress and quiescent structures. AN8614 encodes a single-domain 153 amino-acid cupin showing similarity to germins and oxalate oxidases , some of which surprisingly also show superoxide dismutase activity (Woo et al., 2000). AN8368 encodes a 305 amino-acid protein identified as a spherulin 4like protein. We detected that AN8368 is an ortholog of AN2952, and A. fumigatus and A. clavatus Sph3, which correspond to a glycoside hydrolase essential for the biosynthesis of the exopolysaccharide galactosaminogalactan (Bamford et al., 2015).

The Role of NapA in Carbon Catabolism
Conidia and mycelia from napA mutants are unable to grow on ethanol and arabinose and show reduced growth on fructose, as sole carbon sources. A derepression of alcA during conidiation FIGURE 10 | Glucose starvation induces NapA nuclear localization. NapA does not accumulate in nuclei during growth in arabinose or fructose as sole carbon sources but it does accumulate in nuclei during glucose starvation. Conidia from strain CAM20 (NapA::GFP) were grown for 18 h in MM at 37 • C and then mycelia was shifted to MM with or without glucose or with fructose or glycerol as carbon source for indicated times (0-120 min). Mycelial samples were observed in vivo and photographed every 60 min using Epifluorescence microscopy. Lower panel shows nuclei (DAPI) and NapA::GFP fluorescence in mycelia starved for glucose during 120 min, fixed and photographed. Larger square areas in each picture show enlargements of the areas indicated by smaller squares.
was observed before in conidiation defective mutants carrying alleles of the affected gene fused to the alcA promoter (Arratia-Quijada et al., 2012). The lack of growth on ethanol can be at least partially explained by the fact that NapA is required for full expression of alcR, alcA, and aldA genes during conidiation and presumably also during mycelial growth. Since NapA is needed for normal expression of AlcR and this TF is required for its own transcription, it seems likely that AlcR mediates the regulation of alcP, alcA, alcM, alcS, and aldA exerted by NapA. In this scenario alcR full expression might require CreA derepression, activation by NapA, as well as AlcR autoinduction. The fact that alcR promoter contains putative Ap1 binding sites is consistent with this interpretation. It seems unlikely that CreA represses napA expression because NapA antioxidant function must be required during normal glucose metabolism. In addition, creA d 204 mutants do not show increased resistance to H 2 O 2 and qPCR analysis did not show increased napA mRNA levels in conidia from the creA d 204 mutant, as compared to wild type conidia (not shown). Three lines of evidence indicate that ethanol utilization represents a condition in which NapA antioxidant role becomes more important. First, NapA is required for full expression of most alc genes during conidiation. Second, ethanol growth induces napA expression and NapA nuclear localization. Third, NapA is needed for growth in ethanol. Evidence indicates that ethanol toxicity is mediated by AlcA product acetaldehyde, which itself can cause oxidative stress. Interestingly, it has been considered unlikely that ethanol catabolism alone requires such a high expression levels and subtle regulation of the alc genes (Flipphi et al., 2006). We have added another regulation layer that might help to understand the physiological significance of this complex regulation.
NapA requirement to utilize arabinose and fructose might be explained at different levels. First, NapA is necessary for transcript accumulation of the xylitol/sorbitol dehydrogenase (AN2666) in conidia and possibly in mycelia. AN2666 is 52% identical to A. niger enzyme XdhA, which participates in conversation of xylitol to D-xylulose during arabinose metabolism (de Groot et al., 2007). Xylitol/sorbitol dehydrogenase also shows 77% identity to SdhA, an enzyme involved in reversible transformation between fructose and sorbitol (Koivistoinen et al., 2012). This suggests that AN2666 participates in fructose assimilation and is essential for arabinose assimilation. In addition, arabinose and/or fructose transport might be compromised in napA mutants, considering that NapA regulates several genes encoding proteins with transport activity. Among these, AN10891 and AN2665 encode proteins that are 28 and 26% identical to B. cinerea fructose specific transporter FRT1 (Doehlemann et al., 2005). Second, arabinose and fructose metabolism share some regulatory steps. Within the Aspergilli transcriptional activators AraR, GalR and XilR regulate the metabolic conversion of L-arabinose, D-galactose and D-xylose, respectively. A recent report shows that both AraR and XlnR regulate the pentose catabolism genes, as well as the oxido-reductive D-galactose catabolic pathway (Kowalczyk et al., 2015). However, it is interesting that mutants lacking both AraR and XlnR still show a detectable growth on arabinose, indicating the participation of additional regulatory mechanisms. Also interesting is that A. niger XlnR is somehow related to ROS metabolism, as xlnR null mutants are sensitive to oxidative stress and show a qualitative increase in ROS levels (Raulo et al., 2016). Third, NapA homologs are necessary for expression of enzymes involved in the generation of cellular reducing power (NADPH, GSH and thioredoxins). Although this was not detected during conidiation, during growth A. nidulans napA mutants show a decreased GSH/GSSG ratio and a 1.6 fold decrease in total glutathione content when compared to a WT strain (Thön et al., 2010). This is consistent with the fact that napA mutants are very sensitive to AA/acrolein, known to cause GSH depletion (Kwolek-Mirek et al., 2009). Such a decrease in reducing power might provide additional difficulties in catabolizing carbon sources that are more demanding on these resources. Indeed, ethanol and arabinose catabolism results in an imbalance of redox cofactors (Seiboth and Metz, 2011).
A role for NapA homologs in carbon metabolism appears conserved in fungi. In S. cerevisiae the involvement of Yap1 in responding to carbon stress was proposed after it was recovered in two-hybrid assays as interactor of the Sip2 subunit of the Snf1 kinase, needed for adaptation to glucose limitation, and showing that media shifting from glucose to glycerol or to glucose lacking media induced Yap1 nuclear accumulation (Wiatrowski and Carlson, 2003). Likewise, in S. pombe glucose but not nitrogen starvation induces Pap1 nuclear localization in a process that depends on the stress MAPK Spc1/Sty1 (Madrid et al., 2004). However, a direct involvement of Yap1 or NapA in carbon utilization was not demonstrated in these cases.
In nature, fungal spores germinate on environments containing diverse, often-poor carbon sources as well as drugs produced by competing organisms or that are detritus-derived. In addition, conidia contain chemicals and toxins derived form secondary metabolism as well as auto-inhibitors of germination. Our results support the proposal that oxidative stress is produced during conidiation and that NapA plays a crucial role in the regulation of multiple genes during this process. Although the rates of conidia germination are not decreased in napA mutants under laboratory conditions, the set of genes regulated by NapA during conidiation suggest that they might be important for spores to germinate in natural heterogeneous environments.

AUTHOR CONTRIBUTIONS
JA: Designed experiments, wrote the MS, and obtained funding. AM, FL and OS: Performed and designed experiments, contributed to MS writing.