- Center for Health and Bioresources, AIT Austrian Institute of Technology GmbH, Tulln, Austria
The complex environment of fungi requires a delicate balance between the efforts to acquire nutrition, to reproduce, and to fend off competitors. In Trichoderma reesei, an interrelationship between regulation of enzyme gene expression and secondary metabolism was shown. In this study, we investigated the physiological relevance of the unique YPK1-type kinase USK1 of T. reesei. Usk1 is located in the vicinity of the SOR cluster and is involved in regulation of several genes from this secondary metabolite cluster as well as dihydrotrichotetronine and other secondary metabolites. Moreover, USK1 is required for biosynthesis of normal levels of secondary metabolites in liquid culture. USK1 positively influences cellulase gene regulation, secreted cellulase activity, and biomass formation upon growth in constant darkness on cellulose. Positive effects of USK1 on transcript abundance of the regulator of secondary metabolism, vel1, and the carbon catabolite repressor gene cre1 are in agreement with these functions. In summary, we found that with USK1, T. reesei comprises a unique kinase that adds an additional layer of regulation to the connection of secondary metabolism and enzyme production in fungi.
Introduction
Trichoderma reesei (syn. Hypocrea jecorina) represents one of the most important filamentous fungi for industrial applications (Bischof et al., 2016). T. reesei is a model organism for regulation of plant cell wall degradation due to its efficient machinery for regulation and secretion of carbohydrate active enzymes (Saloheimo and Pakula, 2012; Gupta et al., 2016; Paloheimo et al., 2016). Plant cell wall degrading enzymes are largely co-regulated in T. reesei (Foreman et al., 2003) and their expression is regulated by a network of transcription factors (Benocci et al., 2017), of which XYR1, ACE3, and the carbon catabolite repressor CRE1 are most important. Besides transcription factors, also several signal transduction pathways impact cellulase regulation and especially light as a signal has a profound influence on enzyme expression (Schmoll, 2018).
Trichoderma reesei has a long tradition of application in industry and received generally regarded as safe (GRAS) status (Nevalainen et al., 1994). Recently, the safety of T. reesei as an industrial producer of food and feed enzymes has been re-assessed with a very low level of risk for contamination with harmful toxins (Frisvad et al., 2018). In recent years, secondary metabolism of T. reesei became an interesting field of investigation and also some interconnections with enzyme production became obvious, highlighting the importance of secondary metabolite screening done routinely for industrial production strains.
Sorbicillinoids represent a group of secondary metabolites produced by filamentous fungi, for which anti-inflamatory and antimicrobial functions are reported (Cram, 1948; Maskey et al., 2005) and potential anticancer and anti-HIV applications are discussed (Bringmann et al., 2003; Yao et al., 2015). The biosynthetic genes responsible for production were first elucidated in Penicillium chrysogenum (Salo et al., 2016; Guzman-Chavez et al., 2017). In T. reesei, the production of sorbicillin compounds was investigated at the molecular level. The gene cluster responsible for production of the respective yellow pigments (SOR cluster) is regulated by the two transcription factors YPR1 and YPR2 as well as the carbon catabolite repressor CRE1 and the secondary metabolism regulator LAE1 (Derntl et al., 2016, 2017a; Monroy et al., 2017). The cluster arose in early Hypocreales and lateral gene transfer led to the structure in the T. reesei genome (Druzhinina et al., 2016). The produced sorbicillinoids exert a growth limiting effect on other microbes (Derntl et al., 2017a). Regulation of the SOR cluster is altered in response to the carbon source (Derntl et al., 2016; Dattenböck et al., 2018; Hitzenhammer et al., 2019) and on cellulose a negative feedback cycle is obvious upon growth in light and a positive cycle in darkness (Monroy et al., 2017). Previously, also high level constitutive production of sorbicillinoids was achieved by random mutagenesis of T. reesei RutC30 (Li et al., 2018). However, the regulatory basis for this effect remains to be investigated.
Phosphate residues, attached to and removed from regulatory proteins in a precisely orchestrated manner, represent the currency of signaling processes within the cell. Thereby, biological activity, subcellular localization, half-life, and posttranscriptional modifications as well as interactions with other proteins are regulated (Cohen, 2000; Kosti et al., 2010; Turra et al., 2014). In fungi, diverse physiological functions are known for protein kinases and protein phosphatases including modulation of development, expression of plant cell wall degrading enzymes, secondary metabolism, as well as circadian rhythmicity and light response (Diernfellner and Schafmeier, 2011; Park et al., 2011; Ghosh et al., 2014; de Assis et al., 2015; Winkelströter et al., 2015). Recently, an involvement of different protein phosphatases in many of these mechanisms was shown for T. reesei as well (Rodriguez-Iglesias and Schmoll, 2019). The genome of T. reesei comprises 103 predicted protein kinase genes (Schmoll et al., 2016). Phosphoproteomic analysis of the reaction of T. reesei to inducing conditions provided insights into a complex signaling network, with phosphorylation changes in proteins associated with light-mediated cellulase regulation, carbon sensing, and osmoregulation as well as in the glycolytic pathway (Nguyen et al., 2016).
Here we investigated a kinase with similarity to YPK1 type kinases, which turned out to be unique in Trichoderma spp. With functions in regulation of genes within the SOR cluster, but also in cellulase gene expression, this kinase may represent an interesting evolutionary adaptation of T. reesei to the need of balanced resource distribution between nutrient acquisition and reaction to competitors.
Results
A Unique AGC Family Protein Kinase in the Vicinity of the SOR Cluster
Searching for signaling factors modulating expression of the SOR cluster, we found an AGC family serine/threonine protein kinase represented by the gene model TR_53776 in its vicinity (scaffold 1:478516-479682; 2.3 kb upstream of TR_73618/pks11), which may regulate transcription factor activity or influence enzyme stability by phosphorylation. A recent study on protein phosphorylation in T. reesei (Nguyen et al., 2016) revealed phosphorylation of the transporter within the SOR cluster (TR_43701), which impacts secondary metabolite levels (Monroy et al., 2017). TR_53776 is related to YPK1 type kinases (protein domain cd11651), which are essential for cell growth in yeast and play a role in endocytosis and cell wall integrity. The ypk1 homolog in Neurospora crassa (NCU07280) was found to be essential (Park et al., 2011), but analyses of conditional mutants in Aspergillus nidulans showed functions in growth, hyphal morphogenesis, and delivery of cell membrane and cell wall constituents to the hyphal apex (Colabardini et al., 2013). TR_53776 comprises one potential PEST motif (amino acids 6–24) in its sequence, which indicates regulation of protein stability by phosphorylation (Rechsteiner and Rogers, 1996) as well as eight poor PEST motifs. TargetP 2.0, NetNES 1.1, and WoLF PSORT were used for prediction of localization; however, no targeting sequences to the nucleus or other cellular compartments were detected. Evaluation of available transcriptome data revealed that TR_53776 is regulated in an induction specific manner (Stappler et al., 2017a), but no significant regulation was observed by light, photoreceptors or CRE1 (Tisch et al., 2011b; Tisch and Schmoll, 2013; Monroy et al., 2017).
We further found that the genomic localization of TR_53776 close to the SOR cluster is not syntenic in other fungi such as Aspergillus niger or P. chrysogenum. Accordingly, bidirectional best hit analysis showed that the best blast hits of TR_53776 in these fungi were related to kinases other than TR_53776. Therefore, we performed a phylogenetic analysis of TR_53776 and its most closely related kinases in T. reesei using the best hits of related fungi along with putative homologs from more distantly related fungi. This analysis showed that TR_53776 is not conserved in fungi, not even in closely related Trichoderma spp. like Trichoderma harzianum (Figure 1). Only for Trichoderma atroviride, a homolog was detected. Consequently, we named TR_53776 unique SOR cluster kinase 1 (usk1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/512422/fmicb-11-00974-HTML/image_m/fmicb-11-00974-g001.jpg)
Figure 1. Phylogenetic analysis of usk1 in fungi. The phylogenetic tree was obtained using MEGA4.0 with the maximum parsinomy method. Numbers at branches indicate bootstrap support values. Protein IDs or locus IDs are provided along with species names for sequences originating from Trichoderma reesei, Trichoderma harzianum, Fusarium graminearum, Beauveria bassiana, Magnaporthe grisea, Neurospora crassa, Myceliophthora heterothallica, Penicillium chrysogenum, Aspergillus niger, Aspergillus nidulans, Phlebiopsis gigantea, Phycomyces blakesleeanus, and Postia placenta.
Evaluation of a Relevance of usk1 for Growth
While usk1 is not an ortholog of the well characterized ypk1 kinases, its close relationship still made us evaluate whether there are functional similarities by deleting usk1 in QM6a. In contrast to YPK1 homologs in N. crassa and A. nidulans, T. reesei USK1 is not essential for growth. Strains lacking usk1 are viable and analyses of hyphal extension on solid media with malt extract, cellulose, or glucose did not show any growth defect (Supplementary Figures S1, S2). We also tested the relevance of USK1 for dealing with osmotic stress including ion stress (minimal media with 1 M NaCl or 1 M sorbitol) upon growth in the presence of cellulose or glucose in light and darkness. However, also under these conditions, no significant alterations in growth compared to wild-type was observed (Supplementary Figures S1, S2).
USK1 Positively Impacts Cellulase Gene Expression
Transcript levels of usk1 are upregulated in an induction specific manner (Stappler et al., 2017a). Therefore, we evaluated a relevance of this kinase for regulation of cellulase gene expression. Analysis of transcript levels of cbh1 upon growth on minimal medium with cellulose as carbon source revealed a positive impact of USK1 on cbh1 transcript abundance in darkness (Figure 2A), but not in light (Figure 2B). Accordingly, biomass formation of Δusk1 in liquid media with cellulose as carbon source decreased considerably in darkness (Figure 2C) and specific secreted endo-β-1,4-glucanase activity dropped (Figure 2D). Accordingly, also liberation of reducing sugars from cellulose, representing cellulase activity, decreased (Figure 2E). However, β-glucosidase activity remained largely unchanged in.Δusk1 (Figure 2F). In light only a minor alteration in biomass formation (83 ± 5.4%) was observed between mutant and wildtype. Specific endo-β-1,4-glucanase activity representing cellulase activity remained below detection limits in the deletion strain, as in the wildtype upon growth in constant light.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/512422/fmicb-11-00974-HTML/image_m/fmicb-11-00974-g002.jpg)
Figure 2. Regulatory impact of USK1 on cellulase regulation and on CRE1. Transcript levels of cellobiohydrolase 1 (cbh1) as representative cellulase are shown for growth on microcrystalline cellulose in constant darkness (A; DD) and constant light (B; LL). (C) Biomass formation upon growth on microcrystalline cellulose in constant darkness. (D) Specific endo-1,4 β-D-glucanase activity upon growth on cellulose in darkness. (E) Specific cellulase activity (reducing sugars) upon growth on cellulose in darkness. (F) Specific β-glucosidase activity upon growth on cellulose in darkness. (G,H) Transcript abundance of cre1 upon growth on cellulose in constant darkness (G) or constant light (H). Error bars indicate standard deviations of three biological replicates. Asterisks indicate statistically significant differences compared to wildtype (*p-value < 0.05, **p-value < 0.01, ***p-value < 0.001).
Since CRE1 was previously shown to impact the SOR cluster (Monroy et al., 2017) besides its function as a carbon catabolite repressor (Strauss et al., 1995), we tested an influence of USK1 on transcript abundance of cre1. In darkness, we observed a positive effect of USK1 (about twofold) on transcript levels of cre1 on cellulose (Figure 2G), while no effect was seen in light (Figure 2H).
USK1 Is Involved in Regulation of Secondary Metabolism
Due to its genomic position close to the SOR cluster, we analyzed a potential influence of USK1 on gene regulation of secreted metabolites in general, of those in the SOR cluster and the impact on production of known secondary metabolites produced by T. reesei.
First we tested for a general regulation of secreted metabolite patterns by analysis of secreted compounds during cultivation on minimal medium with cellulose as carbon source by high performance thin layer chromatography (HPTLC) (Hinterdobler et al., 2019). This method provides an informative visual overview on abundance of metabolites in multiple samples and represents a valuable alternative to mass spectrometry, particularly for organisms with only few compounds already known and characterized like T. reesei (Frisvad et al., 2018; Hinterdobler et al., 2019). We found a clear decrease in overall abundance of metabolites upon growth in liquid culture with cellulose as carbon source in darkness in Δusk1 (Figures 3A–D), but no changes in light compared to wildtype under the same conditions (Figures 3E,F). Although biosynthesis of yellow pigments, largely representing sorbicillinoids, is at a low level under these conditions, the positive effect of USK1 on production of individual compounds of this group is still visible (Figures 3C,D). Accordingly, analysis of the supernatant for absorbance at 370 nm, representing part of the sorbicillinoids, is decreased (Figure 3G).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/512422/fmicb-11-00974-HTML/image_m/fmicb-11-00974-g003.jpg)
Figure 3. Impact of USK1 on production of secondary metabolites and their regulators. Wildtype QM6a, parental strain QM6aΔmus53 and Δusk1 were grown on minimal medium in constant light (LL; E,F) or constant darkness (DD; A–D) with cellulose as carbon source. The individual panels show different visualization techniques showing different sets of secondary metabolites of substance classes responsive to the respective technique. (A,E) show visualization at visible light after derivatization with anisaldehyde solution. (B,F) show absorption at 254 nm after development. (C,D) show visualization at visible light after development. For better visibility of the yellow compounds shown in (C) in the saturation was lowered in (D). Analyses were done in biological triplicates and are related to biomass values. (G) Determination of absorbance at 370 nm in the supernatant in relation to the wildtype in constant darkness. Analysis was done in biological triplicates. (H) Quantitative analyses by mass spectrometry with internal standards were performed from strains grown on cellulose in darkness. Abundance of the mentioned metabolites is shown in relation to the wildtype. Values are normalized to biomass.
Quantification of selected secondary metabolites in the supernatant of cellulose grown cultures by mass spectrometry with internal standards revealed decreased levels for the sorbicillin derivative dihydrotrichotetronine (syn. bislongiquinolide or bisorbibutenolide, Shirota et al., 1997; a major product of the SOR cluster, Monroy et al., 2017), alamethicine, orsellinic acid, and paracelsin (Figure 3H). Consequently, USK1 does not only impact sorbicilline compounds associated with the SOR cluster, such as dihydrotrichotetronine, but also other compounds, which is in agreement with HPTLC analysis. Hence, this putative kinase is not specifically targeting sorbicillin production, but presumably impacts regulation of multiple secondary metabolite gene clusters to modulate secondary metabolism more broadly.
In order to gain insight into the regulatory basis of the effect of USK1 on secondary metabolite production, we analyzed its impact on transcript abundance of the general regulator of secondary metabolism, vel1, and the two transcription factor genes associated with the SOR cluster, ypr1 and ypr2. We found that USK1 positively regulates transcript abundance of vel1 in darkness, but slightly negatively in light (Figures 4A,B). For ypr1, we found negative regulation by USK1 in light and darkness (Figures 4C,D), but no effect on ypr2 (Figures 4E,F).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/512422/fmicb-11-00974-HTML/image_m/fmicb-11-00974-g004.jpg)
Figure 4. Influence of USK1 on regulators of secondary metabolism and the SOR cluster. Evaluation of the impact of USK1 on known regulators of secondary metabolism. Transcript abundance of vel1 (A,B), ypr1 (C,D), and ypr2 (E,F) upon growth on cellulose in constant darkness (DD; A,C,E) or constant light (LL; B,D,F). (G–P) Impact of USK1 on SOR cluster genes. Transcript levels of sor1 (TR_73618; G,H), sor2 (TR_73621; I,J), sor5 (TR_73623; K,L), TR_73631 (M,N), and sor4 (TR_43701; O,P) are shown for growth on microcrystalline cellulose in constant darkness (DD; G,I,K,M,O) and constant light (LL; H,J,L,N,P). Error bars indicate standard deviations of three biological replicates. Asterisks indicate statistically significant differences compared to wildtype (*p-value < 0.05, **p-value < 0.01).
Screening for an influence of USK1 on transcript abundance of SOR cluster genes (Figures 4G–P) showed a negative effect USK1 on sor1 (TR_73618), sor5 (TR_73623), and TR_73631 in light and darkness. The gene encoding the transporter SOR4 (TR_43701) is positively regulated in darkness by USK1 (Figure 4O). While USK1 appears to exert a negative effect on transcript levels of several SOR cluster genes, the decreased abundance of secondary metabolites associated with this cluster in the culture filtrate as shown in Figures 3C,D,G may be due to the decreased levels of the transporter TR_43701 (Figure 4O).
Discussion
As one of the most important biotechnological workhorses among filamentous fungi, the physiology of T. reesei and particularly its secondary metabolism are of high interest to industry. It is of utmost importance to ensure the lack of harmful secondary metabolites in food and feed products from fungal fermentations. Therefore, investigation of regulators and the conditions under which they are operative is crucial to understand secondary metabolism in T. reesei.
Protein kinases generally have broad functions in fungi, including regulation of pathogenicity and virulence, stress response, circadian rhythmicity, development, and metabolism (Dickman and Yarden, 1999; Diernfellner and Schafmeier, 2011; Brown et al., 2014; Turra et al., 2014; Manfiolli et al., 2019). The kinomes of fungi share similarity, albeit domain distribution and density of kinomes are suggested to reflect taxonomy (Kosti et al., 2010). In A. nidulans, protein kinases were shown to impact secondary metabolism, particularly the group of mitogen activated protein kinases (De Souza et al., 2013). MAPkinases impact secondary metabolism also in Podospora anserina (Bidard et al., 2012) and N. crassa (Park et al., 2008). USK1 is related to YPK-type protein kinases (particularly to YPK2) which are essential for cell growth and maintenance of cell wall integrity (Chen et al., 1993; Roelants et al., 2002). Saccharomyces cerevisiae YPK2 is subject to phosphorylation by TOR kinases, thereby impacting actin polarization (Kamada et al., 2005). Moreover, YPK2 is implicated in the sphingolipid mediated signaling pathway which targets multiple physiological functions like endocytosis, ubiquitin dependent proteolysis, regulation of nutrient uptake, the cell cycle, translation, and heat stress response (Chung et al., 2001; Cowart and Obeid, 2007). The N. crassa kinase most closely related to USK1 is NCU07280/STK50, which is an essential gene in this fungus (Park et al., 2011) and only shows minor regulation in response to phytosphingosine (Videira et al., 2009). A function in glucose starvation was reported for the fission yeast kinase related to USK1, GAD8 (Saitoh et al., 2015). Also in A. nidulans, the most closely related kinase AN5973 is an essential gene required for polar axis establishment and germling growth (Colabardini et al., 2013; De Souza et al., 2013).
Although USK1 is closely related to these kinases, our analyses revealed that it represents a unique kinase in Trichoderma spp. Accordingly, deletion of usk1 also was not fatal and growth characteristics of Δusk1 were largely normal. Nevertheless, USK1 is highly relevant for expression of cellulolytic enzymes and impacts regulation of genes of the SOR cluster, which is in its genomic vicinity. In Trichoderma parareesei, with the closest homology to T. reesei this localization is conserved, but already in Trichoderma citrinoviride, the locus of usk1 and the SOR cluster is altered as it is in T. atroviride and Trichoderma gamsii. This is in agreement with the SOR cluster being only present in the section Longibrachiatum (Druzhinina et al., 2016). Consequently, the regulatory impact of USK1 on the SOR cluster could have been a beneficial event after acquisition by lateral gene transfer.
Cellulase gene expression is regulated in response to diverse environmental cues, the most important of them being the kind of carbon source and its abundance as well as the impact of light. This regulation is influenced broadly by components of different signal transduction pathways, which are known to be subject to modifications in activity levels by phosphorylation (Bazafkan et al., 2014; Paloheimo et al., 2016; Schmoll, 2018). Transcript levels of usk1 are decreased upon growth in darkness on glucose compared to cellulose in QM6a (Beier, Hinterdobler, and Schmoll, unpublished). Additionally, usk1 showed induction specific regulation (Stappler et al., 2017a). Considering the substantial impact of USK1 on biomass formation and cellulase production upon growth on cellulose, it can be assumed that the major function of USK1 is in nutrient signal transduction and potentially in its coordination with regulation of secondary metabolism. Thereby, the impact of USK1 is not limited to the genes of the cluster and the compounds associated with the SOR cluster. With a relevance of USK1 also to the production of alamethicine and paracelsine, a broader function of USK1 on secondary metabolism is obvious, which is also reflected in considerably decreased levels of secreted metabolites detected in the HPTLC analysis (Figure 3). These findings are in agreement with previous reports, showing that alterations in the levels SOR cluster compounds are often accompanied with modulations of further metabolites associated with other secondary metabolite clusters (Bazafkan et al., 2017; Monroy et al., 2017; Dattenböck et al., 2018). Also for the transcription factor YPR2, which regulates the SOR cluster (Derntl et al., 2016; Monroy et al., 2017), regulation of multiple genes involved in secondary metabolism was shown (Hitzenhammer et al., 2019). For the transcription factor XPP1, which exerts functions in primary and secondary metabolism, a broad impact on both quantity and diversity of secreted metabolites was shown as well (Derntl et al., 2017b). While the actual target(s) of USK1 remain to be elucidated, it can be assumed that in agreement with the results for YPR2 and XPP1, regulation of secondary metabolism is governed by a flat hierarchical network that targets gene regulation broadly rather than specifically, hence establishing a complex network of interdependent regulatory effects. The modulating function of USK1 on the SOR cluster genes and their products as well as on enzyme production moreover add further support for a crosstalk between enzyme expression and secondary metabolism in T. reesei.
Materials and Methods
Strains and Cultivation Conditions
QM6a (ATCC13631) (Martinez et al., 2008) and QM6aΔmus53 (Steiger et al., 2011) was used in this study. For analysis of gene regulation by USK1, strains were cultivated in constant darkness (to exclude interference by circadian rhythmicity) on malt extract agar plates for 14 days to produce conidia for inoculation.
For liquid cultures 108 conidia/L were used for inoculation. Strains grown in Mandels Andreotti minimal medium (Mandels and Andreotti, 1978) with 1% (w/v) cellulose (avicel, SERVA, Alfa Aesar, Karlsruhe, Germany) in constant light (1700 lux, Osram L18W/835 fluorescent bulbs) or in constant darkness at 28°C at 200 r/min. Mycelia and supernatants were harvested after 96 h. Red safety light (darkroom lamp, Philips PF712E, red, E27, 15 W) was applied to avoid random light pulses during harvesting of dark grown cultures.
Hyphal extension was investigated on plates in constant light or constant darkness on the carbon sources mentioned with the experiment and three biological replicates were used. Sexual development was analyzed with strains grown on malt extract media at 22°C in dark-light cycles as described previously (Schmoll, 2013).
Construction of Δusk1
Deletion of usk1 was done in QM6aΔmus53 (deficient in non-homologous end joining) using the streamlined procedure described previously (Schuster et al., 2012) with the hph marker cassette. Transformation was done by protoplasting (Gruber et al., 1990). Successful deletion was confirmed by PCR using primers 53776RTF and 53776RTR. Sequences of all primers used in this study are given in Supplementary Table S1. Copy number determination for integration of the deletion cassette in to the genome of T. reesei was performed as described previously (Tisch and Schmoll, 2013). This analysis showed integration of one single copy of the usk1 deletion cassette into the genome and hence confirmed homologous integration (Supplementary Figure S3).
Nucleic Acid Manipulation and Analysis
Isolation of total RNA, quality control, and RTqPCR were done as described previously (Tisch et al., 2011a; Bazafkan et al., 2017). Sar1 was used as reference gene for RT-qPCR as it was shown to be generally very stable across different conditions (Bazafkan et al., 2017). We considered three biological and three technical replicates for analysis using the CFX Maestro analysis software (Biorad, Hercules, CA, United States). Primers used for analysis are listed in Supplementary Table S1.
Enzyme Assays
Endo-1,4-β-glucanase activity in culture filtrates was determined using azo-CM-cellulose (S-ACMC-L, Megazyme, Wicklow, Ireland) according to the manufacturers instructions.
Cellulase activity as represented by reducing sugars was measured using the pHBAH (p-hydroxybenzoic acid hydrazide) assay as described previously (Mellitzer et al., 2012). Therefore, 1% (w/v) cellulose (Alfa Aesar, Karlsruhe, Germany) served as sole carbon source in the substrate suspension. The absorbance at 410 nm was determined in 96-well flat bottom plates (Greiner Bio One, Kremsmünster, Austria) using a BioTek SynergyTM Mx microplate reader.
Analysis of β-glucosidase was performed using the chromogenic substrate 4-nitrophenyl-β-D-glucopyranosid as described previously (da Silva Delabona et al., 2017).
Biomass for determination of specific enzyme activity was measured in the presence of insoluble cellulose with intracellular protein content as representative criterium as described previously (Schuster et al., 2011; Stappler et al., 2017b). Three biological replicates and technical triplicates were analyzed.
Bioinformatic Analysis
Phylogenetic analysis was performed as described previously (Bazafkan et al., 2015) using the software packages ClustalX1.82 and MEGA4 (Tamura et al., 2007). Therefore, sequences were obtained from JGI Mycocosm. Homology and BLAST searches were performed on the JGI Mycocosm platform1. Screening for PEST domains (Rechsteiner and Rogers, 1996) was done using the online resource ePESTFIND2 with default settings. TargetP 2.0 (Emanuelsson et al., 2007), NetNES 1.1 (la Cour et al., 2004), and WoLF PSORT (Horton et al., 2007) were applied to analyze protein characteristics. For detection of conserved domains, the conserved domain search at NCBI3 (Marchler-Bauer et al., 2017). Statistical significance of values was evaluated by t-test in RStudio (compare_means, ggpubr).
Analysis of Secondary Metabolites Secreted by T. reesei
High performance thin layer chromatography was performed as described previously from culture supernatants of cellulose grown mycelia by extraction with chloroform (Monroy et al., 2017; Hinterdobler et al., 2019). Measurements are related to biomass abundance for every sample.
Abundance of selected secondary metabolites was analyzed by mass spectrometry with internal standards as described previously (Malachova et al., 2014; Monroy et al., 2017).
Data Availability Statement
All datasets relevant for this study are included in the article/Supplementary Material.
Author Contributions
SB contributed to analysis of growth and cellulase production. WH contributed to analysis of secondary metabolite production. AM participated in secondary metabolite analysis and cultivation. HB constructed the recombinant strain. SB and WH participated in drafting the manuscript. MS conceived of the study, participated in analysis and interpretation of results, and wrote the final version of the manuscript. All authors read and agreed to publication of the final manuscript.
Funding
Work of SB, HB, and AM was supported by the Austrian Science Fund (FWF; projects P26935; P24350 to MS). Work of WH was supported by the NFB (Science Fund of Lower Austria; Grant LC16-04 to MS). The funding agencies had no role in study design or publication.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We want to thank Julia Scholda for technical assistance with HPTLC analysis and Marlene Stiegler, Michael Sulyok, and Stefan Böhmdorfer for technical assistance and the possibility to access analytic equipment.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.00974/full#supplementary-material
Footnotes
- ^ https://mycocosm.jgi.doe.gov/mycocosm/home
- ^ http://emboss.bioinformatics.nl/cgi-bin/emboss/epestfind
- ^ https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi?
References
Bazafkan, H., Beier, S., Stappler, E., Böhmdorfer, S., Oberlerchner, J. T., Sulyok, M., et al. (2017). SUB1 has photoreceptor dependent and independent functions in sexual development and secondary metabolism in Trichoderma reesei. Mol. Microbiol. 106, 742–759. doi: 10.1111/mmi.13842
Bazafkan, H., Dattenböck, C., Böhmdorfer, S., Tisch, D., Stappler, E., and Schmoll, M. (2015). Mating type dependent partner sensing as mediated by VEL1 in Trichoderma reesei. Mol. Microbiol. 96, 1103–1118. doi: 10.1111/mmi.12993
Bazafkan, H., Tisch, D., and Schmoll, M. (2014). “Regulation of glycoside hydrolase expression in Trichoderma,” in Biotechnology and Biology of Trichoderma, eds V. K. Gupta, M. Schmoll, A. Herrera-Estrella, I. Druzhinina, and M. G. Tuohy (Oxford: Elsevier), 291–307.
Benocci, T., Aguilar-Pontes, M. V., Zhou, M., Seiboth, B., and de Vries, R. P. (2017). Regulators of plant biomass degradation in ascomycetous fungi. Biotechnol. Biofuels 10:152. doi: 10.1186/s13068-017-0841-x
Bidard, F., Coppin, E., and Silar, P. (2012). The transcriptional response to the inactivation of the PaMpk1 and PaMpk2 MAP kinase pathways in Podospora anserina. Fungal Genet. Biol. 49, 643–652. doi: 10.1016/j.fgb.2012.06.002
Bischof, R. H., Ramoni, J., and Seiboth, B. (2016). Cellulases and beyond: the first 70 years of the enzyme producer Trichoderma reesei. Microb. Cell Fact. 15:106. doi: 10.1186/s12934-016-0507-6
Bringmann, G., Lang, G., Muhlbacher, J., Schaumann, K., Steffens, S., Rytik, P. G., et al. (2003). Sorbicillactone A: a structurally unprecedented bioactive novel-type alkaloid from a sponge-derived fungus. Prog. Mol. Subcell Biol. 37, 231–253. doi: 10.1007/978-3-642-55519-0_9
Brown, N. A., Ries, L. N. A., and Goldman, G. H. (2014). How nutritional status signalling coordinates metabolism and lignocellulolytic enzyme secretion. Fungal Genet. Biol. 72, 48–63. doi: 10.1016/j.fgb.2014.06.012
Chen, P., Lee, K. S., and Levin, D. E. (1993). A pair of putative protein kinase genes (YPK1 and YPK2) is required for cell growth in Saccharomyces cerevisiae. Mol. Gen. Genet. 236, 443–447. doi: 10.1007/bf00277146
Chung, N., Mao, C., Heitman, J., Hannun, Y. A., and Obeid, L. M. (2001). Phytosphingosine as a specific inhibitor of growth and nutrient import in Saccharomyces cerevisiae. J. Biol. Chem. 276, 35614–35621. doi: 10.1074/jbc.M105653200
Cohen, P. (2000). The regulation of protein function by multisite phosphorylation – a 25 year update. Trends Biochem. Sci. 25, 596–601. doi: 10.1016/s0968-0004(00)01712-6
Colabardini, A. C., Brown, N. A., Savoldi, M., Goldman, M. H., and Goldman, G. H. (2013). Functional characterization of Aspergillus nidulans ypkA, a homologue of the mammalian kinase SGK. PLoS ONE 8:e57630. doi: 10.1371/journal.pone.0057630
Cowart, L. A., and Obeid, L. M. (2007). Yeast sphingolipids: recent developments in understanding biosynthesis, regulation, and function. Biochim. Biophys. Acta 1771, 421–431. doi: 10.1016/j.bbalip.2006.08.005
Cram, D. J. (1948). Mold metabolites; the structure of sorbicillin, a pigment produced by the mold Penicillium notatum. J. Am. Chem. Soc. 70, 4240–4243. doi: 10.1021/ja01192a077
da Silva, Delabona, P., Rodrigues, G. N., Zubieta, M. P., Ramoni, J., Codima, C. A., et al. (2017). The relation between xyr1 overexpression in Trichoderma harzianum and sugarcane bagasse saccharification performance. J. Biotechnol. 246, 24–32. doi: 10.1016/j.jbiotec.2017.02.002
Dattenböck, C., Tisch, D., Schuster, A., Monroy, A. A., Hinterdobler, W., and Schmoll, M. (2018). Gene regulation associated with sexual development and female fertility in different isolates of Trichoderma reesei. Fungal Biol. Biotechnol. 5:9. doi: 10.1186/s40694-018-0055-4
de Assis, L. J., Ries, L. N., Savoldi, M., Dinamarco, T. M., Goldman, G. H., and Brown, N. A. (2015). Multiple phosphatases regulate carbon source-dependent germination and primary metabolism in Aspergillus nidulans. G3 (Bethesda) 5, 857–872. doi: 10.1534/g3.115.016667
De Souza, C. P., Hashmi, S. B., Osmani, A. H., Andrews, P., Ringelberg, C. S., Dunlap, J. C., et al. (2013). Functional analysis of the Aspergillus nidulans kinome. PLoS ONE 8:e58008. doi: 10.1371/journal.pone.0058008
Derntl, C., Guzman-Chavez, F., Mello-de-Sousa, T. M., Busse, H. J., Driessen, A. J. M., Mach, R. L., et al. (2017a). In Vivo study of the sorbicillinoid gene cluster in Trichoderma reesei. Front. Microbiol. 8:2037. doi: 10.3389/fmicb.2017.02037
Derntl, C., Kluger, B., Bueschl, C., Schuhmacher, R., Mach, R. L., and Mach-Aigner, A. R. (2017b). Transcription factor Xpp1 is a switch between primary and secondary fungal metabolism. Proc. Natl. Acad. Sci. U.S.A. 114, E560–E569. doi: 10.1073/pnas.1609348114
Derntl, C., Rassinger, A., Srebotnik, E., Mach, R. L., and Mach-Aigner, A. R. (2016). Identification of the main regulator responsible for synthesis of the typical yellow pigment produced by Trichoderma reesei. Appl. Environ. Microbiol. 82, 6247–6257. doi: 10.1128/AEM.01408-16
Dickman, M. B., and Yarden, O. (1999). Serine/threonine protein kinases and phosphatases in filamentious fungi. Fungal Genet. Biol. 26, 99–117. doi: 10.1006/fgbi.1999.1118
Diernfellner, A. C., and Schafmeier, T. (2011). Phosphorylations: making the Neurospora crassa circadian clock tick. FEBS Lett. 585, 1461–1466. doi: 10.1016/j.febslet.2011.03.049
Druzhinina, I. S., Kubicek, E. M., and Kubicek, C. P. (2016). Several steps of lateral gene transfer followed by events of ‘birth-and-death’ evolution shaped a fungal sorbicillinoid biosynthetic gene cluster. BMC Evol. Biol. 16:269. doi: 10.1186/s12862-016-0834-6
Emanuelsson, O., Brunak, S., von Heijne, G., and Nielsen, H. (2007). Locating proteins in the cell using TargetP, SignalP and related tools. Nat. Protoc. 2, 953–971. doi: 10.1038/nprot.2007.131
Foreman, P. K., Brown, D., Dankmeyer, L., Dean, R., Diener, S., Dunn-Coleman, N. S., et al. (2003). Transcriptional regulation of biomass-degrading enzymes in the filamentous fungus Trichoderma reesei. J. Biol. Chem. 278, 31988–31997. doi: 10.1074/jbc.M304750200
Frisvad, J. C., Moller, L. L. H., Larsen, T. O., Kumar, R., and Arnau, J. (2018). Safety of the fungal workhorses of industrial biotechnology: update on the mycotoxin and secondary metabolite potential of Aspergillus niger, Aspergillus oryzae, and Trichoderma reesei. Appl. Microbiol. Biotechnol. 102, 9481–9515. doi: 10.1007/s00253-018-9354-1
Ghosh, A., Servin, J. A., Park, G., and Borkovich, K. A. (2014). Global analysis of serine/threonine and tyrosine protein phosphatase catalytic subunit genes in Neurospora crassa reveals interplay between phosphatases and the p38 mitogen-activated protein kinase. G3 (Bethesda) 4, 349–365. doi: 10.1534/g3.113.008813
Gruber, F., Visser, J., Kubicek, C. P., and de Graaff, L. H. (1990). The development of a heterologous transformation system for the cellulolytic fungus Trichoderma reesei based on a pyrG -negative mutant strain. Curr. Genet. 18, 71–76. doi: 10.1007/bf00321118
Gupta, V. K., Steindorff, A. S., de Paula, R. G., Silva-Rocha, R., Mach-Aigner, A. R., Mach, R. L., et al. (2016). The post-genomic era of Trichoderma reesei: what’s next? Trends Biotechnol. 34, 970–982. doi: 10.1016/j.tibtech.2016.06.003
Guzman-Chavez, F., Salo, O., Nygard, Y., Lankhorst, P. P., Bovenberg, R. A. L., and Driessen, A. J. M. (2017). Mechanism and regulation of sorbicillin biosynthesis by Penicillium chrysogenum. Microb. Biotechnol. 10, 958–968. doi: 10.1111/1751-7915.12736
Hinterdobler, W., Schuster, A., Tisch, D., Ozkan, E., Bazafkan, H., Schinnerl, J., et al. (2019). The role of PKAc1 in gene regulation and trichodimerol production in Trichoderma reesei. Fungal Biol. Biotechnol. 6:12. doi: 10.1186/s40694-019-0075-8
Hitzenhammer, E., Büschl, C., Sulyok, M., Schuhmacher, R., Kluger, B., Wischnitzki, E., et al. (2019). YPR2 is a regulator of light modulated carbon and secondary metabolism in Trichoderma reesei. BMC Genomics 20:211. doi: 10.1186/s12864-019-5574-8
Horton, P., Park, K. J., Obayashi, T., Fujita, N., Harada, H., Adams-Collier, C. J., et al. (2007). WoLF PSORT: protein localization predictor. Nucleic Acids Res. 35, W585–W587. doi: 10.1093/nar/gkm259
Kamada, Y., Fujioka, Y., Suzuki, N. N., Inagaki, F., Wullschleger, S., Loewith, R., et al. (2005). Tor2 directly phosphorylates the AGC kinase Ypk2 to regulate actin polarization. Mol. Cell Biol. 25, 7239–7248. doi: 10.1128/MCB.25.16.7239-7248.2005
Kosti, I., Mandel-Gutfreund, Y., Glaser, F., and Horwitz, B. A. (2010). Comparative analysis of fungal protein kinases and associated domains. BMC Genomics 11:133. doi: 10.1186/1471-2164-11-133
la Cour, T., Kiemer, L., Molgaard, A., Gupta, R., Skriver, K., and Brunak, S. (2004). Analysis and prediction of leucine-rich nuclear export signals. Protein Eng. Des. Sel. 17, 527–536. doi: 10.1093/protein/gzh062
Li, C., Lin, F., Sun, W., Yuan, S., Zhou, Z., Wu, F. G., et al. (2018). Constitutive hyperproduction of sorbicillinoids in Trichoderma reesei ZC121. Biotechnol. Biofuels 11:291. doi: 10.1186/s13068-018-1296-4
Malachova, A., Sulyok, M., Beltran, E., Berthiller, F., and Krska, R. (2014). Optimization and validation of a quantitative liquid chromatography-tandem mass spectrometric method covering 295 bacterial and fungal metabolites including all regulated mycotoxins in four model food matrices. J. Chromatogr. A 1362, 145–156. doi: 10.1016/j.chroma.2014.08.037
Mandels, M., and Andreotti, R. (1978). Problems and challenges in the cellulose to cellulase fermentation. Proc. Biochem. 13, 6–13.
Manfiolli, A. O., Siqueira, F. S., Dos Reis, T. F., Van Dijck, P., Schrevens, S., Hoefgen, S., et al. (2019). Mitogen-activated protein kinase cross-talk interaction modulates the mroduction of Melanins in Aspergillus fumigatus. mBio 10, e215–e219. doi: 10.1128/mBio.00215-19
Marchler-Bauer, A., Bo, Y., Han, L., He, J., Lanczycki, C. J., Lu, S., et al. (2017). CDD/SPARCLE: functional classification of proteins via subfamily domain architectures. Nucleic Acids Res. 45, D200–D203. doi: 10.1093/nar/gkw1129
Martinez, D., Berka, R. M., Henrissat, B., Saloheimo, M., Arvas, M., Baker, S. E., et al. (2008). Genome sequencing and analysis of the biomass-degrading fungus Trichoderma reesei (syn. Hypocrea jecorina). Nat. Biotechnol. 26, 553–560. doi: 10.1038/nbt1403
Maskey, R. P., Grun-Wollny, I., and Laatsch, H. (2005). Sorbicillin analogues and related dimeric compounds from Penicillium notatum. J. Nat. Prod. 68, 865–870. doi: 10.1021/np040137t
Mellitzer, A., Glieder, A., Weis, R., Reisinger, C., and Flicker, K. (2012). Sensitive high-throughput screening for the detection of reducing sugars. Biotechnol. J. 7, 155–162. doi: 10.1002/biot.201100001
Monroy, A. A., Stappler, E., Schuster, A., Sulyok, M., and Schmoll, M. (2017). A CRE1- regulated cluster is responsible for light dependent production of dihydrotrichotetronin in Trichoderma reesei. PLoS ONE 12:e0182530. doi: 10.1371/journal.pone.0182530
Nevalainen, H., Suominen, P., and Taimisto, K. (1994). On the safety of Trichoderma reesei. J. Biotechnol. 37, 193–200.
Nguyen, E. V., Imanishi, S. Y., Haapaniemi, P., Yadav, A., Saloheimo, M., Corthals, G. L., et al. (2016). Quantitative site-specific phosphoproteomics of Trichoderma reesei signaling pathways upon induction of hydrolytic enzyme production. J. Proteome Res. 15, 457–467. doi: 10.1021/acs.jproteome.5b00796
Paloheimo, M., Haarmann, T., Mäkinen, S., and Vehmaanperä, J. (2016). “Production of industrial enzymes in Trichoderma reesei,” in Gene Expression Systems in Fungi: Advancements and Applications, eds M. Schmoll and C. Dattenböck (Heidelberg: Springer International), 23–58.
Park, G., Pan, S., and Borkovich, K. A. (2008). Mitogen-activated protein kinase cascade required for regulation of development and secondary metabolism in Neurospora crassa. Eukaryot. Cell 7, 2113–2122. doi: 10.1128/EC.00466-07
Park, G., Servin, J. A., Turner, G. E., Altamirano, L., Colot, H. V., Collopy, P., et al. (2011). Global analysis of serine-threonine protein kinase genes in Neurospora crassa. Eukaryot. Cell 10, 1553–1564. doi: 10.1128/EC.0514011
Rechsteiner, M., and Rogers, S. W. (1996). PEST sequences and regulation by proteolysis. Trends Biochem. Sci. 21, 267–271.
Rodriguez-Iglesias, A., and Schmoll, M. (2019). Protein phosphatases regulate growth, development, cellulases and secondary metabolism in Trichoderma reesei. Sci. Rep. 9:10995. doi: 10.1038/s41598-019-47421-z
Roelants, F. M., Torrance, P. D., Bezman, N., and Thorner, J. (2002). Pkh1 and Pkh2 differentially phosphorylate and activate Ypk1 and Ykr2 and define protein kinase modules required for maintenance of cell wall integrity. Mol. Biol. Cell 13, 3005–3028. doi: 10.1091/mbc.e02-04-0201
Saitoh, S., Mori, A., Uehara, L., Masuda, F., Soejima, S., and Yanagida, M. (2015). Mechanisms of expression and translocation of major fission yeast glucose transporters regulated by CaMKK/phosphatases, nuclear shuttling, and TOR. Mol. Biol. Cell 26, 373–386. doi: 10.1091/mbc.E14-11-1503
Salo, O., Guzman-Chavez, F., Ries, M. I., Lankhorst, P. P., Bovenberg, R. A., Vreeken, R. J., et al. (2016). Identification of a polyketide synthase involved in sorbicillin biosynthesis by Penicillium chrysogenum. Appl. Environ. Microbiol. 82, 3971–3978. doi: 10.1128/AEM.0035016
Saloheimo, M., and Pakula, T. M. (2012). The cargo and the transport system: secreted proteins and protein secretion in Trichoderma reesei (Hypocrea jecorina). Microbiology 158(Pt 1), 46–57. doi: 10.1099/mic.0.053132-0
Schmoll, M. (2013). “Sexual development in Trichoderma – scrutinizing the aspired phenomenon,” in Trichoderma: Biology and Applications, eds P. K. Mukherjee, B. A. Horwitz, U. S. Singh, M. Mukherjee, and M. Schmoll (Wallingford: CAB International), 67–86.
Schmoll, M. (2018). Regulation of plant cell wall degradation by light in Trichoderma. Fungal Biol. Biotechnol. 5:10. doi: 10.1186/s40694-018-0052-7
Schmoll, M., Dattenböck, C., Carreras-Villasenor, N., Mendoza-Mendoza, A., Tisch, D., Aleman, M. I., et al. (2016). The genomes of three uneven siblings: footprints of the lifestyles of three Trichoderma species. Microbiol. Mol. Biol. Rev. 80, 205–327. doi: 10.1128/MMBR.00040-15
Schuster, A., Bruno, K. S., Collett, J. R., Baker, S. E., Seiboth, B., Kubicek, C. P., et al. (2012). A versatile toolkit for high throughput functional genomics with Trichoderma reesei. Biotechnol. Biofuels 5:1. doi: 10.1186/1754-6834-5-1
Schuster, A., Kubicek, C. P., and Schmoll, M. (2011). Dehydrogenase GRD1 represents a novel component of the cellulase regulon in Trichoderma reesei (Hypocrea jecorina). Appl. Environ. Microbiol. 77, 4553–4563. doi: 10.1128/AEM.00513-11
Shirota, O., Pathak, V., Chowdhury, F. H., Setsuko, S., Takatori, K., and Satake, M. (1997). Structural elucidation of trichotetronines: polyketides possessing a bicyclo [2.2.2] octane skeleton with a tetronic acid moiety isolated from Trichoderma spp. J. Chem. Soc. Perkin Trans. 1, 2961–2964.
Stappler, E., Dattenböck, C., Tisch, D., and Schmoll, M. (2017a). Analysis of light- and carbon-specific transcriptomes implicates a class of G-protein-coupled receptors in cellulose sensing. mSphere 2:e89-17. doi: 10.1128/mSphere.00089-17
Stappler, E., Walton, J. D., and Schmoll, M. (2017b). Abundance of secreted proteins of Trichoderma reesei is regulated by light of different intensities. Front. Microbiol. 8:2586. doi: 10.3389/fmicb.2017.02586
Steiger, M. G., Vitikainen, M., Uskonen, P., Brunner, K., Adam, G., Pakula, T., et al. (2011). Transformation system for Hypocrea jecorina (Trichoderma reesei) that favors homologous integration and employs reusable bidirectionally selectable markers. Appl. Environ. Microbiol. 77, 114–121. doi: 10.1128/AEM.02100-10
Strauss, J., Mach, R. L., Zeilinger, S., Hartler, G., Stoffler, G., Wolschek, M., et al. (1995). Cre1, the carbon catabolite repressor protein from Trichoderma reesei. FEBS Lett. 376, 103–107. doi: 10.1016/0014-5793(95)01255-5
Tamura, K., Dudley, J., Nei, M., and Kumar, S. (2007). MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol. Biol. Evol. 24, 1596–1599. doi: 10.1093/molbev/msm092
Tisch, D., Kubicek, C. P., and Schmoll, M. (2011a). New insights into the mechanism of light modulated signaling by heterotrimeric G-proteins: ENVOY acts on gna1 and gna3 and adjusts cAMP levels in Trichoderma reesei (Hypocrea jecorina). Fungal Genet. Biol. 48, 631–640. doi: 10.1016/j.fgb.2010.12.009
Tisch, D., Kubicek, C. P., and Schmoll, M. (2011b). The phosducin-like protein PhLP1 impacts regulation of glycoside hydrolases and light response in Trichoderma reesei. BMC Genomics 12:613. doi: 10.1186/1471-2164-12-613
Tisch, D., and Schmoll, M. (2013). Targets of light signalling in Trichoderma reesei. BMC Genomics 14:657. doi: 10.1186/1471-2164-14-657
Turra, D., Segorbe, D., and Di Pietro, A. (2014). Protein kinases in plant-pathogenic fungi: conserved regulators of infection. Annu. Rev. Phytopathol. 52, 267–288. doi: 10.1146/annurev-phyto-102313-050143
Videira, A., Kasuga, T., Tian, C., Lemos, C., Castro, A., and Glass, N. L. (2009). Transcriptional analysis of the response of Neurospora crassa to phytosphingosine reveals links to mitochondrial function. Microbiology 155(Pt 9), 3134–3141. doi: 10.1099/mic.0.029710-0
Winkelströter, L. K., Dolan, S. K., Fernanda Dos, Reis, T., Bom, V. L., and Alves et al. (2015). Systematic global analysis of genes encoding protein phosphatases in Aspergillus fumigatus. G3 (Bethesda) 5, 1525–1539. doi: 10.1534/g3.115.016766
Keywords: Trichoderma reesei, Hypocrea jecorina, sorbicillin, dihydrotrichotetronin, secondary metabolism, plant cell wall degradation, signal transduction
Citation: Beier S, Hinterdobler W, Monroy AA, Bazafkan H and Schmoll M (2020) The Kinase USK1 Regulates Cellulase Gene Expression and Secondary Metabolite Biosynthesis in Trichoderma reesei. Front. Microbiol. 11:974. doi: 10.3389/fmicb.2020.00974
Received: 15 November 2019; Accepted: 22 April 2020;
Published: 20 May 2020.
Edited by:
Guodong Liu, Shandong University, ChinaCopyright © 2020 Beier, Hinterdobler, Monroy, Bazafkan and Schmoll. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Monika Schmoll, monika.schmoll@ait.ac.at