Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Endocrine disrupting potency of organic pollutant mixtures isolated from commercial fish oil evaluated in yeast-based bioassays

  • Marek Łukasz Roszko ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Writing – original draft, Writing – review & editing

    roszko@ibprs.pl

    Affiliation Department of Food Analysis, Institute of Agricultural and Food Biotechnology, Rakowiecka, Warsaw, Poland

  • Marta Kamińska,

    Roles Investigation

    Affiliation Department of Food Analysis, Institute of Agricultural and Food Biotechnology, Rakowiecka, Warsaw, Poland

  • Krystyna Szymczyk,

    Roles Formal analysis

    Affiliation Department of Food Analysis, Institute of Agricultural and Food Biotechnology, Rakowiecka, Warsaw, Poland

  • Katarzyna Piasecka-Jóźwiak,

    Roles Writing – review & editing

    Affiliation Department of Fermentation Technology, Institute of Agricultural and Food Biotechnology, Rakowiecka, Warsaw, Poland

  • Beata Chabłowska

    Roles Investigation

    Affiliation Department of Fermentation Technology, Institute of Agricultural and Food Biotechnology, Rakowiecka, Warsaw, Poland

Abstract

The aim of this work was to evaluate the activity of xenobiotic mixtures containing persistent organic pollutants isolated from commercial fish oil samples against sex hormone receptors, including estrogen and androgen. The applied bioassay was based on transgenic yeast strains. The mixtures were extracted from the samples using the semi-permeable membrane dialysis technique and analyzed with gas chromatography/ion trap mass spectrometry. It turned out that mixtures of chemicals isolated from fish oil may interact with human steroid sex hormone receptors in various ways: the tested samples showed both estrogenic and anti-androgenic activity. Calculated 17β-estradiol equivalents for the tested samples ranged between 0.003 and 0.073 pg g–1 (fat). Anti-androgenic activity expressed as the flutamide equivalent concentration was in the 18.58–216.21 ng g–1 (fat) range. Polychlorinated biphenyls and various DDT metabolites were the main fish oil pollutants influencing the receptors. Additivity and/or synergy between chemicals was observed in the ER/AR mediated response.

Introduction

Endocrine disrupting chemicals (EDCs) are compounds that may interfere with human or wildlife endocrine systems and produce adverse developmental, reproductive, neurological and immune effects. A vast number of chemical compounds are known to mimic native human hormones like androgen (AR) and estrogen (ER) and to interact with their nuclear receptors (NRs) [13]. EDCs are able to interact with hormone NRs as agonists or antagonists [4]. They are able not only to activate the receptors in the absence of native hormones but also to increase or inhibit the response of endogenous ligands [5]. Therefore, environmental endocrine disrupters can be grouped into four categories: direct-acting agonists, direct-acting antagonists, indirect-acting agonists, and indirect-acting antagonists [67]. Competitive binding is a required condition for direct-acting compounds, and their affinity can be measured using ligand-binding assays. Some chemicals act indirectly, showing endocrine disrupting potency only after metabolic activation, i.e., after transformation into a related derivative but more native, structurally similar to hormone [89].

The presence of EDCs in the environment has raised concerns about the potentially harmful effects of exposure to such chemicals. The effects include modern civilization diseases such as cancer [1012], and/or wildlife reproductive cycle disturbances [1316]. The general public is chronically exposed to a wide group of EDCs from both natural and anthropogenic sources [17]. EDCs include polychlorinated biphenyls (PCBs), a number of pesticides, polycyclic aromatic hydrocarbons (after metabolic activation), brominated flame retardants, bisphenol A, pharmaceuticals, several classes of phytoestrogens and others; they are quite ubiquitous in the environment [1823]. Several routes of exposure have been identified. However, food is the main exposure path of humans for lipophilic and persistent organic pollutants (POPs) [24]. Most commonly, EDCs in the food/environment are complex mixtures of different compositions, and therefore it is difficult to evaluate their biological potency by means of chemical analyses only. If any conclusions on hormonal activity of a mixed sample are to be meaningfully drawn from data on the potency/effectiveness of individual compounds of the mixture, several conditions must be met: (i) the mechanism of biological activity must be determinable; (ii) the measured activity of the compounds must be correlated with the induced biological effects; and (iii) the effects must be correlated with the combined ED potency. Unfortunately, such assumptions may be not true when ER- and/or AR-interacting compounds are evaluated due to the omnipresence of various compounds in the environment with both agonistic and antagonistic ED activity. In addition, such an approach is truly applicable only if the dose-response curve is linear. However, sigmoidal curves and non-monotonic responses are frequently observed for estrogenic/androgenic agents [2, 2526]. Nevertheless, there is scientific evidence that chemical mixtures of substances sharing the same mode of action elicit predominately additive effects, both in vitro and in vivo [2731]. In addition, synergistic effects could also be possible [3233]. It should, however, be noticed that evaluation of interactions between chemicals is not an easy task, and wrong conclusions might be drawn due to a lack of information on the dose-response relationship [25,32].

The majority of the published studies on EDCs have focused on the ability of individual chemicals to interact with hormone receptors. Quite a few papers dealt with chemical mixtures derived from real-life food samples.

Fish oil is a food constituent. It is a rich source of omega-3 fatty acids and is commonly used as a dietary food supplement. Since fish oil might be contaminated with several classes of persistent chemicals, it has raised concerns on the potential health effects regarding its dietary use. These include potential endocrine disruption caused by the chemical mixtures present in these oils. As reported in our previous papers (e.g., [34]), fish liver oil may be heavily contaminated with chemicals of proven ED potency, including chlorinated biphenyls, brominated diphenyl ethers and its metabolites, or legacy organochlorine pesticides. In this work, estrogen/androgen agonistic and antagonistic activities of fish oil extracts were evaluated using optimized recombinant yeast assays. Since the combined ED potency of a mixture most likely depends on the NR binding potency of the individual mixture compounds, the isolated compounds were fractionated to identify the compounds mostly responsible for the observed activity.

Materials and methods

Chemicals/reagents

HPLC grade organic solvents were exclusively used in this study. N-hexane, cyclohexane, dichloromethane, and toluene were supplied by Sigma-Aldrich (Belefonte, PA, USA). Anhydrous sodium sulfate, potassium hydroxide, and sulfuric acid (96%) of p.a. grade were supplied by Avantor (Gliwice, Poland). Chlorophenol red galactopyranoside (CPRG), basic aluminum oxide, magnesium sulfate, copper sulfate, potassium phosphate, dihydrogen potassium phosphate, ammonium phosphate, iron (II) sulfate, L-threonine, L-leucine, L-histidine, adenine, L-arginine hydrochloride, L-arginine, L-methionine, L-tyrosine, L-isoleucine, L-lysine hydrochloride, L-phenylalanine, L-glutamic acid, L-valine, L-serine, thiamine, pyridoxine, pantothenic acid, inositol, biotin, aspartic acid and glucose were supplied by Sigma-Aldrich (Belefonte, PA, USA). Silica gel 60 (0.063–200 mm) was obtained from Merck (Darmstadt, Germany). YES/YAS yeast strains were obtained from Tigret (Warsaw, Poland). Pure PBDE standards IUPAC #1, 2, 3, 7, 8, 10, 11, 12, 13, 15, 17, 25, 28, 30, 32, 33, 35, 37, 47, 49, 66, 71, 75, 77, 85, 99, 100, 110, 116, 118, 119, 126, 138, 153, 154, 155, 166, 183, and 209 (>97%) were supplied by AccuStandard (New Haven, CT, USA). Aroclor PCB mixtures 1221, 1232, 1242, 1248, 1254, 1260, 1262, and 1268 were obtained from Sigma-Aldrich (Belefonte, PA, USA). o,p-DDT, p,p-DDE, p,p-DDD, p,p-DDT, heptachlor, o,p-DMDT, p,p-DMDT, HCB, trans-/cis-chlordan, and HCH were supplied by IPO (Warsaw, Poland). High purity (>97%) native PAH standards (naphthalene, acenaphtylene, acenaphtene, fluorene, phenanthren, anthracene, fluoranthrene, pyrene, benzo[c]fluorene, benzo[c]phenanthrene, chrysene, benzo[a]anthracene, cyclopenta[c,d]pyrene, benzo[b]fluoranthrene, benzo[j]fluoranthrene, benzo[k]fluoranthrene, 7,12-dimethylbenz[a]anthracene, benzo[a]pyrene, benzo[e]pyrene, 3-methylcholanthrene, indeno[1,2,3-cd]pyrene, dibenzo[a,h]anthracene, benzo[g,h,i]pyrelene, dibenzo[a,l]pyrene, dibenzo[a,i]pyrene, dibenzo[a,e]pyrene, dibenzo[a,h]pyrene, and 5-methylchryzene) were supplied by Dr. Ehrenstorfer (Augsburg, Germany) and Accustandard (New Haven, CT, USA). ER/AR native agonists 17β-estradiol and 17β-hydroxy-testosterone as well as antagonists flutamide and hydroxy-tamoxifen were supplied by Sigma-Aldrich (Belefonte, PA, USA). Streptomycin and ampicillin antibiotics were obtained from Sigma-Aldrich (Belefonte, PA, USA).

Sample preparation

Our previously published sample preparation method with some major modifications was used. The fraction containing organic pollutants was isolated from fish oil samples (4 × 5 g) using semi-permeable dialysis [35]. Dialysates were rotary-evaporated to dryness and de-fatted with gel permeation chromatography. Separations were performed on a 500×10 (Omnifit, Cambridge, UK) column filled with Bio-Beads SX-3 styrene/divinylobenzene resin. The dichloromethane:cyclohexane (1:1) mobile phase eluted at a rate of 1 ml min–1, and a 2 ml sample loop was used. Fractions eluting between 25 and 55 ml were collected. Fractions eluted from individual sub-samples were combined, evaporated to almost dryness and re-dissolved in n-hexane. The phenolic fraction (containing hydroxylated polychlorinated biphenyls/brominated diphenyl ethers) was isolated through alkaline partitioning as described previously [36]. The neutral fraction was fractionated using a 3 g basic alumina column prepared as described previously [37]. The column was eluted with 35 ml of n-hexane, 35 ml of a dichloromethane:n-hexane mixture (1:99), and dichloromethane. The n-hexane fraction containing mainly aliphatic compounds was additionally passed through a glass column containing 1 g of acidic silica gel. Samples were evaporated to dryness and re-dissolved in 500 μl of ethyl acetate. To eliminate possible variations of fraction composition resulting from differences in recovery rates of individual compounds, one half volumes of each fraction were combined into a composite sample, while the second half volumes were analyzed individually. Then, 50 μl of each extract was transferred into a chromatographic vial for GCMS analysis, while the remaining portion was transferred into flat bottom glass chromatographic vials, evaporated to dryness within a gentle stream of nitrogen, and submitted for yeast-based bioassays. Statistical parameters of the applied analytical procedures were discussed in detail previously [3537]. Recovery rates of the tested legacy pesticides from the spiked samples are shown in S1 Table.

GC/MS analyses

Samples were analyzed in a TraceGC Ultra Thermo-Finnigan (Austin, TX, USA) gas chromatograph coupled to a PolarisQ Thermo-Finnigan (Austin, TX, USA) ion trap mass spectrometer equipped with a TriPlus autosampler (Austin, TX, USA) and a PTV injector. Chromatographic separations were performed using Rtx–17 MS medium-polarity fused-silica capillary column (Restek, Bellefonte, PA, USA) and a 5 m × 0.53 mm guard column/retention gap. The following temperature program was used during separations: 40°C 3 min hold, ramp to 200°C (25°C min–1), ramp to 250°C (1.4°C min–1), 2 min hold, ramp to 300°C (25°C min–1), 5 min hold. Then, 40 μl of the extract was introduced into the GC instrument using the solvent split injection method. Initial injector temperature was 40°C and the injection rate was 3 μl s–1. Helium split flow at the injection was set at 200 ml min–1. The split valve was closed after 1 minute and the temperature ramped to 300°C (25°C s–1) and held for 2 min. Subsequently, the split valve was opened, the flow rate set to 100 ml min–1, the temperature ramped to 320°C and held for 5 min. The mass spectrometer ion source was operated at 300°C and the GC transfer line temperature was set at 320°C. Other mass spectrometer operational parameters were as reported previously [36].

Yeast estrogen and androgen bioassays

Yeast-based bioassays were performed according to Schulitz and Metzger [38] with some modifications. Minimal medium was prepared using 13.61 g KH2PO4, 1.98 g (NH4)2SO4, 4.2 g KOH, 0.2 g MgSO4, 0.8 mg Fe2(SO4)3, 50 mg L-leucine, 50 mg L-histidine, 50 mg adenine, 20 mg L-arginine, 20 mg L-methionine, 30 mg L-tyrosine, 30 mg L-isoleucine, 30 mg L-lysine, 25 mg L-phenylalanine, 100 mg L-glutaminic acid, 150 mg L-valine, and 375 mg L-serine dissolved in 1000 ml of deionized water. The solution was sterilized at 121°C for 15 min. A vitamin solution was prepared using 8 mg thiamine, 8 mg pyridoxine, 8 mg pantothenic acid, 40 mg inositol, and 0.4 mg biotin dissolved in 200 ml of deionized water and filtered through a sterile 0.2 mm nylon filter (Whatman PURADISC). The glucose solution (20%), aspartic acid solution (4 mg ml–1), threonine solution (24 mg ml–1) and 20 mM copper sulfate were sterilized at 121°C for 15 min. Combined growth medium was prepared by mixing solutions of glucose (5 ml), aspartic acid (1.25 ml), vitamins (1 ml), threonine (0.5 ml), copper sulfate (0.125 ml) and minimal medium (45 ml). Growth medium was spiked before use with ampicillin and streptomycin to give a final concentration of 1 mg ml–1 each.

Briefly, an overnight yeast (YES/YAS) culture with approx. 1.2 OD@690 nm was centrifuged at 1000 × g and re-suspended in 10 ml of fresh medium. Samples mixed with 120 μl of 4% ethanol solution in water were placed in 2 ml glass vials, inoculated with 60 μl of yeast suspension, shaken by hand, transferred into a tightly closed plastic box with the bottom covered with water (to assure sufficient humidity), and incubated for 8 hours at 36°C. Next, the samples were transferred into a 96-well plastic plate, and their optical density was measured. The culture was terminated with the addition of 20 μl of 0.1% Triton X–100 and 30 μl of lyticase (1 mg ml–1) and mercaptoethanol (50 mM) in 0.1 M phosphate buffer (pH 7.5) solution. Samples were incubated for 45 min at 25°C and afterwards mixed with CPRG solution (1.0 mg ml–1). Evolution of absorbance at 570 and 690 nm was measured six times at 30-minute intervals (total 3 hours).

Test material

Fish (cod and shark liver) oils used as dietary supplements available through retail trade were used for the purposes of this study. Ten different oil brands were used. Pure standards of chemicals found at abundant concentrations in the oil samples including DDE, DDT, PCB mixtures, and several PBDE congeners were tested for estrogenic and androgenic activity. Chemicals were evaluated standalone, in combination and in the presence of a fixed concentration of native hormones (to evaluate antagonistic and combined effects).

Data analysis

Xcalibur 1.2 (GCMS) and Galapagos 2.0 (Plate reader) software packages were used to acquire and analyze the data. The results are given as the means of at least two parallel determinations (±SD, if given). Authentic standards and mass spectra were used to identify chromatographic peaks. Compounds for which standards were not available were identified using mass spectra and quantified on the basis of the nearest available structural analogues. Calibration curves were plotted as peak area vs. analyte concentration.

If some cytotoxic effects inhibited yeast growth, the results of the respective bio-assay were excluded from the calculations. The samples were assessed in terms of agonist or antagonist activity against estrogen or androgen receptors. The results were expressed as mass equivalents of 17β-estradiol, hydroxy-tamoxifen (estrogen agonist), hydroxy-testosterone and flutamide (androgen agonist). Equivalents were calculated from EC 50 values estimated from dose–response curves plotted for samples and native hormones. In cases where EC 50 values were not achieved, EC 20 values were used for calculations. Estimations of response and/or activity for the tested chemicals mixtures were performed using two models: A–simple additive model where Emix is the estimated effect of a mixture of chemicals, and Ei is the effect of individual chemicals at a defined concentration; B–predictive model according to Payne et al. [27], where Emax is the maximal effect (absorbance readings), and the Fi effect is predicted from the dose response curve of the individual chemical at a defined concentration. For the purpose of this study, experimental and measured values were used for calculations instead of using values estimated from the asymmetric Hill functions. Spearman’s rank test was used to test the correlation between the evaluated chemicals concentration and the hormone equivalents. One way ANNOVA was used to evaluate the significance of the differences between measured absorbance values and values estimated from the predictive models.

Results and discussion

Fish oil chemical contamination

POPs identified in the tested fish oil samples are listed in Table 1. Polychlorinated biphenyls (PCB) were found at the highest concentrations, consistent with our previous reports and the reports of other authors [34, 3940]. Relatively high levels of various DDT metabolites and degradation products (including DDMU, DDD, and DDE) were also found, consistent with the reports of Rawn et al. [41], who reported high DDT concentrations in some types of fish oils. The Principal Component Analysis technique was used to compare profiles of PCB congeners found in the tested samples with the composition of some technical mixtures containing those compounds. The analysis showed that the found profiles were most congruent with profiles characteristic of the Aroclor 1260 technical mixture (Fig 1). However, since a number of both abiotic and biological factors affect the profiles, that does not necessarily mean that such a technical mixture was the real contamination source [3435, 37, 42]. Nevertheless, to facilitate interpretation of the obtained results, PCBs concentrations were calculated as Aroclor 1260 equivalents using the authentic Aroclor technical mixture as the standard. Calculations were based on the summary chromatographic peak area of the most abundant congeners.

thumbnail
Fig 1. Principal component analysis scatter plot of PCB profiles found in the tested oil samples and in several PCB-based technical mixtures.

Samples have been indicated with numbers.

https://doi.org/10.1371/journal.pone.0197907.g001

thumbnail
Table 1. Concentrations of organic pollutants identified in the tested oil samples.

https://doi.org/10.1371/journal.pone.0197907.t001

Endocrine disrupting potency of POPs isolated from commercial fish oil

Bioassay results are given in Table 2. Calculated estradiol equivalents for all fractions eluted from the combined samples ranged from values below 0.003 pg g–1 up to 0.073 pg g–1 (fat). None of the tested samples acted as the ER antagonist. Conversely, they produced an additive response and/or show synergy and acted as ER agonists of various activity and potency in the applied assay. The alkaline fraction, i.e., the fraction potentially containing phenolic compounds [36], including metabolites of PCB/PBDE (OH-PCB/PBDE), did not show any hormonal activity in the assay and was not submitted for further evaluations. This was the most likely due to very low concentrations of these compounds in fish oil [34]. Compounds present in the dichloromethane and dichloromethane:n-hexane fractions obtained using basic alumina chromatography were responsible for the observed estrogenic potential. The hexane fraction contained octachlorostyrene, unidentified chlorinated compounds, and heptachlor. Most of the chlorinated biphenyls, DDE, trans-nonachlor, HCB and pentachlorobenzene were found in the hexane:dichloromethane fraction. DDTs, methoxychlor (DMDT), nonachlor, PBDEs and HCHs were found in the dichloromethane fraction. Relative contributions of the chemicals identified in the above fractions are shown in S1 Fig. Since the chemicals eluting in various fractions differed significantly in their spatial structure and polarity, this determined its receptor affinity and matching. In cases of substitution of isomers (like CBs congeners) slight changes in structure may significantly modify the affinity or activity of the compound [43].

thumbnail
Table 2. Calculated equivalent concentrations of (anti-)hormones for the tested oil samples.

LOQ was taken as the concentration of the reference compound causing less than a 10% change in response (<EC10) assuming a 20 g sample mass.

https://doi.org/10.1371/journal.pone.0197907.t002

Hormone equivalents of the combined fractions were higher than the calculated simple sums of the equivalents of the separated fractions. This is not unusual because in most cases, hormone-dependent response curves are exponential rather than linear [26, 4445]. The fractions were combined after chromatography to eliminate possible differences in recovery rates. Interactions between individual ligands may not be excluded at this point. Chemicals present in the fractions may also act in synergy, multiplying the observed response [4647]. Competitiveness of ligands against ER-active sites is also possible. In addition, since the tested mixtures were rather complex, anti-hormonal activity of individual compounds may not be excluded even if such activity of the entire mixture was not indeed observed. In some cases, a subtle and specific disruption in the estrogenic signaling system of the recombinant yeast test might also occur at concentrations that would not result in apparent growth inhibition. This phenomenon has been described as toxic masking [32]. Pliskova et al. [48] reported anti-estrogenic activity of higher-chlorinated CBs while the lower chlorinated ones were weakly estrogenic. Contamination with higher-chlorinated congeners was characteristic of samples tested in this study. Lack of overall anti-estrogenic activity of the CB-containing fraction might very well be caused by some DDT metabolites or other compounds exhibiting weak estrogenic activity, acting in synergy or simply producing an additive response. Anti-estrogenic activity of Aroclor technical mixtures at the tested concentrations against E2 was not proven in this study. In contrast, an additivity or synergistic effect with E2 was observed in all Aroclor mixtures tested, see S2 Fig.

To evaluate the possible interactions between the most abundant chemicals (DDT, DDE, Aroclor 1260) present in oils, some additional experiments were performed using pure standards of those chemicals and its mixtures. These results are shown in Fig 2. The experiments were performed at concentrations similar to those observed in the test samples (taking into account extract pre-concentration). The measured absorbance values for the chemical mixtures have been compared with the effects estimated from the predictive models. As noted previously, DDE, DDT and Aroclor 1260 were found to show weak estrogenicity with DDT found to be the most potent. Higher values were measured for combinations of Aroclor 1260 and E2 mixtures than that predicted from the models. Such observations were also noticed for the tested higher DDE concentrations. The results obtained with the predictive models for DDT were generally in line with the measured values, suggesting a pure additive response (at the tested concentration range). Surprisingly, the applied predicative models have failed when combined mixtures of tested chemicals were evaluated.

thumbnail
Fig 2. Measured and estimated relative response in the yeast estrogen assay for the tested chemicals.

Aroclor 1260, DDT, and DDE concentrations in mixtures were respectively 3.2*10−4, 2.61*10−5, 2.6*10−4; and 4.0*10−6, 6.68*10−7, and 3.2*10−6 (n = 3). Statistically homologous groups have been indicated with letters (significance of the differences was tested between measured values and predictions).

https://doi.org/10.1371/journal.pone.0197907.g002

Hardrup et al. [49] noticed that the prediction of biological effects of hormonally active mixtures containing both chemicals with agonistic as well as chemicals having antagonistic potency with different mathematical models have significant shortcomings and do not always provide sufficient estimation. Brinkmann et al. [50] has also noted that additive models are unable to represent the effects of mixtures of partial agonists due to lesser efficacies compared to native hormones. The maximum achievable effect level and slopes of the dose-response curves are not equal for different compounds. This was also the case for the chemicals tested in our study. In addition, competitive antagonism is also possible. When the most suited model was applied, the predicted estrogenicity for the environmental samples reported by Brinkmann et al. [50] differed in the 5–622% range from the measured value. The results were both over and underestimated depending on the sample evaluated.

In addition, the course of action of the chemical mixtures tested in this study was found to be concentration dependent. At the higher (mixture) concentrations used, the measured combined effect was lower than the predicted one from the effects of individual chemicals at the tested concentrations. When E2 was added to the test mixture a higher response was observed, however not proven statistically different in case of higher mixture concentrations tested. Such observations might indicate that some of the tested chemicals showed anti-estrogenic activity in the presence of weak estrogens. This phenomenon might be explained by the higher affinity of some individual ligands than others to the receptor but not being able to activate the subsequent signaling pathway. Increased response in the presence of E2 might be explained by the extremely high affinity of native hormones to the receptor molecule. However, at lower concentrations, the same tested chemicals mixtures produced significantly higher responses that predicted from the models. This suggests synergy or at least more than additivity between the chemicals. This was also noticed when E2 was additionally added to the tested mixture. The results obtained for tested chemicals mixtures showed significant variability. The mean absorbance measured for the chemical mixture tested at a higher concentration was only proven significantly different from a mixture tested with E2 addition at 6.54*10−11 M. As noted previously, this might be attributed to the differences in the maximum achievable effect of individual chemicals and its receptor binding potency. Such observed alterations might be those most likely related to the chlorinated biphenyls present in the tested samples. As noted previously, they have been reported in the literature to show both estrogenic and anti-estrogenic activity [48]. According to the results of this study, such interactions are significantly concentration dependent.

On the basis of the calculated hormonal activities and determined concentrations of chemicals it may be concluded that the compounds mainly responsible for the observed disrupting potential of the isolated extracts include various DDT metabolites, PCBs, and PBDEs. o,p-DDT is the most potent estrogen among the latter compounds with 6.53 × 10−6 relative estrogen potency (REP) calculated from the EC50 values. Other tested compounds showed REP values within the 1 × 10−7–1 × 10−9 range (S3S5 Figs).

Correlation between hormone equivalents and chemical concentrations

Spearman’s rank test was used to statistically evaluate the results, see Table 3. Individual fractions were taken as cases, while hormone or anti-hormone activity and identified chemicals were taken as variables.

thumbnail
Table 3. Correlation coefficients between calculated concentration equivalents of E2 and flutamide and chemicals identified in the tested samples and fractions (Spearman rank test, α = 0.05).

Statistically significant correlations are bolded.

https://doi.org/10.1371/journal.pone.0197907.t003

It is important to note that CBs were substantively separated in the isolated fractions from PBDE/DDT metabolites except DDE. Some statistically significant monotonic correlations between concentrations and hormone equivalents are seen, including numerous correlations between the E2 equivalent and some CB/BDE/DDT metabolites. As noted earlier, that is not surprising since the latter metabolites were proven to be weak estrogens in the yeast bioassay used in this study. BDE # 47, the most conspicuous congener among all found PBDEs (see Fig 3) turned out to also be the most potent estrogen (S5 and S6 Figs). Its relatively highest share in fish oil PBDE fractions was confirmed in the literature [5152].

thumbnail
Fig 3. Average contribution of individual congeners to the PBDE profile in the tested fish oil samples.

https://doi.org/10.1371/journal.pone.0197907.g003

Other pollutants found in the tested samples include legacy organochlorine insecticides and unidentified chlorine substituted compounds. Concentrations of some of them (e.g., cis-/trans-nonachlor) correlated statistically significantly with the calculated E2 equivalent. Hormonal activity of these pollutants were not evaluated in this study; Klotz et al. [53] reported their weak activity. We noted some HCH/E2 correlations. Additionally, Garg et al. [54] reported HCH estrogenic activity, but Kim et al. [55] reported none.

Statistical tests indicated some chlordane/E2 correlations, but Okoumassoun et al. [56] did not report any. It must be stressed that in respect to the latter compounds the indicated results and correlations may be false. Numerous POPs are frequently correlated with each other in many environmental samples. Those could not be easily isolated or fractionated on a semi-preparative scale (as applied in this study). Since DDT/PCB and HCH/chlordane coexist in the extracts at very different concentrations, it is not possible to undoubtedly evaluate their individual activity/contribution to the observed hormonal potency of the entire mixture. However, as mentioned previously, it was also assumed that no other (unidentified) compound potentially present in the tested samples was hormonally active.

The calculated estrogen equivalent was significantly correlated with the concentration of Aroclor 1260. The estrogen potency of the latter is relatively low but it is an E2 agonist (S2 Fig). This might indicate a possible additivity or synergy between chlorinated biphenyls and other weak estrogens present in food as already demonstrated previously [20,25].

All evaluated oil samples showed anti-androgenic activity in the applied assay. Estimated flutamide equivalents were in the 18.58–216.21 ng g–1 (fat) range. Analysis of fractionated samples indicated a significant correlation between the observed endocrine disrupting potency expressed as the flutamide concentration and the concentration of DDT metabolites and chlorinated biphenyls. In the majority of the tested samples the highest response was observed for dichloromethane and dichloromethane:hexane basic alumina fractions. The n-hexane fractions (composed of mostly aliphatic compounds) showed anti-androgenic activity in two cases (see Table 2, samples 1 and 4). However, compounds responsible for the activity were not identified (different analytical and fractionation methods would be required).

Both DDE and chlorinated biphenyls were found to show anti-androgenic activity using pure standards of those compounds (S7 and S8 Figs). These compounds have also been reported as androgen agonists [5758]. PBDE concentration did not correlate with flutamide equivalent, most likely due to its relatively low concentration compared to CBs/DDE and its low ED potency.

Estrogens and androgens in food

Equivalent concentrations of hormones (E2) estimated in the yeast-based bioassay are low compared to concentrations of endogenous hormones (estradiol, estrone, estriol) naturally occurring in food of animal origin. For example, cow’s milk may contain up to several hundred pg g–1 of endogenous estrogens, including those conjugated with proteins [5960]. These levels are significantly higher than the equivalents calculated in this study. Animal origin food may also be a source of endogenous testosterone [59,61]. Our results are in line with Behr et al., [62] who reported low but detectable estrogenicity of farmed fish meat. The authors suggested possible residues of genistein, i.e., phytoestrogen derived from soy used as aquaculture feed. Estrogenicity of farmed fish muscle demonstrated by Pinto et al. [63] using a YES-based bioassay originated most likely (as suggested by the authors) from the presence of some POPs. Similar to the results of our study, a positive correlation between total PCB levels and the observed estrogen potency was reported by the authors.

There is not much available literature data on anti-androgens equivalents in oil samples or food. The calculated equivalent concentrations seem to be low in the tested samples with the following reservations. The studied activity referred to only one specific element of the biochemical response and only one isoform of the receptor was evaluated. In addition, every compound found in the tested samples may undergo metabolic transformation or activation in humans (including oxidation with various cytochromes), forming some more potent EDCs capable of more seriously affecting human health This has been, for example, well documented for hydroxylated metabolites of chlorinated biphenyls or hydroxylated metabolites of polycyclic aromatic hydrocarbons [1819, 2122, 6467].

Health issues

On the basis of the available literature it is difficult to unambiguously evaluate the health effects of chronic exposure to the chemicals isolated and tested in this study. Scientific communities have not reached a consensus on the potential health effects of low-dose exposure to weak EDCs [26] and consequently on the safe/tolerable exposure/intake levels. The Swedish Chemical Agency [68] identified several obstacles to this, including the complexity of the human endocrine system, variable susceptibility at various human development stages, and long time periods between exposure and occurrence of any observable adverse effects. Both the Swedish Chemical Agency, [68] and Hass et al. [44] pointed out that it was impossible to determine a threshold for a functioning system that under normal conditions was already active. In such a scenario, even the slightest change in hormone levels (hypothetically few molecules) would result in physiological effects. Endocrine disruptor activity estimations are burdened with high uncertainty, especially when experimental results are extrapolated to different species.

Dybing et al. [69] indicated that some level of interaction with or occupation of receptor active sites with endocrine disruptors must be achieved to produce toxic effects in view of the numerous mechanisms of cyto- and homoeostasis protection and the abundance of cellular and molecular elements able to interact with hormonal signaling. Body homeostasis mechanisms are able to counteract any disrupting action produced by xenobiotics at levels below the critical thresholds. However, if the mechanisms are less operational at some stage of organism development, susceptibility to xenobiotics appearing at that stage might be higher.

Regardless of the above, the results of estrogen equivalents calculated for the tested samples may be compared with the acceptable daily intake (ADI) values reported in the literature for the reference compounds. Such an approach has been previously adopted by Caldwell et al. [70]. The WHO [71] proposed an ADI of 0.05 μg E2 kg-1 b.w. day-1 based on a no observed effect level (NOEL). For the purpose of comparison of fish oil related exposure to the toxicity-based benchmark, a conservative scenario was adopted including the following assumptions. An average adult body weight of 60 kg, a maximum concentration (0.073 pg g-1) was used for calculations, and 80 g dietary fat intake (solely from fish oil). Based on such assumptions the average calculated E2 equivalent would amount to 5.84 pg E2 equivalent person-1 day-1. This number is much lower than the 3 μg E2 equivalent person-1 day-1 ADI value. Using the same assumptions, the intake of flutamide equivalent was estimated at 17.28 μg person-1 day-1. The ADI value for flutamide not causing anti-androgenic effects was set at 2.5 μg kg-1 b.w. day-1 which results in 150 μg kg-1 person-1 day-1 [72]. The estimated values in the worst-case scenario are only 10 times lower than the ADI. Despite the fact that the adopted worst-case scenario is unlikely, the calculated intake value might suggest that extensive use of contaminated fish oil in some (sensitive) population groups might be a health concern if anti-androgenic effects are evaluated. In addition, as noted by De Falca et al. [73], a balance between androgens and estrogens may be important in maintaining normal spermatogenesis. Since the tested oil samples showed both estrogenic and anti-androgenic activity, its biological activity might be significantly stronger due to its complex mechanism of action. Exposure to xenoestrogens and anti-androgens during fetal and neonatal development has already been associated with a series of male reproductive disturbances, such as cryptorchidism, hypospadias, impaired fertility (especially due to poor semen quality), and an increased incidence of testicular cancer [73].

As noted previously, limitations of the model and tests applied in this study should be strongly considered. In addition, the mode of action of the chemicals evaluated is only one element in its toxicological profile, which includes carcinogenicity, genotoxicity or interactions with other NRs (including aryl hydrocarbon). This indicates that such contaminated fish oils might produce undesirable health effects, including sex hormone system disruption, despite providing valuable nutrients if used as a food.

Conclusions

The results of this study confirm that residues of POPs present in fish oil may disrupt the human endocrine system by means of agonistic and antagonistic effects on human estrogen and androgen receptors. Equivalent concentrations of the native hormone (E2) estimated using a yeast-based bioassay were low as compared to concentrations of endogenous hormones (estradiol, estrone, estriol) naturally occurring in food of animal origin. However, every tested sample showed some measurable anti-androgenic activity. Residues of polychlorinated biphenyls, DDT, and its metabolites contributed the most to the endocrine disrupting potency in comparison with other POPs present in our fish oil samples.

Supporting information

S1 Fig. Relative contribution of chemicals identified in the basic alumina chromatography fractions.

https://doi.org/10.1371/journal.pone.0197907.s001

(DOCX)

S2 Fig. Relative response in yeast estrogen assay for the tested Aroclor mixtures in the presence of E2 at 2 x 10−10 M (n = 3).

Control–first bar blank, second bar blank + E2. Maximum response from the E2 dose-response curve was used for calculations as a reference maximum effect (Amax).

https://doi.org/10.1371/journal.pone.0197907.s002

(DOCX)

S3 Fig. Dose-response curves for the tested Aroclor mixtures against ER (n = 3).

https://doi.org/10.1371/journal.pone.0197907.s003

(DOCX)

S4 Fig. Dose-response curves for the tested polybrominated biphenyls ethers against ER (n = 3).

https://doi.org/10.1371/journal.pone.0197907.s004

(DOCX)

S5 Fig. Relative response in yeast estrogen assay for the tested PBDEs in the presence of E2 at 2 x 10−10 M (n = 3).

Control–first bar blank, second bar blank + E2.

https://doi.org/10.1371/journal.pone.0197907.s005

(DOCX)

S6 Fig. Dose-response curves for the o,p-DDT and p,p-DDE.

OH-tamoxifen (HT) was used as an anti-estrogen positive control (n = 3).

https://doi.org/10.1371/journal.pone.0197907.s006

(DOCX)

S7 Fig. Antagonist response curves for DDE and Aroclor 1232 against AR.

Flutamide was used as anti-androgen positive control (DHT concentration at 2,78 x 10−9 M). DHT was plotted as agonist positive standard (n = 3).

https://doi.org/10.1371/journal.pone.0197907.s007

(DOCX)

S8 Fig. Antagonist response curves for BDE #47, 99, 100 against AR.

Flutamide was used as anti-androgen positive control (DHT concentration at 2,78 x 10−9 M) (n = 3). DHT was plotted as agonist positive standard.

https://doi.org/10.1371/journal.pone.0197907.s008

(DOCX)

S1 Table. Limits of quantification, limits of determination, recovery rates and recovery relative standard deviation for the evaluated organochlorine insecticides and related compounds.

N/D–not determined.

https://doi.org/10.1371/journal.pone.0197907.s009

(DOCX)

References

  1. 1. Blair RM, Fang H, Branham WS, Hass BS, Dial SL, Moland CL, et al. Estrogenic activity of styrene oligomers after metabolic activation by rat liver microsomes. Toxicol Sci. 2000; 54: 138–153.
  2. 2. Rajapakse N, Ong D, Kortenkamp A. Defining the Impact of Weak Estrogenic Chemicals on the Action of Steroidal Estrogens. Toxicol Sci. 2001; 60: 296–304. pmid:11248142
  3. 3. Bonefeld-Jørgensen EC, Long M, Hofmeister MV, Vinggaard AM. Endocrine-disrupting potential of bisphenol A, bisphenol A dimethacrylate, 4-n-nonylphenol, and 4-n-octylphenol in vitro: new data and a brief review. Environ Health Perspect. 2007; 115 (S1): 69–76.
  4. 4. Diamanti-Kandarakis E, Bourguignon JP, Giudice LC, Hauser R, Prins GS, Soto AM, et al. Endocrine-disrupting chemicals: an Endocrine Society scientific statement. Endocr Rev. 2009; 30 (4): 293–342. pmid:19502515
  5. 5. Marino M, Pellegrini M, La Rosa P, Acconcia F. Susceptibility of estrogen receptor rapid responses to xenoestrogens: Physiological Outcomes. Steroids. 2012; 77: 910–917. pmid:22410438
  6. 6. Cooper RL, Kavlock RJ. Endocrine disruptors and reproductive development: a weight-of-evidence overview. J Endocrinol. 1997; 152 (2): 159–166. pmid:9071972
  7. 7. Layton AC, Sanseverino J, Gregory BW, Easter JP, Sayler GS, Schultz TW. In vitro estrogen receptor binding of PCBs: measured activity and detection of hydroxylated metabolites in a recombinant yeast assay. Toxicol App Pharmacol. 2002; 180: 157–163.
  8. 8. Kitamura S, Sugihara K, Sanoh S, Fujimoto N, Ohta S. Metabolic activation of proestrogens in the environment by cytochrome P450 System. J Health Sci. 2008; 54: 343–355.
  9. 9. Conley JM, Hannas BR, Furr JR, Wilson VS, Gray LE. A demonstration of the uncertainty predicting the estrogenic activity of individual chemicals and mixtures from an in vitro estrogeni receptor transcriptional activation assay (T47D-KBluc) to the in vivo uterothropic assay using oral exposure. Toxicol Sci. 2016; 153:382–395.
  10. 10. Safe SH. Environmental and dietary estrogens and human health: is there a problem? Environ Health Perspect. 1995; 103: 346–351. pmid:7607134
  11. 11. Fucic A, Gamulin M, Ferencic Z, Katic J, Krayer von Krauss M, Bartonova A, et al. Environmental exposure to xenoestrogens and oestrogen related cancers: reproductive system, breast, lung, kidney, pancreas, and brain. Environ Health. 2012; 11 (S8): 1–9.
  12. 12. Hueiwang AJ. Exposure to endocrine disrupting chemicals and male reproductive health. Environ Health. 2014; 2 (55): 1–12.
  13. 13. Montagnani DB, Puddefoot J, Davie TJA, Vinson GP. Environmentally persistent oestrogen-like substances in UK river systems. J Water Environ. 1996; 10: 399–406.
  14. 14. Sonnenschein C, Soto AM. An updated review of environmental estrogen and androgen mimics and antagonists. J Steroid Biochem Mol Biol. 1998; 65: 143–150.
  15. 15. Tyler CR, Jobling S, Sumpter JP. Endocrine disruption in wildlife: A critical review of the evidence. Crit Rev Toxicol. 1998; 28: 319–361. pmid:9711432
  16. 16. Legler J, Zeinstra LM, Schuitemaker F, Lanser PH, Bogerd J, Brouwer A, et al. Comparison of in vivo and in vitro reporter gene assays for short-term screening of estrogenic activity. Environ Sci Technol. 2002; 36: 4410–4415. pmid:12387416
  17. 17. Campbell CG, Borglin SE, Green B, Grayson A, Wozei E, Stringfellow WT. Biologically directed environmental monitoring, fate, and transport of estrogenic endocrine disrupting compounds in water: A review. Chemosphere 2006; 65: 1265–1280. pmid:16979218
  18. 18. Connor K, Ramamoorthy K, Moore M, Mustain M, Chen I, Safe S, Z et al. Hydroxylated polychlorinated biphenyls (PCBs) as estrogens and antiestrogens: structure–activity relationship. Toxicol Appl Pharmacol. 1997; 145: 111–123. pmid:9221830
  19. 19. Elsby R, Ashby J, Sumpter JP, Brooks AN, Pennie WD, Maggs JL, et al. Obstacles to prediction of estrogenicity from chemical structure: assay mediated transformation and the apparent promiscuous nature of the estrogen receptor. Biochem Pharmacol. 2000; 60: 1519–1530. pmid:11020455
  20. 20. Elsby R, Maggs JL, Ashby J, Park BK. Comparison of the modulated effects of human and rat liver microsomal metabolism on the estrogenicity of bisphenol A: implications for extrapolation to humans. J Pharmacol Exp Ther. 2001; 297: 103–113.
  21. 21. Elsby R, Maggs JL, Ashby J, Paton D, Sumpter JP, Park BK. Assessment of the effect of metabolism on the estrogenic activity of xenoestrogens: a two stage approach coupling human liver microsomes and a yeast estrogenic assay. J Pharmacol Exp Ther. 2001; 296: 329–337. pmid:11160614
  22. 22. Fertuck KC, Matthews JB, Zacharewski TR. Hydroxylated benzo[a]pyrene metabolites are responsible for in vitro estrogen receptor–mediated gene expression induced by benzo[a]pyrene, but not elicit urotrophic effects in vivo. Toxicol Sci. 2001; 59: 231–240. pmid:11158716
  23. 23. Bergman A, Heindel JJ, Kasten T, Kidd KA, Jobling S, Neira M, et al. The impact of endocrine disruption: a consensus statement on the state of the science. Environ Health Perspect. 2013; 121: A104–106. pmid:23548368
  24. 24. Liem AKD, Theelen RMC. Dioxins: Chemical Analysis, Exposure and Risk Assessment. CIP-Gegecens Koninklijke Bibliotheek, Den Haag. 1997.
  25. 25. Kortenkamp A, Altenburger R. Synergisms with mixtures of xenoestrogens: A reevaluation using the method of isoboles. Sci Total Environ. 1998; 221: 59–73.
  26. 26. EFSA (European Food Safety Authority). Scientific Opinion on the hazard assessment of endocrine disruptors: Scientific criteria for identification of endocrine disruptors and appropriateness of existing test methods for assessing effects mediated by these substances on human health and the environment. EFSA J. 2013; 11 (3): 3132.
  27. 27. Payne J, Rajapakse N, Wilkins M, Kortenkamp A. Prediction and assessment of the effects of mixtures of four xenoestrogens. Environ Health Perspect. 2000; 108: 983–987. pmid:11049820
  28. 28. Silva E, Rajapakse N, Kortenkamp A. Something from "nothing"—eight weak estrogenic chemicals combined at concentrations below NOECs produce significant mixture effects. Environ Sci Technol. 2002; 36: 1751–1756. pmid:11993873
  29. 29. Rajapakse N., Silva E., Kortenkamp A. Combining xenoestrogens at levels below individual no observed- effect concentrations dramatically enhances steroid hormone action. Environ Health Perspect. 2002; 110: 917–921. pmid:12204827
  30. 30. Hass U, Scholze M, Christiansen S, Dalgaard M, Vinggaard AM, Axelstad M, et al. Combined exposure to anti-androgens exacerbates disruption of sexual differentiation in the rat. Environ Health Perspect. 2007; 115: 122–128.
  31. 31. Seeger B, Klawonn F, Bekale NB, Steinberg P. Mixture Effects of Estrogenic Pesticides at the Human Estrogen Receptor α and β. PLOS ONE. 2016; DOI:10.1371/journal.pone.0147490.
  32. 32. Frische T, Faust M, Meyer W, Backhaus T. Toxic masking and synergistic modulation of the estrogenic activity of chemical mixtures in a yeast estrogen screen (YES). Environ Sci Pollut Res (2009) 16:593–603
  33. 33. Filby AL, Neuparth T, Thorpe KL, Owen R, Galloway TS, Tyler CR. Health Impacts of Estrogens in the Environment, Considering Complex Mixture Effects. Environ Health Perspect. 2007; 115: 1704–1710. pmid:18087587
  34. 34. Roszko M, Szymczyk K, Rzepkowska M, Jędrzejczak R. Preliminary study on brominated dioxins/furans and hydroxylated/methoxylated PBDEs in Baltic cod (Gadus morhua) liver. Comparisonto the levels of analogue chlorinated co-occurring pollutants. Marine Poll Bull. 2015; 96: 165–175.
  35. 35. Roszko M, Rzepkowska M, Szterk A, Szymczyk K, Jędrzejczak R, Bryła M. Application of semi-permeable membrane dialysis/ion trap mass spectrometry technique to determine polybrominated diphenyl ethers and polychlorinated biphenyls in milk fat. Anal Chim Acta. 2012; 748: 9–19.
  36. 36. Roszko M, Szymczyk K, Jędrzejczak R. Simultaneous separation of chlorinated/brominated dioxins, polychlorinated biphenyls, polybrominated diphenyl ethers and their methoxylated derivatives from hydroxylated analogues on molecularly imprinted polymers prior to gas/liquid chromatography and mass spectrometry. Talanta. 2015; 144: 171–183. pmid:26452807
  37. 37. Roszko M, Szymczyk K, Jędrzejczak R. Separation of polychlorinated dibenzo-p-dioxins/furans, non-ortho/mono/di/tri/tetra-ortho-polychlorinated biphenyls, and polybrominated diphenyl ethers groups of compounds prior to their determination with large volume injection gas chromatography-Quadrupole ion storage tandem mass spectrometry. Anal Chim Acta. 2013; 799: 88–98.
  38. 38. Schulitz T, Metzger JW. Determination of estrogenic activity by LYES-assay (yeast estrogen screen-assay assisted by enzymatic digestion with lyticase). Chemosphere. 2004; 57: 1649–1655. pmid:15519410
  39. 39. Struciński P, Piskorska-Pliszczyńska J, Maszewski S, Góralczyk K, Warenik-Bany M, Mikołajczyk S, et al. PCDD/Fs and DLPCBs intake from fish caught in polish fishing grounds in the Baltic Sea–characterizing the risk for consumers. Environ Int. 2013; 56: 32–41. pmid:23567327
  40. 40. Julshamn K, Duinker A, Berntssen M, Nilsen BM, Frantzen S, Nedreaas K, et al. A baseline study on levels of polychlorinated dibenzo-pdioxins, polychlorinated dibenzofurans, non-ortho and mono-ortho PCBs, nondioxin- like PCBs and polybrominated diphenyl ethers in Northeast Arctic cod (Gadus morhua) from different parts of the Barents Sea. Mar Pollut Bull. 2013; 75: 250–258.
  41. 41. Rawn DFK, Breakell K, Vergin V, Nicolidakis H, Sit D, Feeley M. Persistent Organic Pollutants in Fish Oil Supplements on the Canadian Market: Polychlorinated Biphenyls and Organochlorine Insecticides J Food Sci. 2008; 74: 14–19.
  42. 42. Roszko M, Szymczyk K, Jędrzejczak R. Polychlorinated biphenyls (PCBs), polychlorinated diphenyl ethers (PBDEs) and organochlorine pesticides in selected cereals available on the Polish retail market. Sci Total Environ. 2014; 466–467: 136–151.
  43. 43. Preuss TG, Gehrhardt J, Schirmer K, Coors A, Rubach M, Russ A, et al. Nonylphenol isomers differ in estrogenic activity. Environ Sci Technol. 2006; 40:5147–5153. pmid:16955920
  44. 44. Hass U, Christiansen S, Axelstad M, Sřrensen KD, Boberg J. Input for the REACH-Review in 2013 on Endocrine Disruptors. Danish Center on Endocrine Disrupters, DTU Food National Food Institute. 2013.
  45. 45. Wang X, Yang H, Hu X, Zhang X, Zhang Q, Jiang H, et al. Effects of HO-/MeO-PBDEs on Androgen Receptor: In Vitro Investigation and Helix 12-Involved MD Simulation. Environ Sci Technol. 2013; 47: 11802–11809. pmid:24044724
  46. 46. Muhammad AS, Laessig RH, Reed KD. Polybrominated Diphenyl Ethers (PBDEs): New Pollutants-Old Diseases. Clin Med Res. 2003; 1: 281–290. pmid:15931321
  47. 47. Kortenkampa A, Altenburger R. Synergisms with mixtures of xenoestrogens: A reevaluation using the method of isoboles. Sci Total Environ. 1998; 221: 59–73.
  48. 48. Plísková M, Vondrácek J, Canton RF, Nera J, Kocan A, Petrík J, et al. Impact of polychlorinated Biphenyls Contamination on Estrogenic Activity in Human Male Serum Environ Health Persp. 2005; 113 (10): 1277–1284.
  49. 49. Hardrup N, Taxvig C, Pedersen M, Nellemann C, Hass U. Concentration addition, independent action and generalized concentration addition models for mixture effect prediction of sex hormone synthesis In Vitro. Plos ONE. 2013; 8(e70490): 1–13.
  50. 50. Brinkmann M, Hecker M, Giesy JP, Ratte HT, Hollert H, Preuss TG. Generalized concentration addition accurately predicts estrogenic potentials of mixtures and environmental samples containing partial agonists. Toxicol in Vitro. 2018; 46: 294–303. pmid:29066353
  51. 51. Suominen K, Hallikainen A, Ruokojärvi P, Airaksinen R, Koponen J, Rannikko R, et al. Occurrence of PCDD/F, PCB, PBDE, PFAS, and Organotin Compounds in Fish Meal, Fish Oil and Fish Feed Chemosphere. 2011; 85: 300–306.
  52. 52. EFSA (European Food Safety Authority). Scientific Opinion on Polybrominated Diphenyl Ethers (PBDEs) in Food. EFSA Panel on Contaminants in the Food Chain (CONTAM). EFSA J. 2011; 9 (5): 2156.
  53. 53. Klotz DM, Beckman BS, Hill SM, McLachlan JA, Walters MR, Arnold SF. Identification of environmental chemicals with estrogenic activity using a combination of in vitro assays. Environ Health Perspect. 1996; 104: 1084–1089. pmid:8930550
  54. 54. Garg N, Lata P, Jit S, Sangwan N, Singh AK, Dwivedi V, et al. Laboratory and field scale bioremediation of hexachlorocyclohexane (HCH) contaminated soils by means of bioaugmentation and biostimulation. Biodegradation. 2016; 27: 179–193. pmid:27142265
  55. 55. Kim TS, Kim CY, Lee HK, Kang IH, Kim MG, Jung KK, et al. Estrogenic Activity of Persistent Organic Pollutants and Parabens Based on the Stably Transfected Human Estrogen Receptor-α Transcriptional Activation Assay (OECD TG 455). Toxicol Res. 2011; 27: 181–184. pmid:24278570
  56. 56. Okoumassoun LE, Averill BD, Gagne F, Marion M, Denizeau F. Assessing the estrogenic potential of organochlorine pesticides in primary cultures of male rainbow trout (oncorhynchus mykiss) hepatocytes using vitellogenin as a biomarker. Toxicol. 2002; 178:193–207.
  57. 57. Longnecker MP, Gladen BC, Cupul-Uicab LA, RomanoRiquer SP, Weber JP, Chapin RE, et al. In utero exposure to the antiandrogen 1,1-dichloro-2,2-bis(pchlorophenyl)ethylene (DDE) in relation to anogenital distance in male newborns from Chiapas, México. Am J Epidemiol. 2007; 165 (9): 1015–1022.
  58. 58. Yang W, Mu Y, Giesy JP, Zhang A, Yu H. Anti-androgen activity of polybrominated diphenyl ethers determined by comparative molecular similarity indices and molecular docking. Chemosphere. 2009; 75: 1159–1164. pmid:19324393
  59. 59. Hartmann S, Larkus M, Steinhart H. Natural occurrence of steroid hormones in food. Food Chem. 1998; 62: 7–20.
  60. 60. Qin LQ, Wang PY, Kaneko T, Hoshi K, Sato A. Estrogen: one of the risk factors in milk for prostate cancer. Med Hypotheses. 2004; 62: 133–142. pmid:14729019
  61. 61. Malekinejad H, Rezabakhsh A. Hormones in Dairy Foods and Their Impact on Public Health—A Narrative Review. Iran J Public Health.2015; 44: 742–758. pmid:26258087
  62. 62. Behr M, Oehlmann J, Wagner M. Estrogens in the daily diet: In vitro analysis indicates that estrogenic activity is omnipresent in foodstuff and infant formula. Food Chem Toxicol. 2011; 49: 2681–2688. pmid:21801783
  63. 63. Pinto B, Garritano SL, Cristofani R, Ortaggi G, Giuliano A, Amodio-Cocchieri R, et al. Monitoring of polychlorinated biphenyl contamination and estrogenic activity in water, commercial feed and farmed seafood. Environ Monit Assess. 2008; 144: 445–453.
  64. 64. Garner CE, Jefferson WN, Burka LT, Matthews HB, Newbold RR. In Vitro Estrogenicity of the Catechol Metabolites of Selected Polychlorinated Biphenyls. Toxicol Appl Pharmacol. 1999; 154: 188–197. pmid:9925803
  65. 65. Sugihara K, Kitamura S, Sanoh S, Ohta S, Fujimoto N, Maruyama S, et al. Metabolic Activation of the Proestrogens trans-Stilbene and trans-Stilbene Oxide by Rat Liver Microsomes. Toxicol Appl Pharmacol. 2000; 167: 46–54.
  66. 66. Charles GD, Bartels MJ, Zacharewski TR, Gollapudi BB, Freshour NL, Carney EW. Activity of benzo[a]pyrene and its hydroxylated metabolites in an estrogen receptor-α reporter gene assay. Toxicol Sci. 2000; 55: 320–326. pmid:10828263
  67. 67. Nakagawa Y, Suzuki T. Metabolism of 2-hydroxy-4-methoxybenzophenone in isolated rat hepatocytes and xenoestrogenic effects of its metabolites on MFC-7 human breast cancer cells. Chem Biol Interact. 2002; 139: 115–128. pmid:11823001
  68. 68. KemI (Swedish Chemical Agency) The Swedish Chemicals Agency’s position on the possibility to determine threshold levels for endocrine disrupters. Memo. 2013, Reference 4.4b-H13-00416. 2013. pmid:23565126
  69. 69. Dybing E, Doe J, Groten J, Kleiner J, O'Brien J, Renwick AG, et al. Hazard characterization of chemicals in food and diet: dose response, mechanisms and extrapolation issues. Food Chem Toxicol. 2002; 40: 237–282. pmid:11893399
  70. 70. Caldwell DJ, Mastrocco F, Nowak E, Johnston J, Yekel H, Pfeiffer D, et al. An Assessment of Potential Exposure and Risk from Estrogens in Drinking Water. Environ Health Perspect. 2010; 118: 338–344. pmid:20194073
  71. 71. WHO. World Health Organization. Toxicological Evaluation of Certain Veterinary Drug Residues in Food. WHO Food Additive Series: 43. Estradiol-17β, Progesterone, and Testosterone. Geneva:World Health Organization, International Programme on Chemical Safety. Available: http://www.inchem.org/documents/jecfa/ jecmono/v43jec01.htm2000
  72. 72. Schneider S, Fussell KC, Melching-Kollmuss S, Buesen R, Gröters S, Strauss V, et al. Investigations on the dose–response relationship of combined exposure to low doses of three anti-androgens in Wistar rats. Arch Toxicol. 2017; 91(12): 3961–3989. pmid:28879601
  73. 73. De Falca M, Forte M, Laforgia V. Estrogenic and anti-androgenic disrupting chemicals and their impact on the male reproductive system. Front Environ Sci. 2017; 3: 1–12.