# Carbon sources and XlnR-dependent transcriptional landscape of CAZymes in the industrial fungus Talaromyces versatilis: when exception seems to be the rule

Augustina Llanos, Sébastien Dejean, Virginie Neugnot-Roux, Jean Marie

François, Jean-Luc Parrou

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Augustina Llanos, Sébastien Dejean, Virginie Neugnot-Roux, Jean Marie François, Jean-Luc Parrou. Carbon sources and XlnR-dependent transcriptional landscape of CAZymes in the industrial fungus Talaromyces versatilis: when exception seems to be the rule. Microbial Cell Factories, 2019, 18 (1), pp.1-25. 10.1186/s12934-019-1062-8. [https://hal-02177357](https://hal-02177357)

HAL Id: hal-02177357

https://hal.science/hal-02177357

Submitted on 8 Jul 2019

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH

Open Access

<figure><img src="image_3.png" /></figure>

Check for
updates

# Carbon sources and XlnR-dependent transcriptional landscape of CAZymes in the industrial fungus *Talaromyces versatilis*: when exception seems to be the rule

Agustina Llanos<fnref n="1" />,², Sébastien Déjean³, Virginie Neugnot-Roux², Jean M. François<fnref n="1" /> and Jean-Luc Parrou<fnref n="1" />*

**Abstract**

**Background:** Research on filamentous fungi emphasized the remarkable redundancy in genes encoding hydrolytic enzymes, the similarities but also the large differences in their expression, especially through the role of the XlnR/XYR1 transcriptional activator. The purpose of this study was to evaluate the specificities of the industrial fungus *Talaromyces versatilis*, getting clues into the role of XlnR and the importance of glucose repression at the transcriptional level, to provide further levers for cocktail production.

**Results:** By studying a set of 62 redundant genes representative of several categories of enzymes, our results underlined the huge plasticity of transcriptional responses when changing nutritional status. As a general trend, the more heterogeneous the substrate, the more efficient to trigger activation. Genetic modifications of *xlnR* led to significant reorganisation of transcriptional patterns. Just a minimal set of genes actually fitted in a simplistic model of regulation by a transcriptional activator, and this under specific substrates. On the contrary, the diversity of *xlnR*⁺ versus *ΔxlnR* responses illustrated the existence of complex and unpredicted patterns of co-regulated genes that were highly dependent on the culture condition, even between genes that encode members of a functional category of enzymes. They notably revealed a dual, substrate-dependant repressor-activator role of XlnR, with counter-intuitive transcripts regulations that targeted specific genes. About glucose, it appeared as a formal repressive sugar as we observed a massive repression of most genes upon glucose addition to the mycelium grown on wheat straw. However, we also noticed a positive role of this sugar on the basal expression of a few genes, (notably those encoding cellulases), showing again the strong dependence of these regulatory mechanisms upon promoter and nutritional contexts.

**Conclusions:** The diversity of transcriptional patterns appeared to be the rule, while common and stable behaviour, both within gene families and with fungal literature, the exception. The setup of a new biotechnological process to reach optimized, if not customized expression patterns of enzymes, hence appeared tricky just relying on published data that can lead, in the best scenario, to approximate trends. We instead encourage preliminary experimental assays, carried out in the context of interest to reassess gene responses, as a mandatory step before thinking in (genetic) strategies for the improvement of enzyme production in fungi.

**Keywords:** Filamentous fungi, Secretome, Glycosyl hydrolases, Cellulose, Hemicellulose, Biomass deconstruction, Transcription factors, XlnR/XYR1, Glucose repression, High-throughput qPCR

<footnote>*Correspondence: parrou@insa-toulouse.fr</footnote>

<footnote n="1">¹ LISBP, Université de Toulouse, INSA, INRA, CNRS, Toulouse, France
Full list of author information is available at the end of the article</footnote>
## Background

Biomass degradation is an essential yet complex process in nature, particularly for degrading the plant cell wall [1]. This structure is a highly heterogeneous material mainly composed of cellulose, hemicellulose, lignin and pectin, whose degradation is achieved through multiple means, including the use of enzymes with different and complementary activities [2–4]. The secretion of such a wide variety of enzymes is the hallmark of filamentous fungi [1, 5–8]. Part of these enzymes are designed ‘carbohydrate-active enzymes’ (CAZymes) and are compiled into the Carbohydrate Active Enzymes database (http://www.cazy.org/; [9]). They are assigned into five enzyme classes, including Glycoside Hydrolases (GH) and Auxiliary Activities (AA). This classification relies on amino acids sequence similarity, secondary and tertiary fold conservation, and stereochemical architecture of catalytic mechanisms. This classification, however, does not preclude functional redundancy between families, as for instance enzymes bearing xylanase activity that belong to GH10, GH11, GH5, GH7, GH8 and GH43 families [10].

The ever-growing number of sequenced and annotated fungal genomes [11–24] revealed that fungi exhibit tremendous diversity in the number and variety of CAZymes [25–28]. Attempts to correlate growth profiles and gene content in fungal genomes, as for the pectinolytic system [4], indicated that the number of genes could correlate with the ability to degrade specific substrates. Mechanistically however, this correlation appeared drastically more complex when trying to decipher relations between pectin structural elements and the specific subsets of genes used for that deconstruction. From phylogeny viewpoint, the expansion or contraction of gene families is apparently largely lineage specific, and not shared among all fungi of a given lifestyle [29]. In fungi with highly enriched number of CAZymes, it was shown that the increase of GHs is not random across all GH families, but affects specific families to different degrees [18]. Unexpected traits could also be observed, as for plant pathogenic fungi that have in general the highest number of CAZymes while saprophytic fungi contain fewer, even if their high activity in degrading plant biomass would not predict such gene repertoire [28, 30, 31]. Also, the efficient cellulose-degrading fungi as the brown-rot fungi *Postia placenta* faced a significant loss of the gene repertoire typical of cellulolytic microbes, yet developed alternative oxidative mechanisms for biomass deconstruction [32]. Finally, as only a few genes from expanded multi-gene families were actively expressed in a given culture condition, and as it was different in different strains, it has been suggested that family expansion could increase adaptive opportunities rather than activity

under a specific culture condition [33]. This extreme redundancy in gene family members and the diversity of their expression patterns, hence provide foundations for the extraordinary metabolic flexibility reported for filamentous fungi [34].

The filamentous fungus *Talaromyces versatilis*, (formerly *Penicillium funiculosum*), is exploited at the industrial scale to produce an enzymatic cocktail termed *Rovabio®*, which is used as feed additive to enhance the digestibility of cereals-based diets [35, 36] as it is particularly enriched in lignocellulolytic enzymes (i.e. xylanases, arabinofuranosidase, β-xylosidase, etc.) as well as proteases [37]. Independent studies on a *P. funiculosum* strain, also described the richness and hydrolytic potential of the secretome of this fungus when it is grown under different polymeric inducers [38, 39]. However, enzymes recovered in the *Rovabio®* cocktail produced by *T. versatilis* represent only a fraction of the full set of CAZymes that could be in theory encoded by the genome of this fungus [37]. In total, 170 genes were manually annotated as GH or related enzymes, spread into 32 different GH families and relevant for the degradation of the plant polymers (unpublished, proprietary data). This is consistent with other fungal genome versus secretome comparisons, which showed that these microorganisms specifically express and secrete only a subset of CAZymes, generally tailored to adapt to environmental resources [28, 33, 40–43]. This set of secreted proteins is most often species specific, as despite similar genomic potential of certain fungi, their approaches to degrade plant biomass indeed differ markedly in the overall patterns of regulation, activities as well as the specific sets of proteins they employ [32, 41, 44]. A recent, comparative transcriptomic analysis of *A. niger* and *T. reesei* grown on sugarcane bagasse confirmed that these fungi exhibited differences in the precise set of expressed genes and their expression profile, suggesting that they employ different strategies for biomass breakdown, even if a set of CAZymes such as cellulases, hemicellulases and oxidative enzymes was expressed by both fungal strains [45]. It was nevertheless shown that predicted secretomes of fungi with similar lifestyles shared certain characteristics, and not always correlated to gene repertoire and phylogenetic groups [33]. Expression of the secretome was for example largely influenced by the type of substrate or cultivation process [41, 46–48], growth development and pathogenic phase [30]. Also, key genomic mutations could significantly improve enzyme production [49]. All these data point to strict fine tuning of these sets of genes through sophisticated regulatory systems, to achieve optimal growth or adaptation to a given environmental condition [40], as for example in the model of sequential expression
of the genes for an efficient degradation of the lignocellulose [50].

A handful of transcription factors (TF), characterized as activators or repressors of the expression of genes encoding CAZymes, have been isolated in different filamentous fungi (reviewed in [5, 34, 51–53]). Notably, CreA/CRE1 has been documented as the carbon catabolite transcriptional repressor of cellulolytic and xylanolytic genes in response to glucose. However, in filamentous fungi as *A. nidulans*, *N. crassa* and *T. reesei*, if CreA/Cre1/CRE1-dependent derepression occurs in the presence of complex carbohydrates such as lignocellulose, repression yet persists in response to high concentrations of alternative, simple carbon sources such as cellobiose and xylose, which results in CreA/Cre1 nuclear localisation in *A. nidulans* [54]. The transcription factor XlnR/XYR1 was originally reported as a xylanolytic regulator [55]. However, it was emphasized over time that large differences can be found in the sets of genes that are controlled by XlnR in different fungi species [56, 57]. In *M. oryzae* for example, it regulates the expression of genes involved in the pentose catabolic pathway, but not genes encoding hemi-cellulolytic enzymes [58].

The aim of this study was to delineate the role of the transcription factor XlnR in *T. versatilis* by the use of two mutant strains, $\Delta xlnR$ and $xlnR^+$ (deletion and overexpression of the *xlnR* gene, respectively). We explored by high-throughput quantitative PCR the transcriptional landscape of 62 representative genes encoding CAZymes, in the wild type and the two *xlnR* mutant strains. *T. versatilis* was cultivated under a broad range of nutritional conditions, from simple sugars to more heterogeneous polymeric substrates, to evaluate further substrate-dependent transcriptional regulation of these genes, and especially the role of glucose.

## Results

### Balanced choice of GOIs to consider redundancy in the genes coding for similar hydrolytic activities in *T. versatilis*

To better understand the transcriptional network that governs gene expression in the filamentous fungus *Talaromyces versatilis*, (formerly known as *Penicillium funiculosum*), we selected a set of 62 representative genes of interest (GOIs) encoding hydrolytic and auxiliary enzymes in plant biomass degradation (Table 1 and more detailed information in the Additional file 1). A preliminary list of GOIs was selected based on expressed candidates in a RNA-seq study carried out during cultures on glucose and wheat straw [59], and from a proteomic analysis of the Rovabio® cocktail [37]. This list of expressed genes was extended to be more representative of several categories of enzymes. The category ‘cellulases’ gathered

together 3 types of enzymes involved in the degradation of cellulose, namely endoglucanases, cellobiohydrolases and β-glucosidases. The categories ‘backbone hydrolases’, (mainly xylanases), ‘β-xylosidases’, ‘arabinofuranosidases’ and ‘other side chain hydrolases’, (mainly esterases), contained enzymes for hemicellulose degradation. Pectin is another important component of the plant cell wall and its degradation involves, among others, polygalacturonases and pectin lyases, which were included in the study under the category ‘pectinases’. Finally, we considered proteins that have no characterized enzymatic activity but are potentially involved in the degradation process, such as swollenins and hydrophobic binding proteins grouped under the category ‘auxiliary activities’. Figure 1 shows the affiliation of each gene to CAZY GH families. The ‘cellulases’ category notably illustrates the functional redundancy that can be found between GH families, with the six selected endoglucanases distributed in GH5, GH7, GH45 and GH74 families, the GH7 family also containing two of the three cellobiohydrolases.

The 700-bp upstream regions of the coding sequences of these GOIs were scanned for the canonical XlnR binding motif -5'-GGCTAA-3'- [55, 60, 61], as well as for the functional variant, low affinity XlnR binding motif 5'-GGCWAW-3' [60]. We also searched for the CreA binding motif -5'-SCGGRG-3'- [62, 63]. It is worth noticing that almost all selected GOIs exhibit such regulatory sequences in their promoter: 85% possess at least one and up to five CreA motif, 80% present one to four XlnR motif, exactly half of them harbouring the canonical -5'-GGCTAA-3' motif (Additional file 1). From the exhaustive list of the 170 putative genes encoding GH found in the annotated genome of *T. versatilis* (unpublished data), we found that these XlnR and CreA binding motifs present an interesting position bias in the promoter, half of the canonical XlnR binding sites (i.e. interquartile range) being located between -190 and -450 bp upstream the ATG, and between -230 and -480 bp for the CreA motif. These data are in agreement with those reported in *Fusarium* species for which the transcription binding sites mostly located within 600 bp from gene starts [64]. If their mere presence in the promoter regions does not necessarily result in their functionality with respect to transcriptional regulation (see discussion section), this observation suggests that some of these binding sites might be functional in *T. versatilis*.

### A global transcriptional landscape that highlighted the wide plasticity of gene responses, specific of ‘gene-condition’ pairs

We explored 13 different nutritional conditions from simple sugars to complex lignocellulosic material, i.e.
<table><caption>Table 1 GOIs and their affiliation in the different functional categories</caption><thead><tr><th>GOI #</th><th>CAZy family</th><th>Annotation/activity</th><th>Category</th></tr></thead><tbody><tr><td>GOI1</td><td>GH1</td><td>β-Glucosidase</td><td>Cellulase</td></tr><tr><td>GOI2</td><td>GH6</td><td>CbhII</td><td>Cellulase</td></tr><tr><td>GOI3</td><td>GH7</td><td>CbhI</td><td>Cellulase</td></tr><tr><td>GOI4</td><td>GH3</td><td>β-glucosidase</td><td>Cellulase</td></tr><tr><td>GOI5</td><td>GH5</td><td>β-1,4-Endoglucanase</td><td>Cellulase</td></tr><tr><td>GOI6</td><td>GH5</td><td>β-1,4-Endoglucanase</td><td>Cellulase</td></tr><tr><td>GOI7</td><td>GH5</td><td>β-1,4-Endoglucanase</td><td>Cellulase</td></tr><tr><td>GOI8</td><td>GH7</td><td>Cellobiohydrolase (reducing end)</td><td>Cellulase</td></tr><tr><td>GOI9</td><td>GH7</td><td>β-1,4-Endoglucanase</td><td>Cellulase</td></tr><tr><td>GOI10</td><td>GH45</td><td>β-1,4-Endoglucanase</td><td>Cellulase</td></tr><tr><td>GOI11</td><td>GH74</td><td>Endoglucanase</td><td>Cellulase</td></tr><tr><td>GOI12</td><td>GH62</td><td>abf62a</td><td>Arabinofuranosidase</td></tr><tr><td>GOI13</td><td>GH51</td><td>abfA1</td><td>Arabinofuranosidase</td></tr><tr><td>GOI14</td><td>GH62</td><td>Abf62b</td><td>Arabinofuranosidase</td></tr><tr><td>GOI15</td><td>GH62</td><td>Abf62c</td><td>Arabinofuranosidase</td></tr><tr><td>GOI16</td><td>GH54</td><td>Abf-B2</td><td>Arabinofuranosidase</td></tr><tr><td>GOI17</td><td>GH54</td><td>Abf-B1</td><td>Arabinofuranosidase</td></tr><tr><td>GOI18</td><td>GH43</td><td>β-1,4-Xylosidase/α-arabinofuranosidase</td><td>Arabinofuranosidase</td></tr><tr><td>GOI19</td><td>GH54</td><td>α-L-Arabinofuranosidase</td><td>Arabinofuranosidase</td></tr><tr><td>GOI20</td><td>GH54</td><td>α-L-Arabinofuranosidase</td><td>Arabinofuranosidase</td></tr><tr><td>GOI21</td><td>GH54</td><td>α-L-Arabinofuranosidase</td><td>Arabinofuranosidase</td></tr><tr><td>GOI22</td><td>GH11</td><td>XynB</td><td>Backbone hydrolase</td></tr><tr><td>GOI23</td><td>GH11</td><td>XynC</td><td>Backbone hydrolase</td></tr><tr><td>GOI24</td><td>GH10</td><td>XynD</td><td>Backbone hydrolase</td></tr><tr><td>GOI25</td><td>GH11</td><td>XynE</td><td>Backbone hydrolase</td></tr><tr><td>GOI26</td><td>GH11</td><td>XynF</td><td>Backbone hydrolase</td></tr><tr><td>GOI27</td><td>GH11</td><td>XynG</td><td>Backbone hydrolase</td></tr><tr><td>GOI28</td><td>GH11</td><td>XynH</td><td>Backbone hydrolase</td></tr><tr><td>GOI29</td><td>GH11</td><td>XynI</td><td>Backbone hydrolase</td></tr><tr><td>GOI30</td><td>GH5</td><td>β-1,4-Endomannanase</td><td>Backbone hydrolase</td></tr><tr><td>GOI31</td><td>GH3</td><td>β-Xylosidase</td><td>Beta-xylosidase</td></tr><tr><td>GOI32</td><td>GH43</td><td>β-Xylosidase</td><td>Beta-xylosidase</td></tr><tr><td>GOI33</td><td>GH43</td><td>β-Xylosidase</td><td>Beta-xylosidase</td></tr><tr><td>GOI34</td><td>GH43</td><td>GH43</td><td>Beta-xylosidase</td></tr><tr><td>GOI35</td><td>GH43</td><td>GH43</td><td>Beta-xylosidase</td></tr><tr><td>GOI36</td><td>GH3</td><td>β-Xylosidase</td><td>Beta-xylosidase</td></tr><tr><td>GOI37</td><td>GH43</td><td>β-Xylosidase</td><td>Beta-xylosidase</td></tr><tr><td>GOI38</td><td>GH61</td><td>Polysaccharide monooxygenase</td><td>Auxiliary activity</td></tr><tr><td>GOI39</td><td>-</td><td>Hydrophobic binding protein</td><td>Auxiliary activity</td></tr><tr><td>GOI40</td><td>-</td><td>Swollenin</td><td>Auxiliary activity</td></tr><tr><td>GOI41</td><td>-</td><td>Hydrophobic binding protein</td><td>Auxiliary activity</td></tr><tr><td>GOI42</td><td>-</td><td>Hydrophobic binding protein</td><td>Auxiliary activity</td></tr><tr><td>GOI43</td><td>-</td><td>Swollenin</td><td>Auxiliary activity</td></tr><tr><td>GOI44</td><td>GH28</td><td>Polygalacturonase</td><td>Pectinase</td></tr><tr><td>GOI45</td><td>GH53</td><td>Galactanase</td><td>Pectinase</td></tr><tr><td>GOI46</td><td>-</td><td>Pectin lyase</td><td>Pectinase</td></tr><tr><td>GOI47</td><td>GH28</td><td>Endorhamnogalacturonase</td><td>Pectinase</td></tr><tr><td>GOI48</td><td>GH28</td><td>Endopolygalacturonase</td><td>Pectinase</td></tr><tr><td>GOI49</td><td>GH28</td><td>Exopolygalacturonase</td><td>Pectinase</td></tr></tbody></table>
<table><caption>**Table 1 (continued)**</caption><thead><tr><th>GOI #</th><th>CAZy family</th><th>Annotation/activity</th><th>Category</th></tr></thead><tbody><tr><td>GOI50</td><td>GH28</td><td>Exopolygalacturonase</td><td>Pectinase</td></tr><tr><td>GOI51</td><td>GH28</td><td>Xylogalacturonan hydrolase</td><td>Pectinase</td></tr><tr><td>GOI52</td><td>GH78</td><td>α-L-rhamnosidase</td><td>Pectinase</td></tr><tr><td>GOI53</td><td>GH78</td><td>α-L-rhamnosidase</td><td>Pectinase</td></tr><tr><td>GOI54</td><td>GH27</td><td>α-Galactosidase</td><td>Other side chain hydrolase</td></tr><tr><td>GOI55</td><td>GH67</td><td>α-Glucuronidase</td><td>Other side chain hydrolase</td></tr><tr><td>GOI56</td><td>GH67</td><td>α-Glucuronidase</td><td>Other side chain hydrolase</td></tr><tr><td>GOI57</td><td>–</td><td>Acetylxylan esterase</td><td>Other side chain hydrolase</td></tr><tr><td>GOI58</td><td>–</td><td>Acetylxylan esterase</td><td>Other side chain hydrolase</td></tr><tr><td>GOI59</td><td>–</td><td>Carboxylesterase</td><td>Other side chain hydrolase</td></tr><tr><td>GOI60</td><td>–</td><td>Carboxylesterase</td><td>Other side chain hydrolase</td></tr><tr><td>GOI61</td><td>–</td><td>Pectinesterase</td><td>Other side chain hydrolase</td></tr><tr><td>GOI62</td><td>–</td><td>Acetylxylan esterase</td><td>Other side chain hydrolase</td></tr></tbody></table>

more or less heterogeneous polymeric substrates, including also carbon and nitrogen starvation. Normalized fold-change (FC) values of the 62 GOIs during the shift from glucose to the condition of interest, were calculated using glucose as the calibrator condition. The two-dimension hierarchical clustering presented in Fig. 2 identified five clusters of genes (horizontal axis), and four clusters of nutritional conditions (vertical axis) (see also Additional file 2: Figure S1 for quantitative visualisation of gene responses in the five clusters). Gene cluster 1 collected most of the genes encoding cellulases as well as

few proteins with auxiliary activities. These genes were strongly induced in media containing carbon sources as straw (mainly hemicellulose), or avicel (mainly cellulose), with average $log_2$ (FC) values higher than 10 (see also grey boxes in Additional file 3: Figure S2a). Conversely, they were downregulated by simple C5-sugars including xylose, as well as by carbon and nitrogen starvation, with an average 4-fold drop ($log_2$ (FC)) as compared to glucose. The downregulation observed in xylose hence contrasted with their activation by xylan, which indicates that residues released by xylan can trigger induction

<figure><img src="image_4.png" /></figure>
signals that prevail over the downregulation due to the
sole presence of xylose. A quite different situation was
observed for cellobiose, a disaccharide made of glucose
and arising from cellulose degradation, which led to a
slight activation of most of the genes from this gene clus-
ter 1 (average log₂ (FC) of *approx.* 2 as compared to glu-
cose, see Additional file 3: Figure S2a).

As for gene cluster 1, expression of GOIs from gene clusters 4 and 5 was remarkably upregulated in culture media containing polysaccharides (see also Additional file 3: Figure S2d, e), but they did not exhibit any downregulation in the other conditions. A specificity of cluster 5 genes was their better activation in the presence of xylan and xylose, particularly for the two xylanase encoding genes, *xynD* (GOI 24) and *xynF* (GOI 26). A good activation by xylan and xylose was also the most clear hallmark of gene cluster 3 (4<log₂(FC)<10, see Additional file 3: Figure S2c). The induction on xylose was however less potent than that observed in response to xylan, specifically for pectinases (see Additional file 4: Figure S3g), which indicates that xylose may act synergistically with other side-products of xylan hydrolysis for the induction of this functional category of genes, as previously indicated for cellulases. The remaining gene cluster, (number 2), which collected the largest set of GOIs and encoded enzymes with a wide spectrum of predicted hydrolytic activities, appeared quite insensitive to environmental stimuli when referring to this heatmap. This may come from the comparison with other gene clusters, which showed extreme responses to some culture conditions. However, the average expression of gene cluster 2 on straw, soybean and during fermentation was far from insignificant (log₂ (FC) around 4, see Additional file 3: Figure S2b). Notably, polymeric substrates as well as xylose 0.2% sporadically gave rise to high FC values for a few genes (log₂ (FC) > 5).

When referring to nutritional conditions, this 2D-cluster
analysis allowed us to distinguish clusters of two
distinct types (Fig. 2). A first type included complex lig-
nocellulosic-based substrates that overall caused strong
upregulation of almost all the GOIs (first and second
nutritional clusters, on top of the heat map). The sec-
ond type included simpler carbon sources, as well as
carbon and nitrogen starvation, with a transcriptional
response mainly characterized by the downregulation of
cluster 1 genes. The singular positioning of the cellobi-
ose condition, which clustered in this group while it did
not lead to downregulation of cellulases, likely resulted
from its overall modest activation of genes from clusters
2 to 5, similar to the activation observed on arabinose
and xylose. Xylose also distinguished from these condi-
tions, being able to trigger the huge activation of a sub-
set of genes, mainly from gene cluster 3. Interestingly,

xylose and arabinose showed a concentration-dependent
response, as indicated by lower FC values of few GOIs
as the concentration of these sugars increased from 0.2
to 1%, an effect that was particularly prevalent for genes
encoding arabinofuranosidases and side chain hydolases
(see Additional file 4: Figure S3c, d).

About the specific effect of carbon starvation, the left boxplot of Fig. 3 confirmed that the expression level of the overall set of GOIs did not significantly change upon removal of this carbon source from the medium, with the exception of the downregulation observed for cluster 1 genes. Conversely, this nutritional starvation moderately induced just a few genes distributed among gene clusters 2 to 5, mainly those encoding pectinases (see Additional file 4: Figure S3g). This phenomenon was therefore too marginal to conclude on a massive release of glucose repression that might have led to more generalized upregulation, as for example in *A. niger* [50]. Glucose repression was nevertheless an active mechanism in *T. versatilis*, as indicated by the strong repression that was observed upon glucose addition to *T. versatilis* mycelium grown on straw for 24 h (right boxplot in Fig. 3; data from our previous RNA-seq study [59]). In this later experiment, we indeed showed a massive extinction of total RPKM corresponding to GH, as it was previously shown in *A. niger* [65] and *T. reesei* [66] by using the same experimental strategy. To conclude, this modest activation of most of the genes in response to simple sugars, (i.e. arabinose or cellobiose, and even xylose, in the absence of glucose), suggested that best induction of most of these GOIs on hemicellulose- and cellulose-based substrates rely on the synergistic action of inducers that come from their hydrolysis.

The loss-of-function of *xlnR* revealed a dual
activator-repressor role of this transcription factor

We previously identified in *T. versatilis* genome a putative homologue to *A. niger* *xlnR*, (62% identity at the protein level), and successfully deleted this gene [67]. The loss-of-function of *xlnR* in *T. versatilis* slightly reduced growth rate of the mycelium on PDA plates and in liquid MM (data not shown), but it did not abolish growth on the different carbon sources of the study, including xylan and xylose, contrary to what was observed in other fungi as *N. crassa*, *A. oryzae* and *Fusarium* species [61, 68–70]. In Fig. 4, FC values illustrate the ratio of normalised transcript levels between the mutant and the wild-type strain, in a condition of interest, and do not consider the shift from glucose to this medium. The glucose condition was therefore a condition of interest as the others, not a calibrator sample as in Fig. 2. Also, to avoid any ambiguity when referring to clustering results in the remaining part of the manuscript, even for description of classification
<figure><img src="image_5.png" /><figcaption>**Fig. 2** Transcriptional landscape in the WT during the shift from glucose to a condition of interest. Double hierarchical clustering of normalized fold change (FC) values, calculated using glucose as the calibrator sample, and expressed as log₂ (FC). Clustered GOIs (top dendogram) and culture conditions (left dendogram), with the corresponding names at the bottom and right of the heatmap, respectively. The white lines in the heat map delineate the main clusters, i.e. 5 and 4 clusters for GOIs and culture conditions, respectively. The coloured bar at the top of the heat map is referring to functional categories affiliation of GOIs (legend at the bottom left of the figure). Red and green colours indicate increased and decreased FC values, respectively (scale at the bottom of the figure, with saturating intensity for absolute values of log₂ (FC) beyond 9). Cq values that exceeded 35 cycles were considered beyond the limit of detection and represented by the grey cases in the heat map. The bottom, light blue heat map depicts the number of each of the 3 binding motifs identified in the 700 bp upstream the ATG start codon of each gene, i.e. the XlnR binding motif -5'-GGCTAA-3', its low affinity variant 5'-GGCWAW-3', and the CreA binding motif -5'-SCGGRG-3'</figcaption></figure>

results presented in the upcoming Figs. 4 and 5, we will
solely give the name 'gene clusters 1-5' to the clusters that
have been presented in Fig. 2 for the WT response.

The clearest hallmark upon deletion of *xlnR* was the
reduced transcript levels for a subset of genes, which
clustered in the left of the heatmap (Fig. 4). This subset
was enriched for genes encoding xylanases, β-xylosidases
and side chain hydrolases, in response to xylose, xylan

and straw (bottom-left corner, Fig. 4). Very interestingly,
while tracking the origin of these genes relative to their
initial classification in the WT strain exposed to these
carbon sources, almost all downregulated genes in the
Δ*xlnR* background belonged to gene clusters 3 and 5 (See
Additional file 3: Figure S2, and illustrative scatter plots
in Additional file 5: Figure S4 for xylose, xylan and straw
conditions). This observation was consistent with their
<figure><img src="image_6.png" /><figcaption>**Fig. 3** Focus on FC values during carbon starvation and after glucose addition. 'Carbon starvation': Shift of the mycelium from glucose medium to medium without carbon source (glucose condition used as calibrator sample; RT-qPCR data from this study); 'Glucose repression': addition of glucose to mycelium grown on wheat straw for 24 h (wheat straw condition as calibrator sample; RNA-seq data from a previous study [59]). The box-plot gathers $log_2$ (FC) values of the 62 GOIs (i.e. those presented in the 'C starvation' line in Fig. 2, and those reminded in the Additional file 1 for the 'glucose repression' condition)</figcaption></figure>

strong activation when shifting the WT mycelium from glucose to xylose or xylan, and suggested a role of this TF in their activation. It has been stated that there is a core set of genes whose expression always depends on XlnR/XYR1, i.e. the genes coding for XlnD, XynB, XynA, BglA and AguA [57]. We identified the putative homologues of some of these genes in *T. versatilis*, as GOI 36 for XlnD, GOI 22 for XynB, GOI 55 and GOI 56 for AguA. These genes also showed in *T. versatilis* a significant fall of their transcript levels in the absence of *xlnR*, both on xylose and xylan (5.5-fold drop ($log_2$ (FC))), sub-cluster on the right of the bottom-left cluster in Fig. 4).

A noticeable observation for the xylan condition, not for xylose, was a counter-intuitive increase in the transcript levels of many GOIs (Fig. 4). These upregulated genes in the mutant strain mainly corresponded to gene clusters 1 and 4 in the WT strain (Additional file 3: Figure S2a, d), and more precisely to genes encoding cellulases and most arabinofuranosidases (Additional file 4: Figure S3c, f). These data indicated that XlnR might have, on xylan specifically, a dual activator-repressor role. This was consistent with the results observed in the WT strain exposed to this carbon source, as we observed in Fig. 2 a much better xylan-dependent activation of GOIs from gene clusters 3 and 5, positively regulated, than GOIs

from gene clusters 1 and 4, negatively regulated by XlnR. Interestingly, these negatively regulated GOIs on xylan were conversely positively regulated during growth on xylose, as well as on other polymeric carbon sources as straw. We can therefore propose that in the wild type strain grown on xylan, the repressor function of XlnR on the transcription of most cellulases and arabinofuranosidases is triggered by specific side-products of xylan hydrolysis, others than xylose, and not present in other polymeric substrates.

Another remarkable feature of $\Delta xlnR$ cells was the general downregulation observed on straw, more pronounced than in response to xylose, both for the amplitude of the drop of transcript levels, and for the number of GOIs affected by the deletion (Additional file 5: Figure S4c). A similar, generalized effect was observed to some extent on other carbon sources such as avicel and soybean. However, for $\Delta xlnR$ cells exposed to soybean, a notable hallmark was the upregulation of GOIs (right gene cluster in Fig. 4). This cluster gathered all pectinases except GOI45, yet similarly upregulated on soybean. This specific, soybean-triggered upregulation of pectinases in the $\Delta xlnR$ strain, contrasted with their significant downexpression in the presence of other carbon sources as straw and avicel (see also Additional file 4: Figure S3g). As previously mentioned for xylan, this feature highlighted another example of dual activator-repressor role of this transcription factor, whose gene-targets are strictly dependent on the nature of the carbon source.

The sole presence of glucose also appeared as a favorable nutritional context to observe a $\Delta xlnR$-dependent, generalized down-regulation of GOIs, with at least 3-fold drop of transcript levels for 25 GOIs. These downregulated genes on glucose encoded, amongst others, 6 out of the 11 cellulases, 5 out of the 10 pectinases, and 4 out of the 6 proteins with auxiliary activities. On other simple carbon sources, the deletion of *xlnR* did not affect the expression of genes significantly, with the notable exception of genes encoding side chain hydrolases on cellobiose (see Additional file 4: Figure S3d).

## The multicopy integration of *xlnR* surprisingly boosted transcript levels of very limited subsets of genes

In the wild type strain, we could notice moderate two-to five-fold activation of *xlnR* itself during the shift from glucose to the different nutritional conditions (data not shown), which may participate in transcription responses depicted in Fig. 2. It has been reported that addition of extra copies of the xylanolytic regulator or its overexpression can increase transcription of glycosyl hydrolases encoding genes, such as xylanases and endoglucanases, with parallel increase of their secretion [69, 71–75]. We therefore investigated the effect of *xlnR* overexpression in
<figure><img src="image_7.png" /></figure>

T. versatilis by using a strain that carries 9 copies of *xlnR*
under the control of its own promoter (*xlnR*⁺ strain), as
determined by a digital droplet PCR assay on genomic
DNA. As in the previous section, the FC values of GOIs
were calculated as the ratio of normalised transcript lev-
els between the *xlnR*⁺ and the wild-type strain, in each
condition. About *xlnR* transcript levels themselves (data
not shown), an *approx.* 8-fold increase of transcript lev-
els could be observed in the *xlnR*⁺ strain incubated on
xylose and xylan, which correlated to the copy number of
*xlnR* in this strain. Its overexpression was still significant
for the remaining polymeric substrates as for cellobiose
(FC around 6), and less pronounced for arabinose, star-
vation conditions and glucose (2 < FC < 4), confirming the
strong dependency of transcript levels upon nature of the
carbon source.

Not surprisingly, data presented in Fig. 5 showed that
the consequence of *xlnR* overexpression was also highly
dependent upon nutritional environment. Globally, it

affected the expression of *approx.* one-third of GOIs
(upper-right cluster), clearly enriched for genes encod-
ing hemicellulases. We observed this positive transcrip-
tional effect under xylose and complex substrates, not on
arabinose and starvation, neither on cellobiose while this
latter carbon source led to a clear accumulation of *xlnR*
transcripts in this *xlnR*⁺ strain. In more details, the best
substrates to promote the effect of *xlnR* overexpression
were xylan and xylose (26 and 20 GOIs with FC > 3, (i.e.
log₂ (FC) > 1.58), respectively, see also Additional file 6:
Figure S5a, b), and to a lesser extent, avicel (15 GOIs),
fermentation (14 GOIs, see the discussion section) and
soybean (10 GOIs). Much to our surprise with respect to
*xlnR* transcript levels themselves, straw was the least effi-
cient condition of this group, as just a few GOIs showed
a positive response to *xlnR* overexpression (only 7 GOIs
with FC > 3, see also Additional file 6: Figure S5c). As
is shown in Additional files 3: Figure S2 and 4: Figure
S3 to better track these genes, xylan and xylose mainly
<figure><img src="image_8.png" /></figure>

triggered overexpression of GOIs originating from gene clusters 4 and 5, and particularly genes encoding hemi-cellulases (i.e. beta-xylosidases, backbone hydrolases, and to a lesser extent, side chain hydrolases and arabinofuranosidases). On the contrary, a relevant hallmark was the insensitivity to xlnR overexpression of GOIs originating from gene cluster 1, i.e. mainly cellulases-encoding genes, whatever the condition (see Additional file 4: Figure S3f).

As it was previously found with the unexpected Δ*xlnR*− triggered upregulation of cellulases on xylan, and of pectinases on soybean, *xlnR*⁺ strain on xylose highlighted a similar dual repressor-activator function of this TF. It was illustrated with the unforeseen downregulation of a dozen of genes showing FC values lower than minus 3, particularly half of the genes encoding arabinofuranosidases, (6 out of the 10 from this category), as well as side chain hydrolases and proteins with auxiliary activities

(Additional file 4: Figure S3). The most unsettling result came from the glucose condition, which also led to a massive downregulation of GOIs in the *xlnR*⁺ strain, (22 genes with FC ≤ 3, see below and Fig. 6e), particularly those encoding cellulases and proteins with auxiliary activities. A similar trend was observed for this specific set of GOIs on cellobiose, with their downregulation. These results emphasized that the possibility for XlnR TF to boost -or repress- target GOIs, is highly dependent upon the nature of the carbon source and on mechanisms other than the regulation of the sole XlnR protein content. This was illustrated for example by avicel and cellobiose conditions, leading to similar *xlnR* changes both in the WT and in the *xlnR*⁺ strain, yet showing highly divergent expression patterns of GOIs, with upregulation of hemicellulases for the former and downregulation of cellulases and auxiliary activities for the latter.
### *xlnR*<sup>+</sup> versus *ΔxlnR* responses: no generalizable shape, only substrate-specific patterns reflecting the complexity of XlnR-dependent transcriptional regulation

When considering the genes most affected by the overexpression of *xlnR* (i.e. upregulated genes of the top-right cluster in Fig. 5) or by its deletion (i.e. down-regulated genes of the bottom-left cluster in Fig. 4), we found that among the 30 collected genes, only 9 were affected both by *xlnR* deletion and by its over-expression. More focused, condition-specific analysis nuanced this conclusion, as illustrated in Fig. 6 that plots log<sub>2</sub> (FC) ratios in the *xlnR*<sup>+</sup> strain as a function of the *ΔxlnR* strain, where genes expected to be downreg- ulated in the *ΔxlnR* strain and upregulated in the *xlnR*<sup>+</sup> strain, will locate in the upper-left sector of these scat- ter plots. The most pertinent pattern was obtained for strains grown on avicel (Fig. 6a). On this carbon source, downregulation observed in the *ΔxlnR* strain was not strong (−3 < log<sub>2</sub> (FC) ≤ 1 for most of GOIs), yet mas- sive. Reciprocally, avicel was also the most efficient car- bon source to trigger quite generalized upregulation of gene expression in the *xlnR*<sup>+</sup> strain, leading to position- ing of most of GOIs in the upper-left sector of the scat- ter plot.

Consistently with results pointing to XlnR being a dual repressor-activator TF, more composite distribu- tions were observed in xylose and xylan scatter plots. On xylose (Fig. 6b), among the 21 hemicellulases that were downregulated in the *ΔxlnR* strain (FC ≤ 3), only 9 were also upregulated in the *xlnR*<sup>+</sup> strain (FC > 3), with an enrichment of backbone hydrolases originat- ing from gene cluster 5. The remaining 12 genes were either weakly upregulated and not considered, or, in agreement with the repressor role of XlnR, downregu- lated in the *xlnR*<sup>+</sup> strain as were most of arabinofura- nosidases. On xylan (Fig. 6c), as on xylose, few genes encoding hemicellulases originating from gene cluster 5 and few others from cluster 3 showed the expected expression pattern, i.e. downregulation in the *ΔxlnR* strain and upregulation in the *xlnR*<sup>+</sup> strain. But this behaviour was marginal, as many upregulated genes in the *xlnR*<sup>+</sup> strain were also upregulated in the *ΔxlnR* mutant, leading to a massive positioning of the genes in the upper-right sector of the graph, particularly for those encoding arabinofuranosidases and originat- ing from gene cluster 4. Another remarkable situation was observed on straw (Fig. 6d), with very few GOIs responding positively to *xlnR* overexpression while this substrate was the most relevant to observe a massive, downregulation of the full set of GOIs in the *ΔxlnR* mutant. Taking into account the data presented in Additional file 4: Figure S3g, we can speculate that the higher the expression in the WT strain, as on straw, the

lower the possibility to further promote transcript lev- els under *xlnR* overexpression.

Finally, the most outstanding pattern was observed when analyzing the effect of XlnR transcription factor on glucose (Fig. 6e), with almost all the genes being down- regulated in both strain. While tracking the identity of the genes, this behavior particularly affected categories as cellulases, (6 out of the 11 GOIs with FC ≤ 3, four oth- ers being nevertheless downregulated in both strain), pectinases, (5 out of the 10 GOIs with FC ≤ 3, two others being also downregulated in both mutant), and enzymes with auxiliary activities (all the 6 members, 4 significantly affected with FC ≤ 3).

### To belong to the same functional category did not necessarily translate into similar regulatory patterns

Data that have been presented in Fig. 2 first illustrated that members of a functional category were generally widespread within the different clusters of the heatmap, revealing the heterogeneity of their response to the dif- ferent nutritional conditions. Genes encoding cellulases appeared as an exception as almost all members gath- ered in gene cluster 1. Going into details however, three additional genes encoding cellulases showed some divergence in their overall response. GOI 8 (cluster 2) was not induced by the polymeric carbon sources, while GOIs 7 and 1 (cluster 2 and 4, respectively) were unaf- fected by the overall downregulation of cluster 1 genes. Similarities and differences between cellulases were best illustrated by pairwise comparisons presented in Fig. 7. We could observe almost perfectly correlated genes (e.g. GOIs 2, 3 and 11; Pearson correlation coef- ficient of 0.99), and more divergent ones (e.g. GOIs 7, 8 and 1; 0.58 < r < 0.88). Second, the relative quantity of transcript levels appeared as another type of heteroge- neity, even in the WT strain (Fig. 8a, x-axis of the left scatter plot). Similar FC values during the shift from glucose to a condition of interest (y-axis) could indeed hide significant differences in the strength of their pro- motor (3- to 4-log). Third, the scattering of the genes in response to the different *xlnR* backgrounds, previously discerned in Figs. 4 and 5 heatmaps, was perfectly illus- trated when plotting log<sub>2</sub> (FC) ratios in the *xlnR*<sup>+</sup> strain as a function of the *ΔxlnR* strain (Fig. 8a right scatter plot). Whether for xylan, xylose or soybean, both the type of regulation (up- or downregulation) and the amplitude of FC values varied significantly from one gene to another. For substrates as straw or avicel, a better clustering of most of the members did not pre- vent the presence of outliers, away from the median response of this functional category.

The heterogeneity of responses in the regulation of the arabinofuranosidase category was another interesting
<figure><img src="image_9.png" /><figcaption>**Fig. 6** Specificities of *xlnR*⁺ and *ΔxlnR* responses under different culture conditions. Scatter plots of log₂(FC) values in the *xlnR*⁺ strain as a function of log₂(FC) values in the *ΔxlnR* strain. One scatter plot per culture condition. The GOIs are gathered according to the five gene clusters defined in Fig. 1</figcaption></figure>

example. As compared to cellulases, we could notice wider differences between members, both for FC values and relative quantities of transcript levels in a condition of interest (Fig. 8b, left), and for the singularity of patterns in the different *xlnR* backgrounds that clearly illustrated the dual activator-repressor function of XlnR (Fig. 8b, right). Finally, we examined coregulations between genes encoding backbone hydrolases when changing *xlnR*

genetic background, as we previously underlined that a few members of this functional category fitted with logic down- and upregulation in the *ΔxlnR* and *xlnR*⁺ strains, respectively, especially on xylose. The correlation matrices drawn for the WT, *ΔxlnR* and *xlnR*⁺ strains (Additional file 7: Figure S6), showed that even for the few genes that exhibited rather good pairwise correlations in the WT strain (Additional file 7: Figure S6a, arbitrary
<figure><img src="image_10.png" /></figure>

green-coloured squares when $r > 0.90$), these correlations were dismissed in the mutant strains. This was illustrated by the clear dispersion of the points and drop of Pearson's coefficients in the $\Delta xlnR$ and $xlnR^+$ matrices (refer to the same, green-coloured squares in Additional file 7: Figure S6b, c). This result indicated that coinduction of two

genes in the WT strain, under the different substrates, did not necessarily go along with quantitative coregulation of these genes under $\Delta xlnR$ and $xlnR^+$ backgrounds. As a whole, all these observations showed that no consensus response could be proposed within a functional category.
Discussion

A complex transcriptional landscape, mainly shaped
by medium composition

In this work, we have analysed the expression patterns
of 62 representative genes encoding hydrolytic enzymes
in the fungus *Talaromyces versatilis*, (formerly known
as *Penicillium funiculosum*), to highlight key features in

their transcriptomic response to different carbon sources.
As previously reported in different fungi [70, 76–81],
our study underlined the huge plasticity of transcrip-
tional responses when changing nutritional status, with
complex patterns of co-regulated groups of genes. As a
general trend, the more heterogeneous the polymeric
substrate, the more efficient to trigger transcriptional

<figure><img src="image_11.png" /><figcaption>**Fig. 8** Heterogeneity of expression data within functional categories. **a** 11 GOs encoding cellulases and **b** 10 genes encoding arabinofuranosidases, in response to five selected carbon sources (coloured code). Left scatter plot: log₂(FC) values, (i.e. shift from glucose to the medium of interest as in Fig. 2), as a function of the relative quantities of transcript levels (RQs), in the wild type strain. Right scatter plot: log₂(FC) values in the xlnR⁺ strain as a function of log₂(FC) values in the ΔxlnR strain. RQs were calculated using as reference the ‘GOI-condition’ pair sample leading to the highest Cq value in the full set of data from this functional category of genes (i.e. lower transcript level). No symbol drawn at RQ of 1 for the arabinofuranosidases category (**b**), as the culture condition leading to the lowest expression within this group of genes was not plotted (i.e. 5 conditions plotted amongst the 12 used for calculation)</figcaption></figure>
activation when *T. versatilis* was shifted from glucose to these rich media, whose biochemical composition further changes dynamically during its assimilation [50]. Especially, growth on wheat straw, or on the fermentation medium used for industrial cocktail production, led to a powerful activation of nearly all GOIs as previously described for *A. niger* and *T. reesei* also grown on wheat straw [65, 66, 78]. However, despite global similarity of the responses to polymeric substrates, differences in the expression levels of individual genes were observed, as it was found in *A. niger* after comparative RNA-seq based transcriptome responses to wheat straw [65] and coppice willow [79].

More simple sugars as arabinose and cellobiose appeared as weak inducers in *T. versatilis*, strengthening the idea that efficient inducing mechanisms also rely in this fungus on the production of small molecular weight compounds that act in a synergistic manner [34, 50, 65]. Another illustrative result was the comparison of pure xylose versus xylan as substrates: these two carbon sources upregulated almost the same set of genes but the induction was less effective on xylose than on xylan, a result that was already stated by van Peij et al. [72]. In our study, this effect was particularly marked for pectinases, indicating that side-products of xylan hydrolysis may act synergistically with xylose to promote expression of these genes. Similarly in *A. niger*, wheat straw pretreatment abolished the expression of a large set of genes encoding pectinolytic enzymes [82]; these authors explained their observation by the removal of parts of the lignocellulose during pretreatment of this substrate, components which were probably important for induction.

## A limited role of XlnR as direct, positive regulator

The transcription factor XlnR/XYR1 [55], appeared over time as a key regulator for activation of both hemicellulolytic and cellulolytic genes. Expression of most *xlnR/xyr1/xlr1* orthologs is not specifically induced (reviewed in [53]), with a few exceptions as in *T. reesei* where *xyr1* expression appeared strongly activated by sophorose, lactose and cellulose [83, 84]. It was also reported that addition of extra copies of this TF or its overexpression can increase transcription of target genes, with parallel increase of enzyme secretion [69, 71–75]. Again, counterexamples could be found, with the report of a limited impact of this overexpression [85, 86].

As *T. versatilis* genome encodes a protein with strong homology to *A. niger* XlnR/XYR1, we therefore analysed the effects of both the deletion and the overexpression of *xlnR* in this fungus. In our study, we could notice moderate two- to five-fold activation of *xlnR* itself during the shift from glucose to the different nutritional conditions (data not shown). Overexpression upon

multicopy integration of *xlnR* further promoted its transcript levels, but overall, whether for the wild type or the overexpression strain, the relatively modest, carbon source-dependant changes in *xlnR* transcript levels themselves contrasted with the extreme transcript levels of some putative XlnR target genes. These results hence supported the idea that, as in other fungi, XlnR-dependent gene expression does not primarily rely on the sole control of its transcript levels, but on post-translational control [87].

Both mutations, either the deletion or the overexpression of *xlnR*, indeed led to significant reorganisation of GOIs' transcriptional patterns, confirming its role in the global regulation of gene network in *T. versatilis*. However, the full set of “gene-condition” pair possibilities, in the three strain backgrounds, highlighted an outstanding level of complexity of coregulations, especially illustrated by the diversity of $xlnR^+$ versus $\Delta xlnR$ patterns. Just a minimal set of GOIs and/or conditions actually fitted in a simplistic model of regulation by a transcriptional activator, i.e. downregulation in the $\Delta xlnR$ strain and upregulation under XlnR overexpression. The avicel was not the most powerful carbon source for transcriptional activation, but almost all GOIs respected this simple regulatory scheme. Xylose and xylan, classically reported as inducers of hemicellulolytic gene expression, actually triggered the expected down- or upregulation of a tiny group of genes from cluster 5, mainly encoding xylanases and beta-xylosidases. As these few genes showed an enrichment for the XlnR binding motif in their promoter region, (average of 2.5 motif/gene if considering both the canonical and low affinity motif, versus 1.5 for the whole selection of GOIs), these specific regulatory responses in the WT and mutant strains may rely on a positive and probably direct binding of XlnR on the promoter of these genes. Targeted ChIP-PCR or more global ChIPseq-based experiments as in [88] might nevertheless be necessary to clarify the specific role of XlnR and define directly regulated genes by this TF in these conditions, but also in all conditions of interest. Indeed, while variable number of genes were downregulated in the $\Delta xlnR$ mutant strain, the moderate and even the absence of effect of XlnR overexpression was conversely very surprising in conditions other than avicel and xylose. On straw for example, the massive, downregulation of the full set of GOIs in the $\Delta xlnR$ mutant indeed contrasted with the ineffectiveness of its overexpression. Such a lack of activation upon XlnR/XYR1 overexpression was also reported for transcriptional activation of xylanases in *F. oxysporum* [86] and might find explanation in the post-translational control of this TF. Specifically, it was proposed that already high and saturating activity of XlnR in the wild type strain, e.g. through various levels
of phosphorylation elicited by such lignocellulolytic substrates [87], could prevent further activation of gene expression.

Beyond these mechanistic considerations and bearing in mind the zootechnical performance of the 'Rovabio® Advance' cocktail that is produced from the $xlnR^+$ strain, (e.g. body weight gain and feed conversion rate, personal communications from ADISSEO FRANCE SAS and [36]), we specifically focused on this surprising, limited number of overexpressed genes in the fermentation condition (i.e. industrial process). Indeed, solely 14 hits with FC values higher than 3 were identified in this strain, i.e. 5 out of 10 genes encoding pectinases and just a few genes encoding hemicellulases (2 backbone hydrolyses and 2 beta-xylosidases, 3 arabinofuranosidases, 2 side chain hydrolases). Proteomic analysis of the secreted proteins by the wild type and the $xlnR^+$ strain at the end of the fermentation process ([37] and personal communication from O. Guais, CINABio-ADISSEO), showed that despite comparable overall protein concentration and similar number of hits in the secretome of these strains, (67 and 65 characterized proteins, respectively), $xlnR^+$ secretome was enriched in pectinases and hemicellulases. These results actually agree with overexpression data and reinforce the idea that targeting just a minimal number, yet adequate set of genes, can enhance the performance of a cocktail.

## Literature-based predictions of gene regulation: a risky strategy

The influence of the carbon source appeared of prime importance to delineate the set of regulated genes in the strains of interest. In that context, our and previous results from fungal literature showed that it definitely seems difficult to predict the effect of XlnR, even when referring to data from closely related fungi [57, 77, 81]. The plethora of conditions found in the reports, (e.g. nature of the substrate and/or inducer, duration and concentration of the exposure, strain background and overexpression strategy), altogether influence the role of the transcriptional regulators and ultimately, conclusions on gene responses. This makes comparisons between studies extremely tricky, and generalizations often biased, particularly when they are drawn irrespective of the carbon source, or based on a minimal number of GOIs and/or non-overlapping sets of genes between the studies.

A focus on genes encoding cellulases illustrates this diversity of regulatory mechanisms between media and/or fungi (see the review by [89]). For example, a statement as *T. reesei* being the only fungus with a clear XYR1 dependent regulation of the cellulolytic system must be

nuanced, in light of observations made in *A. niger* [72] or *A. oryzae* [69] where their transcription was also heavily impaired or abolished. Other studies concluded for the absence or slight effect of XlnR/XYR1 system on cellulase production in *N. crassa* [70], *F. graminearum* [68] or *T. cellulolyticus* [90]; However, slightly different *T. cellulolyticus* strain backgrounds led to different conclusions between Fujii' and Okuda' studies [75, 90], XlnR being from the latter the main regulator for cellulase gene expression and activity in this fungus. Globally,KITAKE's statement that the regulatory system of cellulolytic enzymes is not only partially conserved, but also significantly diverse [89], can probably be applied to most functional systems. It also strengthens the idea that experimental reassessment of gene responses is mandatory if conclusions must be drawn in other experimental contexts, whether nutritional or genetic. Results obtained in *T. versatilis* did not contravene this idea.

Relative to other fungi, the global behaviour of genes encoding cellulases in *T. versatilis* appeared quite reminiscent of *N. crassa* [70], as they presented almost no effect of XlnR, whether for the deletion or for the overexpression of this TF. This overall behaviour also distinguished these genes from those encoding the other categories of enzymes, suggesting that XlnR should not be involved in the control of their expression in *T. versatilis*. This without counting on their unpredicted downregulation upon $xlnR$ overexpression on cellobiose, and their upregulation upon $xlnR$ deletion on xylan. Such indication for a repressor function of XlnR on cellulases has been reported in literature, yet on glucose, with the upregulation of many genes upon *xyr1* deletion in *T. reesei* [77]. In *T. versatilis* grown on glucose, we observed the reverse situation with almost all the genes encoding cellulases being clearly downregulated in the deletion strain (10 out of the 11 candidates), the most outstanding result being nevertheless, still on glucose, their similar downregulation upon $xlnR$ overexpression. Mechanistically, this striking, downregulation in both mutant strain grown on glucose was inexplicable. But these effects of XlnR were most likely indirect, as a decrease of XlnR motif representation was observed in their promoter as compared to the whole set of GOIs (either the canonical or the low affinity consensus motif), which was actually consistent with the overall lack of effect of XlnR on these genes.

As a conclusion, if we can minimize the role of XlnR on cellulase regulation in our whole set of experimental conditions, the strong induction of these genes in the wild type strain grown in the presence of cellulolytic and hemicellulolytic substrates implies the existence of alternative transcriptional control. Together with their
activation by cellobiose, yet modest, as in *N. crassa* [91, 92], all these results point to the possible contribution of CLR-2/ClrB ortholog, also found in *T. versatilis* genome and proposed as a major regulator of the cellulolytic system in filamentous fungi as *N. crassa*, *A. nidulans* and *A. aculeatus* [91, 93].

**The dual repressor-activator role of TFs: not so marginal observation in *T. versatilis*, and beyond**

Besides genes encoding cellulases, this study revealed other examples of counter-intuitive transcripts regulations that targeted very specific sets of genes, and this, in precise nutritional contexts. Relevant situations were the lower expression of genes encoding auxiliary activities under XlnR overexpression, during growth on xylose and cellobiose; The upregulation of few genes encoding arabinofuranosidases in the Δ*xlnR* strain on xylan, and their downregulation in the *xlnR*⁺ strain on xylose; The upregulation of pectinolytic genes in the Δ*xlnR* strain on soybean, or their downregulation in the *xlnR*⁺ strain on arabinose. Such observations of a dual, substrate-dependant repressor-activator role of XlnR are not new or unusual. We previously cited the upregulation of many cellulase encoding genes upon *xyr1* deletion in *T. reesei* grown on glucose [77]; Bischof et al. highlighted a similar upregulation of gene expression in a Δ*xyr1* mutant of *T. reesei*, for wheat-straw induced chitinases and mannosidases [94]. However, based on variable occurrence and even absence of the Xyr1 binding site in the promoter of these genes, these authors proposed that the action of Xyr1 as a repressor was very unlikely, and they rather considered an indirect, growth-related induction mechanism as this deletion mutant was unable to grow on wheat straw. This alternative seems unlikely in *T. versatilis* since this fungus did not exhibit such growth defects. In addition to XlnR, Kunitake et al. reported highly variable effect of ClbR overexpression on cellulolytic gene expression in *A. aculeatus* [95], with both up- and downregulated genes as compared to WT strain, which is quite reminiscent of the very heterogeneous behaviour of our set of arabinofuranosidases in the *xlnR*⁺ strain grown on xylose. Castro et al. also found that various genes encoding cellulolytic and xylanolytic enzymes were downregulated in the *T. reesei* Δ*cre-1* mutant, even in presence of glucose, suggesting that CRE-1 can act positively on these genes in a direct or indirect manner [77]. Similar observations were reported by Sun et al. with a huge cluster of downregulated genes in the *N. crassa* Δ*cre-1* mutant, demonstrating that CRE-1 functions as a global transcription factor affecting both gene repression and activation in *N. crassa*, both directly and indirectly [96]. Still in this fungus, Thieme et al. recently

reported the de-repression of several pectinolytic genes in the absence of PDR-1 [97]. Several lines of evidences allowed them concluding that besides its activator role, this function in repression of transcription most likely occurs indirectly, in concert with one or several TFs as part of a larger regulatory network. All these examples of dual activator-repressor role of transcription factors, including XlnR, strengthen the extreme specificity of transcriptional responses, with strong influence of both the carbon source and the promoter region context on transcription factors properties.

**Glucose effects in the transcriptional control of GH: more complex patterns than for the expected CCR**

The general mechanisms leading to Carbon Catabolite Repression (CCR) is largely conserved in *S. cerevisiae* and filamentous fungi during growth in the presence of glucose (Review by [54]). In filamentous fungi, this mechanism relies on the pivotal transcriptional regulator CreA/Cre1, whose putative homolog was also identified in *T. versatilis* genome. As it was previously shown in *A. niger* and *T. reesei* [65, 66], glucose did appear in *T. versatilis* as a formal repressive sugar if referring to the massive repression of most GOIs, including cellulases, upon glucose addition to the mycelium grown on wheat straw. While studying the role of PfMig1, (i.e. Cre1/CreA homolog), in the *Penicillium funiculosum* NCIM1228 strain, (also found as *Talaromyces funiculosus*), Randhawa et al. similarly described a PfMig1-dependant, repressive effect of glucose when this sugar was added to the fungus grown on avicel, with major drop of cellulase activity [98]. This later result suggests that in *T. versatilis*, the effect of glucose addition on transcription levels also relies on this transcription factor Cre1/CreA/Mig1. Besides glucose repression, we also noticed a dose-dependent inhibitory effect of xylose and arabinose when used as sole carbon sources. A similar effect was observed in other filamentous fungi, which do not only induce CCR in response to glucose but also to high concentrations of alternative carbon sources, as cellobiose and xylose in *A. niger*, *A. nidulans* and *T. reesei* [83, 99–102]. In such conditions of high substrate concentration, CreA/Cre1 predominantly localized in the nucleus, hence leading to transcriptional repression [83, 99, 102]. Such a mechanism should be transposable to *T. versatilis*, however, this dose-dependent inhibitory effect of xylose and arabinose was observed only for arabinofuranosidases and side chain hydrolases, not for xylanases and cellulases as in other fungi.

The massive drop of transcript levels of most GOIs upon addition of glucose could be considered as a definite proof that CCR is effective in *T. versatilis*. However,
another interesting condition was the removal of glucose from the culture medium (i.e. glucose starvation condition), also found as “de-repressed” or “no-Carbon” condition in literature. The (large-scale) studies carried out in *T. reesei*, *A. niger*, *A. nidulans* and *N. crassa*, reported the release of repression of cellulolytic and xylanolytic genes in the absence of glucose [50, 63, 65, 76, 77, 91, 96, 99, 103], whose products may then act as scouting enzymes playing a foraging role on complex substrates, releasing the inducing molecules that trigger the main hydrolytic response [50, 103, 104]. With the exception of just a few derepressed genes, mainly encoding pectinases, no clear derepression has been observed in *T. versatilis* under glucose starvation condition. On the contrary, we surprisingly observed a drop of transcript levels for cluster 1 genes, (mainly cellulases), which contrasted with previous reports showing *cbh1* expression being unchanged [76, 105] or moderately activated [77] during the shift of *T. reesei* from glucose to a starvation condition. As this phenomenon was also observed with simple, alternative sugars as xylose and arabinose, we favoured the idea that these lower transcript levels of genes encoding cellulases actually result from the absence of glucose after the shift to these different media. That is indeed a common denominator in these conditions, with the exception of the N-starvation context that contained glucose but showed moderate downexpression of almost all GOIs in the WT strain. To our knowledge, this is the first report of a negative effect of the absence of glucose on such functional categories of genes, (here those encoding cellulases and most auxiliary activities), highlighting a counter-intuitive, positive role of this sugar on the basal expression of these specific genes. As a reminder, we observed the opposite trend for genes encoding pectinases, demonstrating once again the strong dependence of these regulatory mechanisms on the promoter context.

The atypical regulatory patterns observed for several “gene-condition” pairs hence challenged the presence of CreA binding motifs in the promoter region of these genes, (average of 1.7 CreA binding motif per promoter region for the whole selection of GOIs). Especially for cellulases upon glucose starvation, the opposite regulatory behaviour might have been expected according to fungal literature on the cellulolytic system (see above), and because they contain an average of 2.2 consensus per gene promoter, which statistically favours a situation where genes that contain multiple adjacent motifs in their promoter regions are more likely direct targets of CRE-1 [96]. Also, while genes that belong to clusters 4 and 5 contain, respectively less than 1 and slightly more than two motifs per gene, no clear difference could be noticed in their overall response to this carbon starvation

condition. Our results then confirmed that it is difficult to anticipate regulatory patterns just referring to TFs binding-sites identification that can return lots of hits, and, even after experimentation, to correlate in silico data to condition-specific changes in regulatory patterns. Sun et al. also stated that it was not possible to predict whether a promoter would be bound by a TF solely on the presence of such motifs [96]. Sound conclusions on the real involvement of TFs hence require (global) experimental validations, as in *N. crassa* where ChIP-PCR of putative target genes showed that CRE-1 binds to only some adjacent 5'-SYGGRG-3' motifs, suggesting that additional regulatory factors affect the specificity of CRE-1 binding to the promoter [96].

With such a variety of data and regulatory patterns, even for those obtained from closely related fungi as *P. funiculosum*, it is difficult to reconcile some of these observations. An example is the positive effect of glucose on cellulase basal expression (our study), and conversely, a basal expression that is negatively regulated by PfMig1 in the presence of glucose [98]. Mechanistically, such apparent discrepancies may rely on genomic sequence differences between the strains, either acting indirectly on TFs localisation and stability, affinity to binding sites, or on local modifications of chromatin structure [53]. Hence, it might be interesting to decipher the exact interplay between glucose, CreA/Cre1/PfMig1 and possibly other TFs and chromatin modifiers. Our attempts at deleting CreA in *T. versatilis* were unfortunately unsuccessful, preventing further investigations on growth phenotypes and the analysis of regulatory patterns in this mutant strain. We faced a similar failure trying to delete the sole zinc finger domain of this transcription factor, while this mutation was successfully obtained in the *P. funiculosum* NCIM1228 strain by using a similar HR-based genetic manipulation [98]. Such technical bottleneck should be alleviated by the implementation of the CRISPR/cas9 system in *T. verstilis*, as it recently proved its efficiency in genome editing in many filamentous fungi species [106–108]. Either for new CreA deletion attempts, or for easier relief from CCR by promoter engineering, of either direct targets of interest or key TFs that might be themselves under the control of CCR [100, 101, 109].

**Gene redundancy in functional categories: the heterogeneity of expression patterns seems to be the rule**

A striking feature in the genome of filamentous fungi, including *T. versatilis*, is the remarkable redundancy in genes coding for the same enzymatic activity [25–28]. In our non-exhaustive set of GOIs, 8 encode xylanases, 7 of
them being members of the GH11 CAZY family. Even more noticeable, 10 encode arabinofuranosidases distributed in CAZY families GH43 (1), GH51 (1), GH54 (5) and GH62 (3), and 11 encode cellulases quite evenly distributed in as much as 7 CAZY families. This redundancy likely goes with slightly different catalytic properties of these enzymes, as for GH54 and GH62 type arabinofuranosidases [110], GH10 and GH11 xylanases [111, 112], or for the cellulolytic activities of multi-enzymatic complex system [113] from Penicillium funiculosum.

Data collected in this study clearly confirmed that classification of proteins according to enzymatic activity does not correlate, for the essential, to uniform expression patterns of the different members. Such a variability has been noticed for a long time, studying specific targets as *T. reesei* *xyn1* and *xyn2* as a function of several inducing substrates [114]. In our study, we found an incredible diversity of responses within categories, from variable fold-change values to opposite regulatory patterns between the genes, with strong dependence upon the nutritional environment. Also and despite similar FC values under a specific condition, 3- to 4-log differences in transcript levels were observed, leading possibly to great differences in the levels of expressed proteins, and ultimately in the corresponding hydrolytic activities for a given substrate. Finally, another source of variability came from the contribution of transcription factors, as XlnR, whose impact on gene expression was both gene-specific and culture condition-dependent.

In a similar work, Mikus et al. [115] specifically focused on the expression patterns of different class II hydrophobin genes from *Hypocrea atroviridis*, under various physiological conditions. They also found that the 10 hydrophobins display different patterns of response to these conditions, illustrating both divergent and redundant functions of these family members. Also, comparative studies between eight mycorrhizal species, with 3 up to 27 hydrophobin genes per genome, confirmed that the fine-tuning of the expression of these genes in various conditions differed considerably between genes in a species of interest, and further showed differentiation between species without following patterns that would be explained by species traits or ecology [116]. Taken altogether, if considering the whole set of gene-condition pair possibilities available in this study, even within the most homogeneous categories (e.g. cellulases) and particularly for the least homogeneous ones (e.g. arabino-furanosidases), the diversity of response seems to be the rule, while common and stable behaviour between members, the exception.

## Conclusions

Amongst the numerous microorganisms exhibiting the capability of biomass degradation, the filamentous fungi *A. niger* and *T. reesei* (*H. jecorina*) stand as pillars of hydrolytic enzymes production. In the present study, we confirmed the extraordinary potential of *Talaromyces versatilis* (formerly known as *Penicillium funiculosum*), another efficient enzyme producer used at the industrial scale for feed. With few exceptions of conserved regulatory patterns as compared to other enzyme producers, our results emphasized instead the outstanding complexity and divergence of responses, between gene family members, and beyond, with the corresponding categories of genes in more or less distant fungi. This plasticity of gene expression clearly relied on the strong influence of both the inducer and the promoter context, for which lots of factors should ideally be considered in addition to dedicated TFs binding sites, as epigenetic effects, sequence nucleotide environment or presence of associated motifs (review from [117]). Further, on the evidence that TFs rarely trigger transcriptional responses as lonesome and direct effectors, but rather act in concert with multiple factors [34], all these elements actually preclude any anticipation of the transcriptional response of a gene, or subsets of genes, in a specific nutritional context. The rational engineering of a fungus of interest and/or the setup of a new biotechnological process to reach optimized, if not customized expression patterns of enzymes, hence appear almost impossible just relying on published data that can lead, in the best cases, to approximate trends. Preliminary, experimental assays carried out in the context of interest, seem therefore mandatory before thinking in (genetic) strategies for the improvement of enzyme production in fungi.

## Methods

### Strains and culture conditions

The industrial strain *Talaromyces versatilis* IMI 378536, deposited in the International Mycological Institute Strain Collection, is a property of ADISSEO FRANCE SAS (patent W0 99/57325 and US 2003/0108642 A1). It has been isolated after UV mutagenesis from the wild type IMI 134755 strain, and is used for the industrial production of the Rovabio® Excel enzymatic cocktail that is enriched for endoxylanases and cellulases [35]. Strains with deletion of *xlnR* (*ΔxlnR*) or overexpression of this gene (referred as *xlnR*⁺ in this study, *id.* to DSM26702 used for the industrial production of the Rovabio® Advance cocktail) were obtained from the IMI 378536 strain, following the protoplasts transformation method described in [67]. The analysis of *T. versatilis* *xlnR*⁺ genomic DNA by southern blot and digital droplet PCR, (carried out on the Get platform from Génotoul), showed
that this mutant strain carries nine copies of *xlnR* under the control of its own promoter, with random integration of a *xlnR* PCR fragment in its genome.

Spores were obtained by growing *T. versatilis* strains on Potato Dextrose Agar (PDA) plates. They were used to inoculate a minimal medium (MM) that contained for 1 L: 1.9 g KH₂PO₄, 0.65 g KCl, 0.65 g MgSO₄, 12.5 mg ZnSO₄, 12.5 mg MnCl₂, 12.5 mg FeSO₄, 5 g NH₄Cl. Unless otherwise stated, the pH was adjusted to 6.0 with 50 mM KH₂PO₄ and MM contained glucose at 10 g/L and was supplemented with 10 mM uridine to allow growth of the *T. versatilis* *ΔxlnR* since this strain is deleted for the *pyrG* locus. This medium was inoculated with 2 × 10⁵ spores/mL and the cultures was performed at 30 °C in Erlenmeyer flasks on a rotary shaker set at 150 rpm.

To expose *T. versatilis* strains to different nutritional conditions, the mycelium obtained after 48 h in MM broth with glucose was filtered through Miracloth (Merck), washed with MM without carbon source and then transferred for the indicated time to fresh media containing the desired carbon source (Table 2). The wild-type and *xlnR*⁺ *T. versatilis* strains were also grown under lab-scale industrial condition that is used for the Rovabio® production, as described in [37]. This fermentation medium is a mix of cellulose and hemicellulose specifically designed for the production of the enzymatic cocktail by *T. versatilis* at industrial scale.

## Mycelia samples, RNA extraction and cDNA synthesis

For transcriptomic analysis, samples of about 50 mg of mycelium were collected by filtration through Miracloth and flash frozen in liquid nitrogen. These frozen samples were mechanically disrupted using the TissueLyser II (Qiagen) in the presence of a 5 mm stainless steel bead (Qiagen), with two sequential high-speed shaking, 30 Hz-set cycles of 3 min. Total RNA was isolated from disrupted mycelia samples using the GeneJET Plant RNA Purification Mini Kit (Thermo-Fischer). An on-column DNase I treatment (Thermo) was added to the protocol, applying 100 µL of the DNase I mix (50 µL of DNase I, 10 µL of 10 X buffer and 40 µL of nuclease-free water) to the column after the first wash, for a 30-min incubation at room temperature and final wash with the wash buffer I. The remaining of the protocol, including a second wash and the elution, was performed as recommended by the Supplier.

RNA was quantified using the ND-1000 UV-visible light spectrophotometer (NanoDrop Technologies) and its quality was assessed on a Bioanalyzer 2100 using the RNA 6000 Nano Labchip kit (Agilent). Only RNA samples with 260/280 nm wavelength ratio of approximately 2 and 260/230 nm wavelength ratio greater than 2 were

retained for analysis. One microgram of total RNA from mycelia samples were used for the cDNA synthesis reaction using the PrimeScript First Strand cDNA Synthesis Kit (Takara), according to the Manufacturers' protocol. The cDNA was diluted 1:10 with water and stored at -20 °C.

## Gene selection and primers design

The choice of a total of 62 genes of interest (GOIs) that encode representative glycoside hydrolases of different GH families is described in the results section. This list of GOIs, manually curated (genome annotation from A. Llanos Ph.D. thesis, (http://www.theses.fr/2014ISAT0029), proprietary data from ADISSEO France SAS), is reported in Table 2. Besides GOIs, 6 other genes, (ubc6, sac7, psm1, npl1, spo7 and spt3), were used as reference genes for RT-qPCR data quantification and normalization [59].

Primers were designed using Vector NTI advance v11 (Life Technologies) with melting temperature of 58-60 °C, length of 18-25 bp and GC content of 50-60%. Most of the genes possess one to several introns in their sequence, which allowed designing the primers at the exon-exon junctions to minimise the amplification of contaminant gDNA. Specificity of the primers for each transcript was checked by BLAST analysis against the *T. versatilis* genome. Amplicon sizes ranged between 70 and 230 bp. Reaction efficiency for each pair of primers was tested by the dilution series method using a mix of cDNA samples as the template. The detailed list of the genes, with corresponding sequences and primers, as well as primer efficiencies and amplicon lengths, is available in the Additional file 1.

## Quantitative PCR on the Biomark

High-throughput quantitative PCR (qPCR) was performed using the Fluidigm® Biomark microfluidic system. Every sample-gene combination was analysed using 96.96 Dynamic Array™ Integrated Fluidic Circuits (IFCs), which ensure performing up to 9216 parallel qPCR reactions in nL-scale volumes. Prior to qPCR, the cDNA samples were pre-amplified (activation at 95 °C for 10 min and 14 PCR cycles of 95 °C for 15 s and 60 °C for 4 min) with the PreAmp Master Mix (Fluidigm®) and a primer mix at 200 nM, which contains all the primers targeting all the genes analyzed on the array. The samples were treated with Exonuclease I (New England BioLabs), for 30 min at 37 °C; 15 min at 80 °C; and 4 °C hold, and then diluted by adding 18 µL of TE-low EDTA buffer (10 mM Tris-HCl, 0.1 mM EDTA, pH 8). Two microliters of these pre-amplified cDNA samples were mixed with 6 µl of reagent mix into the 96-wells of the Sample plate. This reagent mix consisted of 440 µL 2X TaqMan Master Mix
<table><caption>**Table 2** Growth conditions</caption><thead><tr><th>Condition</th><th>Sample</th><th>Harvest time (h)</th><th>Medium composition</th><th>Substrate (reference)</th></tr></thead><tbody><tr><td>Glucose</td><td>AL1</td><td>48</td><td>MM + glucose 1%</td><td></td></tr><tr><td>Cellobiose</td><td>AL3</td><td>2</td><td>MM + cellobiose 0.2%</td><td>D-(+)-Cellobiose (Sigma Aldrich—22150)</td></tr><tr><td>Avicel</td><td>AL18</td><td>24</td><td>MM + avicel 1%</td><td>Avicel PH101 (Sigma Aldrich—11365)</td></tr><tr><td>xylan</td><td>AL20</td><td>24</td><td>MM + xylan 1%</td><td>Beechwood xylan (Sigma Aldrich—X4252)</td></tr><tr><td>wheat straw<sup>a</sup></td><td>AL22</td><td>24</td><td>MM + wheat straw 1%</td><td>Ball-milled wheat straw</td></tr><tr><td>N starvation</td><td>AL29</td><td>1</td><td>MM w/o nitrogen + glucose 1%</td><td></td></tr><tr><td>C starvation</td><td>AL30</td><td>1</td><td>MM (w/o carbon source)</td><td></td></tr><tr><td>Xylose 0.2%</td><td>AL38</td><td>2</td><td>MM + xylose 0.2%</td><td>D-(+)-Xylose (Sigma Aldrich no 95729)</td></tr><tr><td>Xylose 1%</td><td>AL39</td><td>1</td><td>MM + xylose 1%</td><td></td></tr><tr><td>Arabinose 0.2%</td><td>AL41</td><td>2</td><td>MM + arabinose 0.2%</td><td>D-(L)-arabinose (Sigma Aldrich—10850)</td></tr><tr><td>Arabinose 1%</td><td>AL42</td><td>1</td><td>MM + arabinose 1%</td><td></td></tr><tr><td>Soybean<sup>b</sup></td><td>AL43</td><td>24</td><td>MM + soybean 1%</td><td>Ball-milled soybean press cake</td></tr><tr><td>Fermentation<sup>c</sup></td><td>AL49</td><td>140</td><td>Cellulose and hemicellulose mix</td><td></td></tr></tbody></table>

(Applied Biosystems - 430976), 44 µL 20 × DNA Binding Dye Sample Loading Reagent (Fluidigm® - 100-3738), 44 µL 20X Evagreen (Biotium - 31000) plus 132 µL TE. The Sample plate was briefly vortexed and centrifuged. Relative to primer sets preparation, each pair of primers was prepared as a 5 µM solution (each primer) in the Assay loading reagent (Fluidigm®) and dispensed into the 96-wells of the Assay plate. Following priming of the Dynamic Array™ IFC in the IFC Controller HX, 5 µL of cDNA samples and primer sets from the Sample and Assay plates, respectively, were dispensed into the inlets of the IFC. After loading the Array using the IFC controller, the IFC was placed into the Biomark for qPCR. The following protocol was used: Thermal mix steps for 2 min at 50 °C; 30 min at 70 °C; 10 min at 25 °C; Hot start for 2 min at 50 °C and 10 min at 95 °C; followed by 35 cycles of 15 s at 95 °C and 1 min at 60 °C. A melting curve at the end of the protocol verified the specificity of the PCR, absence of contamination and of primer dimers. Data was collected with Fluidigm® Real-Time PCR analysis software using the linear baseline correction method and global auto Cq threshold method. Cq values that exceeded 35 cycles were considered beyond the limit of detection. Each array included a negative control (NTC), a positive control for the samples (human gDNA) and a control for the primers (RNAseP).

### Fold-change calculation and statistical analysis

Quantification of the transcripts by RT-qPCR was performed from three independent biological replicates (samples) for each of the 13 culture conditions. A technical duplicate was carried out for each sample, except

for the glucose samples that were analyzed in technical quadruplicates due to the low transcript levels of most of the GOIs under this condition. Raw qPCR signals were analyzed using the Real-Time PCR Analysis software (Fluidigm®) to determine Cq values. These Cq were then used to calculate the Fold-Change (FC) values, taking into consideration each primer pair efficiency. Data normalization was performed using a multiple reference genes procedure, with the geometrical mean of 6 validated reference genes for *T. versatilis*, i.e. *ubc6*, *sac7*, *psm1*, *npl1*, *spo7* and *spt3* [59]. FC statistics (mean ± SD) were determined from the 3 biological replicates. Unless otherwise stated, the values were expressed as log₂ (FC). Statistical analyses and advanced graphics were conducted using R (https://www.r-project.org/). The full set of expression data used in this article, for the three strain backgrounds, is provided in the Additional file 8.

### Upstream regulatory motifs analysis

Upstream regions of the genes of interest, (700 bp upstream of the ATG codon), were analysed for the presence of cis-acting elements using RSAT Fungi (Regulatory Sequence Analysis Tool—http://rsat-tagc.univ-mrs.fr/rsat/). Each complementary strand of the promoter was scanned for the motifs 5'-GGCTAA-3' (XlnR), 5'-GGCWAW-3' (alternative XlnR motif) and 5'-SCG-GRG-3' (CreA) (Additional file 1). The analysis was also performed with the MEME suite (Multiple EM for Motif Elicitation—http://meme-suite.org) to confirm the results.

<footnote>a [67]</footnote>

<footnote>b [118]</footnote>

<footnote>c Lab-scale industrial condition, carried out at 30 °C, 500 RPM, Ph 4 (proprietary medium composition, ADISSEO FRANCE SAS)</footnote>
## Additional files

**Additional file 1.** Detailed information on GOs. Gene annotations, qPCR primers specifications and CDS sequence, number of each of TF binding motifs in the promoter region, FC values for the 'glucose repression' condition (extracted from our previous RNA-seq study [59]).

**Additional file 2.** Expression profiles in gene clusters, in the WT strain, as a function of culture conditions. Radar chart representing the expression level of GOs ($log_2$ (FC) values), in each of the five gene clusters defined in Fig. 2. Each axis corresponds to a culture condition, as defined in the bottom right legend of the figure. The coloured line represents the average $log_2$ (FC) value of GOs gathered in the cluster.

**Additional file 3.** GOs expression data in each gene cluster. Box plot illustration of $log_2$ (FC) values as a function of the culture conditions and in the three strain backgrounds (grey for the WT, green for the $\Delta xlnR$ mutant, and red for the $xlnR^+$ strain). The ordering of conditions (from the left (Straw) to the right (N-starv)) relies on Fig. 2 classification (top-down order).

**Additional file 4.** GOs expression data in each functional category. Box plot illustration of $log_2$ (FC) values as a function of the culture conditions and in the WT, $\Delta xlnR$ and $xlnR^+$ strains. Legend as in Additional file 3.

**Additional file 5.** Specificity of the $\Delta xlnR$ response patterns on xylose, xylan and straw. Substrate-dependent scatter plots of $log_2$ (FC) values in the $\Delta xlnR$ strain as a function of $log_2$ (FC) values in the wild type strain. The GOs are gathered according to the five gene clusters that have been defined in Fig. 2.

**Additional file 6.** Specificity of the $xlnR^+$ response patterns on xylose, xylan and straw. Substrate-dependent scatter plots of $log_2$ (FC) values in the $xlnR^+$ strain as a function of $log_2$ (FC) values in the wild type strain. The GOs are gathered according to the five gene clusters that have been defined in Fig. 2.

**Additional file 7.** Pairwise correlations between genes encoding backbone hydrolases, in the three strain backgrounds. These matrices were drawn using expression values from Fig. 2 (a: WT strain), Fig. 4 (b: $\Delta xlnR$ strain) and Fig. 5 (c: $xlnR^+$ strain). Same legend as in Fig. 7. We arbitrary coloured in green the pairs that exhibited rather good pairwise correlations in the WT strain (r > 0.90).

**Additional file 8.** Data file of expression data from the study. This lxsx file contains all calculated Log2(FC) values from this study, used for graphics preparation, and that can be readily imported in dedicated statistical analysis softwares (e.g. R, Prism, Statgraphics, etc.). Use the \*Read Me\* sheet for detailed information on the organisation of the data in the different sheets of the file.

## Abbreviations

CAZyme: carbohydrate-active enzymes; CCR: carbon catabolite repression; Cq: quantification cycle in qPCR; C-starv: carbon starvation; ddPCR: digital droplet PCR; FC: fold change; GH: glycosyl hydrolase; GOI: Gene Of Interest; HR-based genetic manipulation: Homologous Recombination-based genetic manipulation; $log_2$ (FC): logarithm base 2 of fold change values; MM: minimal medium; N-Starv: nitrogen starvation; qPCR: quantitative PCR; RNA-seq: whole cell RNA sequencing; RPKM: reads per kilobase of exon model per million mapped reads; RQ: relative quantity; TF: transcription factor.

## Authors' contributions

AL, VNR, JMF and JLP participated in the design of the study. AL carried out the experimental work. AL and JLP analysed the data. SD participated in the statistical analysis. AL drafted the manuscript. AL, SD and JLP contributed to graphics preparation. JLP wrote the manuscript. AL and JMF revised the manuscript. All authors read and approved the final manuscript.

## Author details

## Acknowledgements

We thank Olivier Guais (CINABo-Adisseo) for scientific support during this work and for providing unpublished, proprietary proteomics data.

## Competing interests

The authors declare that they have no competing interests.

## Availability of data and materials

All data generated or analysed during this study are included in additional files.

## Consent for publication

Not applicable.

## Ethics approval and consent to participate

Not applicable. The manuscript does not report data collected from humans or animals.

## Funding

This work was funded by grants from Adisseo (Adisseo SAS, Antony, France) to JLP and JMF, and by CIFRE, Grant Number: 1558/2010 to AL.

## Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Received: 9 November 2018 Accepted: 13 January 2019

Published online: 28 January 2019

## References

1. Segato F, Damasio AR, de Lucas RC, Squina FM, Prade RA. Genomics review of holocellulose deconstruction by *aspergilli*. *Microbiol Mol Biol Rev.* 2014;78:588–613.

2. Bomble YJ, Lin CY, Amore A, Wei H, Holwerda EK, Ciesielski PN, Donohoe BS, Decker SR, Lynd LR, Himmel ME. Lignocellulose deconstruction in the biosphere. *Curr Opin Chem Biol.* 2017;41:61–70.

3. Cragg SM, Beckham GT, Bruce NC, Bugg TD, Distel DL, Dupree P, Etxabe AG, Goodell BS, Jellison J, McGeehan JE, et al. Lignocellulose degradation mechanisms across the Tree of Life. *Curr Opin Chem Biol.* 2015;29:108–19.

4. Benoit I, Coutinho PM, Schols HA, Gerlach JP, Henrissat B, de Vries RP. Degradation of different pectins by fungi: correlations and contrasts between the pectinolytic enzyme sets identified in genomes and the growth on pectins of different origin. *BMC Genomics.* 2012;13:321.

5. Culleton H, McKie V, de Vries RP. Physiological and molecular aspects of degradation of plant polysaccharides by fungi: what have we learned from *Aspergillus*? *Biotechnol J.* 2013;8:884–94.

6. de Vries RP, Visser J. *Aspergillus* enzymes involved in degradation of plant cell wall polysaccharides. *Microbiol Mol Biol Rev.* 2001;65:497–522.

7. Makela MR, Donofrio N, de Vries RP. Plant biomass degradation by fungi. *Fungal Genet Biol.* 2014;72:2–9.

8. Eastwood DC, Floudas D, Binder M, Majcherczyk A, Schneider P, Aerts A, Asiegbu FO, Baker SE, Barry K, Bendiksby M, et al. The plant cell wall-decomposing machinery underlies the functional diversity of forest fungi. *Science.* 2011;333:762–5.

9. Lombard V, Golaconda Ramulu H, Drula E, Coutinho PM, Henrissat B. The carbohydrate-active enzymes database (CAZy) in 2013. *Nucleic Acids Res.* 2014;42:D490–5.

10. Collins T, Gerday C, Feller G. Xylanases, xylanase families and extremophilic xylanases. *FEMS Microbiol Rev.* 2005;29:3–23.

11. Martinez D, Larrondo LF, Putnam N, Gelpke MD, Huang K, Chapman J, Helfenbein KG, Ramaiya P, Detter JC, Larimer F, et al. Genome sequence of the lignocellulose degrading fungus *Panerochaete chrysosporium* strain RP78. *Nat Biotechnol.* 2004;22:695–700.

12. Machida M, Asai K, Sano M, Tanaka T, Kumagai T, Terai G, Kusumoto K, Arima T, Akita O, Kashiwagi Y, et al. Genome sequencing and analysis of *Aspergillus oryzae*. *Nature.* 2005;438:1157–61.

<footnote n="1">¹ LISBP, Université de Toulouse, INSA, INRA, CNRS, Toulouse, France. ² Adisseo France S.A.S, 135 Avenue de Rangueil, 31077 Toulouse, France. ³ Institut de Mathématiques de Toulouse, UMR5219-Université de Toulouse; CNRS-UPS, 31062 Toulouse Cedex 9, France.</footnote>
13. Pel HJ, de Winde JH, Archer DB, Dyer PS, Hofmann G, Schaap PJ, Turner G, de Vries RP, Albang R, Albermann K, et al. Genome sequencing and analysis of the versatile cell factory *Aspergillus niger* CBS 513.88. Nat Biotechnol. 2007;25:221–31.

14. Martinez D, Berka RM, Henrissat B, Saloheimo M, Arvas M, Baker SE, Chapman J, Chertkov O, Coutinho PM, Cullen D, et al. Genome sequencing and analysis of the biomass-degrading fungus *Trichoderma reesei* (syn. *Hypocrea jecorina*). Nat Biotechnol. 2008;26:553–60.

15. Yang Y, Chen M, Li Z, Al-Hatmi AM, Ye Q, Chen H, Wang S, Liao W, Wang J. Genome sequence of *Penicillium capsulatum* strain ATCC 48735, a rare *Penicillium* species used in paper manufactories but that recently caused invasive infection. Genome Announc. 2015;3:e00307–15.

16. Hu L, Taujale R, Liu F, Song J, Yin Q, Zhang Y, Guo J, Yin Y. Draft genome sequence of *Talaromyces verruculosus* (*"Penicillium verruculosum"*) strain TS63-9, a fungus with great potential for industrial production of polysaccharide-degrading enzymes. J Biotechnol. 2016;219:5–6.

17. Gaskell J, Kersten P, Larrondo LF, Canessa P, Martinez D, Hibbett D, Schmoll M, Kubicek CP, Martinez AT, Yadav J, et al. Draft genome sequence of a monokaryotic model brown-rot fungus *Postia (Rhodonia) placenta* SB12. Genom Data. 2017;14:21–31.

18. Peng M, Dilokpimol A, Makela MR, Hilden K, Bervoets S, Riley R, Grigoriiev IV, Hainaut M, Henrissat B, de Vries RP, Granchi Z. The draft genome sequence of the ascomycete fungus *Penicillium subrubescens* reveals a highly enriched content of plant biomass related CAZymes compared to related fungi. J Biotechnol. 2017;246:1–3.

19. Verma S, Gazara RK, Nizam S, Parween S, Chattopadhyay D, Verma PK. Draft genome sequencing and secretome analysis of fungal phytopathogen *Ascochyta rabiei* provides insight into the necrotrophic effector repertoire. Sci Rep. 2016;6:24638.

20. Cerqueira GC, Arnaud MB, Inglis DO, Skrzypek MS, Binkley G, Simison M, Miyasato SR, Binkley J, Orvis J, Shah P, et al. The Aspergillus Genome Database: multispecies curation and incorporation of RNA-Seq data to improve structural gene annotations. Nucleic Acids Res. 2014;42:D705–10.

21. Makela MR, Mansouri S, Wiebenga A, Rytioja J, de Vries RP, Hilden KS. *Penicillium subrubescens* is a promising alternative for *Aspergillus niger* in enzymatic plant biomass saccharification. N Biotechnol. 2016;33:834–41.

22. Suzuki H, MacDonald J, Syed K, Salamov A, Hori C, Aerts A, Henrissat B, Wiebenga A, VanKuyk PA, Barry K, et al. Comparative genomics of the white-rot fungi, *Phanerochaete carnosa* and *P. chrysosporium*, to elucidate the genetic basis of the distinct wood types they colonize. BMC Genomics. 2012;13:444.

23. Galagan JE, Calvo SE, Borkovich KA, Selker EU, Read ND, Jaffe D, FitzHugh W, Ma LJ, Smirnov S, Purcell S, et al. The genome sequence of the filamentous fungus *Neurospora crassa*. Nature. 2003;422:859–68.

24. Galagan JE, Henn MR, Ma LJ, Cuomo CA, Birren B. Genomics of the fungal kingdom: insights into eukaryotic biology. Genome Res. 2005;15:1620–31.

25. Coutinho PM, Andersen MR, Kolenova K, vanKuyk PA, Benoit I, Gruben BS, Trejo-Aguilar B, Visser H, van Solingen P, Pakula T, et al. Post-genomic insights into the plant polysaccharide degradation potential of *Aspergillus nidulans* and comparison to *Aspergillus niger* and *Aspergillus oryzae*. Fungal Genet Biol. 2009;46(Suppl 1):S161–9.

26. Lyu X, Shen C, Fu Y, Xie J, Jiang D, Li G, Cheng J. Comparative genomic and transcriptional analyses of the carbohydrate-active enzymes and secretomes of phytopathogenic fungi reveal their significant roles during infection and development. Sci Rep. 2015;5:15565.

27. Zhao Z, Liu H, Wang C, Xu JR. Correction: Comparative analysis of fungal genomes reveals different plant cell wall degrading capacity in fungi. BMC Genomics. 2014;15:6.

28. Couturier M, Navarro D, Olive C, Chevret D, Haon M, Favel A, Lesage-Meessen L, Henrissat B, Coutinho PM, Berrin JG. Post-genomic analyses of fungal lignocellulosic biomass degradation reveal the unexpected potential of the plant pathogen *Ustilago maydis*. BMC Genomics. 2012;13:57.

29. Krijger JJ, Thon MR, Deising HB, Wirsel SG. Compositions of fungal secretomes indicate a greater impact of phylogenetic history than lifestyle adaptation. BMC Genomics. 2014;15:722.

30. Geiser E, Reindl M, Blank LM, Feldbrugge M, Wierck N, Schipper K. Activating intrinsic carbohydrate-active enzymes of the smut fungus

*Ustilago maydis* for the degradation of plant cell wall components. Appl Environ Microbiol. 2016;82:5174–85.

31. Dean RA, Talbot NJ, Ebbole DJ, Farman ML, Mitchell TK, Orbach MJ, Thon M, Kulkarni R, Xu JR, Pan H, et al. The genome sequence of the rice blast fungus *Magnaporthe grisea*. Nature. 2005;434:980–6.

32. Martinez D, Challacombe J, Morgenstern I, Hibbett D, Schmoll M, Kubicek CP, Ferreira P, Ruiz-Duenas FJ, Martinez AT, Kersten P, et al. Genome, transcriptome, and secretome analysis of wood decay fungus *Postia placenta* supports unique mechanisms of lignocellulose conversion. Proc Natl Acad Sci USA. 2009;106:1954–9.

33. Alfaro M, Castanera R, Lavin JL, Grigoriev IV, Oguiza JA, Ramirez L, Pisabarro AG. Comparative and transcriptional analysis of the predicted secretome in the lignocellulose-degrading basidiomycete fungus. Environ Microbiol. 2016;18:4710–26.

34. Tani S, Kawaguchi T, Kobayashi T. Complex regulation of hydrolytic enzyme genes for cellulosic biomass degradation in filamentous fungi. Appl Microbiol Biotechnol. 2014;98:4829–37.

35. Aquilina G, Bampidis V, De Lourdes Bastos M, Costa LG, Flachowsky G, Gralak MA, Hogstrand C, Leng L, López-Puente S, Martelli G, et al. Scientific opinion on Rovabio® Excel (endo-1, 3(4)-beta-glucanase and endo-1, 4-beta-xylanase) as a feed additive for chickens and turkeys for fattening, laying hens, piglets (weaned) and pigs for fattening, ducks, guinea fowls, quails, geese, pheasants and pigeons. EFSA J. 2013;11:26.

36. Cozannet P, Kidd MT, Montanhini Neto R, Geraert PA. Next-generation non-starch polysaccharide-degrading, multi-carbohydrase complex rich in xylanase and arabinofuranosidase to enhance broiler feed digestibility. Poult Sci. 2017;96:2743–50.

37. Guais O, Borderies G, Pichereaux C, Maestracci M, Neugnot V, Ros-signol M, Francois JM. Proteomics analysis of "Rovabiot Excel", a secreted protein cocktail from the filamentous fungus *Penicillium funiculosum* grown under industrial process fermentation. J Ind Microbiol Biotechnol. 2008;35:1659–68.

38. Ogunmolu FE, Kaur I, Gupta M, Bashir Z, Pasari N, Yazdani SS. Proteomics insights into the biomass hydrolysis potentials of a hypercellulolytic fungus *Penicillium funiculosum*. J Proteome Res. 2015;14:4342–58.

39. Ogunmolu FE, Kaur I, Pasari N, Gupta M, Yazdani SS. Quantitative multiplexed profiling of *Penicillium funiculosum* secretome grown on polymeric cellulose inducers and glucose. J Proteomics. 2018;179:150–60.

40. McCotter SW, Horianopoulos LC, Kronstad JW. Regulation of the fungal secretome. Curr Genet. 2016;62:533–45.

41. Benoit I, Culleton H, Zhou M, DiFalco M, Aguilar-Osorio G, Battaglia E, Bouzid O, Brouwer C, El-Bushari HBO, Coutinho PM, et al. Closely related fungi employ diverse enzymatic strategies to degrade plant biomass. Biotechnol Biofuels. 2015;8:107.

42. Cologna NMD, Gomez-Mendoza DP, Zanoelo FF, Giannesi GC, Guima-raes NCA, Moreira LRS, Filho EXF, Ricart CAO. Exploring *Trichoderma* and *Aspergillus* secretomes: proteomics approaches for the identification of enzymes of biotechnological interest. Enzyme Microb Technol. 2018;109:1–10.

43. Wymelenberg AV, Gaskell J, Mozuch M, Sabat G, Ralph J, Skyba O, Mansfield SD, Blanchette RA, Martinez D, Grigoriev I, et al. Comparative transcriptome and secretome analysis of wood decay fungi *Postia placenta* and *Phanerochaete chrysosporium*. Appl Environ Microbiol. 2010;76:3599–610.

44. Gaskell J, Blanchette RA, Stewart PE, BonDurant SS, Adams M, Sabat G, Kersten P, Cullen D. Transcriptome and secretome analyses of the wood decay fungus *Wolfiporia cocos* support alternative mechanisms of lignocellulose conversion. Appl Environ Microbiol. 2016;82:3979–87.

45. Borin GP, Sanchez CC, de Santana ES, Zanini GK, Dos Santos RAC, de Oliveira Pontes A, de Souza AT, Dal'Mas R, Riano-Pachon DM, Goldman GH, Oliveira JVC. Comparative transcriptome analysis reveals different strategies for degradation of steam-exploded sugarcane bagasse by *Aspergillus niger* and *Trichoderma reesei*. BMC Genomics. 2017;18:501.

46. Florencio C, Cunha FM, Badino AC, Farinas CS, Ximenes E, Ladisch MR. Secretome analysis of *Trichoderma reesei* and *Aspergillus niger* cultivated by submerged and sequential fermentation processes: enzyme production for sugarcane bagasse hydrolysis. Enzyme Microb Technol. 2016;90:53–60.
47. Borin GP, Sanchez CC, de Souza AP, de Santana ES, de Souza AT, Paes Leme AF, Squina FM, Buckeridge M, Goldman GH, Oliveira JV. Comparative secretome analysis of *Trichoderma reesei* and *Aspergillus niger* during growth on sugarcane biomass. PLoS ONE. 2015;10:e0129275.

48. Gong W, Zhang H, Liu S, Zhang L, Gao P, Chen G, Wang L. Comparative secretome analysis of *Aspergillus niger*, *Trichoderma reesei*, and *Penicillium oxalicum* during solid-state fermentation. Appl Biochem Biotechnol. 2015;177:1252–71.

49. Schneider WDH, Goncalves TA, Uchima CA, Couger MB, Prade R, Squina FM, Dillon AJP, Camassola M. *Penicillium echinulatum* secretome analysis reveals the fungi potential for degradation of lignocellulosic biomass. Biotechnol Biofuels. 2016;9:66.

50. van Munster JM, Daly P, Delmas S, Pullan ST, Blythe MJ, Malla S, Kokolski M, Noltorp ECM, Wennberg K, Fetherston R, et al. The role of carbon starvation in the induction of enzymes that degrade plant-derived carbohydrates in *Aspergillus niger*. Fungal Genet Biol. 2014;72:34–47.

51. Amore A, Giacobbe S, Faraco V. Regulation of cellulase and hemicellulase gene expression in fungi. Curr Genomics. 2013;14:230–49.

52. Aro N, Pakula T, Penttila M. Transcriptional regulation of plant cell wall degradation by filamentous fungi. FEMS Microbiol Rev. 2005;29:719–39.

53. Benocci T, Aguilar-Pontes MV, Zhou M, Seiboth B, de Vries RP. Regulators of plant biomass degradation in ascomycetous fungi. Biotechnol Biofuels. 2017;10:152.

54. Brown NA, Ries LNA, Goldman GH. How nutritional status signalling coordinates metabolism and lignocellulolytic enzyme secretion. Fungal Genet Biol. 2014;72:48–63.

55. van Peij NN, Visser J, de Graaff LH. Isolation and analysis of *xlnR*, encoding a transcriptional activator co-ordinating xylanolytic expression in *Aspergillus niger*. Mol Microbiol. 1998;27:131–42.

56. Andersen MR, Vongsangnak W, Panagiotou G, Salazar MP, Lehmann L, Nielsen J. A trispecies *Aspergillus* microarray: comparative transcriptomics of three *Aspergillus* species. Proc Natl Acad Sci USA. 2008;105:4387–92.

57. Klaubauf S, Narang HM, Post H, Zhou M, Brunner K, Mach-Aigner AR, Mach RL, Heck AJR, Altelaar AFM, de Vries RP. Similar is not the same: differences in the function of the (hemi-)cellulolytic regulator XlnR (Xlr1/Xyr1) in filamentous fungi. Fungal Genet Biol. 2014;72:73–81.

58. Battaglia E, Klaubauf S, Vallet J, Ribot C, Lebrun MH, de Vries RP. Xlr1 is involved in the transcriptional control of the pentose catabolic pathway, but not hemi-cellulolytic enzymes in *Magnaporthe oryzae*. Fungal Genet Biol. 2013;57:76–84.

59. Llanos A, Francois JM, Parrou JL. Tracking the best reference genes for RT-qPCR data normalization in filamentous fungi. BMC Genomics. 2015;16:71.

60. Furukawa T, Shida Y, Kitagami N, Mori K, Kato M, Kobayashi T, Okada H, Ogasawara W, Morikawa Y. Identification of specific binding sites for XRY1, a transcriptional activator of cellulolytic and xylanolytic genes in *Trichoderma reesei*. Fungal Genet Biol. 2009;46:564–74.

61. Marui J, Tanaka A, Mimura S, de Graaff LH, Visser J, Kitamoto N, Kato M, Kobayashi T, Tsukagoshi N. A transcriptional activator, AoXlnR, controls the expression of genes encoding xylanolytic enzymes in *Aspergillus oryzae*. Fungal Genet Biol. 2002;35:157–69.

62. Cubero B, Scazzocchio C. Two different, adjacent and divergent zinc finger binding sites are necessary for CREA-mediated carbon catabolite repression in the proline gene cluster of *Aspergillus nidulans*. EMBO J. 1994;13:407–15.

63. Ilmen M, Thrane C, Penttila M. The glucose repressor gene *cre1* of *Trichoderma*: isolation and expression of a full-length and a truncated mutant form. Mol Gen Genet. 1996;251:451–60.

64. Kumar L, Breakspear A, Kistler C, Ma LJ, Xie X. Systematic discovery of regulatory motifs in *Fusarium graminearum* by comparing four *Fusarium* genomes. BMC Genomics. 2010;11:208.

65. Delmas S, Pullan ST, Gaddipati S, Kokolski M, Malla S, Blythe MJ, Ibbett R, Campbell M, Liddell S, Aboobaker A, et al. Uncovering the genome-wide transcriptional responses of the filamentous fungus *Aspergillus niger* to lignocellulose using RNA sequencing. PLoS Genet. 2012;8:e1002875.

66. Ries L, Pullan ST, Delmas S, Malla S, Blythe MJ, Archer DB. Genome-wide transcriptional response of *Trichoderma reesei* to lignocellulose using RNA sequencing and comparison with *Aspergillus niger*. BMC Genomics. 2013;14:541.

67. Delmas S, Llanos A, Parrou JL, Kokolski M, Pullan ST, Shunburne L, Archer DB. Development of an unmarked gene deletion system for the filamentous fungi *Aspergillus niger* and *Talaromyces versatilis*. Appl Environ Microbiol. 2014;80:3484–7.

68. Brunner K, Lichtenauer AM, Kratochwill K, Delic M, Mach RL. Xyr1 regulates xylanase but not cellulase formation in the head blight fungus *Fusarium graminearum*. Curr Genet. 2007;52:213–20.

69. Marui J, Kitamoto N, Kato M, Kobayashi T, Tsukagoshi N. Transcriptional activator, AoXlnR, mediates cellulose-inductive expression of the xylanolytic and cellulolytic genes in *Aspergillus oryzae*. FEBS Lett. 2002;528:279–82.

70. Sun J, Tian C, Diamond S, Glass NL. Deciphering transcriptional regulatory mechanisms associated with hemicellulose degradation in *Neurospora crassa*. Eukaryot Cell. 2012;11:482–93.

71. Noguchi Y, Sano M, Kanamaru K, Ko T, Takeuchi M, Kato M, Kobayashi T. Genes regulated by AoXlnR, the xylanolytic and cellulolytic transcriptional regulator, in *Aspergillus oryzae*. Appl Microbiol Biotechnol. 2009;85:141–54.

72. van Peij NN, Gielkens MM, de Vries RP, Visser J, de Graaff LH. The transcriptional activator XlnR regulates both xylanolytic and endoglucanase gene expression in *Aspergillus niger*. Appl Environ Microbiol. 1998;64:3615–9.

73. Wang J, Wu Y, Gong Y, Yu S, Liu G. Enhancing xylanase production in the thermophilic fungus *Myceliophthora thermophila* by homologous overexpression of *Mtxyr1*. J Ind Microbiol Biotechnol. 2015;42:1233–41.

74. Tamayo-Ramos JA, Orejas M. Enhanced glycosyl hydrolase production in *Aspergillus nidulans* using transcription factor engineering approaches. Biotechnol Biofuels. 2014;7:103.

75. Okuda N, Fujii T, Inoue H, Ishikawa K, Hoshino T. Enhancing cellulase production by overexpression of xylanase regulator protein gene, *xlnR*, in *Talaromyces cellulolyticus* cellulase hyperproducing mutant strain. Biosci Biotechnol Biochem. 2016;80:2065–8.

76. Margolles-Clark E, Ihnen M, Penttila M. Expression patterns of ten hemicellulase genes of the filamentous fungus *Trichoderma reesei* on various carbon sources. J Biotechnol. 1997;57:13.

77. dos Santos Castro L, Antonieto AC, Pedersoli WR, Silva-Rocha R, Persi-noti GF, Silva RN. Expression pattern of cellulolytic and xylanolytic genes regulated by transcriptional factors XRY1 and CRE1 are affected by carbon source in *Trichoderma reesei*. Gene Expr Patterns. 2014;14:88–95.

78. Hakkinen M, Arvas M, Oja M, Aro N, Penttila M, Saloheimo M, Pakula TM. Re-annotation of the CAZy genes of *Trichoderma reesei* and transcription in the presence of lignocellulosic substrates. Microb Cell Fact. 2012;11:134.

79. Pullan ST, Daly P, Delmas S, Ibbett R, Kokolski M, Neiteler A, van Munster JM, Wilson R, Blythe MJ, Gaddipati S, et al. RNA-sequencing reveals the complexities of the transcriptional response to lignocellulosic biofuel substrates in *Aspergillus niger*. Fungal Biol Biotechnol. 2014;1:1–14.

80. Liao H, Li S, Wei Z, Shen Q, Xu Y. Insights into high-efficiency lignocellulolytic enzyme production by *Penicillium oxalicum* GZ-2 induced by a complex substrate. Biotechnol Biofuels. 2014;7:162.

81. Stricker AR, Mach RL, de Graaff LH. Regulation of transcription of cellulases- and hemicellulases-encoding genes in *Aspergillus niger* and *Hypocrea jecorina* (*Trichoderma reesei*). Appl Microbiol Biotechnol. 2008;78:211–20.

82. Daly P, van Munster JM, Blythe MJ, Ibbett R, Kokolski M, Gaddipati S, Lindquist E, Singan VR, Barry KW, Lipzen A, et al. Expression of *Aspergillus niger* CAZymes is determined by compositional changes in wheat straw generated by hydrothermal or ionic liquid pretreatments. Biotechnol Biofuels. 2017;10:35.

83. Lichius A, Seidl-Seiboth V, Seiboth B, Kubicek CP. Nucleo-cytoplasmic shuttling dynamics of the transcriptional regulators *XYR1* and *CRE1* under conditions of cellulase and xylanase gene expression in *Trichoderma reesei*. Mol Microbiol. 2014;94:1162–78.

84. Portnoy T, Margeot A, Seidl-Seiboth V, Le Crom S, Ben Chaabane F, Linke R, Seiboth B, Kubicek CP. Differential regulation of the cellulase transcription factors XRY1, ACE2, and ACE1 in *Trichoderma reesei* strains producing high and low levels of cellulase. Eukaryot Cell. 2011;10:262–71.

85. Bouzid O, de Vries RP. Overexpression of XlnR regulated hemicellulase-encoding genes in *Aspergillus vadosensis*. Curr Biotechnol. 2014;3:252–6.
86. Calero-Nieto F, Di Pietro A, Roncero MI, Hera C. Role of the transcriptional activator XlnR of *Fusarium oxysporum* in regulation of xylanase genes and virulence. Mol Plant Microbe Interact. 2007;20:977-85.

87. Noguchi Y, Tanaka H, Kanamaru K, Kato M, Kobayashi T. Xylose triggers reversible phosphorylation of XlnR, the fungal transcriptional activator of xylanolytic and cellulolytic genes in *Aspergillus oryzae*. *Biosci Biotechnol Biochem.* 2011;75:953-9.

88. Craig JP, Coradetti ST, Starr TL, Glass NL. Direct target network of the Neurospora crassa plant cell wall deconstruction regulators CLR-1, CLR-2, and XLR-1. MBio. 2015;6:e01452-515.

89. Kunitake E, Kobayashi T. Conservation and diversity of the regulators of cellulolytic enzyme genes in Ascomycete fungi. Curr Genet. 2017;63:951-8.

90. Fujii T, Inoue H, Ishikawa K. Characterization of the xylanase regulator protein gene, *xlnR*, in *Talaromyces cellulolyticus* (formerly known as *Acremonium cellulolyticus*). *Biosci Biotechnol Biochem.* 2014;78:1564-7.

91. Coradetti ST, Craig JP, Xiong Y, Shock T, Tian C, Glass NL. Conserved and essential transcription factors for cellulase gene expression in ascomycete fungi. Proc Natl Acad Sci USA. 2012;109:7397-402.

92. Znameroski EA, Coradetti ST, Roche CM, Tsai JC, Iavarone AT, Cate JH, Glass NL. Induction of lignocellulose-degrading enzymes in Neurospora crassa by cellodextrins. Proc Natl Acad Sci USA. 2012;109:6012-7.

93. Kunitake E, Tani S, Sumitani J, Kawaguchi T. A novel transcriptional regulator, ClbR, controls the cellobiose- and cellulose-responsive induction of cellulase and xylanase genes regulated by two distinct signaling pathways in *Aspergillus aculeatus*. *Appl Microbiol Biotechnol.* 2013;97:2017-28.

94. Bischof R, Fourlis L, Limbeck A, Gamauf C, Seiboth B, Kubicek CP. Comparative analysis of the *Trichoderma reesei* transcriptome during growth on the cellulase inducing substrates wheat straw and lactose. *Biotechnol Biofuels.* 2013;6:127.

95. Kunitake E, Kawamura A, Tani S, Takenaka S, Ogasawara W, Sumitani J, Kawaguchi T. Effects of *clbR* overexpression on enzyme production in *Aspergillus aculeatus* vary depending on the cellulosic biomass-degrading enzyme species. *Biosci Biotechnol Biochem.* 2015;79:488-95.

96. Sun J, Glass NL. Identification of the CRE-1 cellulolytic regulon in Neurospora crassa. PLoS ONE. 2011;6:e25654.

97. Thieme N, Wu VW, Dietschmann A, Salamov AA, Wang M, Johnson J, Singan VR, Grigoriev IV, Glass NL, Somerville CR, Benz JP. The transcription factor PDR-1 is a multi-functional regulator and key component of pectin deconstruction and catabolism in Neurospora crassa. *Biotechnol Biofuels.* 2017;10:149.

98. Randhawa A, Ogunyewo OA, Eqbal D, Gupta M, Yazdani SS. Disruption of zinc finger DNA binding domain in catabolite repressor Mig1 increases growth rate, hyphal branching, and cellulase expression in hypercellulolytic fungus *Penicillium funiculosum* NCIM1228. *Biotechnol Biofuels.* 2018;11:15.

99. Brown NA, de Gouvea PF, Krohn NG, Savoldi M, Goldman GH. Functional characterisation of the non-essential protein kinases and phosphatases regulating *Aspergillus nidulans* hydrolytic enzyme production. *Biotechnol Biofuels.* 2013;6:91.

100. de Vries RP, Visser J, de Graaff LH. CreA modulates the XlnR-induced expression on xylose of *Aspergillus niger* genes involved in xylan degradation. *Res Microbiol.* 1999;150:281-5.

101. Mach-Aigner AR, Omony J, Jovanovic B, van Boxtel AJ, de Graaff LH. D-Xylose concentration-dependent hydrolase expression profiles and the function of CreA and XlnR in *Aspergillus niger*. *Appl Environ Microbiol.* 2012;78:3145-55.

102. Mach-Aigner AR, Pucher ME, Mach RL. D-Xylose as a repressor or inducer of xylanase expression in *Hypocrea jecorina* (*Trichoderma reesei*). *Appl Environ Microbiol.* 2010;76:1770-6.

103. Foreman PK, Brown D, Dankmeyer L, Dean R, Diener S, Dunn-Coleman NS, Goedegebuur F, Houfek TD, England GJ, Kelley AS, et al. Transcriptional regulation of biomass-degrading enzymes in the filamentous fungus *Trichoderma reesei*. *J Biol Chem.* 2003;278:31988-97.

104. Kubicek CP, Messner R, Gruber F, Mach RL, Kubicek-Pranz EM. The Trichoderma cellulase regulatory puzzle: from the interior life of a secretory fungus. *Enzyme Microb Technol.* 1993;15:90-9.

105. Ilmen M, Saloheimo A, Onnela ML, Penttilä ME. Regulation of cellulase gene expression in the filamentous fungus *Trichoderma reesei*. *Appl Environ Microbiol.* 1997;63:1298-306.

106. Deng H, Gao R, Liao X, Cai Y. CRISPR system in filamentous fungi: current achievements and future directions. *Gene*. 2017;627:212-21.

107. Donohoe PD, Barrangou R, May AP. Advances in industrial biotechnology using CRISPR-Cas systems. *Trends Biotechnol.* 2018;36:134-46.

108. Wakai S, Arazoe T, Ogino C, Kondo A. Future insights in fungal metabolic engineering. *Bioresour Technol.* 2017;245:1314-26.

109. Tamayo EN, Villanueva A, Hasper AA, de Graaff LH, Ramon D, Orejas M. CreA mediates repression of the regulatory gene *xlnR* which controls the production of xylanolytic enzymes in *Aspergillus nidulans*. *Fungal Genet Biol.* 2008;45:984-93.

110. De La Mare M, Guais O, Bonnin E, Weber J, Francois JM. Molecular and biochemical characterization of three GH62 alpha-I-arabinofuranosidases from the soil deuteromycete *Penicillium funiculosum*. *Enzyme Microb Technol.* 2013;53:351-8.

111. Lafond M, Guais O, Maestracci M, Bonnin E, Giardina T. Four GH11 xylanases from the xylanolytic fungus *Talaromyces versatilis* act differently on (arabino)xylans. *Appl Microbiol Biotechnol.* 2014;98:6339-52.

112. Lafond M, Tauzin A, Desseaux V, Bonnin E, el Ajandouz H, Giardina T. GH10 xylanase D from *Penicillium funiculosum*: biochemical studies and xylooligosaccharide production. *Microb Cell Fact.* 2011;10:20.

113. Karboune S, Geraert PA, Kermasha S. Characterization of selected cellulolytic activities of multi-enzymatic complex system from *Penicillium funiculosum*. *J Agric Food Chem.* 2008;56:903-9.

114. Zeilinger S, Mach RL, Schindler M, Herzog P, Kubicek CP. Different inducibility of expression of the two xylanase genes *xyn1* and *xyn2* in *Trichoderma reesei*. *J Biol Chem.* 1996;271:25624-9.

115. Mikus M, Hatvani L, Neuhof T, Komon-Zelazowska M, Dieckmann R, Schwecke T, Druzhinina IS, von Dohren H, Kubicek CP. Differential regulation and posttranslational processing of the class II hydrophobin genes from the biocontrol fungus *Hypocrea atroviridis*. *Appl Environ Microbiol.* 2009;75:3222-9.

116. Rineau F, Lmalem H, Ahren D, Shah F, Johansson T, Coninx L, Ruytinx J, Nguyen H, Grigoriev I, Kuo A, et al. Comparative genomics and expression levels of hydrophobins from eight mycorrhizal genomes. *Mycorrhiza*. 2017;27:383-96.

117. Inukai S, Kock KH, Bulyk ML. Transcription factor-DNA binding: beyond binding site motifs. *Curr Opin Genet Dev.* 2017;43:110-9.

118. Itavo LC, Soares CM, Itavo CC, Dias AM, Petit HV, Leal ES, de Souza AD. Calorimetry, chemical composition and in vitro digestibility of oilseeds. *Food Chem.* 2015;185:219-25.

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

BMC