
==== Front
Environ Sci Technol
Environ Sci Technol
es
esthag
Environmental Science & Technology
0013-936X
1520-5851
American Chemical Society

39236221
10.1021/acs.est.4c05235
Article
Longitudinal Exposomics in a Multiomic Wellness Cohort Reveals Distinctive and Dynamic Environmental Chemical Mixtures in Blood
https://orcid.org/0000-0001-9463-655X
Sdougkou Kalliroi †
https://orcid.org/0000-0003-2538-8702
Papazian Stefano †‡
https://orcid.org/0000-0001-5141-7111
Bonnefille Bénilde †‡
https://orcid.org/0000-0002-2422-0492
Xie Hongyu †
Edfors Fredrik §
Fagerberg Linn §
Uhlén Mathias §
Bergström Göran ∥⊥
Martin Leah J. #
https://orcid.org/0000-0001-6265-4294
Martin Jonathan W. *†‡
† Department of Environmental Science, Stockholm University, Stockholm 106 91, Sweden
‡ National Facility for Exposomics, Metabolomics Platform, Science for Life Laboratory, Stockholm University, Solna 171 65, Sweden
§ Department of Protein Science, Science for Life Laboratory, KTH-Royal Institute of Technology, Stockholm 100 44, Sweden
∥ Department of Molecular and Clinical Medicine, Institute of Medicine, Sahlgrenska Academy, University of Gothenburg, Gothenburg 40530, Sweden
⊥ Department of Clinical Physiology, Sahlgrenska University Hospital, Region Västra Götaland, Gothenburg 413 45, Sweden
# Independent Researcher, Stockholm 112 33, Sweden
* Email: jon.martin@aces.su.se.
05 09 2024
17 09 2024
58 37 1630216315
26 05 2024
23 08 2024
22 08 2024
© 2024 The Authors. Published by American Chemical Society
2024
The Authors
https://creativecommons.org/licenses/by/4.0/ Permits the broadest form of re-use including for commercial purposes, provided that author attribution and integrity are maintained (https://creativecommons.org/licenses/by/4.0/).

Chemical exposomes can now be comprehensively measured in human blood, but knowledge of their variability and longitudinal stability is required for robust application in cohort studies. Here, we applied high-resolution chemical exposomics to plasma of 46 adults, each sampled 6 times over 2 years in a multiomic cohort, resulting in 276 individual exposomes. In addition to quantitative analysis of 83 priority target analytes, we discovered and semiquantified substances that have rarely or never been reported in humans, including personal care products, pesticide transformation products, and polymer additives. Hierarchical cluster analysis for 519 confidently annotated substances revealed unique and distinctive coexposures, including clustered pesticides, poly(ethylene glycols), chlorinated phenols, or natural substances from tea and coffee; interactive heatmaps were publicly deposited to support open exploration of the complex (meta)data. Intraclass correlation coefficients (ICC) for all annotated substances demonstrated the relatively low stability of the exposome compared to that of proteome, microbiome, and endogenous small molecules. Implications are that the chemical exposome must be measured more frequently than other omics in longitudinal studies and four longitudinal exposure types are defined that can be considered in study design. In this small cohort, mixed-effect models nevertheless revealed significant associations between testosterone and perfluoroalkyl substances, demonstrating great potential for longitudinal exposomics in precision health research.

The first cohort-level application of longitudinal exposomics revealed novel and dynamic coexposures in blood of relevance to precision health.

chemical exposome
longitudinal exposomics
high-resolution mass spectrometry
multiclass targeted
untargeted analysis
blood plasma
Svenska ForskningsrÃ¥det Formas 10.13039/501100001862 2018-02268 VetenskapsrÃ¥det 10.13039/501100004359 2018-03409 document-id-old-9es4c05235
document-id-new-14es4c05235
ccc-price
==== Body
pmcIntroduction

The chemical exposome represents the cumulative sum of environmental chemical exposures throughout an individual’s life course. It includes exposure to natural and anthropogenic chemicals from external sources, such as inhalation of polluted air, intake of dietary substances and pharmaceuticals, and ingestion of contaminated food and water, but also includes internal exposure sources, such as the metabolic products of gut microbiota.1,2 With recognition that our environment is dynamic over time and that our susceptibility to disease changes over the life course, the exposome has always been imagined as a longitudinal endeavor that will require multiple measures of exposure.1,2 For example, the external airborne chemical exposome has been shown to be highly dynamic for individuals over time, and variable among people, even for those living in the same geographical area.3

Comparatively little is known about the dynamics and variability of chemical exposomes in human blood. Sensitive methods for chemical exposomics have been described for blood plasma, including by multiclass targeted analysis4,5 and combined targeted/untargeted analysis by liquid chromatography-high-resolution mass spectrometry (LC-HRMS).6,7 Measurement of the chemical exposome in blood is strategically important because the same sample can be further used for clinical testing, or for multiomic profiling of endogenous molecules whose levels may be impacted by the exposome.8,9 Previous longitudinal studies of small molecules in blood have focused on the metabolome, with limited exploration of environmental chemicals in a few individuals for up to a few months.9,10 Studies of the blood chemical exposome have yet to be reported in a longitudinal cohort, and consequently, there are few measures of inter- and intraindividual variability for priority environmental substances. Such data are necessary to differentiate short-term and chronic exposures and to inform the number of blood sampling events to include in statistically powered exposome studies in the future.

Here, we applied chemical exposomics to recurrent plasma samples of 46 healthy Swedish adults, each sampled 6 times over 2 years in a multiomic wellness profiling study. Using LC-HRMS for quantitative analysis of 83 multiclass targeted analytes (78 priority environmental contaminants and 5 steroid hormones) and parallel untargeted discovery of environmental chemicals and endogenous metabolites,7 here we report the inter- and intraindividual variability for hundreds of molecular environmental exposures, including chemicals not previously detected in human blood. Intraclass correlation coefficients (ICCs) revealed that the stability of the chemical exposome was generally low, compared to parallel measures of plasma proteome, metabolome, lipidome, and microbiome in the same participants. Repeated sampling of the same individuals over time permitted new exposure types to be defined, and for rare and common coexposures to be revealed of relevance to precision health. Mixed-effect modeling also revealed statistically significant exposome-metabolome interactions indicative of endocrine disruption.

Methods

Recruitment and Sampling

Plasma samples were from healthy participants in the Swedish SciLifeLab SCAPIS Wellness Profiling (S3WP) study11 who were previously recruited to the Swedish CArdioPulmonary bioImage Study (SCAPIS), a cohort of 30,000 participants aged 50–65 years and representative of the Swedish population.11,12 S3WP included 6 examination visits in two rounds (4 visits scheduled every 3 months in the first round, and 2 visits scheduled every 6 months in the second round), with 94 of 101 enrolled subjects completing all 6 visits between late 2015 and early 2018.11 At each visit and after overnight fasting, samples of blood, urine, and feces were collected.11 In this work, longitudinal sample sets of plasma from 46 participants (6 samples per person, 276 total samples, 50–200 μL) were selected so that 23 males and 23 females were included, with balanced birth years in the ranges 1950–1955 (n = 14; 7 males and 7 females), 1956–1960 (n = 16; 7 males and 9 females), and 1961–1965 (n = 16; 9 males and 7 females). The study was approved by the Ethical Review Board of Göteborg, Sweden, and all participants provided written informed consent.

Sample Preparation and LC-HRMS Analysis

Plasma samples were prepared and analyzed following a previously described combined targeted and untargeted chemical exposomics method.7 The reader is directed to Sdougkou et al.7 for details on sample preparation, analysis, validation, and quantification; further information can be also found in the Supporting Information and Table S1. The chemical exposomics method was validated for 83 targeted analytes, including environmental contaminants, dietary chemicals, tobacco markers, drugs, and endogenous steroid hormones (Table S1).7 Sample preparation followed a phospholipid removal protocol, and measurements were conducted by ultrahigh pressure LC (Ultimate 3000, Thermo Scientific) with HRMS acquisition (Q Exactive Orbitrap HF-X, Thermo Scientific) in positive and negative electrospray ionization mode (ESI+ and ESI−).7 Spectral acquisition was performed with alternating full scan and data-independent MS/MS acquisition (DIA). Data-dependent acquisition (DDA) with an inclusion list of precursor ions was used for analyte confirmations. Quantification of target analytes was by solvent-based calibration curves with internal standards and by reference standardization,13 which was facilitated by injecting pooled Swedish plasma (see the Supporting Information), for retrospective semiquantification of discovered untargeted substances. For data summaries and statistics, when analytes were detected at concentrations lower than the respective MLOQ, concentrations were substituted by MLOQ/2, and when analytes were nondetect, concentrations were substituted by MLOQ/4.

Untargeted Data Processing and Structural Annotations

For untargeted analysis and spectral library matching, raw data were processed in MS-DIAL (v.4.90)14 for feature alignment across samples, MS1 and DIA MS2 spectral deconvolution, and peak integration (Table S2). Each molecular feature was defined by a chromatographic retention time (RT), an MS1 m/z, and a deconvoluted MS2 spectrum. For feature annotation, spectral matches were considered for total identification scores >700 using MassBankEU,15 MassBank of North America,16 and Global Natural Product Social Molecular Networking (GNPS).17 For each annotated feature, a class (endogenous/environmental) and subclass (e.g., bile acids, PFAS) was assigned following searches of PubChem,18 Human Metabolome Database,19 and Food Metabolome Database.20

The MS-DIAL feature lists from ESI+ and ESI– were combined and analyzed in Python (v. 3.8.16)32 using Jupyter Notebook (v. 6.5.2).33 Redundant features across modes were identified based on mass tolerance of 0.002 Da (after adjusting for (de)protonation) and RT tolerance of 0.2 min, and the redundant feature with the lowest average peak area was discarded. Peak areas of the final combined data set were corrected for instrumental and batch variation using the principal component analysis (PCA) scores of isotope-labeled internal standard signal intensities in each sample,21,22 and then normalized by the analyzed sample volume (complete data reduction workflow in the Supporting Information).

Statistical Analysis and Visualization

PCA was performed in SIMCA (v. 17.0, Umetrics) for the targeted data set. ComplexHeatmap (2.14.0)23 was used to visualize the combined targeted and annotated untargeted data set in R (v. 4.2.1) and R Studio (v. 2023.03.1). Hierachical cluster analysis (HCA) dendrograms and heatmaps were generated based on Pearson correlation as clustering distance and average linkage as the clustering similarity method. For other data visualizations, the Python libraries Plotly (v.5.13.0)24 and Seaborn (v.0.12.2)25 were used. The package SciPy (v. 1.7.3)26 was used for Shapiro-Wilk test, Wilcoxon signed-rank test, Student’s t test, one-way ANOVA, and Pearson correlation coefficient calculation. Statsmodels (v. 0.13.5)27 was used for Bonferroni corrections and Tukey’s test, Pingouin (v. 0.5.3)28 was used for the ICC calculations, and umap-learn (v. 0.5.3) was used for uniform manifold approximation and projection (UMAP).29 For UMAP, default parameters were used: number of neighbors (15), minimum distance (0.1), and Euclidean distance. Unit-variance scaling was applied before the UMAP or PCA analyses. Mixed-effect linear regression analyses were conducted with lme4 (v 1.1–33) in R (v 4.3.0) to examine associations between PFAS and testosterone; models included participant-specific random effects with random intercepts and slopes. Associations between log10 transformed testosterone and log10 transformed PFAS concentrations were examined in unadjusted models, as well as in models adjusted for baseline age and body mass index (BMI) as fixed effects, and in both cases the p-values were Bonferroni corrected for multiple hypotheses using ‘p.adjust’ in R.

Results and Discussion

Multiclass Targeted Analysis

Among all plasma samples, 57 of the 83 target analytes were detected and quantified. These substances belonged to 14 diverse chemical classes (Table S3 and Figure 1a,b), including pesticides (organophosphate and neonicotinoid), flame retardants, PFAS, personal care products, pharmaceuticals, plasticizers, dietary substances, polycyclic aromatic compounds, a nicotine metabolite, and endogenous steroid hormones. For statistical and multivariate analysis, only target analytes with DF > 10% were considered (34 analytes from 9 chemical classes, Figure S1). Comparing mean plasma concentrations of males (n = 138) and females (n = 138), significant differences by sex were identified for parabens as well as certain PFAS and hormones (Figure 1a,b). These were identified with Wilcoxon signed-rank tests (Bonferroni corrected p-values), which were applied after showing non-normal distributions for all 34 analytes (Shapiro-Wilk test, p < 0.05). Methylparaben and propylparaben, which are used in cosmetics, had significantly higher levels in females (methylparaben: males 0.60 ng/mL, females 1.5 ng/mL, p < 0.001; propylparaben: males 0.06 ng/mL, females 0.31 ng/mL, p < 0.01, Figure 1a), consistent with a previous study in urine.30 Conversely, females had significantly lower concentrations of several PFAS (Figure 1b), consistent with known routes of PFAS elimination during pregnancy, breastfeeding, and menstruation,31−34 including for linear PFHxS (males 1.4 ng/mL, females 0.90 ng/mL, p < 0.001), linear PFHpS (males 0.25 ng/mL, females 0.16 ng/mL, p < 0.001), branched PFHpS (males 0.033 ng/mL, females 0.017 ng/mL, p < 0.001) and branched PFOS (males 3.11 ng/mL, females 2.40 ng/mL, p < 0.001). As expected, mean testosterone was significantly higher in males (3.2 ng/mL) than females (0.18 ng/mL, p < 0.001) (Figure 1a), with levels in reference ranges for healthy adults in comparable age groups.35

Figure 1 Targeted chemical exposome and individual variation in S3WP participants. (a, b) Violin plots of average analyte concentrations per individual (log10 scale) for targeted multiclass analytes and targeted perfluoroalkyl and polyfluoroalkyl substances (PFAS), respectively, with color coding according to males (green, n = 23) and females (purple, n = 23). Dashed lines on violin plots indicate mean concentration for each analyte in each sex. Asterisks on top of violin plots indicate the analytes with a significantly different concentration between males and females (Bonferroni corrected p < 0.01; Wilcoxon signed-rank test, n = 138 for each sex). (c, d) Six sampling visits of each individual shown in principal component analysis (PCA) scores plots (3 components; R2Xcum = 43.4%, Q2cum = 29.8%) with samples color-coded by participant ID (c) and sex (d). (e) PCA loadings of the targeted analytes. In the PCA plots, only analytes with a detection frequency >10% are shown, and excluding steroid hormones.

After excluding the detected sex hormones (hydrocortisone, corticosterone, progesterone, and testosterone), PCA was applied to broadly visualize and detect sources of intra- and interindividual variability of environmental exposures (30 targets, Figure 1c–e). In the PCA scores plot colored by S3WP participant ID (Figure 1c), the 6 samples from each individual were often grouped together, indicating relative stability of the target exposome over 2 years. Notably, two individuals were consistently outliers relative to the study population, showing particularly high levels of certain PFAS (i.e., C10–C13 perfluoroalkyl carboxylates for female W0043 colored in yellow; C5–C8 perfluoroalkyl sulfonates for male W0074, colored in gray, Figure 1e). Partial separation between males and females along PC2 (10% of variation, Figure 1d) was mostly explained by sex-specific exposures to parabens and PFAS (Figure 1e), consistent with the statistical differences noted above.

Untargeted Molecular Discovery

Among all 276 plasma samples, a total of 129,547 unique untargeted features were detected across ESI+ and ESI– after blank filtering and removal of redundant features; thus, targeted analytes represented only 0.04% of the overall molecular data set. Substantially fewer features were detectable in any individual (mean 73,426 features per individual, range 67,405–87,159 features, Figure S2). Moreover, only 20,520 features (16% of total data set) were detected in all individuals in at least one visit, somewhat lower than for the targeted analytes (i.e., 26% of targets were detected in all individuals for at least 1 time point). These results indicate that a high proportion of “molecular dark matter”36 in plasma is unique to individuals, possibly representing unique environmental exposures, and/or unique biological response at the metabolome level. This finding is consistent with a recent observation that gut microbiota composition (i.e., internal exposome) explains the majority of variance (i.e., 58%) in individual human plasma metabolites,37 although other sources of environmental exposure have yet to be similarly examined. In untargeted studies of the exposome or the metabolome, it is therefore predictable that the complexity of molecular data sets will increase with larger sample sizes.

By automated matching to open access libraries (based on MS1 accurate mass and MS2 spectra) and after manual curation of all annotations, a total of 462 high-confidence structural candidates were assigned (i.e., at least Level 2a confidence38) in ESI+ and ESI–, corresponding to a 0.4% annotation rate; this total does not include the targeted analytes, which are considered confirmed at Level 1. Combining all Level 1 identified chemicals (targeted analytes and untargeted discoveries confirmed by reference standard;38 see next section), and all Level 2 high-confidence annotations, resulted in a total of 519 substances, including 343 environmental chemicals, 162 endogenous metabolites, and 14 substances with ambiguous classification (Table S4). The 343 environmental chemicals were divided into 11 subclasses, including dietary substances, drugs, industrial chemicals, PFAS, plasticizers, and personal care products (Figure S3), and are hereafter referred to as the “chemical exposome”. The 162 endogenous metabolites were divided into 7 subclasses, including fatty acids, bile acids, amino acids, and hormones. We acknowledge that a strict separation between chemical exposome and endogenous substances leaves certain ambiguities, including for the reason that the human metabolome includes the endogenous metabolomes of other species consumed through diet.39

Confirmation of Selected Structural Annotations

To confirm the Level 2 untargeted molecular annotations, authentic standards (>98% purity) were obtained for 25 analytes, including 3 endogenous and 22 environmental substances. Confirmations were successful for 20 analytes (80% success rate, Tables S5, S6 and Figures S4, S21), indicating an effective untargeted acquisition and data-processing workflow for molecular discovery. When a confirmed analyte was also detectable in the pooled Swedish reference plasma, and its concentration could be quantified with a standard addition curve (0–10 ng/mL range), reference standardization13 was used for semiquantification in individual samples (Table S5). Selected examples of confirmed molecular discoveries, including unexpected and widespread environmental exposures, are discussed below.

Rubber Additives

The substance 1,3-diphenyl guanidine (Table S5 and Figure S4) was recently reported to be the most frequently detected tire-derived contaminant in air40 and indoor dust41 globally and was confirmed here in 83% of samples. Reports of 1,3-diphenyl guanidine in human biofluids remain rare.42−44 A related substance, 4-tert-butylcatechol (TBC, Table S5 and Figure S7), was detected in 92% of samples, despite no previous reports in human biofluids. This substance is reported as a contact allergen,45 and is used as an additive polymerization inhibitor in the rubber, paint, and petroleum industry.45,46

Industrial Chemicals

Triphenylphosphine oxide (TPPO) is a widely used synthetic intermediate in pharmaceutical products47 and was detected here in 3% of samples (Table S5 and Figure S5). TPPO has begun to be widely detected in indoor air, dust, and aquatic systems,48,49 but has only rarely been detected in human samples.42,47 The isomers 1- and 2-naphthalenesulfonate (Table S5 and Figure S12), used in textile, pharmaceutical, and agrochemical production,50,51 were confirmed in 13% of samples here. These two substances have low biodegradability50,51 and neither has been previously confirmed in human biofluids (2-naphthalenesulfonate has been reported previously at Level 27).

We also confirmed 2,6-ditert-butyl-4-nitrophenol (DBNP, Table S5 and Figure S6), a known transformation product of 2,6-ditert-butyl-4-methylphenol, a widely used synthetic antioxidant added to polymers, foods, and cosmetics.52 DBNP has been reported in the environment53,54 and in plastic food packaging55 but not in human biofluids. Its detection in 20% of Swedish samples here deserves further attention due to its biological persistence and toxic potential.54 Two benzotriazole isomers, 4-methyl- and 5-methyl-1H-benzotriazole, were also confirmed and semiquantified in 9% of samples (Table S5 and Figure S10). These benzotriazoles are high-production volume chemicals mainly applied as corrosion inhibitors and ultraviolet stabilizers, widely detected in environmental matrices56−58 but also in a few studies in human samples.59,60

Pesticides

Chlorothalonil-4-hydroxy was confirmed and semiquantified in 100% of plasma samples (1.1–11.5 ng/mL, median 3.7 ng/mL; Table S5 and Figure S9). This finding is consistent with a targeted study of pregnant Swedish women (1997–2015) where it was also detected in all samples (median 4.1 ng/mL),61 yet it remains rarely monitored in humans. It is a known transformation product of the fungicide chlorothalonil, which has not been permitted for agricultural use in Sweden since the 1990s,62 and has been banned in the EU since 2019 due to its carcinogenic properties and risk to fish and amphibians.63 Chlorothalonil-4-hydroxy is considered more toxic,64 more persistent, and more mobile in soil65 than the parent pesticide; thus, its widespread presence in blood deserves attention in exposome studies.

Personal Care Products

Sodium lauryl sulfate was confirmed in 45% of samples (Table S5 and Figure S8), and has not previously been reported in human biofluids to our knowledge, despite being a major surfactant in shampoos.66,67 It is known to be absorbed into the bloodstream in animal models,68 and may also be inhalable during shampoo use.69 This surfactant is rarely monitored and remains unregulated due to its biodegradable nature, but its risks have been debated.70

Longitudinal Stability of the Chemical Exposome

The final data set of 519 annotated substances (Levels 1 and 2) was used to examine the intra- and interindividual exposure variation in this 2 year study with six sampling points. For each substance, we calculated the ICC (Table S7), a nondimensional ratio of the interindividual variance to the total variance (i.e., the sum of inter- and intraindividual variance).71 In the context of exposome research, ICCs are instructive for study design as they describe the extent to which individuals retain their rank order in a study population with repeated measurements of exposure over time.72 Higher ICCs correspond to more stable exposures, that can therefore be measured fewer times throughout the life course, whereas lower ICCs correspond to less stable exposures that may need several repeated measurements to accurately classify exposure over the life course. ICCs range from 0 to 1, with values <0.40 corresponding to poor reproducibility of repeated measurements, 0.40 to 0.75 is considered fair to good, and >0.75 indicates excellent reproducibility.71

The majority of annotated substances in plasma (306 of 519 analytes) had ICCs < 0.40 (Figure 2a), and the mean ICC was significantly higher for endogenous metabolites (0.40) than for the chemical exposome (0.30, Student’s t test, two-tailed, p < 0.001) (Figure 2b). Moreover, while the ICCs for endogenous metabolites were normally distributed, the chemical exposome showed denser distributions toward lower and higher ICC values. More specifically, 66% of ICCs were <0.4 for the chemical exposome, compared to only 45% of ICCs for endogenous metabolites, and 10% of the ICCs were >0.75 for the chemical exposome, compared to only 6% of ICCs for endogenous metabolites (Figure 2b). While the plasma chemical exposome and metabolome have highly stable components, these results mean that the majority of small molecules in plasma must be measured more than once over time to adequately represent an individual’s exposure. A similar conclusion was reached for the target analysis of 24 nonpersistent environmental chemicals in urine of pregnant women and children.73

Figure 2 Longitudinal stability of the chemical exposome and comparison to other molecular profiles. (a) Calculated intraclass correlation coefficients (ICC) of the 519 annotated substances (Level 1 and Level 2) along with their detection frequency (DF), for the chemical exposome (orange, n = 343), endogenous metabolites (blue, n = 162) and compounds with ambiguous origin (gray, n = 14). Four regions of interest are highlighted on the plot: top left (Type 1; rare-stable, low DF, high ICC), top right (Type 2; common-stable, high DF, high ICC), bottom right (Type 3; common-unstable, high DF, low ICC), bottom left (Type 4; rare-unstable, low DF, low ICC). (b) ICC data presented in violin plots separately for the chemical exposome and the endogenous metabolites as well as for the proteome, lipidome, metabolome, and microbiome from Tebani et al.,11 including only individuals present in the exposome study. The solid and dashed lines show median and mean ICC values, respectively, and asterisks indicate significantly different ICC means between the chemical exposome and endogenous metabolites (p < 0.001; Student’s t test, two-tailed) and between the chemical exposome and all other molecular profiles (p < 0.001; one-way ANOVA and Tukey’s test). (c, d) Global profiles across visits applying uniform manifold approximation and projection (UMAP) for the chemical exposome and the endogenous metabolites color-coded by participant ID with gray lines connecting samples representing different visits of the same individual (c) or color-coded by sex (d).

To conclusively link environmental exposures to disease etiology, an objective in exposomics is to identify key perturbations in molecular networks at the metabolome, proteome, and gene expression levels.74 In longitudinal exposome studies, this might be addressed through multiomics, but the relative dynamics of each molecular data set must be considered at the study design stages. Thus, here we compare our results to those of Tebani et al., who reported ICCs for several multiomic profiles in the same study participants.11 We recalculated the ICCs of Tebani et al. to only include individuals with paired chemical exposome measurements (see the Supporting Information), and compared to the relatively low stability of the chemical exposome (mean ICC 0.30, median 0.23), significantly higher mean ICCs were evident for the plasma proteome (mean 0.65, median 0.67), lipidome (mean 0.55, median 0.59), metabolome (mean 0.50, median 0.50), and gut microbiota (mean 0.38, median 0.35) (Figure 2b; one-way ANOVA, Tukey’s test p < 0.001). These results emphasize the importance of repeated measures of the chemical exposome in epidemiological studies to minimize exposure misclassification, and in longitudinal multiomic studies, the exposome should be measured as frequently, or more frequently, than other biomolecular profiles.

The gut microbiome had the most similar distribution of ICCs in comparison to the chemical exposome, with overall lower reproducibility and similar high density at lowest values (Figure 2b). Like the chemical exposome, gut microbiota is shaped by life course exposures, including diet, disease history, and medication, as well as by intrinsic factors such as age and host genetic variation.37 As discussed above, the known link between gut microbiota and plasma small molecules37 could partially explain similar ICC distributions for gut microbiome and the plasma chemical exposome.

UMAP analysis allowed the relative variation of individual exposomes over time to be visualized relative to the study population (Figure 2c). The visits of each individual, colored by participant ID and connected by lines, demonstrate that the complex chemical exposome profile of most individuals has stable factors over time, despite the overall low temporal stability described above (Figure 2c). However, individual variability was evident for both the chemical exposome and endogenous metabolites, with some individuals having remarkably unique and stable profiles and others displaying higher variability between visits (Figure 2c). UMAP visualization of metabolome stability in Tebani et al.11 showed similar results. Notably, UMAP projections of the chemical exposome and endogenous metabolites did not indicate any sex differences (Figure 2d).

Longitudinal Exposure Types

A plot of each chemical’s ICC versus sample detection frequency (Figure 2a) allowed for new insights into chemical exposure characterization that would not have been evident without repeated sampling. Four exposure types (insets#1–4, Figure 2a; see Table S7 for the categorization of each compound) are defined and discussed in this context. First, we categorize Type 1 exposures as rare-stable (i.e., approximately DF < 20%, ICC 0.6–1, 2% of the data set). These exposures may be of high relevance to the small fraction of individuals with consistently elevated exposure over the life course, but could be easily overlooked in cross-sectional study designs, due to being rare in the population. These exposures are therefore of high relevance in a precision health context, as they likely arise from unique behaviors, occupation, or lifestyle factors that could be mitigated if identified. An example was the industrial intermediate triphenylphosphine oxide (ICC 0.64, 3.2% DF; Level 1) which was only detected in 9 samples overall, but consistently in all 6 samples of individual W0015. Other examples included NEtFOSAA (ICC 0.93, 7% DF; Level 1) (consistent for individual W0008), a PFOS precursor with exposure sources linked to food packaging75 and the drug citalopram (ICC 0.88, 9% DF; Level 2)75 (consistent for individual W0090), which acts as a selective serotonin reuptake inhibitor.76

Type 2 exposures were categorized as common-stable (i.e., approximately DF > 80%, ICC 0.6–1, 12% of the data set) and include exposures that are widespread in the population, and with stable levels in individuals. A large proportion of these substances were PFAS and chlorinated analytes, which have biological persistence and long pharmacokinetic half-lives. For example, PFHxS had an ICC of 0.94 (100% DF, reported half-life 8.5 yrs)77 and PFOS isomers had ICCs of 0.95 (100% DF, reported half-life 5.4 yrs for linear PFOS).77 Another example was the fungicide transformation product chlorothalonil-4-hydroxy (ICC 0.80, 100% DF; confirmed Level 1). The half-life of chlorothalonil-4-hydroxy in humans has not yet been determined, but the high ICC reported here and in a previous study of Costa Rican women (ICC 0.81)61 emphasizes the importance of further studies for this environmentally persistent chemical. Other Type 2 substances nevertheless included analytes with reported fast elimination half-lives in the range of days or hours. Such examples were 3-hydroxycotinine (ICC 0.91, 88% DF, half-life 6.6 h;78 confirmed Level 1), and the pesticide pentachlorophenol (ICC 0.83, 98% DF, half-life 20 d;79 Level 1). These latter results suggest that, despite relatively rapid elimination pharmacokinetics, consistent lifestyle factors over the course of 2 years (e.g., tobacco use) can result in steady levels of environmental substances in plasma. In fact, the above ICCs were comparable to that of the targeted steroid hormone testosterone (ICC 0.95, 100% DF; Level 1). While substances in the Type 2 category do not necessarily require repeated sampling to accurately classify exposure in health studies, longitudinal studies with untargeted chemical exposomics may be a powerful approach to discovering emerging persistent chemicals in the population. For example, among the thousands of nonannotated molecular features detected here, these could be prioritized for identification based on simultaneously high ICC and DF.

Type 3 exposures were categorized as common-unstable (i.e., approximately DF > 80%, ICC < 0.4, 18% of the data set). Substances in this category are widespread in the population, but have unstable levels in individuals; thus, the health relevance of these can most powerfully be explored in longitudinal studies with repeated measurements. Examples here included substances associated with exposure through diet, and/or having short elimination half-lives, such as the insecticide metabolite 3,5,6-trichloro-pyridinol (ICC 0.25, 97% DF, half-life 27 h;80 Level 1), which is linked to consumption of imported foods in Sweden.81 Other examples are the caffeine metabolite theobromine (ICC 0.31, 100% DF, half-life 7–12 h;82 Level 2), azelaic acid which is found in grains83 (ICC 0.40, 100% DF; confirmed Level 1), the herbicide ioxynil84 (ICC 0.37, 92% DF; Level 2) and the above-discussed rubber additive 1,3-diphenyl guanidine (ICC < 0.01, 83% DF, half-life 9.6 h in rats;41 Level 1).

Finally, Type 4 exposures were categorized as rare-unstable (i.e., approximately DF < 20%, ICC < 0.4, 23% of the data set). These exposures are more difficult to understand with regard to health and could be deprioritized in exposome studies to further reduce the dimensional complexity and multiple testing. Examples of Type 4 exposures are the industrial chemicals 2,6-ditert-butyl-4-nitrophenol (ICC 0.33, 20% DF; Level 1) and the isomers 1- and 2-naphthalenesulfonate (ICC 0.29, 13% DF; Level 1), both discussed above.

Correlated Coexposures

The complexity and high dimensionality of the chemical exposome represent major challenges to its implementation in health studies; however, it has been proposed that reduction of its dimensional complexity should be possible by grouping correlated exposures.74 Such coexposures may arise if two or more substances have common exposure sources and similar pharmacokinetics, resulting in correlated dynamics. Here we highlight several coexposures among the 46 participants that were evident in HCA heatmaps (Figures 3a,and S22), and for discussion we group these into “common” and “rare” coexposures based on whether these were evident among many participants (Figure 3a, groups A–F, zoom-in Figure S23), or only in individuals or a small fraction of the population (Figure 3a, groups G–M, zoom-in Figure S24), respectively. Both HCA heatmaps are based on all 46 participants and 519 annotated substances (Levels 1 and 2), but for simplicity, we focus on Figure 3a, which is based on the average response of each substance for each individual’s six visits. Figure S22 shows the chemical exposome at each visit to simultaneously visualize intra- and interindividual patterns. To guide interpretation and hypothesis generation, both heatmaps were made interactive (available at https://s3wp-exposomics.serve.scilifelab.se) and metadata are displayed for each substance (i.e., chemical class, subclass, confidence level, ionization mode, DF) and participant (sex, age, BMI).

Figure 3 Clustered exposome profiles of S3WP participants and examples of correlated exposures. (a) Hierarchical cluster analysis heatmap with dendrograms showing the exposome profiles of 46 individuals, each averaged across the 6 clinical visits for 519 annotated substances (Level 1 and Level 2). Color coding of chemical substances (“molecular features”) is according to class, subclass, confidence level of identification, ionization, and detection frequency (DF). Color coding of individuals is according to sex, age, and BMI. Groups of interest are highlighted on the heatmap (in red A–F: common coexposures; in blue G–M: rare coexposures; in black N–P: endogenous metabolites, see zoomed-in versions in Figures S23, S24, and S27). Common coexposures were detected in many participants while rare coexposures were detected only in individuals or a small population fraction. (b) Linear regression plots with 95% confidence intervals across all 6 visits (n = 276) between features that clustered together in the heatmap. The Pearson correlation coefficient (r) is shown for each regression and all analyte pairs showed a significant correlation (p-value <0.001).

One common coexposure (Group A, Figure 3a, zoom-in Figure S23) included the targeted fungicide pentachlorophenol (ICC 0.83, 98% DF) and the dietary fatty acid, pentadecanoic acid, which is a biomarker of dairy consumption85 (Level 2; ICC 0.29, 93% DF; r 0.43, p < 0.005). Pentachlorophenol is also a metabolite of the organochlorine pesticide hexachlorobenzene86 which has been associated with dairy consumption in Sweden.87 Another common coexposure (Group B, Figure S23) included homologues of poly(ethylene glycol) (Level 2, ICC 0–0.31, 36–100% DF), which are approved by the US Food and Drug Administration88 and have been previously detected in environmental,89,90 but not in human samples. Group C coexposures (Figure S23) included sodium lauryl sulfate (ICC 0.09, 45% DF, confirmed Level 1), strongly correlated with tetradecyl sulfate (Level 2; ICC 0.07, 36% DF; r 0.92, p < 0.001, n = 276, Figure 3b) and lauryl diethanolamide (Level 2; ICC 0.01, 32% DF; r 0.66, p < 0.001, n = 276), all of which are ingredients in household cleaning and personal care products.66,67,91 Such correlations between confirmed analytes (targeted/Level 1) and one or more Level 2 substances in a related chemical class provide further confidence in the Level 2 annotations, and more examples are discussed below.

In Group D (Figure S23), the targeted analyte caffeine (ICC 0.50, 99% DF) correlated strongly with the Level 2 caffeine metabolites92 theobromine (ICC 0.31, 99.6% DF; r 0.60, p < 0.001, n = 276 individual samples) and theophylline (ICC 0.52, 100% DF; r 0.80, p < 0.001, n = 276, Figure 3b). While these three alkaloids are typical of coffee plants, they also clustered with 4-methylcatechol (Level 2; ICC 0.60, 99% DF), a compound generated during the roasting of coffee beans.93 Similarly, in Group E (Figure S23) the targeted nicotine metabolite, cotinine (ICC 0.94, 20% DF), was strongly correlated with a secondary nicotine metabolite, 3-hydroxycotinine94 (ICC 0.91, 88% DF; r 0.97, p < 0.001, n = 276, Figure 3b), confirmed here at Level 1. Group F (Figure S23) shows the prominent PFAS class as a coexposure cluster (i.e., many PFAS were highly correlated, Figure S25). Additionally, the HCA heatmap with all 6 visits (Figure S22) provides visual evidence for the high ICCs of these substances, whereby the majority of individuals showed consistently high (red) or low (blue) responses across all visits (also see zoom-in subset in Figure 4a).

Figure 4 Stable PFAS exposome profiles and examples of exposome-metabolome interactions. (a) Zoomed-in subset of the hierarchical cluster analysis heatmap (19 out of 46 participants) with 6 visits shown per individual where a stable group of PFAS is observed. Color coding of features is according to class, subclass, confidence level of identification, ionization, and detection frequency (DF). Color coding of individuals is according to sex, age, and BMI. (b) Linear regression plots with 95% confidence intervals between testosterone and PFAS targets for female individuals across all 6 visits (n = 138). All analyte pairs showed significant correlation (p-value <0.001). The plots also show the analyte distributions in histograms.

Examples of rare coexposures involved groups of chemicals occurring at low DF in the population (approximately <20%), but also include distinctively higher levels of more commonly detected substances in individuals. In Group G (Figure 3a), one distinctive individual had multiple rare coexposures (zoom-in Figure S24), including triphenylphosphine oxide (ICC 0.64, 3.2% DF; Level 1) which was detected in all 6 samples of this individual. As noted above, this substance is used as an intermediate in pharmaceutical products;47 thus, it was interesting to note correlated coexposures in this individual to 2 drugs, trimethoprim95 (Level 2, ICC < 0.01, 14% DF) and ketoprofen96 (Level 2, ICC 0.34, 57% DF), as well as to the targeted flame-retardant, bis(1,3-dichloro-2-propyl) phosphate (ICC < 0.01, 12% DF) and the long-chain PFAS, perfluorotetradecanoate (PFTeDA; ICC 0.31, 9% DF). The reasons for these shared rare exposures are not currently understood, but deserve follow-up study, and would not have been uncovered without multiclass targeted and untargeted chemical exposomics.

In Group J (Figure S24), two rare exposures, namely, the targeted flame-retardant and plasticizer tris(2-butoxyethyl) phosphate (ICC 0.01, 9% DF) and the plasticizer diethylene glycol dibenzoate (Level 2; ICC < 0.01, 9% DF), were strongly correlated (r 0.82, p < 0.001, n = 276, Figure 3b). In another rare coexposure (Group K, Figure S24) the targeted personal care product triclosan (ICC 0.45, 16% DF) clustered with several other chlorinated substances, such as the herbicide 4-chlorophenoxyacetic acid97 (Level 2; ICC 0.50, 48% DF), the pesticide chloroxynil98 (Level 2; ICC 0.01, 33% DF), and an isomer of chlorophenol (Level 2; ICC 0.16, 2.2% DF), the latter discovered based on exact mass, chlorine isotopic pattern (Figure S11) and RT adjacent to 4-chlorophenol (confirmed Level 1). In Group L, the targeted herbicide diuron (ICC < 0.01, 20% DF) clustered with the targeted fungicide carbendazim (ICC < 0.01, 4% DF) and the herbicide ioxynil84 (Level 2; ICC 0.37, 92% DF) (Figure S24).

Chemical Exposomics as a Tool for Precision Health

An individual’s drug metabolism capacity is influenced by both genetic and environmental factors, and understanding these factors is a key objective of personalized medicine, for example, toward precision dosing of drugs.99 We hypothesized that untargeted chemical exposomics could be a useful tool for monitoring drug efficacy (or toxicity) at the individual level by monitoring drugs and detecting their known or unknown metabolites in plasma. We first observed that groups of rare coexposures often consisted of drugs and their known metabolites (Figure S24). For instance, the anti-inflammatory drug diclofenac (Level 1; ICC 0.19, 6.9% DF) and 5-hydroxy diclofenac100 (Level 2; ICC 0.37, 9.8% DF) in Group H, the antihistamine cetirizine101 (Level 2; ICC 0.88, 21% DF) and hydroxyzine (Level 2; ICC0.83, 68% DF) in Group I, and the antidepressant sertraline (Level 2; ICC 0.32, 2.9% DF) and norsertraline102 (Level 2; ICC 0.33, 2.5% DF) in Group M. To further explore the potential of chemical exposomics to identify drug metabolites whose MS2 spectra are unknown, or not publicly available, we suspect-screened all 38,499 nonannotated untargeted features detected in ESI- for: (i) theoretical m/z corresponding to hydroxyl, sulfate or glucuronide metabolites of the above-mentioned pairs of parent drugs and metabolites, and (ii) correlation with the respective parent drug or metabolites (i.e., p < 0.001, n = 276 individual samples, Figure S26). By these combined criteria, a suspected additional hydroxy metabolite of diclofenac (m/z 310.0049, RT 10.74 min, 2 ppm error, ESI-) was found that correlated both with diclofenac (r = 0.61, m/z 294.0094, RT 13.16 min, ESI-) and 5-hydroxy diclofenac (r = 0.92, m/z 312.0190, RT 11.42 min, ESI+), and the corresponding suspect sulfate conjugate (m/z 389.9616, RT 9.50 min, 1.3 ppm error, ESI-) correlated moderately with diclofenac (r = 0.63), and strongly with the new hydroxy metabolite (r = 0.98). Finally, a suspected glucuronide of sertraline (sertraline carbamoyl-O-glucuronide, m/z 524.0884, RT 13.75 min, 0.1 ppm error, ESI-) correlated strongly with both sertraline (r = 0.97, m/z 306.0812, RT 14.09 min, ESI+) and norsertraline (r = 0.98, m/z 292.0653, RT 14.31 min, ESI+). These results demonstrate that untargeted chemical exposomics is a feasible tool to inform the pharmaceutical treatment of individuals.

Hierarchical clustering of natural compounds and endogenous metabolites (Level 2, Figure 3a, groups N–P, zoom-in Figure S27) was also observed. These may be useful indicators of behavior (e.g., dietary habits) or biomarkers of disorders or disease risk. For example, in group N, a single individual deviated from the rest of the population due to the high relative response of several bile acids, including cholic, glycochenodeoxycholic, and ursocholic acid (Figures S27 and S28a), which can indicate hepatic impairment.103 Moreover, 12 participants clustered with relatively high responses of several lipids, including the fatty acids cis-4,10,13,16-docosatetraenoic acid (ICC 0.55, 100% DF), 9Z, 12Z-linoleic acid (ICC 0.44, 100% DF), and arachidonic acid (ICC 0.42, 100% DF), with high correlations among these (Group O, Figures S27 and S28b). High correlations were observed between the microbial metabolite phenylacetylglutamine (ICC 0.65, 100% DF) and the two uremic toxins p-cresol sulfate (ICC 0.62, 100% DF; r 0.80, p < 0.001, n = 276) and 3-indoxyl sulfate (ICC 0.51, 100% DF; r 0.49, p < 0.001, n = 276) (Group P, Figures S27 and S28c). This microbial metabolite and the associated uremic solutes accumulate in patients with chronic kidney disease, and their levels have been associated with adverse outcomes including cardiovascular disease.37,104

Exposome-Metabolome Interaction

Correlations between environmental chemicals and endogenous metabolites can be suggestive of exposome-metabolome interactions, whereby exposure(s) induces a metabolic response. Although the current study is relatively small, we examined for evidence of endocrine disruption between testosterone and PFAS, which are targeted analytes that were detected in most samples. Testosterone and PFAS did not cluster together in the HCA heatmaps because men have much higher testosterone levels than women, and any correlations may be sex-specific. Among all females, positive correlations were observed between testosterone and several PFAS (Figure 4b), including linear PFOS, perfluorononanoate (PFNA), and perfluorodecanoate (PFDA). These three associations were statistically significant in mixed-effect models, which account for the repeated sampling of individuals, and PFOS and PFDA remained significant in adjusted models controlling for baseline age and BMI (Table S8; p < 0.02 after Bonferroni correction). Positive associations between PFOS, PFOA, and PFHxS with total and free testosterone have been reported in postmenopausal women (median 63 yrs),105 similar to the age of women here (median 57 yrs). However, the same study reported no similar effect for PFNA and PFDA on any reproductive hormones.105 These findings for PFAS deserve further study considering that a meta-analysis of prospective studies has reported that higher blood testosterone was associated with an increased risk of breast cancer in postmenopausal women.106

Advances, Limitations, and Future Directions

This study represents the first application of longitudinal chemical exposome profiling in the human blood of multiple individuals. The combined targeted and untargeted chemical exposomics workflow enabled sensitive and precise quantification of priority targeted analytes with simultaneous potential to discover unexpected exposures of potential health relevance. In this longitudinal study of only 46 individuals, we also report unique coexposures, and statistically significant exposome-metabolome interactions, demonstrating promise for chemical exposomics in precision health applications.107 Nevertheless, the vast majority of untargeted molecular features remain unidentified, and the underlying laboratory methods and data analyses must be scaled up for larger studies. While it has been envisioned that unraveling the effects of the chemical exposome on the progression of disease will require network science and systems biology approaches to integrate chemical and biomolecular interactions across many omic levels,74 the relatively low stability of the chemical exposome stands as a major practical challenge to this approach. The wide range of ICCs for hundreds of endogenous and exogenous small molecules reported here can therefore assist in the design of future longitudinal cohorts focused on the chemical exposome.

Data Availability Statement

Interactive HCA heatmaps are available at https://s3wp-exposomics.serve.scilifelab.se. Access to the MS data and participant-level data sets will be made available upon reasonable request by contacting the corresponding author

Supporting Information Available

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.est.4c05235.Methods; detected target analytes (Figure S1); number of annotated features (Figure S2); total number of targeted and annotated untargeted analytes (Figure S3); extracted ion chromatogram and DDA spectrum (Figures S4–S21); hierarchical cluster analysis heatmap (Figure S22); zoomed-in sections of the hierarchical cluster analysis heatmap (Figures S23 and S24); correlations observed for PFAS targeted analytes (Figure S25); correlations observed for Level 1 and Level 2 drugs and drug metabolites (Figure S26); zoomed-in sections of hierarchical cluster analysis heatmap (Figure S27); correlations observed for Level 2 analytes (Figure S28); correlation observed for BADGE-H2O annotated at Level 2 (Figure S29); untargeted molecular discoveries (Table S5)(PDF)

Targeted analytes, method validation, and corresponding native standards (Table S1); MS-DIAL parameters (Table S2); detected and quantified targeted analytes (Table S3); Levels 1 and 2 annotated substances (Table S4); analytes selected for confirmation (Table S6); intraclass correlation coefficient (Table S7); and mixed-effect model parameters (Table S8) (XLSX)

Supplementary Material

es4c05235_si_001.pdf

es4c05235_si_002.xlsx

The authors declare no competing financial interest.

Acknowledgments

This research was supported through a grant from the Swedish Research Council to J.W.M. (VR, 2018-03409) and from the Swedish Research Council for Sustainable Development to J.W.M., F.E., L.F., and M.U. (Formas, 2018-02268). S.P., B.B., and J.W.M. acknowledge Swedish national funding to the National Facility for Exposomics (SciLifeLab). The authors thank Arnold Kochari from the SciLifeLab Data Centre for the guidance in publicly depositing the interactive HCA heatmaps.
==== Refs
References

Wild C. P. Complementing the Genome with an “Exposome”: The Outstanding Challenge of Environmental Exposure Measurement in Molecular Epidemiology. Cancer Epidemiol., Biomarkers Prev. 2005, 14 (8 ), 1847–1850. 10.1158/1055-9965.EPI-05-0456.16103423
Wild C. P. The Exposome: From Concept to Utility. Int. J. Epidemiol. 2012, 41 (1 ), 24–32. 10.1093/ije/dyr236.22296988
Jiang C. ; Wang X. ; Li X. ; Inlora J. ; Wang T. ; Liu Q. ; Snyder M. Dynamic Human Environmental Exposome Revealed by Longitudinal Personal Monitoring. Cell 2018, 175 (1 ), 277–291.e31. 10.1016/j.cell.2018.08.060.30241608
Jamnik T. ; Flasch M. ; Braun D. ; Fareed Y. ; Wasinger D. ; Seki D. ; Berry D. ; Berger A. ; Wisgrill L. ; Warth B. Next-Generation Biomonitoring of the Early-Life Chemical Exposome in Neonatal and Infant Development. Nat. Commun. 2022, 13 (1 ), 2653 10.1038/s41467-022-30204-y.35550507
González-Domínguez R. ; Jáuregui O. ; Queipo-Ortuño M. I. ; Andrés-Lacueva C. Characterization of the Human Exposome by a Comprehensive and Quantitative Large-Scale Multianalyte Metabolomics Platform. Anal. Chem. 2020, 92 (20 ), 13767–13775. 10.1021/acs.analchem.0c02008.32966057
Flasch M. ; Fitz V. ; Rampler E. ; Ezekiel C. N. ; Koellensperger G. ; Warth B. Integrated Exposomics/Metabolomics for Rapid Exposure and Effect Analyses. JACS Au 2022, 2 (11 ), 2548–2560. 10.1021/jacsau.2c00433.36465551
Sdougkou K. ; Xie H. ; Papazian S. ; Bonnefille B. ; Bergdahl I. A. ; Martin J. W. Phospholipid Removal for Enhanced Chemical Exposomics in Human Plasma. Environ. Sci. Technol. 2023, 57 , 10173–10184. 10.1021/acs.est.3c00663.37394749
Rappaport S. M. ; Barupal D. K. ; Wishart D. ; Vineis P. ; Scalbert A. The Blood Exposome and Its Role in Discovering Causes of Disease. Environ. Health Perspect. 2014, 122 (8 ), 769–774. 10.1289/ehp.1308015.24659601
Gao P. ; Shen X. ; Zhang X. ; Jiang C. ; Zhang S. ; Zhou X. ; Rose S. M. S.-F. ; Snyder M. Precision Environmental Health Monitoring by Longitudinal Exposome and Multi-Omics Profiling. Genome Res. 2022, 32 , 1199–1214. 10.1101/gr.276521.121.35667843
Oosterwegel M. J. ; Ibi D. ; Portengen L. ; Probst-Hensch N. ; Tarallo S. ; Naccarati A. ; Imboden M. ; Jeong A. ; Robinot N. ; Scalbert A. ; Amaral A. F. S. ; van Nunen E. ; Gulliver J. ; Chadeau-Hyam M. ; Vineis P. ; Vermeulen R. ; Keski-Rahkonen P. ; Vlaanderen J. Variability of the Human Serum Metabolome over 3 Months in the EXPOsOMICS Personal Exposure Monitoring Study. Environ. Sci. Technol. 2023, 57 , 12752–12759. 10.1021/acs.est.3c03233.37582220
Tebani A. ; Gummesson A. ; Zhong W. ; Koistinen I. S. ; Lakshmikanth T. ; Olsson L. M. ; Boulund F. ; Neiman M. ; Stenlund H. ; Hellström C. ; Karlsson M. J. ; Arif M. ; Dodig-Crnković T. ; Mardinoglu A. ; Lee S. ; Zhang C. ; Chen Y. ; Olin A. ; Mikes J. ; Danielsson H. ; von Feilitzen K. ; Jansson P.-A. ; Angerås O. ; Huss M. ; Kjellqvist S. ; Odeberg J. ; Edfors F. ; Tremaroli V. ; Forsström B. ; Schwenk J. M. ; Nilsson P. ; Moritz T. ; Bäckhed F. ; Engstrand L. ; Brodin P. ; Bergström G. ; Uhlen M. ; Fagerberg L. Integration of Molecular Profiles in a Longitudinal Wellness Profiling Cohort. Nat. Commun. 2020, 11 (1 ), 4487 10.1038/s41467-020-18148-7.32900998
Bergström G. ; Berglund G. ; Blomberg A. ; Brandberg J. ; Engström G. ; Engvall J. ; Eriksson M. ; Faire U. ; de Flinck A. ; Hansson M. G. ; Hedblad B. ; Hjelmgren O. ; Janson C. ; Jernberg T. ; Johnsson Å. ; Johansson L. ; Lind L. ; Löfdahl C.-G. ; Melander O. ; Östgren C. J. ; Persson A. ; Persson M. ; Sandström A. ; Schmidt C. ; Söderberg S. ; Sundström J. ; Toren K. ; Waldenström A. ; Wedel H. ; Vikgren J. ; Fagerberg B. ; Rosengren A. The Swedish CArdioPulmonary BioImage Study: Objectives and Design. J. Intern. Med. 2015, 278 (6 ), 645–659. 10.1111/joim.12384.26096600
Go Y.-M. ; Walker D. I. ; Liang Y. ; Uppal K. ; Soltow Q. A. ; Tran V. ; Strobel F. ; Quyyumi A. A. ; Ziegler T. R. ; Pennell K. D. ; Miller G. W. ; Jones D. P. Reference Standardization for Mass Spectrometry and High-Resolution Metabolomics Applications to Exposome Research. Toxicol. Sci. 2015, 148 (2 ), 531–543. 10.1093/toxsci/kfv198.26358001
Tsugawa H. ; Cajka T. ; Kind T. ; Ma Y. ; Higgins B. ; Ikeda K. ; Kanazawa M. ; VanderGheynst J. ; Fiehn O. ; Arita M. MS-DIAL: Data-Independent MS/MS Deconvolution for Comprehensive Metabolome Analysis. Nat. Methods 2015, 12 (6 ), 523–526. 10.1038/nmeth.3393.25938372
MassBank | Database | Search. https://massbank.eu/MassBank/Search (accessed May 15, 2024).
MassBank of North America. https://mona.fiehnlab.ucdavis.edu/ (accessed March 13, 2021).
GNPS - Analyze, Connect, and Network with your Mass Spectrometry Data. https://gnps.ucsd.edu/ProteoSAFe/static/gnps-splash.jsp (accessed May 14, 2021).
Kim S. ; Chen J. ; Cheng T. ; Gindulyte A. ; He J. ; He S. ; Li Q. ; Shoemaker B. A. ; Thiessen P. A. ; Yu B. ; Zaslavsky L. ; Zhang J. ; Bolton E. E. PubChem 2023 Update. Nucleic Acids Res. 2023, 51 (D1 ), D1373–D1380. 10.1093/nar/gkac956.36305812
Wishart D. S. ; Guo A. ; Oler E. ; Wang F. ; Anjum A. ; Peters H. ; Dizon R. ; Sayeeda Z. ; Tian S. ; Lee B. L. ; Berjanskii M. ; Mah R. ; Yamamoto M. ; Jovel J. ; Torres-Calzada C. ; Hiebert-Giesbrecht M. ; Lui V. W. ; Varshavi D. ; Varshavi D. ; Allen D. ; Arndt D. ; Khetarpal N. ; Sivakumaran A. ; Harford K. ; Sanford S. ; Yee K. ; Cao X. ; Budinski Z. ; Liigand J. ; Zhang L. ; Zheng J. ; Mandal R. ; Karu N. ; Dambrova M. ; Schiöth H. B. ; Greiner R. ; Gautam V. HMDB 5.0: The Human Metabolome Database for 2022. Nucleic Acids Res. 2022, 50 (D1 ), D622–D631. 10.1093/nar/gkab1062.34986597
FooDB. https://foodb.ca/ (accessed February 15, 2023).
Surowiec I. ; Johansson E. ; Torell F. ; Idborg H. ; Gunnarsson I. ; Svenungsson E. ; Jakobsson P.-J. ; Trygg J. Multivariate Strategy for the Sample Selection and Integration of Multi-Batch Data in Metabolomics. Metabolomics 2017, 13 (10 ), 114 10.1007/s11306-017-1248-1.28890672
Redestig H. ; Fukushima A. ; Stenlund H. ; Moritz T. ; Arita M. ; Saito K. ; Kusano M. Compensation for Systematic Cross-Contribution Improves Normalization of Mass Spectrometry Based Metabolomics Data. Anal. Chem. 2009, 81 (19 ), 7974–7980. 10.1021/ac901143w.19743813
Gu Z. ; Eils R. ; Schlesner M. Complex Heatmaps Reveal Patterns and Correlations in Multidimensional Genomic Data. Bioinformatics 2016, 32 (18 ), 2847–2849. 10.1093/bioinformatics/btw313.27207943
Collaborative Data Science; Montréal, QC, Plotly Technologies Inc. 2015.
Waskom M. ; Botvinnik O. ; Hobson P. ; Cole J. B. ; Halchenko Y. ; Hoyer S. ; Miles A. ; Augspurger T. ; Yarkoni T. ; Megies T. ; Coelho L. P. ; Wehner D. ; cynddl ; Ziegler E. ; diego0020 ; Zaytsev Y. V. ; Hoppe T. ; Seabold S. ; Cloud P. ; Koskinen M. ; Meyer K. ; Qalieh A. ; Allan D. Seaborn: V0.5.0 (November 2014), Zenodo 2014.
Virtanen P. ; Gommers R. ; Oliphant T. E. ; Haberland M. ; Reddy T. ; Cournapeau D. ; Burovski E. ; Peterson P. ; Weckesser W. ; Bright J. ; van der Walt S. J. ; Brett M. ; Wilson J. ; Millman K. J. ; Mayorov N. ; Nelson A. R. J. ; Jones E. ; Kern R. ; Larson E. ; Carey C. J. ; Polat İ. ; Feng Y. ; Moore E. W. ; VanderPlas J. ; Laxalde D. ; Perktold J. ; Cimrman R. ; Henriksen I. ; Quintero E. A. ; Harris C. R. ; Archibald A. M. ; Ribeiro A. H. ; Pedregosa F. ; van Mulbregt P. ; et al. SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python. Nat. Methods 2020, 17 (3 ), 261–272. 10.1038/s41592-019-0686-2.32015543
Seabold S. ; Perktold J. Statsmodels: Econometric and Statistical Modeling with Python, Proceedings of the 9th Python in Science Conference; 2010; pp 92–96.
Vallat R. Pingouin: Statistics in Python. J. Open Source Software 2018, 3 (31 ), 1026 10.21105/joss.01026.
McInnes L. ; Healy J. ; Melville J. UMAP: Uniform Manifold Approximation and Projection for Dimension Reduction arXiv 2022. arXiv:1802.03426. arXiv.org e-Print archive. https://arxiv.org/abs/1802.03426.
Smith K. W. ; Braun J. M. ; Williams P. L. ; Ehrlich S. ; Correia K. F. ; Calafat A. M. ; Ye X. ; Ford J. ; Keller M. ; Meeker J. D. ; Hauser R. Predictors and Variability of Urinary Paraben Concentrations in Men and Women, Including before and during Pregnancy. Environ. Health Perspect. 2012, 120 (11 ), 1538–1543. 10.1289/ehp.1104614.22721761
Shu H. ; Lindh C. H. ; Wikström S. ; Bornehag C.-G. Temporal Trends and Predictors of Perfluoroalkyl Substances Serum Levels in Swedish Pregnant Women in the SELMA Study. PLoS One 2018, 13 (12 ), e0209255 10.1371/journal.pone.0209255.30596681
Wong F. ; MacLeod M. ; Mueller J. F. ; Cousins I. T. Enhanced Elimination of Perfluorooctane Sulfonic Acid by Menstruating Women: Evidence from Population-Based Pharmacokinetic Modeling. Environ. Sci. Technol. 2014, 48 (15 ), 8807–8814. 10.1021/es500796y.24943117
Berg V. ; Nøst T. H. ; Huber S. ; Rylander C. ; Hansen S. ; Veyhe A. S. ; Fuskevåg O. M. ; Odland JØ. ; Sandanger T. M. Maternal Serum Concentrations of Per- and Polyfluoroalkyl Substances and Their Predictors in Years with Reduced Production and Use. Environ. Int. 2014, 69 , 58–66. 10.1016/j.envint.2014.04.010.24815340
Varsi K. ; Huber S. ; Averina M. ; Brox J. ; Bjørke-Monsen A.-L. Quantitation of Linear and Branched Perfluoroalkane Sulfonic Acids (PFSAs) in Women and Infants during Pregnancy and Lactation. Environ. Int. 2022, 160 , 107065 10.1016/j.envint.2021.107065.34959199
Kushnir M. M. ; Rockwood A. L. ; Roberts W. L. ; Pattison E. G. ; Bunker A. M. ; Fitzgerald R. L. ; Meikle A. W. Performance Characteristics of a Novel Tandem Mass Spectrometry Assay for Serum Testosterone. Clin Chem. 2006, 52 (1 ), 120–128. 10.1373/clinchem.2005.052167.16299050
da Silva R. R. ; Dorrestein P. C. ; Quinn R. A. Illuminating the Dark Matter in Metabolomics. Proc. Natl. Acad. Sci. U.S.A. 2015, 112 (41 ), 12549–12550. 10.1073/pnas.1516878112.26430243
Dekkers K. F. ; Sayols-Baixeras S. ; Baldanzi G. ; Nowak C. ; Hammar U. ; Nguyen D. ; Varotsis G. ; Brunkwall L. ; Nielsen N. ; Eklund A. C. ; Holm J. B. ; Nielsen H. B. ; Ottosson F. ; Lin Y.-T. ; Ahmad S. ; Lind L. ; Sundström J. ; Engström G. ; Smith J. G. ; Ärnlöv J. ; Orho-Melander M. ; Fall T. An Online Atlas of Human Plasma Metabolite Signatures of Gut Microbiome Composition. Nat. Commun. 2022, 13 (1 ), 5370 10.1038/s41467-022-33050-0.36151114
Schymanski E. L. ; Jeon J. ; Gulde R. ; Fenner K. ; Ruff M. ; Singer H. P. ; Hollender J. Identifying Small Molecules via High Resolution Mass Spectrometry: Communicating Confidence. Environ. Sci. Technol. 2014, 48 (4 ), 2097–2098. 10.1021/es5002105.24476540
Koulman A. ; Volmer D. A. Perspectives for Metabolomics in Human Nutrition: An Overview. Nutr. Bull. 2008, 33 (4 ), 324–330. 10.1111/j.1467-3010.2008.00733.x.21994483
Johannessen C. ; Saini A. ; Zhang X. ; Harner T. Air Monitoring of Tire-Derived Chemicals in Global Megacities Using Passive Samplers. Environ. Pollut. 2022, 314 , 120206 10.1016/j.envpol.2022.120206.36152723
Li Z.-M. ; Kannan K. Occurrence of 1,3-Diphenylguanidine, 1,3-Di-o-Tolylguanidine, and 1,2,3-Triphenylguanidine in Indoor Dust from 11 Countries: Implications for Human Exposure. Environ. Sci. Technol. 2023, 57 (15 ), 6129–6138. 10.1021/acs.est.3c00836.37010350
Gil-Solsona R. ; Nika M.-C. ; Bustamante M. ; Villanueva C. M. ; Foraster M. ; Cosin-Tomás M. ; Alygizakis N. ; Gómez-Roig M. D. ; Llurba-Olive E. ; Sunyer J. ; Thomaidis N. S. ; Dadvand P. ; Gago-Ferrero P. The Potential of Sewage Sludge to Predict and Evaluate the Human Chemical Exposome. Environ. Sci. Technol. Lett. 2021, 8 (12 ), 1077–1084. 10.1021/acs.estlett.1c00848.35647215
Li J. ; Zeng X. ; Liang X. ; Tang S. ; Covaci A. ; Ma X. ; Yang Y. ; Chen D. Gestational Exposure to Plastic Additives and Associations with Placental Function-Related Genes. Environ. Sci. Technol. Lett. 2023, 10 (1 ), 86–92. 10.1021/acs.estlett.2c00870.
Tang S. ; Sun X. ; Qiao X. ; Cui W. ; Yu F. ; Zeng X. ; Covaci A. ; Chen D. Prenatal Exposure to Emerging Plasticizers and Synthetic Antioxidants and Their Potency to Cross Human Placenta. Environ. Sci. Technol. 2022, 56 (12 ), 8507–8517. 10.1021/acs.est.2c01141.35674357
Estlander T. ; Kostiainen M. ; Jolanki R. ; Kanerva L. Active Sensitization and Occupational Allergic Contact Dermatitis Caused by Para-Tertiary-Butylcatechol. Contact Dermatitis 1998, 38 (2 ), 96–100. 10.1111/j.1600-0536.1998.tb05659.x.9506222
Du P. ; Melián-Cabrera I. ; van der Linden B. ; Jansma H. ; Moulijn J. A. Reactant Additive-Triggered Deactivation of Pd/γ-Alumina-Catalyzed Hydrogenation Reactions. A Reactivity and Adsorption Study. Ind. Eng. Chem. Res. 2020, 59 (40 ), 17762–17768. 10.1021/acs.iecr.0c03238.
Liu Y. ; Zhu T. ; Xie Z. ; Deng C. ; Qi X. ; Hu R. ; Wang J. ; Chen J. Human Exposure to Chlorinated Organophosphate Ester Flame Retardants and Plasticizers in an Industrial Area of Shenzhen, China. Int. J. Environ. Res. Public Health 2022, 19 (5 ), 3126 10.3390/ijerph19053126.35270815
Faiz Y. ; Zhao W. ; Feng J. ; Sun C. ; He H. ; Zhu J. Occurrence of Triphenylphosphine Oxide and Other Organophosphorus Compounds in Indoor Air and Settled Dust of an Institute Building. Build. Environ. 2016, 106 , 196–204. 10.1016/j.buildenv.2016.06.022.
Besis A. ; Avgenikou A. ; Pantelaki I. ; Serafeim E. ; Georgiadou E. ; Voutsa D. ; Samara C. Hazardous Organic Pollutants in Indoor Dust from Elementary Schools and Kindergartens in Greece: Implications for Children’s Health. Chemosphere 2023, 310 , 136750 10.1016/j.chemosphere.2022.136750.36241110
Mehra S. ; Chadha P. Bioaccumulation and Toxicity of 2-Naphthalene Sulfonate: An Intermediate Compound Used in Textile Industry. Toxicol. Res. 2020, 9 (2 ), 127–136. 10.1093/toxres/tfaa008.
Hashemi S. H. ; Kaykhaii M. Developments in Methods of Analysis for Naphthalene Sulfonates. Crit. Rev. Anal. Chem. 2017, 47 (2 ), 127–137. 10.1080/10408347.2016.1225256.27541207
Bangee O. D. ; Wilson V. H. ; East G. C. ; Holme I. Antioxidant-Induced Yellowing of Textiles. Polym. Degrad. Stab. 1995, 50 (3 ), 313–317. 10.1016/0141-3910(95)00156-5.
Grigoriadou A. ; Schwarzbauer J. Non-Target Screening of Organic Contaminants in Sediments from the Industrial Coastal Area of Kavala City (NE Greece). Water, Air, Soil Pollut. 2011, 214 (1 ), 623–643. 10.1007/s11270-010-0451-8.
Still K. R. ; Jung A. E. ; Ritchie G. D. ; Jederberg W. W. ; Wilfong E. R. ; Briggs G. B. ; Arfsten D. P. Disposition of 2,6-Di-Tert-Butyl-4-Nitrophenol (DBNP), a Submarine Atmosphere Contaminant, in Male Sprague–Dawley Rats. Environ. Res. 2005, 98 (3 ), 363–367. 10.1016/j.envres.2004.08.009.15910791
Lestido-Cardama A. ; de Quirós A. R. B. ; Bustos J. ; Lomo M. L. ; Losada P. P. ; Sendón R. Estimation of Dietary Exposure to Contaminants Transferred from the Packaging in Fatty Dry Foods Based on Cereals. Foods 2020, 9 (8 ), 1038 10.3390/foods9081038.32752291
Zhu Q. ; Liao C. ; Jiang G. Occurrence of Human Exposure to Benzothiazoles and Benzotriazoles in Indoor Dust in Suizhou and Beijing, China. Chem. Res. Chin. Univ. 2023, 39 (3 ), 508–515. 10.1007/s40242-023-3062-9.
Zhao J.-H. ; Hu L.-X. ; Xiao S. ; Zhao J.-L. ; Liu Y.-S. ; Yang B. ; Zhang Q.-Q. ; Ying G.-G. Screening and Prioritization of Organic Chemicals in a Large River Basin by Suspect and Non-Target Analysis. Environ. Pollut. 2023, 333 , 122098 10.1016/j.envpol.2023.122098.37352960
Bonnefille B. ; Karlsson O. ; Rian M. B. ; Raqib R. ; Parvez F. ; Papazian S. ; Islam M. S. ; Martin J. W. Nontarget Analysis of Polluted Surface Waters in Bangladesh Using Open Science Workflows. Environ. Sci. Technol. 2023, 57 (17 ), 6808–6824. 10.1021/acs.est.2c08200.37083417
Li J. ; Zhao H. ; Zhou Y. ; Xu S. ; Cai Z. Determination of Benzotriazoles and Benzothiazoles in Human Urine by UHPLC-TQMS. J. Chromatogr. B 2017, 1070 , 70–75. 10.1016/j.jchromb.2017.10.045.
Li X. ; Wang L. ; Asimakopoulos A. G. ; Sun H. ; Zhao Z. ; Zhang J. ; Zhang L. ; Wang Q. Benzotriazoles and Benzothiazoles in Paired Maternal Urine and Amniotic Fluid Samples from Tianjin, China. Chemosphere 2018, 199 , 524–530. 10.1016/j.chemosphere.2018.02.076.29455122
Krais A. M. ; van Wendel de Joode B. ; Liljedahl E. R. ; Blomberg A. J. ; Rönnholm A. ; Bengtsson M. ; Cano J. C. ; Hoppin J. A. ; Littorin M. ; Nielsen C. ; Lindh C. H. Detection of the Fungicide Transformation Product 4-Hydroxychlorothalonil in Serum of Pregnant Women from Sweden and Costa Rica. J. Exposure Sci. Environ. Epidemiol. 2023, 34 , 270–277. 10.1038/s41370-023-00580-8.
Wällstedt E. Temporal Trends of 4-Hydroxychlorothalonil in Maternal Serum Samples. 2023.
Kiefer K. ; Müller A. ; Singer H. ; Hollender J. New Relevant Pesticide Transformation Products in Groundwater Detected Using Target and Suspect Screening for Agricultural and Urban Micropollutants with LC-HRMS. Water Res. 2019, 165 , 114972 10.1016/j.watres.2019.114972.31450217
Ashley-Martin J. ; VanLeeuwen J. ; Cribb A. ; Andreou P. ; Guernsey J. R. Breast Cancer Risk, Fungicide Exposure and CYP1A1*2A Gene-Environment Interactions in a Province-Wide Case Control Study in Prince Edward Island, Canada. Int. J. Environ. Res. Public Health 2012, 9 (5 ), 1846–1858. 10.3390/ijerph9051846.22754477
Chaves A. ; Shea D. ; Cope W. G. Environmental Fate of Chlorothalonil in a Costa Rican Banana Plantation. Chemosphere 2007, 69 (7 ), 1166–1174. 10.1016/j.chemosphere.2007.03.048.17498771
Fatimi A. Sodium Tetradecyl Sulfate Molecule: Patent Analysis Based on Chemical Compounds Search, The 1st International Electronic Conference on Processes: Processes System Innovation; MDPI, 2022; p 4.
Atrux-Tallau N. ; Romagny C. ; Padois K. ; Denis A. ; Haftek M. ; Falson F. ; Pirot F. ; Maibach H. I. Effects of Glycerol on Human Skin Damaged by Acute Sodium Lauryl Sulphate Treatment. Arch. Dermatol. Res. 2010, 302 (6 ), 435–441. 10.1007/s00403-009-1021-z.20043170
Bondi C. A. M. ; Marks J. L. ; Wroblewski L. B. ; Raatikainen H. S. ; Lenox S. R. ; Gebhardt K. E. Human and Environmental Toxicity of Sodium Lauryl Sulfate (SLS): Evidence for Safe Use in Household Cleaning Products. Environ. Health Insights 2015, 9 , EHI-S31765 10.4137/EHI.S31765.
Kim Y.-H. ; Lee K. Characterization of Aerosols Produced during Shampoo Use and Harmful Chemicals in Shampoo Aerosols. Environ. Res. 2022, 204 , 111957 10.1016/j.envres.2021.111957.34478728
Asio J. R. G. ; Garcia J. S. ; Antonatos C. ; Sevilla-Nastor J. B. ; Trinidad L. C. Sodium Lauryl Sulfate and Its Potential Impacts on Organisms and the Environment: A Thematic Analysis. Emerging Contam. 2023, 9 (1 ), 100205 10.1016/j.emcon.2023.100205.
Roggeman M. ; Gys C. ; Klimowska A. ; Bastiaensen M. ; Wielgomas B. ; Bamai Y. A. ; Covaci A. Reviewing the Variability in Urinary Concentrations of Non-Persistent Organic Chemicals: Evaluation across Classes, Sampling Strategies and Dilution Corrections. Environ. Res. 2022, 215 , 114332 10.1016/j.envres.2022.114332.36116496
Liljequist D. ; Elfving B. ; Roaldsen K. S. Intraclass Correlation – A Discussion and Demonstration of Basic Features. PLoS One 2019, 14 (7 ), e0219854 10.1371/journal.pone.0219854.31329615
Casas M. ; Basagaña X. ; Sakhi A. K. ; Haug L. S. ; Philippat C. ; Granum B. ; Manzano-Salgado C. B. ; Brochot C. ; Zeman F. ; de Bont J. ; Andrusaityte S. ; Chatzi L. ; Donaire-Gonzalez D. ; Giorgis-Allemand L. ; Gonzalez J. R. ; Gracia-Lavedan E. ; Grazuleviciene R. ; Kampouri M. ; Lyon-Caen S. ; Pañella P. ; Petraviciene I. ; Robinson O. ; Urquiza J. ; Vafeiadi M. ; Vernet C. ; Waiblinger D. ; Wright J. ; Thomsen C. ; Slama R. ; Vrijheid M. Variability of Urinary Concentrations of Non-Persistent Chemicals in Pregnant Women and School-Aged Children. Environ. Int. 2018, 121 , 561–573. 10.1016/j.envint.2018.09.046.30300814
Vermeulen R. ; Schymanski E. L. ; Barabási A.-L. ; Miller G. W. The Exposome and Health: Where Chemistry Meets Biology. Science 2020, 367 (6476 ), 392–396. 10.1126/science.aay3164.31974245
Yang L. ; Wang Z. ; Shi Y. ; Li J. ; Wang Y. ; Zhao Y. ; Wu Y. ; Cai Z. Human Placental Transfer of Perfluoroalkyl Acid Precursors: Levels and Profiles in Paired Maternal and Cord Serum. Chemosphere 2016, 144 , 1631–1638. 10.1016/j.chemosphere.2015.10.063.26517392
Bezchlibnyk-Butler K. ; Aleksic I. ; Kennedy S. H. Citalopram--a Review of Pharmacological and Clinical Effects. J. Psychiatry Neurosci. 2000, 25 (3 ), 241–254.10863884
Olsen G. W. ; Burris J. M. ; Ehresman D. J. ; Froehlich J. W. ; Seacat A. M. ; Butenhoff J. L. ; Zobel L. R. Half-Life of Serum Elimination of Perfluorooctanesulfonate,Perfluorohexanesulfonate, and Perfluorooctanoate in Retired Fluorochemical Production Workers. Environ. Health Perspect. 2007, 115 (9 ), 1298–1305. 10.1289/ehp.10009.17805419
Benowitz N. L. ; Jacob P. Trans-3′-Hydroxycotinine: Disposition Kinetics, Effects and Plasma Levels during Cigarette Smoking. Br. J. Clin. Pharmacol. 2001, 51 (1 ), 53–59. 10.1046/j.1365-2125.2001.01309.x.11167665
Uhl S. ; Schmid P. ; Schlatter C. Pharmacokinetics of Pentachlorophenol in Man. Arch. Toxicol. 1986, 58 (3 ), 182–186. 10.1007/BF00340979.3964082
Institute of Medicine (US) Committee on Military Nutrition Research. Pharmacology of Caffeine National Academies Press: (US); 2001.
Norén E. ; Lindh C. ; Rylander L. ; Glynn A. ; Axelsson J. ; Littorin M. ; Faniband M. ; Larsson E. ; Nielsen C. Concentrations and Temporal Trends in Pesticide Biomarkers in Urine of Swedish Adolescents, 2000–2017. J. Exposure Sci. Environ. Epidemiol. 2020, 30 (4 ), 756–767. 10.1038/s41370-020-0212-8.
Baggott M. J. ; Childs E. ; Hart A. B. ; de Bruin E. ; Palmer A. A. ; Wilkinson J. E. ; de Wit H. Psychopharmacology of Theobromine in Healthy Volunteers. Psychopharmacology 2013, 228 (1 ), 109–118. 10.1007/s00213-013-3021-0.23420115
Abafe O. A. ; Lawal M. A. ; Chokwe T. B. Non-Targeted Screening of Emerging Contaminants in South African Surface and Wastewater. Emerging Contam. 2023, 9 (4 ), 100246 10.1016/j.emcon.2023.100246.
Holtze M. S. ; Sørensen S. R. ; Sørensen J. ; Aamand J. Microbial Degradation of the Benzonitrile Herbicides Dichlobenil, Bromoxynil and Ioxynil in Soil and Subsurface Environments – Insights into Degradation Pathways, Persistent Metabolites and Involved Degrader Organisms. Environ. Pollut. 2008, 154 (2 ), 155–168. 10.1016/j.envpol.2007.09.020.17988770
Jenkins B. ; West J. A. ; Koulman A. A Review of Odd-Chain Fatty Acid Metabolism and the Role of Pentadecanoic Acid (C15:0) and Heptadecanoic Acid (C17:0) in Health and Disease. Molecules 2015, 20 (2 ), 2425–2444. 10.3390/molecules20022425.25647578
Carrizo D. ; Grimalt J. O. ; Ribas-Fito N. ; Torrent M. ; Sunyer J. Pentachlorobenzene, Hexachlorobenzene, and Pentachlorophenol in Children’s Serum from Industrial and Rural Populations after Restricted Use. Ecotoxicol. Environ. Saf. 2008, 71 (1 ), 260–266. 10.1016/j.ecoenv.2007.08.021.17935782
Zamaratskaia G. ; Bjermo H. ; Lindroos A. K. ; Ruokojärvi P. ; Rantakokko P. ; Kiviranta H. ; Lignell S. Breastfeeding during Infancy and Consumption of Fish and Dairy Products Are Associated with Chlorinated Persistent Organic Pollutants in Serum from Swedish Adolescents. Environ. Adv. 2022, 8 , 100210 10.1016/j.envadv.2022.100210.
Shi M. ; Zheng X. ; Ge Y. ; Zhang N. ; Yu L. ; Duan X. ; Liu Y. ; Xue H. ; You J. ; Yin L. Unraveling the Cytotoxicity and Cellular Uptake of Low, Medium and High Molecular Weight Polyethylene Glycol Polymers in MCF-7 Cells by Green UPLC-MS/MS Methods. J. Pharm. Biomed. Anal. 2024, 238 , 115868 10.1016/j.jpba.2023.115868.38000190
Pozza S. A. ; Gonçalves P. B. ; Wouters F. C. ; Vendemiatti J. A. S. ; Nogarotto D. C. ; Pereira-Filho E. R. ; Osório D. M. M. ; Romualdo L. L. ; Godoi J. R. ; Hoinaski L. ; Urban R. C. Particulate Matter Pollution and Non-Targeted Analysis of Polar Compounds in Three Regions of Brazil. Chemosphere 2023, 341 , 139839 10.1016/j.chemosphere.2023.139839.37604347
Eysseric E. ; Beaudry F. ; Gagnon C. ; Segura P. A. Non-Targeted Screening of Trace Organic Contaminants in Surface Waters by a Multi-Tool Approach Based on Combinatorial Analysis of Tandem Mass Spectra and Open Access Databases. Talanta 2021, 230 , 122293 10.1016/j.talanta.2021.122293.33934765
Arvaniti O. S. ; Gkotsis G. ; Nika M.-C. ; Gyparakis S. ; Manios T. ; Thomaidis N. S. ; Fountoulakis M. S. ; Stasinakis A. S. Study on the Occurrence of Artificial Sweeteners, Parabens, and Other Emerging Contaminants in Hospital Wastewater Using LC-QToF-MS Target Screening Approach. Water 2023, 15 (5 ), 936 10.3390/w15050936.
Lelo A. ; Miners J. O. ; Robson R. ; Birkett D. J. Assessment of Caffeine Exposure: Caffeine Content of Beverages, Caffeine Intake, and Plasma Concentrations of Methylxanthines. Clin. Pharmacol. Ther. 1986, 39 (1 ), 54–59. 10.1038/clpt.1986.10.3943270
Müller C. ; Lang R. ; Hofmann T. Quantitative Precursor Studies on Di- and Trihydroxybenzene Formation during Coffee Roasting Using “In Bean” Model Experiments and Stable Isotope Dilution Analysis. J. Agric. Food Chem. 2006, 54 (26 ), 10086–10091. 10.1021/jf062727y.17177545
Jacob P. ; Yu L. ; Duan M. ; Ramos L. ; Yturralde O. ; Benowitz N. L. Determination of the Nicotine Metabolites Cotinine and Trans-3′-Hydroxycotinine in Biologic Fluids of Smokers and Non-Smokers Using Liquid Chromatography - Tandem Mass Spectrometry: Biomarkers for Tobacco Smoke Exposure and for Phenotyping Cytochrome P450 2A6 Activity. J. Chromatogr. B 2011, 879 (3–4 ), 267–276. 10.1016/j.jchromb.2010.12.012.
Sköld O. Resistance to Trimethoprim and Sulfonamides. Vet. Res. 2001, 32 (3–4 ), 261–273. 10.1051/vetres:2001123.11432417
Kantor T. G. Ketoprofen: A Review of Its Pharmacologic and Clinical Properties. Pharmacother., J. Hum. Pharmacol. Drug Ther. 1986, 6 (3 ), 93–102. 10.1002/j.1875-9114.1986.tb03459.x.
Boye B. ; Dieng M. M. ; Brillas E. Degradation of Herbicide 4-Chlorophenoxyacetic Acid by Advanced Electrochemical Oxidation Methods. Environ. Sci. Technol. 2002, 36 (13 ), 3030–3035. 10.1021/es0103391.12144282
Nolte J. ; Heimlich F. ; Graß B. ; Zullei-Seibert N. ; Preuss G. Studies on the Behaviour of Dihalogenated Hydroxybenzonitriles in Water. Fresenius J. Anal Chem. 1995, 351 (1 ), 88–91. 10.1007/BF00324296.
Jackson K. D. ; Achour B. ; Lee J. ; Geffert R. M. ; Beers J. L. ; Latham B. D. Novel Approaches to Characterize Individual Drug Metabolism and Advance Precision Medicine. Drug Metab. Dispos. 2023, 51 (10 ), 1238–1253. 10.1124/dmd.122.001066.37419681
Stülten D. ; Zühlke S. ; Lamshöft M. ; Spiteller M. Occurrence of Diclofenac and Selected Metabolites in Sewage Effluents. Sci. Total Environ. 2008, 405 (1 ), 310–316. 10.1016/j.scitotenv.2008.05.036.18640705
Snowman A. M. ; Snyder S. H. Cetirizine: Actions on Neurotransmitter Receptors. J. Allergy Clin. Immunol. 1990, 86 (6, Part 2 ), 1025–1028. 10.1016/S0091-6749(05)80248-9.1979798
Patel J. ; Spencer E. P. ; Flanagan R. J. HPLC of Sertraline and Norsertraline in Plasma or Serum. Biomed. Chromatogr. 1996, 10 (6 ), 351–354. 10.1002/(SICI)1099-0801(199611)10:6<351::AID-BMC620>3.0.CO;2-O.8949919
Luo L. ; Aubrecht J. ; Li D. ; Warner R. L. ; Johnson K. J. ; Kenny J. ; Colangelo J. L. Assessment of Serum Bile Acid Profiles as Biomarkers of Liver Injury and Liver Disease in Humans. PLoS One 2018, 13 (3 ), e0193824 10.1371/journal.pone.0193824.29513725
Mair R. D. ; Sirich T. L. ; Meyer T. W. Uremic Toxin Clearance and Cardiovascular Toxicities. Toxins 2018, 10 (6 ), 226 10.3390/toxins10060226.29865226
Wang Y. ; Aimuzi R. ; Nian M. ; Zhang Y. ; Luo K. ; Zhang J. Perfluoroalkyl Substances and Sex Hormones in Postmenopausal Women: NHANES 2013–2016. Environ. Int. 2021, 149 , 106408 10.1016/j.envint.2021.106408.33548847
Endogenous Sex Hormones and Breast Cancer in Postmenopausal Women: Reanalysis of Nine Prospective Studies. J. Natl. Cancer Inst. 2002, 94 (8 ), 606–616. 10.1093/jnci/94.8.606.11959894
Zhang P. ; Carlsten C. ; Chaleckis R. ; Hanhineva K. ; Huang M. ; Isobe T. ; Koistinen V. M. ; Meister I. ; Papazian S. ; Sdougkou K. ; Xie H. ; Martin J. W. ; Rappaport S. M. ; Tsugawa H. ; Walker D. I. ; Woodruff T. J. ; Wright R. O. ; Wheelock C. E. Defining the Scope of Exposome Studies and Research Needs from a Multidisciplinary Perspective. Environ. Sci. Technol. Lett. 2021, 8 (10 ), 839–852. 10.1021/acs.estlett.1c00648.34660833
