Influence of Rhizoctonia solani and Trichoderma spp. in growth of bean (Phaseolus vulgaris L.) and in the induction of plant defense-related genes

Many Trichoderma species are well-known for their ability to promote plant growth and defense. We study how the interaction of bean plants with R. solani and/or Trichoderma affect the plants growth and the level of expression of defense-related genes. Trichoderma isolates were evaluated in vitro for their potential to antagonize R. solani. Bioassays were performed in climatic chambers and development of the plants was evaluated. The effect of Trichoderma treatment and/or R. solani infection on the expression of bean defense-related genes was analyzed by real-time PCR and the production of ergosterol and squalene was quantified. In vitro growth inhibition of R. solani was between 86 and 58%. In in vivo assays, the bean plants treated with Trichoderma harzianum T019 always had an increased size respect to control and the plants treated with this isolate did not decrease their size in presence of R. solani. The interaction of plants with R. solani and/or Trichoderma affects the level of expression of seven defense-related genes. Squalene and ergosterol production differences were found among the Trichoderma isolates, T019 showing the highest values for both compounds. T. harzianum T019 shows a positive effect on the level of resistance of bean plants to R. solani. This strain induces the expression of plant defense-related genes and produces a higher level of ergosterol, indicating its ability to grow at a higher rate in the soil, which would explain its positive effects on plant growth and defense in the presence of the pathogen.

Influence of Rhizoctonia solani and Trichoderma spp. in growth of bean (Phaseolus vulgaris L.) and in the induction of plant defense-related genes Introduction The common bean (Phaseolus vulgaris L.) is the third most important food legume crop worldwide, surpassed only by the soybean [Glycine max (L.) Merr.] and peanut (Arachis hypogea L.). Among the southern countries of the European Union, Spain together with Italy and Greece are the main common bean producers. León, a province located at the northwest of Spain, is the main producer province by quantity and quality, with almost 45% of Spanish production in 2014. Socio-economic conditions of León province enabled possible the maintenance of local varieties in traditional cropping systems, which are based in small-scale farms . The high quality of this legume has been awarded with a Protected Geographic Indication (PGI) (EC Reg. n.256/2010 published on 26 March 2010, OJEU L880/17). In the last few years, however dry bean production has gone through difficulties due to relatively low yields (mainly caused by fungus, virus, and bacteria) and insufficient income for growers.
Root rots are the main diseases caused by soil fungi having their incidence on bean yield. Rhizoctonia solani JG Kühn [Teleomorph: Thanatephorus cucumeris (AB Frank) Donk] is the main root rot in León, being detected in 91.8% of affected plants in an evaluation of its occurrence in bean plants (Valenciano et al., 2006). Plant infection occurs through wounds or by a coating of an organ with mycelium, which tears the cuticle and penetrates the epidermis. This pathogen is more aggressive at temperatures between 15 and 18 • C and in moist soils. It is a necrotrophic pathogen, distributed worldwide (Guerrero-González et al., 2011). R. solani is one of the root and hypocotyl pathogen that causes most economic losses worldwide.
Trichoderma (Teleomorph: Hypocrea) is a fungal genus that is found in the soil. It is a secondary fast growing opportunistic invasive, which produces large numbers of spores, enzymes able to degrade the fungal cell wall (chitinases, glucanases, and proteases) and compounds with antimicrobial activity. Many Trichoderma species are also well known as biocontrol agents (BCA) of important phytopathogenic fungi. The primary mechanisms of biocontrol used by Trichoderma in direct confrontation with pathogenic fungi are the mycoparasitism (Papavizas, 1985) antibiosis, and competition for nutrients with the pathogen (Harman and Kubicek, 1998).
Many Trichoderma species colonize the root surface and cause substantial changes in plant metabolism (Harman et al., 2004). The physical interaction between Trichoderma and plants is limited to the first cell layer of the epidermis and the root bark. This symbiotic relationship would thus protect plants against pathogens. Trichoderma induces the expression of genes involved in defense response and promotes plant growth, root development and nutrient availability (Hermosa et al., 2012). During the Trichoderma-plant interaction various classes of metabolites could induce resistance such as proteins with enzymatic activity, low molecular weight compounds, related to the fungal or the plant cell wall, originated by the enzymatic activity of Trichoderma (Woo et al., 2006;Woo and Lorito, 2007) and other secondary metabolites that trigger plant defense mechanisms against the pathogen (Hermosa et al., 2012;Malmierca et al., 2014), by inducing the expression of pathogenesis-related (PR) proteins that reduce the diseases symptoms. Thus, when the plant contacts with a pathogen it is activated a mechanism of systemic acquired resistance (SAR). However, when they interact with a non-pathogen organism the plants activated a mechanism for induced systemic resistance (ISR) Mukherjee et al., 2013).
Squalene is a polyunsaturated terpene that is an intermediate in the ergosterol biosynthetic pathway, which has an essential function in the fungal cell structure. The levels of squalene will influence the level of ergosterol biosynthesis (Garaiová et al., 2013). In addition to its structural function, and as a result of its importance in fungal development, ergosterol is also able to activate the expression of a number of defense genes and could increase the resistance of plants against pathogens (Lochman and Mikeš, 2005).
In this work 23 Trichoderma isolates were collected from bean fields. These isolates were used to study their effect on the growth of bean plants, and also in the defense response of plants against the phytopathogen R. solani. Thus, parameters as plant growth in the presence of the pathogen and/or the different Trichoderma isolates were evaluated, also analyzing the level of expression of defense-related genes in plants treated with the selected Trichoderma isolate.

Trichoderma and R. solani Isolates and Culture Collections
The present study was conducted with twenty-three isolates of Trichoderma (Table 1) collected from the production area of the Protected Geographical Indication (PGI), called "Alubia La Bañeza-León, " without any genetic manipulation and three isolates from other collections. The Trichoderma isolates were stored in the collection "Pathogens and Antagonists of the Laboratory Diagnosis of Pests and Diseases" (PALDPD, University of León, León, Spain). R. solani R43 was also collected from plants of the same PGI and selected by its high virulence ( Table 1).

In vitro Antifungal Assays
Trichoderma isolates were evaluated for their in vitro potential to antagonize the plant pathogenic fungus R. solani using two different tests. For all tests, plugs of 7 mm diameter collected from the edge of growing fungal colonies were used to inoculate potato dextrose agar medium (PDA) in sterile Petri dishes of 9 cm diameter. The dishes were incubated in the dark at 22 • C for 7 days.
The aim of these tests was to study the percentage of R. solani growth inhibition caused by the different Trichoderma isolates.
The antifungal assay on membranes was used to quantify the ability of the Trichoderma isolates to produce metabolites and/or enzymes with inhibitory activity against R. solani. The surface of Petri dishes containing PDA medium was overlaid with a sterile cellophane membrane. Trichoderma plugs, extracted from PDA dishes grown for 7 days at 22 • C, were placed in the center of the dish with the cellophane sheet, containing PDA medium, and incubated for 48 h at 22 • C. Then, the cellophane membranes along with the mycelia of Trichoderma isolates were removed and R. solani plugs were placed in the same plates. Growth of R. solani was recorded after 72 h to calculate the percentage of pathogen growth inhibition ( Figure 1A). Control PDA plates of R. solani, where Trichoderma spp. had not been previously grown, were also prepared in the same conditions as above.  The percentage of inhibition (IM) was calculated after 3 days of growth of R. solani in this medium using the formula %IM = [(C-T)/C]x100 (C: diameter of the R. solani control; T: diameter of R. solani after being exposed to the metabolites of Trichoderma spp.). Experiments were performed with four replicates. The results were compared by analysis of variance (ANOVA) and Fisher least significant difference (LSD) tests using SAS (SAS Institute Inc., 2004, Cary, NC, USA). The direct confrontation assays were used to verify the ability of Trichoderma spp. to overgrow the pathogen. Each Trichoderma isolate was grown in dual culture with R. solani R43. The isolates were placed 5.5 cm apart on the same plate and incubated at 22 • C for 5 days. Experiments were performed with four replicates. The parameters were measured after 5 days: r1 (distance between the pathogen sowing point and furthest point of the colony) and r2 (distance between the pathogen sowing point and the edge of the colony) from where R. solani and Trichoderma mycelia came into contact. Thus, the percentage of inhibition in the direct confrontation assay (ID) was calculated by the formula: %ID =[(r1-r2)/r1]x100 ( Figure 1B). Inhibition of R. solani growth was compared by analysis of variance (ANOVA) and Fisher least significant difference (LSD) tests using SAS. (SAS Institute Inc., 2004, Cary, NC, USA).

In vivo Assay of the Antifungal Activity
These assays were only performed with those Trichoderma isolates that gave percentages of inhibition greater that 40% in membrane assays or 20% in direct confrontation assays, and that were able to sporulate on PDA medium.
The bioassays were performed in climatic chambers with 32 treatments as follow: 15 in vitro selected Trichoderma isolates against R. solani (R43) in order to test the antagonistic activity (RT0-number of Trichoderma isolate); 15 in vitro selected Trichoderma isolates in order to test their effect on plant (CT0number of Trichoderma isolate); one control with R. solani (RC); and one control without fungi (CC). Thirty pots were used per treatment, with polypropylene pots (1 liter capacity) with substrate (80% white peat, 20% black peat and 5.5 pH). Each pot was watered with 250 ml of water prior inoculation. R. solani R43 was inoculated by surface irrigation with 50 ml per pot of a suspension of triturated micromicete culture of this pathogen using five Petri dishes (18 ml of PDA per dish) per liter of water. For control inoculation, only PDA medium was used without any pathogen. Pots were kept in a growth chamber for 8 days at 25 • C (16 h) and 16 • C (8 h), 60% relative humidity (RH) in the dark.
Trichoderma isolates were inoculated on PDA medium to grow in dark conditions (25 • C) for 1 week. After that, they were exposed to light in order to induce the formation of spores. Spores suspensions were prepared at a final concentration of 2 × 10 7 spores/ml. Bean seeds of "Canela" variety were surface sterilized (sodium hypochlorite 1% for 3 min and distilled water for 6 min). Then, they were coated with a spore suspension of each Trichoderma isolate. The seeds were submerged in the spore suspension (45 seeds per 20 ml spores suspension) and they were dried in a flow chamber for 12 h. Coated seeds were sown after 8 days of the inoculation of R. solani R43. The culture was maintained for 45 days with a photoperiod of 16 h light, 25 • C/16 • C (day/night), 60% RH and brightness of 3500 lux. Irrigations were performed every 4 days with tap water (about 250 ml/pot). On the 2nd-4th week a nutrient solution was added (Rigaud and Puppo, 1975). Plants were removed after 45 days from sowing, tissues with symptoms were placed in PDA medium, incubating the plates at 22 • C for 5 days and identifying the fungus for fulfill Koch's postulates. The next parameters were evaluated in removed plants after 45 days from sowing: wet weight and dry weight (72 h in an oven, 82 • C) of the aerial part and root system. The data were transformed by the formula √ x + 0.5 and they were compared by analysis of variance (ANOVA) and Fisher least significant difference (LSD) tests using SAS (SAS Institute Inc., 2004, Cary, NC, USA).

Nucleic Acid Extraction and Manipulation
Genomic DNA from those isolates showing positive and negative phenotypic effects on bean plants, were extracted by growing the Trichoderma isolates in PDB medium (potato dextrose broth). Mycelia were then recovered by filtration, washed with 0.9% NaCl and dried on absorbent filter paper. The procedure for fungal genomic DNA isolation was performed as previously described (Cardoza et al., 2006).

PCR Amplification, Sequencing and DNA Analysis
The amplification of the ITS regions of the nuclear rDNA gene cluster and an approximately 0.56 kb fragment of the tef1 (translation elongation factor 1-α) gene were carried out with the primer pairs ITS1/ITS4 and EF1-728F/EF1-LLErev, respectively, as described previously (Hermosa et al., 2004). The PCR products were purified from agarose gels using the NucleoSpin Extract II Kit (Macherey-Nagel, Düren, Germany), according to the manufacturer's protocol. PCR fragments were sequenced in an ABI 377 Prism Sequencer (Applied Biosystems, Foster City, CA). The Trichoderma ITS and tef1 sequences obtained in this work were analyzed using the online interactive key (available from http://www.isth.info/tools/blast/index.php) (Druzhinina et al., 2005).

Analysis of Expression of Bean Defense-related Genes
Three bean leaves from 45 day-old plants of each treatment were randomly collected and stored at −80 • C until use. Leaves were detached from plants inoculated with Trichoderma isolate showing positive phenotypic results in the in vivo test. Leaves were then reduced to a fine powder in a mortar under liquid nitrogen. Plant RNA isolation were performed as previously described (Malmierca et al., 2013).
cDNA were synthesized using 1 µg total RNA and a Reverse Transcription System with an Oligo(dT) 15 as the primer (Promega, Madison, WT). cDNA were quantified using a Nanodrop 2000 (Thermo Scientific, Wilmington, DE) and used for further studies.

Real Time-PCR Analysis
In order to analyze the effect of Trichoderma treatment and/or R. solani infection of bean plants, oligonucleotides corresponding to seven defense-related genes were designed based on their available sequences (Table 2). PR1, PR2, PR3, and PR4, which encode for pathogenesis related proteins related to the salicylate (SA) pathway; CH5b, CH1 encoding for related to the jasmonate/ethylene pathway (JA/ET), and PAL involved in the phenylpropanoid pathway, were selected to be analyzed in the present study. α-actin (Upchurch and Ramirez, 2010;Guerrero-González et al., 2011) and PvEF1α (this work) encoding genes were used as reference (housekeeping genes) for comparative analysis. The qPCR reactions were carried out using Step One Plus ™ (Applied Biosystems, Foster City, CA). The reactions were performed in a total volume of 20 µl: 10 µl Power SYBR R Green PCR Master Mix (Applied Biosystems, USA), 0.4 µl Forward Primer 10 µM, 0.4 µl Reverse Primer 10 µM, 5 µl cDNA, and H 2 O to 20 µl. The REST 2009 © software (Pfaffl et al., 2002) was used to calculate the relative expression ratio and the significance of the differences between the gene expression levels. For each primer pair used in this work, we performed a standard curve with 320, 160, 80, 40, 20, and 10 ng cDNA to determine the PCR amplification efficiency (E value). Each measurement was made in triplicate (Malmierca et al., 2013).

Quantification of Ergosterol and Squalene
The Trichoderma selected strain was inoculated in 100 ml of CM medium (0.5% malt extract, 0.5% yeast extract, and 0.5% glucose) with 10 6 spores/ml, and incubated 24 h at 28 • C. Then 20 ml from the previous cultures were inoculated on 100 ml potato dextrose broth (PDB medium) and were incubated as before during 24-96 h. The mycelia were filtered through nytal filters (30 µm diameter) and the liquid removed by drying between filter papers. The dry weight of the fungal pellet was calculated. Total intracellular sterols were extracted and ergosterol and squalene content were quantified as previously reported other authors (Cardoza et al., 2007;Ghimire et al., 2009). All measurements were made in duplicate in the Trichoderma selected isolate and a strain used as control of the same species of the selected isolate. The results were compared by analysis of variance (ANOVA) and Fisher least significant difference (LSD) tests using SAS (SAS Institute Inc., 2004, Cary, NC, USA).

Analysis of the in vivo Antagonistic Activity of Trichoderma Isolates with R. solani
The results of the in vitro membrane assays and direct confrontation assays against R. solani, showed above, were used to select the isolates that would be used for the in vivo analysis. Thus, fifteen Trichoderma isolates (marked with asterisk in Tables 1, 3) were tested since they were able to sporulate in PDA medium and showed a percentage of inhibition of R. solani growth higher than 40% in the membrane assays, and/or 20% in the direct confrontation assays.
On plants removed after 45 days from sowing for Koch's postulates, R. solani was found and Trichoderma isolates were also present in the medium.
When dry aerial parts (Figure 2) were analyzed, plants treated with CT019 had the greatest weight of the aerial part, being significantly different from the control (CC). In the RT019 treatment, plants did not show significant differences in comparison with the control (CC), what is indicative of a biocontrol effect by T019 isolate.
In the case of the root system (Figure 3), the situation was similar to that observed in the aerial parts. CT019 treated plants were not significantly different in dry weight, although it was always greater than the weight of the control plants (CC). If the pathogen was present in the soil with same Trichoderma isolate, RT019 treatment, in the case of wet weight, there were no significant differences between control (CC) and control pathogen (RC) (data not shown).
Based on these results, the Trichoderma isolate T019 was selected for further studies since it showed the best positive effects on plant phenotype among all the analyzed isolates (Figure 4).

Molecular Identification of Trichoderma Isolates
Those Trichoderma isolates able to sporulate in PDA medium, showing a percentage of inhibition higher than (1) Growth assay on membranes (effect of Trichoderma isolates on percent growth inhibition of R. solani. using cellophane membranes. (2) in direct confrontation assays (percentage of inhibition of R43 growth when grown confronted with Trichoderma strains during 5 days in PDA medium plates). *Trichoderma isolates selected for in vivo experiments because they were able to sporulate in PDA medium and showed a percentage of inhibition higher than 40% in the membrane assays. and/or 20% in the direct confrontation assays.
40% in the membrane assays, and/or 20% in the direct confrontation assays, were identified. ITS1 region of rDNA and a fragment of the translation elongation factor 1 (tef1) were amplified and sequenced for 12 isolates, and both nucleotide sequences were used for identification at species level. Eight out of these 12 isolates were identified as T. harzianum, and T. atroviride, T. gamsii, T. longibrachiatum and T. citrinoviride species were represented by one isolate ( Table 1).

Effect of Trichoderma Treatment and/or R. solani Infection in the Expression of Bean Defense Related Genes
The amplification efficiencies of the oligo-pairs ( α-actin and PvEF1α were used as housekeeping genes to determine the relative expression level of the other genes analyzed in the present work. Trichoderma T019 strain was selected, based on its positive effects on bean phenotype with and without R. solani infection. The results included in Figure 5 showed that: (i) R. solani down-regulated the expression of all the P. vulgaris defense-related genes analyzed (Figure 5A), raising values of expression ranging between 0.099 (p = 0.029) for PR2 and 0.397 (p = 0.045) for CH1; (ii) In plants treated with T019 compared with control plants (bean plants not treated with Trichoderma nor infected with R. solani) only the CH5b and PR2 were significantly up-regulated, raising comparative expression values of 1.495 (p = 0.000) and 24.492 (p = 0.000), respectively ( Figure 5B). Finally, (iii) treatment with T019 and infection with R. solani significantly up-regulated the expression ratio of the analyzed genes, except PAL, in comparison with plants only infected with R. solani, with values ranging from 1.420 (p = 0.000) to 42.975 (p = 0.000) for CH1 and PR4 respectively ( Figure 5C).
A non-remarkable significant up-regulation of the defenserelated gene expression was observed when expression ratios of the selected genes in plants treated with T. harzianum T019 and infected with R. solani were compared with those from plants untreated or uninfected (data not shown).

Production of Ergosterol and Squalene by the Trichoderma Selected Strains
Strain T019 (T. harzianum) was selected for this analysis because it was those producing the most remarkable positive phenotypic effect on beans (see above). T. harzianum T34 ( Table 1) was used as control strain for comparative purposes because it is a well-known strain, widely characterized (Kullnig et al., 2001).
Production of ergosterol at 24 h of growth did not show significant differences between the strains T. harzianum (T34 and T019) ( Table 4). When the production of ergosterol was analyzed at 96 h, T019 produced significantly higher amounts than its respective control (T34).
Regarding the squalene production at 24 h, T. harzianum (T34 and T019) strains did not show significant differences. However, at 96 h of growth, T019 produced amounts significantly higher than its control.

Discussion
Most of the Trichoderma isolates studied in the present work have been able to control the growth of R. solani. The percentage of growth inhibition in the direct confrontation assays raised values up to 72.77% for T021, but most of the isolates were between 47 and 24%. By contrast, in the antifungal assays on membranes, those percentages ranged mostly between 86 and 58%. In the assays conducted by Campelo et al. (2010) the percentage of inhibition in membrane assays using T. virens T59 (NBT59) and T. atroviride T11 (IMI352941) range between 100 and 84.7%, respectively. The differences observed in the different in vitro assays might be due to the variability of genotypes, with differences in growth, sporulation and in their environmental  biological interactions as a consequence of the prevalence of different modes of action on each of the analyzed strains (Ruano-Rosa et al., 2010).
The in vitro antagonistic activity of the Trichoderma isolates is an indication of their in vivo biological activity against phytopathogenic fungi. However, it has been found that the in vitro antifungal activity of Trichoderma does not correlate in a direct way with its in vivo ability to control diseases caused by phytopathogenic fungi, since many other factors influence that activity (Anees et al., 2010).
Regarding the in vivo assays, bean plants that had been in contact with some Trichoderma isolate had an increased size when the pathogen was not present. Pereira et al. (2014) also observed that T. harzianum was able to promote the growth of common bean plants in comparison to plants grown in its absence. However, in the present work only T019 treated plants did not decrease their size in the presence of R. solani. These data are in agreement with those from Yedidia and coworkers reporting a much stronger effect on cucumber plants treated with T. harzianum, which increased by 75% the length of the root, 95% aerial parts, 80% dry weight and 80% the size of the blade relative to the untreated control (Yedidia et al., 2001). T. harzianum was able to promote the growth of common bean plants in comparison to plants grown in its absence (Pereira et al., 2014).
However, in the present work opposite results were observed to those previously reported (Tello et al., 1985), indicating that infection with R. solani not always resulted in a reduction of the bean plant size, which was explained as a result of the activation of the natural plant-defense mechanisms (Cardoso and Echandi, 1987) that would lead to an enhanced development of the plants when this pathogen was present in the soil. In the present study we have observed that treatment with R. solani resulted in a down-regulation of all the bean defense-related genes analyzed, indicating a certain compensation effect between the intensity of plant defense response and plant growth . Little is known about the effect of Trichoderma treatment and/or infection with R. solani in the expression of bean defenserelated gene expression. As indicated above, interaction of the plant with R. solani during 45 days of growth caused the repression of the seven defense-related genes studied (CH5b, CH1, PR1, PR2, PR3, PR4, and PAL1), as a mechanism to overcome the plant defense response and thus facilitating the progression of the infection process within the plant. Other authors have observed similar results in assays with tomato plants infected with Pseudomonas syringae pv. tomato, where a repression of PR1 and PR4 was shown, as an indication of the decrease in the plant self-defense mechanism, then facilitating disease progression (Zhao et al., 2003).
Expression of P. vulgaris defense-related genes was analyzed in leaves, even when the interaction with the pathogen R. solani is initially produced at the root level, to determine if the signals generated in roots as result of this interaction are able to systematically stimulate the bean defense along long distance from seed to leaf. When the interaction of plant with T019 was analyzed, PR2 was up-regulated, PR4 was slightly restrained but PR1 and PR3 were not affected. Pereira and coworkers observed that T. harzianum also seems to potentiate common bean response against the phytopathogenic fungus R. solani, as shown by the increase in the levels of glu1 and pod3 in the double treatment in comparison to that obtained for plants in the presence of R. solani alone (Pereira et al., 2014). Other studies showed that there was an increase in expression PR1 at 16 h of interaction (Guerrero-González et al., 2011). Rivière and coworkers suggested that PR1 was down-regulated by the β-1,3glucanases (PR2 and PR3) (Rivière et al., 2008). In the present study, after 45 days of the inoculation of the biocontrol isolate we observed a higher expression of PR2, which could be due to an alteration of plant defense responses at these longer times in comparison with shorter ones assayed in the previous reported works. This would result in a higher β-1,3-glucanase activity in the cell wall that would increase oligosaccharides released by its action, then acting as elicitors of the plant defense response and/or of the fungal secondary metabolism (Druzhinina et al., 2011).
In the case of plants inoculated with Trichoderma and seeded in soil infected with R. solani, when compared with plants only infected with R. solani, thus excluding the effect due exclusively to this pathogen, a significant up-regulation of all the analyzed genes was observed, except PAL, indicating that in presence of R. solani several mechanisms are induced in Trichoderma that potentiate its ability to elicit plant defense-responses. In previous works, during the interaction of Solanum tuberosum with R. solani it was observed that after 120 h from the initial infection and damage to the first floor, new outbreaks produce less symptoms, suggesting that the plant might be prepared to defend and avoid disease progression (Lehtonen et al., 2008). Other studies have shown that the PR4 was overexpressed in eggplants treated with BCAs suggesting that these agents promoted defensive reactions within the plant (Angelopoulou et al., 2014). In other assays with pepper plants an increase in the accumulation of PR1 and chitinases was also observed when these plants were previously inoculated with a strain of Fusarium oxysporum and subsequently exposed to Verticillium dahliae (Veloso and Díaz, 2012), emphasizing their importance in the response to pathogen infection and to abiotic stresses (Jung and Hwang, 2000). Thus, it could be inferred that the effect of Trichoderma would therefore modulate the response of the plant and prevent the suppression of defense genes caused by R. solani. However, the effect of the different Trichoderma strains in the different plants would indicate the existence of great differences between them.
In the case of CH5b, related to JA pathway, when R. solani was in contact with the plant for 45 days of development, the expression of this gene was down-regulated in comparison with non-infected plants. When bean plants were inoculated with Trichoderma, its expression was slightly induced. However, if the antagonist and the pathogen were present in the medium, the effect caused by R. solani was the overexpression of this gene. In previous works (Broglie et al., 1991) it was shown that CH5b, responsible for chitinase production, was overexpressed in Nicotiana tabacum and Brassica napus infected with R. solani, and their disease symptoms were reduced. Benhamou et al. (1993) found that the expression of this gene caused a reduction of disease in B. napus infected with the same pathogen, which correlates with similar results observed in strawberry plants inoculated with Botrytis cinerea (Vellicce et al., 2006). In the present study we observed that the contact of the plant with the pathogen during a period of 45 days resulted in a significant repression of this gene, similarly to what it was observed for the expression of CH1, which thereby facilitates disease progression. However, in the presence of Trichoderma a similar effect to that previously described for SA-related genes was observed, so facilitating a prevention of disease progression.
In the case of the phenylalanine ammonia lyase (PAL), this is a key enzyme in the metabolic pathway of phenylpropanoid compounds by catalyzing the amino acid L-phenylalanine  deamination, giving rise to trans-cinnamic acid and ammonia. The trans-cinnamic acid is used for the synthesis of various phenolic compounds, which are precursors in the synthesis of esters, coumarins, flavonoids, and lignin. The production of this enzyme is controlled during plant growth, but is also induced in cells neighboring the infection site and, besides infection, various environmental stimuli such as injury, heavy metal contamination, light and growth regulators (Rahman and Punja, 2005). In the present work, we observed that this gene was down-regulated or not affected in all the conditions assayed, indicating that in the analyzed strain, PAL activity would be more focused in the plant response against environmental stress than in the response against pathogens, which contrast with the role of this enzyme described in pepper, where it is involved in the positive regulation of SA-dependent defense signaling (Kim and Hwang, 2014). Regarding squalene and ergosterol production, in this study it was found that there were differences among the isolates. Ergosterol is a sterol found in the fungal membrane, which although considered by the plant as a PAMP (pathogen associated molecular patterns) (Nürnberger et al., 2004) triggers a series of reactions (Cervone et al., 1997) which would cause an activation of the genes of plant defense (Rossard et al., 2010). Squalene is a precursor of the ergosterol biosynthetic pathway (Malmierca et al., 2013) located in the cellular membranes or accumulated as droplets in the cytoplasm whose physiological function, apart from acting as an ergosterol precursor, still remains unclear, but it might play some role in the elicitation of plant defense responses. Thus, in this assay, an increased production of these compounds by Trichoderma would result in the induction of defense genes in the bean, then the plant could grow better under a pathogen presence in the soil.

Conclusions
Trichoderma isolates inhibit the development of R. solani under in vitro conditions, when they grow in Petri dishes. In vivo conditions, the bean plants that had been in contact with the Trichoderma isolates always had an increased size when the pathogen was not present. When R. solani was present in the soil, the development of the bean plants was significantly reduced. However, bean plants treated with strain T019 did not decrease their size in the presence of R. solani.
The interaction of bean plants with R. solani caused, after 45 days of growth, the down-regulation of the seven defenserelated genes studied (CH5b, CH1, PR1, PR2, PR3, PR4, PAL) as a mechanism to overcome the plant defense response and thus facilitating the progression of the infection process within the plant.
An increased production of ergosterol and squalene by Trichoderma resulted in the induction of defense genes in the bean plants. In this way, the plant would grow better under a pathogen presence in the soil.
T. harzianum T019 shows a positive effect on the level of resistance of bean plants to R. solani. This strain induces the expression of plant defense-related genes and produces a higher level of ergosterol, indicating its ability to grow at a higher rate in the soil, which would explain its positive effects on plant growth and defense in the presence of the pathogen.