ORIGINAL RESEARCH article

Front. Toxicol., 21 May 2025

Sec. Developmental and Reproductive Toxicology

Volume 7 - 2025 | https://doi.org/10.3389/ftox.2025.1579952

This article is part of the Research TopicFrom North to South: A Global Perspective on the Impact of Endocrine Disruptors on Reproductive Development and FunctionView all articles

Residential proximity to agricultural fields, urinary glyphosate levels and breast cancer risk: a case-control study in Argentina

  • 1Instituto de Salud y Ambiente del Litoral (ISAL), Facultad de Bioquímica y Ciencias Biológicas, Universidad Nacional del Litoral (UNL) - Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Santa Fe, Argentina
  • 2Cátedra de Fisiología Humana, Facultad de Bioquímica y Ciencias Biológicas, Universidad Nacional del Litoral, Santa Fe, Argentina
  • 3Instituto de Investigaciones en Señales, Sistemas e Inteligencia Computacional, sinc(i), Facultad de Ingeniería y Ciencias Hídricas, UNL, CONICET, Santa Fe, Argentina
  • 4Departamento de Investigación Epidemiológica, Instituto Nacional de Epidemiología “Dr. Juan H. Jara”, ANLIS Malbrán, Mar del Plata, Argentina
  • 5Departamento de Matemática, Facultad de Bioquímica y Ciencias Biológicas, UNL, Santa Fe, Argentina
  • 6Servicio de Ginecología, Hospital José María Cullen, Santa Fe, Argentina
  • 7Programa de Investigación y Análisis de Residuos y Contaminantes Químicos, Facultad de Ingeniería Química, UNL, Santa Fe, Argentina

Despite accumulated evidence indicating glyphosate herbicide (GLY) presents endocrine disrupting properties, there are still discrepancies. Moreover, few epidemiological studies have focused on hormone-related pathologies. This work aimed to investigate the associations between urinary GLY levels and breast cancer (BC) in women from a region of intense agricultural activity in Argentina, exploring residential proximity to agricultural fields as a potential risk factor for BC. This was a case-control study that involved 90 women from different populations in the Province of Santa Fe, Argentina. Demographic data, lifestyle factors, and residential history were obtained through a questionnaire, while medical outcomes and reproductive history were acquired from medical records. Spot urine samples were collected and the concentrations of GLY and its primary metabolite, aminomethylphosphonic acid (AMPA) were quantified by ultra-high-performance liquid chromatography–mass spectrometry. Odds ratios were estimated to assess the strength of the association between the case/control type and each predictor. GLY concentrations were above the limit of detection (LOD) in 86.1% of samples, with a range of 0.37–10.07 µg GLY/g creatinine. AMPA was not detected in any of the samples analyzed. Although urinary GLY concentrations showed no differences between the case and control groups, women residing near agricultural fields showed an increased risk of BC (OR: 7.38, 95% CI: 2.74–21.90). These original findings show the ubiquitous presence of GLY in adult women from Argentina. Interestingly, women living near agricultural fields have a higher risk of BC, suggesting that exposure not only to GLY but also to agrochemicals in general, could predispose to the development of BC in Argentina. While this study provides valuable insights, further and broader assessments of BC distribution in relation to agrochemical exposure acroos different regions of Argentina are needed.

1 Introduction

In 2022, 2.3 million women worldwide were diagnosed with Breast Cancer (BC), resulting in 670,000 deaths (World Health Organization, 2024). In Argentina, BC is the leading cause of death in women, with 6,100 deaths annually and an estimated 22,000 new cases per year, representing 32.1% of the total cancer incidence (Programa Nacional de Control de Cáncer de Mama (PNCM, 2024). The development of BC is associated with several risk factors, both modifiable and non-modifiable. Among the modifiable risk factors are environmental factors, including chemical agents such as pesticides (Łukasiewicz et al., 2021).

N-(phosphonomethyl)glycine, also known as Glyphosate (GLY), is the active ingredient of commercial formulations named as glyphosate-based herbicides (GBHs). GLY herbicide is a systemic, broad-spectrum and post-emergent agrochemical that is used worldwide for weed control in agriculture, urban and gardening (Acquavella et al., 2004; Andreotti et al., 2018). GLY and its primary metabolite, aminomethylphosphonic acid (AMPA), have been detected in a variety of sources, including water (Aparicio et al., 2013; Bonansea et al., 2017; Mac Loughlin et al., 2017; Singh et al., 2024), soil (Primost et al., 2017; Ayoola et al., 2023), dust resulting from field erosion (Mendez et al., 2017), and food for both human and animal consumption (Rodrigues and De Souza, 2018; Vicini et al., 2021).

Biomonitoring studies have reported the presence of GLY and/or AMPA in urine samples from individuals residing in rural and/or urban areas (Connolly et al., 2022; Chang et al., 2024; Ben Khadda et al., 2025) and, even more worrying, from pregnant and lactating women (Camiccia et al., 2022; Lesseur et al., 2022; Ashley-Martin et al., 2023) and children (Grau et al., 2022; Berni et al., 2023). In addition, GLY concentrations were determined in maternal and umbilical cord serum (Kongtip et al., 2017) and breast milk (Camiccia et al., 2022; Galli et al., 2024). Specifically, in a country with high pesticide use like Argentina, the available data on human exposure to GLY herbicide are limited. To date, concentrations of GLY were documentated only in two regions of Argentina, in the urine of a small rural population from the Province of Chaco (Bressán et al., 2021) and in the urine, plasma, and saliva of an occupationally and non-occupationally exposed male population from the Province of Córdoba (Filippi et al., 2024).

The use of GLY in Argentina began to increase significantly in the mid-1990s, following the introduction of genetically modified GLY resistant soybeans (Benbrook, 2016). This agricultural transition resulted in a substantial increase in GLY application, thereby transforming Argentina into one of the leading consumers of GLY worldwide. Regulatory measures have historically been minimal (Blois and Rendón, 2023); however, growing public concern over health and environmental risks has prompted local-level restrictions and increasing calls for national regulatory reforms, although comprehensive bans or national regulations remain absent (Schmidt et al., 2022).

The utilization of GLY has been extensively implemented throughout the South American continent. In Brazil, a study reported significant contamination of drinking water with GLY and AMPA, particularly in agricultural regions like Paraná, where 100% of the municipalities analyzed exceeded the maximum limits for GLY-AMPA, correlating with increased cancer risks, including BC (Panis et al., 2022). Furthermore, research conducted in Southeast Brazil has revealed the presence of GLY concentrations in drinking water sources reaching up to 8.70 µg/L, thereby exceeding the established safety limits set by national and international regulatory bodies (Lima et al., 2022). In Colombia, in areas dedicated to agricultural activities, the presence of GLY was detected in different water sources near the crops with concentrations ranging from 2.01 to 2.77 µg/L (Alza-Camacho et al., 2016). A recent multicompartmental monitoring study in Uruguay revealed the presence of GLY and other pesticides in the surface water, sediments, and biota of Laguna del Cisne, a subtropical lake utilized as a drinking water source (Rodríguez-Bolaña et al., 2023). The present findings underscore the pervasive environmental occurrence of GLY and underscore the necessity for enhanced surveillance, thereby situating the Argentine context within a more expansive South American framework of pesticide exposure.

Association between GLY exposure and different types of cancer has been researched. Epidemiological studies conducted in Sweden, Canada, and the United States have determined a positive correlation between GLY exposure and the development of non-Hodgkin’s lymphoma (NHL) (Weisenburger, 2021) and multiple myeloma (De Roos et al., 2005). More recently, some authors in a large pooled study did not detect a relationship between GLY (active ingredient) and all types of NHL or multiple myeloma separately, but did find an association with follicular lymphoma, a subtype of NHL (De Roos et al., 2022). Indeed, the carcinogenic potential of GLY has been extensive reviewed and debated by several authoritative and regulatory bodies. The International Agency for Research on Cancer (IARC) classified GLY as ‘probably carcinogenic to humans’ (Group 2A) (EFSA, 2015) based on a specific association with non-Hodgkin’s lymphoma. In contrast to the IARC assessment, the European Food Safety Authority (EFSA) and the Environmental Protection Agency (EPA) concluded shortly after that there was insufficient scientific evidence to consider GLY as possibly carcinogen to humans (Portier et al., 2016; Benbrook, 2019). Despite all this controversy, European Union decided to extend the use of GLY until 2033 (Carrasco Cabrera et al., 2024).

While the debate on its carcinogenic potential continues, increasing evidence indicates GLY as an endocrine disruptor, a chemical with the ability to interfere with hormonal signaling pathways (Muñoz et al., 2021). For instance, GLY has been shown to modulate the activity of sex hormone receptors, particularly the estrogen receptor alpha, by enhancing its transcriptional activation in BC cell lines (Muñoz et al., 2021). Endocrine disruptors, such as pesticides, are of significant concern in BC due to their potential influence on the development and progression of the disease (Wan et al., 2022). In relation to that, a study by Franke et al. (2021) detected higher levels of AMPA in urine samples of women diagnosed with BC compared to controls, suggesting that AMPA exposure may be associated with an increased risk of BC. However, the association between cancer development and GLY exposure in epidemiological studies remains inconclusive. Therefore, further investigations in populations with high GLY exposure, such as those in Argentina, are of great importance.

In the present study, we conducted a case-control study to investigate the associations between urinary GLY levels and BC in women from a region of intense agricultural activity in Argentina, and also exploring residential proximity to agricultural fields as a potential risk factor for BC.

2 Methods

2.1 Population and data collection

This is a case-control study conducted from January to December 2021. The research protocol was approved by The Research Safety and Ethics Advisory Committee (Code CE2018-50) of the Facultad de Bioquímica y Ciencias Biológicas (Universidad Nacional del Litoral, Santa Fe, Argentina). In the current study, 90 women who attended the gynecology service of José María Cullen Provincial Hospital in Santa Fe were recruited to participate in the study. The study included participants aged 18–60 years old who had lived in various regions of the province of Santa Fe for at least 5 years prior to the interview.

The participants were either diagnosed with BC (cases, n = 30) or were cancer-free patients (controls, n = 60). Prior to their participation, the women were provided with comprehensive written information and were asked to provide their signed informed consent. The data of the participants were obtained from their medical records and survey responses. Since not all women completed the questionnaire, sample size for some variables was less than 30 in the case group and less than 60 in the control group. The questionnaire included items pertaining to demographic data, health information, lifestyle factors, reproductive history, family history of BC, and occupational risk factors. The participants’ addresses were obtained from medical records and georeferenced using Google Maps (Online Geocoder). The resulting coordinates were then classified based on the proximity to agricultural fields (as residential proximity). Addresses located within 4,000 m were considered ‘near’, while those located beyond 4,000 m were considered ‘far’ (Thompson et al., 2022). Out of 90 women surveyed, urine samples were collected from 70 patients (16 BC cases and 54 controls). For the remaining 20 patients, the conditions for urine sampling could not be met since women were menstruating or lacked the minimum retention time. The samples were transported to the Instituto de Salud y Ambiente del Litoral (UNL-CONICET) for processing on the day of collection. All samples were de-identified, stored at −80°C, and then shipped to the laboratory of the Programa de Investigación y Análisis de Residuos y Contaminantes Químicos (PRINARC) for analysis of GLY and AMPA concentrations.

2.2 Creatinine analysis

The urine samples were analyzed for creatinine by a colorimetric method using the Wiener Lab Creatinine Kit (code 1260001, Rosario, Argentina). The reaction of creatinine with alkaline picrate in a buffered medium results in the formation of a chromogen, which is then measured at 510 nm. The purpose of this analysis is to ascertain the validity of the samples and to normalize for urine concentration. According to the World Health Organization (World Health Organization, 1996), samples with creatinine levels below 30 mg/dL or above 300 mg/dL should be discarded. Samples with very low creatinine concentrations may interfere with the detection of low-level toxicants, while those with very high creatinine concentrations may indicate dehydration, which may affect the renal secretion, excretion, and/or reabsorption of target chemicals (Barr et al., 2005).

2.3 Urinary GLY and AMPA analyses

2.3.1 Standards and reagents

Crystalline standards of GLY (97%), AMPA (98%), GLY-9-fluorenylmethylchloroformate (GLY-FMOC) (92%) and AMPA-FMOC (98%) were from Dr. Ehrenstorfer (Augsburg, Germany), and the GLY Isotope-Labeled Internal Standard 1,2-13C2 15N (98%) was obtained from Toronto Research Chemicals (Toronto, Canada). Solutions of the derivatizing reagent FMOC-Cl (Sigma, St. Louis, MO, United States) and sodium tetraborate buffer (Anedra, San Fernando, Buenos Aires, Argentina) were separately prepared by dissolving the reagents in acetonitrile and water, respectively. Ultra-high performance liquid chromatography (UHPLC)-grade acetonitrile and methanol (OptimaTM, Fisher Scientific, NJ, United States) and deionized water produced with a Milli-Q System (Millipore, Bedford, MA, United States) were utilized for mobile phase preparation. Additionally, 5 mM ammonium acetate (Anedra, San Fernando, Buenos Aires, Argentina) was employed as modifier to promote ionization. Pesticide grade dichloromethane (Sintorgan, Buenos Aires, Argentina) was used for clean-up purposes.

2.3.2 Analytical method

Human urine samples were analysed for the presence of GLY and AMPA. Sample preparation was adapted from Bernal et al. (2010), with derivatisation based on the method outlined by Demonte et al. (2018), as described below. In an eppendorf tube, 500 μL of human urine was spiked with 20 μL of GLY 1,2- 13C2 15N, followed by the addition of 250 μL of acetonitrile. The mixture was vortexed for 1 min, subjected to 10 min of ultrasound treatment, and centrifuged at 15,000 rpm for 15 min at room temperature to precipitate proteins. Subsequently, 500 μL of the supernatant was transferred to another eppendorf tube, and the protein precipitation step was repeated once more. The supernatant (500 μL) underwent derivatisation by adding 84 μL of borate buffer (40 mM, pH 9) and 84 μL of FMOC-Cl. The concentration of FMOC-Cl was adjusted based on the creatinine content of the original sample. The reaction was allowed to proceed for 2 h at room temperature. Following derivatisation, the extracts were cleaned by liquid-liquid partition with 500 μL of dichloromethane. Finally, a fraction of the aqueous phase was injected into the UHPLC-MS/MS system.

Quality assurance and quality control (QAQC) procedures were implemented as follows. Due the unavailability of blank samples, samples tested without GLY were used as blank samples and each urine sample was spiked with internal standard. Duplicates of each sample were spiked with GLY-FMOC and AMPA-FMOC standard solution to achieve three different concentration levels: one-third of the samples were spiked at 0.50 μg/L, one-third at 1.00 μg/L, and the remaining at 5.00 μg/L.

Percentage recoveries ranged from 62% to 126%, with a relative standard deviation (RSD) lower than 25% for all samples and both compounds. The limit of detection (LOD) and the limit of quantitation (LOQ) were determinate using S/N ratios of 3 and 10, respectively, from 1.00 μg/L spiked samples chromatogramas. LOD were 0.10 μg/L for GLY and AMPA. The LOQs (0.50 μg/L for GLY and AMPA) were experimentally verified by analyzing spiked samples at LOQ level, with recoveries and RSD within acceptable ranges.

To assess linearity, isotopically labelled GLY and AMPA standard working solution was added to urine samples with varying creatinine contents: 31.30 mg/dL (low level), 153.05 mg/dL (medium level), and 300.00 mg/dL (high level). Since AMPA was not detected in samples, the matrix effect was assessed solely by analyzing six standard solutions of GLY 1,2-13C2 15N in triplicate within the range of 0.10–25.00 μg/L (0.10, 0.50, 1.00, 5.00, 10.00, 25.00 μg/L). As it was demonstrated that the matrix effect depends on the level of creatinine in the sample, the concentration of GLY was calculated with the internal standard added to each sample.

2.3.3 Chromatographic system and operating conditions

Liquid chromatography with tandem mass spectrometry (LC-MS/MS) analyses were performed using an Acquity UPLCTM liquid chromatograph (Waters, Milford, MA, United States) coupled to a triple quadrupole mass spectrometer equipped with an electrospray ionization source able to operate in positive and negative-ion mode (TQD, Waters Micromass, United Kingdom).

Chromatographic separation was evaluated on an X-Select CSH-C18 column (3.50 μm particle size, 100.00 × 4.60 mm i.d) at 40°C. Aliquots of 10 μL of standard and/or sample extracts were introduced by means of an auto-sampler (Waters, Milford, MA, United States).

The mobile phase consisted of water and acetonitrile (98:2) + 5 mM NH4Ac (solvent A) and acetonitrile (solvent B). Chromatographic and mass spectrometry data were handled using MassLynx software v 4.1 (Waters, Manchester, United Kingdom).

All these procedures were performed in the PRINARC at the Facultad de Ingeniería Química, Universidad Nacional del Litoral, Santa Fe, Argentina.

2.4 Data analysis

2.4.1 Statistical analysis

Statistical analysis was performed using R software (version 4.2.0). Urinary herbicide concentrations were reported as µg/g creatinine (creatinine-adjusted concentrations). A descriptive analysis of the surveyed variables was performed showing medians and interquartile ranges (IQR) for continuous variables, or frequencies and percentages for categorical variables, as appropriate. The continuous variables included age, weight, height, and body mass index (BMI) calculated from the latter two. Categorical variables were grouped into different dimensions: pesticide exposure (e.g., living near agricultural fields, childhood in rural areas, rural work, use of agrochemicals, family members working with pesticides, pesticide application at home), lifestyle factors (e.g., smoking, alcohol consumption, physical activity, consumption of vegetables, fruits, dairy products, processed meats, red meat, white meat, drinking water source), and socioeconomic status (e.g., educational level, employment status).

Due to the non-normal distribution of the data, non-parametric bivariate tests were used. Each BC patient was matched with two controls based on age, residence proximity classified based on the proximity to agricultural fields and sample provision. Associations between BC and each covariate were initially assessed using the Mann-Whitney U test or Kruskal-Wallis test for continuous variables, and Pearson’s χ2 or Fisher’s exact tests for categorical variables, to identify possible candidate variables for multivariate logistic regression. Variables with a p-value < 0.10 were included in the multivariable logistic regression model. In this model, the response variable was of the case/control type, whereby cases were assigned the value of ‘1’ and controls the value of ‘0’. Odds ratios were estimated with a 95% confidence interval (CI) to evaluate the strength of the association between the response variable and each predictor.

2.4.2 Supervised machine learning

Supervised learning is a machine learning (ML) paradigm in which the data set comprises labeled examples. Each data point contains features and an associated output label. The goal of supervised learning algorithms is to learn a function that maps feature numerical vectors (inputs) to labels (desired outputs) based on example input-output pairs (Russell and Norvig, 2016) also known as training examples (Bishop and Nasrabadi, 2006). A supervised learning algorithm analyzes the training data and produces an inferred function that can be used to map new testing examples, generalizing from the training data to unseen situations.

In this study, several ML classifiers were trained, including Multilayer Perceptron (MLP), Random Forest (RF), Gradient Boosting (Gboost), Bagging and K nearest neighbors (KNN) (Haykin, 1999; James et al., 2013). We have also included a logistic regression linear model as a benchmark for comparison with the machine learning models. The optimal model was selected using random cross-validation, with 1-fold completely random partition having 80% of the complete dataset allocated for training and 20% of the full dataset for testing and performance evaluation. For the MLP model, hyperparameters grid-search was performed with a small subset of the training data. The MLP architecture is hidden_layer_size = 100, activation function = reLU, automatic batch_size, learning_rate = 0.001, and maximum iteration number = 1,000. The optimizer employed was Adam. The other models were used with default parameters: RF, GBoost and Bagging with 100 estimators, KNN with k = 5.

The prediction quality of each model was evaluated using the F1 classical classification measure:

F1=2s+ps++p,s+TPTP+FN,pTPTP+FP

where s+ (recall) measures how good is a classification method at recognizing (and not missing) the true positives; the precision p measures the relation between true positives and false positives; F1, a harmonic score between precision and recall, is used to compare prediction methods. TP, FP and FN represent the number of true positive, false positive and false negative classifications, respectively.

Model explainability, named Explainable Artificial Intelligence (XAI) is crucial in any ML pipeline. In this work, XAI with SHAP (Shapley Additive Explanations) was employed to achieve model explainability. This work used the SHAP method (Lundberg and Lee, 2017), which is based on the game theoretically optimal Shapley values (Shapley, 1953), to explain individual predictions by computing the contribution of each feature to the prediction.

The SHAP explanation method utilizes coalitional game theory to compute Shapley values. In this method, the feature values of a data instance act as players in a coalition. Shapley values explain how to distribute the prediction playout fairly among the features. In the case of tabular data, a player can be an individual feature value. SHAP specifies the explanation as:

gx=0+j=1Mj

where g represents the explanation model, x’ represents an instance o data point, M represents the number of features and ϕj∈R represents the feature attribution for a feature j, also known as the values. The Shapley value for feature j indicates the value contributed by the j-th feature to the prediction output of this particular instance compared to the average prediction for the dataset.

In this study, we used the SHAP Summary Plot to show the importance of each feature for the RF model trained to predict cancer. In this plot, the effect of a feature on the classes is stacked to create the feature importance plot. The summary plot for multiclass classification can show what the model learns from the features. Each point on the summary plot represents a Shapley value for a feature and an instance. The y-axis position corresponds to the feature, and the x-axis position is determined by the Shapley value. The color represents the value of the feature from lowest (red) to highest (green) value. The features are ordered according to their importance. Additionaly, we obtained a detailed SHAP Summary plot for each class.

3 Results

3.1 Population characteristics

The questionnaire was completed by 30 BC cases and 60 controls. Figure 1 indicates the geographic locations of the women’s residences. General characteristics of case and control groups are shown in Table 1 and their dietary habits are displayed in Supplementary Table S1 (Supplementary Material). The control group had a median age of 42 years old (interquartile range (IQR): 36–50 years), while the BC group had a median age of 48 years old (IQR: 41–54 years). Body mass index (BMI) was not significantly different between the groups, with a median of 28.6 (IQR: 23.8–33.5) kg/m2 for controls and a median of 26.7 (IQR: 24.3–30.3) kg/m2 for cases (p = 0.409). Educational level was not significantly different between BC cases and controls (p = 0.467), with only 37.9% of BC cases and 40.7% of controls having completed secondary school. Menopausal status was significantly different between BC cases and controls (p = 0.028). The majority of participants, both cases (56.7%) and controls (61.7%) reported never having had a mammogram. Also, no significant differences were detected between the cases and controls with regard to the age at menarche, family history of BC, oral contraceptive use, hormone replacement therapy, breastfeeding, physical activity, smoking, alcohol consumption, and meat and vegetable intake. Table 2 shows the data on the participants’ exposure to pesticides. In this study, 66.7% of the BC cases and 23.3% of controls reported residing near an agricultural fiel (defined as residential proximity) (p < 0.001).

Figure 1
www.frontiersin.org

Figure 1. (A) Map of the province of Santa Fe showing the geographic locations of the women’s residences. Color represents both BC cases (pink) and the control group (blue). (B) The expanded area of the district of Santa Fe city. José María Cullen Provincial Hospital, where the patients attended for medical attention, is indicated with a red cross. (C) Plot sizes are proportional to the number of participants. The map figure was created using QGIS software version 3.22.15 (QGIS Development Teaikm). Vector layers were provided by the Instituto Geográfico Nacional (IGN) and the Gobierno de la Provincia de Santa Fe.

Table 1
www.frontiersin.org

Table 1. Characteristics of case and control women living in the city of Santa Fe and surroundings.

Table 2
www.frontiersin.org

Table 2. Data on pesticide exposure of case and control women in the city of Santa Fe and surroundings.

Urinary creatinine analysis showed 65 out of 70 samples were valid (14 BC cases and 51 controls) with creatinine values ranging 30–300 mg/dL. GLY was detected in 86.1% of the 65 urine samples tested. Concentrations ranged 0.37–10.07 µg GLY/g creatinine and no significant differences were observed between the BC cases and controls. AMPA was not detected in the samples tested.

The results of the multivariate logistic regression analysis showed a significant increase in the risk of BC for women residing near agricultural fields (OR: 7.38, 95%CI: 2.74-21.90) and for postmenopausal status (OR: 3.80, 95% CI: 1.32-11.70). None of the associations found were related to urinary GLY levels.

3.2 ML analysis

The dataset comprises 90 women, with 30 cases of BC and 60 controls. The target class, ‘Group’, indicates whether the participant has been diagnosed with BC (case group) or not (control group). Table 1 presents the description and details of each dataset feature used in this study. Binary variables (the ‘Yes/No’ type) were assigned values of 0 (No) and 1 (Yes). Additionally, variables with more than two levels were coded incrementally to indicate their respective levels.

In the cross-validation test, each model achieved the following F1 scores at the 20% data test partition: logistic regression 66.7%, MLP 61.1%, RF 77.8%, GBoost 72.2%, Bagging 61.1% and KNN 61.1%. Therefore, RF was selected as the best model for cancer prediction to perform the final training with the complete dataset.

After the cross-validation that was performed in order to select the best classifer, a RF classifier of 100 trees was trained with the complete dataset (no partitions) to classify instances into two classes: ‘no-cancer’ (Class NC) and ‘BC’ (Class BC cases), resulting in F1 = 68.6% for the complete dataset. The SHAP method was then employed to explain the results of the RF model. The features were ranked based on their importance, as shown in Figure 2. The SHAP bar graph shows the mean absolute SHAP values, which represents the average impact of each feature on the model’s predictions for each class. The most important feature for classification into these classes is living near agricultural fields. Figure 3 show the SHAP plots illustrating the impact of each feature on the model output for each separate class. It was observed that for class BC, the variable of ‘living near agricultural fields’ has a very high impact on the model output.

Figure 2
www.frontiersin.org

Figure 2. The SHAP summary bar graph shows the ranking of variable importance based on the mean absolute value (|SHAP value|) for both controls (Class NC) (blue) and BC cases (Class BC cases) (pink). FM: Family members.

Figure 3
www.frontiersin.org

Figure 3. The SHAP summary plot for each class: control cases (left) and BC cases (right) in detail. The y-axis position corresponds to the feature, and the x-axis position is determined by the Shapley value. The color indicates the value of the feature, ranging from lowest (red) to highest (green) value. SHAP values above 0 indicate a positive association with the outcome (for example, regarding BC cases, high values of the feature “living near agricultural fields” have a positive impact on the output of the model for this class). SHAP values below 0 indicate a negative association with the outcome (for example, regarding BC cases, low values of the feature “living near agricultural fields” have a negative impact on the output of the model for this class).

4 Discussion

The present study, similar to others conducted in different countries, demonstrated the ubiquitous presence of urinary GLY in women from Santa Fe, Argentina. The widespread presence of GLY in human samples reflects its pervasive occurrence in the environment, as GLY have been reported in various food matrices (Zoller et al., 2018; Fagan et al., 2020; Cruz and Murray, 2021) and, particularly in Argentina, in different environmental samples (Primost et al., 2017; Lutri et al., 2020; Okada et al., 2020). Moreover, a recent study examining the occurrence of pesticides in 64 small water bodies across regions with intense agricultural activity in Argentina and ten European countries revealed that GLY herbicide exhibited the highest median concentration (Navarro et al., 2024). Even more worrying, Argentina showed the highest overall concentration of pesticides in the water bodies analyzed (Navarro et al., 2024). This evidence suggests a significant burden of GLY exposure in the Argentine population.

In our study, the median GLY concentration, including both healthy subjects and those diagnosed with BC, was 1.23 μg/L, with a maximum of 3.50 µg/L observed in a healthy woman. These urinary GLY levels were higher compared to other studies carried out in Argentina in non-occupationally exposed men (median 0.191 ng/mL) (Filippi et al., 2024) and in a population from a rural village (range < 0,50–3,03 μg/L) with 19.2% of quantifiable samples (Bressán et al., 2021). Moreover, the GLY levels reported in the current work were also higher than those in healthy postmenopausal women in Southern California, United States (median 0.10 μg/L; maximum 3.01 μg/L) (Lucia et al., 2023), healthy lactating women in the United States (mean 0.28 μg/L; range 0.02–1.93 μg/L) (McGuire et al., 2016), and adults in Portugal (median 0.13 μg/L) (Nova et al., 2020). On the other hand, some authors detected similar maximum urinary GLY concentrations of 3.22 μg/L in Danish mothers (Knudsen et al., 2017) and 3.39 μg/L in Swedish young adults (Faniband et al., 2021). To better contextualize GLY and AMPA exposure levels across regions, Table 3 presents a comparative summary of GLY and AMPA levels reported in biological and environmental samples from Argentina, neighboring South American countries, and Europe.

Table 3
www.frontiersin.org

Table 3. GLY and AMPA concentrations detected in biological and environmental samples from South American and European countries.

In our study, the metabolite AMPA was not detected in any of the urine samples analyzed, which differs from other studies that have reported its presence. This discrepancy could be attributed to a number of factors including limited exposure, removal efficiency or limitations in the detection method (Aris and Leblanc, 2011). Our method achieved a LOD of 0.10 µg/L for AMPA, which is comparable to those reported in other biomonitoring study (Faniband et al., 2021). In that study, AMPA was detected in a higher percentage of urine samples than GLY; however, no significant correlation was observed between urinary concentrations of GLY and AMPA. This finding suggests that AMPA levels in urine may not necessarily reflect internal GLY metabolism but could instead be influenced by independent environmental exposure to AMPA. Futhermore, only about 1% of an ingested GLY dose is excreted unchanged in urine, with reported excretion rates for AMPA being even lower. These biological factors may explain the absence of AMPA in urine, particularly when exposure originates primarily from GLY.

Some studies have shown a correlation between pesticide exposure and increased BC risk in vulnerable populations (Franke et al., 2021; De Rezende et al., 2023; Panis et al., 2024). For instance, a case-control study in Paraná, Brazil, a region with extensive pesticide use, examined the impact of pesticide exposure on BC risk in rural women who performed cleaning tasks on pesticide-contaminated equipment and clothing (Panis et al., 2024). The authors found evidence of pesticide exposure, including GLY in urine samples, and revealed that women exposed to pesticides exhibited an elevated risk of BC and lymph node metastasis (Panis et al., 2024). In contrast, in our current study we did not find a direct association between the urinary GLY levels and BC risk. This finding suggests that the relationship between urinary GLY concentration and BC risk may be complex and influenced by additional factors beyond just GLY exposure. Furthermore, taking into account the ubiquitous occurrence of GLY in our healthy and unhealthy population, make it more difficult to establish an association between the levels of this herbicide and the risk of BC and other pathologies.

Importantly, the evaluation of residential proximity revelated that women residing near agricultural fields had increased risk of BC, a finding consistent across both clasiccal statistics and ML approach. This result aligns with previous studies reporting increased BC risk associated with proximity to agricultural areas and surrounding greenness in Spain (O’Callaghan-Gordo et al., 2018). Additionally, a case-control study conducted also in Spain revelead that children residing near agricultural areas had a higher risk of developing various types of cancer (Gómez-Barroso et al., 2016). In France, Carles et al. (2017) found an increased risk of meningioma among adults residing near open field crops. Despite the absence of an association between BC and urinary GLY levels in our study, the positive correlation detected between BC and proximity to cultivated fields suggests that exposure to agrochemicals in general, and not just GLY, may contribute to BC development. While our results are consistent with these previous studies, it is important to note that the existing literature is not yet sufficiently conslusive. This is due, in part, to the heterogeneity in study desings, with some focusing on different cancer types or employing various methods to assess pesticide exposure.

Another interesting finding of this work was the detection of an increased risk of BC in postmenopausal women. This result is consistent with previous data indicating that the majority of BC diagnoses occur in postmenopausal women (Kharb et al., 2020). The postmenopausal period is considered a particularly susceptible phase to the influence of endocrine-disrupting chemicals (EDCs), and growing evidence suggests that EDCs exposure may contribute to BC development, especially during this physiological period (Kass et al., 2020; Gouesse and Plante, 2022).

Notable, GLY has shown estrogenic properties and the ability to induce the proliferation of estrogen-dependent MCF-7 human BC cells (Thongprakaisang et al., 2013; Mesnage and Antoniou, 2017) and ERalpha (ESR1)-positive cholangiocarcinoma cells (Sritana et al., 2018). In vivo studies have also demonstrated GLY’s potential to induce hyperplastic changes in the mammary gland of aged female rats (Zanardi et al., 2020). Based on this evidence, GLY herbicide might be a contribuiting factor in the development of BC. However, further epidemiological research is necessary to address this point, with particular attention to menopausal status.

This study provides evidence for the widespread presence of GLY in urine samples from women residing in Santa Fe, Argentina. Despite the absence of a direct association between GLY levels and BC, we detected a significant association between BC risk and living near agricultural fields. Our findings highlight the importance of considering environmental factors, including pesticides, when assessing BC risk. Finally, this work underscores the need to develop strategies to reduce pesticide exposure and protect the health mainly of those populations living near agricultural fields.

4.1 Limitations and strengths of the study

Our study has limitations, mainly because of the small number of participants, which requires further research with a larger population to increase the power of the study and to better explore these relationships while accounting for relevant confounding factors such as, age, menopause status and family members working with pesticides. Another limitation is that our case-control ratio does not reflect the actual disease incidence in the general population due to recruitment constraints, which could potentially introduce selection bias in our analyses. The limited sample size prevented us from properly stratifying participants by menopausal status, which would be important in future studies given its established association with BC risk. Additionally, the utilization of a single urine sample may not be sufficient to accurately assess GLY exposure levels, given the short half-life of this compound. Nevertheles, our findings indicate that GLY is present in most of the samples analyzed, suggesting the widespread occurrence in the population. Furthermore, the results provide substancial insights into potential risk factors for BC in our population, particularly related with the place of residence. In addition, menopause may represent a vulnerable physiological period associated with increased BC risk in our population. We consider the data from the present study to be of significant importance for future decision-making in Argentina.

Data availability statement

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

Ethics statement

The studies involving humans were approved by The Research Safety and Ethics Advisory Committee of the Facultad de Bioquímica y Ciencias Biológicas, Universidad Nacional del Litoral, Santa Fe, Argentina. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

FD: Formal Analysis, Investigation, Methodology, Writing – original draft. VL: Writing – review and editing, Visualization. GS: Writing – review and editing, Formal Analysis, Software. TR: Writing – review and editing, Formal Analysis. SD’: Formal Analysis, Writing – review and editing. FP: Resources, Writing – review and editing. MR: Methodology, Writing – review and editing. LD: Methodology, Writing – review and editing. MM: Supervision, Writing – review and editing. JV: Conceptualization, Investigation, Supervision, Writing – review and editing, Funding acquisition.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research received funding through grants from the Argentine National Agency of Scientific and Technological Promotion (ANPCyT): PICT 2018 N° 1760, CONICET: PIP 2015 N° 11220150100338CO, and Universidad Nacional del Litoral: CAI+D 2020 PI TIPO II. All grants were awarded to Jorgelina Varayoud.

Acknowledgments

The authors express gratitude to the female participants for their collaboration in the study.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/ftox.2025.1579952/full#supplementary-material

Abbreviations

AMPA, aminomethylphosphonic acid; BC, Breast Cancer; BMI, body mass index; CI, confidence interval; EFSA, European Food Safety Authority; EPA, Environmental Protection Agency; GLY, glyphosate; IARC, International Agency for Research on Cancer; IQR, interquartile ranges; ML, machine learning; OR, Odds ratios; UHPLC-MS/MS, Ultra high-performance liquid chromatography–Mass spectrometry.

References

Acquavella, J. F., Alexander, B. H., Mandel, J. S., Gustin, C., Baker, B., Chapman, P., et al. (2004). Glyphosate biomonitoring for farmers and their families: results from the Farm Family Exposure Study. Environ. Health Perspect. 112, 321–326. doi:10.1289/ehp.6667

PubMed Abstract | CrossRef Full Text | Google Scholar

Alza-Camacho, W. R., García-Colmenares, J. M., and Chaparro-Acuña, S. P. (2016). Determinación voltamétrica de paraquat y glifosato en aguas superficiales. Cienc. Tecnol. Agropecu. 17, 331–345. doi:10.21930/rcta.vol17_num3_art:510

CrossRef Full Text | Google Scholar

Andreotti, G., Koutros, S., Hofmann, J. N., Sandler, D. P., Lubin, J. H., Lynch, C. F., et al. (2018). Glyphosate use and cancer incidence in the agricultural health study. JNCI J. Natl. Cancer Inst. 110, 509–516. doi:10.1093/jnci/djx233

PubMed Abstract | CrossRef Full Text | Google Scholar

Aparicio, V. C., De Gerónimo, E., Marino, D., Primost, J., Carriquiriborde, P., and Costa, J. L. (2013). Environmental fate of glyphosate and aminomethylphosphonic acid in surface waters and soil of agricultural basins. Chemosphere 93, 1866–1873. doi:10.1016/j.chemosphere.2013.06.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Aris, A., and Leblanc, S. (2011). Maternal and fetal exposure to pesticides associated to genetically modified foods in Eastern Townships of Quebec, Canada. Reprod. Toxicol. 31, 528–533. doi:10.1016/j.reprotox.2011.02.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashley-Martin, J., Huang, R., MacPherson, S., Brion, O., Owen, J., Gaudreau, E., et al. (2023). Urinary concentrations and determinants of glyphosate and glufosinate in pregnant Canadian participants in the MIREC study. Environ. Res. 217, 114842. doi:10.1016/j.envres.2022.114842

PubMed Abstract | CrossRef Full Text | Google Scholar

Ayoola, R. T., Olujimi, O. O., Bada, B. S., and Dedeke, G. A. (2023). Seasonal variations in the levels of glyphosate in soil, water and crops from three farm settlements in Oyo state, Nigeria. Heliyon 9, e20324. doi:10.1016/j.heliyon.2023.e20324

PubMed Abstract | CrossRef Full Text | Google Scholar

Barr, D. B., Wilder, L. C., Caudill, S. P., Gonzalez, A. J., Needham, L. L., and Pirkle, J. L. (2005). Urinary creatinine concentrations in the U.S. Population: implications for urinary biologic monitoring measurements. Environ. Health Perspect. 113, 192–200. doi:10.1289/ehp.7337

PubMed Abstract | CrossRef Full Text | Google Scholar

Benbrook, C. M. (2016). Trends in glyphosate herbicide use in the United States and globally. Environ. Sci. Eur. 28, 3. doi:10.1186/s12302-016-0070-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Benbrook, C. M. (2019). How did the US EPA and IARC reach diametrically opposed conclusions on the genotoxicity of glyphosate-based herbicides? Environ. Sci. Eur. 31, 2. doi:10.1186/s12302-018-0184-7

CrossRef Full Text | Google Scholar

Ben Khadda, Z., Bungau, S. G., El Balkhi, S., Ezrari, S., Radu, A.-F., Houssaini, T. S., et al. (2025). Urinary biomonitoring of exposure to glyphosate and its metabolite amino-methyl phosphonic acid among farmers and non-farmers in Morocco. Environ. Toxicol. Pharmacol. 113, 104620. doi:10.1016/j.etap.2024.104620

PubMed Abstract | CrossRef Full Text | Google Scholar

Bernal, J., Bernal, J. L., Martin, M. T., Nozal, M. J., Anadón, A., Martínez-Larrañaga, M. R., et al. (2010). Development and validation of a liquid chromatography–fluorescence–mass spectrometry method to measure glyphosate and aminomethylphosphonic acid in rat plasma. J. Chromatogr. B 878, 3290–3296. doi:10.1016/j.jchromb.2010.10.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Berni, I., Menouni, A., Creta, M., El Ghazi, I., Duca, R.-C., Godderis, L., et al. (2023). Exposure of children to glyphosate in Morocco: urinary levels and predictors of exposure. Environ. Res. 217, 114868. doi:10.1016/j.envres.2022.114868

PubMed Abstract | CrossRef Full Text | Google Scholar

Bishop, C. M., and Nasrabadi, N. M. (2006). Pattern recognition and machine learning. Springer.

Google Scholar

Blois, M. P., and Rendón, C. A. (2023). Un experimento que se hace público: Los roles del conocimiento científico y el Estado en las disputas en torno al glifosato en Argentina.

Google Scholar

Bonansea, R., Filippi, I., Wunderlin, D., Marino, D., and Amé, M. (2017). The fate of glyphosate and AMPA in a freshwater endorheic basin: an ecotoxicological risk assessment. Toxics 6, 3. doi:10.3390/toxics6010003

PubMed Abstract | CrossRef Full Text | Google Scholar

Bressán, I. G., Llesuy, S. F., Rodriguez, C., Ferloni, A., Dawidowski, A. R., Figar, S. B., et al. (2021). Optimization and validation of a liquid chromatography-tandem mass spectrometry method for the determination of glyphosate in human urine after pre-column derivatization with 9-fluorenylmethoxycarbonyl chloride. J. Chromatogr. B 1171, 122616. doi:10.1016/j.jchromb.2021.122616

PubMed Abstract | CrossRef Full Text | Google Scholar

Camiccia, M., Candiotto, L. Z. P., Gaboardi, S. C., Panis, C., and Kottiwitz, L. B. M. (2022). Determination of glyphosate in breast milk of lactating women in a rural area from Paraná state, Brazil. Braz. J. Med. Biol. Res. 55, e12194. doi:10.1590/1414-431x2022e12194

PubMed Abstract | CrossRef Full Text | Google Scholar

Carles, C., Bouvier, G., Esquirol, Y., Piel, C., Migault, L., Pouchieu, C., et al. (2017). Residential proximity to agricultural land and risk of brain tumor in the general population. Environ. Res. 159, 321–330. doi:10.1016/j.envres.2017.08.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Carrasco Cabrera, L., European Food Safety Authority (EFSA), , Di Piazza, G., Dujardin, B., Marchese, E., and Medina Pastor, P. (2024). The 2022 European Union report on pesticide residues in food. EFSA J. 22, e8753. doi:10.2903/j.efsa.2024.8753

PubMed Abstract | CrossRef Full Text | Google Scholar

Chang, V. C., Ospina, M., Xie, S., Andreotti, G., Parks, C. G., Liu, D., et al. (2024). Urinary biomonitoring of glyphosate exposure among male farmers and nonfarmers in the Biomarkers of Exposure and Effect in Agriculture (BEEA) study. Environ. Int. 187, 108644. doi:10.1016/j.envint.2024.108644

PubMed Abstract | CrossRef Full Text | Google Scholar

Chronister, B. N. C., Yang, K., Yang, A. R., Lin, T., Tu, X. M., Lopez-Paredes, D., et al. (2023). Urinary glyphosate, 2,4-D and DEET biomarkers in relation to neurobehavioral performance in ecuadorian adolescents in the ESPINA Cohort. Environ. Health Perspect. 131, 107007. doi:10.1289/EHP11383

PubMed Abstract | CrossRef Full Text | Google Scholar

Connolly, A., Koch, H. M., Bury, D., Koslitz, S., Kolossa-Gehring, M., Conrad, A., et al. (2022). A human biomonitoring study assessing glyphosate and aminomethylphosphonic acid (AMPA) exposures among farm and non-farm families. Toxics 10, 690. doi:10.3390/toxics10110690

PubMed Abstract | CrossRef Full Text | Google Scholar

Cruz, J. M., and Murray, J. A. (2021). Determination of glyphosate and AMPA in oat products for the selection of candidate reference materials. Food Chem. 342, 128213. doi:10.1016/j.foodchem.2020.128213

PubMed Abstract | CrossRef Full Text | Google Scholar

de Melo, K. G., Jacobucci, S. R. F., Garlipp, C. R., Trape, Â. Z., and Rosa, P. C. P. (2020). Determination of glyphosate in human urine from farmers in mato grosso-BR. Interam. J. Med. Health 3.

Google Scholar

Demonte, L. D., Michlig, N., Gaggiotti, M., Adam, C. G., Beldoménico, H. R., and Repetti, M. R. (2018). Determination of glyphosate, AMPA and glufosinate in dairy farm water from Argentina using a simplified UHPLC-MS/MS method. Sci. Total Environ. 645, 34–43. doi:10.1016/j.scitotenv.2018.06.340

PubMed Abstract | CrossRef Full Text | Google Scholar

De Rezende, L. M., Da Silva Santos, S., and Monteiro, G. T. R. (2023). Exposure to pesticides and breast cancer in the city of Petrópolis, Brazil. Environ. Sci. Pollut. Res. 30, 56534–56541. doi:10.1007/s11356-023-26420-8

PubMed Abstract | CrossRef Full Text | Google Scholar

De Roos, A. J., Blair, A., Rusiecki, J. A., Hoppin, J. A., Svec, M., Dosemeci, M., et al. (2005). Cancer incidence among glyphosate-exposed pesticide applicators in the Agricultural Health Study. Environ. Health Perspect. 113, 49–54. doi:10.1289/ehp.7340

PubMed Abstract | CrossRef Full Text | Google Scholar

De Roos, A. J., Fritschi, L., Ward, M. H., Monnereau, A., Hofmann, J., Bernstein, L., et al. (2022). Herbicide use in farming and other jobs in relation to non-Hodgkin’s lymphoma (NHL) risk. Occup. Environ. Med. 79, 795–806. doi:10.1136/oemed-2022-108371

PubMed Abstract | CrossRef Full Text | Google Scholar

European Food Safety Authority (EFSA) (2015). Scientific opinion on an alternative method for the hygienic treatment of bovine colostrum through a series of filtration steps. EFSA J. 13, 4139. doi:10.2903/j.efsa.2015.4139

PubMed Abstract | CrossRef Full Text | Google Scholar

Fagan, J., Bohlen, L., Patton, S., and Klein, K. (2020). Organic diet intervention significantly reduces urinary glyphosate levels in U.S. children and adults. Environ. Res. 189, 109898. doi:10.1016/j.envres.2020.109898

PubMed Abstract | CrossRef Full Text | Google Scholar

Faniband, M. H., Norén, E., Littorin, M., and Lindh, C. H. (2021). Human experimental exposure to glyphosate and biomonitoring of young Swedish adults. Int. J. Hyg. Environ. Health 231, 113657. doi:10.1016/j.ijheh.2020.113657

PubMed Abstract | CrossRef Full Text | Google Scholar

Filippi, I., Fernández, P., Grimalt, J. O., Butinof, M., Amé, M. V., and Muñoz, S. E. (2024). Glyphosate and AMPA in saliva and other traditional human matrices. New findings for less invasive biomonitoring to the exposure to pesticides. Environ. Adv. 15, 100474. doi:10.1016/j.envadv.2023.100474

CrossRef Full Text | Google Scholar

Franke, A. A., Li, X., Shvetsov, Y. B., and Lai, J. F. (2021). Pilot study on the urinary excretion of the glyphosate metabolite aminomethylphosphonic acid and breast cancer risk: the Multiethnic Cohort study. Environ. Pollut. 277, 116848. doi:10.1016/j.envpol.2021.116848

PubMed Abstract | CrossRef Full Text | Google Scholar

Frontera, P., Somma, A., Pérez Parada, A., and Manta, E. (2024). Monitoreo preliminar de glifosato y su principal metabolito en aguas superficiales de la cuenca del Río Santa Lucía y en agua potable. INNOTEC, e647. doi:10.26461/27.01

CrossRef Full Text | Google Scholar

Galli, F. S., Mollari, M., Tassinari, V., Alimonti, C., Ubaldi, A., Cuva, C., et al. (2024). Overview of human health effects related to glyphosate exposure. Front. Toxicol. 6, 1474792. doi:10.3389/ftox.2024.1474792

PubMed Abstract | CrossRef Full Text | Google Scholar

Gómez-Barroso, D., García-Pérez, J., López-Abente, G., Tamayo-Uria, I., Morales-Piga, A., Pardo Romaguera, E., et al. (2016). Agricultural crop exposure and risk of childhood cancer: new findings from a case–control study in Spain. Int. J. Health Geogr. 15, 18. doi:10.1186/s12942-016-0047-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Gouesse, R.-J., and Plante, I. (2022). Environmental exposure to brominated flame retardants: unraveling endocrine and mammary gland effects that may increase disease risk. Toxicol. Sci. 186, 190–207. doi:10.1093/toxsci/kfac006

PubMed Abstract | CrossRef Full Text | Google Scholar

Grau, D., Grau, N., Gascuel, Q., Paroissin, C., Stratonovitch, C., Lairon, D., et al. (2022). Quantifiable urine glyphosate levels detected in 99% of the French population, with higher values in men, in younger people, and in farmers. Environ. Sci. Pollut. Res. 29, 32882–32893. doi:10.1007/s11356-021-18110-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Haykin, S. S. (1999). Neural networks: a comprehensive foundation, 2. Upper Saddle River, NJ: Prentice Hall.

Google Scholar

James, G., Witten, D., Hastie, T., and Tibshirani, R. (2013). An introduction to statistical learning. Springer.

Google Scholar

Kass, L., Gomez, A. L., and Altamirano, G. A. (2020). Relationship between agrochemical compounds and mammary gland development and breast cancer. Mol. Cell. Endocrinol. 508, 110789. doi:10.1016/j.mce.2020.110789

PubMed Abstract | CrossRef Full Text | Google Scholar

Kharb, R., Haider, K., Neha, K., and Yar, M. S. (2020). Aromatase inhibitors: Role in postmenopausal breast cancer. Arch. Pharm. (Weinheim) 353, 2000081. doi:10.1002/ardp.202000081

PubMed Abstract | CrossRef Full Text | Google Scholar

Knudsen, L. E., Hansen, P. W., Mizrak, S., Hansen, H. K., Mørck, T. A., Nielsen, F., et al. (2017). Biomonitoring of Danish school children and mothers including biomarkers of PBDE and glyphosate. Rev. Environ. Health 32, 279–290. doi:10.1515/reveh-2016-0067

PubMed Abstract | CrossRef Full Text | Google Scholar

Kongtip, P., Nankongnab, N., Phupancharoensuk, R., Palarach, C., Sujirarat, D., Sangprasert, S., et al. (2017). Glyphosate and paraquat in maternal and fetal serums in Thai women. J. Agromedicine 22, 282–289. doi:10.1080/1059924X.2017.1319315

PubMed Abstract | CrossRef Full Text | Google Scholar

Lesseur, C., Pathak, K. V., Pirrotte, P., Martinez, M. N., Ferguson, K. K., Barrett, E. S., et al. (2022). Urinary glyphosate concentration in pregnant women in relation to length of gestation. Environ. Res. 203, 111811. doi:10.1016/j.envres.2021.111811

PubMed Abstract | CrossRef Full Text | Google Scholar

Lima, I. B., Boëchat, I. G., Fernandes, M. D., Monteiro, J. A. F., Rivaroli, L., and Gücker, B. (2022). Glyphosate pollution of surface runoff, stream water, and drinking water resources in Southeast Brazil. Environ. Sci. Pollut. Res. 30, 27030–27040. doi:10.1007/s11356-022-24167-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Lucia, R. M., Liao, X., Huang, W.-L., Forman, D., Kim, A., Ziogas, A., et al. (2023). Urinary glyphosate and AMPA levels in a cross-sectional study of postmenopausal women: associations with organic eating behavior and dietary intake. Int. J. Hyg. Environ. Health 252, 114211. doi:10.1016/j.ijheh.2023.114211

PubMed Abstract | CrossRef Full Text | Google Scholar

Łukasiewicz, S., Czeczelewski, M., Forma, A., Baj, J., Sitarz, R., and Stanisławek, A. (2021). Breast cancer—epidemiology, risk factors, classification, prognostic markers, and current treatment strategies—an updated review. Cancers 13, 4287. doi:10.3390/cancers13174287

PubMed Abstract | CrossRef Full Text | Google Scholar

Lundberg, S., and Lee, S. I. (2017). A unified approach to interpreting model predictions. United States: arXiv.

Google Scholar

Lutri, V. F., Matteoda, E., Blarasin, M., Aparicio, V., Giacobone, D., Maldonado, L., et al. (2020). Hydrogeological features affecting spatial distribution of glyphosate and AMPA in groundwater and surface water in an agroecosystem. Córdoba, Argentina. Sci. Total Environ. 711, 134557. doi:10.1016/j.scitotenv.2019.134557

PubMed Abstract | CrossRef Full Text | Google Scholar

Mac Loughlin, T. M., de Castro, M. C., López Aca, V., Orofino, A., Davidovich, I., Bernasconi, C., et al. (2017). La química ambiental como herramienta de transformación social: glifosato en espacios públicos urbanos y propuestas para su prohibición. Available online at: https://host170.sedici.unlp.edu.ar/server/api/core/bitstreams/e37f2069-9534-4322-acf8-43716f887d7b/content.

Google Scholar

McGuire, M. K., McGuire, M. A., Price, W. J., Shafii, B., Carrothers, J. M., Lackey, K. A., et al. (2016). Glyphosate and aminomethylphosphonic acid are not detectable in human milk. Am. J. Clin. Nutr. 103, 1285–1290. doi:10.3945/ajcn.115.126854

PubMed Abstract | CrossRef Full Text | Google Scholar

Mendez, M. J., Aimar, S. B., Aparicio, V. C., Ramirez Haberkon, N. B., Buschiazzo, D. E., De Gerónimo, E., et al. (2017). Glyphosate and Aminomethylphosphonic acid (AMPA) contents in the respirable dust emitted by an agricultural soil of the central semiarid region of Argentina. Aeolian Res. 29, 23–29. doi:10.1016/j.aeolia.2017.09.004

CrossRef Full Text | Google Scholar

Mesnage, R., and Antoniou, M. N. (2017). Facts and fallacies in the debate on glyphosate toxicity. Front. Public Health 5, 316. doi:10.3389/fpubh.2017.00316

PubMed Abstract | CrossRef Full Text | Google Scholar

Muñoz, J. P., Bleak, T. C., and Calaf, G. M. (2021). Glyphosate and the key characteristics of an endocrine disruptor: a review. Chemosphere 270, 128619. doi:10.1016/j.chemosphere.2020.128619

PubMed Abstract | CrossRef Full Text | Google Scholar

Navarro, I., De La Torre, A., Sanz, P., Abrantes, N., Campos, I., Alaoui, A., et al. (2024). Assessing pesticide residues occurrence and risks in water systems: a Pan-European and Argentina perspective. Water Res. 254, 121419. doi:10.1016/j.watres.2024.121419

PubMed Abstract | CrossRef Full Text | Google Scholar

Nova, P., Calheiros, C. S. C., and Silva, M. (2020). Glyphosate in Portuguese adults – a pilot study. Environ. Toxicol. Pharmacol. 80, 103462. doi:10.1016/j.etap.2020.103462

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Callaghan-Gordo, C., Kogevinas, M., Cirach, M., Castaño-Vinyals, G., Aragonés, N., Delfrade, J., et al. (2018). Residential proximity to green spaces and breast cancer risk: the multicase-control study in Spain (MCC-Spain). Int. J. Hyg. Environ. Health 221, 1097–1106. doi:10.1016/j.ijheh.2018.07.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Okada, E., Allinson, M., Barral, M. P., Clarke, B., and Allinson, G. (2020). Glyphosate and aminomethylphosphonic acid (AMPA) are commonly found in urban streams and wetlands of Melbourne, Australia. Water Res. 168, 115139. doi:10.1016/j.watres.2019.115139

PubMed Abstract | CrossRef Full Text | Google Scholar

Panis, C., Candiotto, L. Z. P., Gaboardi, S. C., Gurzenda, S., Cruz, J., Castro, M., et al. (2022). Widespread pesticide contamination of drinking water and impact on cancer risk in Brazil. Environ. Int. 165, 107321. doi:10.1016/j.envint.2022.107321

PubMed Abstract | CrossRef Full Text | Google Scholar

Panis, C., Candiotto, L. Z. P., Gaboardi, S. C., Teixeira, G. T., Alves, F. M., Da Silva, J. C., et al. (2024). Exposure to pesticides and breast cancer in an agricultural region in Brazil. Environ. Sci. Technol. 58, 10470–10481. doi:10.1021/acs.est.3c08695

PubMed Abstract | CrossRef Full Text | Google Scholar

PNCM (2024). Programa Nacional de Control de Cáncer de Mama (PNCM). Available online at: https://www.argentina.gob.ar/salud/instituto-nacional-del-cancer/institucional/pncm.

Google Scholar

Portier, C. J., Armstrong, B. K., Baguley, B. C., Baur, X., Belyaev, I., Bellé, R., et al. (2016). Differences in the carcinogenic evaluation of glyphosate between the international agency for research on cancer (IARC) and the European food safety authority (EFSA). J. Epidemiol. Community Health 70, 741–745. doi:10.1136/jech-2015-207005

PubMed Abstract | CrossRef Full Text | Google Scholar

Primost, J. E., Marino, D. J. G., Aparicio, V. C., Costa, J. L., and Carriquiriborde, P. (2017). Glyphosate and AMPA, “pseudo-persistent” pollutants under real-world agricultural management practices in the Mesopotamic Pampas agroecosystem, Argentina. Environ. Pollut. 229, 771–779. doi:10.1016/j.envpol.2017.06.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodrigues, N. R., and De Souza, A. P. F. (2018). Occurrence of glyphosate and AMPA residues in soy-based infant formula sold in Brazil. Food Addit. Contam. Part A 35, 723–730. doi:10.1080/19440049.2017.1419286

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodríguez-Bolaña, C., Pérez-Parada, A., Tesitore, G., Goyenola, G., Kröger, A., Pacheco, M., et al. (2023). Multicompartmental monitoring of legacy and currently used pesticides in a subtropical lake used as a drinking water source (Laguna del Cisne, Uruguay). Sci. Total Environ. 874, 162310. doi:10.1016/j.scitotenv.2023.162310

PubMed Abstract | CrossRef Full Text | Google Scholar

Russell, S. J., and Norvig, P. (2016). Artificial intelligence: a modern approach. Third edition, Global edition. Boston Columbus Indianapolis: Pearson.

Google Scholar

Schmidt, M., López, V. T., Tobías, M., Grinberg, E., and Merlinsky, G. (2022). Social and environmental conflicts caused by agrochemical use in Salta, Santiago del Estero and Santa Fe, Argentina. Ciênc. Saúde Coletiva 27, 1061–1072. doi:10.1590/1413-81232022273.04852021

PubMed Abstract | CrossRef Full Text | Google Scholar

Shapley, L. S. (1953). A value for n-person games. Contrib. Theory Games 2. doi:10.1515/9781400881970-018

CrossRef Full Text | Google Scholar

Singh, R., Shukla, A., Kaur, G., Girdhar, M., Malik, T., and Mohan, A. (2024). Systemic analysis of glyphosate impact on environment and human health. ACS Omega 9, 6165–6183. doi:10.1021/acsomega.3c08080

PubMed Abstract | CrossRef Full Text | Google Scholar

Sritana, N., Suriyo, T., Kanitwithayanun, J., Songvasin, B. H., Thiantanawat, A., and Satayavivad, J. (2018). Glyphosate induces growth of estrogen receptor alpha positive cholangiocarcinoma cells via non-genomic estrogen receptor/ERK1/2 signaling pathway. Food Chem. Toxicol. 118, 595–607. doi:10.1016/j.fct.2018.06.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, T., Boënne, W., Desmet, N., Seuntjens, P., Bronders, J., and Van Griensven, A. (2015). Quantification and characterization of glyphosate use and loss in a residential area. Sci. Total Environ. 517, 207–214. doi:10.1016/j.scitotenv.2015.02.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Thompson, S., Ritz, B., Cockburn, M., and Heck, J. E. (2022). Prenatal ambient pesticide exposure and childhood retinoblastoma. Int. J. Hyg. Environ. Health 245, 114025. doi:10.1016/j.ijheh.2022.114025

PubMed Abstract | CrossRef Full Text | Google Scholar

Thongprakaisang, S., Thiantanawat, A., Rangkadilok, N., Suriyo, T., and Satayavivad, J. (2013). Glyphosate induces human breast cancer cells growth via estrogen receptors. Food Chem. Toxicol. 59, 129–136. doi:10.1016/j.fct.2013.05.057

PubMed Abstract | CrossRef Full Text | Google Scholar

Vicini, J. L., Jensen, P. K., Young, B. M., and Swarthout, J. T. (2021). Residues of glyphosate in food and dietary exposure. Compr. Rev. Food Sci. Food Saf. 20, 5226–5257. doi:10.1111/1541-4337.12822

PubMed Abstract | CrossRef Full Text | Google Scholar

Wan, M. L. Y., Co, V. A., and El-Nezami, H. (2022). Endocrine disrupting chemicals and breast cancer: a systematic review of epidemiological studies. Crit. Rev. Food Sci. Nutr. 62, 6549–6576. doi:10.1080/10408398.2021.1903382

PubMed Abstract | CrossRef Full Text | Google Scholar

Weisenburger, D. D. (2021). A review and update with perspective of evidence that the herbicide glyphosate (roundup) is a cause of non-hodgkin lymphoma. Clin. Lymphoma Myeloma Leuk. 21, 621–630. doi:10.1016/j.clml.2021.04.009

PubMed Abstract | CrossRef Full Text | Google Scholar

World Health Organization (1996). “Biological monitoring of chemical exposure in the workplace: guidelines,” in International programme on chemical safety. Geneva, Switzerland: World Health Organization.

Google Scholar

World Health Organization (2024). Breast cancer fact sheet. Available online at: https://www.who.int/news-room/fact-sheets/detail/breast-cancer.

Google Scholar

Zanardi, M. V., Schimpf, M. G., Gastiazoro, M. P., Milesi, M. M., Muñoz-de-Toro, M., Varayoud, J., et al. (2020). Glyphosate-based herbicide induces hyperplastic ducts in the mammary gland of aging Wistar rats. Mol. Cell. Endocrinol. 501, 110658. doi:10.1016/j.mce.2019.110658

PubMed Abstract | CrossRef Full Text | Google Scholar

Zoller, O., Rhyn, P., Rupp, H., Zarn, J. A., and Geiser, C. (2018). Glyphosate residues in Swiss market foods: monitoring and risk evaluation. Food Addit. Contam. Part B 11, 83–91. doi:10.1080/19393210.2017.1419509

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: herbicide exposure, urine, breast cancer, case-control study, risk factors

Citation: Doná F, Lorenz V, Stegmayer G, Ricardo T, D’Iorio S, Ponzo F, Repetti MR, Demonte LD, Milesi MM and Varayoud J (2025) Residential proximity to agricultural fields, urinary glyphosate levels and breast cancer risk: a case-control study in Argentina. Front. Toxicol. 7:1579952. doi: 10.3389/ftox.2025.1579952

Received: 19 February 2025; Accepted: 12 May 2025;
Published: 21 May 2025.

Edited by:

Terje Svingen, Technical University of Denmark, Denmark

Reviewed by:

Hannah Lui Park, University of California, Irvine, United States
Leandro Goulart de Araujo, UMR5256 Institut de Recherches sur la Catalyse et l'Environnement de Lyon (IRCELYON), France

Copyright © 2025 Doná, Lorenz, Stegmayer, Ricardo, D’Iorio, Ponzo, Repetti, Demonte, Milesi and Varayoud. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jorgelina Varayoud, dmFyYXlvdWRAZmJjYi51bmwuZWR1LmFy

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.