Emergence and spread of SARS-CoV-2 variants of concern in Canada: a retrospective analysis from clinical and wastewater data.

BACKGROUND: The spread of SARS-CoV-2 has been studied at unprecedented levels worldwide. In jurisdictions where molecular analysis was performed on large scales, the emergence and competition of numerous SARS-CoV-2lineages have been observed in near real-time. Lineage identification, traditionally performed from clinical samples, can also be determined by sampling wastewater from sewersheds serving populations of interest. Variants of concern (VOCs) and SARS-CoV-2 lineages associated with increased transmissibility and/or severity are of particular interest. METHOD: Here, we consider clinical and wastewater data sources to assess the emergence and spread of VOCs in Canada retrospectively. RESULTS: We show that, overall, wastewater-based VOC identification provides similar insights to the surveillance based on clinical samples. Based on clinical data, we observed synchrony in VOC introduction as well as similar emergence speeds across most Canadian provinces despite the large geographical size of the country and differences in provincial public health measures. CONCLUSION: In particular, it took approximately four months for VOC Alpha and Delta to contribute to half of the incidence. In contrast, VOC Omicron achieved the same contribution in less than one month. This study provides significant benchmarks to enhance planning for future VOCs, and to some extent for future pandemics caused by other pathogens, by quantifying the rate of SARS-CoV-2 VOCs invasion in Canada.

Acute Exposure of Zebrafish (Danio rerio) to the Next-Generation Perfluoroalkyl Substance, Perfluoroethylcyclohexanesulfonate, Shows Similar Effects as Legacy Substances.

Perfluoroethylcyclohexanesulfonate (PFECHS) is an emerging perfluoroalkyl substance (PFAS) that has been considered a potential replacement for perfluorooctanesulfonic acid (PFOS). However, there is little information characterizing the toxic potency of PFECHS to zebrafish embryos and its potential for effects in aquatic environments. This study assessed toxic potency of PFECHS in vivo during both acute (96-hour postfertilization) and chronic (21-day posthatch) exposures and tested concentrations of PFECHS from 500 ng/L to 2 mg/L. PFECHS was less likely to cause mortalities than PFOS for both the acute and chronic experiments based on previously published values for PFOS exposure, but exposure resulted in a similar incidence of deformities. Exposure to PFECHS also resulted in significantly increased abundance of transcripts of peroxisome proliferator activated receptor alpha (pparα), cytochrome p450 1a1 (cyp1a1), and apolipoprotein IV (apoaIV) at concentrations nearing those of environmental relevance. Overall, these results provide further insight into the safety of an emerging PFAS alternative in the aquatic environment and raise awareness that previously considered "safer" alternatives may show similar effects as legacy PFASs.

Interspecies Differences in 6PPD-Quinone Toxicity Across Seven Fish Species: Metabolite Identification and Semiquantification.

N-(1,3-Dimethylbutyl)-N'-phenyl-p-phenylenediamine-quinone (6PPD-Q) is a recently identified contaminant that originates from the oxidation of the tire antidegradant 6PPD. 6PPD-Q is acutely toxic to select salmonids at environmentally relevant concentrations, while other fish species display tolerance to concentrations that surpass those measured in the environment. The reasons for these marked differences in sensitivity are presently unknown. The objective of this research was to explore potential toxicokinetic drivers of species sensitivity by characterizing biliary metabolites of 6PPD-Q in sensitive and tolerant fishes. For the first time, we identified an O-glucuronide metabolite of 6PPD-Q using high-resolution mass spectrometry. The semiquantified levels of this metabolite in tolerant species or life stages, including white sturgeon (Acipenser transmontanus), chinook salmon (Oncorhynchus tshawytscha), westslope cutthroat trout (Oncorhynchus clarkii lewisi), and nonfry life stages of Atlantic salmon (Salmo salar), were greater than those in sensitive species, including coho salmon (Oncorhynchus kisutch), brook trout (Salvelinus fontinalis), and rainbow trout (Oncorhynchus mykiss), suggesting that tolerant species might detoxify 6PPD-Q more effectively. Thus, we hypothesize that differences in species sensitivity are a result of differences in basal expression of biotransformation enzyme across various fish species. Moreover, the semiquantification of 6PPD-Q metabolites in bile extracted from wild-caught fish might be a useful biomarker of exposure to 6PPD-Q, thereby being valuable to environmental monitoring and risk assessment.

Validation of Methods for in Vitro-in Vivo Extrapolation Using Hepatic Clearance Measurements in Isolated Perfused Fish Livers.

In vitro biotransformation assays using hepatocytes or liver subcellular fractions, combined with in vitro-in vivo extrapolation (IVIVE) models, have been proposed as an alternative to live fish bioconcentration studies. The uncertainty associated with IVIVE approaches to date has been attributed to assay protocols, model assumptions, or variability of in vivo data. An isolated perfused trout liver model that measures hepatic clearance has been proposed for validating IVIVE predictions in the absence of other confounding factors. Here, we investigated the hepatic clearances of five chemicals (pyrene, phenanthrene, 4-n-nonlyphenol, deltamethrin, and methoxychlor) in this model and compared measured rates to values predicted from published in vitro intrinsic clearances for validation of IVIVE models. Additionally, we varied protein concentrations in perfusates to test binding assumptions of these models. We found that measured and predicted hepatic clearances were in very good agreement (root mean squared error 16.8 mL h-1 g-1) across three levels of protein binding and across a more diverse chemical space than previously studied within this system. Our results show that current IVIVE methods can reliably predict in vivo clearance rates and indicate that discrepancies from measured bioconcentration factors might be driven by other processes, such as extrahepatic biotransformation, etc., and help streamline optimization efforts to the processes that truly matter.

Inflammation of Gill Epithelia in Fish Causes Increased Permeation of Petrogenic Polar Organic Chemicals via Disruption of Tight Junctions.

The epithelial cell layer that lines the gills of fish controls paracellular permeation of chemicals through tight junctions. The integrity of tight junctions can be affected by inflammation, which likely affects the bioavailability of chemicals. Here, the inflammation of the rainbow trout gill cell line RTgill-W1 was induced via exposure to bacterial lipopolysaccharides (LPS). Cells were then coexposed to extracts of oil sands process-affected water (OSPW), which contain complex mixtures of chemicals. After 24 h of exposure, cells exposed to LPS showed a reduction in transepithelial electrical resistance, an indicator of tight junction integrity. Quantitative reverse-transcription polymerase chain reaction (RT-PCR) analysis determined that abundances of transcripts of genes coding for tight junction proteins were significantly less in cells exposed to 20, 50, or 100 mg L-1 LPS. Chemical analysis revealed increased permeation of constituents of OSPW across epithelia at all studied LPS concentrations. These in vitro findings were confirmed in vivo in rainbow trout exposed to LPS and OSPW for 48 h, which resulted in greater accumulation of chemicals relative to that for fish exposed to OSPW alone. Our results demonstrated that inflammation and disruption of tight junctions could lead to greater uptake of potentially harmful chemicals from the environment, which has implications for mixture risk assessment.

Predicting Hepatic Clearance of Psychotropic Drugs in Isolated Perfused Fish Livers Using a Combination of Two In Vitro Assays.

In vitro biotransformation assays with primary trout hepatocytes (RT-HEP) or liver subcellular fractions (RT-S9) have been proposed as valuable tools to help scientists and regulators better understand the toxicokinetics of chemicals. While both assays have been applied successfully to a diversity of neutral organic chemicals, only the RT-S9 assay has been applied to a large number of ionizable organic chemicals. Here, a combination of an in vitro biotransformation assay with RT-HEP with an active transport assay based on the permanent rainbow trout liver cell line RTL-W1 was used to qualitatively predict the potential hepatic clearance of nine psychotropic drugs with various degrees of ionization. Predictions were compared with rates of clearance measured in isolated perfused rainbow trout livers, and the importance of active transport was verified in the presence of the active transport inhibitor cyclosporin A. For the first time, it was demonstrated that a combination of biotransformation and active transport assays is powerful for the prediction of rates of hepatic clearance of ionizable chemicals. Ultimately, it is expected that this approach will allow for use of fewer animals while at the same time improving our confidence in the use of data from in vitro assays in chemical risk assessment.

Effects of in situ experimental selenium exposure on finescale dace (Phoxinus neogaeus) gut microbiome.

Selenium (Se) is an environmental contaminant of global concern that can cause adverse effects in fish at elevated levels. Fish gut microbiome play essential roles in gastrointestinal function and host health and can be perturbed by environmental contaminants, including metals and metalloids. Here, an in-situ Se exposure of female finescale dace (Phoxinus neogaeus) using mesocosms was conducted to determine the impacts of Se accumulation on the gut microbiome and morphometric endpoints. Prior to this study, the gut microbiome of finescale dace, a widespread Cyprinid throughout North America, had not been characterized. Exposure to Se caused a hormetic response of alpha diversity of the gut microbiome, with greater diversity at the lesser concentration of 1.6 µg Se/L, relative to that of fish exposed to the greater concentration of 5.6 µg Se/L. Select gut microbiome taxa of fish were differentially abundant between aqueous exposure concentrations and significantly correlated with liver-somatic index (LSI). The potential effects of gut microbiome dysbiosis on condition of wild fish might be a consideration when assessing adverse effects of Se in aquatic environments. More research regarding effects of Se on field-collected fish gut microbiome and the potential adverse effects or benefits on the host is warranted.

Bone resorption and body reorganization during maturation induce maternal transfer of toxic metals in anguillid eels.

During their once-in-a-lifetime transoceanic spawning migration, anguillid eels do not feed, instead rely on energy stores to fuel the demands of locomotion and reproduction while they reorganize their bodies by depleting body reserves and building up gonadal tissue. Here we show how the European eel (Anguilla anguilla) breaks down its skeleton to redistribute phosphorus and calcium from hard to soft tissues during its sexual development. Using multiple analytical and imaging techniques, we characterize the spatial and temporal degradation of the skeletal framework from initial to final gonadal maturation and use elemental mass ratios in bone, muscle, liver, and gonadal tissue to determine the fluxes and fates of selected minerals and metals in the eels' bodies. We find that bone loss is more pronounced in females than in males and eventually may reach a point at which the mechanical stability of the skeleton is challenged. P and Ca are released and translocated from skeletal tissues to muscle and gonads, leaving both elements in constant proportion in remaining bone structures. The depletion of internal stores from hard and soft tissues during maturation-induced body reorganization is accompanied by the recirculation, translocation, and maternal transfer of potentially toxic metals from bone and muscle to the ovaries in gravid females, which may have direct deleterious effects on health and hinder the reproductive success of individuals of this critically endangered species.

Toxicity of Weathered Sediment-Bound Dilbit to Early Life Stages of Zebrafish (Danio rerio).

Due to high viscosity, bitumen extracted from the Alberta oil sands is diluted with natural gas condensates to form diluted bitumen (dilbit) to facilitate transport through pipelines. Dilbit that is spilled into or near a waterbody is subject to environmental weathering processes such as evaporation and interaction with sediments. This is the first study that assessed the toxicity of weathered sediment-bound dilbit (WSD) to fish early life stages. Exposure of zebrafish (Danio rerio) embryos to water-soluble fractions (WSFs) or water-accommodated fractions (WAFs) of WSD from 30 min to 120 h postfertilization resulted in pericardial edema, yolk sac edema, and incidences of uninflated swim bladder. The presence of oil-mineral aggregates (OMAs) in the WAFs greatly increased toxicity, despite all fractions having similar concentrations of dissolved polycyclic aromatic hydrocarbons (PAHs). There were greater cyp1a mRNA abundances in larvae exposed to WAFs, suggesting that there were differences in bioavailability of PAHs between fractions. However, there was little evidence that embryotoxicity was caused by oxidative stress. Results suggest that evaporation and sediment interaction do not completely attenuate toxicity of dilbit to zebrafish early life stages, and OMAs in exposures exacerbate toxicity.

A Novel Multispecies Toxicokinetic Modeling Approach in Support of Chemical Risk Assessment.

Standardized laboratory tests with a limited number of model species are a key component of chemical risk assessments. These surrogate species cannot represent the entire diversity of native species, but there are practical and ethical objections against testing chemicals in a large variety of species. In previous research, we have developed a multispecies toxicokinetic model to extrapolate chemical bioconcentration across species by combining single-species physiologically based toxicokinetic (PBTK) models. This "top-down" approach was limited, however, by the availability of fully parameterized single-species models. Here, we present a "bottom-up" multispecies PBTK model based on available data from 69 freshwater fishes found in Canada. Monte Carlo-like simulations were performed using statistical distributions of model parameters derived from these data to predict steady-state bioconcentration factors (BCFs) for a set of well-studied chemicals. The distributions of predicted BCFs for 1,4-dichlorobenzene and dichlorodiphenyltrichloroethane largely overlapped those of empirical data, although a tendency existed toward overestimation of measured values. When expressed as means, predicted BCFs for 26 of 34 chemicals (82%) deviated by less than 10-fold from measured data, indicating an accuracy similar to that of previously published single-species models. This new model potentially enables more environmentally relevant predictions of bioconcentration in support of chemical risk assessments.

Toxicokinetic Models for Bioconcentration of Organic Contaminants in Two Life Stages of White Sturgeon (Acipenser transmontanus).

The white sturgeon (Acipenser transmontanus) is an endangered ancient fish species that is known to be particularly sensitive to certain environmental contaminants, partly because of the uptake and subsequent toxicity of lipophilic pollutants prone to bioconcentration as a result of their high lipid content. To better understand the bioconcentration of organic contaminants in this species, toxicokinetic (TK) models were developed for the embryo-larval and subadult life stages. The embryo-larval model was designed as a one-compartment model and validated using whole-body measurements of benzo[a]pyrene (B[a]P) metabolites from a waterborne exposure to B[a]P. A physiologically based TK (PBTK) model was used for the subadult model. The predictive power of the subadult model was validated with an experimental data set of four chemicals. Results showed that the TK models could accurately predict the bioconcentration of organic contaminants for both life stages of white sturgeon within 1 order of magnitude of measured values. These models provide a tool to better understand the impact of environmental contaminants on the health and the survival of endangered white sturgeon populations.

Development of a Comprehensive Toxicity Pathway Model for 17α-Ethinylestradiol in Early Life Stage Fathead Minnows (Pimephales promelas).

There is increasing pressure to develop alternative ecotoxicological risk assessment approaches that do not rely on expensive, time-consuming, and ethically questionable live animal testing. This study aimed to develop a comprehensive early life stage toxicity pathway model for the exposure of fish to estrogenic chemicals that is rooted in mechanistic toxicology. Embryo-larval fathead minnows (FHM; Pimephales promelas) were exposed to graded concentrations of 17α-ethinylestradiol (water control, 0.01% DMSO, 4, 20, and 100 ng/L) for 32 days. Fish were assessed for transcriptomic and proteomic responses at 4 days post-hatch (dph), and for histological and apical end points at 28 dph. Molecular analyses revealed core responses that were indicative of observed apical outcomes, including biological processes resulting in overproduction of vitellogenin and impairment of visual development. Histological observations indicated accumulation of proteinaceous fluid in liver and kidney tissues, energy depletion, and delayed or suppressed gonad development. Additionally, fish in the 100 ng/L treatment group were smaller than controls. Integration of omics data improved the interpretation of perturbations in early life stage FHM, providing evidence of conservation of toxicity pathways across levels of biological organization. Overall, the mechanism-based embryo-larval FHM model showed promise as a replacement for standard adult live animal tests.

Mechanisms of pH-Dependent Uptake of Ionizable Organic Chemicals by Fish from Oil Sands Process-Affected Water (OSPW).

Uptake and effects of ionizable organic chemicals (IOCs) that are weak acids in aqueous solution by fish can differ as a function of pH. While the pH-dependent behavior of select IOCs is well-understood, complex mixtures of IOCs, e.g., from oil sands process-affected water (OSPW), have not yet been studied systematically. Here, we established an in vitro screening method using the rainbow trout gill cell line, RTgill-W1, to investigate pH-dependent cytotoxicity and permeation of IOCs across cultured epithelia using ultra-high-performance liquid chromatography with high-resolution mass spectrometry (UPLC-HRMS). The assay was benchmarked using model chemicals and technical mixtures, and then used to characterize fractions and reconstituted extracts of field-collected OSPW. Significant pH-dependent cytotoxicity of individual IOCs, acidic fractions, and reconstituted extracts of OSPW was observed. In vitro data were in good agreement with data from a 96 h in vivo exposure experiment with juvenile rainbow trout. Permeation of some IOCs from OSPW was mediated by active transport, as revealed by studies in which inhibitors of these active transport mechanisms were applied. We conclude that the RTgill-W1 in vitro assay is useful for the screening of pH-dependent uptake of IOCs in fish, and has applications for in vitro-in vivo extrapolation, and prioritization of chemicals in nontarget screenings.

Bioactivation of Quinolines in a Recombinant Estrogen Receptor Transactivation Assay Is Catalyzed by N-Methyltransferases.

Hydroxylation of polyaromatic compounds through cytochromes P450 (CYPs) is known to result in potentially estrogenic transformation products. Recently, there has been an increasing awareness of the importance of alternative pathways such as aldehyde oxidases (AOX) or N-methyltransferases (NMT) in bioactivation of small molecules, particularly N-heterocycles. Therefore, this study investigated the biotransformation and activity of methylated quinolines, a class of environmentally relevant N-heterocycles that are no native ligands of the estrogen receptor (ER), in the estrogen-responsive cell line ERα CALUX. We found that this widely used cell line overexpresses AOXs and NMTs while having low expression of CYP enzymes. Exposure of ERα CALUX cells to quinolines resulted in estrogenic effects, which could be mitigated using an inhibitor of AOX/NMTs. No such mitigation occurred after coexposure to a CYP1A inhibitor. A number of N-methylated but no hydroxylated transformation products were detected using liquid chromatography-mass spectrometry, which indicated that biotransformations to estrogenic metabolites were likely catalyzed by NMTs. Compared to the natural ER ligand 17ß-estradiol, the products formed during the metabolization of quinolines were weak to moderate agonists of the human ERα. Our findings have potential implications for the risk assessment of these compounds and indicate that care must be taken when using in vitro estrogenicity assays, for example, ERα CALUX, for the characterization of N-heterocycles or environmental samples that may contain them.

A rapid and sensitive fluorometric method for determination of aldehyde oxidase activity.

Previous research has characterized the important role of aldehyde oxidases (AOX) in biotransformation of N-heterocyclic therapeutic drugs and environmental contaminants in mammals. Research pertaining to AOX activity in non-mammalian vertebrates, however, is scarce, despite its biological role as a potentially important metabolic pathway for xenobiotics. One of the limiting factors of research on AOX is that available photometric methods are relatively insensitive, limited in throughput, and prone to cross-reactivity from other enzymes. Therefore, this study aimed to develop a novel and improved fluorometric AOX assay. This assay is based on the conversion of the exogenous aldehyde substrate 4-(dimethyl)amino cinnamaldehyde to its corresponding fluorescent acid by AOX, and was evaluated using partially purified hepatic cytosol from rat, human, and rainbow trout. Purification of native cytosol by heat treatment and ammonium sulfate precipitation resulted in increased specific activity of AOX. Michaelis-Menten kinetic parameters (Kmand Vmax) were comparable to values previously generated by photometric methods. Furthermore, effects of the inhibitor hydralazine on AOX activity revealed half maximal inhibitory concentrations comparable to those generated using conventional methods. Product identity was confirmed by liquid chromatography and mass spectrometry. In summary, this study successfully developed a rapid and sensitive assay for determination of AOX activity in across different vertebrate species that is 4- to 10-fold more sensitive compared to conventional absorbance-based methods. It can be applied in environmental, toxicological, and pharmacological studies relating to identification of AOX substrates, as well as the induction of AOX expression through drugs and environmental contaminants.

Toward Streamlined Identification of Dioxin-like Compounds in Environmental Samples through Integration of Suspension Bioassay.

Effect-directed analysis (EDA) is a powerful strategy to identify biologically active compounds in environmental samples. However, in current EDA studies, fractionation and handling procedures are laborious, consist of multiple evaporation steps, and thus bear the risk of contamination and decreased recoveries of the target compounds. The low resulting throughput has been one of the major bottlenecks of EDA. Here, we propose a high-throughput EDA (HT-EDA) work-flow combining reversed phase high-performance liquid chromatography fractionation of samples into 96-well microplates, followed by toxicity assessment in the micro-EROD bioassay with the wild-type rat hepatoma H4IIE cells, and chemical analysis of bioactive fractions. The approach was evaluated using single substances, binary mixtures, and extracts of sediment samples collected at the Three Gorges Reservoir, Yangtze River, China, as well as the rivers Rhine and Elbe, Germany. Selected bioactive fractions were analyzed by highly sensitive gas chromatography-atmospheric pressure laser ionization-time-of-flight-mass spectrometry. In addition, we optimized the work-flow by seeding previously adapted suspension-cultured H4IIE cells directly into the microplate used for fractionation, which makes any transfers of fractionated samples unnecessary. The proposed HT-EDA work-flow simplifies the procedure for wider application in ecotoxicology and environmental routine programs.

Cross-Species Extrapolation of Uptake and Disposition of Neutral Organic Chemicals in Fish Using a Multispecies Physiologically-Based Toxicokinetic Model Framework.

The potential to bioconcentrate is generally considered to be an unwanted property of a substance. Consequently, chemical legislation, including the European REACH regulations, requires the chemical industry to provide bioconcentration data for chemicals that are produced or imported at volumes exceeding 100 tons per annum or if there is a concern that a substance is persistent, bioaccumulative, and toxic. For the filling of the existing data gap for chemicals produced or imported at levels that are below this stipulated volume, without the need for additional animal experiments, physiologically-based toxicokinetic (PBTK) models can be used to predict whole-body and tissue concentrations of neutral organic chemicals in fish. PBTK models have been developed for many different fish species with promising results. In this study, we developed PBTK models for zebrafish (Danio rerio) and roach (Rutilus rutilus) and combined them with existing models for rainbow trout (Onchorhynchus mykiss), lake trout (Salvelinus namaycush), and fathead minnow (Pimephales promelas). The resulting multispecies model framework allows for cross-species extrapolation of the bioaccumulative potential of neutral organic compounds. Predictions were compared with experimental data and were accurate for most substances. Our model can be used for probabilistic risk assessment of chemical bioaccumulation, with particular emphasis on cross-species evaluations.

Effect-Directed Analysis of Aryl Hydrocarbon Receptor Agonists in Sediments from the Three Gorges Reservoir, China.

The construction of the Three Gorges Dam (TGD) in the Yangtze River raises great concern in ecotoxicological research since large amounts of pollutants enter the Three Gorges Reservoir (TGR) water bodies after TGD impoundment. In this work, effect-directed analysis (EDA), combining effect assessment, fractionation procedure, and target and nontarget analyses, was used to characterize aryl hydrocarbon receptor (AhR) agonists in sediments of the TGR. Priority polycyclic aromatic hydrocarbons (PAHs) containing four to five aromatic rings were found to contribute significantly to the overall observed effects in the area of Chongqing. The relatively high potency fractions in the Kaixian area were characterized by PAHs and methylated derivatives thereof and heterocyclic polycyclic aromatic compounds (PACs) such as dinaphthofurans. Benzothiazole and derivatives were identified as possible AhR agonists in the Kaixian area based on nontarget liquid chromatography-high resolution mass spectrometry (LC-HRMS). To our knowledge, this study is the first one applying the EDA approach and identifying potential AhR agonists in TGR.

Understanding receptor-mediated effects in rainbow trout: in vitro-in vivo extrapolation using physiologically based toxicokinetic models.

The European REACH regulation requires the use of animal experimentation to assess the risk of industrial chemicals. However, the 3R principle (reduction, replacement, refinement) demands the use of suitable alternative test methods. Many dossiers submitted for the authorization of chemicals have attempted to provide the required data without performing new experiments, relying heavily on in silico methods; in vitro assays were scarcely used. We propose a methodology that uses physiologically based toxicokinetic (PBTK) models to extrapolate in vitro data to the in vivo level. We collected experimental results for in vitro and in vivo ethoxyresorufin-O-deethylase and vitellogenin induction following chemical exposure and compared those results with model predictions. We found that the predictive power of aqueous chemical concentrations was limited; median effect concentrations (EC50s) based on internal concentrations in fish correlated better with in vitro EC50s. Our data show that in vitro assays could offer a substitute for fish studies when combined with PBTK models.

Heterocyclic aromatic hydrocarbons show estrogenic activity upon metabolization in a recombinant transactivation assay.

Heterocyclic aromatic hydrocarbons (hetero-PAHs) are increasingly studied at contaminated sites; especially at former industrial facilities where coal tar-oil was handled, e.g., wood treatment plants, high concentrations of hetero-PAHs are frequently detected in groundwater plumes. In previous studies, fractions of groundwater with high estrogenic activity contained hetero-PAHs and their hydroxylated metabolites. To evaluate this preliminary evidence, selected hetero-PAHs were screened for their estrogenic activity in lyticase yeast estrogen screen (LYES) and ER CALUX. All tested substances were inactive in the LYES. Hetero-PAHs such as acridine, xanthene, indole, 2-methylbenzofuran, 2,3-dimethylbenzofuran, dibenzofuran, dibenzothiophene, quinoline, and 6-methylquinoline were positive in the ER CALUX, with estradiol equivalence factors (EEFs) from 2.85 × 10(-7) to 3.18 × 10(-5). The EEF values of these substances were comparable to those of other xenoestrogens (e.g., alkylphenols or bisphenol A) that are sometimes found in surface water. Chemical analyses revealed that T47Dluc cells could metabolize most of the substances. Among the metabolites (tentatively) identified by liquid chromatography-tandem mass spectrometry (LC-MS/MS) were hydroxides and their keto tautomers, sulfates, sulfoxides, and N-oxides. Because of their high concentrations measured in groundwater, we conclude that hetero-PAHs and metabolites may be a potential risk and should be the subject of further research.