Open-access Multivariate relationships in strawberry cultivated with native communities of arbuscular mycorrhizal fungi

ABSTRACT.

The mechanisms underlying the interactions between native mycorrhizal fungal communities and strawberry plants remain unclear. However, the identification of specific associations among variables and their influence on the total experimental variability when using inoculants based on arbuscular mycorrhizal fungi should enable the identification of the most relevant ones. Herein, our objective was to identify and characterize variables related to each other and to the total experimental variability among strawberry plants inoculated with native mycorrhizal communities. Experimental treatments included an uninoculated control and eight multi-specific inoculants from cultivated soils and native forests from reference strawberry-cultivation sites (Bom Princípio, Flores da Cunha, Ipê, and São José of Hortêncio) in Rio Grande do Sul State, Brazil. Morphological, productivity, and quality traits were evaluated. Inoculants obtained from agricultural ecosystems of Bom Princípio and Ipê did not influence the horticultural performance of strawberries, while those from Flores da Cunha largely explained total experimental variability, and therefore, should be considered when selecting the location to obtain inoculants for use on strawberry plants. Number of fruits, fruit flavor, chlorophyll a, and total chlorophyll contents, and, most importantly, root variables, should be included for experimental analysis of ‘Albion’ strawberry responses to multi-specific mycorrhizal inoculants from different locations.

Keywords:
X ananassa Duch.; mycorrhizal biotechnology; principal component analysis; Pearson correlation; canonical correlation.

Introduction

Planning and establishing sustainable agri-food systems are two of the greatest challenges currently faced by horticultural production, which traditionally requires an excess of potentially contaminating chemical inputs (fertilizers and biocides), thus posing a threat to agroecosystems and the environment at large. As this scenario is the opposite of sustainable development goals envisaged on a global scale, there is an urgent need to promote agroecology as a science and driver of environmentally friendly crop management practices. Among the bio-tools that can contribute to enhancing sustainability in horticulture are arbuscular mycorrhizal fungi (AMF), which are naturally present in the soil and establish mutualistic associations with the roots of more than 80% of terrestrial flora, thereby improving the acquisition of water and nutrients by the plant host.

A case in point, strawberry (Fragaria X ananassa Duch.) is an AMF-responsive horticultural crop (Chiomento et al., 2019a). Whereas hydroponic cultivation of strawberries involves large-scale use of fertilizers that contribute to environmental contamination hazards in open drainage systems, mycorrhizal biotechnology can help minimize these inconveniences (Chiomento et al., 2021a). Therefore, researchers have focused on understanding how these microorganisms benefit strawberry plants. The evidence indicates that AMF improve the development of the strawberry root system (Chiomento et al., 2021a) and activate plant defense metabolism, such that berries actually acquire higher levels of phytochemicals (Chiomento et al., 2019a).

Overall, AMF do not show strict host specificity, indicating that plant root systems are often co-colonized by multiple fungal species (Öpik et al., 2009). Therefore, the use of mycorrhizal communities appears to be a promising strategy towards more sustainable horticultural models. Furthermore, plant development is enhanced by indigenous or native AMF (Oliveira et al., 2017) because these microorganisms are physiologically and genetically better adapted to the edaphoclimatic and biogeographical conditions of the local agroecosystem than exotic AMF (Faye et al., 2013).

Indeed, the use of native AMF reportedly improves the adaptation of native hosts to nursery and field conditions (Maltz & Treseder, 2015). However, how mycorrhizal communities interact with strawberry plants in relation to their horticultural potential remains unclear. Principal component analysis (PCA) is a convenient way to determine the influence of the relationships among all the variables studied. This multivariate-covariance structure modeling technique aims to identify latent variables that represent linear combinations of a group of related variables.

The identification of specifically strong associations between variables and their influence on total experimental variability makes it possible to highlight those associations that are most relevance in each case. Thus, for example, while studying strawberry plants, Chiomento et al. (2021b) and Chiomento et al. (2023) adopted statistical analysis techniques that enabled them to describe the relationships among variables and generate other important information to guide reliable practical recommendations. However, information regarding the use of multivariate techniques to understand the effects of native AMF communities on strawberry horticultural potential is scarce at best. Thus, here, we aimed to identify and characterize variables closely related to each other and to the total variability in an experiment in which the strawberry, cultivar ‘Albion’, was grown together with native AMF communities.

Material and methods

Plant material

The study was conducted under hydroponic cultivation in a greenhouse in the municipality of Passo Fundo (28º15'46" S; 52º24'24" W), Rio Grande do Sul State, Brazil.

The day-neutral strawberry cultivar ‘Albion’ was used. The strawberry plants used in the experiment were produced by mother plants obtained from the nursery Llahuén/Chilean Patagonia (33º50'15.41" S; 70º40'03.06" W). In June 2016, the mother plants were transplanted into containers (5 L) filled with the commercial substrate Horta 2® and kept on 1.20 m high benches. In September 2016, the ends of the runners produced were removed from the matrices and transferred to a 72-cell tray (50 cm3 per cell) filled with sterilized sand (120ºC) for 20 min.

Experimental design

The experiment was laid in a completely randomized design with three replicates of a single plant per treatment. Treatments included eight multi-specific mycorrhizal inoculants (AMF communities) (Table 1) collected from soils at reference sites for strawberry cultivation in Rio Grande do Sul State, Brazil (Chiomento et al., 2019b), and a control (uninoculated) treatment. At each reference site, soils were collected from strawberry cultivation sites (SC) and native forest soils (NF). The collected inoculants were identified based on spore morphology and ontogeny.

Table 1
Description of AMF communities identified in soils of reference sites for strawberry cultivation in Rio Grande do Sul State, Brazil.

Procedures

In October 2016, strawberry daughter plants were transplanted into 9 L polyethylene pots filled with sterilized sand. For the treatments inoculated with AMF, 10 g of inoculant soil collected at the reference sites for strawberry cultivation was added to the planting hole of daughter plants at the time of transplanting. A localized irrigation system was used to water the plants individually via dripper rods (2.4 L h-1 per emitter). Irrigation was applied four times a day (total wetting for 10 min.). Additionally, fert-irrigation was performed weekly (Furlani & Fernandes Júnior, 2004).

Measurements

The evaluated attributes included, 1) morphology, 2) productivity, and 3) quality characteristics.

The morphological traits were:

Crown diameter (CD, cm): measured with a digital caliper;

Height of the aerial plant body (HAP, cm) was measured using a ruler;

Crown number (CN);

Fresh mass of the aerial plant body (FMAP, g) and fresh mass of the root system (FMRS, g) were measured on an electronic analytical scale;

Dry mass of the aerial plant body (DMAP, g) and dry mass of the root system (DMRS, g). Tissues were oven-dried at 65ºC for 48h until constant mass and weighed on an electronic analytical scale;

Total root length (TRL, cm), root surface area (RSA, cm2), and root volume (RV, cm3) were analyzed through images using the WinRHIZO® software;

For measurement of TRL, roots were subdivided into the following diameter classes: very fine roots (VFR, Ø < 0.5 mm), fine roots (FR, Ø from 0.5 to 2 mm), and coarse roots (CR, Ø > 2 mm). TRL was determined according to the methodology of Böhm (1979);

Accumulation of dry mass in the aerial plant body (ADMAP, %) was determined according to Atif et al. (2016) using the following equation: ADMAP= (DMAP/FMAP)×100 (1), where DMAP and FMAP are the dry and fresh masses of the aerial plant body, respectively;

Accumulation of dry mass of the root system (ADMRS, %) was determined according to Atif et al. (2016), using the following equation: ADMRS= (DMRS/FMRS)×100 (2), where DMRS is the dry mass of the root system and FMRS is the fresh mass of the root system.

Productivity attributes were measured and calculated using the means of all harvests. The fruits were harvested when they were 85% red. The measured traits were:

Number of fruits (NF);

Fruit weight (FW, g): fruits were weighed on an electronic analytical scale.

As for quality attributes, these included total soluble solids content (SSC), total titratable acidity (TTA), and total anthocyanin content (TAC). One-hundred-gram fruits samples were used. The quality characteristics were determined as follows:

The Dickson quality index (DQI) was determined according to Dickson et al. (1960) using the equation DQI =TDMHAPCD+DMAPDMRS (3), where TDM = total dry mass, HAP = height of the aerial plant body, CD = crown diameter, DMAP = dry mass of the aerial plant body, and DMRS = dry mass of the root system;

Mycorrhizal colonization (MC, %) was determined using the following equation: MC = TNFMR/TNF×100 (4), following the methodologies of Phillips and Hayman (1970), and Trouvelot et al. (1986), where TNFMR is the total number of fragments with mycorrhizal roots, and TNF is the total number of fragments;

Chlorophyll a (CLA), b (CLB), and total (CLT) chlorophyll were measured with an electronic ClorofiLOG 1030 (Falker, Garches, France) chlorophyll-o-meter; 15 measurements were made in the central third of the leaf blade on each sampled plant;

Total soluble solids content (SSC, %) were determined with an analog refractometer;

Total titratable acidity (TTA, % citric acid): determined according to Zenebon et al. (2008);

Fruit flavor (FF): FF= SSC/TTA (5), where SSC is the total soluble solids and TTA is the total titratable acidity.

Total anthocyanin content (TAC, mg of pelargonidin-3-O-glycoside equivalent per 100 g of fresh fruit [mg PE 100 g-1 FF]) was determined using the differential pH method (Lee et al., 2005).

Statistical analysis

Statistical analysis began with the estimation of the means and standard deviations for each variable in each treatment for descriptive characterization of the observations. Subsequently, PCA was applied to identify the variables most related to the treatments, in addition to identifying those that accounted the most for the experimental variability.

After applying PCA, Pearson’s correlation analysis was performed, generating a matrix of correlation coefficients among the treatments and PCA dimensions, and among the measured variables and PCA dimensions. Subsequently, the number of conditions was obtained by the ratio between the highest and lowest eigenvalues of the X’X correlation matrix. A number of conditions ≤100, indicated weak multicollinearity occurrence; a number of conditions between 100 and 1,000 indicated moderate to severe multicollinearity; and a number of conditions ≥1,000, indicated severe multicollinearity (Montgomery et al., 2012). The variance inflation factor was obtained for each variable on the inverse diagonal of the correlation matrix X’X. Severe multicollinearity resulted when the variance inflation-factor value was >10 (Hair et al., 2009). The occurrence of multicollinearity among the explanatory variables was defined by obtaining values for the number of conditions ≥1,000 and variance inflation-factor values >10. When moderate, strong, or high multicollinearity was detected, the variables that caused these results were removed from the dataset and a subsequent diagnosis was made to prove the effectiveness of removing variables from the dataset submitted to statistical analysis.

Finally, a canonical correlation analysis was performed to identify the existence of significant associations among the three groups of variables after identifying those with the lowest contributions to the main components, namely: 1) morphological variables (CN, FMAP, DMAP, FMRS, DMRS, ADMRS, TRL, RSA, VFR, FR, and CR); 2) productive variables (NF and FW); and 3) quality variables (SSC, TTA, FF, TAC, DQI, CLA, CLB, and CLT).

All statistical analyzes were performed in the R software (R Core Team, 2022) with a 5% error probability.

Results

The greatest variability, characterized by the highest standard deviations, occurred for FMRS, DMRS, VFR, FR, FW, TRL, RSA, and RV (Table 2). These results were expected because the evaluations of root systems typically show inherent high variability owing to the proper preparation of the roots for further evaluation. When checking variable means, considering the two mycorrhizal community ecosystems of origin i.e., cropping and forestland, the means of inoculants from agricultural ecosystems (SC) tended towards lower means than those from natural ecosystems (NF) for morphological variables and quality. The opposite was observed for production variables and those of the root system (Table 2).

Considering all the variables measured in our experiment, the results of PCA indicated that the first two components accounted for 41.30% of the total variability observed (Figure 1A), with FMRS, FF, ADMRS, NF, and DMRS accounting for the greater positive contributions, and VFR, FR, TRL, CLA, and CLT accounting for the greater negative contributions in the first two dimensions (Table 3, Figure 2A and B). These findings indicate that these variables were the most relevant in the analysis.

Table 2
Mean and standard deviation of morphological, productive, and quality variables of strawberry cultivar ‘Albion’, inoculated with mycorrhizal communities.

Figure 1
Biplot of principal component analysis with all variables (A) and with the removal of variables ADMAP, MC, RV, HAP, CD (B), considering the first two principal components. CN: crown number; CD: crown diameter (cm); HAP: shoot height (cm); FMAP: shoot fresh mass (g); DMAP: shoot dry mass (g); FMRS: fresh mass of the root system (g); DMRS: dry mass of the root system (g); ADMAP: accumulation of shoot dry mass (%); ADMRS: accumulation of dry mass of the root system (%); VFR: very fine roots (cm); FR: fine roots (cm); CR: coarse roots (cm), NF: number of fruits; FW: fruit weight (g); SSC: total soluble solids content (%); TTA: titratable total acidity (% of citric acid); FF: fruit flavor; TAC: total anthocyanin content (mg PE 100 g-1 FF); TRL: total root length (cm); RSA: root surface area (cm²); RV: root volume (cm³); MC: mycorrhizal colonization (%); DQI: Dickson quality index; CLA: chlorophyll A; CLB: chlorophyll B; CLT: total chlorophyll content.

With this first PCA, we identified that ADMAP, MC, RV, HAP, and CD had greater contributions (positive or negative) to the last three components (Table 3 and Figure 2A). Thus, these five variables were removed from the database and PCA was repeated in their absence (Table 4), which increased the contribution of the first two components to 55.10% (Figure 1B). This finding further supported the contribution of the variables to the first two components (Table 4 and Figure 2B). This approach underlines to importance of highlighting and interpreting FMRS, FF, ADMRS, NF, DMRS, VFR, FR, TRL, CLA, and CLT to further explain the responses of strawberry cultivar ‘Albion’ to multispecific mycorrhizal inoculants from different locations and ecosystems.

By correlating the nine treatments with the PCA dimensions, the multispecific inoculants from Flores da Cunha in an agricultural ecosystem, Flores da Cunha in a natural forest ecosystem, and São José do Hortêncio in an agricultural ecosystem were most related to the first two dimensions, explaining 80.70% of the total variability (Figure 2C). Inoculants from Bom Princípio and Ipê in natural forest ecosystems were the most related to the third and fourth PCA dimensions, which explained a small proportion (6.20%) of the total variability. Similarly, neither Bom Princípio or Ipê inoculants in agricultural ecosystems, São José do Hortêncio in a natural forest ecosystem, or the AMF-uninoculated control exhibited relevant relationships with any PCA dimensions (Figure 2C). These results indicate that the AMF communities from Flores da Cunha stood out as the greatest explanation for the total experimental variability and should therefore be considered relevant for selecting the place of origin of those communities, regardless of the ecosystem.

Table 3
Contribution of all analyzed variables in the experiment in each dimension of the principal components.

Further, the biplot analysis showed that the AMF communities from Flores da Cunha were directly related to a greater extent to TRL, VFR, and FR, and to a lesser extent to RSA, CR, TTA, and TAC (Figure 3). The AMF community from São José do Hortêncio in an agricultural ecosystem, which correlated with the second dimension of the PCA (Figure 2C), was directly related to a larger extent to SSC, DQI, and DMRS and, to a lesser extent, to ADMRS, NF, FW, and CN.

The inoculants formed by the mycorrhizal communities of Bom Princípio in a natural forest ecosystem, Ipê in a natural ecosystem, São José do Hortêncio in a natural ecosystem, and the uninoculated control showed close associations with FF, FMAP, DMAP, CLA, CLB, and CLT (Figure 3). Consistently, mycorrhizal inoculants obtained from the agricultural ecosystems of Bom Princípio and Ipê, which did not show significant relationships with any PCA dimensions (Figure 2C) did not show a close relationship with any measured variable either (Figure 3), indicating that neither inoculant influenced the horticultural performance of strawberries.

Figure 2
Pearson’s linear correlation between the variables and dimensions of the PCA with all variables (A), with the removal of variables ADMAP, MC, RV, HAP, CD (B) and between the nine treatments (1: Bom Princípio in agricultural ecosystem; 2: Bom Princípio in a natural ecosystem; 3: Flores da Cunha in an agricultural ecosystem; 4: Flores da Cunha in a natural ecosystem; 5: Ipê in an agricultural ecosystem; 6: Ipê in a natural ecosystem; 7: São José do Hortêncio in an agricultural ecosystem 8: São José do Hortêncio in a natural ecosystem; 9: control without AMF inoculation) and the PCA dimensions (C). CN: crown number; CD: crown diameter (cm); HAP: shoot height (cm); FMAP: shoot fresh mass (g); DMAP: shoot dry mass (g); FMRS: fresh mass of the root system (g); DMRS: dry mass of the root system (g); ADMAP: accumulation of shoot dry mass (%); ADMRS: accumulation of dry mass of the root system (%); VFR: very fine roots (cm); FR: fine roots (cm); CR: coarse roots (cm), NF: number of fruits; FW: fruit weight (g); SSC: total soluble solids content (%); TTA: titratable total acidity (% of citric acid); FF: fruit flavor; TAC: total anthocyanin content (mg PE 100 g-1 FF); TRL: total root length (cm); RSA: root surface area (cm²); RV: root volume (cm³); MC: mycorrhizal colonization (%); DQI: Dickson quality index; CLA: chlorophyll A; CLB: chlorophyll B; CLT: total chlorophyll content.

Table 4
Contribution of selected variables in each principal component dimension.

The number of conditions was classified as severe multicollinearity, with values of 1.5e17 and 3.6e17 for the groups of morphological and quality variables, respectively, and a variance inflation factor greater than 30 for the same groups. Based on this result, we removed VFR, FR, CR, FF, TTA, RSA, and CLT from the database, in addition to those already indicated by PCA (ADMAP, MC, RV, HAP, and CD). Thus, the number of conditions ranged from 9.72 to 40.81 (weak multicollinearity) and the variance inflation factor from 1.70 to 7.90, within the maximum limit of 10.00.

After the creation of the three groups, a new Pearson correlation analysis was conducted among the variables of each group (Figure 4). Thus, only trivial relationships were significant for morphological variables such as shoot dry mass and shoot fresh mass (Figure 4A). This result was also obtained by correlating the number and weight of fruits in the productive variable group (Figure 4B). No significant correlations were observed in the group of quality traits (Figure 4C).

Figure 3
Principal component analysis for the variables and their association with the nine treatments studied (1: Bom Princípio in an agricultural ecosystem; 2: Bom Princípio in a natural ecosystem; 3: Flores da Cunha in an agricultural ecosystem; 4: Flores da Cunha in a natural ecosystem; 5: Ipê in an agricultural ecosystem; 6: Ipê in a natural ecosystem; 7: São José do Hortêncio in an agricultural ecosystem; 8: São José do Hortêncio in a natural ecosystem; 9: control without AMF inoculation). CN: crown number; CD: crown diameter (cm); HAP: shoot height (cm); FMAP: shoot fresh mass (g); DMAP: shoot dry mass (g); FMRS: fresh mass of the root system (g); DMRS: dry mass of the root system (g); ADMAP: accumulation of shoot dry mass (%); ADMRS: accumulation of dry mass of the root system (%); VFR: very fine roots (cm); FR: fine roots (cm); CR: coarse roots (cm), NF: number of fruits; FW: fruit weight (g); SSC: total soluble solids content (%); TTA: titratable total acidity (% of citric acid); FF: fruit flavor; TAC: total anthocyanin content (mg PE 100 g-1 FF); TRL: total root length (cm); RSA: root surface area (cm²); RV: root volume (cm³); MC: mycorrhizal colonization (%); DQI: Dickson quality index; CLA: chlorophyll A; CLB: chlorophyll B; CLT: total chlorophyll content.

Figure 4
Pearson’s linear correlation within groups of morphological (CN, FMAP, DMAP, FMRS, DMRS, ADMRS, and TRL), productive (NF and FW), and quality (SSC, TTA, DQI, CLA, and CLB) variables. CN: crown number; FMAP: shoot fresh mass (g); DMAP: shoot dry mass (g); FMRS: fresh mass of the root system (g); DMRS: dry mass of the root system (g); ADMRS: root system dry mass accumulation (%); NF: number of fruits; FW: fruit weight (g); SSC: total soluble solids content (%); TTA: titratable total acidity (% of citric acid); TRL: total root length (cm); DQI: Dickson quality index; CLA: chlorophyll A; CLB: chlorophyll B.

To confirm the possible relationships of the variables among the three groups, canonical correlation analysis was performed, which showed no significant effects when relating the groups of morphological variables with those of the productivity variables (Table 5) or the latter with those of the quality variables (Table 6).

Table 5
Crossed canonical loads between morphological and productive variables.
Table 6
Crossed canonical loads between productive and quality variables.

These results show that the number and weight of fruits were not correlated with the morphological and quality variables of strawberries. However, when correlating the groups of morphological variables with those of quality variables (Table 7), the results were significant (p-value < 0.05) for the first three canonical pairs. Therefore, we found that higher DMAP provided fruits with lower TTA (1st pair), while higher root fresh mass provided fruits with higher TTA and lower TRL (2nd pair), and higher root dry mass and root dry mass accumulation provided higher DQI and fruits with lower chlorophyll A content.

Additionally, we observed that a higher DMRS will provide a higher DQI (1st and 3rd pair), and when a reduction in DMRS occurs, there will be a reduction in DQI (2nd pair). Finally, the relationship between a decrease in CN and a lower DQI was verified in the second canonical pair (Table 7).

Table 7
Cross canonical loads between morphological and quality variables.

Discussion

Before and after the selection of variables, the PCA approach underlined the need to analyze and interpret FMRS, DMRS, ADMRS, FF, NF, FR, VFR, TRL, CLA, and CLT to explain the responses of the strawberry cultivar ‘Albion’ to multi-specific mycorrhizal inoculants from different locations and ecosystems. Specifically, PCA makes it possible to simultaneously reduce the dimensionality of a dataset and preserve its variability (Jolliffe & Cadima, 2016) because it is statistically coherent, computationally fast, and scalable (Price et al., 2010). The positive effect of mycorrhizal inoculation on the root system of the plant host (Chiomento et al., 2019a; Chiomento et al., 2021a) is attributed to molecular signaling between symbionts through the release of lipo-chito-oligosaccharides by AMF, which stimulates root formation in the host plant (Oláh et al., 2005).

The correlations between experimental treatments and PCA dimensions allowed us to verify that the multi-specific inoculants from Flores da Cunha in NF and SC ecosystems were directly related to TRL, VFR, FR, RSA, CR, TTA, and TAC. In contrast, the inoculant from the SC in São José was directly related to SSC, DQI, DMRS, ADMRS, NF, FW, and CN. For these ecosystems, 80.70% of the total variability was explained by the first two components, which agreed well with the general rule that at least 70.00% of the total variance in the PCA must be explained by the first two PCs (Rencher, 2002).

By and large, the contribution of mycorrhizal communities from Flores da Cunha to accounting for the total experimental variability observed was the greatest, and should be considered relevant for the selection of the place to obtain the inoculant, especially in the case of agricultural ecosystems. Strawberry crops were established at this soil-inoculant collection site from 2009 to 2016, using cultivar ‘Albion’ (Chiomento et al., 2019b), the same cultivar used in our study. This indicated that the fungal species collected from the NF and SC communities (Table 1) had the highest affinity for ‘Albion’ Because AMF are associated with a wide range of plants, clear host specificity is uncommon, although preferences for plant symbionts arising from effector proteins secreted by AMF to manipulate host cells and facilitate successful infection have been reported (van der Heijden et al., 2015). This alters the structure of the host, suppresses innate immune responses, and alters plant metabolism (Zeng et al., 2018).

Mycorrhizal inoculants obtained from agricultural ecosystems in Bom Princípio and Ipê did not show close relationships with any measured variable. Therefore, they did not influence the horticultural performance of strawberries. As the agricultural ecosystems of Ipê and Bom Princípio were affected by anthropization four and six years ago, respectively, and FC SC and SH SC had been modified for a longer time (Chiomento et al., 2019b), we believe that the mycorrhizal communities BP SC and IP SC were less adapted to strawberry monoculture. Despite the lack of AMF-host specificity, a substantial functional diversity can modulate the benefits generated by microorganisms, which include the plant symbiont, fungal species, and their relationships with environmental conditions (Koch et al., 2017).

When we inoculated strawberry plants of cultivar ‘Albion’ with AMF communities from natural ecosystems (NF), we observed a better performance of this horticultural crop in relation to morphological variables and fruit quality. Thus, we suggest that F. mosseae improves plant growth and fruit quality, as it is a coincident species in communities from NF. The best performance of AMF inoculants from SC was observed for root system and fruit production variables. In this case, C. etunicatum and Glomus sp. (caesaris like) seem to promote strawberry root development and productivity, because they were the two species with equal abundance in the communities of SC.

In Brazil, only one AMF-based commercial inoculant is available, consisting of a single species [Rhizophagus intraradices (formerly Glomus intraradices)]. Traditionally, commercial inoculants are composed of only one fungal species (monospecific) that may not adapt to the conditions of a given agroecosystem. As the use of mycorrhizal species compatible with the plant symbiont intensifies mutualism, we emphasize that the multi-specific inoculants used here originated from soil adapted to strawberry cultivation (Table 1) (Chiomento et al., 2019b). In grassland ecosystems, the mycorrhizal plant root system is colonized by a community of several fungal species (Maherali & Klironomos, 2012). Researchers and industry must focus on studying the potential of multispecific inoculants composed of native fungal species as bio-tools for the sustainable plant cropping.

The inoculation of plants with AMF native to soils under the cultivation of this plant species contributes to improving its performance and adaptation to nursery and field conditions (Maltz & Treseder, 2015), presumably because indigenous AMF are genetically and physiologically adapted to their native hosts (Oliveira et al., 2017). Thus, for example, the restoration of Picconia azorica (Tutin) knobl, native forests was facilitated by the use of AMF, which was more appropriate to the ecological conditions of the environment and the host (Melo et al., 2019). As for strawberries, the addition of native mycorrhizal communities associated with the use of biochar reportedly increases mycorrhizal colonization and improves the development of the root system (Chiomento et al., 2021a).

In the Brazilian subtropics, indigenous AMF communities native to soils at the reference sites for strawberry cultivation (Chiomento et al., 2019b) can be used as bioinput to improve the horticultural potential of this crop in the field as well as under greenhouse conditions. Our findings indicate that the selection of multi-specific mycorrhizal inoculants enhances their positive effects on plant hosts through the maximum benefit of the symbiotic association. Prior to developing commercial inoculants, mycorrhizal genera, species, and strains must be studied for their potential to improve host development. As we verified the dynamic behavior of inoculants depending on their origin, such as those from Flores da Cunha, we believe that careful pre-selection of this bio-tool is indispensable for high-yielding strawberry crop production.

Conclusion

Principal component analysis applied before and after selection of variables, revealed that the number of fruits, fruit flavor, chlorophyll a and total chlorophyll contents, and, mainly, the characteristics of the root system, must be included in the experimental analyzes for further explanation of the responses of the strawberry cultivar ‘Albion’. Particularly, AMF communities obtained from Flores da Cunha, Rio Grande do Sul State, Brazil, explained the greater proportion of total experimental variability and, therefore, should be relevant for selecting the place of origin of the AMF inoculants for use in strawberry crops, as they contain the fungal species with greater affinity for cultivar ‘Albion’ In contrast, the mycorrhizal inoculants obtained from the agricultural ecosystems of Bom Princípio and Ipê, Rio Grande do Sul State, Brazil, did not influence the horticultural performance of strawberry plants, probably owing to the recent intensive anthropization of these agroecosystems.

Acknowledgements

To Bioagro Comercial Agropecuária Ltda. for supplying bare-root strawberry mother plants. AMF community used in this study is regulated by Sistema Nacional de Gestão do Patrimônio Genético e do Conhecimento Tradicional Associado (SisGen) of the Ministry of the Environment, Brazil, according to the registration number A198F50.

References

  • Atif, M. J., & Jellani, G., & Humair, M., & Ahmed, H., & Saleem, N., & Ullah, H., & Khan, M. Z., & Ikram, S. (2016). Different growth media effect the germination and growth of tomato seedlings. Science, Technology & Development, 35(3), 123-127. https://doi.org/10.3923/std.2016.123.127
    » https://doi.org/https://doi.org/10.3923/std.2016.123.127
  • Böhm, W. (1979). Methods of studying root systems. Springer-Verlag.
  • Chiomento, J. L. T., & Costa, R. C., & De Nardi, F. S., & Trentin, N. S., & Nienow, A. A., & Calvete, E. O. (2019a). Arbuscular mycorrhizal fungi communities improve the phytochemical quality of strawberry. The Journal of Horticultural Science and Biotechnology, 94(5), 653-663. https://doi.org/10.1080/14620316.2019.1599699
    » https://doi.org/https://doi.org/10.1080/14620316.2019.1599699
  • Chiomento, J. L. T., & Stürmer, S. L., & Carrenho, R., & Costa, R. C., & Scheffer-Basso, S. M., & Antunes, L. E. C., & Nienow, A. A., & Calvete, E. O. (2019b). Composition of arbuscular mycorrhizal fungi communities signals generalist species in soils cultivated with strawberry. Scientia Horticulturae, 253, 286-294. https://doi.org/10.1016/j.scienta.2019.04.029
    » https://doi.org/https://doi.org/10.1016/j.scienta.2019.04.029
  • Chiomento, J. L. T., & De Nardi, F. S., & Filippi, D., & Trentin, T. S., & Dornelles, A. G., & Fornari, M., & Nienow, A. A., & Calvete, E. O. (2021a). Morpho-horticultural performance of strawberry cultivated on substrate with arbuscular mycorrhizal fungi and biochar. Scientia Horticulturae, 282, https://doi.org/10.1016/j.scienta.2021.110053
    » https://doi.org/https://doi.org/10.1016/j.scienta.2021.110053
  • Chiomento, J. L. T., & De Nardi, F. S., & Kujawa, S. C., & Deggerone, Y. S., & Fante, R., & Kaspary, I. J., & Dornelles, A. G., & Huzar-Novakowiski, J., & Trentin, T. S. (2023). Multivariate contrasts of seven strawberry cultivars in soilless cultivation and greenhouse in southern Brazil. Advanced Chemicobiology Research, 2(1), 62-78. https://doi.org/10.37256/acbr.2120232332
    » https://doi.org/https://doi.org/10.37256/acbr.2120232332
  • Chiomento, J. L. T., & Lima Júnior, E. P., & D’Agostini, M., & De Nardi, F. S., & Trentin, T. S., & Dornelles, A. G., & Huzar-Novakowiski, J., & Calvete, E. O. (2021b). Horticultural potential of nine strawberry cultivars by greenhouse production in Brazil: A view through multivariate analysis. Scientia Horticulturae, 279, https://doi.org/10.1016/j.scienta.2020.109738
    » https://doi.org/https://doi.org/10.1016/j.scienta.2020.109738
  • Dickson, A., & Leaf, A. L., & Hosner, J. F. (1960). Quality appraisal of white spruce and white pine seedling stock in nurseries. The Forestry Chronicle, 36(1), 10-13. https://doi.org/10.5558/tfc36010-1
    » https://doi.org/https://doi.org/10.5558/tfc36010-1
  • Faye, A., & Dalpé, Y., & Ndung’u-Magiroi, K., & Jefwa, J., & Ndoye, I., & Diouf, M., & Lesueur, D. (2013). Evaluation of commercial arbuscular mycorrhizal inoculants. Canadian Journal of Plant Science, 93(6), 1201-1208. https://doi.org/10.4141/cjps2013-326
    » https://doi.org/https://doi.org/10.4141/cjps2013-326
  • Furlani, P. R., & Fernandez Júnior, F. (2004). Cultivo hidropônico de morango em ambiente protegido. Embrapa Clima Temperado.
  • Hair, J. F., & Black, W. C., & Babin, B., & Anderson, R., & Tatham, R. L. (2009). Análise multivariada de dados. Bookman.
  • Jolliffe, I. T., & Cadima, J. (2016). Principal component analysis: A review and recent developments. Philosophical Transactions of the Royal Society A, 374(2065), 1-16. https://doi.org/10.1098/rsta.2015.0202
    » https://doi.org/https://doi.org/10.1098/rsta.2015.0202
  • Koch, A. M., & Antunes, P. M., & Maherali, H., & Hart, M. M., & Klironomos, J. N. (2017). Evolutionary asymmetry in the arbuscular mycorrhizal symbiosis: Conservatism in fungal morphology does not predict host plant growth. New Phytologist, 214(3), 1330-1337. https://doi.org/10.1111/nph.14465
    » https://doi.org/https://doi.org/10.1111/nph.14465
  • Lee, J., & Durst, R. W., & Wrolstad, R. E. (2005). Determination of total monomeric anthocyanin pigment content of fruit juices, beverages, natural colorants, and wines by the pH differential method: collaborative study. Journal of AOAC International, 88(5), 1269-1278. https://doi.org/10.1093/jaoac/88.5.1269
    » https://doi.org/https://doi.org/10.1093/jaoac/88.5.1269
  • Maherali, H., & Klironomos, J. N. (2012). Phylogenetic and trait-based assembly of arbuscular mycorrhizal fungal communities. PLoS One, 7(5), 1-9. https://doi.org/10.1371/journal.pone.0036695
    » https://doi.org/https://doi.org/10.1371/journal.pone.0036695
  • Maltz, M. R., & Treseder, K. K. (2015). Sources of inocula influence mycorrhizal colonization of plants in restoration projects: A meta-analysis. Restoration Ecology, 23(5), 625-634. https://doi.org/10.1111/rec.12231
    » https://doi.org/https://doi.org/10.1111/rec.12231
  • Melo, C. D., & Walker, C., & Krüger, C., & Borges, P. A. V., & Luna, S., & Mendonça, D., & Fonseca, H. M. A. C., & Machado, A. C. (2019). Environmental factors driving arbuscular mycorrhizal fungal communities associated with endemic woody plant *Picconia azorica* on native forest of Azores. Annals of Microbiology, 69(11), 1309-1327. https://doi.org/10.1007/s13213-019-01535-x
    » https://doi.org/https://doi.org/10.1007/s13213-019-01535-x
  • Montgomery, D. C., & Peck, E. A., & Vining, G. G. (2012). Introduction to Linear Regression Analysis. Wiley.
  • Oláh, B., & Brière, C., & Bécard, G., & Dénarié, J., & Gough, C. (2005). Nod factors and a diffusible factor from arbuscular mycorrhizal fungi stimulate lateral root formation in *Medicago truncatula* via the DMI1/DMI2 signalling pathway. The Plant Journal, 44(2), 195-207. https://doi.org/10.1111/j.1365-313X.2005.02522.x
    » https://doi.org/https://doi.org/10.1111/j.1365-313X.2005.02522.x
  • Oliveira, J. R. G., & Resende, G. M., & Melo, N. F., & Yano-Melo, A. M. (2017). Symbiotic compatibility between arbuscular mycorrhizal fungi (autoctone or exotic) and three native species of the Caatinga in different phosphorus levels. Acta Scientiarum. Biological Sciences, 39(1), 59-69. https://doi.org/10.4025/actascibiolsci.v39i1.33486
    » https://doi.org/https://doi.org/10.4025/actascibiolsci.v39i1.33486
  • Öpik, M., & Metsis, M., & Daniell, T. J., & Zobel, M., & Moora, M. (2009). Large-scale parallel 454 sequencing reveals host ecological group specificity of arbuscular mycorrhizal fungi in a boreonemoral forest. New Phytologist, 184(2), 424-437. https://doi.org/10.1111/j.1469-8137.2009.02920.x
    » https://doi.org/https://doi.org/10.1111/j.1469-8137.2009.02920.x
  • Phillips, J. M., & Hayman, D. S. (1970). Improved procedures for clearing roots and staining parasitic and vesicular-arbuscular mycorrhizal fungi for rapid assessment of infection. Transactions of the British Mycological Society, 55(1), 158-161. https://doi.org/10.1016/S0007-1536(70)80110-3
    » https://doi.org/https://doi.org/10.1016/S0007-1536(70)80110-3
  • Price, A., & Zaitlen, N., & Reich, D., & Patterson, N. (2010). New approaches to population stratification in genome-wide association studies. Nature Reviews Genetics, 11, 459-463. https://doi.org/10.1038/nrg2813
    » https://doi.org/https://doi.org/10.1038/nrg2813
  • R Core Team (2022). R: A language and environment for statistical computing. R Foundation for Statistical Computing. https://www.R-project.org/
    » https://www.R-project.org/
  • Redecker, D., & Schubler, A., & Stockinger, H., & Stürmer, S. L., & Morton, J. B., & Walker, C. (2013). An evidence-based consensus for the classification of arbuscular mycorrhizal fungi (*Glomeromycota*). Mycorrhiza, 23, 515-531. https://doi.org/10.1007/s00572-013-0486-y
    » https://doi.org/https://doi.org/10.1007/s00572-013-0486-y
  • Rencher, A. C. (2002). Methods of multivariate analysis John Wiley & Sons.
  • Trouvelot, A., Kouch, J., & Gianinazzi-Pearson, V. (1986). Mesure du taux de mycorhization VA d’un système radiculaire: recherche of method d’estimation ayant une signification fonctionelle. In V. Gianinazzi-Pearson, & S. Gianinazzi (Eds.), Aspects physiologiques et génétiques des mycorhizes, (pp217-221). Inra Press.
  • van der Heijden, M. G. A., & Martin, F. M., & Selosse, M. A., & Sanders, I. R. (2015). Mycorrhizal ecology and evolution: The past, the present, and the future. New Phytologist, 205(4), 1406-1423. https://doi.org/10.1111/nph.13288
    » https://doi.org/https://doi.org/10.1111/nph.13288
  • Zenebon, O., & Pascuet, N. S., & Tiglea, P. (2008). Métodos físico-químicos para análise de alimentos. Instituto Adolfo Lutz.
  • Zeng, T., & Holmer, R., & Hontelez, J., & Lintel-Hekkert, B., & Marufu, L., & Zeeuw, T., & Wu, F., & Schijlen, E., & Bisseling, T., & Limpens, E. (2018). Host- and stage-dependent secretome of the arbuscular mycorrhizal fungus Rhizophagus irregularis The Plant Journal, 94(3), 411-425. https://doi.org/10.1111/tpj.13908
    » https://doi.org/https://doi.org/10.1111/tpj.13908

Publication Dates

  • Publication in this collection
    02 June 2025
  • Date of issue
    2025

History

  • Received
    15 Dec 2023
  • Accepted
    12 Apr 2024
location_on
Editora da Universidade Estadual de Maringá - EDUEM Av. Colombo, 5790, bloco 40, 87020-900 - Maringá PR/ Brasil, Tel.: (55 44) 3011-4253, Fax: (55 44) 3011-1392 - Maringá - PR - Brazil
E-mail: actaagron@uem.br
rss_feed Acompanhe os números deste periódico no seu leitor de RSS
Reportar erro