tcplfit2: an R-language general purpose concentration-response modeling package.

SUMMARY: Many applications of chemical screening are performed in concentration or dose-response mode, and it is necessary to extract appropriate parameters, including whether the chemical/assay pair is active and if so, what are concentrations where activity is seen. Typically, multiple mathematical models or curve shapes are tested against the data to assess the best fit. There are several commercial programs used for this purpose as well as open-source libraries. A widely used system for managing high-throughput screening (HTS) concentration-response data is tcpl (ToxCast Pipeline). The current implementation of tcpl has the concentration-response modeling code tightly integrated with the data management and databasing aspects of HTS data processing. Tcplfit2 is a stand-alone version of the curve-fitting and hitcalling core of tcpl that has been extended to include a large number of standard curve classes and to use benchmark dose modeling. This package will be useful for HTS concentration-response data such as high-throughput whole genome transcriptomics. AVAILABILITY AND IMPLEMENTATION: tcplfit2 is written in R and is available from CRAN.

Application of Cell Painting for chemical hazard evaluation in support of screening-level chemical assessments.

'Cell Painting' is an imaging-based high-throughput phenotypic profiling (HTPP) method in which cultured cells are fluorescently labeled to visualize subcellular structures (i.e., nucleus, nucleoli, endoplasmic reticulum, cytoskeleton, Golgi apparatus / plasma membrane and mitochondria) and to quantify morphological changes in response to chemicals or other perturbagens. HTPP is a high-throughput and cost-effective bioactivity screening method that detects effects associated with many different molecular mechanisms in an untargeted manner, enabling rapid in vitro hazard assessment for thousands of chemicals. Here, 1201 chemicals from the ToxCast library were screened in concentration-response up to ∼100 µM in human U-2 OS cells using HTPP. A phenotype altering concentration (PAC) was estimated for chemicals active in the tested range. PACs tended to be higher than lower bound potency values estimated from a broad collection of targeted high-throughput assays, but lower than the threshold for cytotoxicity. In vitro to in vivo extrapolation (IVIVE) was used to estimate administered equivalent doses (AEDs) based on PACs for comparison to human exposure predictions. AEDs for 18/412 chemicals overlapped with predicted human exposures. Phenotypic profile information was also leveraged to identify putative mechanisms of action and group chemicals. Of 58 known nuclear receptor modulators, only glucocorticoids and retinoids produced characteristic profiles; and both receptor types are expressed in U-2 OS cells. Thirteen chemicals with profile similarity to glucocorticoids were tested in a secondary screen and one chemical, pyrene, was confirmed by an orthogonal gene expression assay as a novel putative GR modulating chemical. Most active chemicals demonstrated profiles not associated with a known mechanism-of-action. However, many structurally related chemicals produced similar profiles, with exceptions such as diniconazole, whose profile differed from other active conazoles. Overall, the present study demonstrates how HTPP can be applied in screening-level chemical assessments through a series of examples and brief case studies.

Assessing utility of thyroid in vitro screening assays through comparisons to observed impacts in vivo.

To better understand endocrine disruption, the U.S. Environmental Protection Agency's (USEPA) Endocrine Disruptor Screening Program (EDSP) utilizes a two-tiered approach to investigate the potential of a chemical to interact with the estrogen, androgen, or thyroid systems. As in vivo testing lacks the throughput to address data gaps on endocrine bioactivity for thousands of chemicals, in vitro high-throughput screening (HTS) methods are being developed to screen larger chemical libraries. The primary objective of this work was to investigate for how many of the 52 chemicals with weight-of-evidence (WoE) determinations from EDSP Tier 1 screening there are available in vitro HTS data supporting a thyroid impact. HTS data from the USEPA ToxCast program and the EDSP WoE were collected for this analysis. Considering the complexity of endocrine disruption and interpreting HTS data, concordance between in vitro activity and in vivo effects ranges from 58 to 78%. Based on this evaluation, we conclude that the current suite of HTS assays is beneficial for prioritizing chemicals for further inquiry; however, without a more detailed analysis, one cannot conclude whether HTS results are the primary mode-of-action. Furthermore, development of in vitro assays for additional thyroid-relevant molecular initiating events is required to effectively predict in vivo thyroid impacts.

Integration of toxicodynamic and toxicokinetic new approach methods into a weight-of-evidence analysis for pesticide developmental neurotoxicity assessment: A case-study with DL- and L-glufosinate.

DL-glufosinate ammonium (DL-GLF) is a registered herbicide for which a guideline Developmental Neurotoxicity (DNT) study has been conducted. Offspring effects included altered brain morphometrics, decreased body weight, and increased motor activity. Guideline DNT studies are not available for its enriched isomers L-GLF acid and L-GLF ammonium; conducting one would be time consuming, resource-intensive, and possibly redundant given the existing DL-GLF DNT. To support deciding whether to request a guideline DNT study for the L-GLF isomers, DL-GLF and the L-GLF isomers were screened using in vitro assays for network formation and neurite outgrowth. DL-GLF and L-GLF isomers were without effects in both assays. DL-GLF and L-GLF (1-100 µM) isomers increased mean firing rate of mature networks to 120-140% of baseline. In vitro toxicokinetic assessments were used to derive administered equivalent doses (AEDs) for the in vitro testing concentrations. The AED for L-GLF was ∼3X higher than the NOAEL from the DL-GLF DNT indicating that the available guideline study would be protective of potential DNT due to L-GLF exposure. Based in part on the results of these in vitro studies, EPA is not requiring L-GLF isomer guideline DNT studies, thereby providing a case study for a useful application of DNT screening assays.

An exposure:activity profiling method for interpreting high-throughput screening data for estrogenic activity--proof of concept.

Rapid high throughput in vitro screening (HTS) assays are now available for characterizing dose-responses in assays that have been selected for their sensitivity in detecting estrogen-related endpoints. For example, EPA's ToxCast™ program recently released endocrine assay results for more than 1800 substances and the interagency Tox21 consortium is in the process of releasing data for approximately 10,000 chemicals. But such activity measurements alone fall short for the purposes of priority setting or screening because the relevant exposure context is not considered. Here, we extend the method of exposure:activity profiling by calculating the exposure:activity ratios (EARs) using human exposure estimates and AC50 values for a range of chemicals tested in a suite of seven estrogenic assays in ToxCast™ and Tox21. To provide additional context, relative estrogenic exposure:activity quotients (REEAQ) were derived by comparing chemical-specific EARs to the EAR of the ubiquitous dietary phytoestrogen, genistein (GEN). Although the activity of a substance in HTS-endocrine assays is not a measure of health hazard or risk, understanding how such a dose compares to human exposures provides a valuable additional metric that can be used in decision-making; substances with small EARs and REEAQs would indicate low priority for further endocrine screening or testing.

A high throughput screening assay for human Thyroperoxidase inhibitors.

Rapid, human relevant assays are needed to assess potential hazards of the many chemicals in commerce. An assay of thyroid peroxidase (TPO) inhibition, using the substrate Amplex Ultra Red, was recently adapted for human TPO (AUR-hTPO). We tested a large number (788) of chemicals through this AUR-hTPO assay and compared performance with published results from an assay using enzyme from rat thyroid microsomes (AUR-rTPO). Coded chemicals, from the US EPA ToxCast Inventory, were tested in a tiered approach: 1) Initial screening at a single concentration; 2) Potency estimation for active chemicals with multiple concentrations; 3) Screening active chemicals for the non-specific activity. The assay gave consistent results for positive chemical methimazole and several positive and negative reference chemicals. hTPO inhibition was observed for 190 chemicals reported as positive in rTPO. Of these, 158 showed no confounding activity (interference due to fluorescence or non-specific protein inhibition). Comparison of all result with rTPO data and with evidence of TPO inhibition found in the literature suggest that the current assay has a higher rate of false negative but a much lower rate of false positive compared with the rTPO screen. These findings underscore the effectiveness of the AUR assay, using hTPO enzyme from engineered cell lines, to identify moderate to strong inhibitors but some improvements may be needed to detect weak TPO inhibitors.

Exploring the effects of experimental parameters and data modeling approaches on in vitro transcriptomic point-of-departure estimates.

Multiple new approach methods (NAMs) are being developed to rapidly screen large numbers of chemicals to aid in hazard evaluation and risk assessments. High-throughput transcriptomics (HTTr) in human cell lines has been proposed as a first-tier screening approach for determining the types of bioactivity a chemical can cause (activation of specific targets vs. generalized cell stress) and for calculating transcriptional points of departure (tPODs) based on changes in gene expression. In the present study, we examine a range of computational methods to calculate tPODs from HTTr data, using six data sets in which MCF7 cells cultured in two different media formulations were treated with a panel of 44 chemicals for 3 different exposure durations (6, 12, 24 hr). The tPOD calculation methods use data at the level of individual genes and gene set signatures, and compare data processed using the ToxCast Pipeline 2 (tcplfit2), BMDExpress and PLIER (Pathway Level Information ExtractoR). Methods were evaluated by comparing to in vitro PODs from a validated set of high-throughput screening (HTS) assays for a set of estrogenic compounds. Key findings include: (1) for a given chemical and set of experimental conditions, tPODs calculated by different methods can vary by several orders of magnitude; (2) tPODs are at least as sensitive to computational methods as to experimental conditions; (3) in comparison to an external reference set of PODs, some methods give generally higher values, principally PLIER and BMDExpress; and (4) the tPODs from HTTr in this one cell type are mostly higher than the overall PODs from a broad battery of targeted in vitro ToxCast assays, reflecting the need to test chemicals in multiple cell types and readout technologies for in vitro hazard screening.

Evaluation of a high-throughput H295R homogenous time resolved fluorescence assay for androgen and estrogen steroidogenesis screening.

The H295R test guideline assay evaluates the effect of test substances on synthesis of 17ß-estradiol (E2) and testosterone (T). The objective of this study was to leverage commercial immunoassay technology to develop a more efficient H295R assay to measure E2 and T levels in 384-well format. The resulting Homogenous Time Resolved Fluorescence assay platform (H295R-HTRF) was evaluated against a training set of 36 chemicals derived from the OECD inter-laboratory validation study, EPA guideline 890.1200 aromatase assay, and azole fungicides active in the HT-H295R assay. Quality control performance criteria were met for all conditions except E2 synthesis inhibition where low basal hormone synthesis was observed. Five proficiency chemicals were active for both the E2 and T endpoints, consistent with guideline classifications. Of the nine OECD core reference chemicals, 9/9 were concordant with outcomes for E2 and 7/9 for T. Likewise, 9/13 and 11/13 OECD supplemental chemicals were concordant with anticipated effects for E2 and T, respectively. Of the 10 azole fungicides screened, 7/10 for E2 and 8/10 for T exhibited concordant outcomes for inhibition. Generally, all active chemicals in the training set demonstrated equivalent or greater potency in the H295R-HTRF assay, supporting the sensitivity of the platform. The adaptation of HTRF technology to the H295R model provides an efficient way to evaluate E2 and T modulators in accordance with guideline specifications.

Reproducibility of organ-level effects in repeat dose animal studies.

This work estimates benchmarks for new approach method (NAM) performance in predicting organ-level effects in repeat dose studies of adult animals based on variability in replicate animal studies. Treatment-related effect values from the Toxicity Reference database (v2.1) for weight, gross, or histopathological changes in the adrenal gland, liver, kidney, spleen, stomach, and thyroid were used. Rates of chemical concordance among organ-level findings in replicate studies, defined by repeated chemical only, chemical and species, or chemical and study type, were calculated. Concordance was 39 - 88%, depending on organ, and was highest within species. Variance in treatment-related effect values, including lowest effect level (LEL) values and benchmark dose (BMD) values when available, was calculated by organ. Multilinear regression modeling, using study descriptors of organ-level effect values as covariates, was used to estimate total variance, mean square error (MSE), and root residual mean square error (RMSE). MSE values, interpreted as estimates of unexplained variance, suggest study descriptors accounted for 52-69% of total variance in organ-level LELs. RMSE ranged from 0.41 - 0.68 log10-mg/kg/day. Differences between organ-level effects from chronic (CHR) and subchronic (SUB) dosing regimens were also quantified. Odds ratios indicated CHR organ effects were unlikely if the SUB study was negative. Mean differences of CHR - SUB organ-level LELs ranged from -0.38 to -0.19 log10 mg/kg/day; the magnitudes of these mean differences were less than RMSE for replicate studies. Finally, in vitro to in vivo extrapolation (IVIVE) was employed to compare bioactive concentrations from in vitro NAMs for kidney and liver to LELs. The observed mean difference between LELs and mean IVIVE dose predictions approached 0.5 log10-mg/kg/day, but differences by chemical ranged widely. Overall, variability in repeat dose organ-level effects suggests expectations for quantitative accuracy of NAM prediction of LELs should be at least ± 1 log10-mg/kg/day, with qualitative accuracy not exceeding 70%.

Evaluation of 147 perfluoroalkyl substances for immunotoxic and other (patho)physiological activities through phenotypic screening of human primary cells.

A structurally diverse set of 147 per- and polyfluoroalkyl substances (PFAS) was screened in a panel of 12 human primary cell systems by measuring 148 biomarkers relevant to (patho)physiological pathways to inform hypotheses about potential mechanistic effects of data-poor PFAS in human model systems. This analysis focused on immunosuppressive activity, which was previously reported as an in vivo effect of perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS), by comparing PFAS responses to four pharmacological immunosuppressants. The PFOS response profile had little correlation with reference immunosuppressants, suggesting in vivo activity does not occur by similar mechanisms. The PFOA response profile did share features with the profile of dexamethasone, although some distinct features were lacking. Other PFAS, including 2,2,3,3-tetrafluoropropyl acrylate, demonstrated more similarity to the reference immunosuppressants but with additional activities not found in the reference immunosuppressive drugs. Correlation of PFAS profiles with a database of environmental chemical responses and pharmacological probes identified potential mechanisms of bioactivity for some PFAS, including responses similar to ubiquitin ligase inhibitors, deubiquitylating enzyme (DUB) inhibitors, and thioredoxin reductase inhibitors. Approximately 21% of the 147 PFAS with confirmed sample quality were bioactive at nominal testing concentrations in the 1-60 micromolar range in these human primary cell systems. These data provide new hypotheses for mechanisms of action for a subset of PFAS and may further aid in development of a PFAS categorization strategy useful in safety assessment.

Environmental mixtures and breast cancer: identifying co-exposure patterns between understudied vs breast cancer-associated chemicals using chemical inventory informatics.

BACKGROUND: Although evidence linking environmental chemicals to breast cancer is growing, mixtures-based exposure evaluations are lacking. OBJECTIVE: This study aimed to identify environmental chemicals in use inventories that co-occur and share properties with chemicals that have association with breast cancer, highlighting exposure combinations that may alter disease risk. METHODS: The occurrence of chemicals within chemical use categories was characterized using the Chemical and Products Database. Co-exposure patterns were evaluated for chemicals that have an association with breast cancer (BC), no known association (NBC), and understudied chemicals (UC) identified through query of the Silent Spring Institute's Mammary Carcinogens Review Database and the U.S. Environmental Protection Agency's Toxicity Reference Database. UCs were ranked based on structure and physicochemical similarities and co-occurrence patterns with BCs within environmentally relevant exposure sources. RESULTS: A total of 6793 chemicals had data available for exposure source occurrence analyses. 50 top-ranking UCs spanning five clusters of co-occurring chemicals were prioritized, based on shared properties with co-occuring BCs, including chemicals used in food production and consumer/personal care products, as well as potential endocrine system modulators. SIGNIFICANCE: Results highlight important co-exposure conditions that are likely prevalent within our everyday environments that warrant further evaluation for possible breast cancer risk. IMPACT STATEMENT: Most environmental studies on breast cancer have focused on evaluating relationships between individual, well-known chemicals and breast cancer risk. This study set out to expand this research field by identifying understudied chemicals and mixtures that may occur in everyday environments due to their patterns of commercial use. Analyses focused on those that co-occur alongside chemicals associated with breast cancer, based upon in silico chemical database querying and analysis. Particularly in instances when understudied chemicals share physicochemical properties and structural features with carcinogens, these chemical mixtures represent conditions that should be studied in future clinical, epidemiological, and toxicological studies.

An expert-driven literature review of "negative" chemicals for developmental neurotoxicity (DNT) in vitro assay evaluation.

To date, approximately 200 chemicals have been tested in US Environmental Protection Agency (EPA) or Organization for Economic Co-operation and Development (OECD) developmental neurotoxicity (DNT) guideline studies, leaving thousands of chemicals without traditional animal information on DNT hazard potential. To address this data gap, a battery of in vitro DNT new approach methodologies (NAMs) has been proposed. Evaluation of the performance of this battery will increase the confidence in its use to determine DNT chemical hazards. One approach to evaluate DNT NAM performance is to use a set of chemicals to evaluate sensitivity and specificity. Since a list of chemicals with potential evidence of in vivo DNT has been established, this study aims to develop a curated list of "negative" chemicals for inclusion in a "DNT NAM evaluation set". A workflow, including a literature search followed by an expert-driven literature review, was used to systematically screen 39 chemicals for lack of DNT effect. Expert panel members evaluated the scientific robustness of relevant studies to inform chemical categorizations. Following review, the panel discussed each chemical and made categorical determinations of "Favorable", "Not Favorable", or "Indeterminate" reflecting acceptance, lack of suitability, or uncertainty given specific limitations and considerations, respectively. The panel determined that 10, 22, and 7 chemicals met the criteria for "Favorable", "Not Favorable", and "Indeterminate", for use as negatives in a DNT NAM evaluation set. Ultimately, this approach not only supports DNT NAM performance evaluation but also highlights challenges in identifying large numbers of negative DNT chemicals.

Structure-based QSAR Models to Predict Repeat Dose Toxicity Points of Departure.

Human health risk assessment for environmental chemical exposure is limited by a vast majority of chemicals with little or no experimental in vivo toxicity data. Data gap filling techniques, such as quantitative structure activity relationship (QSAR) models based on chemical structure information, can predict hazard in the absence of experimental data. Risk assessment requires identification of a quantitative point-of-departure (POD) value, the point on the dose-response curve that marks the beginning of a low-dose extrapolation. This study presents two sets of QSAR models to predict POD values (PODQSAR) for repeat dose toxicity. For training and validation, a publicly available in vivo toxicity dataset for 3592 chemicals was compiled using the U.S. Environmental Protection Agency's Toxicity Value database (ToxValDB). The first set of QSAR models predict point-estimates of POD values (PODQSAR) using structural and physicochemical descriptors for repeat dose study types and species combinations. A random forest QSAR model using study type and species as descriptors showed the best performance, with an external test set root mean square error (RMSE) of 0.71 log10-mg/kg/day and coefficient of determination (R2) of 0.53. The second set of QSAR models predict the 95% confidence intervals for PODQSAR using a constructed POD distribution with a mean equal to the median POD value and a standard deviation of 0.5 log10-mg/kg/day, based on previously published typical study-to-study variability that may lead to uncertainty in model predictions. Bootstrap resampling of the pre-generated POD distribution was used to derive point-estimates and 95% confidence intervals for each POD prediction. Enrichment analysis to evaluate the accuracy of PODQSAR showed that 80% of the 5% most potent chemicals were found in the top 20% of the most potent chemical predictions, suggesting that the repeat dose POD QSAR models presented here may help inform screening level human health risk assessments in the absence of other data.

Python BMDS: A Python interface library and web application for the canonical EPA dose-response modeling software.

Several primary sources of publicly available, quantitative dose-response data from traditional toxicology study designs relevant to predictive toxicology applications are now available, including the redeveloped U.S. Environmental Protection Agency's Toxicity Reference Database (ToxRefDB v2.0), the Health Assessment Workspace Collaborative (HAWC), and the National Toxicology Program's Chemical Program's Chemical Effects in Biological Systems (CEBS). These resources provide effect level information but modeling these data to a curve may be more informative for predictive toxicology applications. Benchmark Dose Software (BMDS) has been recognized broadly and used for regulatory applications at multiple agencies. However, the current BMDS software was not amenable to modeling large datasets. Herein we describe development and use of a Python package that implements a wrapper around BMDS, a software that requires manual input in the dose-response modeling process (i.e., best-fitting model-selection, reporting, and dose-dropping). In the Python BMDS, users can select the BMDS version, customize model recommendation logic, and export summaries of the resultant BMDS output. Further, using the Python interface, a web-based application programming interface (API) has been developed for easy integration into other software systems, pipelines, or databases. Software utility was demonstrated via modeling nearly 28,000 datasets in ToxRefDB v2.0, re-creation of an existing, published large-scale analysis, and demonstration of usage in software such as CEBS and HAWC. Python BMDS enables rapid-batch processing of dose-response datasets using a modeling software with broad acceptance in the toxicology community, thereby providing an important tool for leveraging the publicly available quantitative toxicology data in a reproducible manner.

Evaluation of Chemical Effects on Network Formation in Cortical Neurons Grown on Microelectrode Arrays.

Thousands of chemicals to which humans are potentially exposed have not been evaluated for potential developmental neurotoxicity (DNT), driving efforts to develop a battery of in vitro screening approaches for DNT hazard. Here, 136 unique chemicals were evaluated for potential DNT hazard using a network formation assay (NFA) in cortical cells grown on microelectrode arrays. The effects of chemical exposure from 2 h postplating through 12 days in vitro (DIV) on network formation were evaluated at DIV 5, 7, 9, and 12, with cell viability assessed at DIV 12. Only 82 chemicals altered at least 1 network development parameter. Assay results were reproducible; 10 chemicals tested as biological replicates yielded qualitative results that were 100% concordant, with consistent potency values. Toxicological tipping points were determined for 58 chemicals and were similar to or lower than the lowest 50% effect concentrations (EC50) for all parameters. When EC50 and tipping point values from the NFA were compared to the range of potencies observed in ToxCast assays, the NFA EC50 values were less than the lower quartile for ToxCast assay potencies for a subset of chemicals, many of which are acutely neurotoxic in vivo. For 13 chemicals with available in vivo DNT data, estimated administered equivalent doses based on NFA results were similar to or lower than administered doses in vivo. Collectively, these results indicate that the NFA is sensitive to chemicals acting on nervous system function and will be a valuable contribution to an in vitro DNT screening battery.

ToxRefDB version 2.0: Improved utility for predictive and retrospective toxicology analyses.

The Toxicity Reference Database (ToxRefDB) structures information from over 5000 in vivo toxicity studies, conducted largely to guidelines or specifications from the US Environmental Protection Agency and the National Toxicology Program, into a public resource for training and validation of predictive models. Herein, ToxRefDB version 2.0 (ToxRefDBv2) development is described. Endpoints were annotated (e.g. required, not required) according to guidelines for subacute, subchronic, chronic, developmental, and multigenerational reproductive designs, distinguishing negative responses from untested. Quantitative data were extracted, and dose-response modeling for nearly 28,000 datasets from nearly 400 endpoints using Benchmark Dose (BMD) Modeling Software were generated and stored. Implementation of controlled vocabulary improved data quality; standardization to guideline requirements and cross-referencing with United Medical Language System (UMLS) connects ToxRefDBv2 observations to vocabularies linked to UMLS, including PubMed medical subject headings. ToxRefDBv2 allows for increased connections to other resources and has greatly enhanced quantitative and qualitative utility for predictive toxicology.

Evaluating Chemicals for Thyroid Disruption: Opportunities and Challenges with in Vitro Testing and Adverse Outcome Pathway Approaches.

BACKGROUND: Extensive clinical and experimental research documents the potential for chemical disruption of thyroid hormone (TH) signaling through multiple molecular targets. Perturbation of TH signaling can lead to abnormal brain development, cognitive impairments, and other adverse outcomes in humans and wildlife. To increase chemical safety screening efficiency and reduce vertebrate animal testing, in vitro assays that identify chemical interactions with molecular targets of the thyroid system have been developed and implemented. OBJECTIVES: We present an adverse outcome pathway (AOP) network to link data derived from in vitro assays that measure chemical interactions with thyroid molecular targets to downstream events and adverse outcomes traditionally derived from in vivo testing. We examine the role of new in vitro technologies, in the context of the AOP network, in facilitating consideration of several important regulatory and biological challenges in characterizing chemicals that exert effects through a thyroid mechanism. DISCUSSION: There is a substantial body of knowledge describing chemical effects on molecular and physiological regulation of TH signaling and associated adverse outcomes. Until recently, few alternative nonanimal assays were available to interrogate chemical effects on TH signaling. With the development of these new tools, screening large libraries of chemicals for interactions with molecular targets of the thyroid is now possible. Measuring early chemical interactions with targets in the thyroid pathway provides a means of linking adverse outcomes, which may be influenced by many biological processes, to a thyroid mechanism. However, the use of in vitro assays beyond chemical screening is complicated by continuing limits in our knowledge of TH signaling in important life stages and tissues, such as during fetal brain development. Nonetheless, the thyroid AOP network provides an ideal tool for defining causal linkages of a chemical exerting thyroid-dependent effects and identifying research needs to quantify these effects in support of regulatory decision making. https://doi.org/10.1289/EHP5297.

Novel application of normalized pointwise mutual information (NPMI) to mine biomedical literature for gene sets associated with disease: use case in breast carcinogenesis.

Advances in technology within biomedical sciences have led to an inundation of data across many fields, raising new challenges in how best to integrate and analyze these resources. For example, rapid chemical screening programs like the US Environmental Protection Agency's ToxCast and the collaborative effort, Tox21, have produced massive amounts of information on putative chemical mechanisms where assay targets are identified as genes; however, systematically linking these hypothesized mechanisms with in vivo toxicity endpoints like disease outcomes remains problematic. Herein we present a novel use of normalized pointwise mutual information (NPMI) to mine biomedical literature for gene associations with biological concepts as represented by Medical Subject Headings (MeSH terms) in PubMed. Resources that tag genes to articles were integrated, then cross-species orthologs were identified using UniRef50 clusters. MeSH term frequency was normalized to reflect the MeSH tree structure, and then the resulting GeneID-MeSH associations were ranked using NPMI. The resulting network, called Entity MeSH Co-occurrence Network (EMCON), is a scalable resource for the identification and ranking of genes for a given topic of interest. The utility of EMCON was evaluated with the use case of breast carcinogenesis. Topics relevant to breast carcinogenesis were used to query EMCON and retrieve genes important to each topic. A breast cancer gene set was compiled through expert literature review (ELR) to assess performance of the search results. We found that the results from EMCON ranked the breast cancer genes from ELR higher than randomly selected genes with a recall of 0.98. Precision of the top five genes for selected topics was calculated as 0.87. This work demonstrates that EMCON can be used to link in vitro results to possible biological outcomes, thus aiding in generation of testable hypotheses for furthering understanding of biological function and the contribution of chemical exposures to disease.