Determination of in vitro hepatotoxic potencies of a series of perfluoroalkyl substances (PFASs) based on gene expression changes in HepaRG liver cells

Comparative study on gene expression profiles in the liver of male neonatal mice prenatally exposed to PFOA and its alternative HFPO-DA

Hexafluoropropylene oxide dimer acid (HFPO-DA), which belongs to the class of perfluoroalkyl ether carboxylic acid (PFECA), is a new alternative to perfluorooctanoic acid (PFOA). However, whether HFPO-DA is a safer alternative to PFOA in neonates remains unclear. In this study, we evaluated neonatal hepatic toxicity on postnatal days 9-10 by orally exposing pregnant CD-1 mice to 0.3 or 3.0 mg/kg/day (low or high doses) of HFPO-DA or PFOA from gestation days 15-17. The results showed that exposure of pregnant mice to HFPO-DA and PFOA induced similar phenotypic effects, including significant decreases in neonatal body weight (BW) and significant increases in liver weight relative to BW in the high-dose. Notably, HFPO-DA exposure significantly decreased in neonatal BW in the low-dose group, whereas PFOA did not. Comprehensive gene expression analysis revealed significant alterations in 408 and 1402 differentially expressed genes (DEGs) in the liver of neonates from the low- and high-dose HFPO-DA groups, respectively, while PFOA significantly altered 0 and 292 DEGs in the corresponding groups. Gene set enrichment analysis indicated that the DEGs induced by HFPO-DA and PFOA were enriched in pathway related to "PPAR signaling", "fatty acid metabolism", and "biological oxidations". In addition, transactivation assays revealed that mouse (m)PPARα and mPPARγ activity of HFPO-DA exceeds that of PFOA and molecular docking simulations analysis predicted that the binding conformation differ between PFOA and HFPO-DA. Overall, our findings demonstrate that HFPO-DA consistently affected neonatal phenotypes, liver gene expression and the molecular initiating events involving PPARα/γ, at lower concentrations than PFOA.

Investigating the Mechanism of Neurotoxic Effects of PFAS in Differentiated Neuronal Cells through Transcriptomics and Lipidomics Analysis

Per- and polyfluorinated alkyl substances (PFAS) are pervasive environmental contaminants that bioaccumulate in tissues and pose risks to human health. Increasing evidence links PFAS to neurodegenerative and behavioral disorders, yet the underlying mechanisms of their effects on neuronal function remain largely unexplored. In this study, we utilized SH-SY5Y neuroblastoma cells, differentiated into neuronal-like cells, to investigate the impact of six PFAS compounds─perfluorooctanoic acid (PFOA), perfluorooctanesulfonic acid (PFOS), perfluorodecanoic acid (PFDA), perfluorodecanesulfonic acid (PFDS), 8:2 fluorotelomer sulfonate (8:2 FTS), and 8:2 fluorotelomer alcohol (8:2 FTOH)─on neuronal health. Following a 30 μM exposure for 24 h, PFAS accumulation ranged from 40-6500 ng/mg of protein. Transcriptomic analysis revealed 721 differentially expressed genes (DEGs) across treatments (padj < 0.05), with 11 DEGs shared among all PFAS exposures, indicating potential biomarkers for neuronal PFAS toxicity. PFOA-treated cells showed downregulation of genes involved in synaptic growth and neural function, while PFOS, PFDS, 8:2 FTS, and 8:2 FTOH exposures resulted in the upregulation of genes related to hypoxia response and amino acid metabolism. Lipidomic profiling further demonstrated significant increases in fatty acid levels with PFDA, PFDS, and 8:2 FTS and depletion of triacylglycerols with 8:2 FTOH treatments. These findings suggest that the neurotoxic effects of PFAS are structurally dependent, offering insights into the molecular processes that may drive PFAS-induced neuronal dysfunction.

Time-course characterization of whole-transcriptome dynamics of HepG2/C3A spheroids and its toxicological implications

Physiologically relevant in vitro models are a priority in predictive toxicology to replace and/or reduce animal experiments. The compromised toxicant metabolism of many immortalized human liver cell lines grown as monolayers as compared to in vivo metabolism limits their physiological relevance. However, recent efforts to culture liver cells in a 3D environment, such as spheroids, to better mimic the in vivo conditions, may enhance the toxicant metabolism of human liver cell lines. In this study, we characterized the dynamic changes in the transcriptome of HepG2/C3A hepatocarcinoma cell spheroids maintained in a clinostat system (CelVivo) to gain insight into the metabolic capacity of this model as a function of spheroid size and culture time. We assessed morphological changes (size, necrotic core), cell health, and proliferation rate from initial spheroid seeding to 35 days of continuous culture in conjunction with a time-course (0, 3, 7, 10, 14, 21, 28 days) of the transcriptome (TempO-Seq, BioSpyder). The phenotypic characteristics of HepG2/C3A growing in spheroids were comparable to monolayer growth until ∼Day 12 (Day 10-14) when a significant decrease in cell doubling rate was noted which was concurrent with down-regulation of cell proliferation and cell cycle pathways over this time period. Principal component analysis of the transcriptome data suggests that the Day 3, 7, and 10 spheroids are pronouncedly different from the Day 14, 21, and 28 spheroids in support of a biological transition time point during the long-term 3D spheroid cultures. The expression of genes encoding cellular components involved in toxicant metabolism and transport rapidly increased during the early time points of spheroids to peak at Day 7 or Day 10 as compared to monolayer cultures with a gradual decrease in expression with further culture, suggesting the most metabolically responsive time window for exposure studies. Overall, we provide baseline information on the cellular and molecular characterization, with a particular focus on toxicant metabolic capacity dynamics and cell growth, of HepG2/C3A 3D spheroid cultures over time.

Metabolic and phenotypic changes induced by PFAS exposure in two human hepatocyte cell models

PFAS are ubiquitous industrial chemicals with known adverse health effects, particularly on the liver. The liver, being a vital metabolic organ, is susceptible to PFAS-induced metabolic dysregulation, leading to conditions such as hepatotoxicity and metabolic disturbances. In this study, we investigated the phenotypic and metabolic responses of PFAS exposure using two hepatocyte models, HepG2 (male cell line) and HepaRG (female cell line), aiming to define phenotypic alterations, and metabolic disturbances at the metabolite and pathway levels. The PFAS mixture composition was selected based on epidemiological data, covering a broad concentration spectrum observed in diverse human populations. Phenotypic profiling by Cell Painting assay disclosed predominant effects of PFAS exposure on mitochondrial structure and function in both cell models as well as effects on F-actin, Golgi apparatus, and plasma membrane-associated measures. We employed comprehensive metabolic characterization using liquid chromatography combined with high-resolution mass spectrometry (LC-HRMS). We observed dose-dependent changes in the metabolic profiles, particularly in lipid, steroid, amino acid and sugar and carbohydrate metabolism in both cells as well as in cell media, with HepaRG cell line showing a stronger metabolic response. In cells, most of the bile acids, acylcarnitines and free fatty acids showed downregulation, while medium-chain fatty acids and carnosine were upregulated, while the cell media showed different response especially in relation to the bile acids in HepaRG cell media. Importantly, we observed also nonmonotonic response for several phenotypic features and metabolites. On the pathway level, PFAS exposure was also associated with pathways indicating oxidative stress and inflammatory responses. Taken together, our findings on PFAS-induced phenotypic and metabolic disruptions in hepatocytes shed light on potential mechanisms contributing to the broader comprehension of PFAS-related health risks.

Exposure to PFAS chemicals induces sex-dependent alterations in key rate-limiting steps of lipid metabolism in liver steatosis

Toxicants with the potential to bioaccumulate in humans and animals have long been a cause for concern, particularly due to their association with multiple diseases and organ injuries. Per- and polyfluoro alkyl substances (PFAS) and polycyclic aromatic hydrocarbons (PAH) are two such classes of chemicals that bioaccumulate and have been associated with steatosis in the liver. Although PFAS and PAH are classified as chemicals of concern, their molecular mechanisms of toxicity remain to be explored in detail. In this study, we aimed to identify potential mechanisms by which an acute exposure to PFAS and PAH chemicals can induce lipid accumulation and whether the responses depend on chemical class, dose, and sex. To this end, we analyzed mechanisms beginning with the binding of the chemical to a molecular initiating event (MIE) and the consequent transcriptomic alterations. We collated potential MIEs using predictions from our previously developed ToxProfiler tool and from published steatosis adverse outcome pathways. Most of the MIEs are transcription factors, and we collected their target genes by mining the TRRUST database. To analyze the effects of PFAS and PAH on the steatosis mechanisms, we performed a computational MIE-target gene analysis on high-throughput transcriptomic measurements of liver tissue from male and female rats exposed to either a PFAS or PAH. The results showed peroxisome proliferator-activated receptor (PPAR)-α targets to be the most dysregulated, with most of the genes being upregulated. Furthermore, PFAS exposure disrupted several lipid metabolism genes, including upregulation of fatty acid oxidation genes (Acadm, Acox1, Cpt2, Cyp4a1-3) and downregulation of lipid transport genes (Apoa1, Apoa5, Pltp). We also identified multiple genes with sex-specific behavior. Notably, the rate-limiting genes of gluconeogenesis (Pck1) and bile acid synthesis (Cyp7a1) were specifically downregulated in male rats compared to female rats, while the rate-limiting gene of lipid synthesis (Scd) showed a PFAS-specific upregulation. The results suggest that the PPAR signaling pathway plays a major role in PFAS-induced lipid accumulation in rats. Together, these results show that PFAS exposure induces a sex-specific multi-factorial mechanism involving rate-limiting genes of gluconeogenesis and bile acid synthesis that could lead to activation of an adverse outcome pathway for steatosis.

DNA adductomics aided rapid screening of genotoxic impurities using nucleosides and 3D bioprinted human liver organoids

Current genotoxicity assessment methods are mainly employed to verify the genotoxic safety of drugs, but do not allow for rapid screening of specific genotoxic impurities (GTIs). In this study, a new approach for the recognition of GTIs has been proposed. It is to expose the complex samples to an in vitro nucleoside incubation model, and then draw complete DNA adduct profiles to infer the structures of potential genotoxic impurities (PGIs). Subsequently, the genotoxicity is confirmed in human by 3D bioprinted human liver organoids. To verify the feasibility of the approach, lansoprazole chloride compound (Lanchlor), a PGI during the synthesis of lansoprazole, was selected as the model drug. After confirming genotoxicity by Comet assay, it was exposed to different models to map and compare the DNA adduct profiles by LC-MS/MS. The results showed Lanchlor could generate diverse DNA adducts, revealing firstly its genotoxicity at molecular mechanism of action. Furthermore, the largest variety and content of DNA adducts were observed in the nucleoside incubation model, while the human liver organoids exhibited similar results with rats. The results showed that the combination of DNA adductomics and 3D bioprinted organoids were useful for the rapid screening of GTIs.

Comparative steroidogenic effects of hexafluoropropylene oxide trimer acid (HFPO-TA) and perfluorooctanoic acid (PFOA): Regulation of histone modifications

As a widely used alternative to perfluorooctanoic acid (PFOA), hexafluoropropylene oxide trimer acid (HFPO-TA) has been detected in the environment and humans; however, little is known regarding its male reproductive toxicity. To compare the effects of HFPO-TA on steroid hormone synthesis with PFOA, we exposed Leydig cells (MLTC-1) to non-lethal doses (0.1, 1, and 10 μM) of PFOA and HFPO-TA for 48 h. It was found that the levels of steroid hormones, 17α-hydroxyprogesterone (OHP), androstenedione (ASD), and testosterone (T) were significantly increased in 1 and 10 μM of PFOA and HFPO-TA groups, with greater elevation being observed in the HFPO-TA groups than in the PFOA groups at 10 μM. We further showed that the two rate-limiting steroidogenic genes (Star and Cyp11a1) were up-regulated, while Hsd3b, Cyp17a1, and Hsd17b were down-regulated or unchanged after PFOA/HFPO-TA exposure. Moreover, PFOA exposure significantly up-regulated histone H3K4me1/3 and H3K9me1, while down-regulated H3K4me2 and H3K9me2/3 levels. By contrast, H3K4me2/3 and H3K9me2/3 were enhanced, while H3K4me1 and H3K9me1 were repressed after HFPO-TA treatment. It was further confirmed that H3K4me1/3 were increased and H3K9me2 was decreased in Star and Cyp11a1 promoters by PFOA, while HFPO-TA increased H3K4me2/3 and decreased H3K9me1 in the two gene promoters. Therefore, we propose that low levels of PFOA/HFPO-TA enhance the expression of Star and Cyp11a1 by regulating H3K4 and H3K9 methylation, thus stimulating the production of steroid hormones in MLTC-1 cells. Collectively, HFPO-TA exhibits stronger effects on steroidogenesis compared to PFOA, which may be ascribed to the distinct regulation of histone modifications. These data suggest that HFPO-TA does not appear to be a safer alternative to PFOA on the aspect of male reproductive toxicity.

Hepatotoxicity and lipid metabolism disorders of 8:2 polyfluoroalkyl phosphate diester in zebrafish: In vivo and in silico evidence

8:2 polyfluoroalkyl phosphate diester (8:2 diPAP) has been shown to accumulate in the liver, but whether it induces hepatotoxicity and lipid metabolism disorders remains largely unknown. In this study, zebrafish embryos were exposed to 8:2 diPAP for 7 d. Hepatocellular hypertrophy and karyolysis were noted after exposure to 0.5 ng/L 8:2 diPAP, suggesting suppressed liver development. Compared to the water control, 8:2 diPAP led to significantly higher triglyceride and total cholesterol levels, but markedly lower levels of low-density lipoprotein, implying disturbed lipid homeostasis. The levels of two peroxisome proliferator activated receptor (PPAR) subtypes (pparα and pparγ) involved in hepatotoxicity and lipid metabolism were significantly upregulated by 8:2 diPAP, consistent with their overexpression as determined by immunohistochemistry. In silico results showed that 8:2 diPAP formed hydrogen bonds with PPARα and PPARγ. Among seven machine learning models, Adaptive Boosting performed the best in predicting the binding affinities of PPARα and PPARγ on the test set. The predicted binding affinity of 8:2 diPAP to PPARα (7.12) was higher than that to PPARγ (6.97) by Adaptive Boosting, which matched well with the experimental results. Our results revealed PPAR - mediated adverse effects of 8:2 diPAP on the liver and lipid metabolism of zebrafish larvae.

Exposure to environmental toxicants is associated with gut microbiome dysbiosis, insulin resistance and obesity

Environmental toxicants (ETs) are associated with adverse health outcomes. Here we hypothesized that exposures to ETs are linked with obesity and insulin resistance partly through a dysbiotic gut microbiota and changes in the serum levels of secondary bile acids (BAs). Serum BAs, per- and polyfluoroalkyl substances (PFAS) and additional twenty-seven ETs were measured by mass spectrometry in 264 Danes (121 men and 143 women, aged 56.6 ± 7.3 years, BMI 29.7 ± 6.0 kg/m2) using a combination of targeted and suspect screening approaches. Bacterial species were identified based on whole-genome shotgun sequencing (WGS) of DNA extracted from stool samples. Personalized genome-scale metabolic models (GEMs) of gut microbial communities were developed to elucidate regulation of BA pathways. Subsequently, we compared findings from the human study with metabolic implications of exposure to perfluorooctanoic acid (PFOA) in PPARα-humanized mice. Serum levels of twelve ETs were associated with obesity and insulin resistance. High chemical exposure was associated with increased abundance of several bacterial species (spp.) of genus (Anaerotruncus, Alistipes, Bacteroides, Bifidobacterium, Clostridium, Dorea, Eubacterium, Escherichia, Prevotella, Ruminococcus, Roseburia, Subdoligranulum, and Veillonella), particularly in men. Conversely, females in the higher exposure group, showed a decrease abundance of Prevotella copri. High concentrations of ETs were correlated with increased levels of secondary BAs including lithocholic acid (LCA), and decreased levels of ursodeoxycholic acid (UDCA). In silico causal inference analyses suggested that microbiome-derived secondary BAs may act as mediators between ETs and obesity or insulin resistance. Furthermore, these findings were substantiated by the outcome of the murine exposure study. Our combined epidemiological and mechanistic studies suggest that multiple ETs may play a role in the etiology of obesity and insulin resistance. These effects may arise from disruptions in the microbial biosynthesis of secondary BAs.

Risk-based prioritization of PFAS using phenotypic and transcriptomic data from human induced pluripotent stem cell-derived hepatocytes and cardiomyocytes

Per- and polyfluoroalkyl substances (PFAS) are chemicals with important applications; they are persistent in the environment and may pose human health hazards. Regulatory agencies are con­sidering restrictions and bans of PFAS; however, little data exists for informed decisions. Several prioritization strategies were proposed for evaluation of potential hazards of PFAS. Structure-based grouping could expedite the selection of PFAS for testing; still, the hypothesis that structure-effect relationships exist for PFAS requires confirmation. We tested 26 structurally diverse PFAS from 8 groups using human induced pluripotent stem cell-derived hepatocytes and cardiomyocytes, and tested concentration-response effects on cell function and gene expression. Few phenotypic effects were observed in hepatocytes, but negative chronotropy was observed in cardiomyocytes for 8 PFAS. Substance- and cell type-dependent transcriptomic changes were more prominent but lacked substantial group-specific effects. In hepatocytes, we found upregulation of stress-related and extracellular matrix organization pathways, and down-regulation of fat metabolism. In car­diomyocytes, contractility-related pathways were most affected. We derived phenotypic and transcriptomic points of departure and compared them to predicted PFAS exposures. Conservative estimates for bioactivity and exposure were used to derive a bioactivity-to-exposure ratio (BER) for each PFAS; 23 of 26 PFAS had BER > 1. Overall, these data suggest that structure-based PFAS grouping may not be sufficient to predict their biological effects. Testing of individual PFAS may be needed for scientifically-supported decision-making. Our proposed strategy of using two human cell types and considering phenotypic and transcriptomic effects, combined with dose-response analysis and calculation of BER, may be used for PFAS prioritization.

Impact of perfluoroalkyl substances (PFAS) and PFAS mixtures on lipid metabolism in differentiated HepaRG cells as a model for human hepatocytes

Per- and polyfluoroalkyl substances (PFAS) are environmental contaminants with various adverse health effects in humans including disruption of lipid metabolism. Aim of the present study was to elucidate the molecular mechanisms of PFAS-mediated effects on lipid metabolism in human cells. Here, we examined the impact of a number of PFAS (PFOS, PFOA, PFNA, PFDA, PFHxA, PFBA, PFHxS, PFBS, HFPO-DA, and PMPP) and of some exposure-relevant PFAS mixtures being composed of PFOS, PFOA, PFNA and PFHxS on lipid metabolism in human HepaRG cells, an in vitro model for human hepatocytes. At near cytotoxic concentrations, the selected PFAS and PFAS mixtures induced triglyceride accumulation in HepaRG cells and consistently affected the expression of marker genes for steatosis, as well as PPARα target genes and genes related to lipid and cholesterol metabolism, pointing to common molecular mechanisms of PFAS in disrupting cellular lipid and cholesterol homeostasis. PPARα activation was examined by a transactivation assay in HEK293T cells, and synergistic effects were observed for the selected PFAS mixtures at sum concentrations higher than 25 µM, whereas additivity was observed at sum concentrations lower than 25 µM. Of note, any effect observed in the in vitro assays occurred at PFAS concentrations that were at least four to five magnitudes above real-life internal exposure levels of the general population.

Environmental exposure to per- and perfluoroalkyl substances in early pregnancy and newborn anogenital distance: A prospective cohort study

Per- and polyfluoroalkyl substances (PFAS) are persistent, ubiquitous pollutants, and the current epidemiological evidence regarding the impact of in utero exposure to PFAS on anogenital distance (AGD) is limited and inconclusive. The primary aim of this study was to investigate the potential associations between maternal exposure to PFAS during pregnancy and AGD in newborns. A total of 2273 mother-child pairs were recruited for this study, and both PFAS levels and AGD were measured. Multiple linear regression models were utilized to explore the relationships between individual PFAS and AGD. Additionally, quantile-based g-computation (QGC) was employed to assess the joint effects of mixtures of PFAS on AGD. Our findings showed that maternal exposure to PFOS (β = 0.518, 95% CI: 0.093, 0.942), PFNA (β = 0.487, 95% CI: 0.037, 0.937), PFDA (β = 0.443, 95% CI: 0.048, 0.838), PFUA (β = 0.434, 95% CI: 0.031, 0.838), and PFBS (β = 0.444, 95% CI: 0.124, 0.763) during early pregnancy had a significant positive association with AGD in boys. Similarly, in girls, maternal exposure to PFOS (β = 0.423, 95% CI: 0.006, 0.841), PFNA (β = 0.641, 95% CI: 0.207, 1.074), PFDA (β = 0.670, 95% CI: 0.306, 1.033), PFUA (β = 0.895, 95% CI: 0.509, 1.281), and PFBS (β = 0.474, 95% CI: 0.178, 0.770) had a positive association with AGD, while PFOA (β = -1.254, 95% CI: -1.786, -0.723) had a negative association. QGC models further confirmed that PFAS mixtures were positively associated with AGD. Moreover, PFBS was the primary contributor to the joint effects of PFAS mixtures on AGD. In summary, our study has provided further corroboration for the possibility that PFAS exposure can have an impact on AGD in both boys and girls. The use of AGD as a promising biomarker for endocrine disruption highlights the significance of our findings, which may have valuable clinical implications for reproductive diseases.

Effects of perfluorooctanoic acid (PFOA) on gene expression profiles via nuclear receptors in HepaRG cells: Comparative study with in vitro transactivation assays

Perfluorooctanoic acid (PFOA), a synthetic perfluorinated eight-carbon organic chemical, has been reported to induce hepatotoxicity, including increased liver weight, hepatocellular hypertrophy, necrosis, and increased peroxisome proliferation in rodents. Epidemiological studies have demonstrated associations between serum PFOA levels and various adverse effects. In this study, we investigated the gene expression profiles of human HepaRG cells exposed to 10 and 100 μM PFOA for 24h. Treatment with 10 and 100 μM PFOA significantly modulated the expression of 190 genes and 996 genes, respectively. In particular, genes upregulated or downregulated by 100µM PFOA included peroxisome proliferator-activated receptor (PPAR) signaling genes related to lipid metabolism, adipocyte differentiation, and gluconeogenesis. In addition, we identified the "Nuclear receptors-meta pathways" following the activation of other nuclear receptors: constitutive androstane receptor (CAR), pregnane X receptor (PXR) and farnesoid X receptor (FXR), and the transcription factor, nuclear factor E2-related factor 2 (Nrf2). The expression levels of some target genes (CYP4A11, CYP2B6, CYP3A4, CYP7A1, and GPX2) of these nuclear receptors and Nrf2 were confirmed using quantitative reverse transcription polymerase chain reaction. Next, we performed transactivation assays using COS-7 or HEK293 cells to investigate whether these signaling-pathways were activated by the direct effects of PFOA on human PPARα, CAR, PXR, FXR and Nrf2. PFOA activated PPARα in a concentration-dependent manner, but did not activate CAR, PXR, FXR, or Nrf2. Taken together, these results suggest that PFOA affects the hepatic transcriptomic responses of HepaRG cells through direct activation of PPARα and indirect activation of CAR, PXR FXR and Nrf2. Our finding indicates that PPARα activation found in the "Nuclear receptors-meta pathways" functions as a molecular initiating event for PFOA, and indirect activation of alternative nuclear receptors and Nrf2 also provide important molecular mechanisms in PFOA-induced human hepatotoxicity.