A Chemical Category-Based Prioritization Approach for Selecting 75 Per- and Polyfluoroalkyl Substances (PFAS) for Tiered Toxicity and Toxicokinetic Testing.

Per- and polyfluoroalkyl substances (PFASs) are a group of fluorinated substances of interest to researchers, regulators, and the public due to their widespread presence in the environment. A few PFASs have comparatively extensive amounts of human epidemiological, exposure, and experimental animal toxicity data (e.g., perfluorooctanoic acid), whereas little toxicity and exposure information exists for much of the broader set of PFASs. Given that traditional approaches to generate toxicity information are resource intensive, new approach methods, including in vitro high-throughput toxicity (HTT) testing, are being employed to inform PFAS hazard characterization and further (in vivo) testing. The U.S. Environmental Protection Agency (EPA) and the National Toxicology Program (NTP) are collaborating to develop a risk-based approach for conducting PFAS toxicity testing to facilitate PFAS human health assessments. This article describes the construction of a PFAS screening library and the process by which a targeted subset of 75 PFASs were selected. Multiple factors were considered, including interest to the U.S. EPA, compounds within targeted categories, structural diversity, exposure considerations, procurability and testability, and availability of existing toxicity data. Generating targeted HTT data for PFASs represents a new frontier for informing priority setting. https://doi.org/10.1289/EHP4555.

Introduction

Per- and polyfluoroalkyl substances (PFASs) are a group of fluorinated substances that have generated increased public attention due to their potential health hazard and widespread presence in the environment (Wang et al. 2017; Xiao 2017; Ross et al. 2018). The U.S. Environmental Protection Agency (EPA) Office of Research and Development (ORD) in partnership with the National Toxicology Program (NTP) are currently engaged in producing toxicity information to facilitate human health assessments for PFASs. A few PFASs have comparatively extensive amounts of toxicity data (e.g., perfluorooctanoic acid), but little toxicity information exists for much of the broader set of PFASs identified from preliminary exposure studies that capture potential occurrence in the environment. The hundreds of untested PFASs provide a scenario in which traditional one-by-one toxicity testing would consume tremendous resources and useful toxicity information would not be available for decades. The U.S. EPA’s ToxCast program and the multi-federal agency Tox21 program (which includes the NTP and the U.S. EPA as major partners) have developed the capacity to screen hundreds to thousands of chemicals for bioactivity through in vitro high-throughput toxicity (HTT) testing. Data generated from these assays are already being used to inform hazard identification and prioritize chemicals for further in vivo testing (U.S. EPA 2012, 2014a, 2014b, 2015; Judson et al. 2010, 2015; Kleinstreuer et al. 2017). Within the U.S. EPA, generating such data to inform agency and partner decision making regarding potential human health hazard and risk across the broad landscape of PFASs represents a real-world challenge that HTT coupled with cheminformatic approaches is uniquely designed to address.

This article describes, in brief, the development of the PFAS screening library and the process by which a subset of 75 PFAS substances were selected for HTT screening and tiered toxicity testing, along with mention of the toxicity and toxicokinetic experiments currently underway.

Discussion

Development of the PFAS Screening Library

Since there are no specific chemical catalogs for PFASs, an initial scoping for potentially procurable PFAS substances relied on the use of candidate PFAS structure lists generated from the U.S. EPA’s Distributed Structure–Searchable Toxicity (DSSTox) chemical database. DSSTox currently exceeds 760,000 substances, each of which has undergone some level of chemical structure curation prior to registration (Williams et al. 2017). The largest registered list of PFAS chemicals available at the time this study was initiated was the KEMI PFAS list in DSSTox (named PFASKEMI and available for download at https://comptox.epa.gov/dashboard/chemical_lists/pfaskemi). Approximately 1,200 structures from this list were provided to the chemical contractor for scoping purposes, from which approximately 600 substances were identified as potentially procurable but likely to require on-demand synthesis and exceed standard costs. Based on this preliminary scoping, U.S. EPA funds were secured for the purchase and processing of approximately 400 substances to create a PFAS testing library.

The first procurement phase considered the feasibility of procuring substances of interest to the U.S. EPA. A U.S. EPA workgroup was formed to identify PFASs of interest to U.S. EPA programs and regions and to include PFASs with associated toxicity data that would inform human health risk assessment. The final set of 31 PFASs recommended for further study by this workgroup (list denoted as EPA PFAS WG 31) identified PFASs whose review may support risk evaluation. Also included in the request list for the first phase of PFAS procurements were PFASs that spanned a wider range of U.S. EPA research activities. This larger list of potentially procurable PFASs initially consisted of 89 unique substances (inclusive of EPA PFAS WG 31), which we denote here as EPA-PFAS (note, an updated, expanded U.S. EPA research list, inclusive of salts and anions, is titled EPAPFASRL and is available for download at https://comptox.epa.gov/dashboard/chemical_lists/epapfasrl).

The second phase of PFAS procurements considered a query of the expanded contents of the DSSTox database for chemicals satisfying a range of PFAS-defining criteria, including: fluorines, no aromatic rings (i.e., aliphatic), and molecular weight (MW) . This list was reviewed manually and additional filters were applied to reduce the size of the list from substances to PFAS candidates for possible procurement, for example, excluding heavy metals, halogen salts, low-MW compounds () and compounds for which the ratio of F to C was less than 2:1. The initial set of compounds from this sub-list, for which procurement sources were identified, underwent manual expert DSSTox curation review. The resulting compound list, after confirmation of procurement feasibility, formed the remainder of the initial structure library considered in the present prioritization exercise. This final set of 271 DSSTox-registered substances is referred to herein as the PFAS-Landscape. This set of substances bounded the range of PFASs considered in the below analysis, which was used to identify the candidate 75 subset for tiered toxicity testing.

Categorization of the PFAS Screening Library

Although there are many ways to systematically select a representative subset of structures from a library using different structure-based cheminformatic approaches, in this study predefined, expert-based structural categories were relied upon to characterize the PFAS screening library. The structural categories initially proposed were informed by the work by Buck et al. (2011), who described a systematic terminology for naming and categorizing PFAS substances. The structural categories were used to manually assign each substance in the 271 PFAS library into a respective structural category. To maintain a practical and pragmatic number of structural categories, after this initial assignment was completed, some of the structural categories were combined, for example, n:1 fluorotelomer alcohols and n:2 fluorotelomer alcohols were collapsed into a single category of fluorotelomers. The linkages between the general categories and more specific categories (i.e., subcategories) were retained to offer additional flexibility in the selection of substances for testing. Retaining this layered category information could be particularly useful should specific activity trends within categories be identified and subsequently investigated. Overall, 53 unique structural categories were assigned. Some of the categories contained many more substances (members) than others. Categories containing only one member were referred to as singletons.

Process for Selection of a Subset of 75 PFASs

Using the expert-assigned structural categories described above, we constructed a step-wise workflow to guide efforts to prioritize and weigh various factors for chemical selection within categories. The workflow is graphically illustrated in Figure 1 and consists of an initial PFAS-Landscape characterization step (0), followed by a series of five steps, described in more detail below, to balance the somewhat competing goals of creating a data set that would support read-across within categories while also capturing structural diversity aspects of the PFAS landscape.

Figure 1.

Workflow for selection of structural categories to identify the subset of 75 per- and polyfluoroalkyl substances (PFAS).

Characterizing the PFAS Library (Step 0)

Structural diversity of the full PFAS-Landscape can be represented both in terms of overall chemical counts within categories and labeled by chemical membership in one of the three main groupings, ordered by level of U.S. EPA interest, that is, the 31 PFASs with associated data to inform human health risk assessment (EPA PFAS WG 31), additional PFASs of interest to U.S. EPA researchers (EPA-PFAS), and the remaining PFAS Landscape (PFAS-Landscape) initially identified as procurable. This is illustrated in Figure 2a.

Figure 2.

Per- and polyfluoroalkyl substances (PFAS) Library characteristics by category aligned with workflow steps shown in Figure 1: (a) Total chemical count plot showing the total number of PFASs in a specific structural category color-coded by “status,” that is, membership on a single list (EPA PFAS WG ); (b) ToxVal record count plot showing the number of PFASs with in vivo toxicity studies available within the category, as sourced from the U.S. EPA ToxVal database; (c) Same count plot as in a, showing only the final selected 75 PFAS substances by list membership and category representation. Color figures are available at https://doi.org/10.23645/epacomptox.7479866.

The U.S. EPA’s ToxVal database is a database of source-referenced human health reference or toxicity values collected from in vivo studies that are available through the U.S. EPA CompTox Chemicals Dashboard (https://comptox.epa.gov/dashboard). Figure 2b provides a representation of record(s) from ToxVal as surrogates for in vivo toxicity information per chemical (1, yes; 0, no). These records are mapped onto the same category legend as shown in Figure 2a and demonstrate that toxicity data are available for only a limited subset of categories of interest to the U.S. EPA.

Selection of Structural Categories (Steps 1–5)

Strategies for selecting a subset of PFAS substances for HTT from the PFAS-Landscape focused on two main objectives: 1) maximizing information to support read-across within structure-based groupings, and 2) capturing the structural diversity of the PFAS landscape of interest to the U.S. EPA. Steps 1–3 of the workflow in Figure 1 identified structural categories addressing Objective 1 (maximizing read-across), whereas Steps 4–5 identified structural categories addressing Objective 2 (capturing structural diversity).

In Step 1 of the workflow, structural categories were identified that were both of high interest to the U.S. EPA (EPA PFAS WG 31) and had a record in ToxVal indicating availability of in vivo data. Step 2 identified structural categories that were of broader interest to the U.S. EPA (EPA-PFAS and PFAS-Landscape) and had a record in ToxVal. Structural categories that had a record in ToxVal but were not captured in Steps 1–2, that is, of lesser interest to the agency, were identified in Step 3. Step 4 identified structural categories of interest to the U.S. EPA irrespective of having a record in ToxVal. The remaining structural categories in PFAS-Landscape not captured in Steps 1–4 were considered as part of Step 5.

Substances to address read-across (Objective 1) were drawn from the 10 categories identified in Steps 1–3. Substances for capturing structural diversity were drawn from the remaining 43 categories.

Selection of PFASs within the Prioritized Structural Categories

Test substances were selected on a structural category basis. The process was guided by a quantitative scoring scheme that aimed to capture and rank considerations such as structural variation and physical property information (logKow, vapor pressure), as well as availability of a record in the ToxCast and/or ToxVal databases. The scheme comprised seven aspects described in Table 1.

Table 1

Considerations for selection of the 75 PFAS.

Aspect name	Scoring
1) Structural diversity within a category	Approximated by category size, with score ranging from 1 (≥20 members) to 0 (1 member).
2) Data availability	Availability of in vitro ToxCast data (score=0.5) or ToxVal in vivo data (score=0.75) or both (score=1).
3) Data quantity	Number of ToxVal records for a substance indicating a stronger source–analog for read-across, with scores ranging from 0.15 (for 1 record) to 1 (for ≥20 records).
4) Read-across category-level weight	Value of substance for anchoring read-across trends within a category (e.g., chain length), serving as a source analog (score=0.5) or target analog (score=0.25), or as a target analog for capturing structural diversity (score=0.15).
5) Numerical indicator of U.S. EPA interest	EPA PFAS WG31 (score=1), other U.S. EPA-PFAS (score=0.75), only in PFAS-Landscape (score=0.5).
6) Physicochemical indicators of testability	Both logKow and vapor pressure favorable (score=0.75), one favorable (score=0.5), both unfavorable (score=0). For example, logKow <4.5, vapor pressure <10^3mmHg considered favorable. LogKow and vapor pressure properties relied on OPERA model predictions as available from the U.S. EPA CompTox Chemicals Dashboard. It is recognized that there are issues surrounding the validity of predictions for PFAS substances. The predictions here were used in relative terms within a structural category.
7) Figure 1 workflow step	Step 1 (score=1), Step 2 (score=0.75), Step 3 (score=0.5), Step 4 (score=0.25), Step 5 (score=0).
Total score	Summation of scores from the preceding considerations used to rank each PFAS substance.

Note: OPERA, OPEn structure-activity/property Relationship App; PFAS, per- and polyfluoroalkyl substances.

An initial selection of 75 substances was made based on these scoring considerations. Backup alternatives were also selected based on the same scoring considerations. However, technical challenges were encountered as procurement orders were processed. Of specific note was the physical form of the test substance received (i.e., gas vs. solid), hazmat considerations (e.g., flammability), evidence for volatility/sublimation of stored neat samples, and insolubility in dimethyl sulfoxide (DMSO). These technical considerations resulted in further adjustment of the selection of substances from specific categories, which in turn impacted the degree of structural diversity reflected in the final PFAS procured library as well as the extent to which categories for maximizing read-across were represented. The final set of 75 substances for which DMSO solutions were prepared, and test samples submitted for HTT, is represented graphically by category and status assignment in Figure 2c. The final DSSTox list of 75 PFAS substances, with associated structural information, is labeled EPAPFAS75S1 and is available for download at https://comptox.epa.gov/dashboard/chemical_lists/epapfas75S1.

HTT: Tiered Toxicity and Toxicokinetic Testing

The final set of 75 PFAS substances is currently undergoing targeted and tiered HTT in partnership with NTP. Tier 1 HTT includes in vitro assays focused on multiple end points, including hepatotoxicity, immunotoxicity, developmental toxicity, mitochondrial toxicity, and developmental neurotoxicity along with assays to estimate in vivo toxicokinetics. The assays selected for the Tier 1 toxicity and toxicokinetic characterization were based on both the known in vivo adverse responses of previously tested PFASs and the anticipated effects of a broader range of PFAS. In general, the proposed strategy is to utilize data generated from new approach methods (e.g., HTT toxicological and toxicokinetic assays) in combination with human exposure information (measured and/or predicted) to derive a biological exposure ratio (BER; Thomas et al. 2013). BERs will serve as a measure of potential risk and will be used to prioritize subsequent Tier 2 in vivo testing and inform human health risk assessment. Data generated from in vitro assays will also be used with existing in vivo information to support read-across efforts (Patlewicz et al. 2013).

Summary

A PFAS screening library was constructed and categorized by expert review into 53 structural categories. The final PFAS 75 (EPAPFAS75S1) comprised 46 substances representing 10 of the structural categories with some existing in vivo toxicity information and 29 substances covering a further 24 structural categories (Figure 2c). New PFAS candidates will be selected from the U.S. EPA’s complete, DMSO-solubilized PFAS inventory (EPAPFASINV), now totaling 430 unique DSSTox substances. This complete inventory and the list of 43 PFASs that were procured but found to be DMSO-insoluble (EPAPFASINSOL) as well as the final EPAPFAS75S1 list are available for download at https://comptox.epa.gov/dashboard/chemical_lists/pfasmaster.