Confirming the presence of selected antibiotics and steroids in Norwegian biogas digestate

Farms utilizing sewage sludge and manure in their agronomic plant production are recognized as potential hotspots for environmental release of antibiotics and the resulting promotion of antibiotic resistance. As part of the circular economy, the use of biogas digestates for soil fertilizing is steadily increasing, but their potential contribution to the spreading of pharmaceutical residues is largely unknown. Digestates can be produced from a variety of biowaste resources, including sewage sludge, manure, food waste, and fish ensilage. We developed a method for the detection of 17 antibiotics and 2 steroid hormones and applied the method to detect pharmaceutical residues in digestates from most municipal biogas plants in Norway, covering a variety of feedstocks. The detection frequency and measured levels were overall low for most compounds, except a few incidents which cause concern. Specifically, relatively high levels of amoxicillin, penicillin G, ciprofloxacin, and prednisolone were detected in different digestates. Further, ipronidazole was detected in four digestates, although no commercial pharmaceutical products containing ipronidazole are currently registered in Norway. A simplified risk assessment showed a high risk for soil microorganisms and indicates the tendency for antibiotic-resistant bacteria for penicillin G and amoxicillin. For prednisolone and ipronidazole; however, no toxicity data is available for reliable risk assessments. Graphical abstract Supplementary Information The online version contains supplementary material available at 10.1007/s11356-022-21479-1.


Introduction
Biogas production is increasingly used in circular bioeconomic strategies as a sustainable strategy for organic waste management. Currently, more than 18,000 biogas plants are registered in Europe (Cesaro, 2021). A variety of organic wastes such as sewage sludge and manure, as well as food and household wastes, are anaerobically digested for the production of biogas, containing mainly methane and carbon dioxide (the methane is thereafter used as an energy carrier). The nutrient-rich organic residue from that process is known as biogas digestate. While the use of biomethane for transport is well-known to reduce greenhouse gas emissions, Lyng et al. (2018) showed that the replacement of mineral fertilizer with biogas digestate is equally important for a reduced carbon footprint. Furthermore, the mineral phosphorus resources are finite with a proposed peak around 2030, and today, only 20% of all mined phosphorus is consumed in food (Childers et al. 2011). Thus, increased reuse of nutrients is important.
Unfortunately, in addition to the valuable nutrients, digestates have been shown to contain residues of legacy and emerging organic pollutants which are not fully degraded during the digestion process (Ali et al. 2019;Lindberg et al. 2005;Spielmeyer et al. 2014;Suominen et al. 2014;Widyasari-Metha et al. 2016). These anthropogenic contaminants may ultimately enter the agricultural production system if the digestates are spread on agricultural land (Chen et al. 2019).
Residues of antibiotics and steroids have been detected in sewage sludge, sewage effluents, manure, and other environmental matrices (Chang et al. 2007;Clarke and Smith 2011;Spielmeyer 2018;Verlicchi and Zambello 2015). Both antibiotics and steroid hormones, such as oestrogens and glucocorticoids, have wellknown biological activity and may cause adverse environmental effects. For instance, continuous exposure of soil microbial communities to antibiotic agents by, e.g. manure, spread on agricultural land may lead to elevated levels of antibiotic resistance genes and resistant bacteria (Heuer et al. 2011). Further, the presence of antibiotics can disrupt the natural soil microbial flora and thereby adversely affect biogeochemical processes, such as nitrification, denitrification, and iron reduction (Grenni et al. 2018;Roose-Amsaleg and Laverman 2016). The exposure of aquatic organisms to glucocorticoids can lead to changes in behaviour and immunological responses, as demonstrated in previous studies (e.g. Bal et al. 2017;McNeil et al. 2016).
However, the fate of these contaminants in biogas digestate has only been sparsely investigated, due to lack of appropriate methods. Hence, we developed and validated a multi-compound quantitative trace level analytical method for the investigation of 16 antibiotics and steroid hormone residues in digestates. The choice of compounds was based on a former screening programme of the Norwegian Food Safety Authority. The list includes compounds considered important for food safety in Norway either due to frequent application in Norwegian agriculture or because they are banned for use in Norway according the EU regulation 37/2010/EC (NFSA 2015). The here developed method was applied to representative samples from 12 centralized municipal biogas plants (i.e. most of such plants in Norway) as well as two experimental reactors. Based on the detected pharmaceutical levels, a simplified risk assessment was performed by calculating expected soil concentration caused by the application of digestate and comparing those to ecotoxicity and antibiotic resistance development data.

Collection of biogas digestates and information about the operating conditions
Biogas digestates were collected from twelve major municipal biogas plants in Norway (plants A to L, Tables 1 and S10), as well as two experimental units connected to research stations (plant M and plant I exp , connected to plant I). Each plant sampled about 1 L of their digestate and sent it to the Norwegian University of Life Sciences for quantitative analysis. Detailed information about the operating parameters, as supplied by the biogas plants, can be found in Table S10. The feedstocks used were food waste (E, G, and K), sewage sludge (D, H, J, and I), food waste mixed with sewage sludge (A, B, F, and L), manure mixed with food waste (C), manure (M), and a manure/fish silage combination (I exp ). The digestates are referred to as liquid (< 5% dry matter, subscript L) or solid (20-50% dry matter, subscript S). Some plants produce both fractions; thus, the total number of digestates analysed was 18. Samples were quantified in duplicate for each digestate batch to account for method uncertainty and for the heterogeneity of the material. The same samples were previously analysed for pharmaceuticals and personal care products (PPCPs) by Ali et al. (2019).

Target analytes
For the method validation process, 26 compounds from 8 classes of antibiotics and steroids were chosen as target analytes (Table 1). The target analytes have different linear concentration ranges as determined by the individual calibration curves, as well as different shelf life as standard solution. Based on the individual detector sensitivity and preservability, the compounds were divided into three groups: X, Y, and Z (Table 1). Separate solutions were prepared for each group, for both the native 12 C standards (STDs) and the isotopelabelled internal standards (ISTDs). A complete list of target compounds including structure information and CAS registry numbers is available in Table S1. The stock solution concentrations and compositions are summarized in Table S2.

Chemicals and solutions
Methanol (MeOH, HPLC-grade) and acetonitrile (ACN, HPLC-grade) were purchased from VWR (West Chester, PA, USA). Formic acid, ammonium acetate, disodium ethylene diamine tetra acetate (Na 2 EDTA), citric acid, sodium phosphate dibasic, and phosphoric acid were purchased from Sigma Aldrich (Oslo, Norway). Only grade 1 purified water from Milli-Q water purification systems (Millipore, Bedford, MA, USA) was used for the sample preparation and analysis. All standards and internal standards were purchased according to Table S1 in the Supplementary information. Separate solutions were prepared for the X-, Y-, and Z-compound groups both for the native 12 C standards (STDs) and the isotope-labelled internal standards (ISTDs) ( Table S2).

Sample preparation and clean-up
Based on a previously published method (Hu et al. 2010), a comprehensive sample clean-up protocol was developed and validated ensuring minimum of matrix disturbances and coelution in the final LC-MS/MS quantification (Fig. 1). An aliquot of 2 ± 0.03 g (wet weight, ww) digestate sample was weighed into 15 mL polypropylene tubes (Fig. 1). Internal standards were added, corresponding to 7.5 ng X-compound, 75 ng Y-compounds, and 150 ng Z-compounds (Table S2). Subsequently, 3 mL of extraction solution (Table S3) was added before the sample was vortexed for 10 s, ultrasonicated for 10 min, and centrifuged at 3500 rpm for 15 min. The supernatant was transferred to a glass tube, and the extraction   Table 2 Validation results of all successfully validated compounds, with recovery, method repeatability, method accuracy, matrix effect, and efficiency of the extraction method (EEM) at level 3 of the calibration curve (i.e. 7.5 ng mL −1 for X-compounds, 75 ng mL −1 for Y-compounds, and 150 ng mL −1 for the Z-compounds, see Table S5). procedure was repeated twice without ultrasonication. The combined supernatants were dried to approximately 1 mL with controlled heating at 37 °C under a stream of compressed air (analytical quality, AGA, Porsgrunn, Norway) using a Reacti-Therm III evaporator (Thermo Fisher Scientific Inc., Rockford, USA). The samples were subsequently shaken for 10 s with 4 mL of added grade 1 water. The extract, as well as 1 mL of water used to rinse the sample tube, was passed through a Bond Elut SAX column (500 mg, 3 mL) (Agilent, Santa Clara, USA) preconditioned with 2.5 mL MeOH and 2.5 mL grade 1 water, using light vacuum (water jet). The eluted extract was added to an Oasis HLB solid-phase extraction (SPE) column (200 mg, 6 mL) (Waters, Milford, USA) which was preconditioned with 5 mL MeOH and 5 mL water.
The column was washed with 3 mL of SPE buffer (Table S3) and 3 mL of 5% MeOH in water. The sample was eluted with 3 mL MeOH and then evaporated to dryness at 37 °C. The sample was then reconstituted in 1 mL 20% MeOH in water and vortexed and filtered through a 0.2-μm microcentrifuge filter (Spin-X, Costar, Corning Inc. NY, USA) before the samples were transferred to 2-mL glass vials for quantitative analysis using LC-MS/MS. Triple quadrupole (QqQ) dynamic multiple reaction monitoring (dMRM) with electrospray ionization (MRM-ESI) was applied for quantitative analysis (Fig. 1).

Analysis
Compound-specific chromatographic separation and quantitative detection of the cleaned digestate extracts were conducted on an Agilent 1260 ultra-high performance liquid chromatograph (UHPLC; Agilent Technologies, Waldbronn, Germany), and the detection and quantification were done on an Agilent 6490 triple quadrupole mass spectrometer (Agilent Technologies, Santa Clara, CA, USA) with an Agilent Jet Stream electrospray ion source, using dynamic multiple reaction monitoring (dMRM). For instrument control and method validation and quantification, the Agilent MassHunter software (V B.07.00/Build 7.0.457.0, 2008) was used. All MRM transitions and details from the MassHunter method, as well as the ion source parameters, are described in Tables S3 and S4 in the Supplementary information section.

Method validation
For the method validation, detection limits, recovery, repeatability, accuracy, matrix effect, and efficiency of the extraction method were calculated. The details are given in the Supplementary information. Digestate are a complex matrix that, if not properly treated in advance, would affect the instrumental sensitivity and selectivity of each target analyte differently. To compensate for these variations, a matrix-matched calibration curve was chosen. The calibration curve was prepared with five concentration levels spanning the ranges given in Table 1, as well as one level with no added native compounds (level 0, see Table S5). Compounds were approved if the recovery rate was 40-115%, if the relative coefficient of variation (CV%) for the repeatability was < 15%, if the accuracy was < ± 15% with a CV% < 15%, and if the determination coefficient (R 2 ) of the matrix-matched calibration curve was above 0.985 (Table S6).

Simplified risk assessment
For each of the quantified pharmaceuticals, the highest concentration found in any of the biogas digestates (C DIG ) was used to calculate a predicted environmental concentration (PEC) in soil immediately after application of biogas digestate (Eq. 1). Application rates of digestate vary, due to, e.g. differences in nutrient and heavy metal content. Reasonable estimates according to the Norwegian Agricultural Extension Service is however 8 tonnes of solid or 45 tonnes of liquid digestates (fresh weight) per hectare (personal communication). Assuming an incorporation to a soil depth of 20 cm, and a dry soil bulk density of 1.3 kg l −1 (RHO SOIL ), this corresponded to an application rate of 3.1 or 17.3 g fresh digestate kg −1 dry soil, respectively.
A predicted no-effect concentration (PNEC) for toxicity towards soil organisms was calculated following the standard approach of risk assessments from the European Chemicals Bureau (European Commission 2003): the lowest available EC50 or NOAEL value for soil organisms was divided by an assessment factor of 10-1000, depending on the availability of toxicity data. Tables S11 and S12 summarize the ecotoxicity data used for the PNEC calculations. For prednisolone, no soil toxicity data was available, and the PNEC SOIL was derived by multiplying the PNEC AQUATIC with the soil partitioning coefficient (K d ); see Supplementary information for details. PNEC values for selection of antibiotic-resistant bacteria (ARB) were derived from Menz et al. (2019) by Eq. 2: where PNEC ARB,PW is the PNEC value in soil pore water, K d is the soil-water distribution coefficient, and f(pw) is the fraction of pore water in the soil (0.25). (1) (2) The risk quotient (RQ) was calculated as where a RQ < 0.1 indicates low risk, RQ = 0.1-1 indicates medium risk, and RQ > 1 indicates high risk.

Method validation
Out of a total list of 26 selected target contaminants, 16 compounds were found to be valid for quantification (see Table 2). The detection and quantification thresholds for the validated compounds were highly to medium sensitive (in the pg g −1 to ng g −1 range), and their individual response curves were ranging over three orders of magnitude. The compounds are from six different compound classes, making this a multiclass quantification method appropriate for the screening in biogas digestate. Some of the compounds have lower recoveries than preferred, but since the quantification at ultra-trace level of these target chemicals (pg-range) is associated with an estimated overall method uncertainty between 40 and 50%, recovery rates between 40 and 115% were accepted. The results for all 26 compounds can be found in Table S6. When the matrix effect is ± 20%, no significant matrix effect is assumed. However, for most of the compounds, there was significant matrix effect confirmed, either signal enhancement (> 20%) or ion suppression (< − 20%). This demonstrates the need for a matrix-matched calibration curve. The efficiency of the extraction method was below 40% and above 115% for some compounds, demonstrating the need of well-defined internal standards. This also illustrates a considerable challenge for multi-compound quantification methods for trace level analysis. The deviating extraction efficiency of the compounds is all acceptable as long as the recovery is within the acceptable range, i.e. 40-115%.

Levels of antibiotics and steroids in biogas digestate
Our survey revealed overall low levels of the target substances in the digestates. From the 16 target analytes, 8 were detected above the method quantification limit (MQL) ( Table 3). Four were not detected in any of the samples, namely, ronidazole, enrofloxacin, sarafloxacin, and tiamulin. The antibiotics metronidazole, norfloxacin, difloxacin, and sulfadoxine and the glucocorticoid dexamethasone were detected in at (3) RQ = PEC∕PNEC least one biogas digestate, however, each below their respective MQLs. Despite overall low levels, a few findings were of concern. Amoxicillin, penicillin G, ciprofloxacin, and prednisolone were found at levels above 400 µg kg −1 dw. Also, ipronidazole was found in trace amounts in several digestates, even though the pharmaceutical is not registered for use in Norway.
β-lactams Amoxicillin and penicillin G were found at 460-960 µg kg −1 dw in the food waste digestates E S and E L and the manure digestate M, respectively (Table 3). The detection of both antibiotics was surprising, as β-lactams are expected to rapidly degrade during biological processes by, e.g. hydrolysis (Braschi et al., 2013). The hydrolytic half-life of penicillin G is 60 h at 37 °C (Chadha et al., 2003), which is comparable to the temperature of reactor M (i.e. 35 °C). Consequently, β-lactams are rarely reported in manure despite their common use in animal husbandry. Hence, none of the recent reviews on antibiotics in manure by Spielmeyer (2018) and Wohde et al. (2016) reported the presence of either amoxicillin or penicillin G. Furthermore, the detection of amoxicillin in a digestate from food waste only was unexpected. Globally, antibiotics are frequently detected in animal products and vegetables, even as high as 1500-3000 μg kg −1 fw (in cultivated fish, China and Turkey, reviewed by Chen et al. (2019) Penicillin G is associated with a high risk quotient both towards soil bacteria (2.2, Table 4) and for development of antibiotic-resistant bacteria (ARB, 14.6), while the evaluation of amoxicillin resulted in lower risk quotients. The difference in risk quotients for ARB selection for amoxicillin and penicillin G is mainly due to the large difference in K d values used in Eq. 2. These were estimated from their organic carbon partitioning constants (K OC ) reported in a review by Cycoń et al. (2019); as for our study, no experimental K d was available (details in Supplementary information). The K OC used were 865.5 l kg −1 and 2.68 l kg −1 for amoxicillin and penicillin G, respectively. Consequently, penicillin G is predicted to sorb less to the soil and be as more bioavailable, as reflected in a lower PNEC value for penicillin G than for amoxicillin. Fluoroquinolones Ciprofloxacin is one of the most examined and detected antibiotics worldwide (e.g. Verlicchi and Zambello 2015) and accounts for approximately 90% of the human consumption of quinolones in Norway Table 3 Concentration of selected antibiotics and steroid hormones in the digestates [μg kg −1 wet weight] (μg kg −1 dry weight is given in parenthesis, n = 2 except for E L with n = 4). Operating parameters of the biogas plants A-N can be found in Table S10 , and their degradation during anaerobic digestion is limited (Golet et al. 2003;Lindberg et al. 2006;Zhang and Li 2018). In fact, the concentration may even increase during thermal hydrolysis, probably due to release of intracellular antibiotics (Zhang and Li 2018). Ciprofloxacin can therefore be expected to be found in sewage sludge digestate, and the amount of 430 μg kg −1 dw found in digestate D S in our study corresponds well with the range detected in sewage sludge by TemaNord (2012). A positive result was thus the lack of ciprofloxacin in 7 out of 8 biogas plants receiving sewage sludge, perhaps explained partly also by a decline in ciprofloxacin prescriptions of approximately 50% since 2012 as a strategy to prevent ciprofloxacin resistance (Sakshaug et al. 2017;Sommerschild et al. 2020).
In soils, fluoroquinolones are strongly sorbed to clay particles and organic matter, limiting their bioavailability and thus reducing their impact on soil biota and processes such as nitrogen transformation (Rosendahl et al., 2012). On the other side, the low availability combined with low hydrolysis and thermal and biological degradation (Al-Ahmad et al. 1999;Alexy et al. 2004;Thiele-Bruhn 2003) leads to their persistence in soil. Thus, they can accumulate when repeatedly added, as shown by Dalkmann et al. (2012), who found accumulation of ciprofloxacin in soils irrigated with wastewater. Further, Girardi et al. (2011) found that ciprofloxacin could inhibit soil respiration despite the formation of nonextractable residues. The risk quotient for ciprofloxacin was close to 1, indicating a moderate risk towards soil organisms.
It should be noted that ciprofloxacin was measured at a concentration more than twice as high as the upper boundary of the calibration curve. As the measured concentration of 430 µg kg −1 fw is uncertain, the upper boundary of 205 µg kg −1 fw was used for PEC SOIL calculations. Thus, the risk quotient is probably underestimated. Considering its effect on soil organisms, persistence in soil, and the frequent detection in other studies, further investigations on the levels of ciprofloxacin in biogas digestates are needed.
Sulfonamides Sulfonamides, represented by sulfadiazine and sulfamethazine, were found in several digestates. The concentrations were low, leading to low risk quotients. Neither sulfadiazine nor sulfamethazine was found in sewage sludge in the Nordic countries in 2012 (TemaNord 2012), probably because of their high hydrophobicity and low solid-liquid partition coefficients (K d ) and because their negative charge at high pH hinders electrostatic sorption to the negatively charged surfaces of sludge (Göbel et al. 2005;Zhang and Li 2018). Sulfonamides are also used in veterinary medicine, but the detection frequency in German digestates from swine manure has been low (Spielmeyer et al. 2014).
Nitroimidazoles Ipronidazole was found at trace levels in four digestates from three different plants. Its presence is nevertheless surprising and concerning, as there are no registered pharmaceuticals in Norway for humans or animals containing ipronidazole (Østensen, H., personal communication, 16.07.2021). The digestates were based on sewage sludge alone (J S ) or in combination with food waste (A S , A L , and B S ) suggesting that human use is the origin of the pharmaceutical. For ipronidazole, there was no toxicological data available, and the risk quotient was not calculated.

Steroids
The glucocorticoid hormone prednisolone was found in one digestate sample only, while dexamethasone was found in trace amounts in almost all digestates. Presently, no toxicological information on prednisolone Table 4 Predicted environmental concentration (PEC), predicted no-effect concentration (PNEC) for soil organisms, and the corresponding risk quotients (RQ). ARB, antibiotic-resistant bacteria. PNEC values are based on Table S11, S12, and S15 in the Supplementary information 1 Estimated from a simulated Koc value and from aquatic toxicity data using the partitioning coefficient method 2 Derived from PNEC values in soil pore water estimated by Menz et al. (2019)  is reported for soil organisms, but studies on aquatic organisms have confirmed endocrine disrupting effects (Bal et al. 2017). Based on available toxicity data and K d calculated from a modelled K OC (see details in Supplementary information), a PNEC SOIL of 0.09 μg kg −1 was calculated for prednisolone, yielding a risk quotient of 2. This number is only indicative but confirms the need of conducting in-depth toxicity tests for glucocorticoids and steroids on soil organisms as well as aquatic organisms.
Dexamethasone has received increased attention as the World Health Organization is recommending dexamethasone to treat severe to critical COVID-19 cases (WHO 2020). This synthetic steroid is more effective compared to natural steroids but also more persistent in the environment due to the fluorine moiety in the molecule. As dexamethasone only was found in trace amounts, it was not done a risk evaluation of this compound.

Concluding remarks
Food waste is an understudied matrix with regard to pharmaceutical pollution, as manure and sewage sludge are assumed to be more important entry routes to the soil. However, our results indicate that food waste can be an important entry route to the environment as well. Two other studies on pollutants in digestate from a variety of feedstocks (including food waste) also failed to identify a clear relationship between biogas feedstock and the level of pharmaceuticals and personal care products (PPCP) (Ali et al. 2019, using the same sample set as here) and POPs (Suominen et al. 2014). As confirmed by our results, significant levels of pharmaceutical residues were detected in Norwegian biogas digestates despite well-established retainment technologies and national regulations for pharmaceuticals in veterinary and human medicine. Hence, it is likely that pharmaceutical residues are common in many biogas digestates in other countries as well. If untreated, the presence of antibiotics in biogas digestates may lead to increased antibiotic resistance, harm towards soil organisms, leaching to water bodies, and potential exposure of human consumers, when digestates are applied on agricultural land.
Author contribution Astrid S. Nesse collected digestate, prepared the original draft, performed data analysis, prepared the graphical abstract, and edited the Supplementary information. Stine G. Aanrud developed and adjusted the method for analysis of antibiotics and hormones in biogas digestate, wrote the original draft for part S1 to S8 in the Supplementary information, and edited the manuscript. Jan L. Lyche advised on how to evaluate the ecotoxicity data and edited the manuscript. Trine A. Sogn contributed to the study conception and design and edited the manuscript. Roland Kallenborn was project leader, quality assured the method for analysis of antibiotics and hormones, and edited the manuscript. All authors read and approved the final manuscript.
Funding Open access funding provided by Norwegian University of Life Sciences The research received funding by the Research Council of Norway (RCN) Project 268214/E50: Novel organic pollutants from recycling of organic waste as risk factors for human exposure (NovelPol).
Data availability All the research data used in the study are available in Table 1.

Declarations
Ethics approval The research did not involve human or animal participants, and there was no release of harmful substance to the environment as a result of the study. The authors followed the rules for good scientific practice, as described in the author guidelines.
Consent to participate Not applicable, as there were no human participants in the study.

Consent for publication
All authors gave their explicit consent to publish the manuscript before it was uploaded to ESPR.

Conflict of interest The authors declare no competing interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.