Model systems and organisms for addressing inter- and intra-species variability in risk assessment

The effects of a glyphosate-based herbicide on the bovine gametes during an in vitro embryo production model

Roundup® (R), while it is the most used herbicide globally, and its residues are ubiquitous in urban and suburban areas, its impact on vertebrates' safety remains highly debated. Here, in three in vitro experiments, we investigated the effects of a very low dose (1 ppm) of R on the fertilization capacity and embryo development in cattle. In the first experiment, frozen-thawed bull semen exposed to R for 1 h exhibited reduced motility parameters but unaffected fertilization ability. However, after in vitro fertilization, the rates of embryo formation were significantly lower compared to the untreated controls. In the second experiment, oocytes exposed to R during in vitro maturation showed reduced cleavage rates, and the embryo yield on days 7, 8, and 9 of embryo culture was significantly lower than that of the controls. In the third experiment, oocytes were matured in the presence of R and in a medium containing both R and Zinc, chosen to offer antioxidant protection to the oocytes. Day-7 blastocysts were analyzed for the expression of genes associated with oxidative stress, apoptosis, and epigenetic reprogramming. Exposure to R markedly suppressed embryo formation rates compared to the controls. The combination of R with Zinc restored the blastocyst yield, which on days 8 and 9 was comparable to that of the controls and higher than the groups exposed only to R on all days. The gene expression analysis revealed that R promotes oxidative stress development, triggers apoptosis, and induces epigenetic changes in developing embryos, while zinc presence alleviates these adverse effects of R. These findings imply that even at very low doses, R could be highly toxic, leading to functional abnormalities in both gametes, potentially affecting fertility in both genders.

Advancing probabilistic risk assessment by integrating human biomonitoring, new approach methods, and Bayesian modeling: A case study with the mycotoxin deoxynivalenol

Deoxynivalenol (DON) is a mycotoxin frequently observed in cereals and cereal-based foods, with reported toxicological effects including reduced body weight, immunotoxicity and reproductive defects. The European Food Safety Authority used traditional risk assessment approaches to derive a deterministic Tolerable Daily Intake (TDI) of 1 μg/kg-day, however data from human biomarkers studies indicate widespread and variable exposure worldwide, necessitating more sophisticated and advanced methods to quantify population risk. The World Health Organization/International Programme on Chemical Safety (WHO/IPCS) has previously used DON as a case example in replacing the TDI with a probabilistic toxicity value, using default uncertainty and variability distributions to derive the Human Dose corresponding to an effect size M in the Ith percentile of the population (HDMI) for M = 5 % decrease in body weight and I = 1 %. In this study, we extend this case study by incorporating (1) Bayesian modeling approaches, (2) using both in vivo data and in vitro population new approach methods to replace default distributions for interspecies toxicokinetic (TK) differences and intraspecies TK and toxicodynamic (TD) variability, and (3) integrating biomonitoring data and probabilistic dose-response functions to characterize population risk distributions. We first derive an HDMI of 5.5 [1.4-24] μg/kg-day, also using TK modeling to converted the HDMI to Biomonitoring Equivalents, BEMI for comparison with biomonitoring data, with a blood BEMI of 0.53 [0.17-1.6] μg/L and a urinary excretion BEMI of 3.9 [1.0-16] μg/kg-day. We then illustrate how this integrative approach can advance quantitative risk characterization using two human biomonitoring datasets, estimating both the fraction of population with an effect size M ≥ 5 % as well as the distribution of effect sizes. Overall, we demonstrate that integration of Bayesian modeling, human biomonitoring data, and in vitro population-based TD data within the WHO/IPCS probabilistic framework yields more accurate, precise, and comprehensive risk characterization.

Seminar: Extracellular Vesicles as Mediators of Environmental Stress in Human Disease

BACKGROUND

Extracellular vesicles (EVs), membrane-bound particles containing a variety of RNA types, DNA, proteins, and other macromolecules, are now appreciated as an important means of communication between cells and tissues, both in normal cellular physiology and as a potential indicator of cellular stress, environmental exposures, and early disease pathogenesis. Extracellular signaling through EVs is a growing field of research for understanding fundamental mechanisms of health and disease and for the potential for biomarker discovery and therapy development. EVs are also known to play important roles in mediating the effects of exposure to environmental stress.

OBJECTIVES

This seminar addresses the application of new tools and approaches for EV research, developed in part through the National Institutes of Health (NIH) Extracellular RNA Communication Program, and reflects presentations and discussions from a workshop held 27-28 September 2021 by the National Institute of Environmental Health Sciences (NIEHS) and the National Center for Advancing Translational Sciences (NCATS) on "Extracellular Vesicles, Exosomes, and Cell-Cell Signaling in Response to Environmental Stress." The panel of experts discussed current research on EVs and environmental exposures, highlighted recent advances in EV isolation and characterization, and considered research gaps and opportunities toward identifying and characterizing the roles for EVs in environmentally related diseases, as well as the current challenges and opportunities in this field.

DISCUSSION

The authors discuss the application of new experimental models, particularly organ-on-chip (OOC) systems and in vitro approaches and how these have the potential to extend findings in population-based studies of EVs in exposure-related diseases. Given the complex challenges of identifying cell-specific EVs related to environmental exposures, as well as the general heterogeneity and variability in EVs in blood and other accessible biological samples, there is a critical need for rigorous reporting of experimental methods and validation studies. The authors note that these efforts, combined with cross-disciplinary approaches, would ensure that future research efforts in environmental health studies on EV biomarkers are rigorous and reproducible. https://doi.org/10.1289/EHP12980.

Refining risk estimates for lead in drinking water based on the impact of genetics and diet on blood lead levels using the Collaborative Cross mouse population

Blood lead (Pb) level (BLL) is a commonly used biomarker to evaluate associations with health effects. However, interventions to reduce the adverse effects of Pb require relating BLL to external exposure. Moreover, risk mitigation actions need to ensure protection of more susceptible individuals with a greater tendency to accumulate Pb. Because little data is available to quantify inter-individual variability in biokinetics of Pb, we investigated the influence of genetics and diet on BLL in the genetically diverse Collaborative Cross (CC) mouse population. Adult female mice from 49 CC strains received either a standard mouse chow or a chow mimicking the American diet while being provided water ad libitum with 1000 ppm Pb for 4 weeks. In both arms of the study, inter-strain variability was observed; however, in American diet-fed animals, the BLL was greater and more variable. Importantly, the degree of variation in BLL among strains on the American diet was greater (2.3) than the default variability estimate (1.6) used in setting the regulatory standards. Genetic analysis identified suggestive diet-associated haplotypes that were associated with variation in BLL, largely contributed by the PWK/PhJ strain. This study quantified the variation in BLL that is due to genetic background, diet, and their interactions, and observed that it may be greater than that assumed for current regulatory standards for Pb in drinking water. Moreover, this work highlights the need of characterizing inter-individual variation in BLL to ensure adequate public health interventions aimed at reducing human health risks from Pb.

G × E interactions as a basis for toxicological uncertainty

To transfer toxicological findings from model systems, e.g. animals, to humans, standardized safety factors are applied to account for intra-species and inter-species variabilities. An alternative approach would be to measure and model the actual compound-specific uncertainties. This biological concept assumes that all observed toxicities depend not only on the exposure situation (environment = E), but also on the genetic (G) background of the model (G × E). As a quantitative discipline, toxicology needs to move beyond merely qualitative G × E concepts. Research programs are required that determine the major biological variabilities affecting toxicity and categorize their relative weights and contributions. In a complementary approach, detailed case studies need to explore the role of genetic backgrounds in the adverse effects of defined chemicals. In addition, current understanding of the selection and propagation of adverse outcome pathways (AOP) in different biological environments is very limited. To improve understanding, a particular focus is required on modulatory and counter-regulatory steps. For quantitative approaches to address uncertainties, the concept of "genetic" influence needs a more precise definition. What is usually meant by this term in the context of G × E are the protein functions encoded by the genes. Besides the gene sequence, the regulation of the gene expression and function should also be accounted for. The widened concept of past and present "gene expression" influences is summarized here as Ge. Also, the concept of "environment" needs some re-consideration in situations where exposure timing (Et) is pivotal: prolonged or repeated exposure to the insult (chemical, physical, life style) affects Ge. This implies that it changes the model system. The interaction of Ge with Et might be denoted as Ge × Et. We provide here general explanations and specific examples for this concept and show how it could be applied in the context of New Approach Methodologies (NAM).

Evaluation of Sex Differences in Preclinical Pharmacology Research: How Far Is Left to Go?

Until the last quarter of the 20th century, sex was not recognized as a variable in health research, nor was it believed to be a factor that could affect health and illness. Researchers preferred studying male models for a variety of reasons, such as simplicity, lower costs, hormone confounding effects, and fear of liability from perinatal exposure in case of pregnancy. Equitable representation is imperative for determining the safety, effectiveness, and tolerance of therapeutic agents for all consumers. Decades of female models' underrepresentation in preclinical studies has resulted in inequality in the understanding, diagnosis, and treatment of disease between the sexes. Sex bias has been highlighted as one of the contributing factors to the poor translation and replicability of preclinical research. There have been multiple calls for action, and the inclusion of sex as a biological variable is increasingly supported. However, although there has been substantial progress in the efforts to include more female models in preclinical studies, disparities today remain. In the present review, we consider the current standard practice of the preclinical research setting, why the sex bias exists, why there is the need to include female models, and what risks may arise from continuing this exclusion from experimental design.

Risk Assessment of Transgender People: Development of Rodent Models Mimicking Gender-Affirming Hormone Therapies and Identification of Sex-Dimorphic Liver Genes as Novel Biomarkers of Sex Transition

Transgender (TG) describes individuals whose gender identity differs from the social norms. TG people undergoing gender-affirming hormone therapy (HT) may be considered a sub-group of the population susceptible to environmental contaminants for their targets and modes of action. The aim of this study is to set appropriate HT doses and identify specific biomarkers to implement TG animal models. Four adult rats/group/sex were subcutaneously exposed to three doses of HT (plus control) selected starting from available data. The demasculinizing-feminizing models (dMF) were β-estradiol plus cyproterone acetate, at 0.09 + 0.33, 0.09 + 0.93 and 0.18 + 0.33 mg, respectively, five times/week. The defeminizing-masculinizing models (dFM) were testosterone (T) at 0.45, 0.95 and 2.05 mg, two times/week. Clitoral gain and sperm count, histopathological analysis of reproductive organs and liver, hormone serum levels and gene expression of sex-dimorphic CYP450 were evaluated. In the dMF model, the selected doses-leading to T serum levels at the range of the corresponding cisgender-induced strong general toxicity and cannot be used in long-term studies. In the dFM model, 0.45 mg of T represents the correct dose. In addition, the endpoints selected are considered suitable and reliable to implement the animal model. The sex-specific CYP expression is a suitable biomarker to set proper (de)masculinizing/(de)feminizing HT and to implement TG animal models.

Cumulative Risk Meets Inter-Individual Variability: Probabilistic Concentration Addition of Complex Mixture Exposures in a Population-Based Human In Vitro Model

Although humans are continuously exposed to complex chemical mixtures in the environment, it has been extremely challenging to investigate the resulting cumulative risks and impacts. Recent studies proposed the use of “new approach methods,” in particular in vitro assays, for hazard and dose−response evaluation of mixtures. We previously found, using five human cell-based assays, that concentration addition (CA), the usual default approach to calculate cumulative risk, is mostly accurate to within an order of magnitude. Here, we extend these findings to further investigate how cell-based data can be used to quantify inter-individual variability in CA. Utilizing data from testing 42 Superfund priority chemicals separately and in 8 defined mixtures in a human cell-based population-wide in vitro model, we applied CA to predict effective concentrations for cytotoxicity for each individual, for “typical” (median) and “sensitive” (first percentile) members of the population, and for the median-to-sensitive individual ratio (defined as the toxicodynamic variability factor, TDVF). We quantified the accuracy of CA with the Loewe Additivity Index (LAI). We found that LAI varies more between different mixtures than between different individuals, and that predictions of the population median are generally more accurate than predictions for the “sensitive” individual or the TDVF. Moreover, LAI values were generally <1, indicating that the mixtures were more potent than predicted by CA. Together with our previous studies, we posit that new approach methods data from human cell-based in vitro assays, including multiple phenotypes in diverse cell types and studies in a population-wide model, can fill critical data gaps in cumulative risk assessment, but more sophisticated models of in vitro mixture additivity and bioavailability may be needed. In the meantime, because simple CA models may underestimate potency by an order of magnitude or more, either whole-mixture testing in vitro or, alternatively, more stringent benchmarks of cumulative risk indices (e.g., lower hazard index) may be needed to ensure public health protection.

A Population-Based Human In Vitro Approach to Quantify Inter-Individual Variability in Responses to Chemical Mixtures

Human cell-based population-wide in vitro models have been proposed as a strategy to derive chemical-specific estimates of inter-individual variability; however, the utility of this approach has not yet been tested for cumulative exposures in mixtures. This study aimed to test defined mixtures and their individual components and determine whether adverse effects of the mixtures were likely to be more variable in a population than those of the individual chemicals. The in vitro model comprised 146 human lymphoblastoid cell lines from four diverse subpopulations of European and African descent. Cells were exposed, in concentration−response, to 42 chemicals from diverse classes of environmental pollutants; in addition, eight defined mixtures were prepared from these chemicals using several exposure- or hazard-based scenarios. Points of departure for cytotoxicity were derived using Bayesian concentration−response modeling and population variability was quantified in the form of a toxicodynamic variability factor (TDVF). We found that 28 chemicals and all mixtures exhibited concentration−response cytotoxicity, enabling calculation of the TDVF. The median TDVF across test substances, for both individual chemicals or defined mixtures, ranged from a default assumption (101/2) of toxicodynamic variability in human population to >10. The data also provide a proof of principle for single-variant genome-wide association mapping for toxicity of the chemicals and mixtures, although replication would be necessary due to statistical power limitations with the current sample size. This study demonstrates the feasibility of using a set of human lymphoblastoid cell lines as an in vitro model to quantify the extent of inter-individual variability in hazardous properties of both individual chemicals and mixtures. The data show that population variability of the mixtures is unlikely to exceed that of the most variable component, and that similarity in genome-wide associations among components may be used to accrue additional evidence for grouping of constituents in a mixture for cumulative assessments.