Evaluation of fluoride exposure using disability-adjusted life years and health risk assessment in south-western Iran: A novel Monte Carlo simulation

Consumption of fluoride-contaminated water is a worldwide concern, especially in developing countries, including Iran. However, there are restricted studies of non-single-value health risk assessment and the disease burden regarding fluoride intake nationwide. Prolonged exposure to excessive fluoride has been linked to adverse health effects such as dental and skeletal fluorosis. This can lead to under-mineralization of hard tissues, causing aesthetic concerns for teeth and changes in bone structure, increasing the risk of fractures. As such, we aimed to implement probability-based frameworks using Monte Carlo methods to explore the potential adverse effects of fluoride via the ingestion route. This platform consists of two sectors: 1) health risk assessment of various age categories coupled with a variance decomposition technique to measure the contributions of predictor variables in the outcome of the health risk model, and 2) implementing Monte Carlo methods in doseresponse curves to explore the fluoride-induced burden of diseases of dental fluorosis and skeletal fractures in terms of disability-adjusted life years (DALYs). For this purpose, total water samples of 8053 (N=8053) from 57 sites were analyzed in Fars and Bushehr Provinces. The mean fluoride concentrations were 0.75 mg/L and 1.09 mg/L, with maximum fluoride contents of 6.5 mg/L and 3.22 mg/L for the Fars and Bushehr provinces, respectively. The hazard quotient of the 95th percentile (HQ>1) revealed that all infants and children in the study area were potentially vulnerable to over-receiving fluoride. Sobol' sensitivity analysis indices, including first-order, second-order, and total order, disclosed that fluoride concentration (Cw), ingestion rate (IRw), and their mutual interactions were the most influential factors in the health risk model. DALYs rate of dental fluorosis was as high as 981.45 (uncertainty interval: UI 95 % 353.23-1618.40) in Lamerd, and maximum DALYs of skeletal fractures occurred in Mohr 71.61(49.75-92.71), in Fars Province, indicated severe dental fluorosis but mild hazard regarding fractures. Residents of the Tang-e Eram in Bushehr Province with a DALYs rate of 3609.40 (1296.68-5993.73) for dental fluorosis and a DALYs rate of 284.67 (199.11-367.99) for skeletal fractures were the most potentially endangered population. By evaluating the outputs of the DALYs model, the gap in scenarios of central tendency exposure and reasonable maximum exposure highlights the role of food source intake in overreceiving fluoride. This research insists on implementing defluoridation programs in fluoride-endemic zones to combat the undesirable effects of fluoride. The global measures presented in this research aim to address the root causes of contamination and help policymakers and authorities mitigate fluoride's harmful impacts on the environment and public health.