Mining of Consumer Product Ingredient and Purchasing Data to Identify Potential Chemical Coexposures

Time-dependent hormesis transfer from five high-frequency personal care product components to mixtures

There is potential for personal care products (PCPs) components and mixtures to induce hormesis. How hormesis is related to time and transmitted from components to mixtures are not clear. In this paper, we conducted determination of components in 16 PCP products and then ran frequent itemset mining on the component data. Five high-frequency components (HFCs), betaine (BET), 1,3-butanediol (BUT), ethylenediaminetetraacetic acid disodium salt (EDTA), glycerol (GLO), and phenoxyethanol (POE), and 14 mixtures were identified. For each mixture system, one mixture ray with the actual mixture ratios in the products was selected. Time-dependent microplate toxicity analysis was used to test the luminescence inhibition toxicity of five HFCs and 14 mixture rays to Vibrio qinghaiensis sp.-Q67 at 12 concentration gradients and eight exposure times. It is showed that BET, EDTA, POE, and 13 mixture rays containing at least one J-type component showed time-dependent hormesis. Characteristic parameters used to describe hormesis revealed that the absolute value of the maximum stimulatory effect (|Emin|) generally increased with time. Notably, mixtures composed of POE and S-type components showed greater |Emin| than POE alone at the same time. Importantly, the maximum stimulatory effective concentration, NOEC/the zero effective concentration point, and EC50 remained relatively stable. Nine hormesis transmission phenomena were observed in different mixture rays. While all mixtures primarily exhibited additive action, varying degrees of synergism and antagonism were noted in binary mixtures, with no strong synergism or antagonism observed in ternary and quaternary mixtures. These findings offer valuable insights for the screening of HFCs and their mixtures, as well as the study of hormesis transmission in personal care products.

Characterizing Chemical Exposure Trends from NHANES Urinary Biomonitoring Data

BACKGROUND

Xenobiotic metabolites are widely present in human urine and can indicate recent exposure to environmental chemicals. Proper inference of which chemicals contribute to these metabolites can inform human exposure and risk. Furthermore, longitudinal biomonitoring studies provide insight into how chemical exposures change over time.

OBJECTIVES

We constructed an exposure landscape for as many human-exposure relevant chemicals over as large a time span as possible to characterize exposure trends across demographic groups and chemical types.

METHODS

We analyzed urine data of nine 2-y cohorts (1999-2016) from the National Health and Nutrition Examination Survey (NHANES). Chemical daily intake rates (in milligrams per kilogram bodyweight per day) were inferred, using the R package bayesmarker, from metabolite concentrations in each cohort individually to identify exposure trends. Trends for metabolites and parents were clustered to find chemicals with similar exposure patterns. Exposure variation by age, gender, and body mass index were also assessed.

RESULTS

Intake rates for 179 parent chemicals were inferred from 151 metabolites (96 measured in five or more cohorts). Seventeen metabolites and 44 parent chemicals exhibited fold-changes ≥ 10 between any two cohorts (deltamethrin, di-n -octyl phthalate, and di-isononyl phthalate had the greatest exposure increases). Di-2-ethylhexyl phthalate intake began decreasing in 2007, whereas both di-isobutyl and di-isononyl phthalate began increasing shortly before. Intake for four parabens was markedly higher in females, especially reproductive-age females, compared with males and children. Cadmium and arsenobetaine exhibited higher exposure for individuals > 65 years of age and lower for individuals < 20 years of age.

DISCUSSION

With appropriate analysis, NHANES indicates trends in chemical exposures over the past two decades. Decreases in exposure are observable as the result of regulatory action, with some being accompanied by increases in replacement chemicals. Age- and gender-specific variations in exposure were observed for multiple chemicals. Continued estimation of demographic-specific exposures is needed to both monitor and identify potential vulnerable populations. https://doi.org/10.1289/EHP12188.

Time-dependent hormetic dose responses of skin care product mixtures to Vibrio qinghaiensis sp.-Q67: Appearance and quantification

Hormesis is a widely recognized and extensively studied phenomenon. However, few studies have described the quantitative characteristics of hormesis required for appropriate risk assessment. Although skin care product (SCP) mixtures and their active ingredients can induce the hormesis of Vibrio qinghaiensis sp.-Q67 (Q67), the quantitative characteristics of time-dependent hormetic dose responses in SCPs have not yet been investigated. In this study, 28 SCP mixtures were tested for luminescence toxicity against Q67 after five exposure durations (0.25, 3, 6, 9, and 12 h). With increasing exposure duration, the concentration response curves (CRCs) were classified as constant monotonic nonlinear (S-shaped) for four SCPs, S- to hormetic (J-shaped) for 13 SCPs, and constant J-shaped for 11 SCPs. Of 140 CRCs, 98 were J-shaped. An increased frequency of SCPs inducing hormesis was observed. The toxicity (pEC50) of the SCPs was independent of the exposure duration and product type. The maximum stimulatory effect (Emin) of the 12 SCPs increased with exposure duration. We proposed a modified parameter, the width of inhibition dose zone (WIDZ; EC50/EC10), to depict the width of inhibition dose zone. The WIDZ of S-shaped CRCs were significantly larger than that of J-shaped CRCs. In addition, the characteristic parameters reported in the general literature were analyzed. The good linear relationship between EC50 and the maximum stimulatory effective concentration (ECmin) indicated that toxicity may be transformed into stimulatory effects over exposure durations. The width of stimulation dose zone (WSDZ) and Emin of the seven SCPs had the same increasing trends with increasing exposure duration. The combination of WIDZ with other characteristic parameters (e.g., zero effective concentration point, ECmin, etc.) could better depict hormesis with low-dose stimulation and high-dose inhibition. The quantitative characteristics of the dose-responses of hormesis-inducing SCPs could provide reference basis for the risk assessment of SCP mixtures.

Assessing household lifestyle exposures from consumer purchases, the My Purchases cohort

Consumer purchase data (CPD) is a promising instrument to assess the impact of purchases on health, but is limited by the need for manual scanning, a lack of access to data from multiple retailers, and limited information on product data and health outcomes. Here we describe the My Purchases cohort, a web-app enabled, prospective collection of CPD, covering several large retail chains in Denmark, that enables linkage to health outcomes. The cohort included 459 participants as of July 03, 2023. Up to eight years of CPD have been collected, with 2,225,010 products purchased, comprising 223,440 unique products. We matched 88.5% of all products by product name or item number to one generic food database and three product databases. Combined, the databases enable analysis of key exposures such as nutrients, ingredients, or additives. We found that increasing the number of retailers that provide CPD for each consumer improved the stability of individual CPD profiles and when we compared kilojoule information from generic and specific product matches, we found a median modified relative difference of 0.23. Combined with extensive product databases and health outcomes, CPD could provide the basis for extensive investigations of how what we buy affects our health.

Cutting-edge computational chemical exposure research at the U.S. Environmental Protection Agency

Exposure science is evolving from its traditional "after the fact" and "one chemical at a time" approach to forecasting chemical exposures rapidly enough to keep pace with the constantly expanding landscape of chemicals and exposures. In this article, we provide an overview of the approaches, accomplishments, and plans for advancing computational exposure science within the U.S. Environmental Protection Agency's Office of Research and Development (EPA/ORD). First, to characterize the universe of chemicals in commerce and the environment, a carefully curated, web-accessible chemical resource has been created. This DSSTox database unambiguously identifies >1.2 million unique substances reflecting potential environmental and human exposures and includes computationally accessible links to each compound's corresponding data resources. Next, EPA is developing, applying, and evaluating predictive exposure models. These models increasingly rely on data, computational tools like quantitative structure activity relationship (QSAR) models, and machine learning/artificial intelligence to provide timely and efficient prediction of chemical exposure (and associated uncertainty) for thousands of chemicals at a time. Integral to this modeling effort, EPA is developing data resources across the exposure continuum that includes application of high-resolution mass spectrometry (HRMS) non-targeted analysis (NTA) methods providing measurement capability at scale with the number of chemicals in commerce. These research efforts are integrated and well-tailored to support population exposure assessment to prioritize chemicals for exposure as a critical input to risk management. In addition, the exposure forecasts will allow a wide variety of stakeholders to explore sustainable initiatives like green chemistry to achieve economic, social, and environmental prosperity and protection of future generations.

Online coupling of matrix solid-phase dispersion to direct analysis in real time mass spectrometry for high-throughput analysis of regulated chemicals in consumer products

The current work is the first study on online coupling of matrix solid-phase dispersion (MSPD) to direct analysis in real time mass spectrometry (DART-MS) bridging with solid-phase analytical derivatization (SPAD) based on a graphene oxide nanosheets (GONs)-coated cotton swab. Proof-of-concept demonstrations were explored for high-throughput analysis of a diversity of regulated chemicals in consumer products such as textiles, toys, and cosmetics. On-demand sorbent combinations were blended with samples, packed into MSPD columns, and mounted on a homemade 3D-printed rack module for automated sample feeding. To achieve good synergy between MSPD and DART-MS, a cotton swab with a conical tip deposited with GONs was attached to the bottom of the MSPD column. The swabs serve as a solid-phase microextraction probe for convenient enrichment of the eluted analytes from MSPD, thermal desorption of the enriched analytes by DART, and sensitive detection by a hybrid quadrupole-Orbitrap mass spectrometer. Furthermore, the utility of an on-swab SPAD strategy was demonstrated for the detection of formaldehyde by use of the derivatizing reagent of dansyl hydrazine, contributing to improved ionization efficiency without compromising the overall coherence of the analytical workflow. The MSPD-DART-MS methodology was systematically optimized and validated, obtaining acceptable recovery (71.7-110.3%), repeatability (11.8-19.3%), and sensitivity (limits of detection and quantitation in the ranges of 6.2-19.5 and 23.7-75.9 μg/kg) for 32 target analytes. The developed protocol streamlined sample extraction, clean-up, desorption, ionization, and detection, highlighting the appealing potential for high-throughput analysis of samples with complex matrices.

Prioritization of mixtures of neurotoxic chemicals for biomonitoring using high-throughput toxicokinetics and mixture toxicity modeling

Modern society continues to pollute the environment with larger quantities of chemicals that have also become more structurally and functionally diverse. Risk assessment of chemicals can hardly keep up with the sheer numbers that lead to complex mixtures of increasing chemical diversity including new chemicals, substitution products on top of still abundant legacy compounds. Fortunately, over the last years computational tools have helped us to identify and prioritize chemicals of concern. These include toxicokinetic models to predict exposure to chemicals as well as new approach methodologies such as in-vitro bioassays to address toxicodynamic effects. Combined, they allow for a prediction of mixtures and their respective effects and help overcome the lack of data we face for many chemicals. In this study we propose a high-throughput approach using experimental and predicted exposure, toxicokinetic and toxicodynamic data to simulate mixtures, to which a virtual population is exposed to and predict their mixture effects. The general workflow is adaptable for any type of toxicity, but we demonstrated its applicability with a case study on neurotoxicity. If no experimental data for neurotoxicity were available, we used baseline toxicity predictions as a surrogate. Baseline toxicity is the minimal toxicity any chemical has and might underestimate the true contribution to the mixture effect but many neurotoxicants are not by orders of magnitude more potent than baseline toxicity. Therefore, including baseline-toxic effects in mixture simulations yields a more realistic picture than excluding them in mixture simulations. This workflow did not only correctly identify and prioritize known chemicals of concern like benzothiazoles, organochlorine pesticides and plasticizers but we were also able to identify new potential neurotoxicants that we recommend to include in future biomonitoring studies and if found in humans, to also include in neurotoxicity screening.

Exposure forecasting - ExpoCast - for data-poor chemicals in commerce and the environment

Estimates of exposure are critical to prioritize and assess chemicals based on risk posed to public health and the environment. The U.S. Environmental Protection Agency (EPA) is responsible for regulating thousands of chemicals in commerce and the environment for which exposure data are limited. Since 2009 the EPA's ExpoCast ("Exposure Forecasting") project has sought to develop the data, tools, and evaluation approaches required to generate rapid and scientifically defensible exposure predictions for the full universe of existing and proposed commercial chemicals. This review article aims to summarize issues in exposure science that have been addressed through initiatives affiliated with ExpoCast. ExpoCast research has generally focused on chemical exposure as a statistical systems problem intended to inform thousands of chemicals. The project exists as a companion to EPA's ToxCast ("Toxicity Forecasting") project which has used in vitro high-throughput screening technologies to characterize potential hazard posed by thousands of chemicals for which there are limited toxicity data. Rapid prediction of chemical exposures and in vitro-in vivo extrapolation (IVIVE) of ToxCast data allow for prioritization based upon risk of adverse outcomes due to environmental chemical exposure. ExpoCast has developed (1) integrated modeling approaches to reliably predict exposure and IVIVE dose, (2) highly efficient screening tools for chemical prioritization, (3) efficient and affordable tools for generating new exposure and dose data, and (4) easily accessible exposure databases. The development of new exposure models and databases along with the application of technologies like non-targeted analysis and machine learning have transformed exposure science for data-poor chemicals. By developing high-throughput tools for chemical exposure analytics and translating those tools into public health decisions ExpoCast research has served as a crucible for identifying and addressing exposure science knowledge gaps.

Principal Component Pursuit for Pattern Identification in Environmental Mixtures

BACKGROUND

Environmental health researchers often aim to identify sources or behaviors that give rise to potentially harmful environmental exposures.

OBJECTIVE

We adapted principal component pursuit (PCP)-a robust and well-established technique for dimensionality reduction in computer vision and signal processing-to identify patterns in environmental mixtures. PCP decomposes the exposure mixture into a low-rank matrix containing consistent patterns of exposure across pollutants and a sparse matrix isolating unique or extreme exposure events.

METHODS

We adapted PCP to accommodate nonnegative data, missing data, and values below a given limit of detection (LOD). We simulated data to represent environmental mixtures of two sizes with increasing proportions < LOD and three noise structures. We applied PCP-LOD to evaluate its performance in comparison with principal component analysis (PCA). We next applied principal component pursuit with limit of detection (PCP-LOD) to an exposure mixture of 21 persistent organic pollutants (POPs) measured in 1,000 U.S. adults from the 2001-2002 National Health and Nutrition Examination Survey (NHANES). We applied singular value decomposition to the estimated low-rank matrix to characterize the patterns.

RESULTS

PCP-LOD recovered the true number of patterns through cross-validation for all simulations; based on an a priori specified criterion, PCA recovered the true number of patterns in 32% of simulations. PCP-LOD achieved lower relative predictive error than PCA for all simulated data sets with up to 50% of the data < LOD . When 75% of values were < LOD , PCP-LOD outperformed PCA only when noise was low. In the POP mixture, PCP-LOD identified a rank-three underlying structure and separated 6% of values as extreme events. One pattern represented comprehensive exposure to all POPs. The other patterns grouped chemicals based on known structure and toxicity.

DISCUSSION

PCP-LOD serves as a useful tool to express multidimensional exposures as consistent patterns that, if found to be related to adverse health, are amenable to targeted public health messaging. https://doi.org/10.1289/EHP10479.

Morphological profiling of environmental chemicals enables efficient and untargeted exploration of combination effects

Environmental chemicals are commonly studied one at a time, and there is a need to advance our understanding of the effect of exposure to their combinations. Here we apply high-content microscopy imaging of cells stained with multiplexed dyes (Cell Painting) to profile the effects of Cetyltrimethylammonium bromide (CTAB), Bisphenol A (BPA), and Dibutyltin dilaurate (DBTDL) exposure on four human cell lines; both individually and in all combinations. We show that morphological features can be used with multivariate data analysis to discern between exposures from individual compounds, concentrations, and combinations. CTAB and DBTDL induced concentration-dependent morphological changes across the four cell lines, and BPA exacerbated morphological effects when combined with CTAB and DBTDL. Combined exposure to CTAB and BPA induced changes in the ER, Golgi apparatus, nucleoli and cytoplasmic RNA in one of the cell lines. Different responses between cell lines indicate that multiple cell types are needed when assessing combination effects. The rapid and relatively low-cost experiments combined with high information content make Cell Painting an attractive methodology for future studies of combination effects. All data in the study is made publicly available on Figshare.

CCDB: A database for exploring inter-chemical correlations in metabolomics and exposomics datasets

Inter-chemical correlations in metabolomics and exposomics datasets provide valuable information for studying relationships among chemicals reported for human specimens. With an increase in the number of compounds for these datasets, a network graph analysis and visualization of the correlation structure is difficult to interpret. We have developed the Chemical Correlation Database (CCDB), as a systematic catalogue of inter-chemical correlation in publicly available metabolomics and exposomics studies. The database has been provided via an online interface to create single compound-centric views. We have demonstrated various applications of the database to explore: 1) the chemicals from a chemical class such as Per- and Polyfluoroalkyl Substances (PFAS), polycyclic aromatic hydrocarbons (PAHs), polychlorinated biphenyls (PCBs), phthalates and tobacco smoke related metabolites; 2) xenobiotic metabolites such as caffeine and acetaminophen; 3) endogenous metabolites (acyl-carnitines); and 4) unannotated peaks for PFAS. The database has a rich collection of 35 human studies, including the National Health and Nutrition Examination Survey (NHANES) and high-quality untargeted metabolomics datasets. CCDB is supported by a simple, interactive and user-friendly web-interface to retrieve and visualize the inter-chemical correlation data. The CCDB has the potential to be a key computational resource in metabolomics and exposomics facilitating the expansion of our understanding about biological and chemical relationships among metabolites and chemical exposures in the human body. The database is available at www.ccdb.idsl.me site.

Exposure to Contemporary and Emerging Chemicals in Commerce among Pregnant Women in the United States: The Environmental influences on Child Health Outcome (ECHO) Program

Prenatal chemical exposures can influence maternal and child health; however, few industrial chemicals are routinely biomonitored. We assessed an extensive panel of contemporary and emerging chemicals in 171 pregnant women across the United States (U.S.) and Puerto Rico in the Environmental influences on Child Health Outcomes (ECHO) Program. We simultaneously measured urinary concentrations of 89 analytes (103 total chemicals representing 73 parent compounds) in nine chemical groups: bactericides, benzophenones, bisphenols, fungicides and herbicides, insecticides, organophosphate esters (OPEs), parabens, phthalates/alternative plasticizers, and polycyclic aromatic hydrocarbons (PAHs). We estimated associations of creatinine-adjusted concentrations with sociodemographic and specimen characteristics. Among our diverse prenatal population (60% non-Hispanic Black or Hispanic), we detected 73 of 89 analytes in ≥1 participant and 36 in >50% of participants. Five analytes not currently included in the U.S. biomonitoring were detected in ≥90% of samples: benzophenone-1, thiamethoxam, mono-2-(propyl-6-carboxy-hexyl) phthalate, monocarboxy isooctyl phthalate, and monohydroxy-iso-decyl phthalate. Many analyte concentrations were higher among women of Hispanic ethnicity compared to those of non-Hispanic White women. Concentrations of certain chemicals decreased with the calendar year, whereas concentrations of their replacements increased. Our largest study to date identified widespread exposures to prevalent and understudied chemicals in a diverse sample of pregnant women in the U.S.

Chemical Mixtures in Household Environments: In Silico Predictions and In Vitro Testing of Potential Joint Action on PPARγ in Human Liver Cells

There are thousands of chemicals that humans can be exposed to in their everyday environments, the majority of which are currently understudied and lack substantial testing for potential exposure and toxicity. This study aimed to implement in silico methods to characterize the chemicals that co-occur across chemical and product uses in our everyday household environments that also target a common molecular mediator, thus representing understudied mixtures that may exacerbate toxicity in humans. To detail, the Chemical and Products Database (CPDat) was queried to identify which chemicals co-occur across common exposure sources. Chemicals were preselected to include those that target an important mediator of cell health and toxicity, the peroxisome proliferator activated receptor gamma (PPARγ), in liver cells that were identified through query of the ToxCast/Tox21 database. These co-occurring chemicals were thus hypothesized to exert potential joint effects on PPARγ. To test this hypothesis, five commonly co-occurring chemicals (namely, benzyl cinnamate, butyl paraben, decanoic acid, eugenol, and sodium dodecyl sulfate) were tested individually and in combination for changes in the expression of PPARγ and its downstream target, insulin receptor (INSR), in human liver HepG2 cells. Results showed that these likely co-occurring chemicals in household environments increased both PPARγ and INSR expression more significantly when the exposures occurred as mixtures vs. as individual chemicals. Future studies will evaluate such chemical combinations across more doses, allowing for further quantification of the types of joint action while leveraging this method of chemical combination prioritization. This study demonstrates the utility of in silico-based methods to identify chemicals that co-occur in the environment for mixtures toxicity testing and highlights relationships between understudied chemicals and changes in PPARγ-associated signaling.

Application of Pattern Mining Methods to Assess Exposures to Multiple Airborne Chemical Agents in Two Large Occupational Exposure Databases from France

Surveys of the French working population estimate that approximately 15% of all workers may be exposed to at least three different chemical agents, but the most prevalent coexposure situations and their associated health risks remain relatively understudied. To characterize occupational coexposure situations in France, we extracted personal measurement data from COLCHIC and SCOLA, two large administrative occupation exposure databases. We selected 118 chemical agents that had ≥100 measurements with detected concentrations over the period 2010-2019, including 31 carcinogens (IARC groups 1, 2A, and 2B). We grouped measurements by work situations (WS, combination of sector, occupation, task, and year). We characterized the mixtures across WS using frequent itemset mining and association rules mining. The 275,213 measurements extracted came from 32,670 WS and encompassing 4692 unique mixtures. Workers in 32% of all WS were exposed to ≥2 agents (median 3 agents/WS) and 13% of all WS contained ≥2 carcinogens (median 2 carcinogens/WS). The most frequent coexposures were ethylbenzene-xylene (1550 WS), quartz-cristobalite (1417 WS), and toluene-xylene (1305 WS). Prevalent combinations of carcinogens also included hexavalent chromium-lead (368 WS) and benzene-ethylbenzene (314 WS). Wood dust (6% of WS exposed to at least one other agent) and asbestos (8%) had the least amount of WS coexposed with other agents. Tasks with the highest proportions of coexposure to carcinogens include electric arc welding (37% of WS with coexposure), polymerization and distillation (34%), and construction drilling and excavating (34%). Overall, the coexposure to multiple chemical agents, including carcinogens, was highly prevalent in the databases, and should be taken into account when assessing exposure risks in the workplace.

Human Chemical Exposure from Background Emissions in the United States and the Implication for Quantifying Risks from Marginal Emission Increase

The linear dose-response relationship has long been assumed in assessments of health risk from an incremental chemical emission relative to background emissions. In this study, we systematically examine the relevancy of such an assumption with real-world data. We used the reported emission data, as background emissions, from the 2017 U.S. National Emission Inventory for 95 organic chemicals to estimate the central tendencies of exposures of the general U.S. population. Previously published nonlinear dose-response relationships for chemicals were used to estimate health risk from exposure. We also explored and identified four intervals of exposure in which the nonlinear dose-response relationship may be linearly approximated with fixed slopes. Predicted rates of exposure to these 95 chemicals are all within the lowest of the four intervals and associated with low health risk. The health risk may be overestimated if a slope on the dose-response relationship extrapolated from toxicological assays based on high response rates is used for a marginal increase in emission not substantially higher than background emissions. To improve the confidence of human health risk estimates for chemicals, future efforts should focus on deriving a more accurate dose-response relationship at lower response rates and interface it with exposure assessments.