Dose Addition in the Induction of Craniofacial Malformations in Zebrafish Embryos Exposed to a Complex Mixture of Food-Relevant Chemicals with Dissimilar Modes of Action

A QSAR study for predicting malformation in zebrafish embryo

BACKGROUND

Developmental toxicity tests are extremely expensive, require a large number of animals, and are time-consuming. It is necessary to develop a new approach to simplify the analysis of developmental endpoints. One of these endpoints is malformation, and one group of ongoing methods for simplifying is in silico models. In this study, we aim to develop a quantitative structure-activity relationship (QSAR) model and identify the best algorithm for predicting malformations, as well as the most important and effective physicochemical properties associated with malformation.

METHODS

The dataset was extracted from a reliable database called COMPTOX. Physicochemical properties (descriptors) were calculated using Mordred and RDKit chemoinformatics software. The data were cleaned, preprocessed, and then split into training and testing sets. Machine learning algorithms, such as gradient boosting model (GBM) and logistic regression (LR), as well as deep learning models, including multilayer perceptron (MLP) and neural networks (NNs) trained with train set data and different sets of descriptors. The models were then validated with test set and various statistical parameters, such as Matthew's correlation coefficient (MCC) and balanced accuracy (BAC) score, were used to compare the models.

RESULTS

A set of descriptors containing with 78% AUC was identified as the best set of descriptors. Gradient boosting was determined to be the best algorithm with 78% predictive power.

CONCLUSIONS

The descriptors that were the most effective for developing models directly impact the mechanism of malformation, and GBM is the best model due to its MCC and BAC.

Dose addition in mixtures of compounds with dissimilar endocrine modes of action in in vitro receptor activation assays and the zebrafish sexual development test

BACKGROUND

Human exposure to pesticides is being associated with feminisation for which a decrease of the anogenital distance (AGD) is a sensitive endpoint. Dose addition for the cumulative risk assessment of pesticides in food is considered sufficiently conservative for combinations of compounds with both similar and dissimilar modes of action (MoA).

OBJECTIVE

The present study was designed to test the dose addition hypothesis in a binary mixture of endocrine active compounds with a dissimilar mode of action for the endpoint feminisation.

METHODS

Compounds were selected from a list of chemicals of which exposure is related to a decrease of the AGD in rats and completed with reference compounds. These chemicals were characterised using specific in vitro transcriptional activation (TA) assays for estrogenic and androgenic properties, leading to a final selection of dienestrol as an ER-agonist and flutamide, linuron, and deltamethrin as AR-antagonists. These compounds were then tested in an in vivo model, i.e. in zebrafish (Danio rerio), using sex ratio in the population as an endpoint in order to confirm their feminising effect and MoA. Ultimately, the fish model was used to test a binary mixture of flutamide and dienestrol.

RESULTS

Statistical analysis of the binary mixture of flutamide and dienestrol in the fish sexual development tests (FSDT) with zebrafish supported dose addition.

Basic concepts of mixture toxicity and relevance for risk evaluation and regulation

Exposure to multiple substances is a challenge for risk evaluation. Currently, there is an ongoing debate if generic "mixture assessment/allocation factors" (MAF) should be introduced to increase public health protection. Here, we explore concepts of mixture toxicity and the potential influence of mixture regulation concepts for human health protection. Based on this analysis, we provide recommendations for research and risk assessment. One of the concepts of mixture toxicity is additivity. Substances may act additively by affecting the same molecular mechanism within a common target cell, for example, dioxin-like substances. In a second concept, an "enhancer substance" may act by increasing the target site concentration and aggravating the adverse effect of a "driver substance". For both concepts, adequate risk management of individual substances can reliably prevent adverse effects to humans. Furthermore, we discuss the hypothesis that the large number of substances to which humans are exposed at very low and individually safe doses may interact to cause adverse effects. This commentary identifies knowledge gaps, such as the lack of a comprehensive overview of substances regulated under different silos, including food, environmentally and occupationally relevant substances, the absence of reliable human exposure data and the missing accessibility of ratios of current human exposure to threshold values, which are considered safe for individual substances. Moreover, a comprehensive overview of the molecular mechanisms and most susceptible target cells is required. We conclude that, currently, there is no scientific evidence supporting the need for a generic MAF. Rather, we recommend taking more specific measures, which focus on compounds with relatively small ratios between human exposure and doses, at which adverse effects can be expected.

Incorporating Primer Amplification Efficiencies in Quantitative Reverse Transcription Polymerase Chain Reaction Experiments; Considerations for Differential Gene Expression Analyses in Zebrafish

Quantitative reverse transcription polymerase chain reaction (RT-qPCR) is commonly used to measure the mRNA expression of target genes in zebrafish. Gene expression values from RT-qPCR are typically reported as relative fold-changes, and relative quantification of RT-qPCR data incorporates primer amplification efficiency values for each target gene. We describe the influence of the primer amplification efficiency analysis method on RT-qPCR gene expression fold-change calculations. This report describes (1) a sample analysis demonstrating incorporation of primer amplification efficiency into RT-qPCR analysis for comparing gene expression of a gene of interest between two groups when normalized to multiple reference genes, (2) the influence of differences in primer amplification efficiencies between measured genes on gene expression differences calculated from theoretical delta-Cq (dCq) values, and (3) an empirical comparison of the influence of three methods of defining primer amplification efficiency in gene expression analyses (delta-delta-Cq [ddCq], standard curve, LinRegPCR) using mRNA measurements of a set of genes in zebrafish embryonic development. Given the need to account for the influence of primer amplification efficiency along with the simplicity of using software programs (LinRegPCR) to measure primer amplification efficiency from RT-qPCR data, we encourage using empirical measurements of primer amplification efficiency for RT-qPCR analysis of differential gene expression in zebrafish.

Dose additive maternal and offspring effects of oral maternal exposure to a mixture of three PFAS (HFPO-DA, NBP2, PFOS) during pregnancy in the Sprague-Dawley rat

Simultaneous exposure to multiple per- and polyfluoroalkyl substances (PFAS) is common in humans across the globe. Individual PFAS are associated with adverse health effects, yet the nature of mixture effects after exposure to two or more PFAS remains unclear. Previously we reported that oral administration of hexafluoropropylene oxide-dimer acid (HFPO-DA, or GenX), Nafion byproduct 2 (NBP2), or perfluorooctane sulfonate (PFOS) individually during pregnancy produced maternal and F1 effects. Here, we hypothesized that responses to the combined exposure to these three PFAS would be dose additive. Pregnant Sprague-Dawley rats were exposed to a fixed-ratio equipotent mixture where the top dose contained each PFAS at their ED50 for neonatal mortality (100 % dose = PFOS 3 mg/kg; NBP2 10 mg/kg; HFPO-DA 110 mg/kg), followed by a dilution series (33.3, 10, 3.3, and 1 %) and vehicle controls (0 % dose). Consistent with the single chemical studies, dams were exposed from gestation day (GD)14-18 or from GD8-postnatal day (PND2). Fetal and maternal livers on GD18 displayed multiple significantly upregulated genes associated with lipid and carbohydrate metabolism at all dose levels, while dams displayed significantly increased liver weight (≥3.3 % dose) and reduced serum thyroid hormones (≥33.3 % dose). Maternal exposure from GD8-PND2 significantly reduced pup bodyweights at birth (≥33.3 % dose) and PND2 (all doses), increased neonatal liver weights (≥3.3 % dose), increased pup mortality (≥3.3 % dose), and reduced maternal bodyweights and weight gain at the top dose. Echocardiography of adult F1 males and females identified significantly increased left ventricular anterior wall thickness (~10 % increase), whereas other cardiac morphological, functional, and transcriptomic measures were unaffected. Mixture effects in maternal and neonatal animals conformed to dose addition using a relative potency factor (RPF) analysis. Results support dose addition-based cumulative assessment approaches for estimating combined effects of PFAS co-exposure.

Additive cardiotoxicity of a bisphenol mixture in zebrafish embryos: The involvement of calcium channel and pump

Bisphenol A (BPA) and its analogs, such as bisphenol F (BPF), bisphenol AF (BPAF), and bisphenol B (BPB), are often simultaneously detected in environmental and human specimens. Thus, assessing the toxicity of bisphenol (BP) mixtures is more relevant than assessing that of each BP type. Here, we found that BPs, individually or in a mixture, concentration-dependently and additively increased the mortality of zebrafish embryos (ZFEs) at 96 h post fertilization (hpf) and induced bradycardia (i.e., reduced heart rate) at 48 hpf, indicating their cardiotoxic potency. BPAF was the most potent, followed by BPB, BPA, and BPF. We then explored the mechanism underlying BP-induced bradycardia in ZFEs. Although BPs increased the mRNA expression of the estrogen-responsive gene, treatment with the estrogen receptor inhibitor ICI 182780 did not prevent BP-induced bradycardia. Because they did not change cardiomyocyte counts or heart development-related gene expression, BPs might not affect cardiomyocyte development. By contrast, BPs might impair calcium homeostasis during cardiac contraction and relaxation through the downregulation of the expression of the mRNAs for the pore-forming subunit of L-type Ca2+ channel (LTCC; cacna1c) and sarco/endoplasmic reticulum Ca2+-ATPase (SERCA; atp2a2a). BPs reduced SERCA activity significantly. BPs also potentiated the cardiotoxicity induced by the LTCC blocker nisoldipine, conceivably by inhibiting SERCA activity. In conclusion, BPs additively induced bradycardia in ZFEs, possibly by impeding calcium homeostasis during cardiac contraction and relaxation. BPs also potentiated the cardiotoxicity of calcium channel blockers.

New approach methodologies to facilitate and improve the hazard assessment of non-genotoxic carcinogens-a PARC project

Carcinogenic chemicals, or their metabolites, can be classified as genotoxic or non-genotoxic carcinogens (NGTxCs). Genotoxic compounds induce DNA damage, which can be detected by an established in vitro and in vivo battery of genotoxicity assays. For NGTxCs, DNA is not the primary target, and the possible modes of action (MoA) of NGTxCs are much more diverse than those of genotoxic compounds, and there is no specific in vitro assay for detecting NGTxCs. Therefore, the evaluation of the carcinogenic potential is still dependent on long-term studies in rodents. This 2-year bioassay, mainly applied for testing agrochemicals and pharmaceuticals, is time-consuming, costly and requires very high numbers of animals. More importantly, its relevance for human risk assessment is questionable due to the limited predictivity for human cancer risk, especially with regard to NGTxCs. Thus, there is an urgent need for a transition to new approach methodologies (NAMs), integrating human-relevant in vitro assays and in silico tools that better exploit the current knowledge of the multiple processes involved in carcinogenesis into a modern safety assessment toolbox. Here, we describe an integrative project that aims to use a variety of novel approaches to detect the carcinogenic potential of NGTxCs based on different mechanisms and pathways involved in carcinogenesis. The aim of this project is to contribute suitable assays for the safety assessment toolbox for an efficient and improved, internationally recognized hazard assessment of NGTxCs, and ultimately to contribute to reliable mechanism-based next-generation risk assessment for chemical carcinogens.

Determination of in vitro hepatotoxic potencies of a series of perfluoroalkyl substances (PFASs) based on gene expression changes in HepaRG liver cells

Per- and polyfluoroalkyl substances (PFASs) are omnipresent and have been shown to induce a wide range of adverse health effects, including hepatotoxicity, developmental toxicity, and immunotoxicity. The aim of the present work was to assess whether human HepaRG liver cells can be used to obtain insight into differences in hepatotoxic potencies of a series of PFASs. Therefore, the effects of 18 PFASs on cellular triglyceride accumulation (AdipoRed assay) and gene expression (DNA microarray for PFOS and RT-qPCR for all 18 PFASs) were studied in HepaRG cells. BMDExpress analysis of the PFOS microarray data indicated that various cellular processes were affected at the gene expression level. From these data, ten genes were selected to assess the concentration-effect relationship of all 18 PFASs using RT-qPCR analysis. The AdipoRed data and the RT-qPCR data were used for the derivation of in vitro relative potencies using PROAST analysis. In vitro relative potency factors (RPFs) could be obtained for 8 PFASs (including index chemical PFOA) based on the AdipoRed data, whereas for the selected genes, in vitro RPFs could be obtained for 11-18 PFASs (including index chemical PFOA). For the readout OAT5 expression, in vitro RPFs were obtained for all PFASs. In vitro RPFs were found to correlate in general well with each other (Spearman correlation) except for the PPAR target genes ANGPTL4 and PDK4. Comparison of in vitro RPFs with RPFs obtained from in vivo studies in rats indicate that best correlations (Spearman correlation) were obtained for in vitro RPFs based on OAT5 and CXCL10 expression changes and external in vivo RPFs. HFPO-TA was found to be the most potent PFAS tested, being around tenfold more potent than PFOA. Altogether, it may be concluded that the HepaRG model may provide relevant data to provide insight into which PFASs are relevant regarding their hepatotoxic effects and that it can be applied as a screening tool to prioritize other PFASs for further hazard and risk assessment.

Human risk associated with exposure to mixtures of antiandrogenic chemicals evaluated using in vitro hazard and human biomonitoring data

BACKGROUND

Scientific evidence for underestimated toxicity from unintentional exposure to chemical mixtures is mounting. Yet, harmonized approaches on how to assess the actual risk of mixtures is lacking. As part of the European Joint programme 'Human Biomonitoring for Europe' we explored a novel methodology for mixture risk assessment of chemicals affecting male reproductive function.

METHODOLOGY

We explored a methodology for chemical mixture risk assessment based on human in vitro data combined with human exposure data, thereby circumventing the drawbacks of using hazard data from rodents and estimated exposure intake levels. Human androgen receptor (hAR) antagonism was selected as the most important molecular initiating event linked to adverse outcomes on male reproductive health.

RESULTS

Our work identified 231 chemicals able to interfere with hAR activity. Among these were 61 finally identified as having both reliable hAR antagonist and human biomonitoring data. Calculation of risk quotients indicated that PCBs (118, 138, 157), phthalates (BBP, DBP, DIBP), benzophenone-3, PFOS, methylparaben, triclosan, some pesticides (i.e cypermethrin, β-endosulfan, methylparathion, p,p-DDE), and a PAH metabolite (1-hydroxypyrene) contributed to the mixture effect. The major chemical mixture drivers were PCB 118, BBP, PFOS, DBP, and the UV filter benzophenone-3, together contributing with 75% of the total mixture effect that was primarily driven by high exposure values.

CONCLUSIONS

This viable way forward for mixture risk assessment of chemicals has the advantages of (1) being a more comprehensive mixture risk assessment also covering data-poor chemicals, and (2) including human data only. However, the approach is subjected to uncertainties in terms of in vitro to in vivo extrapolation, it is not ready for decision making, and needs further development. Still, the results indicate a concern for adverse effects on reproductive function in highly exposed boys, especially when considering additional exposure to data-poor chemicals and chemicals acting by other mechanisms of action.

Understanding the Chemical Exposome During Fetal Development and Early Childhood: A Review

Early human life is considered a critical window of susceptibility to external exposures. Infants are exposed to a multitude of environmental factors, collectively referred to as the exposome. The chemical exposome can be summarized as the sum of all xenobiotics that humans are exposed to throughout a lifetime. We review different exposure classes and routes that impact fetal and infant metabolism and the potential toxicological role of mixture effects. We also discuss the progress in human biomonitoring and present possiblemodels for studying maternal-fetal transfer. Data gaps on prenatal and infant exposure to xenobiotic mixtures are identified and include natural biotoxins, in addition to commonly reported synthetic toxicants, to obtain a more holistic assessment of the chemical exposome. We highlight the lack of large-scale studies covering a broad range of xenobiotics. Several recommendations to advance our understanding of the early-life chemical exposome and the subsequent impact on health outcomes are proposed.

Application of AOPs to assist regulatory assessment of chemical risks - Case studies, needs and recommendations

While human regulatory risk assessment (RA) still largely relies on animal studies, new approach methodologies (NAMs) based on in vitro, in silico or non-mammalian alternative models are increasingly used to evaluate chemical hazards. Moreover, human epidemiological studies with biomarkers of effect (BoE) also play an invaluable role in identifying health effects associated with chemical exposures. To move towards the next generation risk assessment (NGRA), it is therefore crucial to establish bridges between NAMs and standard approaches, and to establish processes for increasing mechanistically-based biological plausibility in human studies. The Adverse Outcome Pathway (AOP) framework constitutes an important tool to address these needs but, despite a significant increase in knowledge and awareness, the use of AOPs in chemical RA remains limited. The objective of this paper is to address issues related to using AOPs in a regulatory context from various perspectives as it was discussed in a workshop organized within the European Union partnerships HBM4EU and PARC in spring 2022. The paper presents examples where the AOP framework has been proven useful for the human RA process, particularly in hazard prioritization and characterization, in integrated approaches to testing and assessment (IATA), and in the identification and validation of BoE in epidemiological studies. Nevertheless, several limitations were identified that hinder the optimal usability and acceptance of AOPs by the regulatory community including the lack of quantitative information on response-response relationships and of efficient ways to map chemical data (exposure and toxicity) onto AOPs. The paper summarizes suggestions, ongoing initiatives and third-party tools that may help to overcome these obstacles and thus assure better implementation of AOPs in the NGRA.

Internal relative potency factors based on immunotoxicity for the risk assessment of mixtures of per- and polyfluoroalkyl substances (PFAS) in human biomonitoring

Relative potency factors (RPFs) for per- and polyfluoroalkyl substances (PFAS) have previously been derived based on liver effects in rodents for the purpose of performing mixture risk assessment with primary input from biomonitoring studies. However, in 2020, EFSA established a tolerable weekly intake for four PFAS assuming equal toxic potency for immune suppressive effects in humans. In this study we explored the possibility of deriving RPFs for immune suppressive effects using available data in rodents and humans. Lymphoid organ weights, differential blood cell counts, and clinical chemistry from 28-day studies in male rats from the National Toxicology Program (NTP) were combined with modeled serum PFAS concentrations to derive internal RPFs by applying dose-response modelling. Identified functional studies used diverse protocols and were not suitable for derivation of RPFs but were used to support immunotoxicity of PFAS in a qualitative manner. Furthermore, a novel approach was used to estimate internal RPFs based on epidemiological data by dose-response curve fitting optimization, looking at serum antibody concentrations and key cell populations from the National Health and Nutrition Examination Survey (NHANES). Internal RPFs were successfully derived for PFAS based on rat thymus weight, spleen weight, and globulin concentration. The available dose-response information for blood cell counts did not show a significant trend. Immunotoxic potency in serum was determined in the order PFDA > PFNA > PFHxA > PFOS > PFBS > PFOA > PFHxS. The epidemiological data showed inverse associations for the sum of PFOA, PFNA, PFHxS, and PFOS with serum antibody concentrations to mumps and rubella, but the data did not allow for deduction of reliable internal RPF estimates. The internal RPFs for PFAS based on decreased rat lymphoid organ weights are similar to those previously established for increased rat liver weight, strengthening the confidence in the overall applicability of these RPFs.

Trichloroethylene: An Invisible Cause of Parkinson's Disease?

The etiologies of Parkinson's disease (PD) remain unclear. Some, such as certain genetic mutations and head trauma, are widely known or easily identified. However, these causes or risk factors do not account for the majority of cases. Other, less visible factors must be at play. Among these is a widely used industrial solvent and common environmental contaminant little recognized for its likely role in PD: trichloroethylene (TCE). TCE is a simple, six-atom molecule that can decaffeinate coffee, degrease metal parts, and dry clean clothes. The colorless chemical was first linked to parkinsonism in 1969. Since then, four case studies involving eight individuals have linked occupational exposure to TCE to PD. In addition, a small epidemiological study found that occupational or hobby exposure to the solvent was associated with a 500% increased risk of developing PD. In multiple animal studies, the chemical reproduces the pathological features of PD.Exposure is not confined to those who work with the chemical. TCE pollutes outdoor air, taints groundwater, and contaminates indoor air. The molecule, like radon, evaporates from underlying soil and groundwater and enters homes, workplaces, or schools, often undetected. Despite widespread contamination and increasing industrial, commercial, and military use, clinical investigations of TCE and PD have been limited. Here, through a literature review and seven illustrative cases, we postulate that this ubiquitous chemical is contributing to the global rise of PD and that TCE is one of its invisible and highly preventable causes. Further research is now necessary to examine this hypothesis.

Characterization of different contaminants and current knowledge for defining chemical mixtures in human milk: A review

Hundreds of xenobiotics, with very diverse origins, have been detected in human milk, including contaminants of emerging concern, personal care products and other current-use substances reflecting lifestyle. The routes of exposure to these chemicals include dermal absorption, ingestion and inhalation. Specific families of chemicals are dominant among human milk monitoring studies (e.g., organochlorine pesticides, bisphenol A, dioxins), even though other understudied families may be equally toxicologically relevant (e.g., food-processing chemicals, current-use plasticizers and flame retardants, mycotoxins). Importantly, the lack of reliable human milk monitoring data for some individual chemicals and, especially, for complex mixtures, is a major factor hindering risk assessment. Non-targeted screening can be used as an effective tool to identify unknown contaminants of concern in human milk. This approach, in combination with novel methods to conduct risk assessments on the chemical mixtures detected in human milk, will assist in elucidating exposures that may have adverse effects on the development of breastfeeding infants.

Cumulative maternal and neonatal effects of combined exposure to a mixture of perfluorooctanoic acid (PFOA) and perfluorooctane sulfonic acid (PFOS) during pregnancy in the Sprague-Dawley rat

Globally, biomonitoring data demonstrate virtually all humans carry residues of multiple per- and polyfluoroalkyl substances (PFAS). Despite pervasive co-exposure, limited mixtures-based in vivo PFAS toxicity research has been conducted. Perfluorooctanoic acid (PFOA) and perfluorooctane sulfonic acid (PFOS) are commonly detected PFAS in human and environmental samples and both produce adverse effects in laboratory animal studies, including maternal and offspring effects when orally administered during pregnancy and lactation. To evaluate the effects of combined exposure to PFOA and PFOS, we orally exposed pregnant Sprague-Dawley rats from gestation day 8 (GD8) to postnatal day 2 (PND2) to PFOA (10-250 mg/kg/d) or PFOS (0.1-5 mg/kg/d) individually to characterize effects and dose response curve parameters, followed by a variable-ratio mixture experiment with a constant dose of PFOS (2 mg/kg/d) mixed with increasing doses of PFOA (3-80 mg/kg/d). The mixture study design was intended to: 1) shift the PFOA dose response curves for endpoints shared with PFOS, 2) allow comparison of dose addition (DA) and response addition (RA) model predictions, 3) conduct relative potency factor (RPF) analysis for multiple endpoints, and 4) avoid overt maternal toxicity. Maternal serum and liver concentrations of PFOA and PFOS were consistent between the individual chemical and mixture experiments. Combined exposure with PFOS significantly shifted the PFOA dose response curves towards effects at lower doses compared to PFOA-only exposure for multiple endpoints and these effects were well predicted by dose addition. For endpoints amenable to mixture model analyses, DA produced equivalent or better estimates of observed data than RA. All endpoints evaluated were accurately predicted by RPF and DA approaches except for maternal gestational weight gain, which produced less-than-additive results in the mixture. Data support the hypothesis of cumulative effects on shared endpoints from PFOA and PFOS co-exposure and dose additive approaches for predictive estimates of mixture effects.

Exploring the relationship between blood toxic metal(oid)s and serum insulin levels through benchmark modelling of human data: Possible role of arsenic as a metabolic disruptor

The major goal of this study was to estimate the correlations and dose-response pattern between the measured blood toxic metals (cadmium (Cd), mercury (Hg), chromium (Cr), nickel (Ni))/metalloid (arsenic (As)) and serum insulin level by conducting Benchmark dose (BMD) analysis of human data. The study involved 435 non-occupationally exposed individuals (217 men and 218 women). The samples were collected at health care institutions in Belgrade, Serbia, from January 2019 to May 2021. Blood sample preparation was conducted by microwave digestion. Cd was measured by graphite furnace atomic absorption spectrophotometry (GF-AAS), while inductively coupled plasma-mass spectrometry (ICP-MS) was used to measure Hg, Ni, Cr and As. BMD analysis of insulin levels represented as quantal data was done using the PROAST software version 70.1 (model averaging methodology, BMD response: 10%). In the male population, there was no correlation between toxic metal/metalloid concentrations and insulin level. However, in the female population/whole population, a high positive correlation for As and Hg, and a strong negative correlation for Ni and measured serum insulin level was established. BMD modelling revealed quantitative associations between blood toxic metal/metalloid concentrations and serum insulin levels. All the estimated BMD intervals were wide except the one for As, reflecting a high degree of confidence in the estimations and possible role of As as a metabolic disruptor. These results indicate that, in the case of As blood concentrations, even values higher than BMD (BMDL): 3.27 (1.26) (male population), 2.79 (0.771) (female population), or 1.18 (2.96) μg/L (whole population) might contribute to a 10% higher risk of insulin level alterations, meaning 10% higher risk of blood insulin increasing from within reference range to above reference range. The obtained results contribute to the current body of knowledge on the use of BMD modelling for analysing human data.

Retrospective cumulative dietary risk assessment of craniofacial alterations by residues of pesticides

EFSA established cumulative assessment groups and conducted retrospective cumulative risk assessments for two types of craniofacial alterations (alterations due to abnormal skeletal development, head soft tissue alterations and brain neural tube defects) for 14 European populations of women in childbearing age. Cumulative acute exposure calculations were performed by probabilistic modelling using monitoring data collected by Member States in 2017, 2018 and 2019. A rigorous uncertainty analysis was performed using expert knowledge elicitation. Considering all sources of uncertainty, their dependencies and differences between populations, it was concluded with varying degrees of certainty that the MOET resulting from cumulative exposure is above 100 for the two types of craniofacial alterations. The threshold for regulatory consideration established by risk managers is therefore not exceeded. Considering the severity of the effects under consideration, it was also assessed whether the MOET is above 500. This was the case with varying levels of certainty for the head soft tissue alterations and brain neural tube defects. However, for the alterations due to abnormal skeletal development, it was found about as likely as not that the MOET is above 500 in most populations. For two populations, it was even found more likely that the MOET is below 500. These results were discussed in the light of the conservatism of the methodological approach.

Benchmark dose approach in investigating the relationship between blood metal levels and reproductive hormones: Data set from human study

The main objective of this research was to conduct a dose-response modeling between the internal dose of measured blood Cd, As, Hg, Ni, and Cr and hormonal response of serum testosterone, luteinizing hormone (LH), and follicle-stimulating hormone (FSH). The study included 207 male participants from subjects of 5 different cohorts (patients with prostate, testicular, and pancreatic cancer, patients suffering from various thyroid and metabolic disorders, as well as healthy volunteers), enrolled from January 2019 to May 2021 at the Clinical Centre of Serbia in Belgrade, Serbia. Benchmark dose-response modeling analysis was performed with the PROAST software version 70.1, showing the hormone levels as quantal data. The averaging technique was applied to compute the Benchmark dose (BMD) interval (BMDI), with benchmark response set at 10%. Dose-response relationships between metal/metalloid blood concentration and serum hormone levels were confirmed for all the investigated metals/metalloid and hormones. The narrowest BMDI was found for Cd-testosterone and Hg-LH pairs, indicative of high confidence in these estimates. Although further research is needed, the observed findings demonstrate that the BMD approach may prove to be significant in the dose-response modeling of human data.