Next Article in Journal
Genetic Ancestry Inference and Its Application for the Genetic Mapping of Human Diseases
Next Article in Special Issue
Role of Thylakoid Protein Phosphorylation in Energy-Dependent Quenching of Chlorophyll Fluorescence in Rice Plants
Previous Article in Journal
Level of Murine DDX3 RNA Helicase Determines Phenotype Changes of Hepatocytes In Vitro and In Vivo
Previous Article in Special Issue
PSII Activity Was Inhibited at Flowering Stage with Developing Black Bracts of Oat
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 13
10.3390/ijms22136961
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Trichoderma harzianum on Photosynthetic Characteristics and Fruit Quality of Tomato Plants

by
Igor D. Vukelić
1,
Ljiljana T. Prokić
2,
Gordana M. Racić
1,
Mirjana B. Pešić
2,
Mirjana M. Bojović
1,
Edyta M. Sierka
3,
Hazem M. Kalaji
4,5 and
Dejana M. Panković
1,*
1
Faculty of Ecological Agriculture, Educons University, Vojvode Putnika 87, 21208 Sremska Kamenica, Serbia
2
Faculty of Agriculture, University of Belgrade, Nemanjina 6, 11080 Belgrade, Serbia
3
Faculty of Natural Sciences, Institute of Biology, Biotechnology and Environmental Protection, University of Silesia in Katowice, 28 Jagiellonska, 40-032 Katowice, Poland
4
Department of Plant Physiology, Institute of Biology, Warsaw University of Life Sciences—SGGW, 02-787 Warsaw, Poland
5
Institute of Technology and Life Sciences, National Research Institute, Falenty, Al. Hrabska 3, 05-090 Raszyn, Poland
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(13), 6961; https://doi.org/10.3390/ijms22136961
Submission received: 27 April 2021 / Revised: 20 June 2021 / Accepted: 22 June 2021 / Published: 28 June 2021
(This article belongs to the Collection Photosynthesis in Plants and Algae: Physiological, Biochemical, and Molecular Perspectives)
Browse Figures
Versions Notes

Abstract

:
The beneficial role of fungi from the Trichoderma genus and its secondary metabolites in promoting plant growth, uptake and use efficiency of macronutrients and oligo/micro-nutrients, activation of plant secondary metabolism and plant protection from diseases makes it interesting for application in environmentally friendly agriculture. However, the literature data on the effect of Trichoderma inoculation on tomato fruit quality is scarce. Commercially used tomato cultivars were chosen in combination with indigenous Trichodrema species previously characterized on molecular and biochemical level, to investigate the effect of Trichoderma on photosynthetic characteristics and fruit quality of plants grown in organic system of production. Examined cultivars differed in the majority of examined parameters. Response of cultivar Gružanski zlatni to Trichoderma application was more significant. As a consequence of increased epidermal flavonols and decreased chlorophyll, the nitrogen balance index in leaves has decreased, indicating a shift from primary to secondary metabolism. The quality of its fruit was altered in the sense of increased total flavonoids content, decreased starch, increased Bioaccumulation Index (BI) for Fe and Cr, and decreased BI for heavy metals Ni and Pb. Higher expression of swolenin gene in tomato roots of more responsive tomato cultivar indicates better root colonization, which correlates with observed positive effects of Trichodrema.

1. Introduction

One of the main goals of modern sustainable agriculture is to support raising demands for food and energy of a growing world population but at the same time to maintain soil health and fertility. Soil microorganisms, soil type, agricultural management and plant genotype play crucial roles in determination of soil functioning and fertility [1].
With a production of about 180 million tons in 2017 (FAOSTAT), tomato (Solanum lycopersicum L.) is one of the most important cultivated vegetable in the world. Tomato fruit is a source of vitamins, carotenoids, carbohydrates, phenolic compounds and nutrients which are of vital role in human diet. Among many phenolic compounds in tomato fruit, flavonoids are contributing to its anti-oxidative, anti-cancer, anti-diabetes and cardiovascular protective effects [2]. Carbohydrate content in tomato fruit is considered to be strongly related to fruit yield and quality, as it has important role in fruit growth, composition and ripening [3]. Containing only about 4% of carbohydrate, it is considered as non-starchy vegetable and advised to be used in weight-loss diets [4].
Many factors affect tomato fruit quality: plant nutrients, climatic factors, physical and chemical soil properties and fertilization [5,6]. The Bioaccumulation Index (BI) is often used to investigate the plant efficiency for uptake of elements from the soil and accumulate it in the edible part of the plant [7]. In general, plant symbionts are affecting the uptake of nutrients from the soil [8]; however, their influence on plant metal uptake are still poorly understood [9]. Trichoderma fungi are rhizosphere microorganisms often in symbiotic interaction with plants. The literature data on the effect of Trichoderma inoculation to the plant uptake of nutrients are not consistent. In some cases, the uptake was increased [10] and in others decreased [11]. The plant uptake of nutrients as influenced by Trichoderma is considered to be a key factor for many positive effects of this fungi to plant status, even more so nutrient availability in soils influences the biocontrol efficacy of Trichoderma species [12].
The beneficial properties of Trichoderma and its secondary metabolites on plants concerning improved plant growth parameters, uptake of macronutrients and oligo/micro-nutrients, water use efficiency, photosynthesis, activation of plant secondary metabolism and accumulation of polyphenols have been observed [13,14]. Increased hormone accumulation, altered sugar partitioning, enhanced photosynthetic efficiency and growth promotion of tomato plants were induced by specific strain of T. harzianum [15]. Improved photosynthetic capabilities of different plant species induced by various endophytic strains of Trichoderma due to the increase of photosynthetic pigments or gene expressions regulating chlorophyll biosynthesis, light harvesting complex proteins or Calvin cycle components were documented in more than twenty papers [16]. The upregulation of plant genes or proteins providing optimized internal redox environment in response to abiotic of biotic stresses when roots are colonized by Trichoderma strains was even more significant. However, information about its effect on tomato fruit quality in organic cultivation is scarce. To our knowledge, only two studies on the effect of different Trichoderma species on tomato fruit quality were conducted [13,17]. In the first one, T. longibrachiatum significantly promoted the growth and yield of tomato, obtaining fruits with higher titratable acidity and lower total soluble solids (TSS) content [17]. In the second one the total yield, contents of lycopene and some amino acids increased in the plum tomato fruits after the application of T. harzianum, which was discussed in terms of improved nitrogen metabolism [13].
The recognition between Trichoderma and its host during symbiotic associations is mediated by effectors [18,19,20]. Trichoderma colonizes root surfaces and penetrates to the first or second layers of cells, which is achieved by secretion of cellulolytic and proteolytic enzymes, for example, swolenin that can disrupt the crystalline cellulose structure of plant cell walls [21,22].
The aim of this paper was to study the long-term effects of T. harzianum application on two tomato cultivars (Narvik and Gružanski zlatni) grown in field conditions under polytunnel. For that purpose, we have selected strain T. harzianum SMZC 22660, as it previously exhibited positive effects on the growth of tomato plants, content of epidermal flavonols and decreased the uptake of toxic elements (Cd, Ni and Cr) in a growth chamber experiment [23]. In addition, all examined extracellular enzyme activities of this strain were high, in particular plant cell wall-degrading enzymes which are important for good root colonization. The effects on photosynthetic characteristics were examined by non-destructive measurements on plant leaves. Moreover, fruit was examined for morphological characteristics, and in order to evaluate fruit quality, total phenolic, total flavonoid, soluble sugars and starch content, and uptake efficiency of nutrients and heavy metals were measured.

2. Results

Prior to application, T. harzianum SMZC 22660 was characterized by measurement of extracellular enzyme activities (Table 1). High activities of N-acetyl-β-glucosaminidase, acid phosphatase and naphtol-AS-BI-phosphohydrolase were observed (classified according to the color scale), while low activities were determined for alkaline phosphatase, esterase C4, esterase lipase C8. The suspension of T. harzianum spores (8 × 106 CFU/mL) was applied in the root zone of tomato plants at the beginning of flowering phase and the measurement were done 60 days after application. The plants that were grown in the polytunnel were optimally irrigated by the drip system as confirmed by the SWC and RWC of the leaves (Table 1).
The colonization of the roots of two examined tomato varieties was checked by the expression analysis of swolenin, one of the Trichoderma effector proteins. The expression was lower in the roots of Narvik (lanes 1 and 2 represent two replicate bulk samples) than in Gružanski zlatni (lanes 3 and 4 represent two replicate bulk samples) (Figure 1).
Several morpho-physiological parameters of plants and fruits as well as some fruit quality parameters were measured.
The content of chlorophyll was higher in cultivar Gružanski zlatni in general. However, Trichoderma treated plants exhibited 13% lower chlorophyll content, which was not the case in cultivar Narvik. The response of the content of flavonols was opposite in examined cultivars: in Narvik it decreased by 12%, while in Gružanski zlatni it increased by 25%, in response to Trichoderma treatment. Consequently, NBI has decreased by 21% with treatment in Gružanski zlatni. Generally, NBI was lower in Narvik than in Gružanski zlatni. In both cultivars the content of anthocyanins was low and did not change after Trichoderma treatment.
According to most of the examined morphological parameters of tomato plants and fruits, there were no significant differences neither between cultivars, nor between Trichoderma treatments and controls (Table 2). The most significant difference between the examined cultivars was the number of fruits, almost doubled in cultivar Gružanski zlatni, but it did not respond to Trichoderma treatment.
Content of total polyphenols in the fruit significantly decreased to the similar extent in both investigated genotypes with addition of Trichoderma, 13.7% for Narvik and 14.5% for Gružanski Zlatni (Table 3). Similarly, significant decrease (26%) of flavonoid content was observed when Trichoderma was applied in genotype Narvik; however, flavonoids increased by 48.9% in Gružanski zlatni. The decrease of starch content (37%) after Trichoderma treatment was observed in the same cultivar. Significantly lower content of soluble sugars was observed in tomato fruits of cultivar Narvik in comparison to Gružanski zlatni; however, there were no significant changes with Trichoderma treatment.
Nitrogen content in Cultivar Narvik was 1.28 ± 0.03% and it did not change with Trichoderma presence. However, in cultivar Gružanski zlatni content was 2.21 ± 0.18% and it decreased with Trichoderma treatment to 1.59 ± 0.08%. Content of phosphorus did not change significantly over the two cultivars and treatments; it ranged from the lowest 0.715 ± 0.045 ppm for Gružanski zlatni to the highest 0.79 ± 0.02 ppm for Narvik. Similarly, potassium content did not differ significantly either between treatments or between genotypes: it ranged from the lowest 2.78 ± 0.28 ppm in Gružanski zlatni to the highest 2.98 ± 0.14 ppm in Narvik (Table 4).
The reaction of two cultivars to the fungal treatment was the opposite in case of Mn. In cultivar Narvik it decreased, while in Gružanski zlatni it increased, but the changes were less than 10% (Table 4). The concentration of Fe and Cr were significantly higher in the Trichoderma-treated plants in comparison to the control in both cultivars. The changes of Fe and Cr were lower in Narvik, 23% and 59%, in comparison to Gružanski zlatni 72% and 100%, respectively. On the other hand, concentrations of Zn, Cu and Co were not affected by Trichoderma treatment and were similar in two cultivars as well.
Concentrations of Ni has decreased with the treatment in fruit of both cultivars, 20% in Narvik and 50% in Gružanski zlatni. However, concentrations of Pb responded significantly to the treatment only in Narvik, where it decreased by 50%.
BI for Mn, Ni and Cr ranged from 1.43 to 5.6. The lowest determined BI was for Fe, Pb and Co (˂1). However, high values of BI for Cu and Zn, ranging from 18.5–33, have been calculated. Values of BI in control conditions were similar in two cultivars. However, cultivar Gružanski zlatni had higher BI for Fe and Cr, and lower BI for Ni in response to Trichoderma treatment. Both cultivars responded to the treatment by lowering BI for Pb (Table 5).

Principal Component Analysis

The effect of Trichoderma application on two tomato cultivars was visualized by performing PCA on parameters measured on plants (Figure 2) and fruits separately (Figure 3). The PCA of plant parameters, explained in total 78% of the variability (PC1 51.8%; PC2 26.2%). PC1 clearly separated two examined cultivars (Figure 2). Among the tested parameters, five out of seven were significantly correlated with the first axis. These parameters include NBI, Chl, Flav, Anth and STB. On the other hand, the RWC, Flav, STB and STM were significantly correlated with the second axis, which separated all treatments except one control sample of cultivar Narvik, which was in the group of samples from plants treated with Trichoderma.
The PCA of fruit parameters, explained in total 62.8% of the variability (PC1 40%; PC2 22.8%). Again, PC1 clearly separated two examined cultivars (Figure 3). Among the tested parameters, ten out of nineteen were significantly correlated with the first axis: contents of Pb, P, N, Mn, Fe, Cu, Zn, TP, SS, FN. Meanwhile, the second axis was determined by six out of nineteen parameters: contents of Ni, Cr, Co, S, FH, FM.

3. Discussion

It is known that some selected strains of Trichoderma promote plant growth and development, increase nutrient uptake and finally increase yields [24,25,26,27]. The efficient root colonization of this symbiont depends on the secreted effector molecules. About 20 Trichoderma effector proteins, belonging to three families (cerato platanins, hydrofobins and glycoside hydrolases) were analyzed so far [20]. The first two groups were the most examined in case of T. harzianum. One of the cerato platanin proteins that increases plant root colonization efficiency of Trichoderma is swolenin [21,28].
For this experiment we have selected strain T.harzianum SMZC 22660, due to high activities of N-acetyl-β-glucosaminidase, acid phosphatase and naphtol-AS-BI-phosphohydrolase. Acid phosphatase plays important role in solubilization of bound phosphates, making them available to plants [29]. N-acetyl-β-glucosaminidase is used as main ingredient of plant protection products isolated from T. harzianum P1 [30]. All three enzymes have also been considered important for intensive biocontrol of soil-borne pathogens because of their ability to degrade fungal cell walls [31]. Indeed we have observed good in vitro antagonistic activities of the T.harzianum SZMS22660 strain. For example, values of Biocontrol Index (BCI) calculated from the confronted cultures against F. solani, R.solani, A. alternata, and P.cucurbitaceareum, ranged from 70% to 100%, respectively (unpublished data).
Positive effects of Trichoderma on different tomato genotypes in terms of increasing fresh and dry root mass, shoot dry mass, and stem height were observed in most cases [32]. Beneficial plant growth effects following T. harzianum inoculation have often been explained by the improved plant nutritional status [33]. However, with some genotypes the treatment had no effect or was even detrimental [32]. In our experiments, performed on two tomato cultivars grown in polytunnel, no statistically significant differences of examined growth parameters, neither between cultivars, nor between Trichoderma treatments and controls were observed, which is not in accordance with our previous results [23].
The non-destructive measurements of leaf chlorophyll have been used in many plant species and different environmental conditions providing accurate estimate of leaf chlorophyll content [34]. In this paper, we have used Dualex 4 Scientific leaf-clip that shows a linear relationship with chlorophyll content [35]. This apparatus simultaneously measures epidermal flavonols on the same leaf area, based on the fluorescence excitation ratio method [36], comparing chlorophyll fluorescence induced by ultraviolet (UV) radiation (375 nm) with that induced by red light (650 nm) [37]. A positive relationship (r2 = 0.97) between measurements of total flavonoid content from laboratory analysis and amounts detected non-destructively has been confirmed in medicinal plants [38], grapevine [39] and white cabbage (r2 = 0.93) [40]. The content of chlorophylls and epidermal flavonols in the leaves of plants is known to be an important indicator of nitrogen status in a plant [41]. The chlorophyll content is positively correlated with the nitrogen content while content of epidermal flavonols is inversely correlated to nitrogen content [42]. In our experiment, the lower chlorophyll content and higher epidermal flavonol content in plants grown in the presence of Trichoderma, was statistically significant in Gružanski zlatni, and it was in correlation with decreased nitrogen content. Others have also observed reduction in chlorophyll content of plants inoculated with Trichoderma asperellum TaspHu1 strain [43] and T. pseudoharzianum T1 or T. afroharzianum T52, that was connected with downregulation of chlorophyll synthesis genes in the youngest leaves [44]. While there are reports where no significant increase of chlorophyll a and chlorophyll b in tomato leaves of plants treated with various Trichoderma strains [45], the increase content of photosynthetic pigments and photosynthetic capability in plants grown in presence of Trichoderma species was observed by others [15,16]. Trichoderma presence significantly affected epidermal flavonol content in the opposite directions in both tomato cultivars in our experiment. Others have shown that flavonoids can be induced by symbionts in roots and suggested their involvement in the regulation of a temporary defense response in the root triggered by the symbiont invasion [46]. The increase of total flavonoids and total polyphenols was also observed in leaves from 7 days up to 60 days after Trichoderma foliar application [14]. Moreover, six flavonoids that were significantly increased in bean plants in the presence of Trichoderma strain with major stimulation of plant growth have been identified [47]. In addition, the detailed metabolome and transcriptome analysis have shown that T. harzianum colonization strongly affects and remodels phenylpropanoid pathway of tomato plants [48]. Our data on decreased NBI in Trichoderma treated tomato plants of cultivar Gružanski zlatni indicates a shift from primary to secondary metabolism, and in accordance with previously cited papers, better resistance to diseases or insects might be expected [47,48].
Tomato fruit contains phenolic compounds among which flavonoids accumulate mainly in the peel and highly contribute to its antioxidant activity [49]. Some authors have observed the increase of total phenolic content in tomato fruit [50], grapes [27], edible onion part [51] and cucumber [52], while we observed the decrease of total phenolic compounds in both cultivars in Trichoderma treatment. Similar observations were reported in the study which evaluated the effects of three Trichoderma bioactive metabolites (BAMs) applied to strawberry plants [53]. They recorded reduced levels of total antioxidant capacity (TAC), total phenolic content (TPC), total and individual anthocyanins, as well as for antioxidant proteins, which appear to generally indicate a ROS-enriched environment in the BAM-treated fruits. As we have observed, the increase in total flavonoid content in fruit of Trichoderma treated cultivar Gružanski zlatni, the decrease of total phenolic content in fruit is probably connected with decreased content of non-flavonoids (phenolic acids) in this cultivar. On the other hand, in the study conducted on plum tomatoes in the greenhouse, it was observed that polyphenols in fruits were not affected by T. harzianum strain T22 [13]. The authors discussed that it might be due to the basal constitutive amount of polyphenols in plum tomatoes that is independent of influences by external stimuli. Variations in polyphenol and flavonoid content of tomato fruit depend on tomato fruit variety and this variation could be due to genetic differences as well as different environmental stress conditions and agricultural practices that affect the chemical composition of plants [54,55]. The different response of the flavonoid content in onion varieties to T. asperellum and T22 inoculation was also observed [51]. Ripening of tomato fruits is connected to sugar accumulation. The breakdown of starch is followed by increased content of hexose sugars during formation of red fruit, but depending on applied treatments the content of soluble sugars can decrease with the reduction of starch content during ripening [8]. Determination of sugar content is important for commercial market as it is directly connected with tomato fruit sweetness [56]. To our knowledge, in the literature there are no data showing the effect of Trichoderma on the starch and soluble sugar content in tomato fruit. According to our results, genotype Gružanski zlatni contains significantly more soluble sugars in comparison to Narvik, and this was not affected by Trichoderma treatment. However, Trichoderma addition decreased starch content in genotype Gružanski zlatni.
Our results on the content of micro-elements and heavy metals, in both soil and fruit, as well as the BI values, resemble values reported previously [6]. Though BI for examined elements was similar in two cultivars in control conditions, they responded differently to Trichoderma treatment. Cultivar Gružanski zlatni had a positive response of BI for Fe and Cr. The literature data on the effect of Trichoderma inoculation to the plant uptake of nutrients are not consistent. For example, some authors reported that in response to inoculation with T. harzianum concentrations of Cu, Fe, Mn and Zn significantly increased in roots, shoots and fruits of tomato plants [10] and cucumber roots [33]. This is probably connected with increased availability of the nutrients to plants, as shown for P and Fe in T. asperellum-treated cucumber plants [33]. In contrast, decreased accumulation of iron in leaves of tomato plants grown in presence of T. harzianum T34 was observed [57]. Moreover, decreased concentrations of Cu, Mn and Zn were observed in wheat plants grown on a calcareous medium and inoculated with T. asperellum [11]. Both authors suggest that the decreased concentrations of these elements in plants were due to the competition between plants and Trichoderma. Likewise, study on tomato plants grown in hydroponics, with specific nutrient deficiency, indicated that the effect of T. harzianum inoculation is depended on the deficient element. In case of either Fe- or Cu-deficiency, inoculation was accompanied with the increased uptake of these elements. However, in case of Zn deficiency, the uptake of Zn was suppressed in inoculated plants, due to competition with Trichoderma [58]. The decreased BI in response to Trichoderma treatment was observed for heavy metals: Ni in Gružanski zlatni and Pb in both cultivars. We have also observed reduced accumulation and translocation of toxic levels of elements in the presence of the same strain of Trichoderma previously [23]. Accumulation of Cu was decreased in onion plants in the presence of T. asperellum in conditions of toxic Cu concentrations [59]. Our results for decreased BI for Pb are in accordance with results of authors who observed decrease in uptake of Pb in four examined plant species Miscanthus giganteus L., Salix spp., Phalaris arundinacea L. and Panicum virgatum L. grown in Trichoderma presence [60].
Principal component analysis (PCA) is often used to distinguish significance of the relationship between many examined plant parameters. This multivariate analysis method aims to explain the correlation between a large set of variables in terms of a small number of underlying independent factors [61]. PCA confirmed that the main plant traits that were different in two cultivars in control conditions were the higher content of epidermal flavonols in leaves of Narvik, and higher NBI in leaves of Gružanski zlatni, as a consequence of higher chlorophyll and but lower epidermal flavonols content (Figure 2; Table 2). Gružanski zlatni reacted to Trichoderma treatment by decreased chlorophyll content, increased flavonol content and consequently decreased NBI (Figure 2; Table 2). However, treatment of cultivar Narvik was connected only with lower flavonol content (Figure 2). The main fruit traits that were different in two cultivars in control conditions were higher content of Ni in fruits of Narvik, and higher contents of TP and N, but lower content of Pb in Gružanski zlatni (Figure 3; Table 4) The important traits connected with Trichoderma treatment were the increased Fe and Cr in fruits of Narvik, and increased Mn and TF in fruits of Gružanski zlatni.
Although T. harzianum is one of the most accepted biopesticides, it remains unclear how the positive effects of this fungi depend on the plant genotype [32]. The same authors suggest that the response to Trichoderma spp. is under genetic control and it can be explained that some tomato genotypes have the ability to regulate colonization by Trichoderma. Different tomato genotypes in different ways perform a perception of effectors molecules to identify and enable colonization of the fungus [20]. In our experiment, better root colonization, determined by higher expression of swolenin gene in tomato roots of more responsive tomato cultivar was in correlation with observed positive effects of Trichodrema.

4. Materials and Methods

4.1. Plant Material and Growth Conditions

The field experiment was conducted in Selenča (45°24′52.1′′ N, 19°18′52.4′′ E), in the Province of Vojvodina, Serbia, on an organic vegetable farm in conditions that have been reported previously [62].
Shortly, tomato was grown in a 20 m length, 4.5 m width and 2.5 m maximal height polytunnel, without additional lighting or heating. There were two plots of the same size (10 × 4.5 m) and the experiment was conducted in a randomized block design with two replicates. Irrigation was done by the drip system, daily for one hour equally for all plots. Copper sulfate was only added as a pesticide at two times intervals, prior to flowering and prior to fruit formation. Two medium late tomato cultivars were grown in experiments. They are often used in production for industrial and market use, due to good fruit quality and resistance to diseases. Gružanski zlatni (G2) is resistant to Fusarium oxysporum f. sp. Lycopersici, Verticillium albo-atrum and it is also tolerant to Phytophthora infestans. Narvik (G1) is resistant to Fusarium oxysporum f. sp. Lycopersici and Verticillium albo-atrum. Seeds were sown in plug trays with a peat/perlite (3:1, v/v) medium in February, and after twenty days, the uniform seedlings were planted in plastic pots. At the beginning of April, 20 plants per treatment (about 20cm high) were transplanted into four rows (spacing 50 cm within the row and 50 cm between rows), within the same day. At the end of the June, ten uniform plants from each treatment were chosen to collect fully light-exposed, red-ripe fruits.
Ten uniform plants per plot were used for the experiments. Conidial suspension of the strain Trichoderma harzianum SZMC 20660 (8 × 106 CFU/mL) was applied near the root zone, once during the growing period, in the first flowering phase.
Morphological and physiological parameters were measured 60 days after the treatment: soil water content (SWC), relative water content (RWC), epidermal flavonols (Flav) and total chlorophyll (Chl) (as described below), stem diameter base (STB), stem diameter middle (SDM), fruit weight (FW), fruit height (FH), fruit mass (FM) and number of fruit (FN) per plant. Four to seven red-ripe fruits [63] of similar size per plant were picked, measured (fruit width (FW) and height (FH)), washed, wiped, cut and stored at −80 °C for biochemical analysis.

4.2. Fungal Suspension Application

Strain T. harzianum SZMC 20660 (deposited in the Szeged Microbiological Collection (SZMC); Department of Microbiology) was isolated from the A horizon (5–30 cm) of agricultural soil used in organic agriculture. Species determination was based on their internal transcribed spacer (ITS) sequences [64].
Prior to preparation of fungal suspensions, T. harzianum isolate was pre-incubated at 25 °C in the dark. Suspensions were prepared as follows: pure culture of T. harzianum isolate was grabbed from a Petri dish, resuspended in 100 mL of tap water, and shaken for 2 h on 50 rpm. Conidial suspension (8 × 106 CFU/mL) was applied near the root zone, once during the growing period, in the first flowering phase.

4.3. Detection of the Enzymatic Activity

Trichoderma strain SZMC22660 were preincubated on PDA medium at 25 °C for 5 days. Mycelia discs of 5 mm diameter were cut from the edge of petri dish and inoculated into 200 mL of sterile potato dextrose broth (PDB) liquid medium. Flasks were shaken for 7 days on 25 °C in the dark on an orbital shaker BIOSAN (Lithuania). Liquid culture was filtrated through sterile gauze and filtrate was used for the detection of extracellular enzymatic activity of 19 enzymes (alkaline phosphatase, esterase [C4], esterase-lipase [C8], lipase [C14], leucine arylamidase, valine arylamidase, cistine arylamidase, trypsin, α-chymotripsin, acid phosphatase, Naphtol-AS-BI-phosphohydrolase; α-galactosidase, β-galactosidase, β-glucuronidase, α-glucosidase, β-glucosidase, N-acetil-β-glucosamidase, α-mannosidase, and α-fucosidase) were determined by using the semi quantitative API ZYM system (bioMerieux, Marcy-l’Étoilem, France), according to the manufacturer protocol.
Sixty-five micro liters of each inoculum was transferred into each of the twenty API ZYM strip microtubes prefilled with 5 mL of distilled water for maintaining humidity and incubated at 37 °C in thermostat for 4 h. After incubation period a drop of ZYM A and ZYM B reagents was added to each of twenty wells.
The color reaction was read after 5 min according to the API ZYM reading color scale. According to this scale enzyme activities were scored as following: 0—no color, 1—low, 2—moderate and 3—high.

4.4. RNA Extraction and Expression Analysis of Swolenin by Real-Time RT-PCR

Six plants per treatment were harvested 60 days after the treatment and frozen in liquid N2. Root samples were harvested from lateral roots after quick wash of soil particles. Total RNA was extracted from ~100 mg of frozen and grinded samples by RNeasy plant mini kit, Qiagen (Hilden, Germany), according to the manufacturer’s procedure. cDNA was isolated from ~1 µg of total RNA using oligo(dT) as a primer, by RevertAid First cDNA Synthesis Kit, Fermentas (Waltham, MA, USA).
PCR conditions for the expression analysis of swolenin were the following: 1xBuffer, 3 mM MgCl2, 0.4 mM dNTP, 0.3 µM Primer, 1U Taq polymerase, 2,5 ng BSA. PCR programe was the following: 1′ 940C; 30 cycles 15′′ 940C; 30′′ 600C; 1′45′′ 720C; final extension 10′ 720C.
Swolenin expression was followed by amplification of a 121-bp fragment (Accession number EU370698; Forward primer GTGGCCAGTGTGGAGGTATT; Reverse Primer GTGAGGGATCGAGGTAGCTG). The β-actin was used as reference gene. It was followed by amplification of a 280-bp fragment (Accession number BT013524; Forward primer GACCTGCTCCAC-CATCTTCC; Reverse Primer CAGTGGAGTTGCCGACAAAG) [21].

4.5. Soil Sample Analysis

For chemical investigations, samples were air dried and sieved through a 0.2-mm sieve prior to analysis. Soil pH in water, humus content, free CaCO3 content, available phosphorus (P2O5), available potassium (K2O) and total N were done as described previously [17]. According to the content of CaCO3 (1.41 ± 0.01%), humus (2.9 ± 0.28%) and total N (0.15 ± 0.01%), the soil was weakly calcareous, moderate humic soil with optimal content of N. Available phosphorus content was high (29 ± 0.49 mg P2O5/100 g DW) and available potassium content (23 ± 1.48 mg K2O/100 g DW) was optimal. The soil was weakly alkaline (pH in H2O 7.9 ± 0.01).

4.6. Determination of Soil water Content (SWC) and Leaf Relative Water (RWC) Content

SWC was measured using Theta probe (Delta-T, Cambridge, UK) at a depth of 6 cm. RWC in leaves was determined from the ratio RWC = [(FW − DW)/(TW − DW)] × 100. FW, DW and TW were fresh, dry and turgid weights of 10 leaf discs (2 r = 1 cm), respectively.

4.7. Measurements of Epidermal Flavonols, Total Chlorophyll, Nitrogen Balance Index and Antocyanine

Indices of chlorophyll (Chl), epidermal flavonols (Flav) and ther ratio, NBI as well as antocyanine (Ant) were measured in vivo non-destructively with Dualex sensor (Force-A, Orsay, France) [65]. Measurements were done on ten uniform, fully developed and sun-exposed leaves of the same plants that were later used for measurements of plant and fruit morpho-physiological parameters.

4.8. Extraction of Phenolics and Soluble Sugars

The 500 mg of dried tomato samples (70 °C) were extracted two times with 5 mL of 80% methanol using ultrasonic bath for 30 min at room temperature. After sonication, the samples were centrifuged at 1000× g for 10 min. at 4 °C. The supernatants were collected and were evaporated to dryness by rotary evaporator (Heidolph, Laborota 4000, Schwabach, Germany) under reduced pressure at 40 °C. The residues after evaporation were dissolved in 50 mL of milliQ water and these solutions were used for analysis of total phenolics (TP), total flavonoids (TF) and soluble sugars (SS).

4.9. Determination of Total Phenolic Content (TP)

TP in the samples was determined using Folin–Ciocalteu’s reagent [66]. Tomato extracts were diluted 1:9 with distilled water. Briefly, an aliquot of diluted sample (70 µL) was mixed with 300 µL of Folin–Ciocalteu reagent (diluted 1:10). The samples were mixed and after 5 min 230 µL of 7.5% Na2CO3 was added, followed by incubation for 30 min at 45 °C in dark. After standing for 1 h at room temperature in dark, the absorbance was measured at 765 nm using a Shimadzu spectrophotometer (UV-1800, Shimadzu USA Manufacturing Inc., Canby, UR, USA).TP was expressed as milligrams of gallic acid equivalents per 100 g of dried sample (mg GAE/100 g).

4.10. Determination of Total Flavonoid Content (TF)

TF was determined using a colorimetric assay with aluminium chloride [67]. Briefly, an aliquot of sample (125 µL) was mixed with 625 µL of milliQ water and 37.5 µL of 5% NaNO2. After 6 min, 75 µL of 10% AlCl3 was added to form a flavonoid–aluminium complex. After 5 min 250 µL of 1 M NaOH and 138 µL of distilled water were added in the mixture to make the total volume of 1.25 mL. The absorbance was measured immediately against the blank at 510 nm. TFC was expressed as milligrams of quercetin equivalents per 100 g of dried sample (mg QE/100 g).

4.11. Determination of Soluble Sugars Content (SS)

SS content was determined by the anthrone colorimetric assay using glucose as the standard [68]. Briefly, 1ml of diluted dried tomato extract was added to 5 mL of anthrone reagent. After 5 min, the mixture was allowed to stand for 10 min at boiling water bath. The absorbance was read at 620 nm after 10 min of the sample cooling. SS was expressed as milligrams per g of dried sample.

4.12. Determination of Starch Content (S)

S content was determined according to the method described previously [69]. Briefly, the pellet after soluble sugars extraction mixed with 3.25 mL of 52% perchloric acid and 2.5 mL of distilled water. After 30 min at 0 °C, the mixture centrifuged at 17,000 g for 15 min. The supernatant was collected and the procedure was repeated. The obtained supernatant after repeated starch extraction was collected and mixed with the previous one. The final volume of supernatants was adjusted to 10 mL with distilled water after neutralisation to pH 7 with 4M NaOH and saturated solution of Na2CO3. The starch content in the obtained extracts was determined by anthrone colorimetric assay already described for determination of soluble sugar contents. The obtained results were multiplied with 0.9 (the glucose conversion factor to starch), [70] and were expressed as milligram per gram of dried sample.

4.13. Determination of Macroelements

Total nitrogen content was determined by Kjeldahl method, total phosphorus content and potassium content were determined by spectrophotometry and flame photometry, respectively, after dry destruction (ashing) at 550 °C [71].

4.14. Determination of Micronutrients and Heavy Metal Content in Soil and Tomato Fruit

Samples were dried at 70 °C and grinded to a fine powder, weighed sample, 1 g, was placed in a flask and treated with 3 mL of concentrated HNO3 for 4–5 h. A mixture of HNO3 and HClO4 (ratio of 2:1) was added (3 mL) and 33% H2O2. The mixture was heated at 120–130 °C for 5–6 h, until fumes stopped and resulting solution was clear. Then, 10 mL of deionized water was added into the flask and the solution was boiled again for 10–15 min until the volume was reduced to the half, cooled to room temperature and filtered using Whatman filter paper No. 42. The entire filtrate was mixed and made the volume up to 50 mL with deionized water. All digestions and measurements were performed in triplicate. Metal content was measured by atomic absorption spectroscopy (AAS) (Varian Spec-trAA 220 FS).

4.15. Bioaccumulation Index (BI)

The plant efficiency for uptake of elements in the soil and transfer to the edible part was estimated by the BI. It was calculated as the percentage of the concentration in the edible part of the plant to that in the cultivation soil [7].
B I = M e t a l   c o n c e n t r a t i o n   i n   f r u i t M e t a l   c o n c e n t r a t i o n   i n   s o i l   ×   100

4.16. Statistical Analysis

Data are means of at least three replicates from four to ten different plants from each cultivar/treatment combinations. The results were analyzed by Tukey test using GraphPad Prism software. Principal component analysis (PCA) followed between groups analysis (BGA) was performed using the R platform (version 2.13.1, R Development Core Team, 2011) and the ade4TkGUI package was used [72]. The Pearson’s product moment correlation coefficient was calculated using the cor.test in order to disclose significant relationships between principal components and the variables analyzed (p < 0.05).

5. Conclusions

In our experiments, Trichoderma presence significantly affected Chl and Flav content and NBI in tomato leaves of cultivar Gružanski zlatni, which indicates a shift from primary to secondary metabolism. The increase in TF content in fruit of Trichoderma-treated plants of the same cultivar might contribute to its improved health promoting properties. The same cultivar responded to Trichoderma addition by decreased starch content, but increased BI for Fe and Cr and decreased BI for Ni. BI for Pb was decreased in both cultivars in response to Trichoderma treatment. Better root colonization, determined by higher expression of swolenin gene in tomato roots of more responsive tomato cultivar, was in correlation with observed positive effects of Trichoderma.

Author Contributions

Writing—original draft preparation, I.D.V., G.M.R. and M.M.B.; investigation, I.D.V., G.M.R., M.M.B., M.B.P. and L.T.P.; formal analysis, I.D.V. and G.M.R., methodology, D.M.P. and L.T.P.; writing—review and editing D.M.P., H.M.K. and E.M.S.; supervision, D.M.P. All authors have read and agreed to the published version of the manuscript.

Funding

The research was supported by the Ministry of Education, Science and Technological development of the Republic of Serbia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets used and/or analyzed during the current study are available from the corresponding authors on reasonable request.

Acknowledgments

The research was co-financed by the Ministry of Education, Science and Technological development of the Republic of Serbia (Project No. III43010, TR 31080 and TR 31005). The authors are grateful to Carla Pinheiro from FCT-NOVA, New University of Lisbon, Portugal, for helpful discussions on the statistical analysis of data presented in this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ampt, E.A.; van Ruijven, J.; Raaijmakers, J.M.; Termorshuizen, A.J.; Mommer, L. Linking ecology and plant pathology to unravel the importance of soil-borne fungal pathogens in species-rich grasslands. Eur. J. Plant Pathol. 2019, 154, 141–156. [Google Scholar] [CrossRef] [Green Version]
  2. Frusciante, L.; Carli, P.; Ercolano, M.R.; Pernice, R.; Di Matteo, A.; Fogliano, V.; Pellegrini, N. Antioxidant nutritional quality of tomato. Mol. Nutr. Food Res. 2007, 51, 609–617. [Google Scholar] [CrossRef] [PubMed]
  3. Kanayama, Y. Sugar metabolism and fruit development in the tomato. Hort. J. 2017, 86, 417–425. [Google Scholar] [CrossRef] [Green Version]
  4. Bhowmik, D.; Kumar, K.S.; Paswan, S.; Srivastava, S. Tomato-a natural medicine and its health benefits. J. Pharmacogn. Phytochem. 2012, 1, 33–43. [Google Scholar]
  5. Wang, X.; Xing, Y. Evaluation of the effects of irrigation and fertilization on tomato fruit yield and quality: A principal component analysis. Sci. Rep. 2017, 7, 350. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Klunklin, W.; Savage, G. Effect on quality characteristics of tomatoes grown under well-watered and drought stress conditions. Foods 2017, 6, 56. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Rodriguez-Iruretagoiena, A.; Trebolazabala, J.; Martinez-Arkarazo, I.; de Diego, A.; Madariaga, J.M. Metals and metalloids in fruits of tomatoes (Solanum lycopersicum) and their cultivation soils in the Basque Country: Concentrations and accumulation trends. Food Chem. 2015, 173, 1083–1089. [Google Scholar] [CrossRef]
  8. Centeno, D.C.; Osorio, S.; Nunes-Nesi, A.; Bertolo, A.L.F.; Carneiro, R.T.; Araújo, W.L.; Steinhauser, M.-C.; Michalska, J.; Rohrmann, J.; Geigenberger, P.; et al. Malate Plays a Crucial Role in Starch Metabolism, Ripening, and Soluble Solid Content of Tomato Fruit and Affects Postharvest Softening. Plant Cell 2011, 23, 162. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Wang, W.; Shi, J.; Xie, Q.; Jiang, Y.; Yu, N.; Wang, E. Nutrient exchange and regulation in arbuscular mycorrhizal symbiosis. Mol. Plant 2017, 10, 1147–1158. [Google Scholar] [CrossRef] [Green Version]
  10. Khan, M.Y.; Haque, M.M.; Molla, A.H.; Rahman, M.M.; Alam, M.Z. Antioxidant compounds and minerals in tomatoes by Trichoderma-enriched biofertilizer and their relationship with the soil environments. J. Integr. Agric. 2017, 16, 691–703. [Google Scholar] [CrossRef]
  11. de Santiago, A.; Quintero, J.M.; Avilés, M.; Delgado, A. Effect of Trichoderma asperellum strain T34 on iron, copper, manganese, and zinc uptake by wheat grown on a calcareous medium. Plant Soil 2011, 342, 97–104. [Google Scholar] [CrossRef]
  12. Li, Y.T.; Hwang, S.G.; Huang, Y.M.; Huang, C.H. Effects of Trichoderma asperellum on nutrient uptake and Fusarium wilt of tomato. Crop Prot. 2018, 110, 275–282. [Google Scholar] [CrossRef]
  13. Carillo, P.; Woo, S.L.; Comite, E.; El-Nakhel, C.; Rouphael, Y.; Fusco, G.M.; Borzacchiello, A.; Lanzuise, S.; Vinale, F. Application of Trichoderma harzianum, 6-pentyl-α-pyrone and plant biopolymer formulations modulate plant metabolism and fruit quality of plum tomatoes. Plants 2020, 9, 771. [Google Scholar] [CrossRef] [PubMed]
  14. Şesan, T.E.; Oancea, A.O.; Ştefan, L.M.; Mănoiu, V.S.; Ghiurea, M.; Răut, I.; Pomohaci, C.M. Effects of foliar treatment with a Trichoderma plant biostimulant consortium on Passiflora caerulea L. yield and quality. Microorganisms 2020, 8, 123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Fiorini, L.; Guglielminetti, L.; Mariotti, L.; Curadi, M.; Picciarelli, P.; Scartazza, A.; Sarrocco, S.; Vannacci, G. Trichoderma harzianum T6776 modulates a complex metabolic network to stimulate tomato cv. Micro-Tom growth. Plant Soil 2016, 400, 351–366. [Google Scholar] [CrossRef]
  16. Harman, G.E.; Doni, F.; Khadka, R.B.; Uphoff, N. Endophytic strains of Trichoderma increase plants’ photosynthetic capability. J. Appl. 2021, 130, 529–546. [Google Scholar] [CrossRef] [PubMed]
  17. Ruiz-Cisneros, M.F.; Ornelas-Paz, J.J.; Olivas-Orozco, G.I.; Acosta-Muñiz, C.H.; Sepúlveda-Ahumada, D.R.; Pérez-Corral, D.A.; Rios-Velasco, C.; Salas-Marina, M.A.; Fernández-Pavía, S.P. Efect of Trichoderma spp. and phytopathogenic fungi on plant growth and tomato fruit quality. Rev. Mex. Fitopatol. 2018, 36, 444–456. [Google Scholar]
  18. Viterbo, A.D.A.; Chet, I. TasHyd1, a new hydrophobin gene from the biocontrol agent Trichoderma asperellum, is involved in plant root colonization. Mol. Plant Pathol. 2006, 7, 249–258. [Google Scholar] [CrossRef]
  19. Samolski, I.; Rinco’ n, A.M.; Pinzo’ n, L.M.; Viterbo, A.; Monte, E. The qid74 gene from Trichoderma harzianum has a role in root architecture and plant biofertilization. Microbiology 2012, 158, 129–138. [Google Scholar] [CrossRef]
  20. Ramírez-Valdespino, C.A.; Casas-Flores, S.; Olmedo-Monfil, V. Trichoderma as a model to study effector-like molecules. Front. Biol. 2019, 10, 1030. [Google Scholar] [CrossRef] [PubMed]
  21. Brotman, Y.; Briff, E.; Viterbo, A.; Chet, I. Role of swollenin, an expansin-like protein from Trichoderma, in plant root colonization. Plant Physiol. 2008, 147, 779–789. [Google Scholar] [CrossRef] [Green Version]
  22. Brotman, Y.; Lisec, J.; Méret, M.; Chet, I.; Willmitzer, L.; Viterbo, A. Transcript and metabolite analysis of the Trichoderma-induced systemic resistance response to Pseudomonas syringae in Arabidopsis thaliana. Microbiology 2012, 158, 139–146. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Vukelić, I.; Racić, G.; Bojović, M.; Ćurčić, N.; Mrkajić, D.; Jovanović, D.; Panković, D. Effect of Trichoderma harzianum on morpho-physiological parametars and metal uptake of tomato plants. Matica Srpska J. Nat. Sci. 2020, 139, 61–71. [Google Scholar] [CrossRef]
  24. Vinale, F.; Sivasithamparam, K.; Ghisalberti, E.L.; Woo, S.L.; Nigro, M.; Marra, R.; Manganiello, G. Trichoderma secondary metabolites active on plants and fungal pathogens. Open Mycol. J. 2014, 8, 127–139. [Google Scholar] [CrossRef] [Green Version]
  25. Mehetre, S.T.; Mukherjee, P.K. Trichoderma improves nutrient use efficiency in crop plants. In Nutrient Use Efficiency: From Basics to Advances, 1st ed.; Rakshit, A., Bahadur, H.B., Singh, A.S., Eds.; Springer: New Delhi, India, 2015; pp. 173–180. [Google Scholar]
  26. Woo, S.L.; Ruocco, M.; Vinale, F.; Nigro, M.; Marra, R.; Lombardi, N.; Lorito, M. Trichoderma-based products and their widespread use in agriculture. Open Mycol. J. 2014, 8, 71–126. [Google Scholar] [CrossRef] [Green Version]
  27. Pascale, A.; Vinale, F.; Manganiello, G.; Nigro, M.; Lanzuise, S.; Ruocco, M.; Lorito, M. Trichoderma and its secondary metabolites improve yield and quality of grapes. Crop Prot. 2017, 92, 176–181. [Google Scholar] [CrossRef] [Green Version]
  28. Brotman, Y.; Landau, U.; Cuadros-Inostroza, Á.; Takayuki, T.; Fernie, A.R.; Chet, I.; Viterbo, A.; Willmitzer, L. Trichoderma-plant root colonization: Escaping early plant defense responses and activation of the antioxidant machinery for saline stress tolerance. PLoS Pathog. 2013, 9, e1003221. [Google Scholar] [CrossRef]
  29. Kapri, A.; Tewari, L. Phosphate solubilization potential and phosphatase activity of rhizospheric Trichoderma spp. Braz. J. Microbiol. 2010, 41, 787–795. [Google Scholar] [CrossRef] [Green Version]
  30. Harman, G.E.; Lorito, M.; Di Pietro, A.; Hayes, C.K. N-acetyl-β-glucosaminidase isolated from Trichoderma harzianum. U.S. Patent 5,474,926, 12 December 1995. [Google Scholar]
  31. Harman, G.E.; Howell, C.R.; Viterbo, A.; Chet, I.; Lorito, M. Trichoderma species—Opportunistic, avirulent plant symbionts. Nat. Rev. Microbiol. 2004, 2, 43–56. [Google Scholar] [CrossRef]
  32. Tucci, M.; Ruocco, M.; De Masi, L.; De Palma, M.; Lorito, M. The beneficial effect of Trichoderma spp. on tomato is modulated by the plant genotype. Mol. Plant Pathol. 2012, 12, 341–354. [Google Scholar] [CrossRef]
  33. Yedidia, I.; Srivastva, A.K.; Kapulnik, Y.; Chet, I. Effect of Trichoderma harzianum on microelement concentrations and increased growth of cucumber plants. Plant Soil. 2001, 235, 235–242. [Google Scholar] [CrossRef]
  34. Kalaji, H.M.; Schansker, G.; Brestic, M.; Bussotti, F.; Calatayud, A.; Ferroni, L.; Goltsev, V.; Guidi, L.; Jajoo, A.; Li, P.; et al. Frequently asked questions about chlorophyll fluorescence, the sequel. Photosynth. Res. 2017, 132, 13–66. [Google Scholar] [CrossRef] [Green Version]
  35. Cerovic, Z.G.; Masdoumier, G.; Ghozlen, N.B.; Latouche, G. A new optical leaf-clip meter for simultaneous non-destructive assessment of leaf chlorophyll and epidermal flavonoids. Physiol. Plant. 2012, 146, 251–260. [Google Scholar] [CrossRef] [PubMed]
  36. Cerovic, Z.G.; Ounis, A.; Cartelat, A.; Latouche, G.; Goulas, Y.; Meyer, S.; Moya, I. The use of chlorophyll fluorescence excitation spectra for the non-destructive in situ assessment of UV-absorbing compounds in leaves. Plant Cell Environ. 2002, 25, 1663–1676. [Google Scholar] [CrossRef]
  37. Goulas, Y.; Cerovic, Z.G.; Cartelat, A.; Moya, I. Dualex: A new instrument for field measurements of epidermal ultraviolet absorbance by chlorophyll fluorescence. Appl. Opt. 2004, 43, 4488–4496. [Google Scholar] [CrossRef]
  38. Ibrahim, M.H.; Jaafar, H.Z.E. Dualex 4 reading values as indicator of total flavonoid and chlorophyll content of three varieties of Labisia pumila Benth.(Kacip Fatimah) seedlings under greenhouse condition. Trans. Malays. Soc. Plant Physiol. 2012, 20, 102–105. [Google Scholar]
  39. Agati, G.; Cerovic, Z.G.; Dalla Marta, A.; Di Stefano, V.; Pinelli, P.; Traversi, M.L.; Orlandini, S. Optically-assessed preformed flavonoids and susceptibility of grapevine to Plasmopara viticola under different light regimes. Funct. Plant Biol. 2008, 35, 77–84. [Google Scholar] [CrossRef] [PubMed]
  40. Agati, G.; Tuccio, L.; Kusznierewicz, B.; Chmiel, T.; Bartoszek, A.; Kowalski, A.; Kaniszewski, S.; Grzegorzewska, M.; Kosson, R.; Kaniszewski, S. Nondestructive optical sensing of flavonols and chlorophyll in white head cabbage (Brassica oleracea L. var. capitata subvar. alba) grown under different nitrogen regimens. J. Agric. Food Chem. 2016, 64, 85–94. [Google Scholar] [CrossRef]
  41. Cartelat, A.; Cerovic, Z.G.; Goulas, Y.; Meyer, S.; Lelarge, C.; Prioul, J.L.; Barbottin, A.; Jeuffroy, M.H.; Gate, P.; Agati, G.; et al. Optically assessed contents of leaf polyphenolics and chlorophyll as indicators of nitrogen deficiency in wheat (Triticum aestivum L.). Field Crops Res. 2005, 91, 35–49. [Google Scholar] [CrossRef]
  42. Guller, N.S.; Pehlivan, N.; Karaoglu, S.A.; Guzel, S.; Bozdeveci, A. Trichoderma atroviride ID20G inoculation ameliorates drought stress-induced damages by improving antioxidant defence in maize seedlings. Acta Physiol. Plant. 2016, 38, 19. [Google Scholar] [CrossRef]
  43. Yu, Z.; Wang, Z.; Zhang, Y.; Wang, Y.; Liu, Z. Biocontrol and growth-promoting effect of Trichoderma asperellum TaspHu1 isolate from Juglans mandshurica rhizosphere soil. Microbiol. Res. 2021, 242, 126596. [Google Scholar] [CrossRef]
  44. Liu, B.; Ji, S.; Zhang, H.; Wang, Y.; Liu, Z. Isolation of Trichoderma in the rhizosphere soil of Syringa oblata from Harbin and their biocontrol and growth promotion function. Microbiol. Res. 2020, 235, 126445. [Google Scholar] [CrossRef] [PubMed]
  45. Alexandru, M.; Lazăr, D.; Ene, M.; Sesan, T.E. Influence of some Trichoderma species on photosynthesis intensity and pigments in tomatoes. Rom. Biotech. Lett. 2013, 18, 8499–8510. [Google Scholar]
  46. Hassan, S.; Mathesius, U. The role of flavonoids in root-rhizosphere signalling: Opportunities and challenges for improving plant-microbe interactions. J. Exp. Bot. 2012, 63, 3429–3444. [Google Scholar] [CrossRef] [Green Version]
  47. Mayo-Prieto, S.; Marra, R.; Vinale, F.; Rodríguez-González, Á.; Woo, S.L.; Lorito, M.; Casquero, P.A. Effect of Trichoderma velutinum and Rhizoctonia solani on the metabolome of bean plants (Phaseolus vulgaris L.). Int. J. Mol. Sci. 2019, 20, 549. [Google Scholar] [CrossRef] [Green Version]
  48. Coppola, M.; Diretto, G.; Digilio, M.C.; Woo, S.L.; Giuliano, G.; Molisso, D.; Rao, R. Transcriptome and metabolome reprogramming in tomato plants by Trichoderma harzianum strain T22 primes and enhances defence responses against aphids. Front. physiol. 2019, 10, 745. [Google Scholar] [CrossRef]
  49. Bhandari, S.R.; Cho, M.C.; Lee, J.G. Genotypic variation in carotenoid, ascorbic acid, total phenolic, and flavonoid contents, and antioxidant activity in selected tomato breeding lines. Hortic. Environ. Biotechnol. 2016, 57, 440–452. [Google Scholar] [CrossRef]
  50. Singh, S.P.; Singh, H.B.; Singh, D.K. Effect of Trichoderma harzianum on mineral component and antioxidant activity of tomato fruits. Vegetos 2013, 26, 237–244. [Google Scholar] [CrossRef]
  51. Ortega-García, J.G.; Montes-Belmont, R.; Rodríguez-Monroy, M.; Ramírez-Trujillo, J.A.; Suárez-Rodríguez, R.; Sepúlveda-Jiménez, G. Effect of Trichoderma asperellum applications and mineral fertilization on growth promotion and the content of phenolic compounds and flavonoids in onions. Sci. Hortic. 2015, 195, 8–16. [Google Scholar] [CrossRef]
  52. NawRocka, J.; Szczech, M.; MałolepSza, U. Trichoderma atroviride enhances phenolic synthesis and cucumber protection against Rhizoctonia solani. Plant Prot. Sci. 2017, 54, 17–23. [Google Scholar]
  53. Lombardi, N.; Salzano, A.M.; Troise, A.D.; Scaloni, A.; Vitaglione, P.; Vinale, F.; Marra, R.; Caira, S.; Lorito, M.; d’Errico, G.; et al. Effect of Trichoderma Bioactive Metabolite Treatments on the Production, Quality, and Protein Profile of Strawberry Fruits. J. Agric. Food Chem. 2020, 68, 7246–7258. [Google Scholar] [CrossRef]
  54. Martinez-Valverde, I.; Periago, M.J.; Provan, G.; Chesson, A. Phenolic compounds, lycopene and antioxidant activity in commercial varieties of tomato (Lycopersicum esculentum). J. Sci. Food. Agric. 2002, 82, 323–330. [Google Scholar] [CrossRef]
  55. Silva-Beltrán, N.P.; Ruiz-Cruz, S.; Cira-Chávez, L.A.; Estrada-Alvarado, M.I.; Ornelas-Paz, J.D.J.; López-Mata, M.A.; Del-Toro-Sánchez, C.L.; Ayala-Zavala, J.F.; Márquez-Ríos, E. Total phenolic, flavonoid, tomatine, and tomatidine contents and antioxidant and antimicrobial activities of extracts of tomato plant. Int. J. Anal. Chem. 2015, 2015, 284071. [Google Scholar] [CrossRef] [Green Version]
  56. Matsukura, C. Sugar accumulation in tomato fruit and its modification using molecular breeding techniques. In Functional Genomics and Biotechnology in Solanaceae and Cucurbitaceae Crops, 1st ed.; Ezura, H., Ariizumi, T., Garcia-Mas, J., Rose, J., Eds.; Springer: Berlin/Heidelberg, Germany, 2016; pp. 141–154. [Google Scholar]
  57. Borrero, C.; Trillas, M.I.; Delgado, A.; Avilés, M. Effect of ammonium/nitrate ratio in nutrient solution on control of Fusarium wilt of tomato by Trichoderma asperellum T34. Plant Pathol. 2012, 61, 132–139. [Google Scholar] [CrossRef]
  58. Li, R.X.; Cai, F.; Pang, G.; Shen, Q.R.; Li, R.; Chen, W. Solubilisation of phosphate and 415 micronutrients by Trichoderma harzianum and its relationship with the promotion of 416 tomato plant growth. PLoS ONE 2015, 10, e0130081. [Google Scholar]
  59. Vargas, J.T.; Rodríguez-Monroy, M.; Meyer, M.L.; Montes-Belmont, R.; Sepúlveda-Jiménez, G. Trichoderma asperellum ameliorates phytotoxic effects of copper in onion (Allium cepa L.). Environ. Exp. Bot. 2017, 136, 85–93. [Google Scholar] [CrossRef]
  60. Kacprzak, M.J.; Rosikon, K.; Fijalkowski, K.; Grobelak, A. The effect of Trichoderma on heavy metal mobility and uptake by Miscanthus giganteus, Salix sp.; Phalarisarundinacea, and Panicum virgatum. Appl. Environ. Soil Sci. 2014, 2014, 506142. [Google Scholar] [CrossRef] [Green Version]
  61. Beheshtizadeh, H.; Rezaie, A.; Rezaie, A.; Ghandi, A. Principal component analysis and determination of the selection criteria in bread wheat (Triticum aestivum L.) genotypes. Int. J. Agric. Crop Sci. 2013, 5, 2024. [Google Scholar]
  62. Zivanovic, B.; Vidovic, M.; Komic, S.M.; Jovanovic, L.; Kolarz, P.; Morina, F.; Jovanovic, S.V. Contents of phenolics and carotenoids in tomato grown under polytunnels with different UV-transmission rates. Turk. J. Agric. For. 2017, 41, 113–120. [Google Scholar] [CrossRef]
  63. Grierson, D.; Kader, A.A. Fruit Ripening and Quality. In The Tomato Crop, 1st ed.; Atherton, J., Rudich, J., Eds.; Springer: Dordrecht, The Netherlands, 1986; pp. 241–280. [Google Scholar]
  64. Racić, G.; Körmöczi, P.; Kredics, L.; Raičević, V.; Mutavdžić, B.; Vrvić, M.M.; Panković, D. Effect of the edaphic factors and metal content in soil on the diversity of Trichoderma spp. Environ. Sci. Pollut. Res. 2017, 24, 3375–3386. [Google Scholar] [CrossRef] [Green Version]
  65. Racić, G.; Vukelić, I.; Prokić, L.; Ćurčić, N.; Zorić, M.; Jovanović, L.; Panković, D. The influence of Trichoderma brevicompactum treatment and drought on physiological parameters, abscisic acid content and signalling pathway marker gene expression in leaves and roots of tomato. Ann. Appl. Biol. 2018, 173, 213–221. [Google Scholar] [CrossRef]
  66. Singleton, V.L.; Orthofer, R.; Lamuela-Raventós, R.M. Analysis of total phenols and other oxidation substrates and antioxidants by means of Folin-Ciocalteu reagent. Methods Enzymol. 1998, 299, 52–178. [Google Scholar]
  67. Dewanto, V.; Wu, X.; Adom, K.K.; Liu, R.H. Thermal Processing Enhances the Nutritional Value of Tomatoes by Increasing Total Antioxidant Activity. J. Agric. Food Chem. 2002, 50, 3010–3014. [Google Scholar] [CrossRef]
  68. Yemm, E.W.; Willis, A. The estimation of carbohydrates in plant extracts by anthrone. Biochem. J. 1954, 57, 508. [Google Scholar] [CrossRef] [Green Version]
  69. Kostić, A.Ž.; Gašić, U.M.; Pešić, M.B.; Stanojević, S.P.; Barać, M.B.; Mačukanović-Jocić, M.P.; Avramov, S.N.; Tešić, Ž.L. Phytochemical analysis and total antioxidant capacity of rhizome, above-ground vegetative parts and flower of three iris species. Chem. Biodivers. 2019, 16, e1800565. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Rose, R.; Rose, C.L.; Omi, S.K.; Forry, K.R.; Durall, D.M.; Bigg, W.L. Starch determination by perchloric acid vs enzymes: Evaluating the accuracy and precision of six colorimetric methods. J. Agric. Food Chem. 1991, 39, 2–11. [Google Scholar] [CrossRef]
  71. Jones, J.B. Laboratory Guide for Conducting Soil Tests and Plant Analysis, 1st ed.; CRC Press: Boca Raton, FL, USA, 2001; pp. 109–121. [Google Scholar]
  72. Thioulouse, J.; Dray, S. ade4TkGUI: ade4 Tcl/Tk Graphical User Interface. R Package Version 0.2-5. 2009. Available online: http://CRAN.R-project.org/package=ade4TkGUI (accessed on 20 June 2021).
Ijms 22 06961 g001 550
Figure 1. Expression of the swolenin gene (A) and β-actin reference gene (B) in the roots of tomato variety Narvik (1, 2) and Gružanski zlatni (3, 4). B refers to blank and M referes to the 100 bp DNA Ladder H3 RTU 50 μg/500 μL (A) and 50 bp DNA Ladder MWD50 50 μg/500 μL (B).
Figure 1. Expression of the swolenin gene (A) and β-actin reference gene (B) in the roots of tomato variety Narvik (1, 2) and Gružanski zlatni (3, 4). B refers to blank and M referes to the 100 bp DNA Ladder H3 RTU 50 μg/500 μL (A) and 50 bp DNA Ladder MWD50 50 μg/500 μL (B).
Ijms 22 06961 g001
Ijms 22 06961 g002 550
Figure 2. (a) PCA analysis of Narvik (G1) and Gružanski zlatni (G2) tomato cultivars, grown in the polytunnel without (C) or with T. harzianum application (T), involving chlorophyll (Chl), epidermal flavonols (Flav), nitrogen balance index (NBI), anthocyanin (Anth), relative water content (RWC (%)), stem diameter base (STB (cm)) and stem diameter middle (STM (cm)). (b) Correlation matrix (Pearson’s product correlation coefficient) with the associated p value between the principal components and the analyzed variables.
Figure 2. (a) PCA analysis of Narvik (G1) and Gružanski zlatni (G2) tomato cultivars, grown in the polytunnel without (C) or with T. harzianum application (T), involving chlorophyll (Chl), epidermal flavonols (Flav), nitrogen balance index (NBI), anthocyanin (Anth), relative water content (RWC (%)), stem diameter base (STB (cm)) and stem diameter middle (STM (cm)). (b) Correlation matrix (Pearson’s product correlation coefficient) with the associated p value between the principal components and the analyzed variables.
Ijms 22 06961 g002
Ijms 22 06961 g003 550
Figure 3. (a) PCA analysis of Narvik (G1) and Gružanski zlatni (G2) tomato cultivars, grown in the polytunnel without (C) or with T. harzianum application (T), involving: fruit weight (FW (cm 10−2)), fruit height (FH (cm 10−2)), number of fruit per plant (FN), fruit mass (FM (g)), total polyphenols (TP (GEA mg/100 gDW)), total flavonoids (TF (mg/100 gDW)), starch (S (mg/g DW)), contents of soluble sugars (SS (mg/g DW)), macroelements N (%), P (ppm), K (ppm) and heavy metals (ppm): Mn, Fe, Cu, Zn, Ni, Cr, Pb, and Co. (b) Correlation matrix (Pearson´s product correlation coefficient) with the associated p value between the principal components and the analyzed variables.
Figure 3. (a) PCA analysis of Narvik (G1) and Gružanski zlatni (G2) tomato cultivars, grown in the polytunnel without (C) or with T. harzianum application (T), involving: fruit weight (FW (cm 10−2)), fruit height (FH (cm 10−2)), number of fruit per plant (FN), fruit mass (FM (g)), total polyphenols (TP (GEA mg/100 gDW)), total flavonoids (TF (mg/100 gDW)), starch (S (mg/g DW)), contents of soluble sugars (SS (mg/g DW)), macroelements N (%), P (ppm), K (ppm) and heavy metals (ppm): Mn, Fe, Cu, Zn, Ni, Cr, Pb, and Co. (b) Correlation matrix (Pearson´s product correlation coefficient) with the associated p value between the principal components and the analyzed variables.
Ijms 22 06961 g003
Table
Table 1. Extracellular enzyme activities as determined with API-ZYM system (alkaline phosphatase (A), esterase C4 (B), esterase lipase C8 (C), lipase C14 (D), leucine arylamidase (E), valine arylamidase (F), cysteine arylamidase (G), trypsin (H), α-chymotripsin (I), acid phosphatase (J), naphthol-AS-BI-phosphohydrolase (K), α-galactosidase (L), β-galactosidase (M), β-glucuronidase (N), α-glucosidase (O), β- glucosidase (P), N-acetyl-β-glucosaminidase (Q), α-mannosidase (R), α-mannosidase (S)). According to this scale, enzyme activities were scored as following: 0—no color, 1—low, 2—moderate and 3—high.
Table 1. Extracellular enzyme activities as determined with API-ZYM system (alkaline phosphatase (A), esterase C4 (B), esterase lipase C8 (C), lipase C14 (D), leucine arylamidase (E), valine arylamidase (F), cysteine arylamidase (G), trypsin (H), α-chymotripsin (I), acid phosphatase (J), naphthol-AS-BI-phosphohydrolase (K), α-galactosidase (L), β-galactosidase (M), β-glucuronidase (N), α-glucosidase (O), β- glucosidase (P), N-acetyl-β-glucosaminidase (Q), α-mannosidase (R), α-mannosidase (S)). According to this scale, enzyme activities were scored as following: 0—no color, 1—low, 2—moderate and 3—high.
ABCDEFGHIJKLMNOPQRS
T. harzianum SZMC 226601110000003300000300
Table
Table 2. Mean values of soil water content (SWC), relative leaf water content (RWC), indices of chlorophyll (Chl), epidermal flavonols (Flav), nitrogen balance index (NBI) and anthocyanin (Anth), as determined on ten plants of two tomato cultivars per treatment: cultivar Narvik (G1) and cultivar Gružanski zlatni (G2), without (C) or with T. harzianum application (T).
Table 2. Mean values of soil water content (SWC), relative leaf water content (RWC), indices of chlorophyll (Chl), epidermal flavonols (Flav), nitrogen balance index (NBI) and anthocyanin (Anth), as determined on ten plants of two tomato cultivars per treatment: cultivar Narvik (G1) and cultivar Gružanski zlatni (G2), without (C) or with T. harzianum application (T).
SWC (%)RWC (%)ChlFlavNBIAnth
G1C24.50 ± 2.4 ab49.32 ± 3.7 a29.1 ± 1.8 c0.67 ± 0.03 a43.15 ± 1.2 c0.11 ± 0.005 a
G1T27.33 ± 3.0 a48.73 ± 4.7 a26.4 ± 1.2 c0.59 ± 0.02 b48.23 ± 3.8 c0.12 ± 0.008 a
G2C30.75 ± 0.8 a50.18 ± 4.3 a41.7 ± 1.1 a0.52 ± 0.03 c77.5 ± 1.5 a0.08 ± 0.008 b
G2T30.95 ± 1.1 a54.67 ± 2.6 a36.53 ± 1.2 b0.65 ± 0.03 ab61.25 ± 4.6 b0.09 ± 0.008 b
Means with the same letter in the same column are not significantly different from each other according to Tukey’s test (p < 0.05). Error lines represent ± standard deviation of the mean. Different letters (a, b, c) indicate statistically significant differences according to Tukey’s test (p < 0.05).
Table
Table 3. Mean values of plant (stem diameter base (STB) and stem diameter middle (SDM)) and fruit (fruit weight (FW), fruit height (FH), number of fruit (FN) and fruit mass (FM)) morphological characteristics, measured on ten plants of two tomato cultivars per treatment: cultivar Narvik (G1) and cultivar Gružanski zlatni (G2), without (C) or with T. harzianum application (T).
Table 3. Mean values of plant (stem diameter base (STB) and stem diameter middle (SDM)) and fruit (fruit weight (FW), fruit height (FH), number of fruit (FN) and fruit mass (FM)) morphological characteristics, measured on ten plants of two tomato cultivars per treatment: cultivar Narvik (G1) and cultivar Gružanski zlatni (G2), without (C) or with T. harzianum application (T).
STB
(cm)		SDM
(cm)		FW
(cm × 10−2)		FH
(cm × 10−2)		FNFM
(g)
G1C16.36 ± 0.76 a10.6 ± 1.11 a60.03 ± 2.44 a57.8 ± 3.87 a12 ± 2 b130 ± 9.29 a
G1T16.7 ± 1.31 a12.36 ± 1.01 a61.96 ± 2.05 a61.3 ± 2.87 a11 ± 1 b 143 ± 3.78 a
G2C13.87 ± 1.45 ab11.6 ± 1.05 a59.5 ± 4.37 a57.1 ± 1.42 a21 ± 3 a 133 ± 4.57 a
G2T11.4 ± 0.87 b10.3 ± 0.81 a61.3 ± 3.61 a60.6 ± 1.57 a20 ± 1 a 141 ± 11.35 a
Means with the same letter in the same column are not significantly different from each other according to Tukey’s test (p < 0.05). Error lines represent ± standard deviation of the mean. Different letters (a, b) indicate statistically significant differences according to Tukey’s test (p < 0.05).
Table
Table 4. Mean values of total polyphenols (TP), total flavonoids (TF), starch (S), contents of soluble sugars (SS), total nitrogen (N), phosphorus (P) and potassium (K) were measured on bulked samples of ten plants of two tomato cultivars per treatment: cultivar Narvik (G1) and cultivar Gružanski zlatni (G2), without (C) or with T. harzianum application (T).
Table 4. Mean values of total polyphenols (TP), total flavonoids (TF), starch (S), contents of soluble sugars (SS), total nitrogen (N), phosphorus (P) and potassium (K) were measured on bulked samples of ten plants of two tomato cultivars per treatment: cultivar Narvik (G1) and cultivar Gružanski zlatni (G2), without (C) or with T. harzianum application (T).
TP
(GEA mg/100 gDW)		TF
(mg/100 gDW)		S
(mg/g DW)		SS
(mg/g DW)		N (%)P (ppm)K (ppm)
G1C1632.41 ± 68.74 b455.56 ± 47.13 b30.61 ± 3.81 ab225.18 ± 14.96 bc1.28 ± 0.03 c0.79 ± 0.02 a2.85 ± 0.28 a
G1T1408.29 ± 45.21 c333.26 ± 9.51 c37.65 ± 2.01 a192 ± 31.57 c1.27 ± 0.02 c0.77 ± 0.015 a2.98 ± 0.14 a
G2C2001.77 ± 52.61 a368.14 ± 29.02 c36.52 ± 4.62 a337.07 ± 34.51 a2.21 ± 0.18 a0.71 ± 0.045 a2.78 ± 0.28 a
G2T1710.63 ± 33.31 b548.41 ± 30.66 a22.74 ± 4.09 b284.09 ± 4.26 ab1.59 ± 0.08 b0.74 ± 0.035 a2.80 ± 0.00 a
Means with the same letter in the same column are not significantly different from each other according to Tukey’s test (p < 0.05). Error lines represent ± standard deviation of the mean. Different letters (a, b, c) indicate statistically significant differences according to Tukey’s test (p < 0.05).
Table
Table 5. Mean values of some micro-elements and heavy metals (mg kg−1), measured on fruits bulked from ten plants of two tomato cultivars per treatment: Narvik (G1) and Gružanski zlatni (G2), without (C) or with T. harzianum application (T). Concentrations of the same elements in the soil (b) and Bioaccumulation Index (BI) (c) are presented.
Table 5. Mean values of some micro-elements and heavy metals (mg kg−1), measured on fruits bulked from ten plants of two tomato cultivars per treatment: Narvik (G1) and Gružanski zlatni (G2), without (C) or with T. harzianum application (T). Concentrations of the same elements in the soil (b) and Bioaccumulation Index (BI) (c) are presented.
MnFeCuZnNiCrPbCo
Fruit (a)
G1C7.5 ± 0.1 c34.6 ± 1.0 b9.33 ± 0.15 a12.61 ± 0.61 a1.32 ± 0.05 a0.56 ± 0.005 b0.032 ± 0.01 a0.05 ± 0.01 a
G1T7.1 ± 0.05 d42.79 ± 0.11 a10.79 ± 1.55 ab12.79 ± 0.45 a1.08 ± 0.0 b0.89 ± 0.09 a0.0176 ± 0.03 b0.069 ± 0.02 a
G2C7.8 ± 0.1 b22.52 ± 1.00 c8.85 ± 0.02 ab11.16 ± 0.04 b0.995 ± 0.005 b0.48 ± 0.09 b0.030 ± 0.004 a0.05 ± 0.0005 a
G2T8.05 ± 0.1 a38.77 ± 2.81 a8.11 ± 0.36 b11.45 ± 0.79 ab0.49 ± 0.1 c0.99 ± 0.005 a0.0257 ± 0.045 a0.051 ± 0 a
Soil (b)
496.2 ± 3.913,685 ± 176.7743.55 ± 1.2139.13 ± 1.2023.69 ± 0.5923.07 ± 0.9913.77 ± 0.949.043 ± 0.01
BI (c)
G1C1.510.2521.4232.235.572.430.230.55
G1T1.430.3224.7832.694.563.860.120.76
G2C1.570.1720.3228.524.202.080.210.55
G2T1.620.2918.6229.262.074.290.190.56
Means with the same letter in the same column are not significantly different from each other according to Tukey’s test (p < 0.05). Error lines represent ± standard deviation of the mean. Different letters (a, b, c, d) indicate statistically significant differences according to Tukey’s test (p < 0.05).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vukelić, I.D.; Prokić, L.T.; Racić, G.M.; Pešić, M.B.; Bojović, M.M.; Sierka, E.M.; Kalaji, H.M.; Panković, D.M. Effects of Trichoderma harzianum on Photosynthetic Characteristics and Fruit Quality of Tomato Plants. Int. J. Mol. Sci. 2021, 22, 6961. https://doi.org/10.3390/ijms22136961

AMA Style

Vukelić ID, Prokić LT, Racić GM, Pešić MB, Bojović MM, Sierka EM, Kalaji HM, Panković DM. Effects of Trichoderma harzianum on Photosynthetic Characteristics and Fruit Quality of Tomato Plants. International Journal of Molecular Sciences. 2021; 22(13):6961. https://doi.org/10.3390/ijms22136961

Chicago/Turabian Style

Vukelić, Igor D., Ljiljana T. Prokić, Gordana M. Racić, Mirjana B. Pešić, Mirjana M. Bojović, Edyta M. Sierka, Hazem M. Kalaji, and Dejana M. Panković. 2021. "Effects of Trichoderma harzianum on Photosynthetic Characteristics and Fruit Quality of Tomato Plants" International Journal of Molecular Sciences 22, no. 13: 6961. https://doi.org/10.3390/ijms22136961

APA Style

Vukelić, I. D., Prokić, L. T., Racić, G. M., Pešić, M. B., Bojović, M. M., Sierka, E. M., Kalaji, H. M., & Panković, D. M. (2021). Effects of Trichoderma harzianum on Photosynthetic Characteristics and Fruit Quality of Tomato Plants. International Journal of Molecular Sciences, 22(13), 6961. https://doi.org/10.3390/ijms22136961

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop