Multiresidue Determination of Endocrine Disrupting Compounds in Sewage Treatment Plants (SPE-HPLC-DAD)

This work focuses on the validation and application of solid phase extraction followed by high-performance liquid chromatography (SPE-HPLC) analysis of important endocrine disruptors compounds (EDC) from different classes in highly complex wastewater matrix. The endocrine disruptors investigated included three different categories: pharmaceuticals (sulfamethoxazole (SMZ), trimethoprim (TMP) and diclofenac (DCF)), hormones (estrone (E1), 17α-ethinylestradiol (EE2) and 17β-estradiol (E2)) and plastic materials (bisphenol A (BPA)). The method involves pre-concentration by SPE using Strata-X extraction cartridges followed by HPLC coupled with diode array detector (DAD). As the assumption of homoscedasticity was met for analytical data, ordinary linear regression procedure was applied to the data. The method was considered validated for 7 EDC after consistent evaluation of the key analytical parameters. Recoveries were ranged from 52.3 to 179.6%. Limits of quantification were in the range 6.0-0.4 μg L. The described method was applied to evaluate the occurrence and removal efficiency of EDC in two biological sewage treatment plants (STPs). The influent mean concentrations of E1, E2 and EE2 hormones were 79.54; 175.09 and 102.19 μg L, respectively, while for SMZ, TMP, DCF and BPA were 215.17; 187.01; 218.97 and 87.12 μg L for STP A, respectively. The efficiencies for the removal of EDCs ranged from 0 to 79.4%.


Introduction
Endocrine disruptors compounds (EDC) are components or active ingredients of pharmaceuticals, personal care products, veterinary medicines, "lifestyle compounds", steroids, metabolic regulators, food additives, and therefore, these chemicals are from multiple sources, and hence they are on the recent list of new contaminant candidates from US Environment Protection Agency and European Commission. [1][2][3][4] Most of technologies used in Brazilian's WWTPs (wastewater treatment plants) are conventional biological type, specially of stabilization ponds, activated sludge systems and anaerobic reactors, such as upflow anaerobic sludge blanket (UASB). However, several investigations have shown that EDC are frequently not removed during conventional biological wastewater treatment and also not biodegraded in the environment. [5][6][7][8] Their insufficient removal from the WWTPs is mentioned as the major source of their release into environment. 9,10 In the state of Ceará, biological processes are widely used, because of the great applicability, low capital and operational costs. However, there are still few studies on monitoring these compounds in sewage treatment plant effluents in northeastern Brazil, therefore, more studies on the determination and quantification of these compounds in these environmental matrices are essential. 5 A wide range of methods, mainly based on liquid chromatography-mass spectrometry (LC-MS), gas chromatography-mass spectrometry (GC-MS), solid phase extraction liquid chromatography tandem mass spectrometry (SPE-LC-MS/MS) are being developed. 5,9,11,12 Gas chromatography coupled to electron ionization (EI) mass spectrometry (MS) is both sensitive and selective for the determination of endocrine disruptors, however, derivatization of non-volatile and thermally labile endocrine disruptors is required prior to analysis, 13 besides when analyzing contaminated samples, suppression of electrospray ionization is most likely to occur. 14 This increases the overall analysis time and it may lead to errors to the analytical technique.
High-performance liquid chromatography (HPLC) instruments are one of the most popular and they are still widely applied to analyze these compounds because of their lower cost and greater robustness with sensitive method which can separate a wide variety of endocrine disruptor in variety matrices, some of which are not amenable to gas chromatography. 6 Liquid chromatography with fluorescence detection has been reported to have a low limit of detection (LOD). However, the technique certainly requires derivatization to improve the fluorescence properties for detection. 15 The use of diode array detector (DAD) as a detector for HPLC has proved to be a powerful tool in the determination and identification of compounds as it allows the online acquisition of their UV spectra. 3,4,[6][7][8] In addition, most of above mentioned methods are for one class of endocrine disruptors only. A challenge is presented in the simultaneous extraction and analysis of multiple classes of compounds.
The studies of selective and sensitive analytical methods are an important step to provide well founded data on the behavior of EDC in WWTPs. In order to optimize a treatment process or pre-treatment step so that the emission of undesired pollutants into the receiving waters is prevented, it is necessary more information obtained from analysis of influents and effluents of WWTPs. Based on the above, the goal of this paper is to develop and validate an SPE HPLC-DAD method for the qualitative and quantitative simultaneous determination of selected EDC in two two biological sewage treatment plants (STPs) as well as the evaluation of the occurrence of those EDC in two sewage treatment plants located in the state of Ceará, Brazil.
Sample preparation and HPLC-DAD analyses SPE extraction was performed following the methodology of Pessoa et al. 5 and qualitative and quantitative analyses of extracts collected were analyzed using HPLC.
The first step of sample preparation involved performing wastewater filtrations under vacuum through 0.45 µm glass-fiber filters to remove suspended particulate matter and avoid SPE cartridge clogging. The pH of each sample was then adjusted to 3.0 by addition of 50% (v/v) HCl, after which the analytes were extracted with a Speed Mate 12-port SPE vacuum manifold (Applied Separations, Allentown, PA, USA). The SPE cartridges were initially preconditioned with 10 mL of methanol, and subsequently with 10 mL of Milli-Q water (pH 3, Millipore, Bedford, MA, USA). The samples, typically 500 mL, were then loaded onto the 500 mg Strata-X cartridges (Phenomenex, Lane Cove, NSW, Australia) at a flow rate lower than 2 mL min −1 . The cartridges were then dried for 30 min under vacuum and eluted with 8 mL of methanol. Extracts collected were analyzed using HPLC.
Qualitative and quantitative analyses were performed using a liquid chromatograph HPLC Shimadzu (20A prominence, Kyoto, Japan) with UV-DAD detector (SPD-M20A) (215 nm), C18, 5 µm column, 250 × 4.6 mm inner diameter (Hichrom 5, Berkshire, UK), acetonitrile/HCl 0.1% as mobile phase. Analyses were carried out in a programmed gradient: increase of 10 to 100% acetonitrile in 10 min, returning to 10% in 4 min. The initial flow was 1.4 mL min −1 and after 5 min the flow was increased to 2.0 mL min −1 . The column temperature was kept at 35 °C.

Selectivity
The interference of matrix compounds was assessed by retention times (t R ), UV spectra, and peak purity tests for all endocrine disruptors compounds in both standard solutions and in the fortified matrix. The parallelism between the calibration curves obtained in both solvent and spiked matrices was also considered.  The range of linearity to the target compounds was evaluated based on the steps reported in Barbosa and coworkers' study. 16,17 Visual inspection of the residual plots was performed and in cases of inconclusive visual inspection, it was performed Hartley's F max test to evaluate homoscedasticity. The estimation of the regression parameters by ordinary least squares (OLS) approach requires that the regression residuals follow a normal distribution. In case of heteroscedasticity, the equations of the calibration curves and correlation coefficient (r) were obtained by the method of weighted least squares (WLS). 16,17 Finally, the evaluation of the adjustment of the calibration equations was conducted according to the Student's t r -test for lack of fit.
The Hartley's F max test was conducted by calculus of the F statistics by equation 1 and compared with the critical value F (f = m − 1; k) of Hartley's F table. 18,19 (1) where and are the largest and the smallest variances of the regression, respectively. F values less than or equal to the F critical value show that the regression residuals are homoscedastics.
The t r -test to verify the linear model fit for the calibration curves was performed by calculus of t statistic by equation 2 and compared with the critical value t crit to (n − 2) degrees of freedom and 95% of level of confidence, where r is the correlation coefficient and n is the number of calibration levels. When the t r values were larger than or equal to the critical values, the equation was considered adequate for the calibration curves. [20][21][22] (2)

Significance of regression parameters
After obtaining the adjusted calibration equations it was finally assessed the statistical significance of the regression parameters using the t-test. The calculations were carried out for the t-values for each regression parameter by equations 3 and 4. 19 where b is the y-intercept; s b is the standard deviation of b-values; a is the slope of the calibration curve; s a is the standard deviation of a-values The equations 5 and 6 were used in the calculations of the standard deviations for the parameters estimated by OLS regression. 19 The term s y/x is the standard error of regression and is the the average of the values.
Limit of detection (LOD) and limit of quantification (LOQ) The estimative of limit of detection (equation 7) and limit of quantification (equation 8) were determined based on the parameters of adjusted calibration curves, according to ICH guideline. 23 where σ is the standard deviation of the response based on the standard deviation of the blank; a is the slope of the calibration curve.

Accuracy and precision
The precision of the method was assessed by repeatability. Accuracy was assessed by recovery experiments running with three replicates of spiked samples in levels 0.8, 1.6 and 16.0 µg L −1 for all compounds. The repeatability was reported by relative standard deviation (RSD, in percentage) and accuracy by recovery percentage (Rec), calculated by equations 9 and 10, respectively. In equation 9, the term s r is the standard deviation of repeatability.
where C 1 is the average concentration determined in the fortified sample; C 2 is the concentration determined in the unfortified sample (blank sample); and C 3 is the added concentration to the sample.

Sample collection
The described method was applied to the determination of endocrine disruptors in wastewater samples from wastewater treatment plants. Duplicate grab wastewater samples were collected between September 2013 and March 2015 from two full-scale WWTPs located in the state of Ceará, a semi-arid zone in Brazil. The two STPs (A and B) comprised upflow anaerobic sludge blanket (UASB) with post-chlorination disinfection. Additional details regarding these STPs are given in Table 2. For each STP, five influent and effluent samples were analyzed in order to determine the estrogens removal.
The samples were collected in glass flasks, transferred to 1 L amber glass bottles and preserved by the addition of 10 mL L −1 formaldehyde. Then, the samples were transported to the laboratory in cooling boxes and were prepared for analysis within 48 h. Formaldehyde was also used in the method validation, and no interference was observed. 5,24

Results and Discussion
Calibration curves

Selectivity
The selectivity was observed by association of the t R , peak purity tests, UV spectra, and the parallelism of the calibration curves. Figure 1 shows the overlap chromatograms from prepared samples obtained in matrix and 160 µg L −1 spiked matrix. It can be observed no interference for the seven studied compounds.
The DAD detects the absorption in ultraviolet (UV) to visible (Vis) region, it has multiple photodiode arrays to obtain information over a wide range of wavelengths at one time. If two analytes have different spectra, it is possible to distinguish them using HPLC-DAD even if they have similar t R . Another common DAD application is in determination of the compound purity, wherein the software calculates the rate of absorption through the peak.
The peak purity tests performed by the HPLC-DAD software revealed that all peaks had their purity levels higher than 99%, independent of the matrix ( Table 3).
The similarity is calculated from the comparison of the compound's absorption spectra in the matrix and in standard solutions. The resulting spectrum is a graph where absorption versus wavelength (nm) is related. Absorption in UV-Vis of organic molecules is restricted to certain functional groups (chromophores), which contain valence electrons with low excitation energy. UV-Vis spectra can be used for the qualitative identification of molecules and atomic species, since they are characteristic of a certain molecular structure.
In this study, it was possible to observe similarity level higher than 96% for all studied compounds. It can be seen in Figure 2 the unambiguous comparisons of the overlapping absorption spectra for the compounds studied in matrix with    the standards. Then, these observations guarantee that the method is selective for the seven endocrine disruptors and can be used for monitoring purposes of wastewater samples.

Sensibility, linearity and range
The residual plots visual inspection and Hartley's F max test indicated lack of fit in linear range from 0.8 to 160 µg L −1 for all studied compounds, however, for the 1.6-80.0 µg L −1 considered range, it was possible to observe that OLS was more appropriate than WLS method for BPA and SMZ compounds ( Table 4). The E2 compound showed homoscedastic behavior in a wider linear range (1.6-128 µg L −1 ), as well as TMP compound (0.8-80 µg L −1 ). For DCF, E1 and EE2 compounds, 8.0-160.0 µg L −1 considered range it was more appropriate. The different behavior indicated the need for careful selection of the most adequate model to express the analytical curve and the need for assuring the reliability of the used method. Figure 3 shows simple residuals plots of two compounds.
The linear models for calibration curves showed good fit for all compounds in the selected range, confirmed by the coefficients of correlation r ranging from 0.9926 to 0.9999 and t r -test, which indicated that t r -values were higher than t crit values for (n − 2) degrees of freedom at 95% confidence level. 21,22 The analytical method showed higher sensibility for TMP and SMZ compounds, which can be seen by the slope values from calibration curves in Table 4. This might happen due to the high amount of chromophore groups (unsaturated organic groups that absorb in the ultraviolet and visible regions) present in the molecular structure of TMP and SMZ compounds compared to the other studied compounds besides the presence of nitrogen and sulfur atoms, which have nonbonding electrons in their structures, giving them a greater capacity for absorption of radiation in the range of 170 to 250 nm, then increasing the UV-Vis detector sensibility for TMP and SMZ compounds. At the same time, the hormones E1, E2 and EE2 showed the lower values, probably due to the presence of only one phenol and hydroxyl groups in their molecular structures with most of its structures consisting of single bonds, C−C and C−H, which do not absorb radiation in the UV-Vis. 14  Despite DCF and BPA compounds also presenting chromophore groups in its molecular structure, the DAD detector showed lower sensibility to them compared to TMP and SMZ. This might happen since DCF and BPA have high logK ow (octanol-water partition coefficient) values, which might imply to them a higher affinity to SPE cartridges in the sample preparation step, decreasing the efficiency of methanol extraction step and then decreasing the amount of these compounds present in the sample, thereby decreasing the intensity of absorption for such compounds since quantity of material is directly proportional to the absorption intensity detector. Table 5 shows the t-test results for significance of regression parameters: slopes and intercepts. It can be seen that the t calculated values for most of the regression parameters were greater than the critical t values, with exception for E1, EE2, DCF and TMP compounds using 95% confidence level and (n − 2) degrees of freedom. These corrected calibration equations were used to predict analyte concentrations in subsequent validation studies as in the analyzed wastewater treatment samples.

Limit of detection (LOD), limit of quantification (LOQ), accuracy and precision
The estimated LOD and LOQ for the analytes are reported in Table 6. The LOQ values were 6.0-0.4 µg L −1 . Despite the high values of LOQ found in the present study, when compared to other studies using more sensitive techniques, such as LC-MS/MS, 9 the values were justified, since the environmental samples analyzed showed high levels of the compounds studied.
The values found were similar to those found by Pessoa et al., 5 which studied the hormones E1, E2 and EE2 using derivatization techniques for sample preparation followed by GC-MS analysis. These authors found values of 0.12; 0.04 and 4.7 µg L −1 , respectively.
Almost all EDCs showed recovery rates of approximately 68-106%, to the lowest fortification level. All analytes showed RSD below 24%.
SMZ showed low recovery levels, such as 18.78 and 23.89% for 16 and 1.6 µg L −1 levels, respectively. According to Vidal et al., 6 the low recoveries of SMZ compound by Strata-X cartridges are related to the weak interaction  between SMZ compounds and the Strata-X cartridges polymeric material. Furthermore, the SMZ has the higher solubility comparing to the other studied compounds, and then this endocrine disruptor seems to have more affinityfor aqueous instead of solid phase. TMP compound had different behavior, showing higher recovery values between 205.53-158.22%. For this compound the matrix effect might be more pronounced. Accuracy, in the form of recovery measures (in percentage), is always considered adequate within certain limits (associated with precision values). In cases where the analyte is found at trace level in the sample, as in this study (0.80, 1.6 and 16 µg L −1 ), and there is analytical complexity, the minimum and maximum recovery values are less restrictive than the traditional 70-120% range, being considered as reasonable recovery values between 50 and 120%, with a precision of up to 15%. [25][26][27] It was not developed specific studies to estimate the matrix effect in the present research. However, factors associated with the nature of the analyte (physicalchemical properties) and the chromatographic system (injection system, column and detector) might influence the increasing of chromatographic response induced by the matrix, explaining recovery levels above 120% for some analytes in the present study. The mechanisms related to the chromatographic system are not clear, there are few studies 25 that deal with the matrix effect in the HPLC-DAD system.

Occurrence of endocrine disruptors in influent and effluent samples
It is well known that most of the endocrine disruptors are introduced to environment from the discharge of wastewater treatment plants since these treatments technologies are not specifically designed to remove these compounds. As a consequence, it is extremely important to investigate and understand the fate of EDC within treatment facilities in order to increase removal efficiency and prevent their disposal in surface waters and the contamination of drinking water and ground water.
The distribution of each endocrine disruptors concentrations which were measured in raw and treated (before and after chlorination step) wastewater from two STP (n = 26) is shown in Figures 4-10.

Removal of hormones
According to Figure 4, it is possible to see that neither STP A nor B is efficient to remove E1 compound. In both studied STPs, the mean in the effluent (102.07 and 26.44 µg L −1 , respectively) were higher than influent (79.54 and 25.27 µg L −1 , respectively), this might due to the conversion of E2 and EE2 into E1 before it could be transformed further. [28][29][30][31] It can also be noticed in STP A, an interesting behavior of E1 concentration level through the treatment steps. After UASB process there was removal in E1 concentration, however, after chlorine disinfection step there was an increase in that concentration. This probably might be due the oxidation-reduction (redox) reaction, since E1 compound under anaerobic conditions (R=O) might be reduced to E2 compound (R−OH), which might return to the oxidized form (E1) after chlorine disinfection process.
In terms of E2 compound, the maximum concentration found in STP A was 613.   Figure 6). In terms of total removal, the STP B proved to be more efficient than the STP A, with a 79.4% removal.
The concentration values for the estrogenic steroids (E1, E2 and EE2) found in the present study is higher than those found in the literature, 5 which evaluated the occurrence and removal efficiency of four estrogenic hormones in five biological STPs, located in the state of Ceará, Brazil. The authors found mean influents concentrations of 0.57, 0.14 and 0.42 µg L −1 for E1, E2 and EE2, respectively.
Out of Brazil, average concentrations of estrogenic steroids in effluent of sewage treatment plants (STPs) are lower than this study. In Italian STPs it was found 0.08, 0.01, and 0.05 µg L −1 , for E1, E2 and EE2, respectively. The concentrations of E2 in influents of Japanese STPs ranged from 0.03 to 0.09 µg L −1 and 0.02 to 0.1 µg L −1 in autumn and summer seasons, respectively. 10

Removal of pharmaceuticals
For the TMP antibiotic, the maximum concentration found in the STP A was 673.1 µg L −1 , with the mean values of the influent, post-UASB and effluent of 187.01, 79.44 and 60.43 µg L −1 , respectively ( Figure 7). In terms of total removal, a value of 67.7% was found. For STP B, the   (Figure 7). In terms of total removal, STP B was not efficient.
For the SMZ antibiotic, the maximum concentration found in the STP A was 619.26 µg L −1 , with the mean values of the influent, post-UASB and effluent of 215.17; 161.53 and 71.98 µg L −1 , respectively (Figure 8). In terms of total removal, a value of 66.6% was found. Different behavior was found in STP B, where there was initially an increase in SMZ concentration after UASB reactor and small removal after chlorination step. Thus, there was a global removal of 32.3%.
For the anti-inflammatory DCF, the maximum concentration found in the STP A influent was 594.27 µg L −1 , with influent, post-UASB and effluent mean concentrations of 218.97; 201.61 and 235.34 µg L −1 , respectively ( Figure 9). In terms of total removal, STP B was not efficient. The DCF maximum value of 153.65 µg L −1 was found for the STP B, with the influent, post-UASB and effluent mean concentrations of 71.50; 18.39 and 28.81 µg L −1 , respectively ( Figure 9). In terms of total removal, the STP B proved to be much more efficient than the STP A, with a removal percentage of 59.7%.
Rivera-Jaimes et al. 32     terms of total removal, a value of 24.9% was found. For STP B, treatment was better than STP A, with removal of 69.3%.

Conclusions
Although the occurrence of endocrine disruptors in the wastewater matrices has been extensively investigated in developed countries, it needs to be further studied in developing countries, where the proportion of non treated wastewaters is still very large. This study investigates the occurrence of seven endocrine disruptors from different chemical classes in sewage of the Fortaleza region, located at state of Ceará (Brazil) during two consecutive years. Seven investigated EDCs were consistently found at relevant concentrations (reaching the µg L −1 level) in all STPs studied.
The removal efficiency behavior of the two STPs were different, despite both have the same biological treatment technology, due the differences between the STPs, such as hydraulic holding time, type of inoculum, and EDC influent concentration. Even after biological treatment, the effluents still present considerable amounts of all EDCs, even in high concentrations when compared to those reported in the literature, requiring a post-treatment step for the elimination of these compounds from the effluents before they are discarded into the water bodies.

Author Contributions
Carla B. Vidal was responsible for the conceptualization, data curation, formal analysis, funding acquisition, investigation, methodology, project administration, resources and writing original draft; Pablo G. A. Barbosa for the data curation, formal analysis and writing review and editing; Germana P. Pessoa for the conceptualization and data curation; Patrícia C. Buarque for the conceptualization and data curation; José G. S. Nascimento for the data curation and writing original draft; Antônio L. Farias Filho for the data curation; Mário S. Paz for the formal analysis and writing original draft; André B. dos Santos for the funding acquisition, project administration, resources and supervision; Rivelino M. Cavalcante for the writing review and editing; Ronaldo F. Nascimento for the conceptualization, project administration, resources, writing original draft, review and editing.