Combining phenotypic profiling and targeted RNA-Seq reveals linkages between transcriptional perturbations and chemical effects on cell morphology: Retinoic acid as an example.

The United States Environmental Protection Agency has proposed a tiered testing strategy for chemical hazard evaluation based on new approach methods (NAMs). The first tier includes in vitro profiling assays applicable to many (human) cell types, such as high-throughput transcriptomics (HTTr) and high-throughput phenotypic profiling (HTPP). The goals of this study were to: (1) harmonize the seeding density of U-2 OS human osteosarcoma cells for use in both assays; (2) compare HTTr- versus HTPP-derived potency estimates for 11 mechanistically diverse chemicals; (3) identify candidate reference chemicals for monitoring assay performance in future screens; and (4) characterize the transcriptional and phenotypic changes in detail for all-trans retinoic acid (ATRA) as a model compound known for its adverse effects on osteoblast differentiation. The results of this evaluation showed that (1) HTPP conducted at low (400 cells/well) and high (3000 cells/well) seeding densities yielded comparable potency estimates and similar phenotypic profiles for the tested chemicals; (2) HTPP and HTTr resulted in comparable potency estimates for changes in cellular morphology and gene expression, respectively; (3) three test chemicals (etoposide, ATRA, dexamethasone) produced concentration-dependent effects on cellular morphology and gene expression that were consistent with known modes-of-action, demonstrating their suitability for use as reference chemicals for monitoring assay performance; and (4) ATRA produced phenotypic changes that were highly similar to other retinoic acid receptor activators (AM580, arotinoid acid) and some retinoid X receptor activators (bexarotene, methoprene acid). This phenotype was observed concurrently with autoregulation of the RARB gene. Both effects were prevented by pre-treating U-2 OS cells with pharmacological antagonists of their respective receptors. Thus, the observed phenotype could be considered characteristic of retinoic acid pathway activation in U-2 OS cells. These findings lay the groundwork for combinatorial screening of chemicals using HTTr and HTPP to generate complementary information for the first tier of a NAM-based chemical hazard evaluation strategy.

Identifying candidate reference chemicals for in vitro testing of the retinoid pathway for predictive developmental toxicity.

Evaluating chemicals for potential in vivo toxicity based on their in vitro bioactivity profile is an important step toward animal- free testing. A compendium of reference chemicals and data describing their bioactivity on specific molecular targets, cellular pathways, and biological processes is needed to bolster confidence in the predictive value of in vitro hazard detection. Endogenous signaling by all-trans retinoic acid (ATRA) is an important pathway in developmental processes and toxicities. Employing data extraction methods and advanced literature extraction tools, we assembled a set of candidate reference chemicals with demonstrated activity on ten protein family targets in the retinoid system. The compendium was culled from Protein Data Bank, ChEMBL, ToxCast/Tox21, and the biomedical literature in PubMed. Finally, we performed a case study on one chemical in our collection, citral, an inhibitor of endogenous ATRA production, to determine whether the literature supports an adverse outcome pathway explaining the compound's developmental toxicity initiated by disruption of the retinoid pathway. We also deliver an updated Abstract Sifter tool populated with these reference compounds and complex search terms designed to query the literature for the downstream consequences to support concordance with targeted retinoid pathway disruption.

Viability of cultured human skin cells treated with 1,6-hexamethylene diisocyanate monomer and its oligomer isocyanurate in different culture media.

The isocyanate monomer 1,6-hexamethylene diisocyanate (HDI) and one of its trimers, HDI isocyanurate, are airway and skin sensitizers contained in polyurethane paint. The toxic response of cultured skin cells to these compounds was measured by evaluating the isocyanate concentrations at which 50% of the cells die (i.e., lethal concentration 50%, LC50) because the relative toxicity of each form of HDI should be considered when exposure limits of HDI-based paints are set. By using a luminescent ATP-viability assay, we compared the cytotoxic effects of HDI monomer and HDI isocyanurate on cultured human skin cells (keratinocytes, fibroblasts, and melanocytes) after 4-h isocyanate exposures using culture media with varying levels of nutrients in order to also determine the effects of media composition on isocyanate toxicity. Before analysis, experimental wells were normalized to controls containing cells that were cultured with the same vehicle and media. The measured mean LC50 values ranged from 5 to 200 µM across the experimental conditions, in which HDI isocyanurate in protein-devoid media was the most toxic to cells, producing the lowest LC50 values. For HDI monomer, keratinocytes were the most resistant to its toxicity and melanocytes were the most susceptible. However, when exposed to HDI isocyanurate, the opposite was observed, with melanocytes being the most resilient and the keratinocytes and fibroblasts were more susceptible. Depending on the type of skin cells, dose-response data indicated that HDI isocyanurate was 2-6 times more toxic than HDI monomer when using protein-devoid media whereas HDI isocyanurate was 4-13 times more toxic than HDI monomer when protein-rich media was used. Therefore, if the protein-devoid saline medium alone were used for these experiments, then a significant under-estimation of their relative toxicities in protein-rich environments would have resulted. This difference is because HDI monomer toxicity was more attenuated by the presence of protein in the culture media than HDI isocyanurate toxicity. Thus, conclusions based on comparative toxicity studies and consequent inference applied to potential human toxicity can be affected by in vitro culture media conditions. The physiochemical difference in reactivity of the two forms of HDI to biological molecules most likely explains the observed toxicity differences and may have implications for skin penetration, adverse effects like skin sensitization, and systemic responses like asthma. Future studies are warranted to investigate differences in the biological availability, cellular toxicity, and immunologic sensitization mechanisms for HDI monomer and HDI isocyanurate.

Epigenetic Markers Are Associated With Differences in Isocyanate Biomarker Levels in Exposed Spray-Painters.

Isocyanates are respiratory and skin sensitizers that are one of the main causes of occupational asthma globally. Genetic and epigenetic markers are associated with isocyanate-induced asthma and, before asthma develops, we have shown that genetic polymorphisms are associated with variation in plasma and urine biomarker levels in exposed workers. Inter-individual epigenetic variance may also have a significant role in the observed biomarker variability following isocyanate exposure. Therefore, we determined the percent methylation for CpG islands from DNA extracted from mononuclear blood cells of 24 male spray-painters exposed to 1,6-hexamethylene diisocyanate (HDI) monomer and HDI isocyanurate. Spray-painters' personal inhalation and skin exposure to these compounds and the respective biomarker levels of 1,6-diaminohexane (HDA) and trisaminohexyl isocyanurate (TAHI) in their plasma and urine were measured during three repeated industrial hygiene monitoring visits. We controlled for inhalation exposure, skin exposure, age, smoking status, and ethnicity as covariates and performed an epigenome-wide association study (EWAS) using likelihood-ratio statistical modeling. We identified 38 CpG markers associated with differences in isocyanate biomarker levels (Bonferroni < 0.05). Annotations for these markers included 18 genes: ALG1, ANKRD11, C16orf89, CHD7, COL27A, FUZ, FZD9, HMGN1, KRT6A, LEPR, MAPK10, MED25, NOSIP, PKD1, SNX19, UNC13A, UROS, and ZFHX3. We explored the functions of the genes that have been published in the literature and used GeneMANIA to investigate gene ontologies and predicted protein-interaction networks. The protein functions of the predicted networks include keratinocyte migration, cell-cell adhesions, calcium transport, neurotransmitter release, nitric oxide production, and apoptosis regulation. Many of the protein pathway functions overlap with previous findings on genetic markers associated with variability both in isocyanate biomarker levels and asthma susceptibility, which suggests there are overlapping protein pathways that contribute to both isocyanate toxicokinetics and toxicodynamics. These predicted protein networks can inform future research on the mechanism of allergic airway sensitization by isocyanates and aid in the development of mitigation strategies to better protect worker health.

Influence of Genetic Variance on Biomarker Levels After Occupational Exposure to 1,6-Hexamethylene Diisocyanate Monomer and 1,6-Hexamethylene Diisocyanate Isocyanurate.

We evaluated the impact of genetic variance on biomarker levels in a population of workers in the automotive repair and refinishing industry who were exposed to respiratory sensitizers 1,6-hexamethylene diisocyanate (HDI) monomer and one of its trimers, HDI isocyanurate. The exposures and respective urine and plasma biomarkers 1,6-diaminohexane (HDA) and trisaminohexyl isocyanurate (TAHI) were measured in 33 workers; and genome-wide microarrays (Affymetrix 6.0) were used to genotype the workers' single-nucleotide polymorphisms (SNPs). Linear mixed model analyses have indicated that interindividual variations in both inhalation and skin exposures influenced these biomarker levels. Using exposure values as covariates and a false discovery rate < 0.10 to assess statistical significance, we observed that seven SNPs were associated with HDA in plasma, five were associated with HDA in urine, none reached significance for TAHI in plasma, and eight were associated with TAHI levels in urine. The different genotypes for the 20 significant SNPs accounted for 4- to 16-fold changes observed in biomarker levels. Associated gene functions include transcription regulation, calcium ion transport, vascular morphogenesis, and transforming growth factor beta signaling pathway, which may impact toxicokinetics indirectly by altering inflammation levels. Additionally, in an expanded analysis using a minor allele cutoff of 0.05 instead of 0.10, there were biomarker-associated SNPs within three genes that have been associated with isocyanate-induced asthma: ALK, DOCK2, and LHPP. We demonstrate that genetic variance impacts the biomarker levels in workers exposed to HDI monomer and HDI isocyanurate and that genetics can be used to refine exposure predictions in small cohorts when quantitative personal exposure and biomarker measurements are included in the models.