Next Article in Journal
Effect of Breastmilk Microbiota and Sialylated Oligosaccharides on the Colonization of Infant Gut Microbial Community and Fecal Metabolome
Next Article in Special Issue
The Importance of Sharing Data in Systems Biology
Previous Article in Journal
Mechanisms and Outcomes of Metabolic Surgery in Type 2 Diabetes
Previous Article in Special Issue
Multi-Omics, an Integrated Approach to Identify Novel Blood Biomarkers of Alzheimer’s Disease
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exposure to Environmentally Relevant Levels of PFAS Causes Metabolic Changes in the Freshwater Amphipod Austrochiltonia subtenuis

by
Georgia M. Sinclair
1,*,
Sara M. Long
2,
Navneet Singh
3,
Timothy L. Coggan
3,
Matthew P. J. Askeland
3 and
Oliver A. H. Jones
1
1
Australian Centre for Research on Separation Science (ACROSS), School of Science, RMIT University, Bundoora West Campus, P.O. Box 71, Bundoora, VIC 3083, Australia
2
Aquatic Environmental Stress Research Group (AQUEST), School of Science, RMIT University, Bundoora West Campus, P.O. Box 71, Bundoora, VIC 3083, Australia
3
ADE Consulting Group, Williamstown North, VIC 3016, Australia
*
Author to whom correspondence should be addressed.
Metabolites 2022, 12(11), 1135; https://doi.org/10.3390/metabo12111135
Submission received: 26 October 2022 / Revised: 15 November 2022 / Accepted: 16 November 2022 / Published: 18 November 2022
(This article belongs to the Special Issue Integrated Systems Biology: Challenges and Future Perspectives)

Abstract

:
Per and polyfluoroalkyl substances (PFAS) are of concern to environmental regulators due to their widespread occurrence, persistence and reported toxicity. However, little data exist on the effects of PFAS at environmentally relevant concentrations. The development of molecular markers for PFAS exposure would therefore be useful to better understand the environmental risks of these compounds. In this study, we assessed if such markers could be developed using Gas Chromatography–Mass Spectrometry-based metabolomics. We exposed the freshwater amphipod Austrochiltonia subtenuis to a range of environmentally relevant concentrations of perfluoro-octane sulfonic acid (PFOS), hexafluoropropylene oxide dimer acid (GenX) and perfluorohexanesulphonic acid (PFHxS) for 7 days at five concentrations. A metabolic response was detected in all concentrations and treatments even though the survival rates only differed significantly at the highest exposure levels. The metabolic response differed between compounds but all three PFAS induced changes in the levels of amino acids, fatty acids, and cholesterol, in line with the literature. PFOS was found to bioaccumulate. Both GenX and PFHxS were eliminated from the amphipods, but PFHxS was eliminated at a slower rate than GenX. This information improves our understanding of the sublethal effects of PFAS as well as their environmental fate and behaviour.

Graphical Abstract

1. Introduction

Per and polyfluoroalkyl substances (PFAS) are a diverse group of human-made compounds that have been in use since the 1940s. They are used in a range of products such as Teflon and Scotchgard, food packaging, cosmetics, waterproof textiles and aqueous film-forming foams (AFFFs) [1]. PFAS are very resistant to degradation and are often labelled ‘forever chemicals’. They have been detected in almost every environmental compartment (e.g., soil, sediment, oceans and freshwater) as well as in wildlife and humans throughout the world [2].
It was previously assumed that PFAS were inert and little regard was given to the environmental fate and the ecological impacts of these compounds to organisms and environmental health. This was later realized to be a mistake and, today, there are serious concerns about their potential effects, causing growth/developmental delays and cancer [1]. While there is a substantial amount of data on PFAS toxicity, environmental exposure studies on these compounds tend to use substantially higher concentrations of PFAS than are generally found in the environment [3]. As such, it is still unclear what the biological effects of PFAS may be at environmentally relevant exposure levels [3,4]. Effective management of these compounds at the government, local and industrial levels is therefore challenging because it is difficult to prove that environmental harm has occurred.
Standard ecotoxicology endpoints, such as growth or reproduction, may not register the effects of chronic, low-dose toxicant exposure. However, recent work by Beale et al. [5] showed that environmental toxicological studies can be greatly enhanced when coupled with molecular biology endpoints, such as the use of metabolomics and other omic techniques. Omic techniques include the detection of genes (genomics), mRNA (transcriptomics), and proteins (proteomics) in a biological sample [5]. Metabolomics is the study of small-molecule metabolites such as fats, sugars and amino acids and how these change in response to external factors such as disease, or toxicant exposure [6]. Small changes in the network of metabolites can indicate the level of an organisms’ overall health and a report from the Expert Health Panel for PFAS noted that “better biomarkers of the net effect of all PFAS would be extremely useful” [7]. Finding the range of doses where PFAS concentrations disrupt the biochemistry of biota is therefore likely to enhance our overall understanding of the mode of action, bioaccumulation, and toxicity of PFAS in ecosystems [8,9].
Out of the thousands of potential PFAS, perfluorooctanoic acid (PFOA) and perfluoro-octane sulfonic acid (PFOS) are the most commonly detected in the environment, closely followed by perfluorohexanesulfonic acid (PFHxS). PFOA has already been listed by the European Chemical Agency as a substance of very high concern and, in 2022, UN countries agreed to a global ban of PFHxS, with no exemptions, by adding it to Annex A of the Stockholm Convention on persistent organic pollutants (POPs) [10]. This means the use of both PFOA and PFHxS is now restricted, even though they are still commonly detected in the environment [11]. The phase-out of these compounds has led to the introduction of shorter-chained compounds such as GenX (hexafluoropropylene oxide dimer acid). GenX has recently been detected in ecosystems in Europe and North America, indicating that it may be just as persistent as other PFAS compounds and, as such, it is of increasing concern worldwide [12,13].
This work focused on PFOS, PFHxS and GenX. Three experiments were carried out in controlled laboratory exposures, using the non-model organism Austrochiltonia subtenuis (an endemic and common freshwater amphipod in Victoria, Australia) [14]. They are small, usually <1 cm, benthic shredders and are a food source to many fish, birds and invertebrate species. This species was chosen as it is crucial to the function and structure of invertebrate communities. Gas Chromatography–Mass Spectrometry (GC–MS) was used to identify metabolites in amphipods and assess how they changed in response to environmentally relevant levels of PFAS. We were particularly interested to see if a common biochemical response would be evident across compounds. An additional objective was to calculate bioconcentration factors (BCFs) for each compound and use these to help understand the metabolomic data. The overall aim was to contribute to the understanding of the potential effects on PFAS at environmentally relevant exposure levels and, ultimately, inform environmental regulation of these compounds.

2. Materials and Methods

2.1. Organism

Live amphipods (Austrochiltonia subtenuis) were originally collected from a wild population in Deep Creek, Bulla Rd, Bulla, Victoria, Australia (S 37°37.919157′ E 144°47.995837′). An amphipod culture has long been established (>5 years) and maintained in house at RMIT University (Melbourne, Australia). The amphipods are maintained in aquaria using a standard artificial media (SAM) modified from Borgmann [15]. The SAM consisted of reverse osmosis (RO) water (18.2 mΩ), obtained in house, with 0.23 mM NaHCO3, 0.061 mM CaCl2, 0.032 mM MgSO4, 0.47 mM NaCl, 0.0087 mM KCl, 0.17 mM MgCl2 and 0.0009 mM NaBr. Cultures were maintained at 21 ± 1 °C under a 16:8 h light: dark photoperiod. The culture was fed powdered fish food (Tetramin®, Tetra Werke, Melle, Germany) and yeast–cerophyll–trout chow (YCT), made in house, every second day.

2.2. Chemical and Concentrations

A perfluoro-octane sulfonic acid (PFOS) standard was purchased from Sigma-Aldrich (Melbourne, VIC, Australia). Hexafluoropropylene oxide dimer acid (common name, GenX) and perfluorohexanesulphonic acid (PFHxS) standards were purchased from Wellington laboratories (Guelph, ON, Canada). Experimental solutions of differing concentrations of each compound were prepared from stock solutions diluted with SAM as necessary.

2.3. PFAS Exposure

Environmental standards of PFAS vary around the world. In Australia, the Heads of the Environmental Protection Agency for Australia and New Zealand (HEPA) have set PFAS guideline values based on the precautionary principle. The concentration levels used in this study were based on the PFAS National Environmental Management Plan (NEMP) (Version 2.0, National Chemicals Working Group of the Heads of EPAs Australia and New Zealand, Department of Climate Change, Energy, the Environment and Water, Australian Government), created by HEPA, which suggests a guideline value of 0.13 µg/L of PFOS to protect 95% aquatic species [2]. This was used as the basis of concentrations for all compounds. We exposed amphipods to five concentrations, both above and below the 0.13 µg/L limit (Table 1). This ensured that a) we could be sure of seeing an effect and b) that we could assess differences between high and low exposure levels. Variations in concentrations for each compound were due to small differences in the initial standard concentrations when purchased and the variation introduced during preparation of standards.
For each study, two days prior to the test, 14-day-old (juvenile) amphipods were collected into 2 L beakers until the experiment from the in-house culture tanks, using a 250 µm sieve. On day one of the test, 20 juvenile amphipods were randomly placed in each 600 mL beaker containing 400 mLs of SAM or PFAS treatment solution, including a substrate of 2 × 2 cm square of gauze cut from a sterilized bandage (Livingstone Triangular bandage—bleached calico) purchased from a local supermarket.
Beakers were randomly placed in an incubator with aeration maintained throughout the test. Amphipods were fed every second day with 1 mL powdered fish food (Tetramin®, Tetra Werke, Germany) solution (90 mg of ground Tetramin/50 mLs reversed osmosis water) and 1 mL YCT. The amphipods were exposed to a 16:8 h light: dark photoperiod and the temperature was maintained at 20 °C (±2) throughout the test. Water quality parameters (dissolved oxygen, pH, conductivity, and ammonia) were measured for each treatment at the conclusion of the experiment (Supplementary Materials—Table S1). A solvent control with laboratory-grade methanol at a concentration of 0.25 mL/L was run alongside the exposure treatments (no significant changes in survival or metabolites were detected from the solvent control). Amphipod survival was also recorded in each treatment (Table S2). On day 7, amphipods were counted from each beaker and immediately placed in microcentrifuge tubes on dry ice, to quench metabolites, and stored at −80 °C until analysis. Amphipod tissue was also collected for PFAS analysis to assess bioconcentration. Water samples were collected for PFAS measurement analysis.

2.4. Sample Preparation and Metabolite Extraction

A 300 µL aliquot of a methanol:chloroform (9:1) mixture was added to each sample. Amphipods were then homogenized by manual grinding using a metal spatula (which was cleaned between samples) within the microcentrifuge tubes and placed on dry ice. Samples were then sonicated for 15 min in an ultrasonic water bath and then centrifuged in an Eppendorf 5424 Centrifuge (Macquarie Park, NSW, Australia) at 5500 rpm for 15 min at room temperature. The supernatant was decanted into glass inserts and allowed to evaporate in a fume hood.
Dried extracts were derivatized using 20 µL of methoxyamine in pyridine solution (20 mg/mL). The samples were then vortexed for 30 s and left for 17 h at room temperature. Next, 20 µL of N-methyl-N-(trimethylsilyl)trifluoroacetamide (MSTFA) was added. Samples were vortexed again and left for one hour to react. Analytical-grade hexane was added to each sample (300 µL) prior to GC–MS analysis. A pooled biological quality control (PBQC) was prepared by pooling 7.5 µL of each sample extract, mixed thoroughly and aliquoted into 5 replicates. PBQCs were analysed along with the samples and controls/blanks to assess repeatability, instrument drift and quality control.

2.5. Gas Chromatography–Mass Spectrometry

Samples were analysed in random order. GC–MS analysis was performed using an Agilent 7890B gas chromatograph coupled to an Agilent 5977B mass spectrometer (Agilent, Santa Clara, CA, USA). The gas chromatograph was operated in the splitless mode, with a purge flow at 2 min. Separation was performed with an DB5–MS column (30 m × 250 μm, 0.25 μm) using helium as carrier gas at a flow rate of 1 mL per min. The injection volume was 1 μL and the injector temperature was 250 °C. The oven temperature rose from 35 to 300 °C at 25 °C/min, held at this temperature for 5 min, then increased to 310 °C (at 5 °C/min) and held for 5 min.
Mass spectra were recorded at 1.5 scans/s over an m/z range of 35–550. Metabolites were putatively annotated as level 2 compounds (based upon spectral similarity with public/commercial spectral libraries) where possible, using the National Institute of Standards and Technology (NIST MS search 2.0) library. Raw area retention times (RT) from each significant peak that occurred was recorded (Table S3). Environmental metabolomic reporting standards as described by Morrison et al. [16] were adhered to.

2.6. PFAS Determination Using Liquid Chromatography—Mass Spectrometry

PFAS was measured in water and amphipod tissue samples using an in-house method (ESA-P-ORG-16) based on Coggan et al. [17]. The methodology uses a triple quad LC–MS/MS Shimadzu 8060 with a C18 column (3 mm I.D. × 30 mm, P/N number 228-41606-91). The samples for biota were extracted using pure LC–MS-grade methanol and were tumbled for an hour at 30 rpm. The tumbled samples were then centrifuged for 20 min at 3000 rpm before they were passed through a 0.45 micron filter. Sample pH was then adjusted to 7 via the addition of acidified MilliQ water. Surrogates and internal standards for all compounds were added prior to analysis, as described in Coggan et al. [17].

2.7. Bioconcentration Factors (BCFs)

A reduction in PFAS concentration in the experiment could be due to compounds sticking to the experimental equipment such as the glass beakers, rather than the organisms. To account for this the concentration of PFAS in each experiment was measured both in the water and in the organisms at the end of the experiment via LC–MS/MS. This enabled us to confirm the actual concentration present, rather than just assuming the nominal experimental concentration was correct. We were also able to use these data to calculate the bioconcentration factor (BCF) of each compound.
BCF describes the readiness of chemicals to concentrate in organisms when the compounds are present in the environment. It is a required ecotoxicological parameter for chemical regulation and is generally taken to be the ratio of the concentration of a substance in a specific organism to the exposure concentration, at equilibrium in the media concerned. In this study, BCF was calculated via Equation (1):
BCF = C/Cw,
where BCF = bioconcentration factor, C = the concentration of the test substance in the organism in µg/g and Cw = the concentration of the substance in water in µg/L.

2.8. Statistical Analysis—Survival

The effect of PFAS exposure on amphipod survival was assessed using Student’s t-test (two samples assuming equal variances). Each treatment was compared to controls and significance was set at p < 0.05 (Table S4).

2.9. Metabolite Data Analysis

A combination of multivariate and univariate statistics were used for identification and selection of biomarkers using the total area of the chromatogram peak for each metabolite [18]. Each PFAS exposure was analysed separately.
Values were normalized to the median value of each treatment to account for differences between samples [19]. Values were then transformed (lognatural) to reduce variance between metabolite abundance to allow easier comparison between compounds of high and low concentrations.
Multivariate statistical analysis was applied as a first step for interrogating the data to observe trends, grouping and/or outliers [18]. Given the organisms used in this study were originally wild caught and then cultivated, rather than a lab-based population, their genetic and metabolic variation was expected to be higher than lab grown organisms. It was therefore possible that natural biological variation might mask the effects of PFAS exposure, particularly at low concentration levels. Orthogonal Partial Least Squares—Discriminant Analysis (OPLS-DA) with appropriate model overview and Variable Importance for Projection (VIP) investigating the metabolites causing the shift between groups (supplied in Supplementary Materials—Figures S1–S17) is a common method for determining separation in metabolomic datasets due to treatments when natural biological variation outweighs treatment variation [20]. Therefore, OPLS-DA was carried out using online MetaboAnalyst Version 5.0 (https://www.metaboanalyst.ca/MetaboAnalyst/) (accessed on 15 January 2021) [21] for each experiment separately.
A univariate Analysis Of VAriance (ANOVA) was also conducted to determine significant differences in metabolite levels between treatments (p < 0.05). Following this, Tukey’s post hoc test determined significance between treatments at each concentration. Metabolites were analysed using one-factor ANOVA with treatment as the main effect with R Version 3.0.3, R Foundation for Statistical Computing, Vienna, Austria [22]. Each p-value was adjusted for false discovery rate using the Benjamini–Hochberg (BH) method. The metabolic data were then graphed using box plots. Finally, pathway enrichment analysis was assessed using MetaboAnalyst, Version 5.0, (https://www.metaboanalyst.ca/MetaboAnalyst/) (accessed 15 January 2021) [21], to highlight potential pathways impacted by changes in metabolite abundance.

3. Results and Discussion

3.1. Survival

Survival varied across all treatments, with a decline in survival in the higher concentrations as shown in Figure 1. There was a decrease in total percent survival in high concentrations in all treatments, 21% (PFOS), 18% (GenX), and 12% (PFHxS), compared to the relevant controls. Only the high-dose group for PFOS had a significant reduction in survival compared to controls (p < 0.00) (Supplementary Materials, Table S4).
The survival rates of the controls and the medium exposure were similar, with the controls a little lower than the medium in GenX. However, there is no statistical difference in the survival rates of each group. The variation is higher than ideal, but when dealing with a wild-caught population, some variation is expected; indeed, from our experience, it is not uncommon in this kind of work. Furthermore, the medium dose of GenX was far lower than that has been shown to cause an effect on survival, so we did not expect it to cause an effect. Crucially, this observation backs up the main point that traditional ecotoxicology tests are not sensitive enough to register the effects of PFAS at environmentally relevant exposures levels. In this case, no statistical difference in the survival rates of each group and so testing for survival would show no effects of PFAS exposure. It takes the metabolomics approach to determine the subtle differences of sublethal exposure levels. Moreover, the current paper is not looking at survival but at the effects of PFAS on metabolism at sublethal exposure levels.
Additionally, at the end of the exposures, amphipods from the medium–high- and high-dose treatments for all compounds had visibly reduced biomass but this in itself was not enough to demonstrate an effect. This means that the Australian PFAS NEMP guideline for 80–99% protection of organisms at those concentrations is technically correct—if one is using survival as a metric. However, it is worth noting that the survival of the amphipods was affected following just 7 days of exposure and longer exposures may decrease the survival rates further. In comparison, the marine amphipod, Melita plumulosa, was described to have a PFOS 10-d LC50 of 3.3 mg/L and a reproduction EC50 of 2.2 mg/L [23]. These concentrations are significantly higher than those used in this study. Similarly, the common freshwater amphipod, Hyalella azteca, was reported to have a PFOS 42-d LC50 of 15 mg/L and an EC50 of 3.8 mg/L for reproduction [24]. This suggests amphipods can be considered tolerant of concentrations detected in the majority of PFAS-contaminated sites when using standard endpoints such as survival.
Boudreau et al. [25] stated that at the known environmental concentrations of PFAS, which range from ng/L to mg/L in ecosystems, there is seemingly no adverse risk of PFOS to freshwater ecosystems. This was after the authors had exposed a number of aquatic flora and invertebrates to short term 48 and 96 h PFOS and concluded that PFOS is highly toxic to freshwater systems at 100 mg/L [25]. These concentrations are far above the levels of PFAS in most systems.

3.2. Metabolomics

A total of 27 metabolites were detected using GC–MS. These included amino acids (e.g., valine; proline), carbohydrates (e.g., glycerol), carboxylic acids (e.g., acetic acid; propanoic acid), lipids (e.g., cholesterol), fatty acids (e.g., palmitic acid) and glycerides (e.g., monostearin). Only metabolite features that were statistically significant between treatments were fully identified and used for the multi- and univariate analysis.

3.3. PFOS Exposure

3.3.1. Multivariate Data Analysis for PFOS Exposure

OPLS-PA provided a visual tool to assess the effects of PFOS on the amphipods at different concentrations. There was an overlap of the metabolomic profile responses in amphipods exposed to the low-dose and control groups (Figure 2a) but separation between controls and all other groups (Figure 2b–e).
The Variable Importance for Projection (VIP) indicated similar metabolites to the univariate analysis were responsible for driving the main separation (see Section 3.3.2 and Supplementary Materials, Figures S1–S5). Briefly, isopropyl propane, valine and propanoic acid were shown to be the cause of most of the separation in the lower concentration of PFOS, whereas palmitic acid was highlighted in the higher concentration.

3.3.2. Univariate Data Analysis for PFOS Exposure

Univariate analysis of significant metabolite features was used to further elucidate the overall effects of exposure to PFOS on amphipods. There were seven metabolites whose levels changed significantly (Table 2).
Box plots were used to display significantly responding metabolite abundance trends across the different concentrations of PFOS. There was an increase in monostearin total abundance compared to the controls in the low, low–medium and medium–high groups (Figure 3a), but a statistically significant response was only detected in the low–medium (p = 0.02) and medium–high groups (p = 0.01; Table 2). Ketobutyric acid had the highest, and statistically significant (p = 0.00; Table 2), increase in abundance in the low concentration group compared to the controls. There was also a trending increase in the abundance of this metabolite compared to the control in all groups (Figure 3b). Palmitic acid showed a statistically significant response (p = 0.00; Table 2) at the highest concentration (Figure 3c). Glycerol decreased in abundance in the low–medium group then increased in the M–H- and high-dose groups (Figure 3d). The abundance of propanoic acid, acetic acid and valine significantly changed in response to PFOS (Table 2), and these box plots are listed in Supplementary Materials (Figure S6).

3.4. GenX Exposure

3.4.1. Multivariate Data Analysis for GenX Exposure

The amphipod metabolomic profiles in response to GenX had a distinct response compared to controls in each concentration (Figure 4a–e). There was a clear separation of groups from low to high concentrations.
Cholesterol, monostearin and glycerol were identified using the VIP for each concentration to be responsible for the separation between control and treatment in the OPLS-DA plots (Supplementary Materials, Figures S7–S11), these were later confirmed to be significantly changed via univariate analysis.

3.4.2. Univariate Data Analysis for GenX Exposure

An ANOVA was run for each metabolite following exposure to GenX to identify any significant responses (Table 3).
The metabolites showing the most change in each group were shown using box plots (Figure 5a–d).
Valine was significantly reduced in abundance only in the low concentration group compared to controls (Figure 5a) (p = 0.00; Table 3). Glycerol was significantly reduced in abundance following exposure to each concentration of GenX (Table 3) compared to the control (Figure 5b). Propanoic acid had similar trends to monostearin, with the abundance of both metabolites generally decreasing compared to controls (Figure 5c). The effect of GenX on monostearin started at the lower concentrations, decreased in the low–medium concentrations, increased at medium concentrations before further declining in the high treatments (Figure 5d). Overall, the abundance of these metabolites reduced in abundance compared to controls. Additionally, cholesterol and palmitic acid were found to significantly respond (Table 3), these box plots are located supplementary material (Figure S12a,b).

3.5. PFHxS Exposure

3.5.1. Multivariate Data Analysis for PFHxS Exposure

OPLS-DA was used to visualise the effect of the different concentrations of PFHxS on the amphipod metabolomic profiles. Similar to the response seen for GenX, there was clear separation in metabolite abundance in amphipods exposed to low, L–M, medium, M–H and high PFHxS concentrations (Figure 6a–e). These results indicate that both the newer compounds (GenX and PFHxS) have a stronger effect on metabolism than PFOS even though the latter is currently thought to be of more concern.
The VIP for the OPLSDA highlighted similar metabolites that were identified for driving the separation in PFOS and GenX OPLSDA (Supplementary Materials, Figures S13–S17). Cholesterol, monostearin and propanoic acid were identified to be responsible for the separation between PFHxS and control treatments. These metabolites were then identified to be significant in the univariate analysis.

3.5.2. Univariate Data Analysis for PFHxS Exposure

There were four metabolites that were identified as being statistically different in total abundance following exposure to PFHxS compared to controls (Table 4). Box plots were used to visualise the overall trend in the abundance of these four metabolites (Figure 7a–d).
Propanoic acid increased in abundance at the low concentration and spiked at the low–medium concentration (p = 0.04, Table 4) for PFHxS before dropping in total abundance at higher concentrations compared to the controls. This resulted in a significant reduction in abundance at a high concentration (p = 0.00, Table 4) (Figure 7a). There was a significant decline in abundance in the glyceride monostearin in the medium and high concentrations of PFHxS (Table 4), whereas there was little effect on monostearin abundance in the remaining concentrations, low, L–M and M–H (Figure 7b). PFHxS had a significant effect on the level of cholesterol, which was reduced in abundance at the low and high concentrations (p = 0.03, p = 0.01, respectively) (Table 4). However, cholesterol remained relatively consistent with control levels at the other concentrations (Figure 7c). Finally, the amino acid proline significantly decreased in abundance compared to the controls in three of the five concentrations of PFHxS (Table 4 and Figure 7d).

3.6. Overall PFAS Effect on Metabolites

Using metabolomics, we can begin to monitor, and rapidly screen, the small biological molecules in an organism at a given point in time and use these to help to elucidate the biological effects of PFAS [4,26]. The univariate analysis for each PFAS exposure highlighted significantly responding metabolites at different concentrations. There are several metabolites that were found to have a significant response in two and three of the PFAS exposures. For example, monostearin, palmitic acid and propanoic acid had a similar response in each PFAS exposure. It is however worth noting that the responses of the metabolites did not line up perfectly for each compound and PFAS concentration. For example, an increase in propanoic acid was detected in the low PFOS exposure, in the L–M group for PFHxS and the medium concentration of GenX. PFOS and GenX had similar effects on palmitic acid at each concentration, slightly dropping in abundance in the first three concentrations but increasing at high concentrations. This indicates that there may be different biological responses to low-dose PFAS exposure compared to high-dose exposure. The metabolites that change are all involved with fatty acid metabolism, potentially indicating an identifiable biological effect of exposure to fluorinated compounds. The result of fatty acid metabolism being impacted by PFAS exposure is consistent with the recent review investigating common biochemical pathways in response to PFAS [5]. This review highlighted those effects on fatty acid metabolism were seen in a range of organisms following PFAS exposure whether assessed by transcriptomics, lipidomics, or metabolomics [5].

3.7. Metabolomic Pathways

The biochemical pathways potentially disrupted due to exposure to the PFAS treatments were amino sugar metabolism; aspartate metabolism; pyruvate metabolism (e.g., acetic acid p = 0.04); steroid biosynthesis; steroidogenesis; bile acid biosynthesis (e.g., cholesterol, p = 0.00; palmitic acid, p = 0.00); glycerolipid metabolism (e.g., glycerol, p = 0.00; palmitic acid, p = 0.00); galactose metabolism (e.g., glycerol, p = 0.00); fatty acid metabolism; fatty acid elongation (e.g., palmitic acid, p = 0.00, monostearin, p = 0.00); propanoate metabolism (e.g., ketobutyric acid, p = 0.00; propanoic acid, p = 0.03); valine, leucine and isoleucine degradation (e.g., valine, p = 0.05); arginine and proline metabolism (e.g., Proline, p = 0.00).
Disruption in the pathways outlined in Figure 8 can impact the organism’s energy balance but this would likely be a general toxicant response and not a response to PFAS specifically. Fatty acid metabolism was affected, affecting the maintenance and development of the cuticle, leading to a delay in molts, growth, and development [27,28]. There is a possibility that reduced growth and size could lead to a decline in total population fitness over time, but more work would be needed to assess this.
There was a change in the metabolite profile of each treatment at every concentration of PFAS compared to controls. This response was dose and contaminant dependent. Therefore, a linear response to PFAS cannot be assumed when organisms are exposed to increasing concentrations of these compounds. Studies investigating the effect of PFAS exposure on humans have made similar observations, where PFAS had a non-monotonic response in relation to kidney disease, cognitive function, diabetes, and fertility [29,30,31,32]. The results shown here challenge the common perception that increasing concentration equals increasing toxicity and accumulation of PFAS in a linear fashion. Instead, the data suggest metabolic responses can vary depending on the level of exposure.
Exposure to PFAS has been shown to increase fatty acid and energy metabolism in a range of aquatic organisms—for example, the water flea Daphnia magna [33] and the eastern school prawn Metapenaeus macleayi [34]. A study of turtles from PFAS-impacted waterways in Queensland, Australia, found an inflammation response and re-routing of central carbon metabolites [4]. With many studies being published, similar metabolites/pathways are emerging across organisms. In a recent review, Beale et al. [5], listed the common biochemical pathways affected following PFAS exposure that were also detected in this study. These alterations could result in metabolic deficiencies or perturbations that might lead to a decline in physical condition, fertility, and survival.

3.8. PFAS Water Concentration

A sample from the high-dose group from each treatment was measured via LC–MS (Table 5). A large reduction in PFOS concentration compared to the nominal concentration was observed, with a slightly smaller reduction for PFHxS. In contrast, there was little reduction in GenX levels in the water between the start and end of the experiment. This could be due to GenX binding with the equipment used, such as the beakers or substrate, rather then accumulating in the amphipods but steps were taken to reduce this during sample preparation and analysis, as per Coggan et al. [17].

3.9. Bioconcentration Factors

PFAS concentrations were below the limit of detection in the amphipods measured from low, low–med and medium concentrations. The BCFs at a high concentration of each compound are shown in Table 6. Only PFOS was found to bioaccumulate—GenX was quickly eliminated from the amphipods, and PFHxS was eliminated but at a slower rate than GenX. Exposure to PFHxS and GenX still induced metabolic changes even though BCFs were low. Only four of the eight metabolites detected to significantly change in response to PFAS did so in correlation with bioaccumulation. This would seem to indicate that bioaccumulation is not necessarily needed for a metabolic response to occur in an organism due to PFAS exposure.
PFOS has been reported to have high bioaccumulation potential in invertebrates, fish, reptiles, birds and mammals [4,13,35,36,37,38]. This was also found in the current study, where the amphipods had accumulated a greater amount of PFOS compared to GenX and PFHxS following the 7-day exposure. It is however important to note that studies on the bioaccumulation of PFAS are variable and not necessarily reported in a standardised fashion (e.g., dry weight versus wet weight), which limits their effectiveness in some cases. McCarthy et al. [39] also highlight the limited data that have been collected for fish, benthic organisms, and terrestrial biota (including reptiles and amphibians) and state that toxicology testing across all classes of organisms has not been consistent [39]. Additionally, Taylor et al. [34] noted that differences among individuals, species and environments can lead to a highly variable bioaccumulation of PFAS in aquatic systems. Despite this, it is still generally acknowledged that toxicity must come with bioaccumulation [40]. For example, Karnjanapiboonwong et al. [41] investigated the response of earthworms (Eisenia fetida) to a range of PFAS. They found that toxicity and accumulation of PFAS increased with increasing chain length, and similar results were reported by Taylor et al. [34].
Chain length is thought to also affect how and which PFAS are metabolised and excreted from organisms [33,37,42,43,44,45,46]. For example, GenX (C6) did not appear to have bioaccumulation potential when exposed to small fish, and it was hypothesised to be eliminated after ingestion in comparison to PFOS (C8) [13]. This result from Hassell et al. [13] is consistent with the GenX and PFOS exposure seen in this study using amphipods. PFOS accumulated in the fish and the amphipods, whereas GenX had the lowest BCF in the amphipods, either due to lack of uptake or rapid elimination. GenX is a relatively new/emerging PFAS, and little research has been carried out on its bioaccumulation potential in organisms. This should be addressed in future studies given these chemicals are expected to accumulate in ecosystems over a period of time [11,47].
PFHxS is also a C6 fluoro-carbon compound and, as such, it is assumed to have a similarly low bioaccumulation potential. In this study, the BCF of PFHxS was in between that of PFOS and GenX, in agreement with studies claiming BCF increases with increasing chain length [48]. However, chain length may not be the only factor determining the BCF of these compounds. It is possible that the shape and overall width of the molecule may also need to be considered.

4. Conclusions

This is the first study investigating the effects of PFAS using metabolomics on the freshwater amphipod Austrochiltonia. subtenuis. The work shows that metabolomics can add value to ecotoxicology studies on PFAS. Metabolic effects were observed at lower concentrations of PFAS than were seen with more traditional ecotoxicological endpoints such as survival or bioaccumulation, indicating that our understanding of the effects of these compounds is not complete. Metabolomic pathways detected to be affected by PFAS exposure have been shown to respond in other biota exposed to PFAS. These pathways include fatty acid, arginine, proline, glycolipid, galactose, and aspartate metabolism. We recommend additional research on the metabolome of invertebrates to achieve a detailed understanding of the entire consequence of these changes. Here, we recommend further investigations using omic techniques to better understand the sublethal effects of PFAS on organisms at environmentally relevant concentrations—specifically investigating the effect on a variety of organisms and investigating the bioaccumulation of PFAS detected in ecosystems as well as investigating whether temperature, salinity and other parameters impact the metabolomic responses seen to be affected by PFAS and bioaccumulation. With a larger variety of studies, confidence in the sublethal effects of PFAS will continue to emerge. This study increases our overall understanding of PFAS toxicity at environmentally relevant concentrations. With further investigations, metabolite responses could potentially help PFAS-monitoring remediation and efforts to be better targeted.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/metabo12111135/s1, Table S1: Water quality for amphipod PFAS exposure (PFOS, GenX and PFHxS). Table S2: Total survival of amphipods following PFAS exposure (PFOS, GenX and PFHxS). Table S3: Retention time and m/z of significant metabolites detected. Table S4: Survival t-test (two samples assuming equal variance). Figures S1–S5: Orthogonal Partial Lease Squares—Discriminant Analysis (OPLS-DA)—amphipod separation PFOS model overview and VIP. Figure S6a–c: Significant metabolite features from amphipods when exposed to PFOS at low (0.04 µg/L), low–medium (0.20 µg/L), medium (1.00 µg/L), medium–high (5.00 µg/L) and high (25.00 µg/L) concentrations. (a) Propanoic acid; b) acetic acid c) valine. Figure S7–S11: Orthogonal Partial Lease Squares—Discriminant Analysis (OPLS-DA)—amphipod separation to GenX model overview and VIP. Figure S12a,b: Significant metabolites from amphipods exposed to GenX at low (0.03 µg/L), low–medium (0.16 µg/L), medium (0.80 µg/L), medium–high (4.00 µg/L) and high (20.00 µg/L) concentrations. (a) Cholesterol and (b) palmitic acid. Figures S13–S17: Orthogonal Partial Lease Squares—Discriminant Analysis (OPLS-DA)—amphipod separation following low concentration to PFHxS model overview and VIP.

Author Contributions

Conceptualization, G.M.S., O.A.H.J. and S.M.L.; methodology, G.M.S., S.M.L. and O.A.H.J.; formal analysis, G.M.S., N.S., T.L.C. and M.P.J.A. writing—original draft preparation, G.M.S.; writing—review and editing, S.M.L., O.A.H.J., N.S., T.L.C. and M.P.J.A. All authors have read and agreed to the published version of the manuscript.

Funding

The lead author thanks the Australian government for the award of an Australian Government Research Training Program (RTP) Scholarship administered by RMIT University Research Stipend.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in article and Supplementary Materials.

Acknowledgments

All authors thank Rebecca Reid and Hung Vu for maintenance and management of the amphipod cultures and the experimental set up.

Conflicts of Interest

The authors declare no conflict of interest, there is no connection between RMIT University and the subject of this manuscript.

References

  1. Pelch, K.E.; Reade, A.; Wolffe, T.A.M.; Kwiatkowski, C.F. PFAS Health Effects Database: Protocol for a Systematic Evidence Map. Environ. Int. 2019, 130, 104851. [Google Scholar] [CrossRef] [PubMed]
  2. Heads of Environmental Protection Agency Australia and New Zealand (HEPA) PFAS National Environmental Management Plan. Available online: https://www.dcceew.gov.au/sites/default/files/documents/pfas-nemp-2.pdf (accessed on 8 August 2019).
  3. Sinclair, G.M.; Long, S.M.; Jones, O.A.H. What Are the Effects of PFAS Exposure at Environmentally Relevant Concentrations? Chemosphere 2020, 258, 127340. [Google Scholar] [CrossRef] [PubMed]
  4. Beale, D.J.; Hillyer, K.; Nilsson, S.; Limpus, D.; Bose, U.; Broadbent, J.A.; Vardy, S. Bioaccumulation and Metabolic Response of PFAS Mixtures in Wild-Caught Freshwater Turtles (Emydura macquarii macquarii) Using Omics-Based Ecosurveillance Techniques. Sci. Total Environ. 2021, 806, 151264. [Google Scholar] [CrossRef]
  5. Beale, D.J.; Sinclair, G.M.; Shah, R.; Paten, A.M.; Kumar, A.; Long, S.M.; Vardy, S.; Jones, O.A.H. A Review of Omics-Based PFAS Exposure Studies Reveals Common Biochemical Response Pathways. Sci. Total Environ. 2022, 845, 157255. [Google Scholar] [CrossRef] [PubMed]
  6. Jones, O.A.H.; Maguire, M.L.; Griffin, J.L.; Dias, D.A.; Spurgeon, D.J.; Svendsen, C. Metabolomics and Its Use in Ecology. Austral Ecol. 2013, 38, 713–720. [Google Scholar] [CrossRef]
  7. Buckley, N.; Sim, M.; Douglas, K.; Håkansson, H. Expert Health Panel for Per- and Poly-Fluoroalkyl Substances (PFAS). Available online: https://www.health.gov.au/internet/main/publishing.nsf/Content/C9734ED6BE238EC0CA2581BD00052C03/$File/expert-panel-report.pdf (accessed on 30 May 2019).
  8. Jones, O.A.H. Illuminating the Dark Metabolome to Advance the Molecular Characterisation of Biological Systems. Metabolomics 2018, 14, 101. [Google Scholar] [CrossRef]
  9. Yao, X.; Cao, D.; Wang, F.; Zhang, W.; Ma, C.; Song, M. An Overview of Omics Approaches to Characterize the Effect of Perfluoroalkyl Substances in Environmental Health. TrAC—Trends Anal. Chem. 2019, 121, 115367. [Google Scholar] [CrossRef]
  10. Hogue, C. Governments Agree to Eliminate PFHxS. C&EN Glob. Enterp. 2022, 100, 14. [Google Scholar]
  11. Brandsma, S.H.; Koekkoek, J.C.; van Velzen, M.J.M.; de Boer, J. The PFOA Substitute GenX Detected in the Environment near a Fluoropolymer Manufacturing Plant in the Netherlands. Chemosphere 2019, 220, 493–500. [Google Scholar] [CrossRef]
  12. Xiao, F. Emerging Poly- and Perfluoroalkyl Substances in the Aquatic Environment: A Review of Current Literature. Water Res. 2017, 124, 482–495. [Google Scholar] [CrossRef]
  13. Hassell, K.L.; Coggan, T.L.; Cresswell, T.; Kolobaric, A.; Berry, K.; Crosbie, N.D.; Blackbeard, J.; Pettigrove, V.J.; Clarke, B.O. Dietary Uptake and Depuration Kinetics of Perfluorooctane Sulfonate, Perfluorooctanoic Acid, and Hexafluoropropylene Oxide Dimer Acid (GenX) in a Benthic Fish. Environ. Toxicol. Chem. 2020, 39, 595–603. [Google Scholar] [CrossRef]
  14. Jeppe, K.J.; Kellar, C.R.; Marshall, S.; Colombo, V.; Sinclair, G.M.; Pettigrove, V. Bifenthrin Causes Toxicity in Urban Stormwater Wetlands: Field and Laboratory Assessment Using Austrochiltonia (Amphipoda). Environ. Sci. Technol. 2017, 51, 7254–7262. [Google Scholar] [CrossRef] [PubMed]
  15. Borgmann, U. Systematic Analysis of Aqueous Ion Requirements of Hyalella Azteca: A Standard Artificial Medium Including the Essential Bromide Ion. Arch. Environ. Contam. Toxicol 1996, 188, 356–363. [Google Scholar] [CrossRef]
  16. Morrison, N.; Dan, A.E.; Ae, B.; Bundy, J.G.; Tim, A.E.; Ae, C.; Currie, F.; Matthew, A.E.; Davey, P.; Nathan, A.E.; et al. Standard Reporting Requirements for Biological Samples in Metabolomics Experiments: Environmental Context. Metabolomics 2007, 3, 203–210. [Google Scholar] [CrossRef]
  17. Coggan, T.L.; Anumol, T.; Pyke, J.; Shimeta, J.; Clarke, B.O. A Single Analytical Method for the Determination of 53 Legacy and Emerging Per- and Polyfluoroalkyl Substances (PFAS) in Aqueous Matrices. Anal. Bioanal. Chem. 2019, 411, 3507–3520. [Google Scholar] [CrossRef]
  18. Lacalle-Bergeron, L.; Izquierdo-Sandoval, D.; Sancho, J.V.; López, F.J.; Hernández, F.; Portolés, T. Chromatography Hyphenated To High Resolution Mass Spectrometry In Untargeted Metabolomics For Investigation Of Food (BIO)Markers. TrAC Trends Anal. Chem. 2020, 135, 116161. [Google Scholar] [CrossRef]
  19. Sinclair, G.M.; O’Brien, A.L.; Keough, M.; De Souza, D.P.; Dayalan, S.; Kanojia, K.; Kouremenos, K.; Tull, D.L.; Coleman, R.A.; Jones, O.A.H.; et al. Using Metabolomics to Assess the Sub-Lethal Effects of Zinc and Boscalid on an Estuarine Polychaete Worm over Time. Metabolomics 2019, 15, 108. [Google Scholar] [CrossRef]
  20. Blasco, H.; Błaszczyński, J.; Billaut, J.C.; Nadal-Desbarats, L.; Pradat, P.F.; Devos, D.; Moreau, C.; Andres, C.R.; Emond, P.; Corcia, P.; et al. Comparative Analysis of Targeted Metabolomics: Dominance-Based Rough Set Approach versus Orthogonal Partial Least Square-Discriminant Analysis. J. Biomed. Inform. 2015, 53, 291–299. [Google Scholar] [CrossRef]
  21. Pang, Z.; Chong, J.; Zhou, G.; De Lima Morais, D.A.; Chang, L.; Barrette, M.; Gauthier, C.; Jacques, P.É.; Li, S.; Xia, J. MetaboAnalyst 5.0: Narrowing the Gap between Raw Spectra and Functional Insights. Nucleic Acids Res. 2021, 49, W388–W396. [Google Scholar] [CrossRef]
  22. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020; Available online: https://www.r-project.or (accessed on 15 January 2021).
  23. Simpson, S.L.; Liu, Y.; Spadaro, D.A.; Wang, X.; Kookana, R.S.; Batley, G.E. Chronic Effects and Thresholds for Estuarine and Marine Benthic Organism Exposure to Perfluorooctane Sulfonic Acid (PFOS )–Contaminated Sediments: Influence of Organic Carbon and Exposure Routes. Sci. Total Environ. 2021, 776, 146008. [Google Scholar] [CrossRef]
  24. Krupa, P.M.; Lotufo, G.R.; Mylroie, E.J.; May, L.K.; Gust, K.A.; Kimble, A.N.; Jung, M.G.; Boyda, J.A.; Moore, W. Chronic Aquatic Toxicity of Perfluorooctane Sulfonic Acid (PFOS) to Ceriodaphnia Dubia, Chironomus Dilutus, Danio Rerio, and Hyalella Azteca. Ecotoxicol. Environ. Saf. 2022, 241, 113838. [Google Scholar] [CrossRef] [PubMed]
  25. Boudreau, T.M.; Sibley, P.K.; Mabury, S.A.; Muir, D.G.C.; Solomon, K.R. Laboratory Evaluation of the Toxicity of Perfluorooctane Sulfonate (PFOS) on Selenastrum Capricornutum, Chlorella Vulgaris, Lemna Gibba, Daphnia Magna, and Daphnia Pulicaria. Arch. Environ. Contam. Toxicol. 2003, 44, 307–313. [Google Scholar] [CrossRef] [PubMed]
  26. Ankley, G.T.; Cureton, P.; Hoke, R.A.; Houde, M.; Kumar, A.; Kurias, J.; Lanno, R.; McCarthy, C.; Newsted, J.; Salice, C.J.; et al. Assessing the Ecological Risks of Per- and Polyfluoroalkyl Substances: Current State-of-the Science and a Proposed Path Forward. Environ. Toxicol. Chem. 2020, 40, 564–605. [Google Scholar] [CrossRef]
  27. Lithodidae, D.A.; Beder, A.M.; Copeman, L.A.; Eckert, G.L. Crustacean Biology. J. Crustac. Biol. J. 2018, 38, 728–738. [Google Scholar] [CrossRef] [Green Version]
  28. Li, W.; Li, S.; Liu, J.; Wang, X.; Chen, H.; Hao, H.; Wang, K. Vital Carbohydrate and Lipid Metabolites in Serum Involved in Energy Metabolism during Pubertal Molt of Mud Crab (Scylla paramamosain). Metabolites 2021, 11, 651. [Google Scholar] [CrossRef]
  29. Kyun, S.; Ding, N.; Han, D. Perfluoroalkyl Substances and Cognitive Function in Older Adults: Should We Consider Non-Monotonic Dose-Responses and Chronic Kidney Disease? Environ. Res. 2021, 192, 110346. [Google Scholar] [CrossRef]
  30. Park, S.; Ding, N.; Han, D. Perfluoroalkyl Substances and Cognitive Function: Is the Association Non-Monotonic in Humans? Environ. Epidemiol. 2019, 3, 305. [Google Scholar]
  31. Lind, L.; Zethelius, B.; Salihovic, S.; Lind, M. Circulating Levels of Perfluoroalkyl Substances and Prevalent Diabetes in the Elderly. Diabetologia 2014, 473–479. [Google Scholar] [CrossRef] [PubMed]
  32. Tarapore, P.; Ouyang, B. Perfluoroalkyl Chemicals and Male Reproductive Health: Do PFOA and PFOS Increase Risk for Male Infertility? Int. J. Environ. Res. Public Health 2021, 18, 3794. [Google Scholar] [CrossRef] [PubMed]
  33. Kariuki, M.N.; Nagato, E.G.; Lankadurai, B.P.; Simpson, A.J.; Simpson, M.J. Analysis of Sub-Lethal Toxicity of Perfluorooctane Sulfonate (PFOS) to Daphnia Magna Using1H Nuclear Magnetic Resonance-Based Metabolomics. Metabolites 2017, 7, 15. [Google Scholar] [CrossRef] [Green Version]
  34. Taylor, M.D.; Bräunig, J.; Mueller, J.F.; Crompton, M.; Dunstan, R.H.; Nilsson, S. Metabolomic Profiles Associated with Exposure to Per- And Polyfluoroalkyl Substances (PFASs) in Aquatic Environments. Environ. Sci. Process. Impacts 2019, 21, 1980–1990. [Google Scholar] [CrossRef] [PubMed]
  35. Ahrens, L.; Bundschuh, M. Fate and Effects of Poly- and Perfluoroalkyl Substances in the Aquatic Environment: A Review. Environ. Toxicol. Chem. 2014, 33, 1921–1929. [Google Scholar] [CrossRef] [PubMed]
  36. Ghisi, R.; Vamerali, T.; Manzetti, S. Accumulation of Perfluorinated Alkyl Substances (PFAS) in Agricultural Plants: A Review. Environ. Res. 2019, 169, 326–341. [Google Scholar] [CrossRef] [PubMed]
  37. Palmer, K.; Bangma, J.T.; Reiner, J.L.; Bonde, R.K.; Korte, J.E.; Boggs, A.S.P.; Bowden, J.A. Per- and Polyfluoroalkyl Substances (PFAS) in Plasma of the West Indian Manatee (Trichechus manatus). Mar. Pollut. Bull. 2019, 140, 610–615. [Google Scholar] [CrossRef] [PubMed]
  38. Wen, B.; Zhang, H.; Li, L.; Hu, X.; Liu, Y.; Shan, X.; Zhang, S. Bioavailability of Perfluorooctane Sulfonate (PFOS) and Perfluorooctanoic Acid (PFOA) in Biosolids-Amended Soils to Earthworms (Eisenia fetida). Chemosphere 2015, 118, 361–366. [Google Scholar] [CrossRef] [PubMed]
  39. McCarthy, C.; Kappleman, W.; DiGuiseppi, W. Ecological Considerations of Per- and Polyfluoroalkyl Substances (PFAS). Curr. Pollut. Rep. 2017, 3, 289–301. [Google Scholar] [CrossRef]
  40. Wang, W.-X. Bioaccumulation and Biomonitoring. In Marine Ecotoxicology; Academic Press: Cambridge, MA, USA, 2016; pp. 99–119. [Google Scholar] [CrossRef]
  41. Karnjanapiboonwong, A.; Deb, S.K.; Subbiah, S.; Wang, D.; Anderson, T.A. Perfluoroalkylsulfonic and Carboxylic Acids in Earthworms (Eisenia fetida): Accumulation and Effects Results from Spiked Soils at PFAS Concentrations Bracketing Environmental Relevance. Chemosphere 2018, 199, 168–173. [Google Scholar] [CrossRef]
  42. Chiesa, L.M.; Nobile, M.; Pasquale, E.; Balzaretti, C.; Cagnardi, P.; Tedesco, D.; Panseri, S.; Arioli, F. Detection of Perfluoroalkyl Acids and Sulphonates in Italian Eel Samples by HPLC-HRMS Orbitrap. Chemosphere 2018, 193, 358–364. [Google Scholar] [CrossRef]
  43. D’eon, J.C.; Mabury, S.A. Exploring Indirect Sources of Human Exposure to Perfluoroalkyl Carboxylates (PFCAs): Evaluating Uptake, Elimination, and Biotransformation of Polyfluoroalkyl Phosphate Esters (PAPs) in the Rat. Environ. Health Perspect. 2011, 119, 344–350. [Google Scholar] [CrossRef]
  44. Land, M.; De Wit, C.A.; Bignert, A.; Cousins, I.T.; Herzke, D.; Johansson, J.H.; Martin, J.W. What Is the Effect of Phasing out Long-Chain per- and Polyfluoroalkyl Substances on the Concentrations of Perfluoroalkyl Acids and Their Precursors in the Environment? A Systematic Review. Environ. Evid. 2018, 7, 4. [Google Scholar] [CrossRef] [Green Version]
  45. Whitacre, D.M. Reviews of Environmental Contamination and Toxicology. Rev. Environ. Contam. Toxicol. 2015, 233, 107–130. [Google Scholar] [CrossRef]
  46. Zhao, S.; Zhu, L.; Liu, L.; Liu, Z.; Zhang, Y. Bioaccumulation of Perfluoroalkyl Carboxylates (PFCAs) and Perfluoroalkane Sulfonates (PFSAs) by Earthworms (Eisenia fetida) in Soil. Environ. Pollut. 2013, 179, 45–52. [Google Scholar] [CrossRef] [PubMed]
  47. Gomis, M.I.; Vestergren, R.; Borg, D.; Cousins, I.T. Comparing the Toxic Potency in Vivo of Long-Chain Perfluoroalkyl Acids and Fluorinated Alternatives. Environ. Int. 2018, 113, 1–9. [Google Scholar] [CrossRef] [PubMed]
  48. Menger, F.; Pohl, J.; Ahrens, L.; Carlsson, G.; Örn, S. Behavioural Effects and Bioconcentration of Per- and Polyfluoroalkyl Substances (PFASs) in Zebrafish (Danio rerio) Embryos. Chemosphere 2020, 245, 125573. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Mean (±SD) survival of amphipods Austrochiltona subtenuis following exposure to different concentrations of PFOS, PFHxS and GenX for 7 days. Two-tailed t-test between control and concentration revealed survival of amphipods were significantly reduced (p = 0.00) * at the highest PFOS concentration. Ctl = control, LM = low–medium, and MH = medium–high.
Figure 1. Mean (±SD) survival of amphipods Austrochiltona subtenuis following exposure to different concentrations of PFOS, PFHxS and GenX for 7 days. Two-tailed t-test between control and concentration revealed survival of amphipods were significantly reduced (p = 0.00) * at the highest PFOS concentration. Ctl = control, LM = low–medium, and MH = medium–high.
Metabolites 12 01135 g001
Figure 2. Orthogonal PLS-DA of polar metabolites from amphipods exposed to PFOS (blue) and control (pink) for 7 days at (a) low (0.04 µg/L), (b) low–medium (0.20 µg/L), (c) medium (1.00 µg/L), (d) medium–high (5.00 µg/L) and (e) high (25.00 µg/L) concentrations. Orthogonal PLS-DA model validation data can be found in Supplementary Materials (Figures S1–S5).
Figure 2. Orthogonal PLS-DA of polar metabolites from amphipods exposed to PFOS (blue) and control (pink) for 7 days at (a) low (0.04 µg/L), (b) low–medium (0.20 µg/L), (c) medium (1.00 µg/L), (d) medium–high (5.00 µg/L) and (e) high (25.00 µg/L) concentrations. Orthogonal PLS-DA model validation data can be found in Supplementary Materials (Figures S1–S5).
Metabolites 12 01135 g002
Figure 3. Top four significant metabolite features from amphipods when exposed to PFOS at low (0.04 µg/L), low–medium (0.20 µg/L), medium (1.00 µg/L), medium–high (5.00 µg/L) and high (25.00 µg/L) concentrations. (a) Monostearin; (b) ketobutyric acid; (c) palmitic acid; and (d) glycerol.
Figure 3. Top four significant metabolite features from amphipods when exposed to PFOS at low (0.04 µg/L), low–medium (0.20 µg/L), medium (1.00 µg/L), medium–high (5.00 µg/L) and high (25.00 µg/L) concentrations. (a) Monostearin; (b) ketobutyric acid; (c) palmitic acid; and (d) glycerol.
Metabolites 12 01135 g003
Figure 4. Orthogonal PLS-DA: separation of polar metabolites from amphipods exposed to GenX for 7 days at (a) low (0.03 µg/L), (b) low–medium (0.16 µg/L), (c) medium (0.80 µg/L), (d) medium–high (4.00 µg/L) and (e) high (20.00 µg/L) concentrations. Orthogonal PLS-DA Model overview can be found in Supplementary Materials (Figures S7–S11).
Figure 4. Orthogonal PLS-DA: separation of polar metabolites from amphipods exposed to GenX for 7 days at (a) low (0.03 µg/L), (b) low–medium (0.16 µg/L), (c) medium (0.80 µg/L), (d) medium–high (4.00 µg/L) and (e) high (20.00 µg/L) concentrations. Orthogonal PLS-DA Model overview can be found in Supplementary Materials (Figures S7–S11).
Metabolites 12 01135 g004
Figure 5. Significant metabolites from amphipods exposed to GenX at low (0.03 µg/L), low–medium (0.16 µg/L), medium (0.80 µg/L), medium–high (4.00 µg/L) and high (20.00 µg/L) concentrations. (a) Valine; (b) glycerol; (c) propanoic acid; (d) monostearin.
Figure 5. Significant metabolites from amphipods exposed to GenX at low (0.03 µg/L), low–medium (0.16 µg/L), medium (0.80 µg/L), medium–high (4.00 µg/L) and high (20.00 µg/L) concentrations. (a) Valine; (b) glycerol; (c) propanoic acid; (d) monostearin.
Metabolites 12 01135 g005
Figure 6. Orthogonal PLS-DA: separation of polar metabolites from amphipods exposed to PFHxS for 7 days at (a) low (0.03 µg/L), (b) low–medium (0.15 µg/L), (c) medium (0.76 µg/L), (d) medium–high (3.78 µg/L) and (e) high (18.92 µg/L) concentrations. Orthogonal PLS-DA model overview can be found in Supplementary Materials (Figures S13–S17).
Figure 6. Orthogonal PLS-DA: separation of polar metabolites from amphipods exposed to PFHxS for 7 days at (a) low (0.03 µg/L), (b) low–medium (0.15 µg/L), (c) medium (0.76 µg/L), (d) medium–high (3.78 µg/L) and (e) high (18.92 µg/L) concentrations. Orthogonal PLS-DA model overview can be found in Supplementary Materials (Figures S13–S17).
Metabolites 12 01135 g006
Figure 7. (ad) Significant metabolites from amphipods exposed to PFHxS at low (0.03 µg/L), low–medium (0.15 µg/L), medium (0.76 µg/L), medium–high (3.78 µg/L) and high (18.92 µg/L) concentrations. (a) Propanoic acid; (b) monostearin; (c) cholesterol; (d) proline.
Figure 7. (ad) Significant metabolites from amphipods exposed to PFHxS at low (0.03 µg/L), low–medium (0.15 µg/L), medium (0.76 µg/L), medium–high (3.78 µg/L) and high (18.92 µg/L) concentrations. (a) Propanoic acid; (b) monostearin; (c) cholesterol; (d) proline.
Metabolites 12 01135 g007
Figure 8. Metabolic pathways affected by PFAS exposure following a 7-day exposure. Red outline indicates the metabolites that changed significantly in abundance (listed in Table 2, Table 3 and Table 4).
Figure 8. Metabolic pathways affected by PFAS exposure following a 7-day exposure. Red outline indicates the metabolites that changed significantly in abundance (listed in Table 2, Table 3 and Table 4).
Metabolites 12 01135 g008
Table 1. Nominal concentrations (µg/L) of PFOS, GenX and PFHxS used in three separate Austrochiltonia subtenuis 7-day exposures using analytical-grade PFAS standards.
Table 1. Nominal concentrations (µg/L) of PFOS, GenX and PFHxS used in three separate Austrochiltonia subtenuis 7-day exposures using analytical-grade PFAS standards.
TreatmentLowLow–MediumMediumMedium–HighHigh
Control00000
PFOS0.040.201.005.0025.00
GenX0.030.160.804.0020.00
PFHxS0.030.150.763.7818.92
Table 2. Significantly changed metabolites following exposure to PFOS identified from a one-way ANOVA. Values in bold represent statistical significance following Tukey’s post hoc test. All values listed are BH-adjusted p values (p < 0.05), using a 95% family-wise confidence level. Degrees of freedom = 24/5.
Table 2. Significantly changed metabolites following exposure to PFOS identified from a one-way ANOVA. Values in bold represent statistical significance following Tukey’s post hoc test. All values listed are BH-adjusted p values (p < 0.05), using a 95% family-wise confidence level. Degrees of freedom = 24/5.
MetaboliteMean Square
Residual
Comparison between Treatment and Control
LOW-CTRLL–M–CTRLMED–CTRLM–H–CTRLHIGH–CTRL
Monostearin3 × 10110.090.021.000.010.85
Ketobutyric acid2 × 10110.001.000.700.960.41
Palmitic acid3 × 10110.171.001.000.080.00
Glycerol1 × 10130.020.010.010.010.93
Propanoic acid2 × 10110.740.950.990.520.50
Acetic acid2 × 10141.000.111.000.990.94
Valine1 × 10130.200.060.720.120.08
Table 3. Significantly changed metabolites in amphipods following exposure to GenX. Values in bold represent statistical significance following Tukey’s post hoc test. All values listed are BH-adjusted p values (p < 0.05) using a 95% family-wise confidence level. Degrees of freedom = 21/5.
Table 3. Significantly changed metabolites in amphipods following exposure to GenX. Values in bold represent statistical significance following Tukey’s post hoc test. All values listed are BH-adjusted p values (p < 0.05) using a 95% family-wise confidence level. Degrees of freedom = 21/5.
MetaboliteMean Square
Residual
Comparison between Treatment and Control
LOW-CTRLL–M–CTRLMED–CTRLM–H–CTRLHIGH–CTRL
Valine9 × 10100.001.000.750.750.56
Glycerol9 × 10100.010.040.000.010.00
Monostearin9 × 10100.440.740.030.870.70
Propanopic acid9 × 10100.010.021.000.040.02
Cholesterol9 × 10100.270.440.000.230.10
Palmitic acid9 × 10100.650.480.190.120.02
Table 4. Significantly changed metabolites identified from a one-way ANOVA following exposure to PFHxS compared to controls. Values in bold represent statistical significance following Tukey’s post hoc test. All values listed are BH-adjusted p values (p < 0.05) using a 95% family-wise confidence level. Degrees of freedom = 20/5.
Table 4. Significantly changed metabolites identified from a one-way ANOVA following exposure to PFHxS compared to controls. Values in bold represent statistical significance following Tukey’s post hoc test. All values listed are BH-adjusted p values (p < 0.05) using a 95% family-wise confidence level. Degrees of freedom = 20/5.
MetaboliteMean Square ResidualComparison between Treatment and Control
LOW-CTRLL–M-CTRLMED-CTRLM–H-CTRLHIGH-CTRL
Propanoic acid0.100.860.041.001.000.00
Monostearin0.431.001.000.010.990.00
Cholesterol0.480.030.991.001.000.01
Proline0.310.020.020.100.000.25
Table 5. Nominal and Measured PFAS concentrations (µg/L)in treatments (PFOS, GenX and PFHxS). Control water < 0.01 µg/L of PFOS and PFHxS and < 0.1 µg/L of GenX (limit of detection).
Table 5. Nominal and Measured PFAS concentrations (µg/L)in treatments (PFOS, GenX and PFHxS). Control water < 0.01 µg/L of PFOS and PFHxS and < 0.1 µg/L of GenX (limit of detection).
TreatmentDay 0
(Nominal)
Day 7
(Measured)
PFOS25.000.11
GenX20.0019.71
PFHxS18.9212.3
Table 6. Bioconcentration factors for each treatment of PFAS at a high concentration (µg/L).
Table 6. Bioconcentration factors for each treatment of PFAS at a high concentration (µg/L).
TreatmentHigh R1High R2High R3High R4High R5
PFOS18.1827.27190.91118.189.09
GenX0.010.010.010.010.01
PFHxS0.050.070.040.080.03
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sinclair, G.M.; Long, S.M.; Singh, N.; Coggan, T.L.; Askeland, M.P.J.; Jones, O.A.H. Exposure to Environmentally Relevant Levels of PFAS Causes Metabolic Changes in the Freshwater Amphipod Austrochiltonia subtenuis. Metabolites 2022, 12, 1135. https://doi.org/10.3390/metabo12111135

AMA Style

Sinclair GM, Long SM, Singh N, Coggan TL, Askeland MPJ, Jones OAH. Exposure to Environmentally Relevant Levels of PFAS Causes Metabolic Changes in the Freshwater Amphipod Austrochiltonia subtenuis. Metabolites. 2022; 12(11):1135. https://doi.org/10.3390/metabo12111135

Chicago/Turabian Style

Sinclair, Georgia M., Sara M. Long, Navneet Singh, Timothy L. Coggan, Matthew P. J. Askeland, and Oliver A. H. Jones. 2022. "Exposure to Environmentally Relevant Levels of PFAS Causes Metabolic Changes in the Freshwater Amphipod Austrochiltonia subtenuis" Metabolites 12, no. 11: 1135. https://doi.org/10.3390/metabo12111135

APA Style

Sinclair, G. M., Long, S. M., Singh, N., Coggan, T. L., Askeland, M. P. J., & Jones, O. A. H. (2022). Exposure to Environmentally Relevant Levels of PFAS Causes Metabolic Changes in the Freshwater Amphipod Austrochiltonia subtenuis. Metabolites, 12(11), 1135. https://doi.org/10.3390/metabo12111135

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop