High-Throughput Screening Data Interpretation in the Context of In Vivo Transcriptomic Responses to Oral Cr(VI) Exposure

Extracellular Vesicles altered by a Per- and Polyfluoroalkyl Substance Mixture: In Vitro Dose-Dependent Release, Chemical Content, and MicroRNA Signatures involved in Liver Health

Per- and polyfluoroalkyl substances (PFAS) have emerged as high priority contaminants due to their ubiquity and pervasiveness in the environment. Numerous PFAS co-occur across sources of drinking water, including areas of North Carolina (NC) with some detected concentrations above the Environmental Protection Agency's health advisory levels. While evidence demonstrates PFAS exposure induces harmful effects in the liver, the involvement of extracellular vesicles (EVs) as potential mediators of these effects has yet to be evaluated. This study set out to evaluate the hypothesis that PFAS mixtures induce dose-dependent release of EVs from liver cells, with exposures causing differential loading of microRNAs (miRNAs) and PFAS chemical signatures. To test this hypothesis, a defined PFAS mixture was prioritized utilizing data collected by the NC PFAS Testing Network. This mixture contained three substances, PFOS, PFOA, and PFHxA, selected based upon co-occurrence patterns and the inclusion of both short-chain (PFHxA) and long-chain (PFOA and PFOS) substances. HepG2 liver cells were exposed to equimolar PFAS, and secreted EVs were isolated from conditioned media and characterized for count and molecular content. Exposures induced a dose-dependent release of EVs carrying miRNAs that were differentially loaded upon exposure. These altered miRNA signatures were predicted to target mRNA pathways involved in hepatic fibrosis and cancer. Chemical concentrations of PFOS, PFOA, and PFHxA were also detected in both parent HepG2 cells and their released EVs, specifically within a 15-fold range after normalizing for protein content. This study therefore established EVs as novel biological responders and measurable endpoints for evaluating PFAS-induced toxicity.

The placenta epigenome-brain axis: placental epigenomic and transcriptomic responses that preprogram cognitive impairment

Aim: The placenta-brain axis reflects a developmental linkage where disrupted placental function is associated with impaired neurodevelopment later in life. Placental gene expression and the expression of epigenetic modifiers such as miRNAs may be tied to these impairments and are understudied. Materials & methods: The expression levels of mRNAs (n = 37,268) and their targeting miRNAs (n = 2083) were assessed within placentas collected from the ELGAN study cohort (n = 386). The ELGAN adolescents were assessed for neurocognitive function at age 10 and the association with placental mRNA/miRNAs was determined. Results: Placental mRNAs related to inflammatory and apoptotic processes are under miRNA control and associated with cognitive impairment at age 10. Conclusion: Findings highlight key placenta epigenome-brain relationships that support the developmental origins of health and disease hypothesis.

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.

Crypt and Villus Transcriptomic Responses in Mouse Small Intestine Following Oral Exposure to Hexavalent Chromium

Oral exposure to hexavalent chromium (Cr(VI)) induces tumors in the mouse duodenum. Previous microarray-based transcriptomic analyses of homogenized mouse duodenal tissue have demonstrated Cr(VI)-induced alterations in various cellular pathways and processes. However, X-ray fluorescence microscopy indicates that chromium localizes primarily to the duodenal villi following exposure to Cr(VI), suggesting that previous transcriptomic analyses of homogenized tissue provide an incomplete picture of transcriptomic responses in the duodenum. Herein, transcriptomic analyses were conducted separately on crypt and villus tissue from formalin-fixed paraffin-embedded transverse duodenal sections from the same study in which microarray-based analyses were previously conducted. A total of 28 groups (7 doses × 2 timepoints × 2 tissue compartments) were analyzed for differential gene expression, dose-response, and gene set enrichment. Tissue compartment isolation was confirmed by differences in expression of typical markers of crypt and villus compartments. Fewer than 21 genes were altered in the crypt compartment of mice exposed to 0.1-5 ppm Cr(VI) for 7 or 90 days, which increased to hundreds or thousands of genes at ≥20 ppm Cr(VI). Consistent with histological evidence for crypt proliferation, a significant, dose-dependent increase in genes that regulate mitotic cell cycle was prominent in the crypt, while subtle in the villus, when compared with samples from time-matched controls. Minimal transcriptomic evidence of DNA damage response in either the crypts or the villi is consistent with published in vivo genotoxicity data. These results are also discussed in the context of modes of action that have been proposed for Cr(VI)-induced small intestine tumors in mice.

Inferred inactivation of the Cftr gene in the duodena of mice exposed to hexavalent chromium (Cr(VI)) in drinking water supports its tumor-suppressor status and implies its potential role in Cr(VI)-induced carcinogenesis of the small intestines

Carcinogenicity of hexavalent chromium [Cr (VI)] has been supported by a number of epidemiological and animal studies; however, its carcinogenic mode of action is still incompletely understood. To identify mechanisms involved in cancer development, we analyzed gene expression data from duodena of mice exposed to Cr(VI) in drinking water. This analysis included (i) identification of upstream regulatory molecules that are likely responsible for the observed gene expression changes, (ii) identification of annotated gene expression data from public repositories that correlate with gene expression changes in duodena of Cr(VI)-exposed mice, and (iii) identification of hallmark and oncogenic signature gene sets relevant to these data. We identified the inactivated CFTR gene among the top scoring upstream regulators, and found positive correlations between the expression data from duodena of Cr(VI)-exposed mice and other datasets in public repositories associated with the inactivation of the CFTR gene. In addition, we found enrichment of signatures for oncogenic signaling, sustained cell proliferation, impaired apoptosis and tissue remodeling. Results of our computational study support the tumor-suppressor role of the CFTR gene. Furthermore, our results support human relevance of the Cr(VI)-mediated carcinogenesis observed in the small intestines of exposed mice and suggest possible groups that may be more vulnerable to the adverse outcomes associated with the inactivation of CFTR by hexavalent chromium or other agents. Lastly, our findings predict, for the first time, the role of CFTR inactivation in chemical carcinogenesis and expand the range of plausible mechanisms that may be operative in Cr(VI)-mediated carcinogenesis of intestinal and possibly other tissues.

Development of the genomic inflammatory index (GII) to assess key maternal antecedents associated with placental inflammation

INTRODUCTION

Placental inflammation is associated with a variety of adverse health outcomes, including poor pregnancy outcomes as well as later in life health. The current clinical methodologies for evaluating placental histology for inflammation are limited in their sensitivity. The objective of this study was to develop a genomic inflammatory index (GII) that can be utilized as a biomarker to effectively quantify and evaluate placental inflammation.

METHODS

RNA-sequencing of n = 386 placentas from the Extremely Low Gestational Age Newborn (ELGAN) cohort was conducted. Transcriptional data for a biologically-targeted set of 14 genes, selected for their established role in pro-inflammatory signaling pathways, were aggregated to construct the GII. Multiple linear regression models were used to examine relationships between 47 perinatal factors and the GII.

RESULTS

The GII demonstrated a nine-fold difference across subjects and displayed positive trends with other indicators of placental inflammation. Significant differences in the GII were observed for race where women who self-identified as Black displayed higher levels of placental inflammation than those who self-identified as White women (p < 0.001). Additionally, married Black women showed reduced placental inflammation compared to those who were unmarried (beta value: 0.828, p-value: 0.032). Placentas from women who were treated with steroids during the delivery of the infant displayed higher GII levels than those who were not (p = 0.023).

DISCUSSION

Overall, the GII demonstrated an association between various perinatal factors and placental inflammation. It is anticipated that the GII will provide a novel genomics tool for quantifying placental inflammation, allowing for further investigation of causes, and ultimately the prevention, of inflammation in the placenta.

Predictive modeling of biological responses in the rat liver using in vitro Tox21 bioactivity: Benefits from high-throughput toxicokinetics

Computational methods are needed to more efficiently leverage data from in vitro cell-based models to predict what occurs within whole body systems after chemical insults. This study set out to test the hypothesis that in vitro high-throughput screening (HTS) data can more effectively predict in vivo biological responses when chemical disposition and toxicokinetic (TK) modeling are employed. In vitro HTS data from the Tox21 consortium were analyzed in concert with chemical disposition modeling to derive nominal, aqueous, and intracellular estimates of concentrations eliciting 50% maximal activity. In vivo biological responses were captured using rat liver transcriptomic data from the DrugMatrix and TG-Gates databases and evaluated for pathway enrichment. In vivo dosing data were translated to equivalent body concentrations using HTTK modeling. Random forest models were then trained and tested to predict in vivo pathway-level activity across 221 chemicals using in vitro bioactivity data and physicochemical properties as predictor variables, incorporating methods to address imbalanced training data resulting from high instances of inactivity. Model performance was quantified using the area under the receiver operator characteristic curve (AUC-ROC) and compared across pathways for different combinations of predictor variables. All models that included toxicokinetics were found to outperform those that excluded toxicokinetics. Biological interpretation of the model features revealed that rather than a direct mapping of in vitro assays to in vivo pathways, unexpected combinations of multiple in vitro assays predicted in vivo pathway-level activities. To demonstrate the utility of these findings, the highest-performing model was leveraged to make new predictions of in vivo biological responses across all biological pathways for remaining chemicals tested in Tox21 with adequate data coverage (n = 6617). These results demonstrate that, when chemical disposition and toxicokinetics are carefully considered, in vitro HT screening data can be used to effectively predict in vivo biological responses to chemicals.

An adverse outcome pathway for small intestinal tumors in mice involving chronic cytotoxicity and regenerative hyperplasia: a case study with hexavalent chromium, captan, and folpet

Small intestinal (SI) tumors are relatively uncommon outcomes in rodent cancer bioassays, and limited information regarding chemical-induced SI tumorigenesis has been reported in the published literature. Herein, we propose a cytotoxicity-mediated adverse outcome pathway (AOP) for SI tumors by leveraging extensive target species- and site-specific molecular, cellular, and histological mode of action (MOA) research for three reference chemicals, the fungicides captan and folpet and the transition metal hexavalent chromium (Cr(VI)). The gut barrier functions through highly efficient homeostatic regulation of SI epithelial cell sloughing, regenerative proliferation, and repair, which involves the replacement of up to 10¹¹ cells per day. This dynamic turnover in the SI provides a unique local environment for a cytotoxicity mediated AOP/MOA. Upon entering the duodenum, cytotoxicity to the villous epithelium is the molecular initiating event, as indicated by crypt elongation, villous atrophy/blunting, and other morphologic changes. Over time, the regenerative capacity of the gut epithelium to compensate declines as epithelial loss accelerates, especially at higher exposures. The first key event (KE), sustained regenerative crypt proliferation/hyperplasia, requires sufficient durations, likely exceeding 6 or 12 months, due to extensive repair capacity, to create more opportunities for the second KE, spontaneous mutation/transformation, ultimately leading to proximal SI tumors. Per OECD guidance, biological plausibility, essentiality, and empirical support were assessed using modified Bradford Hill considerations. The weight-of-evidence also included a lack of induced mutations in the duodenum after up to 90 days of Cr(VI) or captan exposure. The extensive evidence for this AOP, along with the knowledge that human exposures are orders of magnitude below those associated with KEs in this AOP, supports its use for regulatory applications, including hazard identification and risk assessment.

Identifying the Transcriptional Response of Cancer and Inflammation-Related Genes in Lung Cells in Relation to Ambient Air Chemical Mixtures in Houston, Texas

Atmospheric pollution represents a complex mixture of air chemicals that continually interact and transform, making it difficult to accurately evaluate associated toxicity responses representative of real-world exposure. This study leveraged data from a previously published article and reevaluated lung cell transcriptional response induced by outdoor atmospheric pollution mixtures using field-based exposure conditions in the industrialized Houston Ship Channel. The tested hypothesis was that individual and co-occurring chemicals in the atmosphere relate to altered expression of critical genes involved in inflammation and cancer-related processes in lung cells. Human lung cells were exposed at an air-liquid interface to ambient air mixtures for 4 h, with experiments replicated across 5 days. Real-time monitoring of primary and secondary gas-phase pollutants, as well as other atmospheric conditions, was simultaneously conducted. Transcriptional analysis of exposed cells identified critical genes showing differential expression associated with both individual and chemical mixtures. The individual pollutant identified with the largest amount of associated transcriptional response was benzene. Tumor necrosis factor (TNF) and interferon regulatory factor 1 (IRFN1) were identified as key upstream transcription factor regulators of the cellular response to benzene. This study is among the first to measure lung cell transcriptional responses in relation to real-world, gas-phase air mixtures.

A role for microRNAs in the epigenetic control of sexually dimorphic gene expression in the human placenta

Aim: The contribution of miRNAs as epigenetic regulators of sexually dimorphic gene expression in the placenta is unknown. Materials & methods: 382 placentas from the extremely low gestational age newborns (ELGAN) cohort were evaluated for expression levels of 37,268 mRNAs and 2,102 miRNAs using genome-wide RNA-sequencing. Differential expression analysis was used to identify differences in the expression based on the sex of the fetus. Results: Sexually dimorphic expression was observed for 128 mRNAs and 59 miRNAs. A set of 25 miRNA master regulators was identified that likely contribute to the sexual dimorphic mRNA expression. Conclusion: These data highlight sex-dependent miRNA and mRNA patterning in the placenta and provide insight into a potential mechanism for observed sex differences in outcomes.

Comparison of Gene Expression Responses in the Small Intestine of Mice Following Exposure to 3 Carcinogens Using the S1500+ Gene Set Informs a Potential Common Adverse Outcome Pathway

Carcinogenesis of the small intestine is rare in humans and rodents. Oral exposure to hexavalent chromium (Cr(VI)) and the fungicides captan and folpet induce intestinal carcinogenesis in mice. Previously (Toxicol Pathol. 330:48-52), we showed that B6C3F1 mice exposed to carcinogenic concentrations of Cr(VI), captan, or folpet for 28 days exhibited similar histopathological responses including villus enterocyte cytotoxicity and regenerative crypt epithelial hyperplasia. Herein, we analyze transcriptomic responses from formalin-fixed, paraffin-embedded duodenal sections from the aforementioned study. TempO-Seq technology and the S1500+ gene set were used to analyze transcription responses. Transcriptional responses were similar between all 3 agents; gene-level comparison identified 126/546 (23%) differentially expressed genes altered in the same direction, with a total of 25 upregulated pathways. These changes were related to cellular metabolism, stress, inflammatory/immune cell response, and cell proliferation, including upregulation in hypoxia inducible factor 1 (HIF-1) and activator protein 1 (AP1) signaling pathways, which have also been shown to be related to intestinal injury and angiogenesis/carcinogenesis. The similar molecular-, cellular-, and tissue-level changes induced by these 3 carcinogens can be informative for the development of an adverse outcome pathway for intestinal cancer.

Review of transcriptomic responses to hexavalent chromium exposure in lung cells supports a role of epigenetic mediators in carcinogenesis

Inhalation exposure to hexavalent chromium [Cr(VI)] is associated with increased risk of lung cancer with a mode of action (MOA) postulated to involve non-mutagenic key events, yet molecular-level events remain uncertain. Previously-published transcriptomic studies in the lung and lung cells were reviewed to evaluate molecular events in the MOA. This study aimed to (i) identify biological pathways that are consistently modulated by Cr(VI) in the lung through the compilation of transcriptomic-based databases, (ii) predict interactions between epigenetic regulators and transcriptional responses, and (iii) relate findings to previous literature to postulate a mechanism of action underlying Cr(VI)-induced lung cancer involving changes in genomic/epigenomic signatures. This cross-study comparison identified 372 genes with Cr(VI)-induced expression alterations in multiple studies. Pathway enrichment analyses of the commonly modulated genes demonstrated that pathways involved in cytotoxicity / cell proliferation were highly enriched, as well as the general suppression of genes involved in DNA damage repair. These signaling alterations were predicted to be regulated by DNA methylation, histone modifications, and microRNAs; and published evidence substantiates the role of these epigenetic regulators in Cr(VI)-induced carcinogenicity. Findings support the influence of epigenetic alterations on cell signaling related to Cr(VI)-induced cytotoxicity/cell proliferation, and decreases in DNA repair signaling leading to tumorigenesis.

Identifying Attributes That Influence In Vitro-to-In Vivo Concordance by Comparing In Vitro Tox21 Bioactivity Versus In Vivo DrugMatrix Transcriptomic Responses Across 130 Chemicals

Recent efforts aimed at integrating in vitro high-throughput screening (HTS) data into chemical toxicity assessments are necessitating increased understanding of concordance between chemical-induced responses observed in vitro versus in vivo. This investigation set out to (1) measure concordance between in vitro HTS data and transcriptomic responses observed in vivo, focusing on the liver, and (2) identify attributes that can influence concordance. Signal response profiles from 130 substances were compared between in vitro data produced through Tox21 and liver transcriptomic data through DrugMatrix, collected from rats exposed to a chemical for ≤5 days. A global in vitro-to-in vivo comparative analysis based on pathway-level responses resulted in an overall average percent agreement of 79%, ranging on a per-chemical basis between 41% and 100%. Whereas concordance amongst inactive chemicals was high (89%), concordance amongst chemicals showing in vitro activity was only 13%, suggesting that follow-up in vivo and/or orthogonal in vitro assays would improve interpretations of in vitro activity. Attributes identified to influence concordance included experimental design attributes (eg, cell type), target pathways, and physicochemical properties (eg, logP). The attribute that most consistently increased concordance was dose applicability, evaluated by filtering for experimental doses administered to rats that were within 10-fold of those related to likely bioactivity, derived using Tox21 data and high-throughput toxicokinetic modeling. Together, findings suggest that in vitro screening approaches to predict in vivo toxicity are viable particularly when certain attributes are considered, including whether activity versus inactivity is observed, experimental design, chemical properties, and dose applicability.

A review of the genotoxic, mutagenic, and carcinogenic potentials of several lower acrylates

Lower alkyl acrylate monomers include methyl-, ethyl-, n-butyl-, and 2-ethylhexyl acrylate. These acrylates are used in the manufacture of acrylic polymers and copolymers for plastics, food packaging, adhesives, and cosmetic formulations. Although there is limited potential for human environmental exposure, occupational exposure can occur via inhalation and dermal contact. Recently, new genotoxicity data have been generated, along with in silico and in vitro read-cross analyses, for these acrylates. The availability of high-throughput screening (HTS) data through the ToxCast™/Tox21 databases allows for consideration of computational toxicology and organization of these data according to the ten key characteristics of carcinogens. Therefore, we conducted a comprehensive review to evaluate the mechanistic, toxicokinetic, animal, and human data, including HTS data, for characterizing the potential carcinogenicity, mutagenicity, and genotoxicity of these acrylates. Toxicokinetic data demonstrate that these acrylates are metabolized rapidly by carboxylesterase hydrolysis and conjugation with glutathione. HTS data demonstrated an overall lack of bioactivity in cancer-related pathways. Overall, the genotoxicity and mutagenicity data support a cytotoxic, non-genotoxic mechanism for these acrylates. Cancer bioassay studies conducted by the oral, dermal, and inhalation routes in animal models with these acrylates did not show any increase in tumor incidence, with two exceptions. At high doses, and secondary to chronic site-of-contact irritation and corrosion, rodent forestomach tumors were induced by oral gavage dosing with ethyl acrylate, and skin tumors were observed following chronic dermal dosing with 2-ethylhexyl acrylate in C3H/HeJ inbred mice (a strain with deficiencies in wound healing), but not in the outbred NMRI strain. For both dermal and forestomach cancers, tumorigenesis is secondary to high doses and long-term tissue damage, shown to be reversible. With evidence that these chemicals are not genotoxic, and that they cause forestomach and dermal tumors through chronic irritation and regenerative proliferation mechanisms, these acrylates are unlikely to pose a human cancer hazard.

Integration of mechanistic and pharmacokinetic information to derive oral reference dose and margin-of-exposure values for hexavalent chromium

The current US Environmental Protection Agency (EPA) reference dose (RfD) for oral exposure to chromium, 0.003 mg kg-1  day-1 , is based on a no-observable-adverse-effect-level from a 1958 bioassay of rats exposed to ≤25 ppm hexavalent chromium [Cr(VI)] in drinking water. EPA characterizes the confidence in this RfD as "low." A more recent cancer bioassay indicates that Cr(VI) in drinking water is carcinogenic to mice at ≥30 ppm. To assess whether the existing RfD is health protective, neoplastic and non-neoplastic lesions from the 2 year cancer bioassay were modeled in a three-step process. First, a rodent physiological-based pharmacokinetic (PBPK) model was used to estimate internal dose metrics relevant to each lesion. Second, benchmark dose modeling was conducted on each lesion using the internal dose metrics. Third, a human PBPK model was used to estimate the daily mg kg-1 dose that would produce the same internal dose metric in both normal and susceptible humans. Mechanistic research into the mode of action for Cr(VI)-induced intestinal tumors in mice supports a threshold mechanism involving intestinal wounding and chronic regenerative hyperplasia. As such, an RfD was developed using incidence data for the precursor lesion diffuse epithelial hyperplasia. This RfD was compared to RfDs for other non-cancer endpoints; all RfD values ranged 0.003-0.02 mg kg-1  day-1 . The lowest of these values is identical to EPA's existing RfD value. Although the RfD value remains 0.003 mg kg-1  day-1 , the confidence is greatly improved due to the use of a 2-year bioassay, mechanistic data, PBPK models and benchmark dose modeling.

Benchmark Dose Modeling Estimates of the Concentrations of Inorganic Arsenic That Induce Changes to the Neonatal Transcriptome, Proteome, and Epigenome in a Pregnancy Cohort

Prenatal inorganic arsenic (iAs) exposure influences the expression of critical genes and proteins associated with adverse outcomes in newborns, in part through epigenetic mediators. The doses at which these genomic and epigenomic changes occur have yet to be evaluated in the context of dose-response modeling. The goal of the present study was to estimate iAs doses that correspond to changes in transcriptomic, proteomic, epigenomic, and integrated multi-omic signatures in human cord blood through benchmark dose (BMD) modeling. Genome-wide DNA methylation, microRNA expression, mRNA expression, and protein expression levels in cord blood were modeled against total urinary arsenic (U-tAs) levels from pregnant women exposed to varying levels of iAs. Dose-response relationships were modeled in BMDExpress, and BMDs representing 10% response levels were estimated. Overall, DNA methylation changes were estimated to occur at lower exposure concentrations in comparison to other molecular endpoints. Multi-omic module eigengenes were derived through weighted gene co-expression network analysis, representing co-modulated signatures across transcriptomic, proteomic, and epigenomic profiles. One module eigengene was associated with decreased gestational age occurring alongside increased iAs exposure. Genes/proteins within this module eigengene showed enrichment for organismal development, including potassium voltage-gated channel subfamily Q member 1 (KCNQ1), an imprinted gene showing differential methylation and expression in response to iAs. Modeling of this prioritized multi-omic module eigengene resulted in a BMD(BMDL) of 58(45) μg/L U-tAs, which was estimated to correspond to drinking water arsenic concentrations of 51(40) μg/L. Results are in line with epidemiological evidence supporting effects of prenatal iAs occurring at levels <100 μg As/L urine. Together, findings present a variety of BMD measures to estimate doses at which prenatal iAs exposure influences neonatal outcome-relevant transcriptomic, proteomic, and epigenomic profiles.