2.1. Study design and data collection
This is a descriptive, cross-sectional, and quantitative exploratory study with the objective of determining a clinicopathological signature associated with exposure to pesticides in rural women with breast cancer. Figure 1 illustrates the study design.
This study complied with the national and international regulatory standards for research involving humans and was approved by the Research Ethics Committee (CEP) of the State University of Western Paraná (UNIOESTE), under the number CAAE 35524814.4.0000.0107. All volunteers signed a free and informed consent form. Data, including unique materials, documentation, and codes used in the analysis, are available at https://github.com/Laboratorio-de-Analise-de-Dados/Article_pesticide_exposure.
A total of 923 women who attended the Francisco Beltrão Cancer Hospital (Ceonc) from May 2015 to April 2023 with images suggestive of breast lesions identified by mammograms and ultrasound were included. The study comprised the Eighth Paraná Health Region, which comprises 27 municipalities characterized predominantly by rural family work. To obtain the diagnosis of breast cancer, a biopsy of the suspicious lesion was performed by a pathologist, followed by anatomopathological analysis and immunohistochemistry. After excluding patients with benign lesions, 386 patients were included in the study as having a breast cancer diagnosis. Patients were categorized as occupationally exposed (n = 208) or unexposed to pesticides (n = 141).
The Eighth Regional Health Region of Paraná comprises about 500,000 inhabitants living in the following 27 municipalities: Ampére, Barracão, Bela Vista da Caroba, Boa Esperança do Iguaçu, Bom Jesus do Sul, Capanema, Cruzeiro do Iguaçu, Dois Vizinhos, Éneas Marques, Flor da Serra Do Sul, Francisco Beltrão, Manfrinópolis, Marmeleiro, Nova Esperança do Sudoeste, Nova Prata do Iguaçu, Pérola D'oeste, Pinhal de São Bento, Planalto; Pranchita, Realeza, Renascença, Salgado Filho, Salto do Lontra, Santa Izabel do Oeste, Santo Antônio do Sudoeste, São Jorge D'Oeste, and Verê. We chose this region to develop the study because Paraná is the state that sells the fourth highest amount of pesticides in Brazil.26 This reflects an extensive use of pesticides in the state’s agricultural activities, which play a significant role in the Gross Domestic Product (GDP) of the 27 municipalities that make up this health area.26 More than 50% of the region’s inhabitants engage in agricultural activities, with a particular focus on family farming. This population is subject to considerable pesticide exposure, especially glyphosate, atrazine, and 2,4-dichlorophenoxyacetic (2,4-D), which are widely used in soybean, corn, and wheat monocultures in the region.26 Figure 1A depicts the correlation between breast cancer cases and the amount of pesticides used by municipalities in the Eighth Health Region of Paraná.
Data were grouped into three categories: patient characteristics, tumor characteristics, and disease characteristics (Figure 1B).
2.1.1. Patient characteristics
The following data were used as defining parameters of patient characteristics: age in years at diagnosis and menopausal status at diagnosis, dichotomized as presence and absence. The patient's weight at diagnosis was also obtained in kilograms (kg), height in meters (m), and body mass index (BMI) in kg/m2. The profile of occupational pesticide exposure was obtained using a standardized data collection instrument validated for this purpose.27 The exposure criteria are based on continuous, unprotected, and direct handling of pesticides. Thus, included in the group of pesticide-exposed patients were rural women with a history of direct handling of pesticides without use of protective gloves during preparation and/or dilution of the poisonous solution, application of pesticides, and/or decontamination of personal protective equipment and/or washing of clothes used during spraying, who reported living at least 50% of their lives under direct handling of pesticides at least twice a week for every week of the year. The group of women unexposed to pesticides comprises urban workers with no previous or current history of occupational pesticide exposure.
2.1.2. Tumor characteristics
The following parameters were considered to contain the tumor characteristics: estrogen (ER) and progesterone (PR) receptors' expression profile in percentage (%), considering values greater than zero as positive, and zero as negative; amplification of the epidermal human growth factor receptor 2 (HER2) considering values of "3+" and "2+" with the positive FISH amplification test as positive and values zero, "1+" and "2+" without the FISH amplification test as negative; proliferation index Ki67 in %, considering values below/equal to 14% or above as cut-off (≤14% as low and <14% as high proliferation); molecular subtyping of breast tumors considering the classes Luminal A = any positivity for ER and/or PR and ki67 below/equal to 14%, Luminal B = any positivity for ER and/or PR and ki67 above 14%, HER2-amplified = any ER/PR/ki67 value and presence of amplification for HER2, and Triple Negative = ER/PR/HER2 negative and any ki67 value (as described by the St. Gallen Consensus)28; tumor size represented in mm, histological grade categorized as low (grades 1 and 2) and high (grade 3). Lymph node invasion, presence of angiolymphatic emboli, and occurrence of distant metastases were dichotomized as presence or absence.
2.1.3. Disease characteristics
The parameters considered disease characteristics were stratification of death risk and recurrence (stratified into low risk, intermediate risk, and high risk, as described in Joint Ordinance No. 5 of April 18, 2019),29 chemoresistance development, disease recurrence, and death. All were dichotomized as presence or absence.
2.2. Data analysis
All data wore processed in Python version 3.10.12. A descriptive statistical data analysis was performed, including disease, tumor, and patient characteristics. A ratio was calculated (dividing the frequencies of each molecular subtype) and the differences analyzed by Chi-square and Fisher's exact test for samples with n less than 5. A p>0.05 was considered significant. The parameters that showed significance were subjected to the contingency chi-square test to confirm their independence.