Tricking Arthrinium malaysianum into Producing Industrially Important Enzymes Under 2-Deoxy D-Glucose Treatment.

This study catalogs production of industrially important enzymes and changes in transcript expression caused by 2-deoxy D-glucose (2-DG) treatment in Arthrinium malaysianum cultures. Carbon Catabolite Repression (CCR) induced by 2-DG in this species is cAMP independent unlike many other organisms. Higher levels of secreted endoglucanase (EG), β-glucosidase (BGL), β-xylosidase (BXL), and filter paper activity assay (FPase) enzymes under 2-DG treatment can be exploited for commercial purposes. An integrated RNA sequencing and quantitative proteomic analysis was performed to investigate the cellular response to 2-DG in A. malaysianum. Analysis of RNASeq data under 2-DG treated and control condition reveals that 56% of the unigenes do not have any known similarity to proteins in non-redundant database. Gene Ontology IDs were assigned to 36% of the transcripts (13260) and about 5207 (14%) were mapped to Kyoto Encyclopedia of Genes and Genomes pathway (KEGG). About 1711 genes encoding 2691 transcripts were differentially expressed in treated vs. control samples. Out of the 2691 differentially expressed transcripts, only 582 have any known function. The most up regulated genes belonged to Pentose Phosphate Pathways and carbohydrate degradation class as expected. In addition, genes involved in protein folding, binding, catalytic activity, DNA repair, and secondary metabolites were up-regulated under 2-DG treatment. Whereas genes encoding glycosylation pathways, growth, nutrient reservoir activity was repressed. Gene ontology analysis of the differentially expressed genes indicates metabolic process (35%) is the pre-dominant class followed by carbohydrate degradation (11%), protein folding, and trafficking (6.2%) and transport (5.3%) classes. Unlike other organisms, conventional unfolded protein response (UPR) was not activated in either control or treated conditions. Major enzymes secreted by A. malaysianum are those degrading plant polysaccharides, the most dominant ones being β-glucosidase, as demonstrated by the 2D gel analysis. A set of 7 differentially expressed mRNAs were validated by qPCR. Transmission electron microscopy analyses demonstrated that the 2-DG treated cell walls of hyphae showed significant differences in the cell-wall thickness. Overall 2-DG treatment in A. malaysianum induced secretion of large amount of commercially viable enzymes compared to other known species.


INTRODUCTION
Endophytes represent a promising group of organisms, as they are mostly untapped reservoirs of metabolic diversity. They are often able to degrade cellulose and produce an extraordinary diversity of metabolites (Gianoulis et al., 2012). Arthrinium sp. are endophytes grown in diverse ecological niches, also known to produce bioactive compounds that have pharmacological and medicinal applications (Crous and Groenewald, 2013). The first Arthrinium sp., A. saccharicola was isolated from Saccharum officinarum and was later reported from various habitats such as plants, air and seawater (Miao et al., 2006). A. malaysianum is a known contaminant on Termitomyces cultures (Sawhasan et al., 2012) and in the present study, this organism is isolated from Termitomyces mycellial mats. Reports of several other endophytes belonging to Arthrinium sp. (A. arundinis, A. saccharicola) isolated from marine brown algae Sargassum sp. or from a mixed culture of Termitomyces isolated from termite comb, displayed considerable cellulolytic activities along with antifungal and antioxidant activity (Hong et al., 2015). Arthrinium saccharicola, Arthrinium affsacchari, Arthrinium phaeospermum isolated from Miscanthus sp. are known to produce industrially important enzymes (Shrestha et al., 2015).However, protein secretory pathways and the spectrum of proteins natively secreted by this fungus to the extracellular medium still remain virtually unexplored. The production and storage of several secondary metabolites from Arthrinium sp. have been reported earlier. Bostrycin, ergosterol, succinic acid, various phenolic compounds and antibiotics such as 1-threo-β-hydroxiaspartic acid are some of the known pharmaceutically important products from these genera. Papulacandins, with inhibitory and antifungal activity against bacteria, molds, yeasts, and fungi are also some of the most commonly found secreted products from this group (Traxler et al., 1980;Larrondo et al., 1996;Aissaoui et al., 1998;Adelantado et al., 2005;Calvo et al., 2005).
Carbohydrates are the primary and preferred source of metabolic carbon for most organisms and are used for generating energy and producing biomolecules (de Sousa Lima et al., 2014). The gene Mig1/CreA/CRE1, a zinc-finger transcription factor conserved in most fungal species mediates carbon catabolite repression (CCR) in presence of easily metabolized carbohydrates such as glucose (Ruijter and Visser, 1997). Relief of CCR by glucose depletion or deletion of creA was found to activate expression of genes involved in cell wall-degrading enzymes (CWDEs). Besides being required for CWDE gene expression, alleviation of CCR is also needed for the activation of alternative metabolic pathways like ethanol assimilation and virulence. These proteins that are responsive to carbon catabolite (de) repression, are generally soluble. CAZymes are released under CC (de) Repression. Low molecular weight carbohydrates can act as inducers for the expression of other genes encoding CAZymes, either working individually or together (Delmas et al., 2012). In Neurospora crassa, a number of genes that encode cellulolytic and pectinolytic enzymes show increased transcription levels upon carbon starvation as well as during early exposure to cellulose or pectin (Coradetti et al., 2012;Benz et al., 2014).
In recent years, high throughput next generation sequencing technologies like Roche, Illumina, and Applied Biosystems SOLiD, Ion Torrent are increasingly being used for non-model organisms or where reference genome is not available. Several investigations have been successfully carried out using these techniques lately (Zhou et al., 2014). Recently several de novo transcriptome sequencing using Illumina platform has been reported in species such as Safflower (Lulin et al., 2012), Anopheles (Crawford et al., 2010); Eucalyptus (Mizrachi et al., 2010); making it an attractive option for researchers.
In this study, 2-deoxy D-glucose was used in the growth medium as metabolite for down regulating carbon catabolic repression in A. malaysianum followed by induction of CWDE. Glucose is a major structural precursor in filamentous fungi as the dominant monosaccharide component for cell wall βglucans. When glycolysis is inhibited by 2-DG in aerobic condition, cells continue to produce ATP using alternative energy sources, i.e., amino and fatty acids metabolism. Pentose Phosohate Pathway plays an important role in the production of proteins and therefore this pathway is expected to be relevant especially in an organism with an efficient protein production system (Wick et al., 1957;Chen and Gueron, 1992;Limón et al., 2011). This is a first of its kind study in endophytic fungus A. malaysianum cataloging the metabolic and transcriptional changes taking place under 2-DG regulation. This study further substantiates if 2-DG can be used as a precursor for production of commercially applicable enzymes.

Growing Cultures in Control and Treated Conditions
A. malaysianum was isolated as a contaminant from the mat of Termitomyces clypeatus MTCC 5091 and cultured in potato dextrose agar. The mycelia of A. malaysianum was later transferred to 100 ml of growth media in 500 ml shake-flasks containing (%, w/v), Glucose-5, Malt extract-5, Potato extract-20%. Apart from this, other components and micronutrients such as ammonium dihydrogen phosphate 2.5, (%, w/v) of CaCl2.2H2O, 0.037; KH2PO4, 0.087; MgSO4.7H2O, 0.05; Boric acid,0.057; FeSO4.7H2O, 0.025; MnCl2.4H2O, 0.0036; NaMoO4.4H2O, 0.0032; ZnSO4.7H2O 0.03 was added to the culture at pH 5.0, and incubated at 30 • C under shaking condition (250 rpm) in orbital shaker. The mycelia mats were harvested by filtration through a double-layered muslin cloth and washed several times with sterile distilled water. The mycelia mats were then stored in −80 • C (with excess liquid squeezed out), for 18s RNA sequence extraction. For determination of dry weight, samples collected at various time intervals were centrifuged at 9000 × g for 30 min, and the resulting supernatant was filtered through a 0.45 µ (Millipore).The dry weight of mycelia were measured after repeated washing of the mycelial pellets with distilled water and drying overnight at 70 • C. The biomass concentration was measured as grams of dry weight.
For studying 2-DG response on the cultures, mycelium of A. malaysianum was added to the media containing a single carbon source (1% of glucose or cellobiose). pH of the medium was adjusted to 5.0. The cultures were grown under shaking conditions at 30 ± 2 • C for requisite number of days. A separate set containing the above media ingredients along with addition of 2-DG (0.05%) was also done simultaneously.

Free-Radical Scavenging Activity Assay
Each culture was incubated for 7 days and filtered to separate the mycelium from the media. 100 mg wet mycelial mats were extracted and thoroughly washed with distilled water and homogenized on ice, for 5 mins. The mycelial suspension was centrifuged at 15,000 rpm for 20 min at 4 • C. The pellet was resuspended in distilled water and sonicated on ice 4 times for 1 min using ultrasonic sonicator at full amplitude. The suspension was centrifuged at 8000 rpm for 10 min at 4 • C to precipitate the coarse particle. The precipitate was then washed, air dried and stored at 4 • C. The antioxidant activity of aqueous fungal extracts and the standards were assessed on the basis of the radical scavenging effect of the stable DPPH (2, 2-Diphenyl-1-pycrylhydrazyl; Leung et al., 2006). The DPPH, a stable free radical with the characteristic absorption at 517 nm, was used to study the radical scavenging effects on the fungus. The image of hyphae was captured by using light microscopy.
β-Glucosidase (E.C.3.2.1.21) assay was carried out in the reaction mixture containing 2 mM pNPG as substrate in 0.1 M sodium acetate buffer at pH 5.0 and an appropriate amount of enzyme. One unit of β-Glucosidase activity was defined as the amount of enzyme needed to liberate 1 µM of p-nitrophenol (pNP) per minute under the standard assay conditions. β-Xylosidase (E.C.3.2.1.37) was assayed using a similar method with p-nitrophenyl-β-D-xylopyranoside (pNPX) as substrate. Xylanase was assayed using birch wood xylan as substrate. The solution of 0.9 ml xylan (1%, w/v) and the 0.1 ml of enzyme at appropriate dilution were incubated at 50 • C for 10 mins. The reducing sugar was determined by the DNS method at 540 nm with xylose as standard. Filter paper degrading activity (FPA) was assayed according to the IUPAC method using 1 × 6 cm (50 mg) Whatman No.1 filter paper strips as a substrate.

Total RNA Isolation from Control and Experimental Fungus Samples
Total RNA was isolated from 100 mg of tissue using Xcelgen Plant RNA Miniprep kit. The quality of RNA was analyzed on 1% denaturing agarose gel. The paired-end cDNA sequencing libraries were prepared for both samples, separately using Illumina TruSeq RNA Library Preparation Kit. Library preparation was started with mRNA fragmentation followed by reverse transcription, second-strand synthesis, pair-end adapter ligation, and finally ended by index PCR amplification of adaptor-ligated library. Library quantification and qualification was performed on Caliper LabChip GX using HT DNA High Sensitivity Assay Kit.
The expression levels of assembled unigenes/transcripts were measured using fragments per kilobases per million reads (FPKM) method. We mapped treated and control assemblies back into merged assembly for annotation consistency using BLAT (Kent, 2002).

Generating CDS From Unigenes Using Codon Preference
The identification of coding sequences (CDS) from expressed transcripts still remains a challenge. Although many methods exist, a modified log likelihood approach based on codon preference (Staden and McLachlan, 1982) works best. Initially, we generated a codon usage table using CUSP from EMBOSS (Rice et al., 2000), where the entire transcript set was used for generating training dataset. Then we iterated this process several times using the modified outputs from the previous dataset and extracted codon preference data from that and used them for training in the subsequent prediction. Fourth iteration produced the best CDS that were subsequently used for analysis.

qPCR Protocols
Based on their functional importance, seven genes were selected for validating the results of the RNA-Seq analysis. The primers for these genes were designed with Primer Express 3.0 software (Applied Biosystems). Total RNA was extracted from both sets of independent cultures grown in absence (control) and presence (treated) of 2-DG in cellobiose media as described above, and then converted to cDNA by random priming. PCR reactions were run in four replicates. The transcription level of genes was determined according to the 2 − CT method, using actin as a reference gene for the normalization of gene expression levels.

Study of A. malaysianum Ultra-Structure by Transmission Electron Microscopy (TEM)
Ultra structural damage of 2-DG treated mycelia was ascertained through TEM studies carried out with a TECNAI SPIRIT (FEI, The Netherlands) instrument having an accelerating voltage of 60 kV. Culture aliquots were centrifuged at 10,000 rpm for 10 min; the mycelial pellet was fixed with 2.5% glutaraldehyde in 0.25 M phosphate buffer (pH 7.2) for 24 h, washed thrice with the same buffer and post fixed with 1% w/v osmium tetroxide (2 h). This was followed by 2 h incubation with 0.5% w/v aqueous uranyl acetate solution followed by dehydration using a graded alcohol series (from 50 to 100% v/v) and embedding in spurr. Mycelial thin sections (60 nm) were cut with a diamond knife (DUPONT, Germany) in a UC6 Leica Ultra microtome (Germany) and mounted on copper grids before examination by TEM.
Water extracts of mycelia of both treated (cellobiose +2-DG) and control (cellobiose) of A. malaysianum were subjected to possible antioxidant activity by the DPPH free radical scavenging method. The model of scavenging the stable DPPH radical is widely used to evaluate antioxidant activities over a relatively short time compared to other methods. Scavenging effect of mycelia extracts on DPPH radicals was found to be concentration dependent. Free radical scavenging effect of 8.1% at a lower concentration of 2 mg/ml of mycelia extracts was achieved in 2-DG treated mycelium as shown in Table 1B.

Proteome Analysis of β-Glucosidase Secretion
Total proteins secreted to extracellular milieu during cultivation of A. malaysianum in treated (cellobiose in presence of 2-DG) and control (cellobiose in absence of 2-DG) medium was extracted at different time interval (48, 60, 72 h) and was analyzed by 2-DGE analysis, ( Figure 2B). Most of the protein spots were concentrated in pI range 4-7 and between molecular mass of 30 and 120 kDa. It was clear from the gel images that extracellular β-glucosidase in absence of 2-DG was produced after 60 h of cultivation whereas extracellular β-glucosidase in presence of 2-DG was produced within 48 h and secreted into the medium ( Figure 2B).

Transcriptome Assembly and Annotation
We generated deep transcriptome sequences with an approximate coverage of 250X from control and 2-DG treated samples using Illumina platform. Transcript reads from control and treated samples were assembled separately and together (merged assembly) for analyzing differential expression. Out of the several assemblers used, Trinity assembler (Grabherr et al., 2011) generated the best results. Further, downstream data analysis was carried out using merged assembly. We used dinucleotide relative abundance (Gentles and Karlin, 2001) to calculate possible contaminants and removed about 84 sequences from the merged assembly. Final assembly resulted in 21192 genes (termed as unigenes) coding for 27252 transcripts from control samples, 22372 unigenes coding for 32982 transcripts from 2-DG treated samples and 22847 unigenes encoding 36248 transcripts from merged samples ( Table 2). The length distribution of the assembled transcripts are represented in Figure 4A.
Putative functions were assigned to 44% (16,096) transcripts based on their similarities to non redundant database (Supplementary Table 1). Out of which, 152 (0.4%) transcripts matched with ribosomal genes. A significantly large number of orthologs were detected with Ascomycetes fungus Eutypa followed by Bacteria ( Figure 4B).
Taxonomy analysis A. malaysianum unigenes indicated that it had 55% sequence similarity with Eutypa lata UCREL1 (3588) species, followed by Colletotrichum higginsianum (1499)   Carbon content(g/l) 10 ± 0.001 10 ± 0.001 10 ± 0.001 large part of the secretome was classified as hypothetical proteins followed by phosphorylation related proteins, Glycosyl hydrolases and transferase. The major gene ontology processes in A. malaysianum were consistent with those of Eutypa lata, Gaeumannomyces graminis, Thielavia sp., Togninia, and Colletotrichum (O'Connell et al., 2012;Jaramillo et al., 2015). Among the GO classifications, 7351 transcripts were assigned to molecular process; 230 under cellular component and 838 under biological process. Under each of the three main categories (biological process, cellular component and molecular function), "metabolic process, ", "cell, " "cellular process, " "binding, " and "catalytic activity" terms occupied the largest proportion in both the control and experimental samples. We also noticed a high percentage of experimental transcripts from categories of "transporter activity, " "organelle" and "protein binding." "Nutrient reserve activity, " "viral reproduction" and "protease activity" was more prominent in 2-DG treated transcripts.

ABOUT 7.3% GENES ARE DIFFERENTIALLY EXPRESSED UNDER 2-DG TREATED CONDITIONS
We used a fold change value of >=2 for estimating differentially expressed genes. A total of 2691 transcripts encoded by 1711 genes were differentially expressed under control and treated conditions. Out of this, 1150 transcripts were down regulated and 1540 transcripts were up-regulated. The list of genes and their individual functions with their KO id's are provided in Supplementary Table 3. Out of the differentially expressed genes, the number of carbohydrate metabolism (11%) genes was the highest in number followed by secondary metabolites (5.1%), transporters (5.3%), binding proteins (4.8%), catalytic activity (4.6%), and cell rescue and virulence genes (3.7%; Figure 5). Differentially expressed transcripts were analyzed further on several key functions such as (i) carbohydrateactive enzymes(CAZymes) and cell wall biogenesis (ii) Secondary metabolism, unfolded protein response and virulence genes iii) FIGURE 2 | Growth, carbon consumption, and synthesis of β-glucosidase in A. malaysianum in presence of 0.5 mg/ml of 2DG. (A) Time dependent profile of fungal growth, cellobiose utilization, and kinetics of intracellular and extracellular β-glucosidase production, (B) 2D Gel electrophoresis analysis of the secreted proteins of A. malaysianum at three different hours. Extracellular β-glucosidase secreted into the culture medium in absence (A,B) in presence of 2DG (0.5 mg/ml) at 48, 60, and 72 h is highlighted as (Muller et al., 2015). metabolic pathways and enzymes, and (iv) Other functions mainly on cellular transporters and cell repair.

Enhanced Expression of Cazymes and Cell Wall Homeostasis During 2-DG Treatment
Enzymes that break down plant cell wall can be categorized into the following classes: glycoside hydrolases (GH), glycosyl transferases (GT) polysaccharide lyases (PL), carbohydrate esterases (CE), and auxiliary activities (AA). All CAZymes and auxiliary activity enzymes can be attached to carbohydratebinding modules (CBMs). Our transcriptomic analysis revealed a vast majority of carbohydrate-active enzymes encoded by A. malaysianum. About 3447 CAZymes were detected from the merged transcriptome (Supplementary Table 4).
Majority of CAZymes belonged to GH (117)     can be regarded as a hallmark of convergent evolution of A. malaysianum caused due to mixed culture with Termitomyces clypeatus.

2-DG Triggers Expression of Virulence Associated Genes, Unfolded Protein Response And Biosynthesis of Secondary Metabolites
Virulence factors are difficult to define by classical terminology. Intrinsically associated with pathogenicity, such as endotoxin and molecules referred to as modulins (e.g., cell wall polysaccharides and lipids), can contribute to the pathogenic process. Transcriptome analysis A. malaysianum in 2-DG treated condition revealed the presence of potential virulence associated genes (Table 4). Some requisite virulence factors, such as the toxins or polysaccharide capsules proteins, betalactamase like protein were expressed suggesting the pathogenic nature of the organism in stressed environmental condition. Constitutively expressed components of microorganisms, such as endotoxin and molecules referred to as modulins (e.g., cell wall polysaccharides and lipids), that contribute to the pathogenic process were also expressed. Surprisingly, several genes with important functions in protein synthesis and trafficking were also upregulated ( Table 3). Many of them were proteins involved in folding and maturation like pdi-1, ER resident chaperones, bip-1, grp 78, and several heat shock proteins. These genes were induced in response to increased up regulation of CWDE influenced by 2-DG. These data suggested that the secretory pathway in this fungus regulates increased lignocellulase synthesis and secretion by up regulating the ER stress responsive proteins to reduce ER load (Figures 6, 7). Protein disulphide isomerase (PDI) is one of the major ER foldases. The enzyme catalyses the formation and rearrangement of disulphide bridges in proteins during folding. The pdi gene was upregulated in treated condition.and its expression level was qPCR ( Figure 8B). The differential expression levels of several chaperons and binding protein were also upregulated (Figure 6), although the expression of bZIP transcription factor and HAC 1 showed minor increase (Figures 6, 8B).
Secondary metabolites are involved in intracellular, intercellular and interspecific interactions. A. is known to produces a wide variety of secondary metabolites, and this motivated us to find the molecular basis of this production by de novo transcriptomics. We identified 30 core genes related to secondary metabolism including 12 polyketide synthase (PKSs), four putative PKS-like enzymes, seven fattyacid synthases, (Monooxygenase FAD-binding protein etc. Supplementary Table 1; Octave et al., 2006).
Other artifacts of 2-DG treatments in the culture media

2-DG Reduced Cell Viability Followed by Cell Death
A large number cell repair, senescence, p38, mapk, glutathione-s-transferase, pho81 was induced by 2-DG  Table 4). This is in concurrence to the condition induced by carbon starvation in fungi (Emri et al., 2004;Pócsi et al., 2004;Katz et al., 2006;Szilágyi et al., 2013). This data suggest that macrophagy is a part of late response to 2-DG in A. malaysianum as observed in old cultures by gradual pigmentation and mycelia fragmentation. However, a more detailed study is required to further establish this effect. The production of cell wall hydrolyzing enzymes may also trigger cell death process in old carbon starved cultures . In our study, there was up regulation of cell wall integrity pathway genes such as WSC, SNF-1, and other map kinases that are usually present in carbon starvation or depletion like condition ( Table 3). The degradation of glutathione is significant, because glutathione is utilized as an energy source, supports intracellular amino acid transport and protects cells against ROS. Interestingly, 2DG resulted in up regulation of the expression of glutathione-s-transferase but down regulation of superoxide (Figure 6).

Expanded Transporter Gene Families in Treated Condition
In plant pathogens, cellular transporters are responsible not only for export of compounds involved in pathogenesis and virulence, but they also may play an essential role in protection against plant defense compounds (e.g., secondary metabolites) during pathogenesis, possibly by exporting hostderived antimicrobial compounds out of the cell (Pao et al., 1998;Schoonbeek et al., 2002). The transportation system is involved       in uptake of essential nutrients and ions, excretion of metabolic end products, and deleterious substances, and communication between cells and the environment. Transporters comprised of 5.3% of our classified transcripts (the second largest cluster after carbohydrate biosynthesis). Specific transporters elevated such as those of ABC, hexoses, solute, sugar:H+ symporter, MFS transporters. A total of 46 genes encoding transporters were identified in the merged transcriptome highest order of MFS transporters. MFS transporters are involved in the transport of monosaccharides, oligosaccharides, inositols, drugs, amino acids, nucleosides, organophosphate esters, Krebs cycle metabolites, and a large variety of organic and inorganic anions and cations.Compared to the control, a significant increase in MFS transporters was observed in the transcriptome of treated. The gene number of sugar transporter (SP) family of MFS subfamily was higher indicating the uptake of more carbohydrate-produced nutrients. Comparative analysis with other fungi revealed that the MFS alpha-glucoside transporter, MFS multidrug transporter is overrepresented, suggesting export of more metabolism products in the extracellular. Presence of a high affinity nicotinate permease that catalyzes nicotinic acid (vitamin B3) uptake, vitamin H transporter, putative vitamin b6 transporter, bsu1 protein, vitamin H transporter 1/biotin symporter vht1, reflecting that might be dependent from the host plant for vitamin supply (Table 3).

Validating RNAseq Results with qPCR Techniques
Seven differentially expressed transcripts, of which most of them were related to metabolism (Hexokinase, Fructose bisphosphate, Isocitrate lyase, Transaldolase), protein folding (Protein disulfide isomerase, HAC) and cell wall degrading enzyme (β-glucosidase) were selected for validating the results obtained RNAseq analyses using qPCR (Table 5). In this analysis, a comparison has been made for the expression of the chosen genes between the Control (absence of 2-DG) and treated (presence of 2-DG) as illustrated in Figure 8. This study showed that most of the gene related to metabolism and protein folding were found to be upregulated in 2-DG treatment compared to the control, which showed no such up-regulation. We examined transcript levels of β-glucosidase differentially expressed gene through qPCR, which were also found to be differentially expressed at 2-D proteomic analysis ( Figure 2B). The expression levels of two well-known UPR maker genes such as PDI and HACA was determined using qPCR and was found that though there was increased expression of Protein disulfide isomerase the expression of HAC remained unaltered. The results obtained from qPCR data are concordant with our transcriptomics data, confirming the reliability of the results with minimum discrepancies.

Cellulose Degrading Enzyme Production in Response to 2-Deoxyglucose are Results of Shunting of Catabolite Repression
A. malaysianum has an innate ability to secrete lignocellulosedegrading enzymes in the extracellular medium ( Table 1A). The higher level of production of these enzymes makes A. malaysianum an industrially important strain for biofuel production that can be rapidly exploited. In vitro addition of 2-DG in cellobiose medium resulted in further increase of β-glucosidase and endoglucanase activity by several folds (Table 1B). It may be proposed that such rise in β-glucosidase results from the combined action of two mechanisms acting simultaneously; (i) Shunting down of the catabolite repression of the cell by glucose anti-metabolite 2-DG in the medium; (ii) cellobiose induces lignocellulolytic enzymes. The nonmetabolizable glucose analog 2-DG has been used extensively to induce catabolite repression in bacteria and fungi (Kornberg and Lambourne, 1994). Moreover, the use of 2-deoxyglucose (2-DG), as an antimetabolite has been routinely practiced for selection of the carbohydrate degrading enzyme hyper-producing mutants. The addition of 2-DG to cells causes a glucose starvation-like response, inhibiting growth and reducing viability even in the presence of abundant glucose. 2-DG is taken up and converted to 2-deoxyglucose-6-phosphate (2-DG-6P). However, the absence of a hydroxyl group on C-2 prevents the further catabolism of 2DG-6P by phosphogluco-isomerase. Accumulation of 2DG-6P in the cytoplasm result in inhibition of growth by repressing early glycolytic enzymes, thereby inhibiting glycolysis. Whether these are the only means by which 2-DG short-circuits normal glucose utilization needs to be explored further.
Exogeneous addition of 2-DG in a growth media containing a sole carbon source extended the lag phase, and inhibited growth in fungi (Zhu et al., 2015). In presence of preferred carbon source (glucose), sugar uptake and utilization was unaffected and maximum growth was observed in 48 h in absence as well as in presence of 2-DG. However, in cellobiose medium 2-DG caused delayed growth, indicating that 2-deoxyglucose led to CCR of cellobiose utilization (Figures 1A,B). Similar result was observed in Thermoanaerobacterium aotearoense where 0.5 mg/ml of  2-DG strongly inhibited growth (Zhu et al., 2015). The CCR is generally controlled by the second messenger cyclic adenosine 3' ,5'-monophosphate (cAMP), which has a fundamental role in global gene regulation (Lin et al., 2013). Catabolite repression by 2-DG was not reversed by adding exogenous cAMP (10 mM; Figure 1C). Similar studies on another Ascomycetous Aspergillus nidulans also corroborates our observation of no CCR relief in presence of cAMP (Borgia and Sypherd, 1977;Lee et al., 1996). Sample spreadsheet of data analysis using the 2 −∆∆CT method. The fold change in expression of the target gene relative to the internal control gene (-actin). The samples were analyzed using qPCR and the Ct data were imported into Microsoft Excel. The mean fold change in expression of the target calculated where ∆∆C T (C T , Target -C, Actin ) treated-(C T,Target -C, Actin ) control.
Since 2-DG hampers microorganism's growth, it is widely used as an inducer of CCR (Zhu et al., 2015). We added different concentration of 2-DG (0.25-1 mg/ml) to two different carbon source (1% cellobiose and glucose media) to detect the CCR by monitoring the culture growth. 2-DG delayed growth in cellobiose but not in glucose cultures. In glucose medium maximum growth (measured as mycelia dry weight) was observed in 48 h both in absence and presence of 2DG ( Figure 1A). The growth rate measured in absence of 2-DG was 0.31 mg/g and 0.22 mg/g at 1 mg/ml concentration. In cellobiose medium maximum growth in absence of 2DG was observed at 60 h (0.271 mg/g) and was lagged upto 36 h with increasing concentration of 2DG (0.25-1 mg/ml). The growth density was also visibly reduced with increasing 2DG concentration (Figures 3A,B). Addition of 10 mM cAMP showed no effect on the growth delay caused by 2-DG in cellobiose medium (Figure 1C), unlike in E. coli where exogenous addition of cAMP in the medium provided CCR relief (Ullmann and Monod, 1968;Isaacs et al., 1994). This proves that in A. malaysianum, 2-DGinduced CCR is cAMP-independent.

Release of Extracellular β-Glucoside in Response to 2-DG is Involved in Cellular Signaling
In filamentous fungus, signal transduction pathways have been found to be largely associated with cell wall weakening (Wang et al., 2014). In N. crassa addition of 1 mM cellobiose and 2deoxy-glucose inhibited the secretion of β-glucoside. Although the actual mechanism behind this was not discussed it was hypothesized that 2-DG was nearly as effective as glucose in blocking induction. Segregated and pellet like hyphae is probably the consequence of hyphae fragmentation due to production of cell wall hydrolytic enzymes like endochitinase and endo β-1,3 glucanases (Szilágyi et al., 2012(Szilágyi et al., , 2013 (Figures 3A-D).
Since 2-DG is a non-metabolizable glucose analog and also a glycosylation inhibitor, it markedly reduced the amount of cell wall components resulting in declining dry cell mass and progressive fragmentation of the hyphae in A. malaysianum (Figures 3E-H, 10). Cell wall of S. cerevisae grown in presence of 0.2%, 2-DG contained less glucose and mannose compared to the normal cell (McCartney et al., 2014). In the present study, we gained a more complex view of the changes taking place in the cell wall homeostasis. We detected the induction of several genes encoding hydrolytic enzymes such as chitinase, glucan β-1,3glucosidase (Table 1B) and down regulation of chitin synthase in 2-DG treated condition (Figure 6; Table 4). Similar induction of several genes encoding cell wall hydrolyzing enzymes (chitinase, α and β-glucanases and mannanases) has been detected in carbon starved cultures of Aspergillus niger (Nitsche et al., 2012). Up regulation of putative β-glucosidase encoding genes might also contribute to the degradation of the fungal cell wall, as observed by enhanced extracellular β-glucosidase activity in presence of 2-DG (Figure 2; Table 3; Szilágyi et al., 2010).

Distinctive Expansion of Gene Families
Associated with Plant Cell Wall Degradation, Secondary Metabolism, and Nutrient Uptake in the Transcriptome of A. malaysianum The expanded CAZyme arsenals of A. malaysianum were similar to those of Eutypa lata, Thielavia terrestris, and Colletotrichum sp. (Table 6). Interestingly, these fungi are known to be pathogen on host plants and have been isolated as endophytes from other sources also. Glycoside hydrolases  (Zhao et al., 2013), which possesses high numbers of pectolytic enzymes, A. malaysianum showed a wider array of enzymes that target cellulose and hemicellulose. This includes endo-β-1,4-cellulases (GH5), β-glucosidases (GH3), xyloglucan transglucosylase/hydrolases (GH16), and β-xylosidases (GH43). We found significant expansion of genes coding for CAZymes containing the CBM1 domain, a carbohydrate-binding module. CBM1 strongly binds to crystalline cellulose and that is required for full activity of fungal cellulases similar to the findings from E. lata that has a similar CBM1gene family expansion pattern (Klosterman et al., 2011).
P450s have a broad spectrum of functions in fungi, from housekeeping activities, such as synthesis of essential membrane lipids, to synthesis of secondary metabolites, and detoxification of xenobiotic compounds . Differences in the number of genes belonging to the different P450 classes were observed among species, particularly between the Ascomycetes and Basidiomycete (Brakhage, 2013). A. malaysianum harbors high copies of antifungal, antioxidant, monooxygenases and cytochomeP450 superfamily genes including the genes responsible for the production of rare alkaloid-Loleine (Wang et al., 2008; Table 3). Three different types of CYP families were identified in A. malaysianum: cytochrome c, b560, and the P450. The high diversity of secondary metabolites is related to the diversity of the CYP genes. For example, the 219 CYP genes in the Ganoderma lucidum genome resulted in a large number of different secondary metabolites .
Pathways associated with Hexokinase, Transaldolase, F-6-phosphate, and glucose-6 phosphate dehydrogenase and alcohol dehydrogenase production were induced whereas Isocitrate lyase and Phosphogluco-isomerase were suppressed (Figures 8A,B). The activity and expression fold of these marker metabolic enzymes was quantified by enzymatic assay and also determined by qPCR data reinforcing our suggested carbon flow cycle (Livak and Schmittgen, 2001; Figure 10). The accumulation of non metabolizable 2DG-6-P in the cytoplasm forced the cell to recover from the glucose sink by entering into Pentose phosphate Pathway and gluconeogenesis pathway. This leads to up regulation of several metabolic enzymes like transaldolase (talB), transketolase (tktA), fructose bis phosphate (fbp), glucose 1 dehydrogenase (gdh), gluconolactonase (gnl), ribose 5 phosphate isomerase (rpiB; Figure 6).
CreA/Cre1-knock out mutants display severe phenotypic changes such as reduced growth, abnormal hyphal morphology and sporulation. Up regulation of different genes linked to carbohydrate metabolism (e.g., glycoside hydrolases and glycoside transferases), cellular transport, signaling molecules like MAPK and SNF2 as well as certain chaperons and heat shock protein that is induced by several stresses, including glucose limitation (Portnoy et al., 2011). It may be presumed that 2-DG works in a manner homologous to these Cre mutants creating a similar environment and displays similar traits (Figure 10).

CONCLUSION
Tricking A. malaysianum into producing industrially important enzymes with 2-Deoxy D-Glucose treatment is a novel way to increase commercial usability of this fungi. The ability of A. malaysianum to produce hemicellulolytic enzymes renders it an attractive organism for the production of fuels and renewable energy. We investigated 2-DG induced carbon metabolism and carbohydrate degrading enzymes secretion in this organism in order to better understand the regulatory mechanism of this non metabolized sugar and how to manipulate it. Transcriptomics study and our experimental data indicate that A. malaysianum is regulated in a CCR independent mechanism in presence of 2DG. The induced synthesis of secondary metabolites, fatty acid and lipids, and proteins associated with cell wall degrading enzymes demonstrated the benefits of addition of 2-DG in the culture medium. 2-DG induces a carbon starving like environment causing increased production of proteins and extracellular enzymes. As a result, endoplasmic reticulum activity increases considerably leading in further up regulation of transporters, disulfide isomerases, and binding proteins. Plethora of secretory enzymes linked to carbohydrate degradation and

DATA AVAILABILITY
Sequence data is publicly available in Genbank with SRA id SRR746797. The Bioproject id is PRJNA186407 and Biosample id is SAMN01923006.