RNA Sequencing Reveals Differential Gene Expression of Cerrena Unicolor in Response to Variable Lighting Conditions

To elucidate the light-dependent gene expression in Cerrena unicolor FCL139, the transcriptomes of the fungus growing in white, blue, green, and red lighting conditions and darkness were analysed. Among 10,413 all-unigenes detected in C. unicolor, 7762 were found to be expressed in all tested conditions. Transcripts encoding putative fungal photoreceptors in the C. unicolor transcriptome were identified. The number of transcripts uniquely produced by fungus ranged from 20 during its growth in darkness to 112 in the green lighting conditions. We identified numerous genes whose expression differed substantially between the darkness (control) and each of the light variants tested, with the greatest number of differentially expressed genes (DEGs) (454 up- and 457 down-regulated) observed for the white lighting conditions. The DEGs comprised those involved in primary carbohydrate metabolism, amino acid metabolism, autophagy, nucleotide repair systems, signalling pathways, and carotenoid metabolism as defined using Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The analysis of the expression profile of genes coding for lignocellulose-degrading enzymes suggests that the wood-degradation properties of C. unicolor may be independent of the lighting conditions and may result from the overall stimulation of fungal metabolism by daylight.


Introduction
Wood-colonising fungi are known to have unique biochemical pathways enabling them to assimilate a vast array of simple and complex nutrients and to produce a variety of metabolites. Most wood decay fungi are strictly saprotrophic, utilising dead wood as a food base. However, there are commonly known parasitic species of fungi (pathogens) growing and producing sporophores on living trees [1]. Elucidation of the mechanism of biological wood decay is not only important ecologically due to the role in the carbon and nitrogen cycle, but also has economic significance and plays a beneficial role in human life [2]. The environmental impact of saprotrophic fungi arises from their metabolic versatility, which includes the production of a wide range of enzymes directly or indirectly involved in the degradation of organic residues [3]. Individual wood-decomposing fungus produces various enzymes, degrades different plant material, and colonises diverse ecological niches. Therefore, understanding the physiology of wood-degrading fungi and the ability to affect their metabolic potential may have far-reaching impact for not only future research directions but also human life.
In nature, efficient lignocellulose degradation, which occurs during the process of wood decay, is mainly achieved by the activity of white rot basidiomycetes. As known, environmental factors greatly influence decay caused by white rot fungi [4], and light is one of the most important signals for every living cell. Sunlight serves as either a source of energy or information from the environment and can be considered crucial for successful competition, survival, and development in nature [5]. In fungi, light controls developmental decisions, stress response, and physiological adaptations as well as the circadian clock. To sense light, only a few photoreceptor systems have developed during evolution, comprising flavin-based blue-light, retinal-based green-light (such as rhodopsin), and linear tetrapyrrole-based red-light sensors [6]. All these major classes of photoreceptors are found in fungi suggesting that they can detect specific wavelengths by discrete receptor proteins [7].
The emergence of high throughput next-generation sequencing (NGS) techniques has accelerated the discovery of novel genes and has provided an opportunity to further explore and understand the differential gene expression at the genomic and transcriptional level. Since the first published transcriptome study [8], an increasing number of transcriptomes of wood decay basidiomycetes have been reported [9][10][11][12][13][14]. In contrast to the static genome, the transcriptome reflects the differential gene expression in response to developmental or environmental factors, which constitute a dynamic link between an organism's genome and its physical characteristics [15][16][17]. Transcriptome analysis performed recently for Cerrena unicolor, a wood-degrading basidiomycete with great importance in biotechnology and medicine [18][19][20], revealed differential expression of genes coding for proteins engaged in degradation of various kinds of wood [14]. Moreover, the dependence of C. unicolor cultivation in different lighting conditions and synthesis of certain enzymes has also been demonstrated [21].
Here, we report a comparative transcriptome analysis based on high throughput RNA sequencing experiment of C. unicolor FCL139 cultivated in different lighting conditions. Our results provide a comprehensive dataset for analysis of transcriptional profile changes occurring in fungal cells in response to variable lighting conditions, contributing to better understanding of C. unicolor photobiology and light dependent metabolism and behaviour of this environmentally and commercially important fungus.

Transcriptome Sequencing and Identification of Transcripts Coding for Fungal Photoreceptors
An mRNA-Seq approach was employed to gain insight into the gene expression profile of C. unicolor FCL139 during growth on ash sawdust in various lighting conditions comprising white (W), green (G), blue (B) and red (R) light and darkness (D) as a control. Fifteen cDNA samples (three replicates per lighting condition) were prepared from C. unicolor mycelia and sequenced as pair-end reads using the Solid 5500 platform resulting in a total average of 113.6 million (for C. unicolor cDNA samples obtained after growth in green (G) light) to 128.3 million reads (red (R) light) ( Table 1). After filtering off primer-adaptor sequences, low-quality and ambiguous reads, and rRNA mapping reads, the total average number of high-quality reads ranged from 64.5 million for (G) to 93.3 million for the control (darkness, D). These high-quality filtered reads were mapped to the 12,966 predicted genes of C. unicolor using a Subjunc algorithm. Reads from individual sample replicates were concatenated so that there was one mapping per sample (e.g., one for white (W) light, one for green (G) light, etc.). The total number of reads that were uniquely mapped for individual samples ranged from 34.2 million for (G) to over 47.8 million for (W) ( Table 1). Altogether, 10,413 all-unigenes were found for C. unicolor mycelia grown in various lighting conditions, among which 7762 were expressed in all tested culture variants ( Figure 1). Among C. unicolor all-unigenes for which a molecular function was assigned (based on the similarity of the predicted protein sequences to proteins characterised in other organisms or the presence of conserved functional domains within the predicted protein sequence), several genes encoding putative fungal photoreceptors were found ( Table 2). Four of them were annotated as white collar proteins sharing similarity with respective proteins of Neurospora crassa. Two genes encoded proteins similar to the cryptochrome of N. crassa, while putative protein products of two other genes were similar to the phytochrome of Ceratodon purpureus and one to Trametes versicolor opsin ( Table  2).  Among C. unicolor all-unigenes for which a molecular function was assigned (based on the similarity of the predicted protein sequences to proteins characterised in other organisms or the presence of conserved functional domains within the predicted protein sequence), several genes encoding putative fungal photoreceptors were found ( Table 2). Four of them were annotated as white collar proteins sharing similarity with respective proteins of Neurospora crassa. Two genes encoded proteins similar to the cryptochrome of N. crassa, while putative protein products of two other genes were similar to the phytochrome of Ceratodon purpureus and one to Trametes versicolor opsin ( Table 2).

Analysis of Differentially Expressed Genes (DEGs)
For the analysis of DEGs, genes that were expressed with fragments per kilo base per million reads (FPKM) >1 in at least one of the lighting variants were included and considered as "expressed genes". The number of "expressed genes" varied between the tested lighting variants and was 8267 for W, 8539 for G, 8282 for B, 8349 for R, and 8206 for D. The analysis of C. unicolor transcriptomes revealed numerous genes whose expression differed significantly (p-value <0.05) between the dark variant (control) and each of the light variants used (Table 3 and Figure A1). In general, the greatest number of DEGs (454 up-regulated and 457 down-regulated) was observed when the fungus was grown in the white light, compared to the control conditions. In contrast, the lowest number of DEGs (45 up-regulated and 171 down-regulated) was observed during C. unicolor growth in the green light, compared to the control. Only in the white lighting variant, the numbers of up-regulated and down-regulated genes were similar. In the other tested conditions (except for red vs. darkness), down-regulated genes predominated (Table 3). Table 3. Number of C. unicolor genes whose expression differed significantly (at p < 0.05), during fungus growth in the tested lighting conditions in comparison to the control (darkness). Only genes whose transcripts abundance (FPKM) in at least one of the tested conditions was >1 were considered.

Compared Lightning Conditions
Number

KEGG Pathway Enrichment of DEGs
The KEGG pathway analysis was performed for differentially expressed genes ( Figure 2). Generally, most of the DEGs were ascribed into the "metabolic pathways" category in almost all lighting conditions except for the blue light. The genes encoding various metabolic pathways were mostly down-regulated in each of the analysed lighting variants ( Figure 2). The genes representing KEGG pathways related to biosynthesis of secondary metabolites and antibiotics were fairy abundant in all the lighting variants as well. The blue light specifically induced expression of genes potentially engaged in histidine metabolism, while the expression of genes related to retinol metabolism was up-regulated in mycelia grown in the red light. Specific induction of the genes encoding putative proteins belonging to a two-component system, MAPK and RIG-I-like receptor-based signalling pathways, lysine degradation, sphingolipid and pyruvate metabolism, DNA repair pathways, endocytosis and autophagy was observed in the white light conditions. Only a few KEGG categories were found as unique in individual lightning variants when down-regulated genes were considered in the pathway analysis. In the mycelia grown in the blue and green light, a number of down-regulated genes related to ribosome biogenesis were observed. Fungus growth in the red light resulted in repression of putative genes engaged in lysine synthesis, while the green light specifically down-regulated genes involved in the citric cycle. The expression of genes related to the cell cycle, DNA replication, and tyrosine metabolism pathways was found to be down-regulated in C. unicolor grown in the white light, whereas all genes whose expression was influenced by the blue light fell into categories shared by at least one other light variant ( Figure 2).

Light Regulates Primary Metabolism of C. Unicolor
The KEGG pathway analysis revealed that light exerted a global influence on the primary metabolism of C. unicolor. In comparison to darkness (control), during fungus growth in the white lighting conditions, substantially decreased expression of three genes encoding putative enzymes, namely methylenetetrahydrofolate reductase (NADPH) [ genes was observed during fungus growth in the white light culture variant. The genes were potentially engaged in the metabolism of tyrosine leading to its transformation to pyruvate, succinate, or melanin. Most of the genes whose expression was repressed during C. unicolor growth in the white light encoded putative enzymes are involved in purine/pyrimidine metabolism, DNA replication, and homologous recombination or participating in protein processing in the endoplasmic reticulum. On the other hand, the white light induced expression of putative enzymes that degrade lysine, isoleucine, leucine, and valine and enhance degradation of fatty acids. Starch and sucrose metabolism seemed to be influenced in the mycelia growing in the white lighting conditions by the up-regulated expression of genes coding for trehalose 6-

Light Regulates Primary Metabolism of C. Unicolor
The KEGG pathway analysis revealed that light exerted a global influence on the primary metabolism of C. unicolor. In comparison to darkness (control), during fungus growth in the white lighting conditions, substantially decreased expression of three genes encoding putative enzymes, namely methylenetetrahydrofolate reductase (NADPH) [ The genes were potentially engaged in the metabolism of tyrosine leading to its transformation to pyruvate, succinate, or melanin. Most of the genes whose expression was repressed during C. unicolor growth in the white light encoded putative enzymes are involved in purine/pyrimidine metabolism, DNA replication, and homologous recombination or participating in protein processing in the endoplasmic reticulum. On the other hand, the white light induced expression of putative enzymes that degrade lysine, isoleucine, leucine, and valine and enhance degradation of fatty acids. Starch and sucrose metabolism seemed to be influenced in the mycelia growing in the white lighting conditions by the up-regulated expression of genes coding for trehalose 6

Effect of Light on the Expression of Wood-Degrading Enzymes
Since C. unicolor is a white rot fungus and the mycelia were cultivated on ash sawdust during this experiment, analysis of the expression pattern of genes coding for wood degrading enzymes was conducted. Transcripts of genes encoding all groups of enzymes (cellulases, hemicellulases, and lignin-modifying enzymes) were found in the transcriptomes of C. unicolor cultivated in the tested lightning conditions (Table 4). However, the highest number of transcripts of genes related to wood degradation (16, comprising both up-or down-regulated genes) was noted when white light was applied during fungus growth, while the lowest number (only four transcripts) was recorded in the green lighting variant. In all lightning conditions except for the green light variant, the up-regulated genes related to wood degradation predominated. The white light induced mainly the expression of genes encoding lignin-degrading auxiliary (LDA) enzymes comprising alcohol oxidase (XLOC_000665) and five cytochromes P-450 (XLOC_007791, XLOC_005202, XLOC_006061, XLOC_007410, and XLOC_001403). Similarly, the blue light induced the expression of putative LDA enzymes including aryl-alcohol dehydrogenase (XLOC_008717), pyranose 2-oxidase (XLOC_009983), and monooxygenase (XLOC_003435) and the red light influenced the expression of cytochrome P-450 (XLOC_007791, XLOC_00671, XLOC_006683, XLOC_001969, and XLOC_007717). On the other hand, the fungus growth in the red light variant resulted in down-regulation of the laccase coding gene (XLOC_011551), whereas the white light repressed four laccase genes (XLOC_008979, XLOC_008955, XLOC_000669, and XLOC_011551). The expression of genes coding for all hemicellulases and almost all cellulases was induced by the different lightning conditions. Only exoglucanase (XLOC_006587) and endo-β-1,4-glucanase (XLOC_007951) coding genes were down-regulated when the blue and white light variants were applied, respectively.

Signalling Pathways
Besides the above-mentioned genes coding for enzymes that are directly engaged in C. unicolor primary metabolism and wood degradation, the transcripts of several genes coding for proteins of various signalling pathways were detected in the mycelia cultured in the different lightning conditions, indicating that light may be a regulatory factor affecting other cellular processes ( Figure 3 and Table A1). It was observed that the white and blue light induced the expression of the mitogen-activated protein kinase (MAPK) gene, which may be responsible for signal transmission to the pathways related to filamentation, osmoregulation, or cell cycle arrest. Both the white and blue light variants also induced the FoxO signalling pathway (the "Forkhead box" (FOX)). However, the blue light signal seemed to be further transduced via FoxO to the processes of oxidative stress response and DNA repair, whereas the white light may have affected ubiquitin mediated proteolysis. The latter process may also have been affected by the white light via an alternative signalling pathway involving a RIG-I-like receptor, whose expression was induced in the studied conditions. The intracellular calcium level may be additionally regulated by white light via induced expression of cGMP-dependent protein kinase (cGMP-PKG). Several primary metabolism pathways comprising β-oxidation of lipids, glucose uptake, and cholesterol synthesis may be indirectly regulated by white light affecting the AMP-activated protein kinase (AMPK) pathway. The primary metabolic processes mentioned above may also be indirectly regulated by protein phosphatase 2 (PP2A, a known inhibitor of AMPK), whose expression was negatively affected by the red light.

Discussion
Reaction of fungi to illumination is multidirectional, and light can significantly enhance their fitness [5]. Despite the knowledge about the prevalence of photoreceptors in fungi, the influence of light on metabolic processes has not been investigated in detail in most cases. The analyses of fungi response to light is predominantly restricted to white light or daylight effects vs. darkness as the most extreme values [5,7]. In this work, transcriptomic response to the different wavelengths of light was studied in C. unicolor FCL139. From 20 to 112 transcripts were found to be specific to one of the culturing conditions, suggesting that C. unicolor can respond to specific light wavelengths in different manner. This observation is in agreement with the previous studies and the finding that even significant number of transcribed genes of organisms may be engaged in response to light [5,23]. In N. crassa, 31% of the expressed genes were affected two-fold or more by exposure to light [24]. Among all assigned unigenes of C. unicolor, those coding for putative photoreceptor proteins were found. The presence of transcripts coding for putative proteins similar to N. crassa white collar proteins and cryptochromes, C. purpureus phytochromes, and T. versicolor opsins suggests that C. unicolor is able to sense from blue to far-red light. The similarity of C. unicolor receptors to plant phytochromes may be surprising. It should be noticed that no phytochrome-like light sensor protein

Discussion
Reaction of fungi to illumination is multidirectional, and light can significantly enhance their fitness [5]. Despite the knowledge about the prevalence of photoreceptors in fungi, the influence of light on metabolic processes has not been investigated in detail in most cases. The analyses of fungi response to light is predominantly restricted to white light or daylight effects vs. darkness as the most extreme values [5,7]. In this work, transcriptomic response to the different wavelengths of light was studied in C. unicolor FCL139. From 20 to 112 transcripts were found to be specific to one of the culturing conditions, suggesting that C. unicolor can respond to specific light wavelengths in different manner. This observation is in agreement with the previous studies and the finding that even significant number of transcribed genes of organisms may be engaged in response to light [5,23]. In N. crassa, 31% of the expressed genes were affected two-fold or more by exposure to light [24]. Among all assigned unigenes of C. unicolor, those coding for putative photoreceptor proteins were found. The presence of transcripts coding for putative proteins similar to N. crassa white collar proteins and cryptochromes, C. purpureus phytochromes, and T. versicolor opsins suggests that C. unicolor is able to sense from blue to far-red light. The similarity of C. unicolor receptors to plant phytochromes may be surprising. It should be noticed that no phytochrome-like light sensor protein has been characterised in basidiomycetes, and plants and fungi share the same structure of the conserved photosensing region of the receptor [25], which may suggest an unusual structure of the carboxy-terminal domain in C. unicolor.
The analysis of DEGs performed for C. unicolor revealed numerous genes whose expression varies considerably during fungus growth. The KEGG pathway analysis of differentially expressed genes revealed significant enrichment of metabolic pathways indicating that light influenced mainly primary metabolism of C. unicolor. The expression of genes engaged in metabolic pathway was mostly down-regulated in each of the analysed lighting variant except for blue light. Concomitantly, in white light conditions, up-regulation of genes engaged in amino acid degradation, autophagy, nucleotide repair systems, signalling pathways, and carotenoid metabolism was observed. In turn, 319 differentially expressed carbohydrate-active enzymes participated in carbon metabolism were identified in Pleurotus eryngii under blue light stimulation [26]. In N. crassa, light increased cellular metabolism causing simultaneously significant oxidative stress to the organism. In response to this stress, protective photopigments and antioxidants were produced, and genes involved in ribosome biogenesis were transiently repressed [24]. The number of down-regulated genes related to ribosome biogenesis was also observed in the C. unicolor mycelia grown in the blue and green light. In Aspergillus ornatus, glucose metabolism was decreased by light prior to production of conidia [5]. While blue light is a commonly known signal triggering up-regulation of carotenogenesis [27], the red light induction of genes related to respective pathway was observed in C. unicolor. Transcriptome analysis of Cordyceps militaris albino mutant and its sibling strain revealed that the previously identified photoreceptors were expressed similarly in response to light. Both strains shared 235 DEGs after the light treatment and all DEGs were grouped into three functional categories: biological processes, cellular components, and molecular functions [28]. Similar discrepancies in the metabolism in different fungal species were previously reported in two model organisms: N. crassa and Aspergillus nidulans [29]. Light triggers significant expressional changes within the first 24-48 h in A. nidulans and darkness results in a massive transcriptional reprogramming of gene expression [30].
In Trichoderma reesei cultivated in darkness, significant increase in transcripts coding for cellulases was observed [31]. This phenomenon may be regulated by cAMP-dependent protein kinase (PKA), which influences cellulase synthesis positively in light and negatively in dark conditions. In C. unicolor positive changes in the number of transcripts involved in the AMPK pathway was observed in white lighting. However, only one transcript coding for cellulase was found to be up-regulated in these conditions. Tang et al. [32] proved that light-induced brown film formation in Lentinula edodes occurred via MAPK and a two-component system, with possible engagement of white collar proteins or phytochromes. Our results show that both pathways (AMPK and MAPK based) were induced by the white light in C. unicolor. Moreover, up-regulation of genes coding for cytochrome P-450 and monooxygenase by the red and blue light, respectively, which was observed in our analysis, seem to support the role of putative white collar proteins or phytochromes in AMPK and MAPK based signal transduction. Phytochromes as red light receptors were also reported to be involved in regulation of secondary metabolism in A. nidulans [33] and our data suggest that synthesis of several enzymes involved in secondary metabolism was induced by the light of same wavelength. Furthermore, in C. unicolor, the MAPK pathway induction by the white light resulted in enhanced filamentation, which was proved to be a stress response in bacteria to DNA damage [34]. MAPK was suggested to play an important role as part of the DNA repair system in Trichoderma atroviride [35]. In turn, PP2A was proved to be part of the FREQUENCY(FRQ)-based oscillator with circadian properties in N. crassa [36]. PP2A expression was down-regulated by red light in our experiment, raising the possibility that C. unicolor is sensitive to day-night changes.
Recently reported comparative transcriptomic analysis of C. unicolor performed in uncontrolled lighting conditions revealed differential expression of genes engaged in wood degradation during fungus growth on ash, maple, and birch sawdust, and significant differences were observed for cellulase-coding genes [14]. Currently, transcripts of genes coding for cellulases, hemicellulases, and lignin-modifying enzymes were found during fungus growth on ash sawdust medium in all lighting conditions, which may suggest that wood-degradation properties of C. unicolor depend on the medium applied rather than lighting conditions [14,37]. The high number of transcripts related to wood degradation observed during fungus growth in the white light may have resulted from the overall stimulation of its metabolism by daylight [38]. It has been demonstrated recently that blue light efficiently boosted laccase synthesis in C. unicolor and Pycnoporus sanguineus, whereas in Phlebia lindtneri the highest laccase activity were observed in green lighting conditions in LH or cellulose-based medium [21]. Laccase was also up-regulated 10 h after blue light exposure in Coprinopsis cinerea [39]. On the contrary, green light efficiently inhibited laccase in other P. sanguineus strains and maximum enzyme activity was noticed in the darkness [40].
In conclusion, understanding how C. unicolor and other fungi sense the environment is of intellectual, industrial, and medicinal importance. The detection of transcripts encoding for putative fungal photoreceptors in C. unicolor transcriptome strongly suggests that the fungus can sense light and respond to various wavelength of light in a slightly diverse manner, by changing its overall metabolism. The up-regulation of gene expression of several signalling pathways, which can be considered as a hub integrating and transducing various environmental signals, implies that light together with other factors may affect cellular processes of the fungus.

Medium and Growth Conditions
The C. unicolor FCL139 strain was obtained from the culture collection of the Regensberg University. The fungal culture was maintained in 2% (w/v) malt agar slants. As an inoculum, pieces of agar culture were grown in the Lindeberg and Holm (LH) medium [41] in non-agitated conical flasks for 10 days at 28 • C. Next, the mycelial mats were collected and homogenised in a disperser (IKA, Warszawa, Poland). The fragmented mycelial culture (1% v/v) was used as a standard inoculum for further studies. Solid-state lignocellulose C. unicolor cultures were grown at 28 • C on 1 g of sterile ash sawdust (wood particles < 4 mm) soaked with 9 mL of distilled water and solidified with 1.5 % of agar in Petri dishes covered with a nylon membrane filter (0.45 µm). The seeds were then cultivated at 28 • C in incubators equipped with illumination led cassettes (KT 115, Binder, Neckarsulm, Germany) for 8 days. Continuous lighting conditions (20 lux) were provided throughout the entire period of C. unicolor cultivation. The following lighting variants were applied: white (W) (4000-4750 K), green (G) (510-520 nm), blue (B) (465-470 nm), and red (R) (620-625 nm) light and darkness (D) as a control.

RNA Extraction and Sequencing
Fifty milligrams of C. unicolor mycelia were frozen in liquid nitrogen and hand ground to a fine powder with a pre-chilled mortar and pestle. Polyadenylated (polyA) mRNA was directly isolated using oligo (dT) magnetic beads. The Dynabeads mRNA DIRECT Kit (Ambion, Waltham, MA, USA) designed for simple and rapid isolation of polyA mRNA directly from the crude lysate was employed for polyA mRNA extraction following the manufacturer's protocol. After homogenisation of hyphae, the DNA shearing step (allowing reduction of viscosity) was performed by forcing the lysate 3-5 times through a 21-gauge needle using a 2 mL syringe. The mRNA quantity and quality were evaluated using a Qubit 2.0 fluorimeter and a Qubit RNA HS Assay Kit (Thermo Fisher Scientific, Waltham, MA, USA), as well as an Agilent 2100 Bioanalyzer and an RNA 6000 Nano Kit (Agilent Technologies, Santa Clara, CA, USA). The transcriptome library for each sample was prepared using a Solid Total RNA-Seq Kit (Applied Biosystems, Waltham, MA, USA) following the manufacturer protocol for low-input RNA samples. The median size of the mRNA fragments after RNase III fragmentation ranged from 138 to 195 bp. The median size of the amplified cDNA whole transcriptome libraries ranged from 225 to 356 bp. Each of the 15 constructed cDNA libraries (3 replicates per lighting variant) was sequenced using the Solid 5500 platform. The sequencing data were derived from the paired-end reads (50 bp forward and 35 bp reverse) of the template in the Solid templated bead, using forward and reverse ligation chemistry (F3, FWD1 and F5, RNA REV1 Seq. Primers, respectively) in a single run. The processing of original images to sequences and base-calling for the obtained paired-end reads were performed by the LifeScope pipeline (version 1.6) (Life Technologies, Carlsbad, CA, USA).

NGS Data Analysis
The SOLiD sequencing results of paired-end reads in triplicates for each C. unicolor sample were converted into FASTQ format using the solid2fastq program from the BFAST package [42] with default settings. Adapter trimming and removal of low-quality reads was performed with cutadapt [43]. Filtered and processed results of individual cDNA replicates reads of C. unicolor are summarised in Table 1. The C. unicolor rRNAs were identified de novo with Barrnap 0.7 (Available online: https://github.com/Victorian-Bioinformatics-Consortium/barrnap). Based on these predictions, rRNA sequences were fetched from the C. unicolor genome using a custom Python script, and the rRNA-mapping reads were excluded from the sequencing results using Bowtie 2 [44]. The C. unicolor v1.1 (Cerun2) genome assembly, downloaded from the U.S. Department of Energy Joint Genome Institute (DOE JGI, http://jgi.doe.gov), was used as the reference for the mapping of the obtained filtered reads. Reads from each of the sample replicates were concatenated using a Subjunc algorithm of the Subread software package so that there was one mapping per sample to the reference genome [45]. Re-annotation of the C. unicolor genome based on available JGI annotation and high-throughput sequence data of five transcriptomes was performed with the Cufflinks software package [46]. The gene and transcript identifiers assigned by Cufflinks were used in further analyses. The default BLAST E-value cutoff 10 −5 was used during similarity searches.
Transcript abundance was calculated on the basis of the number of mRNA-Seq reads mapped to a given gene model with RSEM, accurate transcript quantification software, from RNA-Seq data [47] and normalised to TPM (transcripts per million) and FPKM (fragments per kilobase of exon per million fragments mapped) values, respectively. Identification of the differentially expressed genes (DEGs) based on the abundance of transcripts and isoforms of the same gene between mycelia grown in darkness (control) and each of the lighting variants was performed with the DESeq 2 package [48].
Individual genes and transcripts were automatically annotated on the basis of predicted sequence similarities to known proteins of the SwissProt and TrEMBL databases as well as the presence of conserved domains (PFAM database).

Nucleotide Sequence Accession Numbers
The raw data from RNAseq experiments have been deposited to the Sequence Read Archive database of NCBI under BioProject number PRJNA401095.

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.