Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 12 October 2023
Sec. Cancer Epidemiology and Prevention
This article is part of the Research Topic Childhood Leukemias in Latin America: Epidemiology, Causality, Novel Predictive Profiles and Therapeutic Strategies View all 22 articles

Childhood acute lymphoblastic leukemia survival and spatial analysis of socio-environmental risks in Mexico

Jaqueline Calderon-Hernandez,Jaqueline Calderon-Hernandez1,2Lizet Jarquin-YaezLizet Jarquin-Yañez3Luis Reyes-ArreguinLuis Reyes-Arreguin1Luis A. Diaz-PadillaLuis A. Diaz-Padilla4Jose Luis Gonzalez-CompeanJose Luis Gonzalez-Compean5Pablo Gonzalez-MontalvoPablo Gonzalez-Montalvo6Rebeca Rivera-Gomez,Rebeca Rivera-Gomez7,8Jairo R. Villanueva-ToledoJairo R. Villanueva-Toledo9Kristal PechKristal Pech10Oscar ArrietaOscar Arrieta11Yelda A. Leal,,*Yelda A. Leal11,12,13*
  • 1Facultad de Medicina/CIAAS, Universidad Autónoma de San Luis Potosi, San Luis Potosi, Mexico
  • 2Global Public Health Program, Boston College, Boston, MA, United States
  • 3Consejo Nacional de Humanidades, Ciencias y Tecnologías (CONAHCYT), Ciudad de México, Mexico
  • 4Servicio de Oncología Pediátrica de la Unidad Médica de Alta Especialidad (UMAE), Centro Médico Nacional “Ignacio García Téllez”, Instituto Mexicano de Seguro Social (IMSS), Mérida, Yucatán, Mexico
  • 5Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional Unidad Tamaulipas, Victoria, Tamaulipas, Mexico
  • 6Servicio de Oncología Pediátrica, Hospital O’Horán, Servicios de Salud de Yucatán/Facultad de Medicina, Universidad Autónoma de Yucatán, Yucatán, Mexico
  • 7Facultad de Ciencias de la Salud, Universidad Autónoma de Baja California, Tijuana, Baja California, Mexico
  • 8Comisión de Salud Fronteriza Mexico-Estados, Registro Poblacional de Cáncer de Tijuana BajaREG, Tijuana, Baja California, Mexico
  • 9Consejo Nacional de Humanidades, Ciencias y Tecnologías (CONAHCYT) – Fundación IMSS, A.C., Ciudad de México, Mexico
  • 10Departamento de Epidemiología del Instituto Nacional de Cancerología, Ciudad de México, Mexico
  • 11Coordinación del Registro Nacional de Cáncer del Instituto Nacional de Cancerología, Ciudad de México, Mexico
  • 12Registro Poblacional de Cáncer Mérida, Unidad Médica de Alta Especialidad (UMAE), Centro Médico Nacional “Ignacio García Téllez”, Instituto Mexicano de Seguro Social (IMSS), Mérida, Yucatán, Mexico
  • 13Centro Institucional de Capacitación y Registro de Cáncer, Coordinación de Investigación en Salud, Instituto Mexicano de Seguro Social (IMSS), Ciudad de México, Mexico

Background: Acute lymphoblastic leukemia (ALL) etiology remains largely unknown; incidence patterns by age, sex, and geographical distribution suggest a potential environmental role.

Aim: To identify ALL clusters from four contrasting urban areas of Mexico and to characterize the sources of environmental carcinogens.

Methods: Hospital-based ALL cases (n = 443) diagnosed in children <19 years old from the Metropolitan Zones of Merida and San Luis Potosi, the State of Mexico, and Tijuana were analyzed (2015–2020). ALL cases were coded according to the International Classification of Diseases for Oncology. ALL clusters were identified by Kernel Density, and excess risk was estimated. Data of particulate matter ≤2.5 µm (PM2.5) concentrations measured by community-monitoring stations were analyzed. Geocoded datasets of benzene, polycyclic aromatic hydrocarbons, and PM2.5 sources were analyzed to characterize patterns of exposure in ALL clusters.

Results: The survival rate for ALL ranged from 61.5% to 78.6%. Seven ALL clusters with excess risk (RR 1.4–2.3, p < 0.05) were identified. The carcinogen sources included artisanal brick kilns, gas stations, cement works, carpentry, paint, and chemical manufacturing establishments. PM2.5 levels ranged from 15 µg/m3 to 37 µg/m3 among study areas.

Conclusion: ALL clusters were identified at the community level; the excess risk could be explained by small-scale carcinogen sources. The levels of PM2.5 in outdoor air ranged from 3 to 6 times above the World Health Organization (WHO) air quality guidelines. Healthcare providers must raise awareness of the increased risk of ALL in children living near sources of environmental carcinogens; cancer control and prevention strategies must be steered from a multi-sectoral and multi-action perspective to protect children’s health.

Introduction

Childhood and adolescent cancer is a significant and relevant public health challenge; each year, an estimated 400,000 children <19 years old are diagnosed with cancer globally. Ninety percent of them live in low- and middle-income countries (LMICs) where treatment is often unavailable or unaffordable. As a result, less than 30% of children with cancer in LMICs survive, compared to more than 80% of those in high-income countries (HICs) (1, 2). The World Health Organization (WHO) Global Initiative on Childhood Cancer proposes reaching at least a 60% cancer survival rate, thereby saving one million lives of children with cancer by 2030 (3). The Global Initiative for Childhood Cancer is part of the response to the World Health Assembly resolution Cancer Prevention and Control through an Integrated Approach, focused on the reduction of premature mortality from non-communicable diseases (NCDs) and the achievement of universal health coverage (4).

Leukemia is the most frequently occurring cancer among children and adolescents; acute lymphoblastic leukemia (ALL) is the most common type, and over 80% of ALL cases are classified as B-lineage (B-ALL). ALL is characterized by the rapid proliferation of poorly differentiated progenitor cells inside the bone marrow (5, 6). Leukemia, among other health conditions, has been included in the group of “new pediatric conditions” (7). Some chemical toxicants in the environment are known and suspected causal factors of the disease (8, 9). It is estimated that 17% (7% to 42%) of all cancer disease burdens in children under five can be attributable to environmental causes; in today’s world, children are surrounded by an estimated 350,000 chemicals (10, 11).

Although the ALL etiology in children and adolescents remains largely unknown, evidence suggests exposure to ionizing radiation, albeit rare, and prenatal/postnatal exposures to pesticides are linked to an increased risk of ALL (1214). In addition, certain behaviors, such as pica (placing hands and objects in one’s mouth) and playing outside, increase children’s odds of exposure (4). Regarding exposures in urban environments, residential proximity to petrol stations and industrial sites and exposure to mixtures of benzene, particulate matter ≤2.5 µm (PM2.5), and polycyclic aromatic hydrocarbons (PAHs) from traffic emission are related to high risk of ALL (1518).

Interpretation of epidemiologic evidence to determine causality is complex and relies on a wide range of supporting data. The study of cancer and environmental exposures often involves collecting and studying large, complex datasets of cases, environmental exposure patterns or biological indexes, and additional variables (19). In the absence of desirable data on cancer cases from population-based cancer registries (PBCRs), hospital-based cancer registries (HBCRs) could be an alternative.

For this reason, the WHO, through the Global Initiative for Cancer Registry (GICR; www.gicr.iarc.fr) and the International Agency for Research on Cancer (IARC), has promoted establishing cancer registry and surveillance strategies, with standards and quality procedures to provide reliable and better data regarding cancer burden and setting priorities for cancer control and prevention, especially in LMICs.

Therefore, to increase our knowledge and awareness of environmental factors and ALL in Mexican children <19 years old, we designed this research to compare ALL survival rates in four contrasting urban areas (northern, central, and southern). Also, we present an overview analyzing the distribution of the disease from a territorial perspective aiming to identify ALL clusters, visualize areas vulnerable to the disease, and characterize fixed small-scale carcinogens hotbeds at the local/neighborhood level.

Materials and methods

Study areas

This study was conducted in four urban areas of Mexico: first, the Metropolitan Zone of Merida (MZ-Merida), located in the southeastern state of Yucatan Mexico (centroid: −89.62166, 20.96777). It includes the municipalities of Conkal, Merida, Kanasin, Ucu, and Uman. Merida City started a PBCR in 2016 to collect data on cancer cases adhering to international standards (20). The MZ-Merida has had accelerated urbanization and economic growth within the last few years; the size of the population is 995,129. Second, the Metropolitan Zone of San Luis Potosi (MZ-SLP) is located in the south-central state of San Luis Potosí (centroid: −100.94282, 22.27046). It includes the San Luis Potosi and Soledad de Graciano Sanchez municipalities, serving 911,908 inhabitants; the area has undergone accelerated industrialization and urban growth within the last decade (21). Third, the State of Mexico (centroid: −99.64537, 19.35596), located in the central part of the country, is a heavily industrialized and densely populated area (1,878,372 inhabitants). It is composed of 125 municipalities, which include Jiquipilco, Otzolotepec, Temoaya, Calimaya, Chapultepec, Lerma, Mexicaltzingo, San Mateo Atenco, Tianguistenco, Toluca, Coatepec Harinas, Ixtapan de la Sal, Tenancingo, and Villa Guerrero. Fourth, the municipality of Tijuana (Tijuana), located in the northwest state of Baja California (centroid: −117.00371, 32.5027), has a population of 922,523. Tijuana borders the southwest region of the United States and is characterized by a high density of gasoline and diesel engine vehicles (22). In 2018, the Tijuana-Cancer Registry was established to monitor and track cancer morbidity in the area (Figure 1).

FIGURE 1
www.frontiersin.org

Figure 1 Location of study areas on Mexico geopolitical map.

Source of patient data

Based on hospital records, ALL cases of children <19 years old from the MZ-Merida and Tijuana were included. ALL cases from the MZ-SLP and the State of Mexico were obtained from the registry of the Mexican Association of Aid to Children with Cancer (AMANC by its acronym in Spanish). These are a subset of hospital-based cases. AMANC is a non-profit organization (NPO) involved in the support of the population without medical insurance (formerly Seguro Popular). The HBCR works under international guidelines and standards developed by the IARC and International Association of Cancer Registries (IACR), including the International Classification of Diseases for Oncology (ICD-O3.2) (23, 24). The registry reports pediatric cancer cases through the International Classification of Childhood Cancer (ICCC) (2). Case-finding information was recorded on the following variables: sex, age, date of birth, place of residence, and basis of diagnosis. The HBCR uses CanReg5 software, which prevents the analysis of duplicates and/or non-existent codes. In addition, IARC-Tools and Link Plus 2.0 software were employed as internal quality control measures (25).

The datasets were validated with the IARC-CanReg5 methodology as previously mentioned. Cases with incomplete addresses were excluded to prevent geo-localization. This project was approved by the Health Research Council from the Mexican Institute of Social Security (IMSS by its acronym in Spanish) R-2016-785-065. Case anonymity was maintained by assigning an identifier to each case and eliminating patient names from the dataset.

Spatial exploratory analysis for ALL cluster identification

A descriptive spatial cluster analysis was calculated by Kernel Density Estimator using QGis©. Areas with the highest density of ALL cases were identified according to a distance radius of 2,310 m (meter) for the MZ-Merida, 1,755 m for the MZ-SLP, 12,259 m for the State of Mexico, and 2,215 m for Tijuana. A 200-m cell size was set for each area of study. After ALL spatial clusters were identified, clusters with the highest number of cases were selected, and the area was clipped to obtain a geographic polygon. To estimate the ALL excess risk, a Poisson distribution model was assumed, and the average population of each polygon was obtained from census data (2010 and 2020) belonging to the National Institute of Statistics and Geography (INEGI by its acronym in Spanish). The analysis was conducted with SatScan© software v 10.1.

Identification of emission sources and geo dataset integration

The latest IARC cancer monographs and National Toxicology Program Reports (US EPA) were consulted, and an exhaustive systematic literature review spanning from 2010 to 2022 was conducted to identify chemical substances associated with an increased risk of ALL in children <19 years old. Benzene, PAHs, and PM2.5 (26, 27) were selected to create an inventory of small, fixed sources. A geocoded dataset with the addresses of chemical and paint manufacturing establishments, artisanal brick kilns, gas stations, carpentry, and cement works was created for each area of study. Data were obtained from the National Statistical Directory of Economic Units (DENUE by its acronym in Spanish).

Scenario characterization by emission sources in the ALL clusters

The polluting emission sources were counted for each Basic Geostatistical Area (BGA) by the vector layer in the UTM Zone 14 N spatial reference system inside each polygon defined by the ALL clusters. Then, hierarchical cluster analysis was carried out to identify BGA clusters for each emission source of pollutants. We used the Euclidean and Manhattan distance matrices and the complete, single, and average linkage among others until the most appropriate number of clusters and the most compact and defined form were obtained. The selected method was the Ward.D2 with the Manhattan distance matrix, and all analyses were conducted with R studio© software (28). The proportion of BGA clusters for emission sources was estimated by each ALL cluster to identify the main emission sources.

Time trends of annual average measurements of PM2.5

PM2.5 hourly data for each area of study were obtained from various sources. The MZ-Merida data were acquired from the Mexican National Air Quality Reports (2011–2022). There is a single air quality monitoring station that began gathering data in 2013, but adequate and verifiable data (average/year) were only available from 2015 to 2016 (29). The MZ-SLP data were obtained from the Secretary of Ecology and Environmental Management of San Luis Potosí (30). The MZ-SLP has four monitoring stations; only two reported data for PM2.5 for 2016 (30). State of Mexico data were acquired from the National System of Air Quality Information (SINAICA, by its acronym in Spanish) (31), air monitoring started in 2011 (32), and eight stations reported PM2.5 (33). Finally, data belonging to Tijuana were also acquired from the Mexican National Air Quality Reports (2011–2022), and only two out of the four monitoring stations in the city reported PM2.5 levels (33). It is important to note that only the State of Mexico reported concentrations of PM2.5 with enough data (at least 75% of the hourly and daily averages) as required by the Official Mexican Standard (NOM) to analyze time trends (27). Only data that met the criteria of the NOM-025-SSA1-2014 were included in the analysis (29).

Statistical analysis

Descriptive statistics of ALL cases in CanReg5 were calculated, grouping by age and sex. Survival analysis was performed on R software (version 4.3.0). A non-parametric test was conducted by analyzing Kaplan–Meier curves for ALL cases in the MZ-Merida, the MZ-SLP, the State of Mexico, and Tijuana, and log-rank tests were used to determine the statistical significance of differences. p < 0.05 was considered to indicate statistical significance. Prevalence rates were calculated in each study area by considering average population data from the 2010 and 2020 censuses; age groups under study at the municipal and BGA level were also taken into account (prevalence rates at the BGA level was calculated for the selected ALL cluster polygons). Spatial analysis methods involving Kernel Density Estimation, discrete Poisson modeling, and hierarchical clustering were performed in Rstudio©, SatScan©, and QGis© software.

Results

ALL cases distributed by age, sex, and area of study

This study analyzed a total of 443 ALL cases diagnosed between 2015 and 2020 in children <19 years old. Fifty-nine percent were male (n = 261). Distribution by residence area was as follows: 95 (21.5%) cases from the MZ-Merida, 39 (8.8%) cases from the MZ-SLP, 192 (43.3%) from the State of Mexico, and 117 (26.4%) from Tijuana. Disease distribution among age groups was as follows: 91 cases (20.5%) between 0 and 4 years old, 139 cases (31.4%) between 5 and 9 years old, 119 cases (26.9%) between 10 and 14 years old, and 94 cases (21.2%) for an adolescent group between 15 and 19 years old; data are shown in Table 1.

TABLE 1
www.frontiersin.org

Table 1 Demographic characteristics of children <19 years old with acute lymphoblastic leukemia (2015–2020).

Survival analysis

In the MZ-Merida, the 5-year survival rate in children <19 years old diagnosed with ALL was 78.6% (95% CI, 70.1–88.2). Furthermore, survival analysis by sex found that male children in this area had a lower survival rate of 75.3% (95% CI, 64.5–87.9) when compared to female children with a rate of 85.4% (95% CI, 73.1–99.7). These results were not statistically significant (p = 0.34). The MZ-SLP documented the lowest overall survival rate of 61.5% (95% CI, 48.0–78.9). Further analysis by sex found that male children in this area also had the lowest survival rate at 56.5% (95% CI, 39.5–80.9) when compared to female children at 68.8% (95% CI, 49.4–85.7). These results were not statistically significant (p = 0.55). Cases from the State of Mexico had a 72.9% (95% CI, 65.4–81.2) survival rate. Male and female children in this area had a similar survival rate distribution at 72.3% (95% CI, 62.6–83.5) and 73.3% (95% CI, 62.2–86.2), respectively; but these results were not statistically significant (p = 0.76). Finally, Tijuana yielded a 5-year survival rate of 70.6% (95% CI, 58.8–84.8) in children <19 years old. Female children had a 66.7% (95% CI, 49.3–90.2) survival rate, while male children had 73.7% (95% CI, 59.2–91.7), although the difference was not statistically significant (p = 0.37); data are shown in Figure 2.

FIGURE 2
www.frontiersin.org

Figure 2 Kaplan–Meier curve for ALL survival in children.

Spatial ALL case distribution

The spatial analysis identified 24 ALL clusters distributed by the following: five in the MZ-Merida, six in the MZ-SLP, seven in the State of Mexico, and six in Tijuana. This preliminary spatial analysis selected 11 ALL clusters for carcinogens’ source characterization (identified with a red circle in Figure 3): three of them in the MZ-Merida identified as Cluster 1, located in the west-central area; Cluster 2, downtown area; and Cluster 3, located in the central-east area (Figure 3). Four cement works, 37 carpentry, 20 gas stations, and three chemical manufacturing establishments were identified in this area.

FIGURE 3
www.frontiersin.org

Figure 3 Density of acute lymphoblastic leukemia by study areas. Clusters with the highest number of cases were selected for characterization, identified by a red circle.

In the MZ-SLP, three clusters were identified: Cluster 1 in the north-central area, Cluster 2 in the downtown area, and Cluster 3 in the central-south area (Figure 3). The analysis geolocated 23 gas stations, 80 brick kilns, 16 paint manufacturing establishments, and six chemical manufacturing establishments in these clusters.

Three clusters were identified within the State of Mexico: Cluster 1 in the north-central region encompassing Jiquipilco, Otzolotepec, and Temoaya municipalities; Cluster 2, in the central region including the Calimaya, Chapultepec, Lerma, Mexicaltzingo, San Mateo Atenco, Tianguistenco, and Toluca municipalities; and Cluster 3 located in the south-central area surrounding the Coatepec Harinas, Ixtapan de la Sal, Tenancingo, and Villa Guerrero municipalities (Figure 3). Sources of carcinogens included two cement works, 411 carpentry, 161 gas stations, 230 brick kilns, 12 paint manufacturing establishments, and 50 chemical manufacturing establishments.

In Tijuana, two clusters were identified: Cluster 1 in the west and Cluster 2 in the northeast of the municipality (Figure 3). Seven carpentry, 52 gas stations, one brick kiln, three paint manufacturing establishments, and 12 chemical manufacturing establishments were pinpointed.

Prevalence of ALL and excess risk

The prevalence rate of ALL within the total identified eleven-clusters was 2.4 per 10,000 inhabitants (173 cases within a total population of 709,571). Excess risk (RR) was calculated in each study area by assuming a discrete Poisson distribution. The MZ-Merida had the highest RR at 2.3 (log-likelihood ratio 7.6; p-value 0.00032) in Clusters 1 and 3. The State of Mexico yielded an RR of 2.2 (log-likelihood ratio 5.5; p-value 0.0046), and Tijuana had an RR of 1.4 (log-likelihood ratio 4.1; p-value 0.03). The MZ-SLP yielded a non-significant excess risk (Table 2).

TABLE 2
www.frontiersin.org

Table 2 Prevalence and excess of risk for acute lymphoblastic leukemia in selected clusters.

Hierarchical clustering for identification of main emission sources by BGA

In the MZ-Merida, the analysis of hierarchical clustering showed a high proportion of BGA clusters of carpentry in the three ALL clusters (13.3%, 16.7%, and 24.1%), gas stations (13.3%, 16.7%, and 24.1%), and cement works (6.7%, 8.3%, and 6.9%) for Cluster 1, Cluster 2, and Cluster 3, respectively. In addition, 13.3% was for chemical manufacturing establishments in Cluster 1. In the MZ-SLP, a high proportion of chemicals manufacturing (22.2%) in Cluster 3, brick kilns (20%) in Cluster 1, and gas stations (19.4%) in Cluster 2 were identified. In addition, painting manufacturing establishments and gas stations were identified in the three clusters. In the State of Mexico, the BGA clusters of carpentry were identified in the three ALL clusters (36.8%, 28%, and 22.2% for Cluster 2, Cluster 1, and Cluster 3, respectively). The BGA cluster for gas stations was 38.9%, 24.4%, and 20% for Cluster 3, Cluster 2, and Cluster 1, respectively. In addition, in Cluster 2, brick kilns were pinpointed. In Tijuana, the high proportion of BGA clusters was gas stations (27.9% and 20.9%) (Figure 4).

FIGURE 4
www.frontiersin.org

Figure 4 Prevalence of Basic Geostatistical Area (BGA) clusters of fixed polluting sources within the area of the selected ALL (acute lymphoblastic leukemia) clusters.

Time trends of annual average measurements of PM2.5

The annual average concentrations of PM2.5 were analyzed in the four areas of study. The residents in these urban and metropolitan areas were exposed to outdoor air levels of annual average PM2.5 that were three to seven times higher than the WHO air quality annual average guideline. In 2019, Tijuana experienced the highest value of PM2.5 with 37 µg/m3, which was 7.4 times higher than the WHO guideline, while the MZ-Merida had the lowest value (15 µg/m3) (Figure 5).

FIGURE 5
www.frontiersin.org

Figure 5 Annual average concentration of PM2.5 in the four areas of study in Mexico (2011–2022). Dotted lines represent the World Health Organization (WHO) and the Official Mexican Standard (NOM) air quality guideline changes over time.

Discussion

The proportion of ALL cases in the four areas of study was higher in male than female; these data are consistent with the report of the global incidence of childhood and adolescent cancer (34). In regard to age, most ALL cases occurred in children aged 0–9 years within the MZ-Merida, the MZ-SLP, and Tijuana. In the State of Mexico, almost 70% of ALL afflicted children aged 5–14 years. This evidence further supports the hypothesis that ALL initiates in utero. Molecular proof of prenatal leukemogenesis is provided by backtracking of leukemic somatic mutations (gene translocations) in cord or blood samples at birth (35). Age, area of residence, and paternal or maternal exposures are responsible for differential patterns of exposure to environmental hazards. Exposure to carcinogens during preconception, prenatal, postnatal, and preadolescent periods is critical for ALL onset.

Survival analysis suggests the MZ-SLP had the lowest 5-year survival rate (61.5%) when compared to the national average of 61.8% as reported by a retrospective cohort study (2005–2015) (36). Although the 5-year survival rates for the State of Mexico, the MZ-Merida, and Tijuana (72.6%, 78.6%, and 70.6%, respectively) were higher than the national average, these values fall below those of HICs (37). Differences in survival rates may be attributed to inequities within healthcare system organizations, for example, differences in early diagnosing or early treatment, which are known to increase survival rates. Quality of medical facilities, health insurance, and state of residence are all factors known to further affect disease outcomes. For instance, uninsured children were almost twice to succumb to disease than insured children (34).

The ALL 5-year survival rate within the Mexican population has remained low throughout the years, with a study suggesting it being as low as 58% (38). This survival rate has remained low despite efforts to improve the healthcare system, suggesting additional external or environmental exposures may be at play (39). The analysis shows that children from the MZ-SLP have the lowest survival rate out of all the different areas analyzed, in addition to an extremely high number of known carcinogen sources. One of the clusters analyzed, for example, was composed of 80 artisanal brick kilns (20% of total BGA sources) known to be major mixture carcinogen sources (Figure 4). Further research has identified high levels of PAHs, benzene, and epigenetic changes in children living within these neighborhoods (40, 41). A recent publication reported that emissions from brick kilns are especially dangerous sources of pollutants since they employ materials such as oils, tires, wood, and industrial/municipal wastes as combustion materials (42).

ALL cases (ICD-O3.2) reported by HBCRs were analyzed to identify disease clusters containing patterns of small-scale carcinogen sources. Clusters of ALL were identified first, followed by an estimate of excess risk in each area of study. Spatial analysis detected 11 different clusters within the regions analyzed. However, only seven of these clusters had a statistically significant excess risk (RR ranged from 1.4 to 2.3). The highest estimated RR values were found in the MZ-Merida region (Table 2). Further analysis of fixed polluting sources in this study area found that 50% of them are related to benzene emissions (gas stations and carpentry), especially in Cluster 1. Chemical manufacturing was also prevalent in this region, accounting for almost 40% of carcinogen sources in Cluster 2 (Figure 4).

PM2.5 data were only available for the city of Merida from 2015 to 2016, with an annual average of 15 µg/m3. This value is three times greater than the WHO’s guideline for outdoor air. In addition to fixed sources of benzene and PAHs, practices such as burning solid wastes, incinerators, and crematories increase levels of carcinogens in the air. Studies have even reported particle-bound PAHs, possibly originating from burning fossil fuels and solid wastes (43).

Drinking contaminated water could be an additional route of carcinogen exposure for children living in the MZ-Merida region since filtration of pollutants by the soil characteristics and karstic aquifer facilitates and degrades drinking water quality in this part of the country (44).

In the State of Mexico, we identified gas stations and carpentry as the major carcinogen contributors (66% according to the type of fixed polluting source identified). The ALL Cluster 3 had an excess risk of 2.2 (p < 0.0046). This cluster comprises the south-central region of the state and includes the Coatepec Harinas, Ixtapan de la Sal, Tenancingo, and Villa Guerrero municipalities, with floriculture being the main economic activity within this area. One of the country’s main flower producers, 60% of the production process in this area is overseen by women, making them and their children vulnerable populations. The flower market has seen strong growth within the past two decades, creating more than 250,000 jobs. The Mexican floricultural market is extremely competitive internationally, with European countries being the main consumers (71% of the world’s production) (45, 46). A highly profitable market, its negative effects on children’s health are often underplayed. Although our research did not focus on pesticide exposure, we identified an area within this region yielding high ALL prevalence. Pesticide exposure is a well-documented, causal factor related to ALL (47). It was very interesting to see a high proportion of gas stations (38.9%) identified within this cluster, in an area where important economic activities demand the transportation of product. Gas stations are a catalyst for carcinogen exposure through gas and diesel combustion. This is a great example of an extremely complex scenario where carcinogen exposure may come from a variety of sources, providing spatial statistical methodologies a critical tool for cancer surveillance.

Disparities in death rates and cancer burdens associated with environmental pollutants are sometimes concealed by aggregated analysis of cancer data at the state, municipal, or medical facility level. In this study, we examined the ALL spatial distribution as a point process to identify highly vulnerable areas for ALL within urban areas. High spatial resolution information on pollution levels and disease rates is necessary to characterize inequities in air pollution, water or economic activities exposures, and related health risks. Hazards inside cities are not distributed evenly, contributing to persistent health disparities between neighborhoods and population sub-groups (48). Thus, spatiotemporal hazard distribution and scenario characterization analysis in small areas must be encouraged to protect children’s health.

Two main clusters were identified in Tijuana, with an ALL excess risk of 1.4 (p < 0.03). Hierarchical cluster analysis demonstrated a significant proportion of carcinogen sources in the form of gas stations (20% and 27.9%). Tijuana is strategically located within the United States–Mexico border and is home to 43 open entry ports. The San Ysidro Port of Entry (POE) is one of the busiest border crossings in the world, overseeing 70,000 vehicles on a daily basis. This area also sees large quantities of industrial imports and exports (e.g., “maquiladoras”), usually transported by large commercial vehicles across the Otay Mesa POE. These two main clusters were located in the northmost region of Tijuana, approximately 5–10 km from the two main POEs. It is estimated that the annual average of PM2.5 levels in this area is 7.4 times greater than the WHO air quality guidelines (Figure 5).

There is evidence of high benzene concentrations in Tijuana’s air (49). Benzene is the second-highest substance released into Mexico’s air vehicle emissions, with diesel engines being the major contributors (50). Additional sources of benzene include fumes of paint, gas stations, and cigarette smoke. The IARC classified benzene as a human carcinogen, highlighting its morbidity (27). ALL onset through benzene’s toxicological mechanisms in the lung and liver is well understood for cancer in these organs. After inhalation, benzene is metabolized to benzene oxide and hydroquinone. These products are further converted into benzoquinone in bone marrow, initiating cytotoxicity and DNA damage, thereby increasing the risk of developing ALL (51).

Several reports indicate that populations living near industrial or urban areas exhibit a higher proportion of cancer cases, with more than 90% of air pollution-related deaths occurring in LMICs (52). A recent study in Colombia identified leukemia cases clustering within a 1-km radius surrounding industrial sites (53). Another study in San Luis Potosi, Mexico, documented high excess risk in ALL cases clustered near brick kilns, downtown neighborhoods, surrounding industrial zones, heavy traffic zones, and areas with a high density of diesel-burning vehicles (Jarquin-Yañez et al, 2023 data in process of publication). In Spain, García-Perez et al. reported ALL excess risk in children living in urban areas (OR = 1.36; 95% CI 1.02–1.80) (54) and excess risk of ALL in children living within a 2.5-km radius from sources emitting pollutants (OR = 1.31; 95% CI 1.03–1.67) (54, 55). Sources of stationary and mobile polluting sources may include service stations, incinerators, industrial facilities, crematories, urban landfills, artisanal brick kilns, vehicles, and areas burning biomass as fuel (56). These sources emit benzene, PAHs, and PM2.5 into the atmosphere. Most recently, these substances have been associated with an elevated risk of pediatric ALL (16, 57, 58). This study was interested in further understanding the geographical distribution of ALL within areas of interest, with the goals of identifying the disease’s spatial patterns and characterizing small-scale carcinogen sources.

Due to the variability in spatial pollutant distribution among cities, cancer surveillance programs should conduct rapid and precise geospatial analysis to identify areas of greatest vulnerability. UNICEF’s 2022 report has urged local, state, and federal governments to prioritize the protection of pediatric environmental health, citing extensive evidence that simultaneous exposure to environmental hazards, even at low doses, can result in adverse health outcomes in exposed populations (59, 60).

The rapid economic and urban growth in urban areas, coupled with weak regulation policies, has increased the likelihood of carcinogen exposure among urban children. International science-based reports of carcinogens now recognize the importance of considering real-world exposure scenarios (61, 62). This approach is especially relevant for countries like Mexico, where 80% of the population resides in urban areas and is exposed to pollutant mixtures on a daily basis (63).

The lack of evidence from epidemiological studies is a critical issue surrounding environmental exposures affecting human health. Most of the evidence comes from studies in countries with solid regulations and bans on pollutants as a control measure. In Mexico, the magnitude of the burden of carcinogens and pollution exposures is poorly documented. Despite well-documented evidence from epidemiological and toxicological studies, pollution is often not recognized as a cause of cancer. Therefore, concrete actions are necessary to translate this evidence into cancer control measures for children and the general population. Children and adolescents who survive cancer are at greater risk of developing other health conditions at early ages (64). Moreover, children are more susceptible to exposure to environmental carcinogens because their organs and systems are developing, and their bodies are smaller compared to adults. In addition, children have specific behaviors such as putting their hands and objects in their mouths and playing outdoors, which increase their exposure to carcinogenic sources (4, 65).

The study has some limitations. The data analyzed were restricted to a period of 5 years (2015 to 2020) for ALL of MZ-SLP and the State of Mexico, and a sub-database of ALL was analyzed. The dataset of sources of exposure was consulted from DENUE, which is a non-exhaustive registry of pollution, and thus, there is a possibility of misclassification bias due to non-registered exposure.

Conclusion

The results suggest that children <19 years old living in neighborhoods with mixed environmental carcinogenic have an excess risk for ALL. Integrated spatial analysis is necessary to identify vulnerable populations in Mexico. Reducing exposure to carcinogens is essential for cancer prevention. However, it is a challenging task that requires sustained and coordinated actions from national and local governmental authorities. Expanding the network of PBCRs in Mexico is a priority, along with novel surveillance systems that incorporate environmental data and geospatial analysis to identify high-vulnerability areas. Multi-sectoral involvement of the healthcare sector, governmental agencies for environmental monitoring and control, urban planners, researchers, industry, and the population is necessary for cancer control and prevention strategies.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by Health Research Council from the Mexican Institute of Social Security (IMSS by its acronym in Spanish) R-2016-785-065. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

LY was the lead author of this paper; she conceptualized, wrote, edited, reviewed all sections, and approved the final version for submission. C-HJ conceptualized, wrote, edited, reviewed all sections, and approved the final version for submission. J-YL geocoded the databases and performed the spatial analysis. J-YL, R-AL, PK, V-TJ, and G-CJL reviewed and validated the database, produced the figures and tables, and approved the final version for submission. D-PLA, G-MP, R-GR, and AO edited and reviewed all sections. All authors contributed to the article and approved the submitted version.

Funding

The study was supported by the Mexican Institute of Social Security Foundation and the Health Research Council of the Mexican Institute of Social Security (FUNIMSS-FIS/IMSS/PROT/PRIO/18/078; PI: LY).

Acknowledgments

The authors are grateful to the cancer registry team and all the physicians and personnel of institutions for their collaboration with the data collection. The authors thank the Mexican Association of Aid to Children with Cancer (AMANC by its acronym in Spanish) for having provided the data. The authors greatly thank Miguel Cid-Rosas from Boston College for his critical review.

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.

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.

References

1. Huang J, Ssentongo P, Sharma R. Editorial: Cancer burden, prevention and treatment in developing countries. Front Public Health (2023) 10:1124473. doi: 10.3389/fpubh.2022.1124473

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Steliarova-Foucher E, Stiller C, Lacour B, Kaatsch P. International Classification of Childhood Cancer, third edition. Cancer (2005) 103:1457–67. doi: 10.1002/cncr.20910

PubMed Abstract | CrossRef Full Text | Google Scholar

3. World Health Organization. CureAll framework: WHO global initiative for childhood cancer: increasing access, advancing quality, saving lives (2021). Available at: https://apps.who.int/iris/bitstream/handle/10665/347370/9789240025271-eng.pdf?sequence=1&isAllowed=y (Accessed March 19, 2023).

Google Scholar

4. World Health Organization. Cancer prevention and control in the context of an integrated approach. World Health Assembly Resolution WHA70, 12 (2017). Available at: https://apps.who.int/gb/ebwha/pdf_files/WHA70/A70_R12-en.pdf (Accessed March 19, 2023).

Google Scholar

5. World Health Organization. Assessing national capacity for the prevention and control of noncommunicable diseases: report of the 2019 global survey (2020). Available at: https://apps.who.int/iris/handle/10665/331452 (Accessed June 4, 2023).

Google Scholar

6. Zhang J, Walsh MF, Wu G, Edmonson MN, Gruber TA, Easton J, et al. Germline mutations in predisposition genes in pediatric cancer. New Engl J Med (2015) 373:2336–46. doi: 10.1056/NEJMoa1508054

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Landrigan PJ. Environmental hazards for children in USA. Int J Occup Med Environ Health (1998) 11:189–94.

PubMed Abstract | Google Scholar

8. Landrigan PJ, Miodovnik A. Children’s health and the environment: an overview. Mount Sinai J Medicine: A J Trans Personalized Med (2011) 78:1–10. doi: 10.1002/msj.20236

CrossRef Full Text | Google Scholar

9. Landrigan PJ, Fuller R, Fisher S, Suk WA, Sly P, Chiles TC, et al. Pollution and children’s health. Sci Total Environ (2019) 650:2389–94. doi: 10.1016/j.scitotenv.2018.09.375

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Lam CG, Howard SC, Bouffet E, Pritchard-Jones K. Science and health for all children with cancer. Sci (1979) (2019) 363:1182–6. doi: 10.1126/science.aaw4892

CrossRef Full Text | Google Scholar

11. Binagwaho A, Laborde A, Landrigan PJ. Safeguarding children’s health in a changing global environment. Lancet (2022) 400:1176–8. doi: 10.1016/S0140-6736(22)01797-4

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Onyije FM, Olsson A, Baaken D, Erdmann F, Stanulla M, Wollschläger D, et al. Environmental risk factors for childhood acute lymphoblastic leukemia: an umbrella review. Cancers (Basel) (2022) 14:382. doi: 10.3390/cancers14020382

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Pech K, Pérez-Herrera N, Vértiz-Hernández ÁA, Lajous M, Farías P. Health risk assessment in children occupationally and para-occupationally exposed to benzene using a reverse-translation PBPK model. Int J Environ Res Public Health (2023) 20:2275. doi: 10.3390/ijerph20032275

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Van Maele-Fabry G, Gamet-Payrastre L, Lison D. Household exposure to pesticides and risk of leukemia in children and adolescents: Updated systematic review and meta-analysis. Int J Hyg Environ Health (2019) 222:49–67. doi: 10.1016/j.ijheh.2018.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Wei T, Jiao R, Nakyeyune R, Zang Z, Shao Y, Shen Y, et al. Exposure to outdoor air pollution at different periods and the risk of leukemia: a meta-analysis. Environ Sci pollut Res (2021) 28:35376–91. doi: 10.1007/s11356-021-14053-8

CrossRef Full Text | Google Scholar

16. Symanski E, Tee Lewis PG, Chen T-Y, Chan W, Lai D, Ma X. Air toxics and early childhood acute lymphocytic leukemia in Texas, a population based case control study. Environ Health (2016) 15:70. doi: 10.1186/s12940-016-0154-8

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Malavolti M, Malagoli C, Filippini T, Wise LA, Bellelli A, Palazzi G, et al. Residential proximity to petrol stations and risk of childhood leukemia. Eur J Epidemiol (2023) 38:771–82. doi: 10.1007/s10654-023-01009-0

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Ortega-García JA, López-Hernández FA, Cárceles-Álvarez A, Fuster-Soler JL, Sotomayor DI, Ramis R. Childhood cancer in small geographical areas and proximity to air-polluting industries. Environ Res (2017) 156:63–73. doi: 10.1016/j.envres.2017.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

19. CDC. Guidelines for examining unusual pattern of cancer and environmental concerns. (2022). Available at: https://www.cdc.gov/nceh/cancer-environment/pdfs/guidelines-for-examining-unusual-patterns-of-cancer-and-environmental-concerns-h.pdf (Accessed July 10, 2023).

Google Scholar

20. Leal YA, Reynoso-Noverón N, Aguilar-Castillejos LF, Meneses-García A, Mohar A, Piñeros M. Implementation of the population-based cancer registry in the city of Mérida, Mexico: process and early results. Salud Publica Mex (2019) 62:96. doi: 10.21149/10501

CrossRef Full Text | Google Scholar

21. Gobierno de Mexico, Data Mexico. Municipio de San Luis Potosí. (2020). Available at: https://dataMexico.org/en/profile/geo/san-luis-potosi (Accessed June 4, 2023).

Google Scholar

22. Instituto Nacional de Ecología y Cambio Climático (INECC), Dirección de Investigación sobre la Contaminación Urbana y Regional (DGICUR), Dirección de Investigación sobre la Calidad del Aire (DICA). Estudio de emisiones y actividad vehicular en Baja California, México. Reporte final (2011). Available at: https://www.gob.mx/cms/uploads/attachment/file/112408/2011_CGCSA_RSD_Baja_California.pdf (Accessed June 4, 2023).

Google Scholar

23. (ENCR) ENoCR. Development/recommendations on cancer registration standards: coding incidence date (2022). France: European Network of Cancer Registries (ENCR). Available at: https://www.encr.eu/sites/default/files/pdf/incideng.pdf (Accessed June 4, 2023).

Google Scholar

24. Fritz APC, Jack A, Shanmugaratnam K, Sobin L, Parkin D, Whelan S. Clasificación Internacional de Enfermedades para Oncología (CIE-O)-3. Washington, DC. USA: Organizacion Panamericana de la Salud/Organizacion Mundial de la Salud (2003).

Google Scholar

25. Ferlay J, Burkhard C, Whelan S, Parkin D. Check and conversion programs for cancer registries (IARC tools for cancer registries). Technical Report 42 (2005). International Agency for Research on Cancer. Available at: https://cri.tums.ac.ir/2/pbcr/References/Check%20and%20Conversion%20Programs%20for%20Cancer%20Registries.pdf (Accessed June 4, 2023).

Google Scholar

26. Lioy PJ, Weisel CP, Millette JR, Eisenreich S, Vallero D, Offenberg J, et al. Characterization of the dust/smoke aerosol that settled east of the World Trade Center (WTC) in lower Manhattan after the collapse of the WTC 11 September 2001. Environ Health Perspect (2002) 110:703–14. doi: 10.1289/ehp.02110703

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Intenational Agency for Research on Cancer (IARC). IARC monographs on the evaluation of carcinogenic risks to humans: benzene. (2018). International Agency for Research on Cancer. Available at: https://publications.iarc.fr/Book-And-Report-Series/Iarc-Monographs-On-The-Identification-Of-Carcinogenic-Hazards-To-Humans/Benzene-2018 (Accessed July 18, 2023).

Google Scholar

28. Jiménez Cuadrillero MÁ. RPubs - Clustering Jerárquico en R. Clustering Jerárquico en R. Available at: https://rpubs.com/mjimcua/clustering-jerarquico-en-r (Accessed June 4, 2023).

Google Scholar

29. NORMA Oficial Mexicana. NOM-025-SSA1-2014, Salud ambiental. Valores límite permisibles para la concentración de partículas suspendidas PM10 y PM2.5 en el aire ambiente y criterios para su evaluación (2014). Available at: https://www.dof.gob.mx/nota_detalle.php?codigo=5357042&fecha=20/08/2014#gsc.tab=0 (Accessed June 4, 2023).

Google Scholar

30. SEGAM. Red de monitoreo de la calidad del aire de la ciudad de San Luis Potosí y zona conurbada. Calidad del aire (2023). Available at: https://slp.gob.mx/segam/Paginas/Calidad%20del%20Aire/Calidad-del-Aire.aspx (Accessed [Accessed June 4, 2023).

Google Scholar

31. Instituto Nacional de Ecología y Cambio CLimático. Sistema Nacional de Información de la Calidad del Aire, SINAICA (2023). Available at: https://sinaica.inecc.gob.mx (Accessed June 4, 2023).

Google Scholar

32. INECC. Informe Nacional de la Calidad del Aire 2013. (2014). Available at: https://sinaica.inecc.gob.mx/archivo/informes/Informe2013.pdf (Accessed June 4, 2023).

Google Scholar

33. INECC. Informe Nacional de la Calidad del Aire 2020. (2021). México. Available at: https://sinaica.inecc.gob.mx/archivo/informes/Informe2020.pdf (Accessed June 4, 2023).

Google Scholar

34. INEGI. Estadísticas a propósito del día mundial contra el cáncer (4 de febrero). Datos nacionales. (2023). Available at: https://www.inegi.org.mx/contenidos/saladeprensa/aproposito/2023/EAP_Cancer.pdf (Accessed June 4, 2023).

Google Scholar

35. Marcotte EL, Spector LG, Mendes-de-Almeida DP, Nelson HH. The prenatal origin of childhood leukemia: potential applications for epidemiology and newborn screening. Front Pediatr (2021) 9:639479. doi: 10.3389/fped.2021.639479

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Muñoz-Aguirre P, Huerta-Gutierrez R, Zamora S, Mohar A, Vega-Vega L, Hernández-Ávila JE, et al. Acute lymphoblastic leukaemia survival in children covered by seguro popular in Mexico: A national comprehensive analysis 2005–2017. Health Syst Reform (2021) 7(1):e1914897. doi: 10.1080/23288604.2021.1914897

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Kandola A. What to know about leukemia survival rates. (2023). Available at: https://www.medicalnewstoday.com/articles/322990 (Accessed June 4, 2023).

Google Scholar

38. Castro-Ríos A, Reyes-Morales H, Pelcastre-Villafuerte BE, Rendón-Macías ME, Fajardo-Gutiérrez A. Socioeconomic inequalities in survival of children with acute lymphoblastic leukemia insured by social security in Mexico: a study of the 2007–2009 cohorts. Int J Equity Health (2019) 18:40. doi: 10.1186/s12939-019-0940-3

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Cabrera-Rivera LT, Sweetser B, Fuster-Soler JL, Ramis R, López-Hernández FA, Pérez-Martínez A, et al. Looking towards 2030: strengthening the environmental health in childhood–adolescent cancer survivor programs. Int J Environ Res Public Health (2022) 20:443. doi: 10.3390/ijerph20010443

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Díaz de León-Martínez L, Flores-Ramírez R, Rodriguez-Aguilar M, Berumen-Rodríguez A, Pérez-Vázquez FJ, Díaz-Barriga F. Analysis of urinary metabolites of polycyclic aromatic hydrocarbons in precarious workers of highly exposed occupational scenarios in Mexico. Environ Sci pollut Res (2021) 28:23087–98. doi: 10.1007/s11356-021-12413-y

CrossRef Full Text | Google Scholar

41. Alegría-Torres JA, Barretta F, Batres-Esquivel LE, Carrizales-Yáñez L, Pérez-Maldonado IN, Baccarelli A, et al. Epigenetic markers of exposure to polycyclic aromatic hydrocarbons in Mexican brickmakers: A pilot study. Chemosphere (2013) 91:475–80. doi: 10.1016/j.chemosphere.2012.11.077

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Barrera V, Contreras C, Mugica-Alvarez V, Galindo G, Flores R, MIranda J. PM2.5 Characterization and Source Apportionment Using Positive Matrix Factorization at San Luis Potosi City, Mexico, during the Years 2017–2018. Atmosphere (Basel) (2023) 14:1160. doi: 10.3390/atmos14071160

CrossRef Full Text | Google Scholar

43. Alvarez-Ospina H, Giordano S, Ladino LA, Raga GB, Muñoz-Salazar JI, Leyte-Lugo M, et al. Carabali G. Particle-bound polycyclic aromatic hydrocarbons (pPAHs) in Merida, Mexico. Aerosol Air Qual Res (2021) 21:200245. doi: 10.4209/aaqr.200245

CrossRef Full Text | Google Scholar

44. Rodriguez G AP, Araujo Leon JA. The glyphosate herbicide in Yucatan, Mexico. MOJ Bioequivalence Bioavailability (2018) 5(6):284–286. doi: 10.15406/mojbb.2018.05.00115

CrossRef Full Text | Google Scholar

45. Estrada Elvires AC, Contreras-Castillo JM, Valdivia-Alcalá R. Production and competitiveness of Mexican floriculture. Agro Productividad (2022) 15(6):11–17. doi: 10.32854/agrop.v14i6.2186

CrossRef Full Text | Google Scholar

46. Pereira PCG, Parente CET, Carvalho GO, Torres JPM, Meire RO, Dorneles PR, et al. A review on pesticides in flower production: A push to reduce human exposure and environmental contamination. Environ pollut (2021) 289:117817. doi: 10.1016/j.envpol.2021.117817

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Hernández A, Menéndez P. Linking pesticide exposure with pediatric leukemia: potential underlying mechanisms. Int J Mol Sci (2016) 17:461. doi: 10.3390/ijms17040461

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Castillo MD, Kinney PL, Southerland V, Arno CA, Crawford K, van Donkelaar A, et al. Estimating intra-urban inequities in PM 2.5 -attributable health impacts: A case study for Washington, DC. Geohealth (2021) 5(11):e22021GH000431. doi: 10.1029/2021GH000431

CrossRef Full Text | Google Scholar

49. Bretón JGC, Bretón RMC, Guarnaccia C, Severino R del CL, Marrón MR, Guzmán AR, et al. Characterization and sources of aromatic hydrocarbons (BTEX) in an urban site of Tijuana, Baja California, Mexico. AIP Conference Proceedings (2018), 1982(1):020020. doi: 10.1063/1.5045426

CrossRef Full Text | Google Scholar

50. Secretary of Environment and Natural Resources (SEMARNAT). Registro de Emisiones y Transferencia de Contaminantes (RETC) (2022). Available at: https://www.gob.mx/semarnat/acciones-y-programas/registro-de-emisiones-y-transferencia-de-contaminantes-retc (Accessed June 4, 2023).

Google Scholar

51. Khalade A, Jaakkola MS, Pukkala E, Jaakkola JJ. Exposure to benzene at work and the risk of leukemia: a systematic review and meta-analysis. Environ Health (2010) 9:31. doi: 10.1186/1476-069X-9-31

PubMed Abstract | CrossRef Full Text | Google Scholar

52. International Bank for Reconstruction and Development/The World Bank. The global health cost of ambient PM2.5 Air pollution. (2020). Available at: https://documents1.worldbank.org/curated/en/202401605153894060/pdf/World-The-Global-Cost-of-Ambient-PM2-5-Air-Pollution.pdf (Accessed June 4, 2023).

Google Scholar

53. Rodriguez-Villamizar LA, Moreno-Corzo FE, Valbuena-Garcia AM, Uribe Pérez CJ, Brome Bohórquez MR, García García HI, et al. Childhood leukemia in small geographical areas and proximity to industrial sources of air pollutants in three Colombian cities. Int J Environ Res Public Health (2020) 17:7925. doi: 10.3390/ijerph17217925

PubMed Abstract | CrossRef Full Text | Google Scholar

54. García-Pérez J, Gómez-Barroso D, Tamayo-Uria I, Ramis R. Methodological approaches to the study of cancer risk in the vicinity of pollution sources: the experience of a population-based case–control study of childhood cancer. Int J Health Geogr (2019) 18:12. doi: 10.1186/s12942-019-0176-x

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Ojeda Sánchez C, García-Pérez J, Gómez-Barroso D, Domínguez-Castillo A, Pardo Romaguera E, Cañete A, et al. Exploring blue spaces’ Effects on childhood leukaemia incidence: A population-based case–control study in Spain. Int J Environ Res Public Health (2022) 19:5232. doi: 10.3390/ijerph19095232

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Mudu P, Pérez Velasco R, Zastenskaya I, Jarosinska D. The importance and challenge of carcinogenic air pollutants for health risk and impact assessment. Eur J Public Health (2020) 30. doi: 10.1093/eurpub/ckaa165.841

CrossRef Full Text | Google Scholar

57. Heck JE, He D, Contreras ZA, Ritz B, Olsen J, Hansen J. Parental occupational exposure to benzene and the risk of childhood and adolescent acute lymphoblastic leukaemia: a population-based study. Occup Environ Med (2019) 76:527–9. doi: 10.1136/oemed-2019-105738

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Asenjo S, Nuñez O, Segú-Tell J, Pardo Romaguera E, Cañete Nieto A, Martín-Méndez I, et al. Cadmium (Cd) and Lead (Pb) topsoil levels and incidence of childhood leukemias. Environ Geochem Health (2022) 44:2341–54. doi: 10.1007/s10653-021-01030-w

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Boyd DR, Orellana M UN, Human Rights Council, Special Rapporteur on the Issue of Human Rights Obligations Relating to the Enjoyment of a Safe CH and SE, UN, Human Rights Council. Special Rapporteur on the Implications for Human Rights of the Environmentally Sound Management and Disposal of Hazardous Substances and Wastes. The right to a clean, healthy and sustainable environment : non-toxic environment : report of the Special Rapporteur on the Issue of Human Rights Obligations Relating to the Enjoyment of a Safe, Clean, Healthy and Sustainable Environment (2022). Available at: https://digitallibrary.un.org/record/3957797?ln=en#record-files-collapse-header (Accessed June 4, 2023).

Google Scholar

60. Kortenkamp A, Faust M. Regulate to reduce chemical mixture risk. Sci (1979) (2018) 361:224–6. doi: 10.1126/science.aat9219

CrossRef Full Text | Google Scholar

61. Madia F, Worth A, Whelan M, Corvi R. Carcinogenicity assessment: Addressing the challenges of cancer and chemicals in the environment. Environ Int (2019) 128:417–29. doi: 10.1016/j.envint.2019.04.067

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Lunn RM, Mehta SS, Jahnke GD, Wang A, Wolfe MS, Berridge BR. Cancer hazard evaluations for contemporary needs: highlights from new national toxicology program evaluations and methodological advancements. JNCI: J Natl Cancer Institute (2022) 114:1441–8. doi: 10.1093/jnci/djac164

CrossRef Full Text | Google Scholar

63. IARC. Air pollution and cancer. (2013). France: International Agency for research Cancer. Available at: https://www.iarc.who.int/wp-content/uploads/2018/07/AirPollutionandCancer161.pdf (Accessed June 4, 2023).

Google Scholar

64. Robison LL, Hudson MM. Survivors of childhood and adolescent cancer: life-long risks and responsibilities. Nat Rev Cancer (2014) 14:61–70. doi: 10.1038/nrc3634

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Landrigan PJ, Fuller R, Acosta NJR, Adeyi O, Arnold R, Basu N(, et al. The Lancet Commission on pollution and health. Lancet (2018) 391:462–512. doi: 10.1016/S0140-6736(17)32345-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: acute lymphoblastic leukemia (ALL), environmental exposure, childhood & adolescent, spatial analysis, cancer registry

Citation: Calderon-Hernandez J, Jarquin-Yañez L, Reyes-Arreguin L, Diaz-Padilla LA, Gonzalez-Compean JL, Gonzalez-Montalvo P, Rivera-Gomez R, Villanueva-Toledo JR, Pech K, Arrieta O and Leal YA (2023) Childhood acute lymphoblastic leukemia survival and spatial analysis of socio-environmental risks in Mexico. Front. Oncol. 13:1236942. doi: 10.3389/fonc.2023.1236942

Received: 08 June 2023; Accepted: 22 August 2023;
Published: 12 October 2023.

Edited by:

Juan Manuel Mejia-Arangure, Universidad Nacional Autonoma de Mexico, Mexico

Reviewed by:

Maria S. Pombo-de-Oliveira, National Cancer Institute (INCA), Brazil
Silvia Jiménez-Morales, National Institute of Genomic Medicine (INMEGEN), Mexico

Copyright © 2023 Calderon-Hernandez, Jarquin-Yañez, Reyes-Arreguin, Diaz-Padilla, Gonzalez-Compean, Gonzalez-Montalvo, Rivera-Gomez, Villanueva-Toledo, Pech, Arrieta and Leal. 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: Yelda A. Leal, eWVsZGFfbGVhbDAzQHlhaG9vLmNvbS5teDs=; eWVsZGEubGVhbEBpbXNzLmdvYi5teA==

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.