Transcriptome and secretome analysis of Aspergillus fumigatus in the presence of sugarcane bagasse

Background Sugarcane bagasse has been proposed as a lignocellulosic residue for second-generation ethanol (2G) produced by breaking down biomass into fermentable sugars. The enzymatic cocktails for biomass degradation are mostly produced by fungi, but low cost and high efficiency can consolidate 2G technologies. A. fumigatus plays an important role in plant biomass degradation capabilities and recycling. To gain more insight into the divergence in gene expression during steam-exploded bagasse (SEB) breakdown, this study profiled the transcriptome of A. fumigatus by RNA sequencing to compare transcriptional profiles of A. fumigatus grown on media containing SEB or fructose as the sole carbon source. Secretome analysis was also performed using SDS-PAGE and LC-MS/MS. Results The maximum activities of cellulases (0.032 U mL-1), endo-1,4-β--xylanase (10.82 U mL-1) and endo-1,3-β glucanases (0.77 U mL-1) showed that functional CAZymes (carbohydrate-active enzymes) were secreted in the SEB culture conditions. Correlations between transcriptome and secretome data identified several CAZymes in A. fumigatus. Particular attention was given to CAZymes related to lignocellulose degradation and sugar transporters. Genes encoding glycoside hydrolase classes commonly expressed during the breakdown of cellulose, such as GH-5, 6, 7, 43, 45, and hemicellulose, such as GH-2, 10, 11, 30, 43, were found to be highly expressed in SEB conditions. Lytic polysaccharide monooxygenases (LPMO) classified as auxiliary activity families AA9 (GH61), CE (1, 4, 8, 15, 16), PL (1, 3, 4, 20) and GT (1, 2, 4, 8, 20, 35, 48) were also differentially expressed in this condition. Similarly, the most important enzymes related to biomass degradation, including endoxylanases, xyloglucanases, β-xylosidases, LPMOs, α-arabinofuranosidases, cellobiohydrolases, endoglucanases and β-glucosidases, were also identified in the secretome. Conclusions This is the first report of a transcriptome and secretome experiment of Aspergillus fumigatus in the degradation of pretreated sugarcane bagasse. The results suggest that this strain employs important strategies for this complex degradation process. It was possible to identify a set of genes and proteins that might be applied in several biotechnology fields. This knowledge can be exploited for the improvement of 2G ethanol production by the rational design of enzymatic cocktails. Electronic supplementary material The online version of this article (10.1186/s12864-018-4627-8) contains supplementary material, which is available to authorized users.


Background
The demand for energy has increased continuously worldwide, which has raised concerns about sustainability and has prompted the search and development of advanced renewable and sustainable sources of energy [1]. Bioethanol has been noted as an alternative fuel to tackle these issues [1][2][3]. In Brazil, ethanol production relies on the fermentation of sucrose from sugarcane to yield the so-called first-generation (1G) bioethanol [4,5]. The current Brazilian production is estimated at 30 billion liters per year, but the growing appeal of this fuel has called for investments in the development of new technologies to produce ethanol [6]. Large amounts of sugarcane straw and bagasse are generated every year in Brazil, so this biomass could be used as a substrate to produce 2G bioethanol, which in a few years will compete with 1G ethanol costs [6][7][8][9][10][11].
Previously, specific cellulose-, hemicellulose-, pectin-, and lignin-degrading enzymes were identified as secreted by A. fumigatus in the presence of different carbon sources (Avicel, cellulose, rice straw, starch, xylan, corn and soybean) that can be used in the lignocellulosic bioenergy industry [22,[24][25][26]. To gain more insight into how efficiently A. fumigatus AF293 can depolymerize the sugarcane bagasse, a complex biomass important for Brazilian 2G ethanol production, and to identify genes and proteins responsible for these lignocellulosic breakdown reactions, we examined the transcriptional response by RNA-Seq and proteomic profile by mass spectrometry (LC-MS/MS) of A. fumigatus that was cultivated in the presence of steam-exploded sugarcane bagasse (SEB).

Strains, media, and growth conditions
A. fumigatus AF293, gently donated by the Prof. Dr. Sérgio Akira Uyemura (University of São Paulo, BR), was grown on YAG medium (2% (w/v) dextrose, 0.5% (w/v) yeast extract, 0.1% (v/v) trace elements and 1.8% (w/v) agar) at 37°C for two days. Spores were harvested and inoculated to a final concentration of 1 × 10 8 per 50 mL of YNB culture with 1% (w/v) fructose as the carbon source at 37°C for 16 h (h) in a rotary shaker with agitation at 200 rpm. Afterward, the mycelia were transferred to 1% (w/v) SEB (47.5% cellulose; 9.0% hemicellulose and 34.3% lignin) or 1% (w/v) fructose as the carbon source for 3, 6, 12, 18, 24, 48 and 72 h. Fructose was used as a control in all experimental conditions [26]. Mycelia were harvested by filtration through Whatman grade 1 filters (GE Healthcare, Grandview Blvd. Waukesha, WI, USA), washed once with sterile cool water and kept at − 80°C until RNA extraction. Supernatants were collected to measure enzymatic activity and determine the secretome. All the experiments described below were performed in three biological replicates.
Enzymatic activities were also measured by the dinitrosalicylic acid (DNS) assay [27]. Cellulase activities were measured with β-glucan and low-viscosity carboxymethylcellulose (CMC) as substrates, and xylanase activities were measured with the xyloglucan. Briefly, 20 μL of the supernatant from the samples grown in presence of 1% SEB for 24, 48 and 72 h were mixed with 30 μL of sodium acetate buffer 100 mM (pH 5.5) and 50 μL of substrate at 0.5% (w/v) final concentration to achieve a final volume of 100 μL. The reactions were incubated at 40°C for 5 min for β-glucan and for 10 min for xyloglucan and CMC substrates. The reaction was stopped by adding 100 μL of DNS. All the reactions were performed in triplicate. The calculation of enzyme activities was based on a corresponding standard containing glucose. One unit (U) of enzymatic activity was defined as the amount of enzyme needed to liberate 1 μmol of reducing sugars per minute.

RNA isolation and cDNA synthesis
Fungal biomass was harvested at different times from SEB or fructose culture conditions, and mycelia were ground in liquid nitrogen using a mortar and pestle. Total RNA was purified by using the "Direct-zol™ RNA MiniPrep" kit according to the manufacturer's instructions (Zymo Research, Irvine, CA, EUA) using the oncolumn DNAse treatment. RNA integrity was confirmed with a bioanalyzer by using the "Agilent RNA 6000 Nano" kit (Agilent Technologies, Santa Clara, CA, EUA) and the "Plant Total RNA Nano" protocol. RNA was quantified on a Qubit® 2.0 fluorimeter (Thermo Fisher Scientific, Waltham, MS, EUA) with the Qubit® RNA BR Assay kit (Thermo Fisher Scientific, Waltham, MS, EUA). cDNA was synthesized from 1 μg of mRNA using SuperScript® II Reverse Transcriptase (Invitrogen, Carlsbad, CA, EUA).

Library preparation and RNA sequencing
RNA sequencing libraries were prepared using the "Tru-Seq Stranded mRNA HT Sample Prep" kit (Illumina, San Diego, CA, EUA), mRNA enrichment was performed using magnetic beads coupled with oligo (dT). Sequencing was carried out in the HiSeq 2500 system (Illumina, San Diego, CA, EUA) at the NGS facility located at the Brazilian Bioethanol Science and Technology Laboratory (CTBE), Campinas, SP, Brazil.

Bioinformatic analysis of RNA-Seq data
FastQC [28] was used to check the quality of the sequencing reads visually. Removal of the remaining adapter sequences and quality trimming with a sliding window of size 4, minimum quality of 20, and length filtering (to keep reads with a length of at least 60 bp) were performed with Trimmomatic v0.32, [29]. Clean reads were screened against a database of ribosomal RNA with the aid of SortMeRNA [30]. High-quality reads were mapped in a strand-specific manner by using TopHat2 [31] against the genome sequence of A. fumigatus Af293 obtained from ASPGD [32,33]. The number of exon-exon junctions at different levels of read subsampling was employed to confirm sequencing saturation with RSeQC [34]. Mapping of the reads to the features of the exons were summarized at the gene level by using the function featureCounts from the Rsubread v1.12.6 package [35] in R v3.0.2 [36] and the annotation file in GFF3 format from ASPGD. Differential gene expression was analyzed with edgeR [37] in R [36]. Briefly, genes with at least one CPM (counts per million) in at least three samples were kept for analysis, which was equivalent to removing genes with low and noisy expression. The expression values were normalized by the trimmed mean of M-values (TMM) method to account for differences in the composition of RNA [38]. After the dispersion was estimated and the biological coefficient of variation was computed, the differentially expressed genes were called by fitting a negative binomial model with generalized linear models (GLS) that included factors for the TMM and the dispersion estimates [37]. A likelihood ratio test was performed to provide a p-value for differential expression. The p-values were adjusted for multiple testing by the method of Benjamini-Hochberg, to control the false discovery rate (FDR) [39]. The full R script used for the analysis and the raw count matrix are available in Additional file 1: Figure S1 and Additional file 2: Table S1, respectively. Genes with FDR values lower than 0.05 and log2-fold changes greater than 1.0 or lower than − 1.0, i.e., a difference of twice the expression level in either direction, were considered differentially expressed.

Supernatant analysis by SDS-PAGE
Supernatants (50 mL) from SEB or fructose culture conditions were lyophilized until completely dry and resuspended in 2 mL of buffer (Tris-HCl 50 mM, pH 6.8; 1 mM DTT; and 1 mM protease inhibitor), and 15 μL was separated by 10% SDS-PAGE (110 V, 90 min). The proteins were visualized by staining with 0.1% Coomassie Brilliant Blue R250 (w/v), which was followed by destaining with 45% methanol and 10% acetic acid solution (v/v). The protein concentration was determined by Bradford's Assay (Bio-Rad Protein Assay Hercules, CA, EUA) [40]. Prior to mass spectrometry, all the bands from the SDS-PAGE gels were manually excised, reduced, alkylated, digested with trypsin, and purified (Promega, Madison, WI, EUA -V5111) according to a previously described method [41].

Identification of proteins by coupled system of the LC-MS/MS type
Peptides were sequenced on a Synapt G2 HDMS (Waters, Milfords, MS, EUA) mass spectrometer coupled to a UPLC NanoAcquity system with 1D technology (Waters, Milfords, MS, EUA) and captured by a C18 Symmetry column (5 μm, 180 μm × 20 mm) (Waters, Milfords, MS, EUA). The peptides were separated by using a 2-90% acetonitrile gradient in 0.1% formic acid and an HSS T3 analytical column (1.8 μm, 75 μm × 100 mm) (Waters, Milfords, MS, EUA) with a flow of 300 μL min − 1 for 120 min. The data were acquired on a Waters Synapt G2S Q-TOF mass spectrometer equipped with a NanoLock-Spray (Waters, Milfords, MS, EUA). The experiments were performed in the HDMSE mode (data-independent analysis). The mass spectra were processed with the Pro-teinLynxGlobalServer (PLGS) software version 3.1. The proteins were identified by comparison to the Aspergillus UNIPROT database (207,966 proteins) [42]. The defined parameters were automatic tolerance for precursors and ion products, minimum of three corresponding ion fragments per peptide, minimum of seven corresponding ion fragments per protein, trypsin missed cleavage, carbamidomethylation as a fixed modification, oxidation of methionine as a variable modification, and 4% FDR peptide.
For CAZy enzyme identification, the proteins in the secretome were screened with a library of hidden Markov models by using HMMER3 [46] of carbohydrateactive enzymes obtained from dbCAN [47]. Hits were considered positive on the basis of the dbCAN recommendations.

qRT-PCR analysis
After RNA-Seq analysis, 4 DEGs (Differentially Expressed Genes) were selected, including sugar transporters and CAZymes, for qRT-PCR analysis. RNA was extracted and purified as previously described. cDNA was synthesized from 5 μg of RNA using SuperScript® II Reverse Transcriptase (Invitrogen, Carlsbad, CA, USA). Quantitative PCR (qPCR) analyses were performed according to Semighini et al. [48]. The abundance of the respective mRNAs was normalized using β-tubulin probes. The primers for the investigated genes are listed in Additional file 3: Table S2.

Functional enrichment
Genes identified as differentially expressed were analyzed by FunCat functional enrichment [49]. The CAZy proteins from the secretome were classified according to the GO-Slim classifications from the AspGD based on the ontology "Molecular function" (GO Categorization Slim Mapper) [50].

Venn diagrams
The area-proportional Venn diagrams were drawn based on images generated with free online software [51].

Enzymatic analysis
To evaluate the activity of enzymes produced by A. fumigatus in the presence of sugarcane bagasse, we performed enzymatic assays using xyloglucan, β-glucan, and CMC as substrates. Specific endo-1,4-β-xylanase and endo-1,4-βglucanase activities, were also investigated using Azo-Xylan (Birchwood) and Azo-CM-Cellulose (both from Megazyme International, Bray, Ireland), respectively. The enzymes from supernatants derived from A. fumigatus cultures were capable of hydrolyzing cellulose (CMC) and hemicelluloses (β-glucan and xylan), but no activities were detected in xyloglucans ( Fig. 1a-b).
We observed that the activities depended on duration of growth, with maximum activities detected within a b Fig. 1 Enzyme activity. Enzymatic activities (U mL −1 ) of supernatants from A. fumigatus culture against different substrates after 24, 48 and 72 h growth on sugarcane bagasse. Each bar represents the mean and the standard deviation of values from three independent experiments. a Azo-xylan, β-glucan and xyloglucan; b 1% CMC and Azo-CMC three days, except for endo-1,4-β-D-xylanase, with peak activity after two days culture and no activities detected prior to 24 h of incubation. Maximum activities of cellulases (0.032 U mL − 1 ), endo-1,4-β-xylanase (10. 82 U mL − 1 ) and endo-1,3-β glucanases (0.77 U mL − 1 ) were detected for A. fumigatus while in A. niger (0.002 U mL − 1 ; 2.3 U mL − 1 and 0.4 U mL − 1 , respectively) and T. reesei RUT-C30 (0.0039 U mL − 1 , 0. 4 U mL − 1 and 0.1 U mL − 1 , respectively) are described on the same biomass [4]. These results indicate that A. fumigatus is an excellent producer of an arsenal of hydrolytic enzymes, with activities superior to the hypercellulolytic strain T reesei RUT30-C.
To gain more insight into the hydrolytic enzymes of A. fumigatus specific for sugarcane bagasse breakdown, we selected the cultivation of 24 h (when we detected enzyme activity) because we were interested in the initial process of SEB breakdown. In this time we determinate the transcriptome and the secretome responses of this strain.

Analysis of the transcriptome of A. fumigatus under the influence of sugarcane bagasse as the substrate
To identify potential new genes involved in SEB breakdown, we analyzed the transcriptome by RNA-Seq after 24 h of cultivation. After RNA sequencing, each sample generated approximately 14 to 16 million paired-end reads. RNAseq data were analysed by comparing the mycelium grown on SEB and that grown on fructose. We observed 2227 genes differentially expressed (FDR < 0.05, |log2FC| > 1) in SEB where 1181 were upregulated, while 1045 were upregulated in fructose conditions (downregulated in SEB) (Additional file 4: Table S3). Gene ontology (GO) and the functional catalogue (Fun-Cat) classified the differentially expressed genes functionally in 18 different enriched categories [40]. Two significant categories among upregulated genes were Metabolic Processes (GO:0008152) and Protein Synthesis (GO:0006412) and in downregulated genes Metabolic Processes (GO:0008152) and Energy (GO:0006112) (Additional file 5: Table S4). Genes related to the regulation of C-compound and carbohydrate metabolic processes represent the two GO terms commonly enriched for both up-and downregulated genes in the SEB condition, including genes encoding carbohydrateactive enzymes (CAZymes) and transporters. Given the importance of CAZymes for the degradation of biomass, we directed our efforts toward a better understanding of the transcriptional profile of these enzymes and sugar transporters.
We found 197 differentially expressed CAZyme genes classified based on the CAZy database (http://www.cazy. org) [52]. Concerning the 566 CAZyme genes predicted in the A. fumigatus AF293 genome categorized into the different classes (247 GHs, 105 GTs, 96 CEs, 59 AAs, 15 PLs, and 44 CBMs) (Additional file 6: Table S5), we concluded that 35% of CAZyme genes were differentially expressed in our data, which highlights the potential that a wide spectrum of hydrolytic enzymes were produced. However, glycosyl transferases appeared only in a very small percentage (~1%), suggesting a secondary role in polysaccharide degradation (Fig. 2a).
Taken together, endoglucanases, cellobiohydrolases and beta-glucosidases were significantly upregulated, suggesting the cellulose degradation potential of this strain, and the abundance of hemicellulases highlights, once again, the great potential of A. fumigatus in complex biomass deconstruction.
To validate RNAseq data and get additional information about the expression over time, we have performed qRT-PCR for 4 DEGs that encode enzymes essential to biomass degradation, Afu4g07850 (LPMO) , Afu1g14710 (β-glucosidase), Afu6g11610 (1,4-β-Dglucan-cellobiohydrolase) and Afu3g02090 (β-xylosidase), during 3, 6, 12, 18 and 24 h of cultivation in SEB and fructose. The expression profiles of these genes behaved in different ways: Afu1g14710 and Afu3g02090 genes were strongly induced at the beginning (6 h) of the growth in SEB, and their expression decreased after 6 h, while Afu4g07850 had an increasing gene expression during the time course, and Afu6g11610 increased at 24 h (Fig. 3).

Sugar transporters identified during RNA sequencing
Approximately 106 genes encoding sugar transporters have been reported in the Aspergillus genome, and only 88 genes were described as encoding sugar transporters in A. fumigatus strain Z5, which are distributed among the SP, FHS, SHS, and GPH families (the SP family includes 79 genes) [48]. Additionally, the genomes of filamentous fungi also encode large numbers of major facilitator superfamily (MFS) transporters. Among them, 25 transporters were differentially expressed on SEB, classified as encoding MFS hexose transporter, MFS and sugar transporter, UDP-Glc/Gal endoplasmic reticulum nucleotide sugar transporter, nucleotide sugar transporter, hexose transporter protein, high affinity glucose/hexose transporter, MFS glucose transporter, MFS lactose transporter, MFS maltose transporter, and xylose transporter (Fig. 4a).  Orthologous genes in Aspergillus and non-Aspergillus species were identified by sequence analysis according to Aspergillus Genome Database (AspGD; http://www. aspgd.org) [32]. Three orthologous genes encoding possible putative xylose transporters (Afu1g03530, Afu4g14610, and Afu6g14442) and five related to cellobiose transporters (Afu3g01670, Afu6g14500, Afu6g14560, Afu7g05100, and Afu8g04480) have been identified in N. crassa, A. oryzae, A. niger, and A. nidulans (Additional file 7: Table S6).

Characterization of the secretome of A. fumigatus in the presence of sugarcane bagasse
Once the transcriptome was characterized, we analyzed the secreted protein profiles of A. fumigatus cultivated in the same condition by SDS-PAGE and LC-MS/MS. The total protein secreted by the fungi was approximately 300 μg mL − 1 in SEB versus 112 μg mL − 1 in the fungi grown on fructose (Additional file 8: Figure S2). In the SEB supernatant, we detected 128 secreted proteins, and only 44 were detected in the fructose supernatant (Additional file 9: Table S7), 27 of which are the same for both conditions (Fig. 5a). As previously described, in the A. fumigatus genome, 566 CAZyme genes (461 proteins, excluding GTs) are predicted [25]; 271 of them are predicted to be secreted (Additional file 6: Table S5). The presence of a signal peptide in these proteins was inferred using Sig-nalP 3.0 (http://www.cbs.dtu.dk/services/SignalP/) [44]. Approximately 78 proteins identified in SEB correspond to CAZymes in the classes GH (46), CE (8), AA (8), PL (2), and CBM (15), which is 18.61% of the total CAZymes encoded by the A. fumigatus genome. The remaining proteins are classified as proteins of unknown functions (17%), hydrolases/peptidases/proteases/binding proteins (14%), oxidoreductases (6%), transferases (2%), and lyases (2%) (Fig. 5b). We have also identified intracellular proteins, suggesting cell lysis or an unknown mechanism.
We have also identified lignin-depolymerizing enzymes such as catalase-peroxidase, cellobiose dehydrogenase, catalase B, FAD-dependent oxidase, laccase, and Cu-Zn superoxide dismutase. Although the sugarcane bagasse employed here was treated by steam explosion, traces of lignin might have remained in the substrate, which justifies the secretion of these enzymes by the fungus.
We can conclude that the most important CAZymes  Proteins identified in this study, which have also been identified in a previous study [23] b Proteins identified in this study, which have also been identified in a previous study [22] c Proteins identified in this study, which have also been identified in a previous study [24] GH93, CE1, CE16 and AA9 (LPMO)) were secreted and play important roles in biomass degradation. For the first time, new proteins such as GH16 (endo-1,4-betaglucanase), GH5 (endoglucanase), LPMO (AA9), swollenin and GH3 (β-glucosidase) were identified in the Aspergillus fumigatus secretome, probably because we used sugarcane bagasse as the source of carbon, and these enzymes can be specific to this complex biomass.

Integration of secretomics and transcriptomics
We observed weak correlations between transcriptome and secretome datasets, mainly because we chose the same time (24 h) to isolate mRNA and proteins. Considering that Aspergillus needs at least a few hours to translate mRNA to protein and to secrete it, these data provide an idea about which proteins are transcribed earlier or produced constitutively. Among the 1181 upregulated genes in transcriptome, 63 encoded proteins were detected in the secretome. As the same way, 16 of secreted proteins were identified as downregulated in RNAseq data (Fig. 6). In addition, the weak correlation observed could be the result of the influence of some factors that alter transcription and translation mechanisms [22,23,54].
The both data revealed a significant upregulation of secreted CAZymes, which is important for the observation of specific alterations triggered by different conditions. Among the 135 upregulated CAZyme genes, 48 encoded proteins were detected in the secretome, mainly cellulolytic and hemicellulolytic enzymes (Table 3). Four enzymes, β-1,3-endoglucanase EglC (Afu3g00270), FADdependent oxygenase (Afu3g00840), FAD/FMN-containing isoamyl alcohol oxidase MreA (Afu6g03620) and oligopeptidase family protein (Afu8g04730), were downregulated when the strain was cultivated in SEB but were detected in the secretome.
The variation in expression during fungal growth in the presence of SEB was observed for some genes in the qRT-PCR data (Fig. 3). We observed that the gene expression depends on the duration of incubation, which explains the low correlation between the data of the secretome and transcriptome and supports the hypothesis that the approach used in this study is able to provide information about the modulation of gene expression of A. fumigatus in the presence of sugarcane bagasse. Genes encoding enzymes that are upregulated in the transcriptome (the highest log 2 FC) and are present in the proteome, as well, are described in Table 3. The main enzymes produced by the fungus are endo-β1,4-xylanase, cellobiohydrolase B, endoglucanase and arabinofuranosidase, which leads us to speculate that A. fumigatus spent its energy on the transcription and production of cellulolytic and xylanolytic compounds, which were probably being secreted for biomass degradation and the release of small sugars.
To gain more insight into new enzymes and to identify new genes specific to sugarcane biomass hydrolysis, we have chosen to investigate A. fumigatus because it is widely distributed in the environment, it can degrade plant biomass, and it is an excellent enzyme producer [21]. To the best of our knowledge, this is the first study on transcriptional response and secretome of A. fumigatus grown on sugarcane steam-exploded bagasse.
There are few omic studies on biomass hydrolysis by A. fumigatus [22][23][24][25], and only one study has analyzed the transcriptome of A. fumigatus when grown Fig. 6 Intersection of differentially expressed genes and secreted proteins. Venn diagram of transcriptome (upregulated genes and downregulated genes) and proteins secreted by A. fumigatus when grown on SEB on polysaccharide substrates. Miao et al. [25] conducted a transcriptional study on the induction of CAZymes by this fungus grown on cellulose, oat spelt xylan, rice straw and sucrose. The authors showed that important genes are differentially expressed in each carbon source.
Here, we found few discrepancies in the number of induced genes when we used SEB as carbon source ( Table 4). The main CAZyme families (GH1, GH3, GH5, GH6, GH7, GH10, GH11, GH12, GH43, GH62, GH67, GH74, AA9, CE3, CE5 and CE16 [85]) were also induced. However, important genes including GH45, GH51, GH54, GH93, GH115, and CE1 were downregulated in SEB. Furthermore, two important genes were exclusively induced in SEB: PL4 (Afu4g03780 and Afu8g00820), which plays a role in pectin breakdown, and CBM35 (Afu8g02510), which is known to bind primarily to xylan and mannans [86]. The distinct gene expression was most likely due to substrate composition and to cultivation time, which was 24 h for our analysis. Shorter cultivation times could point out new pattern to gene induction. In qRT-PCR data, we observed that 6 h  is the beginning of gene expression, which might represent a standard mechanism in which A. fumigatus acts in contact with complex biomass sources and should contain more highly induced enzymes. These results also explain the percentage (35%) of differentially expressed CAZymes identified in the RNA-Seq data at 24 h. Because glycosyltransferases (GTs) contribute to fungal cell remodeling, the percentage of upregulated genes was low (~1%). Likewise, Miao et al. [25] described that genes encoding GTs are downregulated in Z5 strain, which supports the idea that GTs do not directly participate in the hydrolysis of complex biomass.
In this sense, the biomass itself has to be investigated in order to understand sugarcane biomass hydrolysis as well as possible. A similar work performed by Borin et al., 2017 [18] described a transcriptional response of A. niger and T. reesei grown in SEB for different periods. They found 190 upregulated CAZymes from 62 different families in A. niger, and 105 genes of 51 CAZyme families in T. reesei, whereas we detected 135 upregulated CAZymes in the A. fumigatus transcriptome herein. The number of genes induced by each microorganism was different and depended on time. A higher number of DEGs in A. niger and T. reesei was observed in 24 h of culture, and so we then compared our data at this duration, and again small differences in upregulated CAZymes were observed.
After biomass hydrolysis, which breaks down cellulose and hemicellulose into mono-or disaccharides, the released sugars need to be transported into the cells through a large number of sugar transporters, most of which have not been characterized yet [20,21,87]. One of the main challenges concerning biofuel production from lignocellulosic biomass is the inability of organisms to grow on, to transport, and to ferment sugars other than glucose (e.g., xylose and cellobiose). Gaining a better insight into potential xylose and/or cellobiose transporters seems to be a good approach to overcome this challenge [87]. These transporters represent an important industrial tool that can be applied to different industrial processes [88][89][90]. Additionally, the genomes of filamentous fungi also encode large numbers of major facilitator superfamily (MFS) transporters. For example, the T. reesei and A. nidulans genomes have been predicted to encode 164 and 357 proteins belonging to MFS, respectively, although the exact number of proteins involved in sugar transport remains unknown [91][92][93].
We were also interested in new sugar transporters, such as the xylose transporter. Until now, no sugar transporter for A. fumigatus related to biomass breakdown has been described. We verified that the 25 DEG homolog transporters had particular expression profiles, upregulated or downregulated, suggesting that SEB hydrolysis released enough glucose, xylose or cellobiose to regulate sugar transporter gene expression. Another interesting finding was that two sugar transporters (Afu6g14442 and Afu1g03530) were highly induced in SEB (25.5 times and 10 times, respectively) and in the presence of 1% xylose, which revealed that these sugar transporters could be specific xylose transporters in Aspergillus fumigatus. Overexpression of xylose transporters in S. cerevisiae is a fast way to use xylose and may improve ethanol productivity [94].
These data showed that, although the transcriptome data did not reveal potential new enzyme targets for the deconstruction of sugarcane biomass deconstruction, the secretome analyses indicated key enzymes that may be essencial for this hydrolysis and which act synergistically for efficient deconstruction. We clearly observed the need for accessory enzymes secreted as LPMO and swollenin, which have never been described in other A. fumigatus secretome analyses [95], which allowed us to conclude that the other secreted CAZymes together with AA9 identified in A. fumigatus form a potential arsenal of hydrolytic enzymes. LPMOs have recently been implicated in lignocellulosic biomass degradation. Although these enzymes were first classified into the GH61 and CBM33 families, they are currently classified into the AA9 and AA10 families. These enzymes cleave the lignocellulosic biomass glycosidic bonds through an oxidative mechanism that provides new ends for the recognition of cellulases and for action on cellulose [95,96]. In addition, AA9s have been identified in A. niger, M. thermophile, T. asperellum, and T. reesei secretome growth in sugarcane bagasse, which allows us to conclude that they play an important role in lignocellulosic biomass breakdown [4,[95][96][97][98][99].
Many studies have focused on LPMO enzymes, and some works have even characterized them, but no investigations into LPMO enzymes from A. fumigatus are available so far [96]. Hence, the role of most of these enzymes remains unclear, and AA expression during A. fumigatus growth on bagasse suggests that they play an important part in biomass degradation.
Together, the transcriptome and secretome have shown several enzymes that A. fumigatus uses to hydrolyze SEB and which most likely act synergistically to depolymerize cellulose and hemicellulose. In most Aspergillus species, distinct genes encode the same class of enzymes (isoenzymes) [77], as observed by the data regarding the A. fumigatus secretome and transcriptome. Our results suggest that complete hydrolysis of this lignocellulose biomass to simple sugars, such as glucose, xylose, and arabinose, requires the combined actions of several enzymes that have different substrate specificities and act synergistically. The great potential of this species is evident, and its enzymes can contribute to optimization of enzymatic cocktails for use in 2G bioethanol production.

Conclusion
Through these findings, it is suggested that different biomasses require a set of enzymes due its complexity and A. fumigatus Af293 is an excellent CAZymes producer for sugarcane biomass breakdown. The analysis of proteome and transcriptome revealed a set of CAZymes highly expressed and secreted, such as cellulases, hemicellulases, delignification and auxiliary enzymes necessary to SEB breakdown. In addition, from CAZymes proteins, LPMOs, which could contribute to better degradation of cellulose, were also detected in A. fumigatus secretome. Cellobiohydrolases, endoglucanases and LPMOs can act synergistically in cellulose depolymerization and LPMOs can be included in the most advanced enzymatic cocktails. Altogether, the data show that despite the pathogenicity of A. fumigatus, it can produce a wide variety of enzymes, which can be expressed in a nonpathogenic microorganism and may contribute to the optimization of currently marketed enzymatic cocktails for the viable production of 2G bioethanol.

Funding
This study was supported by Fundação de Amparo À Pesquisa do Estado de São Paulo (FAPESP, 2014/10466-0). AVB and LEG were funded by scholarships granted by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq). PFG was funded by a scholarship granted by Coordenação de Aperfeiçoamento de Pessoal de Ensino Superior (CAPES).

Availability of data and materials
The datasets generated and/or analysed during this study are included in this published article and its Additional files 1, 2, 3, 4, 5, 6, 7, 8 and 9. The raw data from A. fumigatus RNAseq are available at the NCBI's Gene Expression Omnibus under the accession code GSE95798 (https://www.ncbi.nlm.nih.gov/geo/query/ acc.cgi?acc=GSE95798), the underlying strand-specific raw short reads were deposited at the NCBI's Short Read Archive and are associated to the BioProject PRJNA376829 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA376829).
Authors' contributions PFG and AVB contributed substantially to the experimental work, data analyses and were involved in writing the manuscript and preparation of figures/tables, and offered intellectual contributions. LEG carried out secretome experiments and intellectual contributions. ESS carried out RNAseq experiments, data analyses and offered intellectual contributions. DMRP performed the analysis of RNAseq data, offered intellectual contributions and wrote the manuscript. SAU offered intellectual contributions at various stages of this work and provided critical reviews of the manuscript. TMD participated in the design and coordination of the study, wrote the manuscript, and conducted the final analysis of the data. All the authors have read and approved the final manuscript.
Ethics approval and consent to participate Not applicable.

Consent for publication
Not applicable.