Next Article in Journal
Research on Multimodal Prediction of E-Commerce Customer Satisfaction Driven by Big Data
Previous Article in Journal
Cryptanalysis of Ateniese–Steiner–Tsudik-Authenticated Group Key Management Protocol
Previous Article in Special Issue
Innovative Functional Lactic Acid Bacteria Fermented Oat Beverages with the Addition of Fruit Extracts and Lyophilisates
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Three-Dimensional Analysis of the Impact of Different Concentrations of Glyphosate on the Growth of Cocoa (Theobroma cacao)

by
Juan Diego Valenzuela-Cobos
1,
Fabricio Guevara-Viejó
1,
Purificación Galindo-Villardón
1,2 and
Purificación Vicente-Galindo
1,2,3,*
1
Centro de Estudios Estadísticos, Universidad Estatal de Milagro (UNEMI), Milagro 091050, Ecuador
2
Department of Statistics, University of Salamanca, 37008 Salamanca, Spain
3
Institute for Biomedical Research of Salamanca (IBSAL), 37007 Salamanca, Spain
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(18), 8180; https://doi.org/10.3390/app14188180
Submission received: 21 May 2024 / Revised: 5 July 2024 / Accepted: 22 July 2024 / Published: 11 September 2024

Abstract

:
Ecuadorian cocoa possesses important organoleptic characteristics such as aroma and flavor, called fine and aromatic cocoa. The objective of this study was to evaluate the initial growth responses of young cocoa seedlings to glyphosate in a dose progression in 45 cocoa plants (5 months old), which were transplanted into pots with substrate adjusted to pH 6.0–6.5. Glyphosate doses (0 to 904 g e.e. ha−1) were applied every two weeks, evaluating the impact at 30 and 60 days post-application. Data on shikimate accumulation parameters, chlorophyll content and PSII quantum efficiency were subjected to multivariate analysis using a three-dimensional scatter plot. The results indicated that high concentrations of glyphosate contributed to higher shikimate concentration and lower PSII quantum efficiency. The findings for the variables crop damage, stem height and stem diameter were evaluated by ANOVA. Similarities were reported between the results of the variables height and diameter, and significant differences (p < 0.05) in the variable crop damage for all treatments were also reported. In terms of phytotoxic reaction and growth parameters, the most efficient treatment was DO4, since the seedlings with this dosage showed a low percentage of damage (10%) and the best indices in terms of height and diameter. The least efficient treatment was D15. The control plants (DO1) showed a crop damage of >50% because the absence of control favored weed proliferation. These indications highlight the need to adjust glyphosate doses according to the specific needs of each crop and the development stage of the plant in order to reduce negative effects and maximize potential benefits.

1. Introduction

Cocoa (Theobroma cacao) is a tree native to countries in America, Africa, Asia and Oceania located in equatorial and tropical latitudes which have a tropical-humid environment [1]. The fruit of the tree consists of a cob containing between 20 and 50 seeds with highly valued organoleptic characteristics, such as aroma and flavor. These kernels are the main ingredient in the chocolate industry and are used in the production of a wide range of manufactured products, including cosmetics, beverages and pharmaceuticals [2]. In economic terms, cocoa is an agricultural product of vital relevance for the world market. According to ICCO [3], during the period 2019–2020 world production totaled 4.7 million tons. Of this annual production, 18.4% (0.9 million tons) was generated in Latin America, with the Ivory Coast (2.1 million tons) and Ghana (0.8 million tons) leading production, followed by Ecuador (0.32 million tons). In Ecuador, cocoa cultivation is the economic mainstay of many rural families that are dedicated solely to this activity, especially in the foothills of the Andes Mountains [4,5]. According to a report issued by the National Finance Corporation (CFN) [6], the harvested area increased by 13% between 2013 and 2016, contributing 12% to the national agricultural Gross Domestic Product (GDP) and 1.50% to the total GDP in 2019.
A limiting factor in agricultural production is the presence of so-called “weeds”. These weeds act at the same tropical level as crops, competing for soil nutrients, water and solar energy, and serve as a refuge for thriving pests. In plantations of Theobroma cacao, allelopathy secreted by weeds has the potential to decrease total bean yield and quality by up to 30%. To control the growth of these unwanted hosts, cultural techniques, such as tillage, and chemical methods, such as glyphosate, have been implemented, of which the use of herbicides has shown greater efficacy [7,8].
Glyphosate, a broad-spectrum herbicide, acts as the predominant chemical component in numerous systemic herbicides intended for the management of a wide range of annual and perennial plant species [9]. This compound regulates weed growth by interfering with the production of aromatic amino acids, which are essential for the synthesis of essential proteins in the primary and secondary metabolic processes of the plants affected by its action [9,10]. Glyphosate is currently the most widely used herbicide worldwide and is considered the safest for human health under the technical conditions indicated: direct spraying on weeds, use in concentrations between 1 and 3% and low doses and low exposure time for the person handling the herbicide [11]. Scientists and regulatory authorities globally, including the US EPA, have evaluated the carcinogenic potential of glyphosate since 1985 in animals and concluded that it does not exhibit genotoxicity within the suggested dose range. For this reason, the use of the pesticide continues to demonstrate safety for the health of farmers, while increasing crop yield and quality, generating higher profits for farmers and the agricultural industry [12,13].
Although pesticides help protect the cocoa plant, their indiscriminate use can cause damage such as stunting or even death. Adu-Yeboah et al. (2023) demonstrated that high doses of glyphosate result in significant accumulation of shikimate and decreased photosynthesis, which affects the health and growth of cocoa [14]. On the other hand, Kayode Olufemi et al. (2020) observed that moderate doses may be less harmful, suggesting that glyphosate, applied with caution, could be part of a weed management regime without significant detrimental effects [15]. In Ecuador, there is little scientific research on the collateral impact of glyphosate application for weed control on cocoa plantations of commercial interest. In this context, it is important to examine the developmental responses of Ecuadorian cocoa seedlings in the presence of different doses of glyphosate, using multivariate statistical analysis tools to reveal complex relationships between the parameters under study, which cannot be visualized with the naked eye.
The multivariate data analysis revolution began to spread in the 1990s. Data mining techniques, such as three-dimensional scattering algorithms, involve the synthesis of large volumes of information in a lower-dimensional space. The graphical representation is carried out in a two-dimensional environment expressing points for observations and vectors for variables. The main advantage of this type of approach, in contrast to traditional statistical methods of classification/sorting, is the localization and understanding of patterns, trends and sequences hidden in the data. The application of three-dimensional scatter plots in agronomic data promises to positively transform decision making in agriculture, leading to significant improvements in the efficiency of weed control treatments and, consequently, in productivity [16].
The main objective of this paper was to implement a three-dimensional dispersion analysis as a statistical visualization tool to evaluate the influence of different glyphosate-based treatments on 45 cocoa plants (five months old).

2. Materials and Methods

2.1. Seedling Selection and Sample Preparation

This experimental study was conducted at the “Esmeraldas” farm, located in the city of Montalvo, province of Los Ríos, Ecuador. Five-month-old Theobroma cacao plants of the CCN-51 variety grown under in situ conditions on this property were used for this purpose. Forty-five seedlings were selected and distributed in a sampling area of approximately 400 m2, all free of visible diseases, nutritional deficiencies or any other physical anomaly. The characteristics considered for selection were similarity in size, height, structure and general physical condition. The selected seedlings were labeled and registered for follow-up and subsequent analysis.
Seedlings were transplanted in 7 L pots with a substrate composed of a mixture of 40% compost, 30% coconut fiber, 20% coarse sand and 10% vermiculite. This substrate was adjusted to a pH of 6.0–6.5 by adding agricultural lime. Samples of the substrate (500 g) were collected for the determination of physicochemical parameters under standardized methodology. The quantification of organic matter was performed following the Walkley–Black method, whose approach is based on oxidation (potassium dichromate) and acidity (sulfuric acid), followed by titration with ferrous sulfate. Phosphorus availability was extracted by Olsen’s method, using sodium bicarbonate, molybdate reagents and ascorbic acid, with absorbance at 800 nm. The number of exchangeable cations (K, Ca, Mg) was measured by the ammonium acetate method, while the exchangeable acidity (H+Al) was calculated by titration after extraction with 1 N KCl. The sum of exchangeable bases (SB) resulted from the total amount of exchangeable basic cations K, Ca and Mg. The effective cation exchange capacity (T) resulted from the sum of SB and H+Al. Finally, base saturation (V) was calculated as the ratio of the sum of exchangeable bases to the effective cation exchange capacity. The report revealed the following contents: organic matter (M.O.) = 25.0 g.∙dm−3, P (resin) = 30.0 mg∙dm−3, K = 2.0 mmolc∙dm−3, Ca = 25.0 mmolc∙dm−3, Mg = 15.0 mmolc∙dm−3, H+Al = 20.0 mmolc∙dm−3, SB = 42.0 mmolc∙dm−3 and T = 62.0 mmolc∙dm−3, with a base saturation (V) of 68% [17].

2.2. Treatments and Experimental Design

The experiment consisted of the direct spraying of different doses of glyphosate (Vargas brand, Quito, Ecuador) on transplanted cocoa seedlings. A total of 15 treatments were used, and the doses of each were assigned randomly, combining amounts implemented in previous studies using glyphosate, ranging from 0 to 904 g e.e. ha−1. They were applied by foliar spraying, simulating an accidental supply and ensuring foliar coverage. The process was repeated with a frequency of 15 days for 2 months. Plant evaluations were carried out in two instances: on Days 30 and 60 after the first herbicide application [18].

2.3. Measurement of Variables

At the end of the experiment (60 days), the following variables were measured: phytotoxic crop damage reactions (%) and plant biochemical reactions such as shikimate accumulation (μg/g). To assess plant photosynthetic health, chlorophyll content (mg/m²) and PSII quantum efficiency (Fv/Fm) were calculated. Subsequently, growth parameters such as plant height (m) and stem diameter (cm) were measured. These variables allowed a comprehensive evaluation of the effect of glyphosate on various physiological and growth aspects of cocoa seedlings. The data were used to confirm the homogeneity of the seedlings collected and to correlate the results of the analyses with the physical characteristics of the plants.

2.3.1. Measurement of Crop Damage

In young cocoa seedlings, visual evaluations of phytotoxicity caused by glyphosate were carried out 60 days after application. A herbicide tolerance scale from 1 to 9 developed by the European Weed Research Council (EWRC) was used for this purpose. The scale ranked plant response as follows: (1) no visible effects; (2) very mild effects, such as stunting or dwarfing and chlorosis; (3) mild effects with stunting and chlorosis, but these effects can be reversed; (4) moderate effects of chlorosis and/or stunting, with most being possibly reversible; (5) severe chlorosis and/or stunting; (6) more severe damage; (7) increased severity of damage; (8) even greater severity of damage and (9) plant death [19].

2.3.2. Photochemical Activity of the PSII

To evaluate quantum efficiency, PSII chlorophyll fluorescence was measured using a portable fluorometer (PAM 2100, WALZ, Effeltrich, Germany). Measurements were performed on fully expanded leaves of the fifteen cocoa seedlings selected for each dose. The slides were adapted to a photon flux density (FFD) of 1000 μmol∙m−2 s−1, and instantaneous maximum (F’m) and stable (F) fluorescence measurements were recorded, following the methodology of Genty et al. [20].

2.3.3. Shikimate Accumulation

To evaluate the accumulation of Shikimate, plant material was crushed in a mortar with liquid nitrogen to a fine powder, then weighed and stored at −80 °C. To determine the amount of shikimic acid by HPLC, 900 μL of 0.25 N HCl was placed in a microcentrifuge tube with 100 mg of crushed plant tissue and allowed to stand for 5 min. Subsequently, the extracts were vortexed for 5 min and ultrasonication was performed for a period of 8 min at 25 °C (2210 Branson Ultrasonicator, Markham, OH, Canada). The cellular components were then separated by centrifuging at maximum speed in a benchtop microcentrifuge for 10 min. The sediment was discarded and the clear liquid was passed through a 0.22 μm syringe filter (MilliporeSigma™ Millex®, San Jose, CA, USA) for direct HPLC analysis of shikimate according to the methodology of Kretzmert et al. [21].

2.3.4. Chlorophyll Content

Leaf tissue samples were collected in the morning between 8:00 and 10:00 am to ensure consistency in moisture content and metabolic activity. The collected leaves were immediately placed in sealed plastic bags and stored on ice for transport to the laboratory. Fifty grams of sample were used. They were crushed with 80% acetone in a biker until a homogeneous mixture was obtained. The mixture was centrifuged to separate the liquid extract (containing chlorophyll) from the solid residues. The extract was analyzed with a spectrophotometer. The absorbance was measured at different specific wavelengths of 666 nm, 649 nm and 470 nm using 95% ethanol as a target.
Chlorophyll content in cocoa seedlings was determined following Arnon’s equations by measuring the absorbance of the sample at different wavelengths [22]. The total chlorophyll content, as seen in Equation (3), is the sum of chlorophyll a (Chl a) content (Equation (1)) and chlorophyll b (Chl b) content (Equation (2)). In the formulas, A666 refers to the absorbance measured at 666 nm and A649 to the absorbance measured at 649 nm:
c a = 13.77 A 666 4.63 A 649
Equation (1). Chlorophyll a (Chl a)
c b = 24.63 A 649 6.66 A 666
Equation (2). Chlorophyll b (Chl b)
c t = 7.11 A 649 + 20 A 666
Equation (3). Total Chlorophyll Content

2.4. Statistical Analysis

A three-dimensional scattering analysis was performed to track the patterns between biochemical reactions (shikimate accumulation) and photosynthetic health parameters (chlorophyll content, PSII quantum efficiency) of cocoa plants after application of different doses of glyphosate. Multivariate analysis was performed with R software, version 4.1.1 (R Core Team, Vienna, Austria). In addition, the results of phytotoxic reaction (crop damage) and growth parameters (height, diameter) were studied by analysis of variance (ANOVA) to evaluate the significance of differences at the p < 0.05 level. Subsequently, Duncan’s test was applied (p < 0.05). This process was carried out in the software (SPSS ver. 26.0.0.0.0).

2.5. Three-Dimensional Scatter Plot

Three-dimensional scatter plots are an extension of two-dimensional scatter plots and are used to identify possible correlations between observations and variables in a data set on a single plot of three mutually perpendicular axes, where each variable is interpreted as a Cartesian coordinate vector. Cartesian coordinate vectors are defined such that each element of the matrix is the inner product of the vectors corresponding to its rows and columns [23].
The analysis of a three-dimensional scatter plot goes beyond data visualization: it also involves the understanding of different statistical tools involved, such as three-dimensional Euclidean space and measures of central tendency. The exploration begins with the definition of a three-dimensional Euclidean space, with a set of the tern (x, y, z) of real numbers, formally defined as R 3 . The observations represent the relationship between the three variables, where f is a function of two variables, z = f(x,y).
The Euclidean distance between two observations (Equation (4)) indicates the distance between the influence of two different treatments P( p 1 , q 1 , r 1 ) and Q ( p 1 , q 1 , r 1 ) [24].
d ( P , Q ) = ( x 2 x 1 ) 2 + ( y 2 y 1 ) 2 + ( z 2 z 1 ) 2
Equation (4). Euclidean distance fitted to a three-dimensional plane.
The calculation of the mean of each variable (Equation (5)) allows us to understand the concentration of the data in each dimension and to detect possible biases in the data [25].
Mean for variable “xa
a ¯ = 1 n   i = 1 n a i
Equation (5). Average of the variables in a three-dimensional plane.
Where:
  • a ¯ represents the mean of the variable.
  • a i represents each individual observation of the variable.
  • n is the total number of observations.
  • a can be x, y, z depending on the variable of interest.
Covariance (Equation (6)) measures the linear relationship between two variables, indicating the degree to which these two change together. Given a positive covariance, the variables tend to increase; given a negative covariance, one variable tends to increase while the other variable tends to decrease and given a covariance close to zero, a weak linear relationship between the variables is inferred.
Covariance between variable “a” and variable “b”.
C o v a , b = 1 n 1   i = 1 n a i a ¯ b i b ¯
Equation (6). Calculation of covariance in a three-dimensional plane.
Where:
  • a ¯ and b ¯ represents the mean of the variable.
  • a i and b i are the individual observations of the variables a and b , respectively.
  • n is the total number of observations.
  • a and b can be any pair of the variables x, y, or z to compute C o v x , y , C o v x , z or C o v y , z .
Pearson’s correlation (Equation (7)) is another measure of association, which allows us to denote how related two variables are to each other [26]. For a three-dimensional plane, the Pearson correlation is calculated in pairs.
Correlation between variable “a” and variable “b”.
r a , b = C o v   ( a , b ) σ a σ b  
Equation (7). Calculation of Pearson’s correlation in a three-dimensional plane.
Where:
  • r a , b represents the correlation coefficient between the variables a and b .
  • C o v ( a , b ) is the covariance between the variables a and b .
  • n are the standard deviations of the variables a and b , respectively.
  • a and b can be any pair of the variables x, y, or z to compute r x , y , r x , z or r y , z .

3. Results and Discussion

The purpose of this study was to expose, by means of a three-dimensional scatter plot and ANOVA, the complex effects of different doses of glyphosate on various parameters of cocoa seedling development. The treatments used in the research, with their respective doses, were the following: DO1 (0 g e.e. ha−1); DO2 (25 g e.e. ha−1); DO3 (50 g a.e. ha−1); DO4 (80 g e.e. ha−1); DO5 (144 g e.e. ha−1); DO6 (288 g e.e. ha−1); DO7 (325 g e.e. ha−1); DO8 (432 g e.e. ha−1); DO9 (535 g e.e. ha−1); DO10 (625 g e.e. ha−1); DO11 (700 g e.e. ha−1); DO12 (768 g e.e. ha−1); DO13 (803 g e.e. ha−1); DO14 (864 g e.e. ha−1) and DO15 (904 g e.e. ha−1).
Three plants were treated for each treatment and the characteristics of each group were evaluated after the application of the doses.

3.1. Multivariate Statistics for the Analysis of Photosynthetic Health and Biochemical Reaction of Seedlings

Figure 1 expresses the three-dimensional scatter plot that relates to the patrons existing between the parameters of photosynthetic health and biochemical reaction. Plane 1, corresponding to the horizontal variable, represents the accumulation of shikimate; Plane 2, belonging to the vertical variable, denotes the quantum efficiency of PSII and Plane 3, corresponding to the depth variable, collects information about chlorophyll content. The Euclidean distance allowed us to quantify how far apart the different treatments are according to the combination of their effects on the plants. The graph reduced the complexity of the multivariate information and reflects the normal distribution of 15 points, concerning the number of treatments, around the three Cartesian coordinates (x, y and z).
In this analysis, we observed that high concentrations of glyphosate, particularly from DO5 treatment, showed higher shikimate accumulation and lower chlorophyll content, contributing to lower quantum efficiency. On the other hand, lower concentrations of glyphosate, from DO1 to DO4, exhibited the highest ranges of PSII quantum efficiency and chlorophyll content. Plant performance responses were not perfectly positively correlated with applied doses, which may be attributed to the influence of other factors, such as genetic variability and the interaction of allelopathy-producing microorganisms, among others. Deterioration in quantum efficiency was found to be associated with a reduction in chlorophyll content, which implies a decrease in the plant’s ability to capture and use light energy effectively. As a cluster analysis, it was determined that the DO4 (80 g e.e. ha−1) treatment was the one that showed the highest balance among the evaluated parameters, emphasizing the higher PSII quantum efficiency and the higher chlorophyll content.
Similar studies showed that glyphosate can interfere with the enzyme EPSPS, which is essential for amino acid synthesis in plants, which could be a contributing factor to the observed effects on cocoa growth [27]. In another study, De Maria [28] indicated that the phytotoxic effect on tissues affects young leaves more due to the degree of increased shikimate in the tissues. The herbicidal activity of glyphosate lies in its ability to inhibit the EPSPS enzyme, blocking the shikimate pathway [29]. This pathway is essential for the biosynthesis of aromatic compounds in plants. In the study of Mueller et al. [30], shikimate was found to accumulate in glyphosate-resistant and glyphosate-sensitive Conyza canadensis plants. This phenomenon suggests that shikimate accumulation may not be directly related to glyphosate resistance, indicating that other mechanisms may be involved in the resistance of some plants to this herbicide [31].
For Maxwell and Johnson [32], the quantum efficiency of photosynthesis, especially the PSII parameter (Φ), is crucial to evaluate the response of plants to abiotic stresses such as drought, salinity and temperature extremes. These environmental conditions are important to determine and control for effective plant development. Quantum efficiency is fundamental for the design of lighting systems in controlled agriculture, such as greenhouses or vertical farming systems, where artificial light is used to maximize photosynthesis. In another study, Mateos-Naranjo et al. [33] indicated that exposure to glyphosate can cause alterations in the ultrastructure of chloroplasts in eucalyptus plants, which could negatively influence quantum efficiency. Glyphosate can inhibit chlorophyll synthesis by interfering with the production of key precursors, leading to a decrease in the ability of plants to capture light energy and effectively photosynthesize. Other studies showed that herbicide application to soybean significantly reduced chlorophyll content, which is related to a lower rate of photosynthesis and compromised plant growth [34]. Chlorophyll is not only essential for light absorption during photosynthesis but also plays a key role in protecting plant tissues against photooxidative damage. As chlorophyll concentration decreases, plants suffer a reduction in biomass and their overall capacity for growth and development.
A similar study by Zobiole et al. [35] observed that soybean plants treated with glyphosate showed a significant reduction in PSII quantum efficiency and a decrease in chlorophyll concentration. On the other hand, the study of Adu-Yeboah et al. [14] indicated that a glyphosate dosage of 360 g a.e. ha−1 or more resulted in >60% leaf damage and a more than 10-fold increase in shikimate accumulation. On the other hand, in the study of Adu-Yeboah et al. [14], whose experimentation was carried out under field conditions in Ghana, the sensitivity of cocoa plantations to spraying with different doses of glyphosate was also indicated. The authors found that glyphosate doses ≥360 g a.e. ha−1 caused leaf damage to more than 60% of plant structures, a greater than 10-fold increase in shikimate accumulation and a significant decrease in chlorophyll content. These findings highlight the need to develop and promote the use of herbicides with more specific active ingredients in order to minimize the collateral impact on non-target plants.

3.2. Analysis of Variance (ANOVA) for the Evaluation of Phytotoxic Reaction and Growth Parameters (Height, Diameter)

Table 1 shows the results of the investigation with respect to phytotoxic reaction and growth parameters. The findings indicated that there were similarities between treatments in terms of crop damage, height and stem diameter. Higher concentrations of the herbicide inhibited plant development, while lower concentrations showed fewer adverse effects. The means obtained allowed establishing that the most efficient treatment was DO4, using a dose of 80 g e.a. ha−1 since the seedlings treated with this dosage showed a low percentage of damage and the best indices in terms of height and diameter. The least efficient treatment was D15 because the seedlings with this dosage showed a lower development and the second-highest percentage of damage to the crop.
Leaf damage was measured in percentages, showing necrosis, chlorosis and wilting of the plant. Cocoa plants belonging to the control treatment (DO1) showed a crop damage of >50% because the absence of control favored the proliferation of weeds. A similar result occurred with the study of Ríos [31]: the control proved to be lower due to the high density of weeds competing for light, nutrients and space. The doses of 803 and 904 g e.e. ha−1 caused leaf scorch in young leaves, as well as changes in the shape of leaf margins and depigmentation along the leaflets. After application of doses of 144 and 288 g a.e. ha−1, cocoa seedlings exhibited depigmentation on older leaves, abnormal coloration, wavy blade-shaped edges and blighted leaves but produced healthy new leaves 60 days after application. Although the symptoms were worrisome, the plants recovered.
The study by Martínez Ramírez et al. [32] indicated that the highest dose of 0.200 kg∙-ha−1 caused significant phytotoxic effects. Initially, phytotoxicity symptoms, such as chlorosis, were severe, reaching a grade 4 but then the symptoms decreased to a grade 3 after 45 days of application and were completely eliminated after 60 days, leaving no residual effects on the treated plant population. This observation is important because it shows how the dose of a chemical can affect plant health in the short term, but also indicates long-term recovery if conditions permit [33].
A similar study was conducted by Konrad et al. [36], in which they observed that low doses of glyphosate did not interfere with plant growth for the BRS-Maués cultivar. In addition, the morphology of plants exposed to high doses of this herbicide was also affected, with a reduction in stem height and diameter. Helander, et al. [37] observed that in trees and shrubs the herbicide can significantly reduce both height and stem diameter, especially in young plants. This is due to inhibition of essential amino acid synthesis and general damage to photosynthetic capacity.
In contrast, Magalhães Farias et al. [38] found that insufficient g a.e. ha−1 of glyphosate on H. brasiliensis caused damage by reducing plant development in height and diameter. However, they recovered the growth of these plants within 180 days after application. In the statistical analysis of Neves et al. [39], they observed a 15% increase in the height variable at doses of 18–33.5 g a.e. ha−1 and maximum responses of 16% at doses of 27 g a.e. ha−1. Glyphosate is a broad-spectrum herbicide that can function as a plant growth regulator, provided that the correct doses are applied with precision. Similarly, Silva et al. [40] reported that, at low doses, the herbicide can stimulate an increase in plant height, coinciding with the results found in this study.
Barbosa et al. [41] found that this herbicide on coffee seedlings can reduce their growth, with more pronounced effects when the herbicide is absorbed by the leaves. These effects include reductions in height, root length and total dry mass. In the study by Gonçalves et al. [42], high doses of glyphosate were found to cause severe intoxication in coffee seedlings, resulting in significant reductions in height, leaf area, and shoot and root biomass. However, the use of foliar fertilizers can mitigate the negative effects of these doses. It is important to know the dose-response effect of the plant to glyphosate because extremely high doses affect foliar and root uptake. This study could be relevant to improve management and dosage practices in agricultural applications to minimize plant damage and control pests or diseases.

4. Conclusions

In terms of photosynthetic health, high rates of shikimate accumulation, specifically from DO5 treatment, contributed to a detrimental effect on seedlings, contributing to lower quantum efficiency. Meanwhile, lower concentrations of glyphosate, from DO1 to D04, exhibited higher PSII quantum efficiency and higher chlorophyll content.
In terms of phytotoxic reaction and growth parameters, the most efficient treatment was DO4 (80 g e.e. ha−1), since the seedlings with this dosage showed a low percentage of damage (10%) and the best indices in terms of height and diameter. The least efficient treatment was D15.
The cocoa plants belonging to the control treatment (DO1) showed a crop damage of >50% due to the fact that the absence of control favored the proliferation of weeds; this allows us to deduce that the application of glyphosate is necessary.

Author Contributions

Conceptualization, J.D.V.-C. and F.G.-V.; formal analysis, J.D.V.-C. and F.G.-V.; investigation, F.G.-V.; methodology, J.D.V.-C. and F.G.-V.; supervision, P.G.-V. and P.V.-G.; writing—original draft, J.D.V.-C., F.G.-V., P.G.-V. and P.V.-G.; writing—review and editing, P.G.-V. and P.V.-G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Universidad Estatal de Milagro (UNEMI) Scholarship.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Acknowledgments

The authors are grateful to the Universidad Estatal de Milagro (UNEMI).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Muller, M.W.; Valle, R.R. Ecofisiologia do cultivo do cacaueiro. In Ciência, Tecnologia e Manejo do Cacaueiro, 2nd ed.; Valle, R.R., Ed.; Centro de Pesquisas do Cacau Ceplac: Brasília, Brasil, 2012; pp. 31–66. [Google Scholar]
  2. Sánchez Mora, F.D.; Garcés Fiallos, F.R. Moniliophthora roreri (Cif y Par) Evans et al. en el cultivo de cacao. Sci. Agropecu. 2012, 3, 249–258, Universidad Nacional de Trujillo. Available online: https://www.redalyc.org/articulo.oa?id=357633703006 (accessed on 30 July 2024). [CrossRef]
  3. ICCO. Producción de Granos de Cacao. La Organización Internacional del Cacao. 2021. Available online: https://www.icco.org/about-us/international-cocoa-agreements/doc_download/4577-production-qbcs-xlvi-no-2.html (accessed on 30 July 2024).
  4. Vargas, E.M.; Molina, X.C.; Cevallos, E.Z. Recorrido Histórico de La Importancia Del Cacao Para La Economía de Ecuador; Sinergias Educativas: Ecuador, Quevedo, 2022. [Google Scholar]
  5. Cervantes, X.; Mendoza, E.; Vásconez, G.; Samaniego, E. Visualización Bibliométrica En El Área Del Cultivo de Cacao Evaluando Competencias Administrativas y Socioeconómicas a Nivel de América Latina. 2020. Available online: https://www.japss.org/upload/33.%20Molina.pdf (accessed on 1 May 2024).
  6. CFN. Ficha Sectorial. 2023. Available online: https://www.cfn.fin.ec/wpcontent/uploads/downloads/biblioteca/2023/fichas-sectoriales-2-trimestre/Ficha-SectorialCacao.pdf (accessed on 6 May 2024).
  7. Marín, C.; Weiner, J. Efectos de la densidad y el patrón de siembra en la supresión de malezas y el rendimiento de grano en tres variedades de maíz bajo alta presión de malezas. Investig. Malezas 2014, 54, 467–474. [Google Scholar] [CrossRef]
  8. Oerke, E.-C. Crop losses to pests. J. Agric. Sci. 2006, 144, 31–43. [Google Scholar] [CrossRef]
  9. Carlisle, S.M.; Trevors, J.T. Glyphosate in the environment. Water Air Soil Pollut. 1988, 39, 409–420. [Google Scholar] [CrossRef]
  10. U.S. Department of Agriculture Forest Service (Servicio Forestal de los Estados Unidos). Glifosato: Perfil de Información Sobre Herbicidas; Boletín; Administración de Medicamentos de los Estados Unidos, Región Noroeste del Pacífico: Washington, DC, USA, 1997; p. 25.
  11. Paz-y-Miño, C.; López Cortés, A. Impactos Genéticos, Sanitarios y Ambientales Causados por el Glifosato en Ecuador. Editorial Nova. 2014. Available online: http://site.ebrary.com/id/10935200 (accessed on 1 May 2024).
  12. Williams, G.M.; Aardema, M.; Acquavella, J.; Berry, S.C.; Brusick, D.; Burns, M.M.; de Camargo, J.L.V.; Garabrant, D.; Greim, H.A.; Kier, L.D.; et al. A Review of the Carcinogenic Potential of Glyphosate by Four Independent Expert Panels and Comparison to the IARC Assessment. Crit. Rev. Toxicol. 2016, 46, 3–20. [Google Scholar] [CrossRef]
  13. Lara, M.; Zubiate, G. Control de plaguicidas en Ecuador: ¿un problema subestimado? Rev. Bionatura 2020. [Google Scholar] [CrossRef]
  14. Adu-Yeboah, P.; Lowor, S.T.; Segbefia, M.A.; Konlan, S.; Pobee, P. Physiological and growth responses of cacao to glyphosate exposure. Rev. Cienc. Ambient. Salud Parte B 2023, 58, 91–99. [Google Scholar] [CrossRef]
  15. Olufemi, A.K.; Babatunde, A.K.; Olufemi, O.O.; Adewale, A.S.; Osasogie, U.; Mohammed, I.; Akanji, A.S. Evaluation of herbicides and manual weed control methods in the establishment of Theobroma cacao L. in Cross River State of Nigeria. Mundo J. Adv. Res. Rev. 2020, 8, 298–306. [Google Scholar] [CrossRef]
  16. Guevara-Viejó, F.; Valenzuela-Cobos, J.D.; Vicente-Galindo, P.; Galindo-Villardón, P. Application of K-Means Clustering Algorithm to Commercial Parameters of Pleurotus spp. Cultivated on Representative Agricultural Wastes from Province of Guayas. J. Fungi 2021, 7, 537. [Google Scholar] [CrossRef]
  17. Carvalho, L.B.D.; Alves, P.L.C.A.; Duke, S.O. Hormesis with glyphosate depends on coffee growth stage. An. Acad. Bras. Ciênc. 2013, 85, 813–822. [Google Scholar] [CrossRef]
  18. Solano Neto, A.; Yamashita, O.M.; Carvalho, M.A.C.D.; Felito, R.A.; Rocha, A.M.D.; Rabelo, H.D.O. Efeito de doses de glyphosate em deriva simulada em plantas jovens de cacau. Principia 2021, 1, 134. [Google Scholar] [CrossRef]
  19. Estimado, B.S.; Sandral, G.A.; Spencer, D.; Khan, M.R.I.; Higgins, T.J.V. The tolerance of three transgenic subterranean clover (Trifolium subterraneum L.) lines with the bxn gene to herbicides containing bromoxynil. Aust. J. Agric. Res. 2003, 54, 203. [Google Scholar] [CrossRef]
  20. Genty, B.; Briantais, J.-M.; Panadero, N.R. The relationship between the quantum yield of photosynthetic electron transport and quenching of chlorophyll fluorescence. Biochim. Biophys. Acta (BBA)-Gen. Subj. 1989, 990, 87–92. [Google Scholar] [CrossRef]
  21. Kretzmert, K.A.; Sammons, D.; Frantz, J.P.; Morris, P.J. A Quick and Colorful Method to Detect and Quantify Shikimic Acid in Whole Plant Tissues or Plant Extracts. In Weed Technology; Weed Science Society of America: Westminster, CO, USA, 2011; pp. 78–83. [Google Scholar]
  22. Del Valle Tascón, S.; Sanz, M.J.; Calatayud, A.; Barreno, E. Coeficientes de Extinción de Clorofilas y Feofitinas (a y b) En DMSO y Ecuaciones Para El Cálculo de Sus Concentraciones. In Stvdia Botánica; Departamento de Biología Vegetal, Universitat de Valencia: Valencia, Spain, 1994; pp. 115–121. [Google Scholar]
  23. Henzi, A.; Derecho, M. A rank-based sequential test of independence. Biometrika 2024, asae023. [Google Scholar] [CrossRef]
  24. Herrera, E. La genética de poblaciones y el origen de la diversidad humana. Rev. Med. Hondureña 2013, 81, 40–45. [Google Scholar]
  25. Walpole, R.; Myers, R.; Myers, S.; Ye, K. Probabilidad y Estadística Para Ingeniería y Ciencias, 9th ed.; Pearson Educación: Naucalpan de Juárez, México, 2012. [Google Scholar]
  26. Vargas, A. Estadística Descriptiva e Inferencial, 2nd ed.; Universidad de Castilla-La Mancha: Murcia, Spain, 1995. [Google Scholar]
  27. Antolinez Sandoval, E.Y.; Almanza Merchán, P.J.; Barona Rodriguez, A.F.; Polanco Díaz, E.; Serrano Cely, P.A. Current State of Cocoa Plantation: A Review of its Main Limitations. Cien. Agri 2020, 17, 1–11. [Google Scholar] [CrossRef]
  28. De Maria, N. Efecto del Glifosato Sobre la Simbiosis “Lupinus albus-Bradyrhizobium” sp. (“Lupinus”); Complutense de Madrid: Madrid, Spain, 2004. [Google Scholar]
  29. García-Villanueva, L.A.; Cuapio-Ortega, V.H.; Henández-Paniagua, I.Y.; Fernández-Villagómez, G.; Rodrigo-Ilarri, J.; Rodrigo-Clavero, M.E.; Andraca-Ayala, G.L.; Hernández-Cruz, G.B.; Banda-Santamaría, S. Effects of Glyphosate on the Environment and Human Health. Nat. Environ. Pollut. Technol. 2024, 23, 17–32. [Google Scholar] [CrossRef]
  30. Mueller, T.C.; Massey, J.H.; Hayes, R.M.; Main, C.L.; Stewart, C.N. Shikimate Accumulates in Both Glyphosate-Sensitive and Glyphosate-Resistant Horseweed (Conyza canadensis L. Cronq.). J. Agric. Food Chem. 2003, 51, 680–684. [Google Scholar] [CrossRef]
  31. Ríos, A. Control de Malezas y Fitotoxicidad al Maíz Transgénico RR Con Dos Formulaciones de Glifosato. Ph.D. Dissertation, Escuela Agrícola Panamericana, Tegucigalpa, Honduras, 2011. [Google Scholar]
  32. Maxwell, K.; Johnson, G.N. Chlorophyll fluorescence—A practical guide. Rev. Botánica Exp. 2000, 51, 659–668. [Google Scholar] [CrossRef]
  33. Mateos-Naranjo, E.; Pérez-Martín, A. Efectos de las concentraciones subletales de glifosato en el crecimiento y el rendimiento fotosintético de especies no objetivo Bolboschoenus Maritimus. Quimiosfera 2013, 93, 2631–2638. [Google Scholar]
  34. Martínez Ramírez, R.; Zuaznábar Zuaznábar, R.; Rodríguez García, J. Efectividad En El Control de Malezas y Fitotoxicidad de Guateque GD 75 En Caña de Azúcar. Rev. Ing. Agrícola 2022, 12, 60–64. Available online: https://www.redalyc.org/articulo.oa?id=586272858008 (accessed on 1 May 2024).
  35. Zobiole, L.H.S.; Oliveira, R.S.; Kremer, R.J.; Constantin, J.; Bonato, C.M.; Muñiz, A.S. Eficiencia En El Uso Del Agua y Fotosíntesis de La Soja Resistente al Glifosato Afectada Por El Glifosato. Bioquímica Fisiol. Pestic. 2010, 97, 182–193. [Google Scholar]
  36. Konrad, B.N.L.; Pinheiro, S.C.; Ferreira, C.C.; Hoffmann, E.K.; Albertino, S.M.F. Using Glyphosate on Guarana Seedlings in the Amazon. Moléculas 2023, 28, 5193. [Google Scholar] [CrossRef] [PubMed]
  37. Helander, M.; Saloniemi, I.; Saikkonen, K. Glyphosate in northern ecosystems. Trends Plant Sci. 2012, 17, 569–574. [Google Scholar] [CrossRef] [PubMed]
  38. Magalhães Farias, C.C.; Rondón Neto, R.M.; Magalhães Farias, C.B. Efeitos de subdoses de glyphosate em plantas jovens de seringueira (Hevea brasiliensis Aubl.). Rev. Bras. Herbic. 2012, 11, 119–125. [Google Scholar]
  39. Neves, D.C.; Furlani Junior, E.; Valter Filho, V.V. Hormese no crescimento do algodoeiro por subdoses de glifosato. In Congresso Brasileiro Do Algodão De Moda Anais; Embrapa Algodón: Brasília, Brazil, 2009; pp. 915–922. [Google Scholar]
  40. Silva, J.C.D.; Arf, O.; Gerlach, G.A.X.; Kuryiama, C.S.; Rodrigues, R.A.F. Efeito hormese de glyphosate em feijoeiro. Pesqui. Agropecu. Demasiado. 2012, 42, 295–302. [Google Scholar] [CrossRef]
  41. Barbosa, E.A.; Silva, I.M.D.; França, A.C.; Silva, E.D.B.; Matos, C.D.C.D. Evaluation of leaf and root absorptions of glyphosate in the growth of coffee plants. Arq. Inst. Biol. 2020, 87, E0762018. [Google Scholar] [CrossRef]
  42. Luisa Lucca Gonçalves, A.; Antônio Mendes Pereira, G.; Aparecida Silva, L.; Patrik Alves, D.; Franchini Pucci, L.; Afonso Reis Gonçalves, V.; Rodrigues dos Reis, M. Use of foliar fertilizers as attenuators of injuries caused by glyphosate drift in young coffee plantations. CS Coffee Sci. 2020, 15, 1–7. [Google Scholar] [CrossRef]
Figure 1. Three-dimensional scatter plot of biochemical reactions and photosynthetic health of cocoa seedlings as a function of different doses of glyphosate. Euclidean distance was calculated considering the following variables: shikimate accumulation, chlorophyll content and PSII quantum efficiency.
Figure 1. Three-dimensional scatter plot of biochemical reactions and photosynthetic health of cocoa seedlings as a function of different doses of glyphosate. Euclidean distance was calculated considering the following variables: shikimate accumulation, chlorophyll content and PSII quantum efficiency.
Applsci 14 08180 g001
Table 1. Effect of glyphosate application on growth parameters of cocoa seedlings.
Table 1. Effect of glyphosate application on growth parameters of cocoa seedlings.
TreatmentCrop Damage (%)Height (m)Diameter (cm)
DO123.00 b ± 2.640.66 ce ± 0.041.10 d ± 0.11
DO215.33 a ± 0.280.70 d ± 0.02 1.13 de ± 0.12
DO312.00 a ± 2.640.75 f ± 0.021.18 fg ± 0.15
DO410.00 a ± 1.13 0.83 h ± 0.051.23 g ± 0.03
DO525.00 b ± 1.560.77 fg ± 0.021.16 fg ± 0.03
DO633.00 c ± 1.730.75 eg ± 0.03 1.10 de ± 0.04
DO750.16 d ± 0.280.78 e ± 0.111.15 ef ± 0.04
DO853.41 de ± 2.960.66 de ± 0.031.05 cd ± 0.06
DO956.33 e ± 1.150.64 cd ± 0.031.00 cd ± 0.13
DO1053.50 de ± 1.560.65 cd ± 0.041.03 cd ± 0.09
DO1155.00 de ± 1.29 0.62 ac ± 0.060.95 bc ± 0.04
DO1272.33 f ± 1.510.60 bc ± 0.070.93 b ± 0.05
DO1375.66 f ± 1.080.57 bc ± 0.050.90 b ± 0.09
DO1471.66 f ± 1.880.53 ab ± 0.11 0.91 b ± 0.07
DO1575.00 f ± 1.210.50 a ± 0.070.85 a ± 0.09
Note: The results of crop damage, stem height and diameter data are presented according to each treatment applied; Different letters within a column indicate significant differences (p < 0.05).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Valenzuela-Cobos, J.D.; Guevara-Viejó, F.; Galindo-Villardón, P.; Vicente-Galindo, P. Three-Dimensional Analysis of the Impact of Different Concentrations of Glyphosate on the Growth of Cocoa (Theobroma cacao). Appl. Sci. 2024, 14, 8180. https://doi.org/10.3390/app14188180

AMA Style

Valenzuela-Cobos JD, Guevara-Viejó F, Galindo-Villardón P, Vicente-Galindo P. Three-Dimensional Analysis of the Impact of Different Concentrations of Glyphosate on the Growth of Cocoa (Theobroma cacao). Applied Sciences. 2024; 14(18):8180. https://doi.org/10.3390/app14188180

Chicago/Turabian Style

Valenzuela-Cobos, Juan Diego, Fabricio Guevara-Viejó, Purificación Galindo-Villardón, and Purificación Vicente-Galindo. 2024. "Three-Dimensional Analysis of the Impact of Different Concentrations of Glyphosate on the Growth of Cocoa (Theobroma cacao)" Applied Sciences 14, no. 18: 8180. https://doi.org/10.3390/app14188180

APA Style

Valenzuela-Cobos, J. D., Guevara-Viejó, F., Galindo-Villardón, P., & Vicente-Galindo, P. (2024). Three-Dimensional Analysis of the Impact of Different Concentrations of Glyphosate on the Growth of Cocoa (Theobroma cacao). Applied Sciences, 14(18), 8180. https://doi.org/10.3390/app14188180

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