# Fungal X-Intrinsic Protein Aquaporin from Trichoderma atroviride: Structural and Functional Considerations

Maroua Ben Amira, Mohamed Faize, Magnus Karlsson, Mukesh Dubey, Magdalena Frąc, Jacek Panek, Boris Fumanal, Aurélie Gousset-Dupont, Jean-Louis Julien, Hatem Chaar, et al.

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

## To cite this version:

Maroua Ben Amira, Mohamed Faize, Magnus Karlsson, Mukesh Dubey, Magdalena Frąc, et al. Fungal X-Intrinsic Protein Aquaporin from Trichoderma atroviride: Structural and Functional Considerations. Biomolecules, 2021, 11 (2), pp.338. 10.3390/biom11020338. hal-03153574

HAL Id: hal-03153574

https://hal.inrae.fr/hal-03153574

Submitted on 26 Feb 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished 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
Article

# Fungal X-Intrinsic Protein Aquaporin from *Trichoderma atroviride*: Structural and Functional Considerations

Maroua BEN AMIRA <fnref n="1" />,<fnref n="2" />, Mohamed FAIZE <fnref n="3" />, Magnus KARLSSON <fnref n="4" />, Mukesh DUBEY <fnref n="4" />, Magdalena FRĄC <fnref n="5" />, Jacek PANEK <fnref n="5" />, Boris FUMANAL <fnref n="1" />, Aurélie GOUSSET-DUPONT <fnref n="1" />, Jean-Louis JULIEN <fnref n="1" />, Hatem CHAAR <fnref n="6" />, Daniel AUGUIN <fnref n="7" />, Robin MOM <fnref n="1" />, Philippe LABEL <fnref n="1" /> and Jean-Stéphane VENISSE <fnref n="1" />,*

**Citation:** Ben Amira, M; Faize, M.; Karlsson, M.; Dubey, M.; Frąc, M.; Panek, J.; Fumanal, B.; Gousset-Dupont, A.; Julien, J.-L.; Chaar, H.; et al. Fungal X-Intrinsic Protein Aquaporin from *Trichoderma atroviride*: Structural and Functional Considerations. *Biomolecules* **2021**, *11*, 338. https://doi.org/10.3390/biom11020338

Received: 11 January 2021
Accepted: 18 February 2021
Published: 23 February 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

**Abstract:** The major intrinsic protein (MIP) superfamily is a key part of the fungal transmembrane transport network. It facilitates the transport of water and low molecular weight solutes across biomembranes. The fungal uncharacterized X-Intrinsic Protein (XIP) subfamily includes the full protein diversity of MIP. Their biological functions still remain fully hypothetical. The aim of this study is still to deepen the diversity and the structure of the XIP subfamily in light of the MIP counterparts—the aquaporins (AQPs) and aquaglyceroporins (AQGPs)—and to describe for the first time their function in the development, biomass accumulation, and mycoparasitic aptitudes of the fungal bioagent *Trichoderma atroviride*. The fungus-XIP clade, with one member (*TriatXIP*), is one of the three clades of MIPs that make up the diversity of *T. atroviride* MIPs, along with the AQPs (three members) and the AQGPs (three members). *TriatXIP* resembles those of strict aquaporins, predicting water diffusion and possibly other small polar solutes due to particularly wider ar/R constriction with a Lysine substitution at the LE2 position. The XIP loss of function in *ΔTriatXIP* mutants slightly delays biomass accumulation but does not impact mycoparasitic activities. *ΔTriatMIP* forms colonies similar to wild type; however, the hyphae are slightly thinner and colonies produce rare chlamydospores in PDA and specific media, most of which are relatively small and exhibit abnormal morphologies. To better understand the molecular causes of these deviant phenotypes, a wide-metabolic survey of the *ΔTriatXIPs* demonstrates that the delayed growth kinetic, correlated to a decrease in respiration rate, is caused by perturbations in the pentose phosphate pathway. Furthermore, the null expression of the *XIP* gene strongly impacts the expression of four expressed *MIP*-encoding genes of *T. atroviride*, a plausible compensating effect which safeguards the physiological integrity and life cycle of the fungus. This paper offers an overview of the fungal XIP family in the biocontrol agent *T. atroviride* which will be useful for further functional analysis of this particular MIP subfamily in vegetative growth and the environmental stress response in fungi. Ultimately, these findings have implications for the ecophysiology of *Trichoderma* spp. in natural, agronomic, and industrial systems.

**Keywords:** aquaporin; uncharacterized X-Intrinsic proteins; *Trichoderma atroviride*; 3D modeling; chlamydospores; pentose phosphate pathway; stress responses

<footnote n="1">¹ Université Clermont Auvergne, INRAE, PIAF, 63000 Clermont-Ferrand, France; marouabenamira@gmail.com (M.B.A.); j-louis.julien@uca.fr (J.-L.J.); boris.fumanal@uca.fr (B.F.); aurelie.gousset@uca.fr (A.G.-D.); Robin.mom@uca.fr (R.M.); philippe.label@inrae.fr (P.L.); j-stephane.venisse@uca.fr (J.S.V.)</footnote>

<footnote n="2">² Faculté des Sciences de Bizerte, Zarzouna 7021, Tunisia</footnote>

<footnote n="3">³ Laboratory of Plant Biotechnology, Ecology and Ecosystem Valorization, Faculty of Sciences, University Chouaib Doukkali, 24000 El Jadida, Morocco; faizemohamed@yahoo.fr</footnote>

<footnote n="4">⁴ Department of Forest Mycology and Plant Pathology, Swedish University of Agricultural Sciences, Uppsala BioCenter, SE-75007, Uppsala, Sweden; Magnus.Karlsson@slu.se (M.K.); Mukesh.Dubey@slu.se (M.D.)</footnote>

<footnote n="5">5 Institute of Agrophysics, Polish Academy of Sciences, Doświadczalna 4, 20-290 Lublin, Poland; m.frac@ipan.lublin.pl (M.F.); j.panek@ipan.lublin.pl (J.P.)</footnote>

<footnote n="6">⁶ Crop Improvement Laboratory, National Institute of Agronomy of Tunisia (INRAT), Ariana 2049, Tunisia; chaarh@yahoo.com</footnote>

<footnote n="7">7 Université d'Orléans, Laboratoire de Biologie des Ligneux et des Grandes Cultures, UPRES EA 1207, INRA-USC1328 Orléans, France; auguin@univ-orleans.fr</footnote>

<footnote>* Correspondence: j-stephane.venisse@uca.fr</footnote>
# 1. Introduction

The *Trichoderma* species are filamentous ascomycete fungi belonging to the order Hypocreales in the class Sordariomycetes. The *Trichoderma* genus is predominantly composed of green-spored and filamentous soil-dwelling fungi found worldwide, which is, above all, cosmopolitan and ubiquitous in the environment. Like any chemo-heterotrophic organisms, they obtain their reduced carbon source from other organisms whether they are living or dead. With the exception of the pathogenic and biotrophic behavior of a few specific *Trichoderma* strains [1–4], *Trichoderma* strains are mainly saprophytes and necrotrophic mycoparasites. Even more promising, and notably with plants, a few strains are part of symbiotic and/or commensal relationships where they establish a robust and long-lasting colonization of organ surfaces with possible asymptomatic penetrations into the epidermis and a few cells below this level [5,6]. As a result, by acting as plant symbionts (endophytism), specific *Trichoderma* strains have been extensively developed as biocontrol agents against different plant diseases due to their ability to successfully antagonize plant pathogenic fungi or bacteria [7–9]. Several indirect and direct coordinated strategies play substantial roles in their ability to compete with other microorganisms [9,10]. In addition, root colonization by *Trichoderma* spp. may mitigate the detrimental effects of abiotic stresses such as salinity and drought by enhancing global plant resistance, uptake and use of nutrients, root growth and development, as well as crop productivity [11,12].

The strategies of *Trichoderma* spp. deployed against antagonized phytopathogens are complex. Our understanding of the mechanisms of biocontrol and the *Trichoderma*-plant-pathogen interaction is incomplete, especially concerning nutrition. *Trichoderma* spp. display a remarkable nutritional versatility. They succeed in various heterotrophic and biotrophic interactions with a very wide range of hosts, being able to exhibit various adapted polyphagical lifestyles by decomposing dead organic matter [13]. They feed by extracellular digestion, secreting into the extracellular environment a cocktail of digestive enzymes that break down biopolymers (cellulose, hemi-cellulose, chitin, lignin, lipids, proteins, etc.) [14]. The resulting external soluble products are subsequently absorbed into the fungal cells by transmembrane diffusions common to all fungi. This route requires an arsenal of integral membrane proteins, channels and transporters, such as a large array of sugar transporters (major facilitator superfamily (MFS)) [15], nitrogen transporters (di/tripeptide transporters (PTR)) [16], and ATP-binding cassette (ABC) transporters [17]. However, most potential systems for membrane transport proteins encoded by filamentous fungi have not yet been characterized, including the major intrinsic protein (MIP) superfamily, also called aquaporins (AQPs) [18]. These channels facilitate the transport of water and a variety of small molecules across biological membranes, the plasmalemic and the endomembrane system, of unicellular and multicellular organisms [19]. In filamentous fungi, preliminary evidence shows that MIPs are involved in the regulation of hyphal growth, sclerotia formation, conidiation, spore germination, virulence, secondary metabolisms, and the transport of various molecules [20–25]. In this respect, the impact of transcriptional and protein regulation of MIPs on the fungus' life cycle remains a fascinating area that deserves further exploration.

The MIP superfamily is present in all kingdoms from Prokaryotes (archaea and eu-bacteria) to Eukaryotes (fungi, insects, animals, and plants) [19]. Most of these clades share two subfamilies: the aquaporins (AQPs, *lato sensu*), and the aquaglyceroporins (AGPs or AQGPs), which distinguish themselves by their primary sequences and their ability to transport various solutes. The aquaporins *lato sensu* encompass both the “orthodox” or “classical” aquaporins (AQPs, *stricto sensu*), which are mainly water channels. The aquaglyceroporins additionally transport small organic compounds such as urea and glycerol due to differential hydrophobic and pore size features. Fungi possess both of these two MIP prototype groups in which four distinct sub-groups can be distinguished based on their primary sequences: the orthodox fungal water channels, the facultative fungal aquaporins, the fungal aquaglyceroporins, and the fungal uncategorized (X)-intrinsic proteins (XIPs) [26]. This fourth group, XIP, resembles orthodox aquaporins. It is
predominantly shared between certain clades of fungi and Viridiplantae, and marginally present in some protozoa [26–31]. So far, XIP has been found to be the only MIP subgroup common to both plants and mycetes, implying a particular evolutionary relationship between these two unrelated phylogenetic lineages that needs to be resolved. In addition, this suggests specific roles in solute transport and/or protein regulation of these channels which deserve better understanding.

Like any other MIP, the XIP protein protomer structure is highly conserved and resembles an “hourglass model” [29,30,32]. In plants, several *XIP* copies per genome constitute the XIP subfamily. XIP proteins are located at the plasma membrane of the epidermal and parenchyma cells and at the endomembrane system. They transport certain polar solutes and water for specific members [30,33,34]. In filamentous fungi, and in a great majority of cases, the fungal XIP subfamily consists of a single member per species [35]. In *Trichoderma harzianum*, *XIP* is constitutively transcribed and up-regulated during interaction with the plant pathogen *Fusarium solani* on olive tree roots [32]. In silico analysis showed that fungal XIP transports water and possibly other small polar molecules like H₂O₂, but intriguingly, not glycerol concerning the plant counterparts. However, unlike some AQPs or AQGs from fungi that have been functionally characterized [21,23,25,36], knowledge of XIP 3D structure, expression and physiological roles is still hypothetical. Thus, a fungal XIP prototype needs to be considered carefully.

In this context, the present study has a triple objective. After having introduced the MIP diversity in the *Trichoderma atroviride* genome, we first investigate the sequence and structural characteristic of the fungal XIP subfamily with a special focus on *TriatXIP*. Secondly, because the functions of the fungal-XIP are yet to be understood, the developmental, anatomical, physiological and biochemical effects of five XIP deletion mutants (Δ*TriatXIP*) were subsequently analyzed under specific nutritional contexts and mycoparasitism situations against different phytopathogen agents. Thirdly, the XIP member must be considered as a unit integrating a network of channels, and profiling of the expression of the six other MIP paralogues present in the *T. atroviride* genome has been recorded in the deviant mutant phenotypes subjected to different contrasting abiotic environments. We make use of a comprehensive and multidisciplinary analytical approach, together with the results of recently published studies, to better understand the roles that the fungal-XIPs play in water transport and, more broadly, in the *Trichoderma* life cycle and its sub-jacent metabolisms.

## 2. Materials and Methods

### 2.1. *Fungal Strains and Culture*

*Trichoderma atroviride* strain IMI 206040 (wild type) was used for this study. The plant pathogenic strains *Botrytis cinerea* B05.10, *Fusarium graminearum* ACCC 36966, and *Rhizoctonia solani* ACCC 36246 were used as fungal preys. Strains were routinely inoculated on potato dextrose agar (PDA; peeled potatoes, 200 g; dextrose, 10 g; distilled water, 1000 mL; pH 6.5) in Petri plates, and incubated in continuous darkness at 23 °C. Chlamydospores were observed after 20 days of incubation on PDA in continuous darkness at 23 °C, and after 14 days on two corn flour media: cornmeal and water (200 g/L), pH 4.2, and corn flour 62.86 g/L, glycerol 7.54 mL/L, pH 4.2 [37].

### 2.2. *Sequence Collection*

The nucleotide and corresponding amino acid sequences of *MIP* genes were retrieved from *Trichoderma* spp. available at the Joint Genome Institute (JGI; http://genome.jgi-psf.org/) and the NCBI GenBank GSS databases (http://www.ncbi.nlm.nih.gov/). Each putative MIP sequence was carefully scrutinized including verification of the expected MIP motifs and the prediction of the transmembrane topology with Interproscan from EMBL (http://www.ebi.ac.uk/Tools/pfa/iprscan/). Protein names and accession numbers used in this work are listed in (Figure S1).
## 2.3. Precision on MIP Naming Used in this Work

In the bibliography and for convenience, the AQP and MIP terminologies (as well as their acronyms) are considered synonymous and are used interchangeably [38]. Because the MIP superfamily from fungi includes “true” AQP with strict water channel ability (which corresponds to a strict functional terminology), we use the acronym MIP for this study in a broader sense to describe the “aquaporin” channels *lato sensu*, and the aquaporin and its acronym AQP according to its strict functional terminology.

## 2.4. In Silico Molecular Characterization of Fungal XIP

Multiple alignments of the XIP protein sequences were computed with ClustalW, and the results were viewed with Jalview (v2.10.5). The phylogenic study was carried out with the coding region of *XIP* sequences. Sequences were aligned with MUSCLE (v3.8.31) [39]. Ambiguous regions (i.e., containing gaps and/or poor alignment) were removed with Gblocks [40]. The phylogenetic tree was reconstructed using the maximum likelihood method implemented in the PhyML program (v3.1/3.0 aLRT) [41]. Reliability for the internal branch was assessed using the bootstrapping method including 1000 bootstrap replicates. Graphical representation and edition of the phylogenetic tree were performed with TreeDyn [42]. MIP from *T. atroviride* (3 AQP: AQP-31598, AQP-43816, AQP-6990; 3 AQGP: Fsp-like_39327, Fsp-like_283564, AQGP-Other_90169; [28] were added to the analysis as XIP intraspecific out-group sequences, and the sequences from *Aspergillus terreus* (AAJN01000055.1) and *Penicillium marneffei* (XP002149425.1) were used as XIP interspecific out-group references [29]. The subcellular localization of *T. atroviride* MIP was predicted using WoLF PSORT, a protein subcellular localization prediction tool available at http://www.genscript.com/wolf-psort.html. In addition, predicted localizations were confirmed with the use of DeepLoc-1.0 (http://www.cbs.dtu.dk/services/DeepLoc/) and the Cello prediction system (http://cello.life.nctu.edu.tw/). Because these three tools generated similar data, only those using WoLF PSORT were shown in this work.

## 2.5. Homology Modeling and Atomic System Building

Homology molecular modeling of *TriatXIP* from *T. atroviride* was performed with the I-TASSER (Iterative Threading ASSEMBly Refinement) program suite [43]. As for the molecular dynamics simulations, the atomic system was built up using the charmm-gui web interface [44]. The homology model of *TriatXIP* was inserted into a 3-palmitoleoyl-2-oleoyl-d-glycero-1-phosphatidylcholine (YOPC) lipids bilayer and solvated with TIP3 water. An isotonic concentration of 150 mM KCl was then generated. The all-atom simulation was performed with Gromacs (v.2018.1) [45] in a charmm36 m force field [46]. A first minimization step was followed by six equilibration steps before the 10 ns production phase. Pressure and temperature were kept constant at 1 bar and 303.15 Kelvin, respectively, using the Berendsen method during equilibration and the Parrinello-Rahman and Nose-Hoover methods during production. The Lennard-Jones interactions threshold was set at 12 ångströms and the long-range electrostatic interactions were calculated according to the particle mesh Ewald method.

Concerning the trajectory analyses, permeability coefficients (pf) were calculated according to the collective coordinate method [47,48] from a 5 ns sub-trajectory after 5 ns of additional equilibration for each protomer leading to four repetitions. Water molecules were monitored through the MDAnalysis library [49]. The pore diameter was calculated using HOLE software [50] from five frames per nanosecond of the same sub-trajectory for each protomer.

For the structure analyses, electrostatic potentials were calculated with APBS [51]. Prior to this operation, PDB2PQR [52] was used with the PARSE forcefield [53] to type the atoms of the structure for the solver. PyMOL was used to analyze and illustrate the models [54].
## 2.6. Construction of the Gene Deletion Vector, Transformation and ΔTriatXIP Mutant Validation

Genomic DNA was isolated using a hexadecyltrimethylammonium bromide-based method [55]. The 5'-flank region (1.05 kb before the start codon) and 3'-flank region (including the stop codon) of *TriatXIP* were amplified from genomic DNA of *T. atroviride* using gene-specific primer pairs (Table S1), and ups F / ups R and ds F / ds R as described previously [56–58]. Three fragment multisite gateway cloning systems (Invitrogen, Carlsbad, CA, USA) were used to construct gene deletion vectors. Gateway entry clones of the purified 5'-flank and 3'-flank PCR fragments were generated by BP reaction following the procedure described by the manufacturer (Invitrogen, Carlsbad, CA, USA). The entry clone for hygromycin resistance cassette (hygB) generated during our previous studies [56–58] from pCT74 vector [59] was used. The gateway LR recombination reaction was performed using the entry plasmid of 5'-flank, 3'-flank, hygB and destination vector pPm43GW [60] to generate the deletion vector. The deletion vector was transformed into *Agrobacterium tumefaciens* strain AGL-1 following a freeze-thaw procedure [61] and positive clones were selected on YEP (10 g/l yeast extract, 10 g/l bacto peptone, 5 g/l NaCl; pH 7.0) plates containing 35 µg/mL rifampicin (Sigma-Aldrich, St. Louis, MO, USA) and 100 µg/mL spectinomycin (Sigma-Aldrich, St. Louis, MO, USA).

*Agrobacterium tumefaciens*-mediated transformation (ATMT) was performed based on a previous protocol for *T. harzianum* [62]. Transformed strains were selected on plates containing 100 µg/mL of hygromycin (Sigma-Aldrich, St. Louis, MO, USA). Validation of homologous integration of the deletion cassettes in putative transformants was carried out using a PCR screening approach with primer combinations targeting the hygB cassette and sequences flanking the deletion cassettes as described previously [56–58]. The PCR-positive transformants were tested for mitotic stability and purified by two rounds of single spore isolation [56–58]. The transcript level of the *TriatXIP* gene in *T. atroviride* wild type and gene deletion strains were determined by RT-PCR using RevertAid premium reverse transcriptase (Fermentas, St. Leon-Rot, Germany) and their respective primer pairs (Table S1). For all analysis, five independent Δ*TriatXIP* deletion strains (named Δ*TriatXIPA* to Δ*TriatXIPE*) were used to confirm that the observed phenotypes were attributed to the deletion of *TriatXIP* but not to ectopic insertions.

## 2.7. Histological Analysis

The morphology of each strain was recorded by optical microscopy (×40 and ×100). After 10 days of cultivation, mycelium was absorbed on clear tape, stained with 1% (w/v) toluidine blue for 1 min, washed with water for 5 s, and then apposed on microscope slides. All images were processed with a Zeiss Axio Observer Z1 microscope, digital camera and Zen imaging software system (Zeiss, Jena, Germany).

## 2.8. Screening for Antagonistic Activity of *T. Atroviride*

In vitro antagonistic activity of the *T. atroviride* strains (wild type and the five Δ*TriatXIP* strains) were independently evaluated against three phytopathogenic fungi: *B. cinerea*, *F. graminearum*, and *R. solani*. A 5-mm-diameter disc of actively growing mycelia of the *T. atroviride* strains was placed on the periphery of a PDA plate. Similarly, phytopathogens were placed on the opposite side of the plate. Assays were incubated at 23 °C for 10 days. Antagonistic activity was measured as zone inhibition and growth reduction for a period of 10 days. All experiments were conducted in independent biological triplicates.

## 2.9. Biolog Phenotype Microarray Experiments

Phenotype MicroArrays™ for Filamentous Fungi (Biolog FF, Biolog Inc., Hayward, CA, USA) were used to compare the metabolic profiles on 95 single compounds of the *T. atroviride* wild type strain and the five Δ*TriatXIPs*. A list of the compounds and their biochemical families is given in Figure S2. The assimilation of compounds was reflected after
96 h of mycelial growth and quantified by measuring the optical density (OD) in the wells at 750 nm, the wavelength at which hyaline mycelium has maximum absorbance. To quantify catabolism, the wells contained a tetrazolium dye that turns into a purple insoluble precipitate when reduced due to mitochondrial activity. This was measured at the maximum absorbance of the reduced tetrazolium salt of 490 nm. Beforehand, conidia of each strain were generated on PDA plates at 23 °C in constant darkness. The FF plate procedure was carried out essentially as per the manufacturer's instructions, with slight modifications [63,64]. The experiment was performed on three replicate plates with separately prepared inocula, and incubated at 23 °C under constant darkness.

The statistical analysis and graphical output were performed using R software (version 3.6.3, R-core Team) complemented by FactoMineR (2.3) [65], FactoExtra (1.0.6) [66], tidyverse (1.3.0) [67], ggplot2 (3.3.0) [68], gridExtra (2.3) [69], egg (0.4.5) [70], lemon (0.4.3) [71], multcomp (1.4.12) [72], multcompView (0.1.8) [73], broom (0.5.5) [74] and RcolorBrewer (1.1.2) [75] libraries. Data was scaled to unit variance for PCA. Data smoothing of Figure S14 was obtained after a linear model fit. Multiple comparisons of Figure S15 were made under the generalized linear model with data fitted to a single stratum analysis of variance and the Tukey honest significant difference test. *p*-values were rounded to the third decimal place.

## 2.10. Screening for MIP Transcriptional Responses of T. atroviride in Various Nitrogen and Carbon Amended Media

A 5 mm disc of the wild type and the five ΔTriatXIP strains were inoculated on PDA medium amended with 0.5 M of D-Glucose, D-Ribose, D-Mannose, NO₃NH₄, Urea, Glycerol, D-Mannitol, L-Asparagine and L-Aspartic acid, and incubated at 23 °C under constant darkness for five days. At the end of the experiment, the mycelia were removed from the plates and stored at -80 °C for further molecular analysis. All experiments were conducted in independent biological triplicates.

## 2.11. RNA Extraction and Quantitative RT-PCR

Total RNA was extracted from mycelia after being ground to a fine powder in liquid nitrogen and treated with a CTAB extraction buffer (Bromide Cetyltrimethylammonium) as previously described [28]. RNA concentrations were determined by spectrophotometry at OD 260/280 (spectrophotometer ND-1000, Nanodrop, France), and quality was checked by using 2% TAE/agarose electrophoresis. Two µg of total RNA were reverse-transcribed with Oligo-dT using the SuperScript® III First-Strand Synthesis System for RT-PCR (Invitrogen, Carlsbad, CA, USA). cDNA was diluted 40-fold with sterile water. The abundance of MIP-related transcripts was determined by real-time qPCR with a MyiQ instrument (Bio-Rad, Hercules, CA, USA). MIP gene expression levels were calculated by the 2^ΔΔCT^ method [76]. PCR amplifications were done in 15 µL of PCR reaction using MESA GREEN qPCR MasterMix Plus (Eurogentec, Liège, Belgium) from 2 µL of 40-fold diluted cDNA template. PCR cycling conditions were 1 cycle for 3 min at 95 °C, followed by 35 cycles for 15 s at 95 °C, 15 s at 54 °C, and 20 s at 72 °C. PCR reactions were ended by generating a dissociation curve of 56 to 94 °C with an increment of 0.5 °C/10 s to check for primer dimers and nonspecific amplification. All PCR technical samples were assayed in duplicate. The geometric mean of Ct of five genes encoding to 18SrRNA, 28SrRNA, Tubulin, Actin, and Gpdh were used as internal references to normalize MIP expression for their stable expressions during fungus treatment. These genes were chosen from a panel of widely used housekeeping genes [28–77] using the software application BestKeeper v1 [78]. Each referrer belongs to protein families involved in different cell processes in order to minimize the risk of co-regulations. Specific primer pairs for each MIP member were designated in consensus zones after the alignment of MIP sequences retrieved from Trichoderma spp. with the Primer3plus application (http://www.bioinformatics.nl/primer3plus). Specific amplification of only one desired band was observed using each primer combination for qRT-PCR analysis, and PCR efficiency was 100 ± 3% for all primer pairs. Primer pairs
are listed in (Table S1). For statistical analysis, the biological replicates correspond to the five *TriatXIP* null mutants and three wild type independent mycelial cultures. Gene expression levels between parental and mutant strains were computed by a one-way analysis of variance (ANOVA) followed by a Tukey's honest significant difference (HSD) post hoc test (*p* < 0.05).

# 3. Results and Discussion

## 3.1. *Tichoderma atroviride* MIPsub-class

### 3.1.1. *T. atroviride* XIP, AQP and AQGP Inventory

In recent years, the physiological relevance of the MIP-facilitated transmembrane diffusion of water and various solutes has been largely acknowledged in Eukaryotes such as plants and animals. Comparatively, such information from fungi is lacking, except for ecotomycorrhizal (EcM) fungi that develop mutualistic associations with the roots of a large range of plants. Fungi exhibit great phenotypic plasticity in their responses to their immediate growth environment, correlated with remarkable nutritional versatility that models a great variety of polyphagical lifestyles. MIP integrates a wide and complex network of transmembrane channels that provides controlled cell internalization of inorganic and organic matter from the immediate environment. This matter, regarded as the molecular building blocks of life, is of prime importance in the context of nutrition (metabolism) as well as for maintaining cell homeostasis, osmoregulation and turgor, and determining the final shape of the cells. The MIP superfamily from the *Trichoderma* genus includes the functional aquaporins (AQPs) (*stricto sensu*) with the "orthodox fungal water channels", the functional aquaglyceroporins (AQGs) with the *Fps*-like aquaglyceroporins and the facultative fungal aquaporins (or "other" aquaglyceroporins), and the presumed fungal (un-)orthodox aquaporins X-intrinsic protein (XIP) (Figure 1, Figures S1,S2,S3) [28–79]. Concerning *T. atroviride*, its genome contains seven *MIP* homologous genes that cover the three MIP sub-families: three *AQGs* including one "*other AQGP*" (JGI Protein id: 90169; then called *TriatOtherAQGP-90169*) and two *Fps*-like *AQGs* (*TriatFpsAQGP-39327*, *TriatFpsAQGP-283564), three *classic AQPs* (*TriatAQP-31598*, *TriatAQP-43816*, *TriatAQP-6990*), and one *XIP* (*TriatXIP-319992*) (Figures S4,S5). Every related protein falls under the MIP phenotype: each member is a small integral protein, ranging from 276 to 344 AA, with theoretical MWs from 30.132 to 37.57 kDa and *pI*s from 6.60 to 9.14, except for *TriatFpsAQGP-283564* with 594 AA, 65.625 KDa, and a *pI* 5.88, respectively; Table 1), for which every 3D prediction is confidently an aquaporin structure (Figure S6).

<table><caption>**Table 1.** Nomenclature and protein properties of *Triat* MIPs from *Trichoderma atroviride*.</caption><thead><tr><th><sup>a</sup>Loci<br>Proposed gene name/Locus</th><th>Size <sup>b</sup>MW<br>(aa) (kDa)</th><th><sup>b</sup>pI</th><th><sup>c</sup>TMH</th><th><sup>d</sup>SubCL</th><th colspan="2"><sup>e</sup>NPA</th><th><sup>f</sup>ar/R SF</th></tr></thead><tbody><tr><td colspan="8">Fungal Uncategorized X-intrinsic proteins (XIPs)</td></tr><tr><td><i>TriatXIP-319992</i></td><td>301 32.103 8.75</td><td>6</td><td>PM</td><td>NPT</td><td>LB</td><td>LE</td><td>NPA N-S-Q-K</td></tr><tr><td colspan="8">Aquaporin (AQGs)</td></tr><tr><td><i>TriatAQGP-31598</i></td><td>276 30.132 7.70</td><td>6</td><td>PM</td><td>NPA NPV</td><td>Y-M-A-R</td><td></td><td>276</td></tr><tr><td><i>TriatAQGP-43816</i></td><td>312 32.729 6.60</td><td>6</td><td>PM</td><td>NPA NPV</td><td>F-H-T-R</td><td></td><td></td></tr><tr><td><i>TriatAQGP-6990</i></td><td>303 32.453 9.14</td><td>6</td><td>PM</td><td>NPA NPA</td><td>F-H-T-R</td><td></td><td></td></tr><tr><td colspan="8">Aquaglyceroporin (AQGs)</td></tr><tr><td><i>TriatOtherAQGP-90169</i></td><td>344 37.57 6.75</td><td>6</td><td>*</td><td>PM</td><td>SPA NMA</td><td>F-A-T-R</td><td></td></tr><tr><td><i>TriatFpsAQGP-39327</i></td><td>299 32.653 7.60</td><td>6</td><td>PM</td><td>NPA NLA</td><td>W-G-Y-R</td><td></td><td></td></tr><tr><td><i>TriatFpsAQGP-283564</i></td><td>594 65.625 5.88</td><td>6</td><td>PM</td><td>NPT NPS</td><td>W-T-A-R</td><td></td><td></td></tr></tbody></table>

<sup>a</sup>Loci, Protein IDs and AQP types are based on JGI assembly. <sup>b</sup>MW, Protein molecular weight; <sup>c</sup>pI, protein isoelectric point.

<sup>d</sup>TMH, Number of transmembrane helices predicted by TMHMM and SOSUI analysis tools; <sup>e</sup>regions adjusted according to alignment with characterized orthologs [28,29]. <sup>dSubCL</sup>, Predicted subcellular localization using Wolfpsort server: PM,
plasma membrane; Details in Figure S7. eNPA, Asparagine, Proline, Alanine. far/R SF, ar/R Selectivity Filter (H2-H5-LE1-LE2).

In this respect, their 3D structure resembles that of an hourglass formed by six trans-membranes α-helices (TMH1 to TMH6) connected by five inter-helical loops (Loop A to E). The N- and C-termini are localized intracellularly. At the center of the pore formed by the six TMH helices, two distinct constriction structures are formed: one helicoidal with two highly conserved "NPA" (Asn-Pro-Ala) motifs, and the second called "ar/R" selectivity filter (ar/R SF). The latter is composed of four specific amino acid residues that spread along the MIP primary structure but whose side chains fit into the channel (Table 1). These two constrictions act as a selective filter, determining the solute permeability of the MIPs [80]. Finally, the in silico prediction of the spatial expression of *T. atroviride* MIPs shows that they are predominately found in plasma membranes, but may also be located in endoplasmic reticulum, peroxisomes, mitochondria and vacuoles (Figure S7).

MIPs in pathogenic fungi may act as attractive targets for antifungal drugs [35]. In this respect, knowledge emerges on certain fungal AQP and AQGP counterparts, which is not the case for the "uncharacterized X-intrinsic protein" subgroup. The fungal *XIP*, with *TriatXIP* (protein ID 319992) as a representative member, is the core of our issue for this present work.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Phylogenetic analysis of *XIP* from the *Trichoderma* genera and selected fungal MIP. AQPs, classical aquaporins; AGPs, aquaglyceroporins; XIP, X-intrinsic protein. Red and blue squares indicate aquaglyceroporin and aquaporin sub-family nodes, respectively. The bootstrap values indicated at the nodes are based on 1000 bootstrap replicates. Branch values lower than 50% are hidden. The distance scale denotes the evolutionary distance expressed in number of nucleic acid substitutions per site. *MIP* sequences from *T. harzianum* (CBS 226.95 v1.0 as reference, JGI) were used as intraspecific out-group references [28], and *XIP* sequences from *Aspergillus terreus* and *Penicillium marneffei* were used as interspecific out-group references [29]. *AQP*, *AQGP* and *XIP* sequences and accession numbers attached after each species name are listed in Figure S1.</figcaption></figure>
### 3.1.2. Sequence and Structural Characteristic of TriatXIP

In light of the available sequenced genomes from fungi, and coinciding with what is commonly observed in the genomes of several filamentous fungi [28–79], one XIP gene is present in the *T. atroviride* genome (named *TriatXIP*). *TriatXIP* is made up of four introns and six exons, resulting in a CDS of 906 nucleotides of length. *TriatXIP* is predicted to encode a protein of 301 amino acids, named *TriatXIP*, with a theoretical molecular mass of about 32 kDa which typically corresponds to an MIP protomer. Like any typical MIP, the *TriatXIP* protein protomeric structure is highly conserved and resembles an “hour-glass model” (Figure S6). It is composed of six transmembrane alpha-helix (TMH1-TMH6) connected by five extramembrane loops (LA-LE). The central pore comprises two additional and short helical segments (HB and HE) which both possess a highly conserved “Asn-Pro-Ala” (NPA) motif that forms a narrow region of the pore (NPT on HB and NPA on HE in *TriatXIP*). Because of the positive electrostatic field generated locally by the hemi-helix dipole, this region is believed to be essential for proton exclusion [81]. The substrate selectivity of MIP mainly depends on another constriction: the aromatic/arginine (ar/R) selectivity filter composed of four residues [82–85].

In order to access the atomic details of water transport through *T. atroviride* XIP, a homology model of its three-dimensional structure was built following the same process as in our previous work [28]. However, this time, the model was integrated into an atomic system mimicking the cellular conditions encountered by this transmembrane protein family (i.e., inserted in a lipidic bilayer and solvated with water and 150 mM KCl ions) and molecular dynamics were performed.

The ar/R constriction residues of the *TriatXIP* are displayed in Figure 2a, S5 and Table 1. The first three residues, N81, S211 and Q225 are typical of XIP ar/R constriction and correlate with previously observed data on *T. harzianum* [28]. However, the conserved Arginine (R) is replaced by a Lysine (K) on *T. atroviride* XIP. While the Lysine side chain is also positively charged, its interactions with water will not be exactly the same as the Arginine guanidinium group, and could change the transport properties of the pore. The surface electrostatic potential of the pore shows the expected radiation of the ‘ar/R Arginine or Lysine’ side chain positive charge in the upper part of the channel (Figure 2b) as well as the positive electrostatic field located at the center of the pore and generated by hemi-helices HB and HE dipoles. Finally, we can also observe the local negative regions created by the oxygens of pore lining residues which correspond to protein-water interaction sites (Figure 2b). The mean (± standard error) pore diameter for *T. atroviride* at the ar/R constriction is 1.86 ± 0.1 Ångströms which is close to the *T. harzianum* XIP pore diameter at this same location of 1.8 Ångströms [28]. This substitution by Lysine seems to be exclusive to these few fungal-XIP [35]. It is not observed in MIPs from plants, animals, insects or Prokaryotes, where Arginine is replaced by Asparagine or more marginally, inter alia, by Valine, Cysteine or Serine [29,30,86]. It will be very interesting to understand what kind of functional effect this residue change might have on the protein and, ultimately, on the fungal behavior.

We can also notice the absence of the salt bridge on the HE hemi-helix which is thought to allow for a small shift of the helix position and, hence, to increase pore size in AQGs [28]. The conserved Aspartate located one position after the ‘ar/R Arginine’ in AQGs is replaced by a Cysteine in *T. atroviride* XIP, however this Aspartate and the Arginine situated one helix turn after are conserved in the three AQGs of *T. atroviride* (Figures S2,S3,S4).

From the 5 ns of sub-trajectory simulation produced and sampled from a 100 ns full simulation result [47], we can observe water permeation as well as a typical single file organization of water molecules inside the pore, especially in the NPA region, the vestibules looking a little more disorganized than in strict AQPs (see movie Figure S9, Figure S8). Using the collective coordinate method [47,48], water permeability coefficients (*pf*) were calculated for each protomer. The mean (± standard error) *pf* is 2.7 ± 0.4·10⁻¹⁴ cm³s⁻¹, which is typical for *pf* values of AQPs computed from molecular dynamics [87–89].
To conclude, based on the similar electrostatic profiles of strict AQPs and *Triat*XIP pores and the similar organization of water continuum and permeability coefficients, it seems very likely that *Triat*XIP is an efficient water channel. However, the wider ar/R constriction and the replacement of its typical Arginine by Lysine might generate functional differences between *Triat*XIP and strict AQPs, and might favour small polar solutes other than water.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Structural analysis of *Triat*XIP. (a) Schematic representation of one protomer of *Triat*XIP as seen from the extra-cellular compartment. The backbone is represented in a diagram and the pore lining residues in sticks. Residues of the ar/R SF constriction are labeled. K230 is unusual as it replaces the very conserved Arginine of the constriction. (b) Schematic representation of the tetramer with a focus on one protomer for which the pore surface is figured. Electrostatic potential (KT/e) was computed with APBS (the atoms were previously typed with the PARSE forcefield in PDB2PQR) and projected on the pore surface. The channel surface has a typical positive potential at its NPA region, where the two helices (represented by cylinders) meet and in the extra-cellular vestibule because of the positively charged lysine 230 side-chain. A water molecule is figured in purple interacting with the two Asparagines of the NPA region. Local electro-negative potentials are located at the interaction sites for water molecules with carbonyl groups of the backbone. The two Asparagines of the NPA region, the Lysine 230 of the ar/R constriction and the carbonyl groups of pore-lining residues are represented with sticks. All representations were made with PyMol software [54].</figcaption></figure>

## 3.2. Functional Characteristics of *T. atroviride* *Triat*XIP

### 3.2.1. General Issue and Objectives of this Item

There is an increasing understanding that the MIP superfamily performs an astonishing variety of physiological functions that define the biological activities of organisms. Among the fungal MIPs, XIPs are not functionally characterized to date, which means their role is still obscure. Promising channel functions predicted by structure modeling reveal that the plasma membrane *Triat*XIP proteins might act as strictly aquaporin-like proteins. However, we cannot exclude that *Triat*XIP proteins can transport in vivo specific solutes as some orthodox aquaporins and/or aquaglyceroporins, possibly due to their ability to interact with potential regulators.
Furthermore, *TriatXIP* is part of the five *TriatMIP* genes that are transcripted in *T. atroviride* (this item will be discussed in greater detail later in this article). It is hazardous to correlate the steady state levels of the transcriptional expression of a gene with the putative functions of its related protein; however, it is very attractive to consider the XIP channel among the potential candidates facilitating the diffusion of matter across the plasma membrane involved in fungal metabolism and growth. Accordingly, it could be expected (but nonetheless questionable) that the lack of a single aquaporin belonging to a superfamily directly impairs a wild phenotype—in this case, *T. atroviride*.

Indeed, should *TriatXIP* be significantly expressed and predicted to transport water, its failure among the fungal MIP network may impact the overall intracellular water content and intracellular turgor pressure. Such a hydric upheaval, even though infinitesimal, can reset the overall metabolism of the organism. To help shed light on some potential XIP roles in fungi, five *TriatXIP* deletion mutants (Δ*TriatXIP*) from *T. atroviride* were designed, and their phenotypes, cell metabolism, and subjacent MIP transcriptional profiles were analyzed in strains developing in stressful biological conditions.

### 3.2.2. *TriatXIP* Disruption Moderately Impacts Fungal Growth but Annihilates Chlamydospore Formation

Despite the transcript accumulations of *TriatXIP* in the wild type, *TriatXIP* deletion does not impact fungal growth during the first 48 h of cultivation under stress-free growing conditions (Figure 3a). Contrasted mycelial growths in Δ*TriatXIP* are significantly distinguishable from the third day of cultivation. Ultimately, all *Trichoderma* strains intensively cover the entire surface of the Petri dishes after four days of cultivation for the wild strain and five days for the transformants. This is consistent with studies carried out on yeast or filamentous organisms where nullified *AQP* or *AQGP* gene expressions result in growth failure and are connected to an assumed decrease in plasma-membrane water permeability [23,90,91]. In this context, considering the pivotal role of water channels in maintaining cellular homeostasis and metabolism integrity, *TriatXIP* disruption may result in osmoregulatory imbalance that significantly impacts *T. atroviride* growth.

At a microscopic level, subtle perceptible differences in the hyphal morphology of transformants compared to the wild type strain are revealed (Figure 3b, Figure S10). Concerning the conidiospores, their appearance is similar between strains. The most notable difference is the very low presence of chlamydospores (or even a complete lack according to the experiment; Figure S10) with the mutants after 20 days of cultivation on a PDA medium and specific corn flour media [92]. Furthermore, most chlamydospores display abnormal morphologies and/or are somewhat smaller (≈6.5 µm) than those observed with the wild strain (≈8.5 µm) (Figure 3c). These data would indicate that *TriatXIP* plays a role in the chlamydosporogenesis in *T. atroviride*, thus directly interfering with one of the asexual reproduction processes of the fungi. Chlamydospores usually form from the mycelium or certain conidia and play a major role as survival structures for various pathogenic and saprophytes fungi. However, even today, very little is known about the molecular networks involved in the formation and maturation of chlamydospores and the natural factors that enhance their formation, especially for the *Trichoderma* genus [92,93]. A genome-wide transcriptional analysis of these mutants will offer an opportunity to further understand the chlamydosporogenesis process in *Trichoderma*.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Phenotypes of *T. atroviride* wild type and the five Δ*TriatXIP* mutants upon growth on PDA medium. (a) The mycelial growth corresponds to the mean of three independent biological experiments, and bars represent the biological standard deviation. *, Data are statistically significantly different between the wild type and each transformant, as verified by one-way ANOVA analysis followed by Tukey's post hoc test (p < 0.05). (b) Light microscopy of hyphae and spores of the parental strain and the Δ*TriatXIPa* mutant for illustration after 20 days of cultivation on PDA media in darkness at 23 °C. The scale bar represents 10 or 20 µm. Samples were colored with toluidine blue (1% *w/v*) for 1 min prior to the observations. ChlSp, Chlamydospores; Sp, spores; Hy, hyphes. Each Δ*TriatXIP* are shown in Figure S10. (c) Example of morphology of Chlamydospores (ChlSp) and conidiospores (Sp) from the parental strain and the five Δ*TriatXIP* mutants after 14 days of cultivation on corn flour medium (corn flour 62.86 g/L, glycerol 7.54 mL/L, pH 4.2) in darkness at 23 °C. The scale bar represents 10 µm. Samples were colored with toluidine blue (1% *w/v*) for 1 min prior to the observations by light microscopy.</figcaption></figure>

### 3.2.3. Evaluation of *T. atroviride* Wild Type and ΔTriatXIP Mutants as Mycoparasites

A great deal of literature describes *T. atroviride* as a mycoparasite, drastically reducing the mycelial growth and conidiation of a panel of economically important aerial and soil-borne phytopathogenic fungi [94]. In this respect, this strain is now increasingly used as a bioprotectant in organic cropping systems around the world. *Trichoderma* biocontrol properties are described as a complex combination of several mechanisms, including nutrient competition and direct mycoparasitism, which involves the production of antifungal metabolites and cell-wall degrading enzymes. Although very little is known about the contribution of MIPs in establishing and deploying mycoparasitic mechanisms of biocontrol agents [91–95], it has recently been suggested that aquaporins are involved in the *Trichoderma*-prey interactions [28,95].

To substantiate the role of *TriatXIP* during antagonism, *T. atroviride* wild type and the five Δ*TriatXIPs* were confronted with *B. cinerea*, *R. solani*, and *F. graminearum*, given that *T. atroviride* stops the growth of these phytopathogens. The results show that, despite a significant slowdown in mycelia growth when wild type, transformants and prey mycelia meet together at four days of cultivation; no specific knock-on effect in growth kinetics for the Δ*TriatXIPs* is observed throughout the confrontation (i.e., before and until the full draping of the mycelia) (Figure 4; Figure S11). Similarly, no differences in growth kinetics of the fungal prey fungi are recorded during the first steps of interaction with the
*ΔTriatXIPs* in comparison with the wild type, nor after six days of dual confrontation. When contact is established between strains, prey growth is definitively stopped. The data clearly show that *ΔTriatXIPs* are capable of annihilating phytopathogen growth as with the wild type. This demonstrates that deleting *TriatXIP* does not cause the apparent failure in transformants in an antagonistic capacity, nor does it alter their ability to confront the potential aggressivity of the preys, at least with the phytopathogen strains used.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Mycoparasitic activity of the five ΔTriatXIP mutants and *T. atroviride* wild strain against *Rhizoctonia solani*, *Botrytis cinerea* and *Fusarium graminearum* as phytopathogenic hosts. (a) The mycelial growths on PDA medium correspond to the mean of three independent biological experiments, and bars represent the biological standard deviation. *, Data of the five ΔTriatXIP growth rates, and each wild strain (antagonist and phytopathogenic hosts) placed in confrontation situations were statistically and significantly different from that of each strain upon stress-free cultivation at the corresponding time point, as verified by one-way ANOVA analysis followed by Tukey's post hoc test (p < 0.05). Statistical calculations carried out included (*T. atroviride* WT + Host WT vs. *T. atroviride* WT), (ΔTriatXIPs + Host WT vs. ΔTriatXIPs), (Host WT + *T. atroviride* WT vs. Host WT), and (Host WT + ΔTriatXIPs vs. Host WT). (b) Plate confrontation assays of the *T. atroviride* parental strain (first plate) and the ΔTriatXIPa mutant (second plate; ΔTriatXIPa is used for illustration, all confrontations implicating the five ΔTriatXIPs are presented in Figure S11, and the phytopathogenic hosts (third plate). Pictures were taken six days after inoculation of the two fungi on opposite sides of the plate.</figcaption></figure>

### 3.2.4. Metabolic Analysis of ΔTriatXIP Mutants

The XIP loss of function in ΔTriatXIP mutants slightly delays biomass accumulation. At the microscopic level, these mutants express discrete differential phenotypic traits on the hyphae and nearly zero chlamydosporogenesis generating mostly chlamydospores with anomalous morphologies. In addition to the fact that fungal cells must continuously regulate the water permeability of their membranes in order to sustain water uptake
during growth and/or osmotic adjustment [20], a phenomena in which *TriatXIP* may play
a role; it is possible to hypothesize that *TriatXIP* can play a wider role by regulating C-
and N-matter flow, notably in metabolic homeostasis. To answer this question, we ana-
lyzed the metabolic profiles of the Δ*TriatXIPs* in comparison with the wild type strain. To
do so, we placed these analyses in a larger context where, in a natural composting envi-
ronment, the microorganisms are capable of metabolizing and competing for a high vari-
ety of vital compounds (or macronutrients). Filamentous fungi, such as the *Trichoderma*
genus, play a vital role in the ecological equilibria between microorganisms, in particular
by being responsible for a significant proportion of the hydrolysis of soil resources [96].
Considering that macronutrients are one of the major determinants of fungal phenotype,
we used the high-throughput method of the Phenotype MicroArray™ (PM) system (Bi-
olog FF) to explain phenotypic changes of the Δ*TriatXIP* strains to contrasted nutrient uti-
lizations. Based on the absorbance measurements at OD<sup>750</sup>, which reflect biomass produc-
tion (i.e., mycelial growth), and OD<sup>490</sup>, which reflects the respiration rate and the substrate
use (i.e., catabolism), this phenotype microarray assay has been introduced as a way to do
high-throughput screening of diverse basic nutrient sources as well as evaluate growth
and metabolic characterizations of various microorganism strains, including assorted
*Trichoderma* spp. [14–97]. However, the interpretation of each between-strain difference in
the speed of substrate use is not straightforward and has to be done with great care: each
change in color could be produced by a redox shift in tetrazolium dye due to cell respira-
tion (NADH production) that can interplay with i) inherent variations in some enzymatic
activities of the fungi during its growth, ii) an increase of turbidity due to fungal body
proliferation, or iii) a change in medium color after producing specific metabolites in re-
sponse to the candidate metabolite of the microplate and/or by the production of meta-
bolic wastes by the fungi. Nevertheless, and despite the methodological uncertainties, a
change in absorbance values can be interpreted as a useful indicator of metabolic activities
that can suggest general metabolic trends. In this respect, this approach has proved effec-
tive by offering many cell-wide perspectives, especially in understanding how a gene in-
itially encoded in the genome and under control of the environment is displayed in fine
at the cellular level, notably throughout the phenotypes of an organism [98,99].

In the present work, the metabolic profiles tied to 95 C- and N-compounds were an-
alyzed (metabolites listed in Figure S12). As observed for several fungal strains (e.g., *Acre-
monium*, *Aspergillus*, *Neurospora*, and so on), and in line with the metabolic profiles rec-
orded for species from the genera *Trichoderma* [100], *T. atroviride* exhibits a wide range of
substrate assimilation (respiration) and rapid growth (biomass) (both being significantly
colinear—Figure S13) which covers the major biochemical classes (Figure 5; Figure S14).
Such physiological aptitudes can be brought into correlation with the great advantage for
this genus of colonizing many worldwide ecological niches. In detail, the most assimilated
carbon sources are carbohydrates (oses and osides) followed by polyols, carboxylic acids,
and amino acids. The nucleic acids and the amino- and amido-compounds are consumed
at a much lower level. Still to date, the metabolic richness structuring the functional di-
versity of microbial communities (in soil and/or within a "host" organism) are far from
being listed.

First of all, this data integrates the paradigm that carbohydrates contribute most to
the overall metabolic activity of microbiota, and that the related carbohydrate regulatory
pathways appear to be modulated mainly through perturbations in the communities
[101,102].
<figure><img src="image_8.png" /></figure>
**Figure 5.** Comparative carbon source utilization profiles of the five ΔTriatXIP mutants and the *T. atroviride* parental strain. Strains were grown on 95 carbon sources (FF-plates) using the Biolog Phenotype Microarray system. The order of the organic sources is the rank of respiration rate and substrate use measured at OD490 nm on 95 organic sources and water after 96 h of incubation. Each organic source is colorized according to biochemical class. The mycelial respiration corresponds to the mean of three independent biological experiments per strain, and bars represent the biological standard deviation. **Respiration rates from the five ΔTriatXIP that were significantly different from that of the *T. atroviride* parental strain as verified by t-test (p < 0.05). *Statistical data not significantly different (p > 0.05) between the five ΔTriatXIP mutants and the parental strain, but reflect interesting metabolic trends. Statistical data are detailed for the 95 carbon sources in Figure S15.

### 3.2.5. Comparative Analysis of Transformants and Parental Strains

Concerning the comparative analysis of cell responses between the transformants and the parental strain, 16 differential metabolic patterns can be significantly identified (p < 0.05) (Figure 5, Figure S15). There is decreased metabolization in seven oses, five amino acids, three carboxylic acids and one polyol, reflecting a possible reduction in the resource uses.

Concerning the oses, disruptant phenotypes are perturbed on D-Mannose (p < 0.0001), D-Ribose (p = 0.007), N-acetyl-D-Glucosamine (p = 0.013), β-Methyl-D-Galactoside (p = 0.013), Glycerol (p = 0.015), and D-Arabinose (p = 0.03). In fungi, Mannose is an interesting candidate because it integrates multifaceted pivotal functions in cell physiology. For example, it is notably present in the composition of a panel of mannose polymers, mannans and soluble peptidomannans, which form an integral part of the cell wall surface of most fungi [103]. The Mannose residues here are crucial for maintaining the cell wall integrity and fungal viability [104]. Furthermore, Mannose is considered potent cell organic protectant, as it may rank among the most abundant polyols (with trehalose and raffinose) in cells exposed to various stress conditions. In this respect, it can constitute up to 10% of the dry weight of certain filamentous fungi [105]. As regards the nutritional context, the Mannose polymers that structure the fungal cell wall surface or those that are in the compositions of the plant cell wall pectic polysaccharides (rhamnogalacturonans) and several plant secondary metabolites (e.g., anthocyanins, flavonoids and triterpenoids) can be metabolized. The releasing D-Mannose can be used as sole carbon sources, then being metabolized to supply various fungal pathways including the pentose-phosphate cycle, glycolysis and the amino sugar and nucleotide sugar metabolisms. Similarly, N-acetylglucosamine (GlcNAc), an amide derivative of the monosaccharide glucose, is the major monomer of the polysaccharide chitin (β-(1,4)GlcNAc). Chitin is an essential structural component located in the inner layer of the fungal cell wall close to the plasma membrane. With various β-(1,3/1,6)-glucans, they constitute the structural scaffold of the cell wall, thus conferring strength and rigidity to the cell wall as a whole in most fungi [106]. Concerning β-Methyl-D-Galactoside, a methyl derivative of the monosaccharide galactose, it integrates the composition of various galacto-containing oligosaccharides and branched polysaccharides such as galactomannans and galactoglucomannans. These polymers are the main groups of plant hemicelluloses, but they have also been isolated from cell walls in several fungi. Interestingly, these two monomers are linked to Mannose metabolism through the amino sugar and nucleotide sugar metabolisms. D-Ribose and D-Arabinose (two major constituents of the plant cell wall lignocelluloses) are two of the intermediates of the pentose catabolic pathway, and in this respect, several fungi can use these aldopentoses as their sole carbon sources. In *Trichoderma*, the pentose sugar pathway involves several enzymatic steps, most of which have been well-characterized, especially those linked to the applications of certain sugars in food and the feed and biofuel industries [107–109]. Fundamentally, it generates NADPH and pentoses, some precursors that are of great importance for the biosynthesis of nucleic acids and amino acids. Interestingly, the mutants show a decrease in metabolic activity when they develop on the pentose pathway intermediates D-Arabinose and D-Ribose (as with D-Mannose and its derivative L-Rhamnose). It is plausible that the pentose catabolic pathway is affected in the transformants, impacting de facto their metabolic activities and redox equilibrium. It may
explain the subtle changes during fungal growth, as previously observed for *Candida albi-cans* [110]. Among the polyoses, interesting statistical trends emerge involving reducing diholosides D-Cellobiose (β-d-Glcp-1,4-d-Glcp) and D-Melibiose (α-d-Galp-1,6-d-Glcp) which noticeably alter ΔTriatXIP growth and respiration in comparison with the parental strain (p = 0.061, p = 0.087, resp.). These diholosides are produced through hydrolysis of the plant cell wall by fungi before being internalized to supply energy to cells.

When considering the group of substrates belonging to the amino acid class, and particularly those that produced the most significant biomass, significant metabolic differences can be found between wild type and transformants such as L-Alanyl-Glycine (p = 0.005), L-Aspartic acid (p = 0.013) and L-Asparagine (p = 0.032). When amino acid use induces very moderate biomass, differential growth is significantly observed on L-Phenylalanine (p < 0.0001) and Glycyl L-Glutamic acid (p = 0.006). Besides being proteinogenic amino acids, several of these amino acids can be diverted to various metabolic pathways such as the Krebs-Henseleit TCA cycle, providing cell energy. For example, L-Asparagine and its derivative L-Aspartic acid are linked to the TCA intermediate deaminated keto acid oxaloacetate. This implies that these amino acids are preferred by *T. atroviride* to supplement the supply of mineral nitrogen and energy, as previously observed with the Basidiomycota *Scleroderma sinnamariense* [111]. However, the amino acids L-Serine, L-Alanine, L-Threonine and L-Succinic acid are metabolized as often as L-Asparagine without inducing contrasting biomass between strains (Figure S15). Consequently, a possible interpretation is that the nitrogen pathways would be impacted, particularly the urea cycle, as it is intricately linked to the TCA cycle, first and foremost. This hypothesis is strengthened by the fact that ΔTriatXIP biomass and respiration give an interesting decreased tendency on L-Ornithine (p = 0.078). However, the power of the experimental design (with three repetitions for each strain) hinders any clear decision. The urea cycle integrates biochemical reactions to convert ammonia (NH₃), a toxic waste originating from the amino acid catabolism, into urea ((NH₂)₂CO), a less toxic nitrogen source for excretion. From there, under conditions of amino acid fluctuation, L-Asparagine and L-Aspartic acid, with their derivative oxaloacetate, are known for regulating various enzymes from the urea cycle, maintaining the flow of nitrogen into the urea cycle via the carbamoylation of the L-Ornithine by the Ornithine transcarbamoylase enzyme. The fact that XIP transporters can be linked to the urea cycle opens up a whole new realm of possibilities that should be explored in future experiments, especially since XIPs are known for facilitating the transport of urea [33]. Lastly, L-Asparagine is one of the alternative forms of storing nitrogen in plants, and consequently, could be a favorable source of organic nitrogen for fungal growth. Interestingly, several amino acids including L-Asparagine and L-Aspartic acid triggered fungal metabolic activity, notably by playing key roles in stimulating the expression of many genes, including those required for the production of extracellular hydrolases (i.e., proteases, chitinases and glucanases) and transporters [112]. These amino acids, as regulators, integrate the nutritional processes whatever the lifestyle of the fungi, parasitism or saprophytism [113].

To conclude, several of the C- and N-compounds mentioned above are quite common in yeast and fungi (including *Trichoderma* spp.). They contribute to many key physiological processes by modulating spore germination, (a-)sexual development and forming fruiting bodies, and enhancing bioactivity of vegetative cells. However, the interconnections between these metabolic pathways are eminently complex and warrant further discussion, especially in fungi.

Our data clearly demonstrate that the direct mutagenesis of *TriatXIP* alters central C- and N-metabolisms in *T. atroviride*. However, it is very tricky to correlate deviant metabolisms with metabolites that would not be transported by XIP a priori, due to the configuration of the canal, which is too narrow. This is all the more important since a plausible alternative mechanism deserves to be highlighted: it cannot exclude that the *TriatXIP* deletion has a direct impact on cell respiration by disrupting gas exchanges, in particular for O₂ and/or CO₂. Recent research has revealed that some human and plant MIPs are
involved in the transport of these hydrophobic gazes [114–116]. Although the underlying
functional and structural mechanisms are not yet clearly elucidated [117,118], this new
evidence opens a very interesting horizon, especially for the XIP subclass, which presents
several structural singularities and for which the transport of gaseous molecules has not
yet been characterized. Clearly, more investigation is required to confirm how *TriatXIP*
participates in the C- and N-metabolisms and energy-releasing pathways, and how this
assembly incorporates Chlamydosporogenesis.

3.2.6. *TriatXIP Transcriptional Regulation in the Wild Strain*

Irrespective of how Δ*TriatXIP* transformants grow, and even though the most signif-
icant variations in carbon and nitrogen source use profile is limited to no more than a few
molecules, the fact is that a lack of *TriatXIP* interferes with both the carbohydrate and
nitrogen pathways. Of the numerous cell functions that could be awarded to these metab-
olites, those for which they act as sensing compounds and/or growth regulators attract
special attention. For example, D-Mannose, D-Ribose, D-Cellobiose and D-Melibiose are
known to be potent natural inducer components that modulate the production and/or ac-
tivity of various classes of hydrolases in *Trichoderma* spp. [108,119,120].

In this regard, because exchanges of matter between compartments are constantly
and finely regulated, there is a strong likelihood that any fluctuating events in the imme-
diate environment of the fungi will modulate the transcriptional steady-state levels of *Tri-
atXIP*.

Consequently, for the transcriptional analysis, we used the assays inherent to some
nutritional complementations with various nitrogen (urea, NH<sub>4</sub>NO<sub>3</sub>) or C- and N-ele-
ments (D-Mannose, D-glucose, D-Ribose, L-Aspartic, L-Asparagine), and two osmopro-
tectants polyols (Glycerol, Mannitol). In the wild type, results show that *TriatXIP* is sig-
nificantly induced when PDA is complemented with NH<sub>4</sub>NO<sub>3</sub>, D-Mannose, L-Aspartic
and L-Asparagine (Figure 6a). *A contrario*, *TriatXIP* is significantly repressed by D-Glucose
and Mannitol. Glycerol, Urea and D-Ribose have no effect on *TriatXIP* expression. Scan-
ning 1.5kb upstream of the start codon of *TriatXIP* using the yeast *Saccharomyces cerevisiae*
dedicated promoter database SCPD showed that these modulations can be attributable to
an over-representation of *cis*-regulated motifs targeted by nucleic acid binding proteins
known to be involved in various carbon and nitrogen metabolic processes (data not
shown). These results support previous analyses conducted on the promoter of *XIP* from
*T. harzianum* that exhibits similar *cis*-element patterns [28].

The sensing functions dependent on specific C- and N-metabolites are possible ways
of cell regulation *TriatXIP* dependent. This hypothesis fuels very exciting fundamental
questions. Because metabolisms are intricate, our data offer a new opportunity to explore
the cellular signaling as non-linear pathways that need to be further studied in fungi, es-
pecially for the *Trichoderma* genus. However, fundamentally, the assimilation of the spe-
cific components mentioned above could be viewed through the multifunctional perme-
ant-channel abilities of the members constituting this diversified MIP superfamily to
which *TriatXIP* belongs.
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Relative transcription ratios of the expressed MIP genes from *Trichoderma atroviride* wild type and the five Δ*Triat*XIP mutant upon growth on PDA medium supplemented with various inorganic and organic compounds. (a) Differential XIP transcript accumulation in the wild strain. (b) Differential AQP and AQGP transcript accumulation in the wild type and the five Δ*Triat*XIP mutants. Differential transcript levels for each gene were estimated using real-time qRT-PCR analyses, and normalized by the expression of five housekeeping genes. Relative transcript abundance rates were obtained by the E~ΔΔCt~ method. Data correspond to means of three biological replicates for the wild strain and the mean of the five Δ*Triat*XIP mutants. Bars represent the biological standard error. *, Data are statistically significantly different between the wild and the five Δ*Triat*XIP strains, as verified by one-way ANOVA analysis followed by Tukey's post hoc test (*p* < 0.05). Expression levels are displayed by the depth of color, where red represents down-regulated regulations and blue denotes up-regulated regulations. Statistical analysis showing the significant difference between wild type and mutants is detailed in Figure S16.</figcaption></figure>

### 3.3. *The Fungal XIP, AQP and AQGP: Three Pivotal MIP Pawn Families on Complex Chessboard Membrane Transport*

In the wild strain, the *TriatXIP* gene is co-transcribed with two AQP genes (*AQP-31598* and *AQP-43816*) and two *Fps-like* genes (*AQGP-39327* and *AQGP-283564*) (Figure 6b). The *AQP-43816* and the “Other”-*AQGP-90169* do not seem to be expressed in our biological conditions. On the nutritional environments described above, the four *Triat*MIPs transcribed are up-regulated in most of the studied conditions, with the exception of AQP-31598 and AQGP-39327 which are repressed in the presence of D-Glucose and/or Mannitol. AQP-43816 and AQP-31598 are not modulated by Mannitol, as for AQGP-39327 with Urea. More importantly, these transcriptional patterns significantly contrasted with those monitored with the Δ*TriatXIP* mutants: the down-regulations of the four *Triat*MIPs expressed observed with the wild type are exacerbated with the transformants, except for AQP-43816 with Mannitol, and AQGP-283564 with D-Glucose and Mannitol, which are up-regulated (Figure 6b, Figure S16).

It is of great interest to assess the MIP gene expression in its entirety. Indeed, we recall i) the representativeness of the three fungal MIP subgroups in the transcriptional network of the fungi with at least one member per subgroup which is expressed and modulated in fluctuating environments, ii) the substantial differences in gene regulation between members from a same subgroup that share relative high sequence identity, and iii) a subcellular localization of all members that is preferentially predicted to plasma membrane. In addition, despite the significant induced expression of *TriatXIP* in the wild type strains under contrasted environmental conditions, the Δ*TriatXIP* mutants show moderate altered phenotypes and seem to conserve their mycoparasitic features. In view of all these circumstances, simplistic and then possibly erroneous conclusions would be to deduce that the *TriatXIP* gene probably does not perform essential physiological functions in the development life cycle and/or in the responses to environmental changes and/or mycoparasitism of *T. atroviride*. That can be reminiscent of the fact that XIPs represent the smallest fungal MIP subgroup and they seem to be frequently missing in certain fungal clades [20–26]. However, this would not take into account the fact that mutants have
significantly deviant C- and N-metabolisms, and they fail to produce abundant and classic chlamydospores under ad hoc conditions.

All these results increase the sense of plausible interconnected compensating effects between *Triat*MIP from a same subgroup and counterparts belonging to the three different subgroups, in particular to facilitating common transport across biomembranes of water and a variety of low molecular weight solutes which need to be identified. As a complement, these compensating events do not exclude the possibility of the divergent function of certain MIP members belonging to each subgroup, but whose combination remains essential for the metabolic integrity of the organism: first by the amount of permeable activity from each protomer, but also for the tetramer structures (homo- and hetero-oligomers) where permeability activity can be potentiated when protomers interact with each other in a cooperative mechanism, as observed for yeast, plant and mammal MIP [121–124]. To the best of our knowledge, there is no homo- or hetero-oligomer complex described in filamentous MIP, opening a vast new field to be further explored.

In terms of solute permeation and foremost with regards to water, the apparent contradiction between the constitutive accumulation of *TriatXIP* transcripts that can be modulated under various biological conditions and the lack of noticeable important growth reduction between *ΔTriatXIP* and *T. atroviride* wild type strains in presence of stressors may be explained by an overlap in water specificity of *TriatXIP* with other MIPs. Because *TriatXIP* and the two co-expressed *TriatAQPs* (three paralogs phylogenetically linked) (Figure 6b) are co-expressed (in particular *TriatAQP-31598* which displays similar transcriptional patterns to *TriatXIP* under stresses), we can conjecture that these genes are jointly required for cell tolerance and mycoparasitism in *T. atroviride*, and that the function of *TriatXIP* protein may be complemented by AQP-31598 in *ΔTriatXIP* strains. In *S. cerevisiae*, some combined activities of orthodox AQPs and the aquaglyceroporin Fps1 have been suspected to be essential for adjusting appropriate osmotic balance, in particular with regard to hypo-osmotic conditions [125,126]. Interestingly, in addition to playing an integral role in water transport, MIP have also been considered to be potential components in sensing osmotic changes in microorganisms, plants or mammals [127–131]. All fluctuating contexts that modulate MIP expression may result in a modification in plasma membrane water permeability concomitantly to a change in membrane surface tension [132]. Ultimately, because each MIP can be considered a piece maintaining the global cellular homeostasis puzzle, the decrease in *ΔTriatXIP* growth (slightly exacerbated under stresses) may result in complex osmoregulatory imbalances cumulating deviances in water availability and membrane surface tension.

In addition to water, XIPs are involved in transporting small solutes in plants, such as glycerol, urea, ammonia, H₂O₂, boron, arsenic, antimony and silicon [30, 33, 34, 83]. In terms of multifunctional permeant-channels, these overlap mechanisms can involve transporting specific solutes which play important roles in nutrient uptake and osmoregulation, such as certain linear polyols (Glycerol, Mannitol, Arabitol, Xylitol, Sorbitol) and nitrogen (Urea, NO₃⁻, NH₄⁺), to name just a few [133, 134], or certain unexpected permeants considered harmful or toxic at high dosage such as arsenite, antimonite and various by-products of the glycolysis and tricarboxylic acid cycle (Ethanol, Acetate, Lactate, Malate) [121, 135–137]. Interestingly, *ΔTriatXIP* exhibits contrasted growth phenotype in the presence of some of these elements such as glycerol, and plausibly with Xylitol (p = 0.151) (Figure 5, Figure S15). Conversely, no impact can be seen with Lactate, Malate, or Erythritol. This might imply the plausible transport of these two polyols, Glycerol and Xylitol. However, these data contrast with the structure modeling simulations that tend to demonstrate that the hourglass of *TriatXIP* presents a pore constriction of 1.86Å in diameter (which corresponds to the diameter of a water molecule). These physicochemical features are too restrictive to enable molecules other than water to pass through. It is therefore reasonable to assume that *TriatXIP* may act in conjunction with specific channels and transporters dedicated to transport these molecules. The example of glycerol diffusion clearly demonstrates this compensating event. Fps1-like and Yfl054-like
aquaglyceroporins are two AQGP subgroups present in yeasts and some filamentous fungi [18–26]. Although the Fps1-like function has not yet been fully understood, it is a glycerol-facilitator significantly regulated by osmotic changes. Also, Fps1-like protein has a constriction region above 3.4 Å wide and is defined by more hydrophobic residues, as observed for several GLP and AQGP from *T. harzianum* [28–138]. In numerous eukaryotic and prokaryotic organisms, glycerol acts like a carbon source and a major protectant compatible solute under low $a_w$ in high osmolarity environments [139]. In response to fluctuating environments where cell osmolarity is impacted, Fps1-like protein can mediate the entrance and release of Glycerol into the cell. However, Fps1-like-deleted cells are able to grow in media with glycerol as the sole carbon source, suggesting a second role of Fps1-like on glycerol uptake [140]. Similarly, GUP1 and GUP2—two active glycerol transporters—play prominent roles in glycerol fluxes in yeasts, complementing Fps1-like protein absence [141]. These compensating occurrences involve a plethora of organic molecules whose transport can be ensured by several transporters and/or channels [142]. In our work, the observed deviant phenotypes of *ΔTriatXIP* transformants may result from indirect consequences of the capture and/or metabolism of these large solutes where the homeostatic regulatory mechanisms inherent to the hydric status had to be significantly altered by the lack of *TriatXIP*. In any case, the loss of XIP activity is not lethal to the organism, possibly due to compensation events, notwithstanding the fact that the chlamydosporogenesis process is deeply impacted.

**4. Conclusion and Perspectives**

Water is the main component of all living cells, and its flow is a fundamental property of life. The MIP superfamily, including the fungal-XIP, is an active candidate in the complex network of multifunctional and/or specific channels and transporters, and this ensemble modulates the high velocity of water exchange across the biological membrane inherent to most organisms’ physiological survival processes. Together with the recently published results [26,29,31,35], our present results complement our previous works on Trichoderma MIP [28] and give a deeper insight into the fungal XIP subgroup and its possible regulation and involvement in the nutritional and mycoparasitic processes of *T. atroviride*, alongside the activities of its paralogue aquaporins and aquaglyceroporins. It is essential to consider that one of the direct consequences of the XIP loss-of-function is its effect on the expression of the other expressed *MIP*-encoding genes of *T. atroviride*, some being up-regulated and others down-regulated. Kinetics in the growth of both the mutant strains and the wild type occurred quite similarly, though delayed for the mutants, and in this respect, the *TriatMIP* transcriptional differences are not the consequence of variable rates of nutrient uptake. We therefore conclude that the fungal *XIP* gene—through the function of its related proteinaceous product—is an important variable in the balance of MIP gene expression in the fungal growth or the cell responses to changing environmental conditions in general, and in particular, during chlamydosporogenesis. However, *Trichoderma* spp. are organisms with specialized lifestyles. An attractive research priority for the future would be the prediction of different XIP ecological functions from different fungal species, supported by extensive structural analyses on the diversity of the sequences involved. Such multidisciplinary analyses will allow a better understanding of the physiological roles of this singular MIP subclass and their functioning in various physiological mechanisms in fungi.

**Supplementary Materials:** The following are available online at www.mdpi.com/2218-273X/11/2/338/s1, Figure S1: Protein MIP sequences from *Trichoderma* genus and XIPs from *Aspergillus terreus* and *Aspergillus terreus*, Figure S2: Multiple sequence alignment of MIP from *Trichoderma* genus, Figure S3: Detailed Figure S2, Figure S4: Aminoacid sequence alignment of the seven MIP from *Trichoderma atroviride*, Figure S5: Aminoacid sequence alignment of the fungal XIP from *Trichoderma* genus, Figure S6: Predicted 3D models of the seven *Trichoderma atroviride* MIP proteins, Figure S7: Details of predicted subcellular localization of *Trichoderma atroviride* MIP, Figure S8: Root mean square displacement of the alpha carbons of the tetrameric modeled structure of *TriatXIP* (without
the extra-cellular loops and intra-cellular extremities) along the whole 10 nanoseconds production trajectory, Figure S9: Movie, Water transport inside *TriatXIP T. atroviride* for 9 nanoseconds. On the left part, the backbone of the monomer is represented with ribbons, on the right part, the lysine 230 of the ar/R constriction, Asparagines 101 and 227 of the NPA region and carbonyl groups of the pore lining residues are represented with sticks. Water molecules coloured in orange or purple correspond to water molecules crossing a 15-angström long section of the pore containing both the ar/R constriction and the NPA region during the trajectory. The MDAnalysis toolkit was used to monitor water molecule displacement [49]. We can observe a single file continuum of water and the typical re-orientation of molecules at the NPA region due to the electropositive potential produced by hemi helices HB and HE dipoles. Figure S10: Light microscopy of hyphae and spores of the *Trichoderma atroviride* wild strain and the five Δ*TriatXIP* mutants, Figure S11: Confrontations of the five Δ*TriatXIP* mutants and the wild strain with *Rhizoctonia solani*, *Botrytis cinerea* and *Fusarium graminearum* on PDA medium, Figure S12: List of the 95 metabolites used for the study, Figure S13: Colinearity between biomass production (i.e mycelial growth) and respiration (i.e. catabolism) per biochemical class from the *T. atroviride* wild strain and the five Δ*TriatXIP* mutants in cultivation on the 95 metabolites using the Phenotype MicroArray™ (PM) system (Biolog FF MicroPlate, MT2 serie), Figure S14: Respiration rate of *T. atroviride* wild strain and the five Δ*TriatXIP* mutants according to the 8 biochemical classes available on the Phenotype MicroArray™ (PM) system (Biolog FF MicroPlate, MT2 serie), Figure S15: Comparative carbon source utilization profiles of the five Δ*TriatXIP* mutants and the *T. atroviride* parental strain, Figure S16: Detailed relative transcription ratios of the expressed MIP genes from *Trichoderma atroviride* wild type and the five Δ*TriatXIP* mutants upon growth on PDA medium supplemented with various inorganic and organic compounds presented in Figure 6B, Table S1: Primers used for the *TriatXIP* mutant construction, and the *Trichoderma atroviride* MIP qPCR amplifications

**Author Contributions:** Conceived and designed experiments, J.S.V.; Performed experiments and analyzed data: M.B.A., M.F., M.D., M.K., M.F., J.P., B.F., R.M. and D.A.; Validation: M.K., M.F., D.A., P.L. and J.S.V.; Formal analysis: B.F., M.K., M.F., D.A., P.L. and J.S.V.; Resources: M.D., M.K., J.P., M.F., H.C., D.A. and P.L.; Data curation: M.F., M.K., M.F., D.A. and J.S.V.; Writing-original draft preparation: J.S.V., R.M., D.A. and M.F.; Review and editing, M.F., M.D., M.K., J.L.J, M.F., B.F., A.G.D., R.M., D.A., P.L. and J.S.V.; Supervision, M.F., M.K. and J.S.V.; Project administration and Funding acquisition, J.S.V. All authors have read and collectively agreed to the published version of the manuscript.

**Funding:** This work was granted access to the French HPC resources of IDRIS under the allocation 2019-A0060710751 and 2020-AP010710612 made by GENCI. It was supported by I-SITE CAP 20-25 (ANR grant 16-IDEX-0001) Emergence 2017 from the University of Clermont-Auvergne, and by the PHC program "Uthic" from Campus France (grant 34861PF) under the joint aegis of the Tunisian Ministry of Foreign Affairs and the French Ministry of Higher Education and Scientific Research.

**Institutional Review Board Statement:** Not applicable

**Informed Consent Statement:** Not applicable

**Data Availability Statement:** All relevant data are within the paper and its Supplementary Materials files.

**Acknowledgments:** The authors wish to thank Céline Sac and Caroline Savel for their technical assistance in fungal cultivation and molecular biology, and Nicole Brunel-Michac for her technical assistance in microscopic analysis. The Biolog analyses were performed using equipment purchased within the EU Operational Programme Eastern Poland 2007–2013—Regional Laboratory of Renewable Energy, IA PAS. http://www.polskawschodnia.gov.pl/. The authors also wish to acknowledge Christie Nielsen-Chaar for generously providing the final linguistic revision of the manuscript. We are also grateful to the anonymous reviewers for their constructive comments.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

**Abbreviations:**

<table><tr><td>AQP</td><td>Aquaporin</td></tr><tr><td>AQGP</td><td>Aquaglyceroporin</td></tr></table>
<table><tr><td>ar/R SF</td><td>Aromatic/arginine" selectivity filter</td></tr><tr><td>H2O2</td><td>Hydrogen peroxide</td></tr><tr><td>MIP</td><td>Major Intrinsic Protein</td></tr><tr><td>NH4NO3</td><td>Ammonium nitrate</td></tr><tr><td>NPA</td><td>Asparagine-proline-alanine" motif</td></tr><tr><td>PIP</td><td>Plasma membrane intrinsic protein</td></tr><tr><td>TMH</td><td>trans-alpha helical transmembranes region</td></tr><tr><td>TIP</td><td>Tonoplast intrinsic protein</td></tr><tr><td>XIP</td><td>X-intrinsic protein</td></tr></table>

References

1. Samuels, G.J.; Dodd, S.L.; Gams, W.; Castlebury, L.A.; Petrini, O. *Trichoderma* species associated with the green mold epidemic of commercially grown *Agaricus bisporus*. *Mycologia* **2002**, *94*, 146–170.

2. Serda Kantarciog, A.; Celkan, T.; Yucel, A.; Mikami, Y.; Kuruglglu, S.; Mitani, H.; Altas, K. Fatal *Trichoderma harzianum* infection in a leukemic pediatric patient. *Med. Mycol.* **2009**, *47*, 207–215.

3. Molnar-Gabor, E.; Doczi, I.; Hatvani, L.; Vagvolgyiand Kredics, C. Isolated sinusitis sphenoidalis caused by *Trichoderma longibrachiatum* in an immunocompetent patient with headache. *J. Med Microbiol.* **2013**, *62*, 1249–1252.

4. Li Destri Nicosia, M.G.; Mosca, S.; Mercurio, R.; Schena, L. Dieback of *Pinus nigra* ceedlings caused by a ctrain of *Trichoderma viride*. *Plant Dis.* **2015**, *99*, 44–49.

5. Zimand, G.; Elad, Y.; Chet, I. Effect of *Trichoderma harzianum* on *Botrytis cinerea* pathogenicity. *Phytopathology* **1996**, *86*, 1255–1260.

6. Gazis, R.; Chaverri, P. Diversity of fungal endophytes in leaves and stems of wild rubber trees (*Hevea brasiliensis*) in Peru. *Fungal Ecol.* **2010**, *3*, 240–254.

7. Singh, A.; Shukla, N.; Kabadwal, B.C.; Tewari, A.K.; Kumar, J. Review on Plant-*Trichoderma*-Pathogen Interaction. *Int. J. Curr. Microbiol. App. Sci.* **2018**, *7*, 2382–2397.

8. Ben Amira, M.; Lopez, D.; Triki, M.A.; Khouaja, A.; Chaar, H.; Fumanal, B.; Gousset-Dupont, A.; Bonhomme, L.; Label, P.; Goupil, P.; et al. Beneficial effect of *Trichoderma harzianum* strain Ths97 in biocontrolling *Fusarium solani* causal agent of root rot disease in olive trees. *Biol. Control* **2017**, *110*, 70–78.

9. Oszust, K.; Cybulska, J.; Frąc, M. How do *Trichoderma* genus fungi win a nutritional competition battle against soft fruit pathogens? A report on niche overlap nutritional potentiates. *Int. J. Mol. Sci.* **2020**, *21*, 4235.

10. Gajera, H.; Rinkal, D.; Sunil, P.; Mansukh, K.; Balubhai, G. Molecular mechanism of *Trichoderma* as biocontrol agents against phytopathogen system—A review. *Curr. Res. Microbiol. Biotechnol.* **2013**, *1*, 133–142.

11. Hermosa, R.; Viterbo, A.; Chet, I.; Monte, E. Plant-beneficial effects of *Trichoderma* and of its genes. *Microbiology* **2012**, *158*, 17–25.

12. López-Bucio, J.; Pelagio-Flores, R.; Herrera-Estrella, A. *Trichoderma* as biostimulant: Exploiting the multilevel properties of a plant beneficial fungus. *Sci. Hortic.* **2015**, *196*, 109–123.

13. Druzhinina, I.S.; Seidl-Seiboth, V.; Herrera-Estrella, A.; Horwitz, B.A.; Kenerley, C.M.; Monte, E.; Mukherjee, P.K.; Zeilinger, S.; Grigoriev I.V.; Kubicek, C.P. *Trichoderma*: The genomics of opportunistic success. *Nat. Rev. Microbiol.* **2011**, *9*, 749–759

14. Oszust, K.; Pawlik, A.; Siczek, A.; Janusz, G.; Gryta, A.; Bilińska-Wielgus, N.; Frąc, M. Efficient cellulases production by *Trichoderma atroviride* G79/11 in submerged culture based on soy flour-cellulose-lactose. *BioResources* **2017**, *12*, 8468–8489.

15. Zhang, W.; Kou, Y.; Xu, J.; Cao, Y.; Zhao, G.; Shao, J.; Wang, H.; Wang, Z.; Bao, X.; Chen, G.; et al. Two major facilitator super-family sugar transporters from *Trichoderma reesei* and their roles in induction of cellulase biosynthesis. *J. Biol. Chem.* **2013**, *288*, 32861–32872.

16. Vizcaínoa, J.A.; Cardoza, R.E.; Hauserc, M.; Hermosa, R.; Rey, M.; Llobell, A.; Becker, J.M.; Gutiérrez, S.; Monte, E. ThPTR2, a di/tri-peptide transporter gene from *Trichoderma harzianum*. *Fungal Genet. Biol.* **2006**, *43*, 234–246.

17. Sun, J.; Zhang, T.; Li, Y.; Wang, X.; Chen, J. Functional characterization of the ABC transporter TaPdr2 in the tolerance of bio-control the fungus *Trichoderma atroviride* T23 to dichlorvos stress. *Biol. Control* **2018**, *10*, 004.

18. Pettersson, N.; Filipsson, C.; Becit, E.; Brive, L.; Hohmann, S. Aquaporins in yeasts and filamentous fungi. *Biol. Cell* **2005**, *97*, 487–500.

19. Engel, A.; Walz, T.; Agre, P. The aquaporin family of membrane water channels. *Curr. Opin. Struct. Biol.* **1994**, *4*, 545–553.

20. Nehls, U.; Dietz, S. Fungal aquaporins: Cellular functions and ecophysiological perspectives. *Appl. Microbiol. Biotechnol.* **2014**, *98*, 8835.

21. Li, T.; Hu, Y.J.; Hao, Z.P.; Li, H.; Chen, B.D. Aquaporin genes GintAQPF1 and GintAQPF2 from *Glomus intraradices* contribute to plant drought tolerance. *Plant Signal Behav.* **2013**, *8*, e24030.

22. Li, T.; Hu, Y.J.; Hao, Z.P.; Li, H.; Wang, Y.S.; Chen, B.D. First cloning and characterization of two functional aquaporin genes from an arbuscular mycorrhizal fungus *Glomus intraradices*. *New Phytol.* **2013**, *197*, 617–630.

23. An, B.; Boqiang, L.; Hu, L.; Zhanquan, Z.; Guozheng, Q.; Shiping, T. Aquaporin 8 regulates cellular development and reactive oxygen species production, a critical component of virulence in Botrytis cinerea. *New Phytol.* **2016**, *209*, 1668–1680.
24. Turgeman, T.; Shatil-Cohen, A.; Moshelion, M.; Teper-Bamnolker, P.; Skory, C.D.; Lichter, A.; Eshel, D. The role of aquaporins in pH-dependent germination of *Rhizopus delemar* spores. *PLoS ONE* **2016**, *11*, e0150543.

25. Ding, M.; Li, J.; Fan, X.; He, F.; Yu, X.; Chen, L.; Zou, S.; Liang, Y.; Yu, J. Aquaporin1 regulates development, secondary metabolism and stress responses in *Fusarium graminearum*. *Curr. Genet.* **2018**, *64*, 1057–1069.

26. Xu, H.; Cooke, J.E.K.; Zwiazek, J.J. Phylogenetic analysis of fungal aquaporins provide insight into their possible role in water transport of mycorrhizal associations. *Botany* **2013**, *91*, 495–504.

27. Danielson, J.A.H.; Johanson, U. Unexpected complexity of the aquaporin gene family in the moss Physcomitrella patens. *BMC Plant Biol.* **2008**, *8*, 45.

28. Ben Amira, M.; Mom, R.; Lopez, D.; Chaar, H.; Ali, K.; Pujade-Renaud, V.; Fumanal, B.; Gousset-Dupont, A.; Bronner, G.; Label, P.; et al. MIP diversity from *Trichoderma*: Structural considerations and transcriptional modulation during mycoparasitic association with *Fusarium solani* olive trees ». *PLoS ONE* **2018**, *13*, e0193760.

29. Gupta, A.B.; Sankararamakrishnan, R. Genome-wide analysis of major intrinsic proteins in the tree plant *Populus trichocarpa*: Characterization of XIP subfamily of aquaporins from evolutionary perspective. *BMC Plant Biol.* **2009**, *9*, 134.

30. Lopez, D.; Bronner, G.; Brunel, N.; Auguin, D.; Bourgerie, S.; Brignolas, F.; Carpin, S.; Tournaire-Roux, C.; Maurel, C.; Fumanal, B.; et al. Insights into Populus XIP aquaporins: Evolutionary expansion, protein functionality, and environmental regulation. *J. Exp. Bot.* **2012**, *63*, 2217–2230.

31. Venkatesh, J.; Yu, J.W.; Gastonb, D.; Park, S.W. Molecular evolution and functional divergence of X-intrinsic protein genes in plants. *Mol. Genet. Genom.* **2015**, *290*, 443–460.

32. Lopez, D.; Ben Amira, M.; Brown, D.; Muries, B.; Brunel-Michac, N.; Bourgerie, S.; Porcheron, B.; Lemoine, R.; Chrestin, H.; Mollison, E.; et al. The *Hevea brasiliensis* XIP aquaporin subfamily: Genomic, structural and functional characterizations with relevance to intensive latex harvesting. *Plant Mol. Biol.* **2016**, *91*, 375–396.

33. Bienert, G.P.; Bienert, M.D.; Jahn, T.P.; Boutry, M.; Chaumont, F. Solanaceae XIP are plasma membrane aquaporins that facilitate the transport of many uncharged substrates. *Plant J.* **2011**, *66*, 306–317.

34. Noronha, H.; Araújo, D.; Conde, C.; Martins, A.P.; Soveral, G.; Chaumont, F.; Delrot, S.; Gerós, H. The grapevine uncharacterized intrinsic protein 1 (VvXIP1) is regulated by drought stress and transports glycerol, hydrogen peroxide, heavy metals but not water. *PLoS ONE* **2016**, *11*, e0160976.

35. Verma, R.K.; Prabh, N.D.; Sankararamakrishnan, R. New subfamilies of major intrinsic proteins in fungi suggest novel transport properties in fungal channels: Implications for the host-fungal interactions. *BMC Evol. Biol.* **2014**, *14*, 173.

36. Dietz, S.; von Bülow, J.; Beitz, E.; Nehls, U. The aquaporin gene family of the ectomycorrhizal fungus *Laccaria bicolor*: Lessons for symbiotic functions. *New Phytol.* **2011**, *190*, 927–940.

37. Li, Y.Q.; Song, K.; Li, Y.C.; Chen, J. Statistical culture-based strategies to enhance chlamydospore production by *Trichoderma harzianum* SH2303 in liquid fermentation. *J. Zhejiang Univ.-Sci. B (Biomed Biotechnol.)* **2016**, *17*, 619–627.

38. Perez Di Giorgio, J.; Soto, G.; Alleva, K.; Jozefkowicz, C.; Amodeo, G.; Muschietti, J.P.; Ayub, N.D. Prediction of aquaporin function by integrating evolutionary and functional analyses. *J. Membr. Biol.* **2014**, *247*, 107–125.

39. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. *Nucleic Acids Res.* **2004**, *32*, 1792–1797.

40. Castresana, J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. *Mol. Biol. Evol.* **2000**, *17*, 540–552.

41. Guindon, S.; Dufayard, J.F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New Algorithms and Methods to Estimate Maximum-Likelihood Phylogenies: Assessing the Performance of PhyML 3.0. *Syst. Biol.* **2010**, *59*, 307–321.

42. Chevenet, F.; Brun, C.; Banuls, A.L.; Jacq, B.; Chisten, R. TreeDyn: Towards dynamic graphics and annotations for analyses of trees. *BMC Bioinform.* **2006**, *10*, 439.

43. Yang, J.; Yan, R.; Roy, A.; Xu, D.; Poisson, J.; Zhang, Y. The I-TASSER Suite: Protein Structure and Function Prediction. *Nat. Methods* **2015**, *12*, 7–8.

44. Wu, E.L.; Cheng, X.; Jo, S.; Rui, H.; Song, K.C.; Dávila-Contreras, E.M.; Qi, Y.; Lee, J.; Monje-Galvan, V.; Venable, R.M.; et al. CHARMM-GUI Membrane Builder toward realistic biological membrane simulations. *J. Comput. Chem.* **2014**, *35*, 1997–2004.

45. Abraham, M.J.; James, M.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindahl, E. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. *SoftwareX* **2015**, *1-2*, 19–25.

46. Huang, J.; Rauscher, S.; Nawrocki, G.; Ran, T.; Feig, M.; de Groot, B.L.; Grubmüller, H.; MacKerell, A.D. CHARMM36m: An Improved Force Field for Folded and Intrinsically Disordered Proteins. *Nat. Methods* **2017**, *14*, 71–73.

47. Zhu, F.; Tajkhorshid, E.; Schulten, K. Collective Diffusion Model for Water Permeation through Microscopic Channels. *Phys. Rev. Lett.* **2004**, *93*, 224501.

48. Mom, R.; Muries, B.; Benoit, P.; Robert-Paganin, J.; Réty, S.; Venisse, J.S.; Padua, A.; Label, P.; Auguin, D. Voltage-gating of aquaporins, a putative conserved safety mechanism during ionic stresses. *FEBS Lett.* **2020**, doi:10.1002/1873-3468.13944.

49. Michaud-Agrawal, N.; Denning, E.J.; Woolf, T.B.; Beckstein, O. MDAnalysis: A toolkit for the analysis of molecular dynamics simulations. *J. Comput. Chem.* **2011**, *32*, 2319–2327.

50. Smart, O.S.; Neduvelil, J.G.; Wang, X.; Wallace, B.A.; Sansom, M.S.P. HOLE: A Program for the Analysis of the Pore Dimensions of Ion Channel Structural Models. *J. Mol. Graph. Theor.* **1996**, *14*, 354–360.

51. Baker, N.A.; Sept, D.; Joseph, S.; Holst, M.J.; McCammon, J.A. Electrostatics of Nanosystems: Application to Microtubules and the Ribosome. *Proc. Natl. Acad. Sci. USA* **2001**, *98*, 10037–10041.
52. Dolinsky, T.J.; Czodrowski, P.; Li, H.; Nielsen, J.E.; Jensen, J.H.; Klebe, G.; Baker, N.A. PDB2PQR: Expanding and upgrading automated preparation of biomolecular structures for molecular simulations. *Nucleic Acids Res.* **2007**, *35*, W522–W525.

53. Tang, C.L.; Alexov, E.; Pyle, A.M.; Honig, B. Calculation of PKas in RNA: On the Structural Origins and Functional Roles of Protonated Nucleotides. *J. Mol. Biol.* **2007**, *366*, 1475–1496.

54. DeLano, W.L. The *PyMOL* Molecular Graphics System, Schrödinger LLC Wwwwpymol.org. Version 1. **2002** http://www.pymol.org. http://www.pymol.org.

55. Nygren, C.M.R.; Eberhardt, U.; Karlsson, M.; Parrent, J.L.; Lindahl, B.D.; Taylor, A.F.S. Growth on nitrate and occurrence of nitrate reductase-encoding genes in a phylogenetically diverse range of ectomycorrhizal fungi. *New Phytol.* **2008**, *180*, 875–889.

56. Dubey, M.K.; Ubhayasekera, W.; Sandgren, M.; Jensen, D.F.; Karlsson, M. Disruption of the Eng18B ENGase gene in the fungal biocontrol agent *Trichoderma atroviride* affects growth, conidiation and antagonistic ability. *PLoS ONE* **2012**, *7*, e36152.

57. Dubey, M.K.; Broberg, A.; Jensen, D.F.; Karlsson, M. Role of the methylcitrate cycle in growth, antagonism and induction of systemic defence responses in the fungal biocontrol agent *Trichoderma atroviride*. *Microbiology* **2013**, *159*, 2492–2500.

58. Dubey, M.K.; Broberg, A.; Sooriyaarachchi, S.; Ubhayasekera, W.; Jensen, D.F.; Karlsson, M. The glyoxylate cycle is involved in pleotropic phenotypes, antagonism and induction of plant defence responses in the fungal biocontrol agent *Trichoderma atroviride*. *Fungal Genet. Biol.* **2013**, *58–59*, 33–41.

59. Lorang, J.M.; Tuori, R.P.; Martinez, J.P.; Sawyer, T.L.; Redman, R.S.; Rollins, J.A.; Wolpert, T.J.; Johnson, K.B.; Rodriguez, R.J.; Dickman, M.B.; et al. Green fluorescent protein is lighting up fungal biology. *Appl. Environ. Microbiol.* **2001**, *67*, 1987–1994.

60. Karimi, M.; De Meyer, B.; Hilson, P. Modular cloning in plant cells. *Trends Plant Sci.* **2005**, *10*, 103–105.

61. Xu, R.; Qingshun, L.Q. Protocol: Streamline cloning of genes into binary vectors in *Agrobacterium* via the Gateway(R) TOPO vector system. *Plant Methods* **2008**, *4*, 4–10.

62. Utermark, J.; Karlovsky, P. Genetic transformation of filamentous fungi by *Agrobacterium tumefaciens*. *Nat. Protoc.* **2008**, doi:10.1038/nprot.2008.83.

63. Frąc, M. Mycological evaluation of dairy sewage sludge and its influence on functional diversity of soil microorganisms (in Polish). *Acta Agrophys. Monogr.* **2012**, *1*, 141–142.

64. Pawlik, A.; Ruminowicz-Stefaniuk, M.; Frąc, M.; Mazur, A.; Wielbo, J.; Janusz, G. The wood decay fungus *Cerrena unicolor* adjusts its metabolism to grow on various types of wood and light conditions. *PLoS ONE* **2019**, *14*, e0211744.

65. Le, S.; Josse, J.; Husson, F. FactoMineR: An R Package for Multivariate Analysis. *J. Stat. Softw.* **2008**, *25*, 1–18.

66. Kassambara, A.; Mundt, F. Factoextra: Extract and Visualize the Results of Multivariate Data Analyses. R package version 1.0.6. CRAN.R-project.org/package=factoextra. 2019. https://cran.r-project.org/web/packages/factoextra/index.html (accessed on 1 January 2021)

67. Wickham, H.; Averick, M.; Bryan, J.; Chang, W.; D'Agostino McGowan, L.; François, R.; Grolemund, G.; Hayes, A.; Henry, L.; Hester, J.; et al. Welcome to the tidyverse. *J. Open Source Softw.* **2019**, *4*, 1686.

68. Wickham, H. *ggplot2: Elegant Graphics for Data Analysis*; Springer: New York, NY, USA, 2016. ISBN 978-3-319-24277-4.

69. Auguie, B. gridExtra: Miscellaneous Functions for "Grid" Graphics. R package version 2.3. CRAN.R-project.org/package=gridExtra. 2017. https://cran.r-project.org/web/packages gridExtra/index.html (accessed on 1 January 2021)

70. Auguie, B. Egg: Extensions for 'ggplot2': Custom Geom, Custom Themes, Plot Alignment, Labelled Panels, Symmetric Scales, and Fixed Panel Size. R package version 0.4.5. CRAN.R-project.org/package=egg. 2019. https://cran.r-project.org/web/packages/egg/index.html https://rdrr.io/cran/egg/

71. McKinnon Edwards, S. lemon: Freshing Up your 'ggplot2' Plots. R package version 0.4.3. CRAN.R-project.org/package=lemon. 2019. https://cran.r-project.org/web/packages/lemon/index.html (accessed on 1 January 2021)

72. Hothorn, T.; Bretz, F.; Westfall, P. Simultaneous Inference in General Parametric Models. *Biom.* **2008**, *50*, 346–363.

73. Graves, S.; Piepho, H.P.; Selzer, L.; Dorai-Raj, S. multcompView: Visualizations of Paired Comparisons. R package version 0.1-8. CRAN.R-project.org/package=multcompView. 2019. https://cran.r-project.org/web/packages/multcompView/index.html (accessed on 1 January 2021)

74. Robinson, D.; Hayes, A. Broom: Convert Statistical Analysis Objects into Tidy Tibbles. R package version 0.5.5. CRAN.R-project.org/package=broom. 2020. https://cran.r-project.org/web/packages/RColorBrewer/index.html (accessed on 1 January 2021)

75. Neuwirth, E. RColorBrewer: ColorBrewer Palettes. R Package Version 1.1-2. 2014. Available online: https://CRAN.R-project.org/package=RColorBrewer.

76. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. *Methods* **2001**, *25*, 402–408.

77. Seidl, V.; Song, L.; Lindquist, E.; Gruber, S.; Koptchinskiy, A.; Zeilinger, S.; Schmoll, M.; Martínez, P.; Sun, J.; Grigoriev, I.; et al. Transcriptomic response of the mycoparasitic fungus *Trichoderma atroviride* to the presence of a fungal prey. *BMC Genom.* **2009**, *10*, 567.

78. Pfaffl, M.W.; Tichopad, A.; Prgomet, C.; Neuvians, T.P. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: bestKeeper-excel- based tool using pair-wise correlations. *Biotechnol. Lett.* **2004**, *26*, 509–515.

79. Abascal, F.; Irisarri, I.; Zardoya, R. Diversity and evolution of membrane intrinsic proteins. *Biochim. Biophys. Acta* **2014**, *1840*, 1468–1481.

80. Törnroth-Horsefield, S.; Wang, Y.; Hedfalk, K.; Johanson, U.; Karlsson, M.; Tajkhorshid, E.; Neutze, R.; Kjellbom, P. Structural mechanism of plant aquaporin gating. *Nature* **2006**, *439*, 688–694.
81. De Groot, B.L.; Grubmueller, H. The Dynamics and Energetics of Water Permeation and Proton Exclusion in Aquaporins. *Curr. Opin. Struct. Biol.* **2005**, *15*, 176-183.

82. Shivaraj, S.M.; Deshmukh, R.K.; Rai, R.; Bélanger, R.; Agrawal, P.K.; Dash, P.K. Genome-wide identification, characterization, and expression profile of aquaporin gene family in flax (*Linum usitatissimum*). *Sci. Rep.* **2017**, *27*, 46137.

83. Beitz, E.; Wu, B.; Holm, L.M.; Schultz, J.E.; Zeuthen, T. Point Mutations in the Aromatic/Arginine Region in Aquaporin 1 Allow Passage of Urea, Glycerol, Ammonia, and Protons. *Proc. Natl. Acad. Sci. USA* **2006**, *103*, 269-274.

84. Wang, Y.i.; Schulten, K.; Tajkhorshid, E. What Makes an Aquaporin a Glycerol Channel? A comparative study of AqpZ and GlpF. *Structure* **2005**, *13*, 1107-1118.

85. Hub, J.S.; de Groot, B.L. Mechanism of selectivity in aquaporins and aquaglyceroporins. *Proc. Natl. Acad. Sci. USA* **2008**, *105*, 1198-1203.

86. Faize, M.; Fumanal, B.; Luque, F.; Ramírez-Tejero, J.A.; Zou, Z.; Qiao, X.; Faize, L.; Gousset-Dupont, A.; Roeckel-Drevet, P.; Label, P.; et al. Genome Wild Analysis and Molecular understanding of the Aquaporin Diversity in Olive Trees (Olea Europaea L.). *Int. J. Mol. Sci.* **2020**, *21*, 4183.

87. Hub, J.S.; Aponte-Santamaría, C.; Grubmüller, H.; de Groot, B.L. Voltage-Regulated Water Flux through Aquaporin Channels In Silico. *Biophys. J.* **2010**, *99*, L97-L99.

88. De Groot, B.L.; Grubmüller, H. Water Permeation across Biological Membranes: Mechanism and Dynamics of Aquaporin-1 and GlpF. *Science* **2001**, *294*, 2353-2357.

89. Hashido, M.; Ikeguchi, M.; Kidera, A. Comparative simulations of aquaporin family: AQP1, AQPZ, AQP0 and GlpF. *FEBS Lett.* **2005**, *579*, 5549-5552.

90. Sidoux-Walter, F.; Pettersson, N.; Hohmann, S. The *Saccharomyces cerevisiae* aquaporin Aqy1 is involved in sporulation. *Proc. Natl. Acad. Sci. USA* **2004**, *101*, 17422-17427.

91. Vieira, P.M.; Zeilinger, S.; Brandão, R.S.; Vianna, G.R.; Georg, R.C.; Gruber, S.; Aragão, F.J.L.; Ulhoa, C.J. Overexpression of an aquaglyceroporin gene in Trichoderma harzianum affects stress tolerance, pathogen antagonism and Phaseolus vulgaris development. *Biol. Control* **2018**, *126*, 185-191.

92. Li, L.; Qu, Q.; Tian, B.; Zhang, K.Q. Induction of Chlamydospores in *Trichoderma harzianum* and *Gliocladium roseum* by Antifungal Compounds Produced by *Bacillus subtilis* C2. *J. Phytopathol.* **2005**, *153*, 686-693.

93. Mukherjee, P.K.; Kenerley, C.M. Regulation of Morphogenesis and Biocontrol Properties in *Trichoderma virens* by a VELVET Protein, Vel1. *Appl. Environ. Microbiol.* **2010**, *76*, 2345-2352.

94. Zachow, C.; Berg, C.; Müller, H.; Monk, J.; Berg, G. Endemic plants harbour specific *Trichoderma* communities with an exceptional potential for biocontrol of phytopathogens. *J. Biotechnol.* **2016**, *235*, 162-170.

95. Vieira, P.M.; Siqueira Guedes Coelho, A.; Steindorff, S.A.; Linhares de Siqueira, S.J.; do Nascimento Silva, R.; Ulhoa, C.J. Identification of differentially expressed genes from Trichoderma harzianum during growth on cell wall of Fusarium solani as a tool for biotechnological application. *BMC Genom.* **2013**, *14*, 177.

96. Sood, M.; Kapoor, D.; Kumar, V.; Sheteiwy, M.S.; Ramakrishnan, M.; Landi, M.; Araniti, F.; Sharma, A. *Trichoderma*: The "Secrets" of a Multitalented Biocontrol Agent. *Plants* **2020**, *9*, 762.

97. Seidl, V.; Druzhinina, I.S.; Kubicek, C.P. A screening system for carbon sources enhancing β-*N*-acetylglucosaminidase formation in *Hypocrea atroviridis* (*Trichoderma atroviride*). *Microbiology* **2006**, *152*, 2003-2012.

98. Frąc, M.; Gryta, A.; Oszust, K.; Kotowicz, N. Fast and Accurate Microplate Method (Biolog MT2) for Detection of Fusarium Fungicides Resistance/Sensitivity. *Front. Microbiol.* **2016**, *7*, 489.

99. Pertile, G.; Panek, J.; Oszust, K.; Siczek, A.; Oleszek, M.; Gryta, A.; Frąc, M. Effect of different organic waste on cellulose-degrading enzymes secreted by *Petriella setifera* in the presence of cellobiose and glucose. *Cellulose* **2019**, *26*, 7905-7922.

100. Kubicek, C.P.; Bissett, J.; Kullnig-Gradinger, C.M.; Druzhinina, I.S.; Szakacs, G. Genetic and metabolic diversity of *Trichoderma*: A case study on South-East Asian isolates. *Fungal Genet Biol.* **2003**, *38*, 310-317.

101. Zak, J.; Willig, M.; Moorhead, D.; Wildman, H. Functional diversity of microbial communities, a quantitative approach. *Soil Biol. Biochem.* **1994**, *26*, 1101-1108.

102. Garland, J. Analysis and interpretation of community-level physiological profiles in microbial ecology. *FEMS Microbiol. Ecol.* **1997**, *24*, 289-300.

103. Melton, L.D.; Smith, B.G.; Roshita, I.; Schröder, R. Mannans in primary and secondary plant cell walls. *N. Z. J. For. Sci.* **2009**, *39*, 153-160.

104. Lesage, G.; Bussey, H. Cell Wall Assembly in *Saccharomyces cerevisiae*. *Microbiol. Mol. Biol. Rev.* **2006**, *70*, 317-343.

105. Ruijter, G.J.; Bax, M.; Patel, H.; Flitter, S.J.; van de Vondervoort, P.J.; de Vries, R.P.; vanKuyk, P.A.; Visser, J. Mannitol is required for stress tolerance in *Aspergillus niger* conidiospores. *Eukaryot Cell.* **2003**, *2*, 690-698.

106. Latgé, J.P. The cell wall: A carbohydrate armour for the fungal cell. *Mol. Microbiol.* **2007**, *66*, 279-290.

107. Pail, M.; Peterbauer, T.; Seiboth, B.; Hametner, C.; Druzhinina, I.; Kubicek, C.P. The metabolic role and evolution of l-arabinitol 4-dehydrogenase of *Hypocrea jecorina*. *Eur. J. Biochem.* **2004**, *271*, 1864-1872.

108. Mach-Aigner, A.R.; Gudynaite-Savitch, L.; Mach, R.L. L-Arabitol Is the Actual Inducer of Xylanase Expression in *Hypocrea jecorina* (*Trichoderma reesei*). *Appl. Environ. Microbiol.* **2011**, *77*, 5988-5994.

109. Seiboth, B.; Herold, S.; Kubicek, C.P. Metabolic engineering of inducer formation for cellulase and hemicellulase gene expression in *Trichoderma reesei*. *Subcell Biochem.* **2012**, *64*, 367-390.
110. Trejo-Hernández, A.; Andrade-Domínguez, A.; Hernández, M.; Encarnación, S. Interspecies competition triggers virulence and mutability in Candida albicans-*Pseudomonas aeruginosa* mixed biofilms. *ISME J.* **2014**, *8*, 1974–1988.

111. Bechem, E.E.T. Utilisation of organic and inorganic nitrogen sources by *Scleroderma sinnamariense* mont. Isolated from *Gnetum africanum* welw. *Afr. J. Biotechnol.* **2012**, *11*, 9205–9213.

112. Malmberg, S.E.; Adams, C.M. Insulin signalling and the general amino acid control response. Two distinct pathways to amino acid synthesis and uptake. *J. Biol. Chem.* **2008**, *283*, 19229–19234.

113. Mallikharjuna Rao, K.L.N.; Raju, S.; Ravisankarc, H. Cultural conditions on the production of extracellular enzymes by *Trichoderma* isolates from tobacco rhizosphere. *Braz J Microbiol.* **2016**, *47*, 25–32.

114. Musa-Aziz, R.; Chen, L.M.; Pelletier, M.F.; Boron, W.F. Relative CO2/NH3 selectivities of AQP1, AQP4, AQP5, AmtB, and RhAG. *Proc. Natl. Acad. Sci. USA* **2009**, *106*, 5406–5411.

115. Heckwolf, M.; Pater, D.; Hanson, D.T.; Kaldenhoff, R. The Arabidopsis thaliana aquaporin AtPIP1:2 is a physiologically relevant CO₂ transport facilitator. *Plant J.* **2011**, *67*, 795–804.

116. Zwiazek, J.J.; Xu, H.; Tan, X.; Navarro-Ródenas, A.; Morte, A. Significance of oxygen transport through aquaporins. *Sci. Rep.* **2017**, *7*, 40411.

117. Tournaire-Roux, C.; Sutka, M.; Javot, H.; Gout, E.; Gerbeau, P.; Luu, D.T.; Bligny, R.; Maurel, C. Cytosolic pH regulates root water transport during anoxic stress through gating of aquaporins. *Nature* **2003**, *425*, 393–397.

118. Tan, X.; Xu, H.; Khan, S.; Equiza, M.A.; Lee, S.H.; Vaziriyeganeh, M.; Zwiazek, J.J. Plant water transport and aquaporins in oxygen-deprived environments. *J. Plant Physiol.* **2018**, *227*, 20–30.

119. Parrou, J.L.; Jules, M.; Beltran, G.; François, J. Acid trehalase in yeasts and filamentous fungi: Localization, regulation and physiological function. *FEMS Yeast Res.* **2005**, *5*, 503–511.

120. Jourdier, E.; Cohen, C.; Poughon, L.; Larroche, C.; Monot, F.; Ben Chaabane, F. Cellulase activity mapping of *Trichoderma reesei* cultivated in sugar mixtures under fed-batch conditions. *Biotechnol. Biofuels* **2013**, *6*, 79.

121. Yaneff, A.; Sigaut, L.; Marquez, M.; Alleva, K.; Pietrasanta, L.I.; Amodeo, G. Heteromerization of PIP aquaporins affects their intrinsic permeability. *Proc. Natl. Acad. Sci. USA* **2014**, *111*, 231–236.

122. Berny, M.C.; Gilis, D.; Rooman, M.; Chaumont, F. Single mutations in the transmembrane domains of maize plasma membrane aquaporins affect the activity of monomers within a heterotetramer. *Mol. Plant.* **2016**, *9*, 986–1003.

123. El Tarazi, A.; Lussier, Y.; Da Cal, S.; Bissonnette, P.; Bichet, D.G. Functional recovery of AQP2 recessive mutations through hetero-oligomerization with wild-type counterpart. *Sci. Rep.* **2016**, *6*, 33298.

124. Sjohann, J.; Bath, P.; Neutze, R.; Hedfalk, K. Applying bimolecular fluorescence complementation to screen and purify aquaporin protein: Protein complexes. *Protein Sci.* **2016**, *25*, 2196–2208.

125. Hohmann, S.; Krantz, M.; Nordlander, B. Yeast osmoregulation. In *Methods Enzymol*; Dieter, H., Helmut, S., Eds.; Academic Press: Cambridge, MA, USA, **2007**; Volume 428, pp. 29–45.

126. Ahmadpour, D.; Geijer, C.; Tamás, M.J.; Lindkvist-Petersson, K.; Hohmann, S. Yeast reveals unexpected roles and regulatory features of aquaporins and aquaglyceroporins. *Biochim. Biophys. Acta* **2014**, *1840*, 1482–1491.

127. Wood, J.M. Bacterial responses to osmotic challenges. *J. Gen. Physiol.* **2015**, *145*, 381.

128. Soveral, G.; Madeira, A.; Loureiro-Dias, M.C.; Moura, T.F. Membrane tension regulates water transport in yeast. *Biochim. Biophys. Acta* **2008**, *1778*, 2573–2579.

129. Ozu, M.; Dorr, R.A.; Gutiérrez, F.; Politi, M.T.; Toriano, R. Human AQP1 is a constitutively open channel that closes by a membrane-tension-mediated mechanism. *Biophys. J.* **2013**, *104*, 85–95.

130. Hill, A.E.; Shachar-Hill, Y. Are Aquaporins the Missing Transmembrane Osmosensors? *J. Membr. Biol.* **2015**, *248*, 753–765.

131. Muries, B.; Mom, R.; Benoit, P.; Brunel-Michac, N.; Cochard, H.; Drevet, P.; Petel, G.; Badel, E.; Fumanal, B.; Gousset-Dupont, A.; et al. Aquaporins and water control in drought-stressed poplar leaves: A glimpse into the extraxylem vascular territories. *Environ. Exp. Bot.* **2019**, *162*, 25–37.

132. Ozu, M.; Galizia, L.; Acuña, C.; Amodeo, G. Aquaporins: More Than Functional Monomers in a Tetrameric Arrangement. *Cells* **2018**, *7*, E209.

133. Ampah-Korsah, H.; Anderberg, H.I.; Engfors, A.; Kirscht, A.; Norden, K.; Kjellstrom, S.; Kjellbom, P.; Johanson, U. The aquaporin splice variant NbXIP1:1α is permeable to boric acid and is phosphorylated in the N-terminal domain. *Front. Plant Sci.* **2016**, *7*, 862.

134. Wallace, I.S.; Shakesby, A.J.; Hwang, J.H.; Choi, W.G.; Martínková, N.; Douglas, A.E.; Roberts, D.M. *Acyrthosiphon pisum* AQP2: A multifunctional insect aquaglyceroporin. *Biochim. Biophys. Acta* **2012**, *1*, 627–635.

135. Wysocki, R.; Chéry, C.C.; Wawrzycka, D.; Van Hulle, M.; Cornelis, R.; Thevelein, J.M.; Tamás, M.J. The glycerol channel Fps1p mediates the uptake of arsenite and antimonite in *Saccharomyces cerevisiae*. *Molec. Microbiol.* **2001**, *40*, 1391–1401.

136. Mollapour, M.; Piper, P.W. Hog1 mitogen-activated protein kinase phosphorylation targets the yeast Fps1 aquaglyceroporin for endocytosis, thereby rendering cells resistant to acetic acid. *Mol. Cell. Biol.* **2007**, *27*, 6446–6456.

137. Bienert, G.P.; Desguin, B.; Chaumont, F.; Hols, P. Channel-mediated lactic acid transport: A novel function for aquaglyceroporins in bacteria. *Biochem J.* **2013**, *454*, 559–570.

138. Daxiong, F.; Libson, A.; Miercke Larry, J.W.; Weitzman, C.; Nollert, P.; Krucin-ski, J.; Stroud, R.M. Structure of a Glycerol-Conducting Channel and the Basis for Its Selectivity. *Science* **2000**, *290*, 481–486.

139. Hohmann, S. Osmotic Stress Signaling and Osmoadaptation in Yeasts. *Microbiol. Mol. Biol. Rev.* **2002**, *66*, 300–372.
140. Tamas, M.J.; Luyten, K.; Sutherland, F.C.; Hernandez, A.; Albertyn, J.; Valadi, H.; Li, H.; Prior, B.A.; Kilian, S.G.; Ramos, J.; et al. Fps1p controls the accumulation and release of the compatible solute glycerolin yeast osmoregulation. *Mol. Microbiol.* **1999**, *31*, 1087–1104.

141. Holst, B.; Lunde, C.; Lages, F.; Oliveira, R.; Lucas, C.; Kielland-Brandt, M. GUP1 and its close homologue GUP2, encoding multimembrane-spanning proteins involved in active glycerol uptake in *Saccharomyces cerevisiae*. *Mol. Microbiol.* **2000**, *37*, 108–124.

142. Saier, M.H. A Functional-Phylogenetic Classification System for Transmembrane Solute Transporters. *Microbiol. Mol. Biol. Rev.* **2000**, *64*, 354–411.