A Method Scope Extension for the Simultaneous Analysis of POPs, Current-Use and Banned Pesticides, Rodenticides, and Pharmaceuticals in Liver. Application to Food Safety and Biomonitoring

The screening of hundreds of substances belonging to multiple chemical classes in liver is required in areas such as food safety or biomonitoring. We adapted a previous QuEChERS-based method in blood to the liver matrix and applied to these fields of study. The validation of the method allowed the inclusion of 351 contaminants, 80% with a LOQ < 2 ng/g. In the analysis of 42 consumer liver samples, we detected trace levels of 29 different contaminants. The most frequent and concentrated was 4,4’-DDE. POPs accounted for 66% of the compounds detected. In no case was the MRL reached for any of the contaminants detected. We also applied the method to 151 livers of wild birds to perform a biomonitoring pilot study in the Canary Islands. We detected 52 contaminants in 15 bird species. These were also mostly POPs, although high frequencies and concentrations of anticoagulant rodenticides (AR) and some other agricultural pesticides also stand out. POPs and AR contamination levels were significantly higher in terrestrial birds, raptors and particularly in nocturnal birds. Pesticide contamination levels were also higher in terrestrial birds, as well as in non-raptors and diurnal birds. The validated method is simple, robust, and sensitive and performs well in a variety of practical scenarios, where it can be carried out relatively quickly and inexpensively.


Introduction
Animal liver is a common matrix for the search of chemical residues for several reasons. On the one hand, the use of veterinary drugs, which is widely accepted in veterinary practice to treat bacterial infections, parasitism (internal or external), inflammation, and other animal diseases or their symptoms in livestock practice [1,2], may condition the appearance of residues of these chemicals in foods of animal origin, especially in the liver [3][4][5]. This is particularly relevant, since the presence of antimicrobial compounds can induce the spread of drug-resistant pathogenic bacterial strains or produce allergic reactions in humans [3,6]. In addition, pesticide products used in agriculture can leave The main objective of this research is to demonstrate the applicability of the validated method in the two fields of application mentioned above: (a) verification of compliance with maximum residue limits in livers intended for human consumption; and (b) biomonitoring of contaminants in wildlife. For this purpose, a set of samples was collected for each of the two independent studies. The samples are described in the following subsections.

Sampling for the Food Safety Study
To verify the applicability of this method for the control of residues subject to MRLs in livers intended for human consumption, 46 liver samples from butcheries, supermarkets, and the general slaughterhouse of Gran Canaria were acquired: 34 samples of beef liver and 12 samples of chicken liver. All the samples, as they were acquired, were transferred to the laboratory and frozen at −20 • C until they were processed.

Sampling for the Biomonitoring Study
The validated method was applied to real samples of wildlife specimens that were received in our laboratory for forensic analyses in the period between September 2020 and Toxics 2021, 9, 238 6 of 62 May 2021. Thus, we studied a series of 151 fresh liver samples belonging to 15 different species of birds. All the specimens were sent by environmental agents or by the Tafira Fauna Rehabilitation Centre, within the framework of the Strategy for the Prevention and Control of Poisoning in the Canary Islands [36]. All the birds included in this part of the study died from different classes of trauma, and there was no suspicion that they died of poisoning. The species included in this study were: Accipiter nisus (n = 5); Ardea cinerea (n = 12); Asio otus canariensis (n = 34); Burrhinus oedicnemus (n = 10); Buteo buteo insularum (n = 12); Calonectris diomedea (n = 8); Ciconia ciconia (n = 2); Corvus corax canariensis (n = 16); Egretta garzetta (n = 4); Falco eleanorae (n = 2); Falco pelegrinoides (n = 6); Falco tinnunculus canariensis (n = 14); Larus michaellis (n = 14); Turdus merula (n = 4); and Tyto alba (n = 8). The livers, received or extracted at in situ necropsy, were kept at −24ºC until the moment of their processing for analysis. No animals were sacrificed for the purposes of this work.

Statistical Analyses
All statistical analyses were performed with GraphPad Prism v9.2 software (Graph-Pad Software, CA, USA). The distribution of the variables included in this study was evaluated using the Kolmogorov-Smirnov test. The concentration of most of the contaminants detected did not follow a normal distribution, so the results are expressed in terms of median and range. For this same reason nonparametric tests to check for statistical differences between groups were employed, as these evaluate the median rather than the mean, which is appropriate given the relatively high number of undetected values in some groups. Homogeneity of variance (homoscedasticity) was previously tested using Levene's test. The Kruskal-Wallis and Mann-Whitney U tests were used as nonparametric tests for overall and pairwise comparisons, respectively. However, as an additional check, pairwise comparisons were also performed using Student's t-test after logarithmic transformation of the data. A P-value of less than 0.05 (two-tailed) was considered statistically significant. The prevalence of exposure to each contaminant for each species was calculated as the percentage of animals with that residue detected in the liver over the total number of individuals of that species studied. For the study of determinants in the series used for biomonitoring, the response variables considered for comparisons were the amount in the liver of (a) the sum of non-persistent pesticides; (b) the sum of persistent organic pollutants; and (c) the sum of rodenticides.

Method Scope Extension Optimization
In our previous research we optimized and validated a multi-residue method for the analysis of 360 substances (pharmaceuticals, pesticides, rodenticides, and POPs) in blood for biomonitoring purposes [16,28]. Therefore, this is not an ex-novo methodological development, but an extension of the scope of our previously published method to include a new matrix, the liver. However, for a better method performance in this more complex matrix, we considered optimizing the previously established chromatographic conditions, including recalculation of RTs, as well as optimization of MRM transitions to allow for higher sensitivity, as well as adjusting qualifiers and qualifier ratios, and identifying possible interferences with matrix components. The compounds are shown in alphabetical order in Appendix A along with their retention time, transitions, and their collision energies. As we did with the original method in blood, we decided to directly inject the extracts obtained in acetonitrile for LC-MS/MS and GC-MS/MS analyses, without using evaporation and solvent change, to avoid the loss of the more volatile compounds. Several authors, including our group [30,37,38], have shown that ACN, although not the most commonly used solvent in GC-MS/MS, is an appropriate solvent for this type of analysis.
The final number of validated compounds in this scope extension counts 351 chemicals and metabolites compared to 360 in the previous work. With respect to the original method, there are 18 compounds that met the validation criteria in whole blood, which do not meet the validation criteria when the method is applied to liver samples: acetaminophen, chlor-Toxics 2021, 9,238 7 of 62 fenapyr, corticosterone 21 acetate, phenbutatin oxide, iprodione, isocarbophos, leptophos, malaoxon, malathion, marbofloxacin, methomyl oxime, N,N,-dimethyl-N-tolylsulfamide, paraoxon ethyl, parathion ethyl, penicillin G, phosmet oxon, piperacillin, and trichlorfon. On the other hand, the opposite occurred with 9 compounds. Dichlorvos, doramectin, metalaxyl, methiocarb-sulfoxide, moxidectin, oxime, pthalimide, pyrimicarb-desmethyl and spirotetramat met the validation criteria in the presence of liver matrix and could therefore be included in the method in liver, whereas in blood they did not and had to be left out.

Validation Parameters
For confirmation of compound identity and selectivity, 2 MRM transitions were used, one for quantification (Q) and one for confirmation (q). A maximum deviation of ±30% was tolerated for the ion ratio. Similarly, a maximum deviation of ±0.1 min was established for the retention time.
We first studied the linearity of the response by injecting blank liver extract samples spiked with all analytes at 11 levels (range 0.4-40 ng/g) and processed in quintuplicate with the method described in Section 2.3. The linearity study on the response (R 2 ), indicated that this was higher than 0.98 for all analytes in the range studied.
To determine accuracy and precision, % recovery and % relative standard deviation (RSD) was calculated. A recovery within the range 70-120% and RSD values ≤ 20% was considered acceptable, as specified in the guidelines employed [33,34]. Recovery and RSD experiments were performed with blank liver samples fortified at least in four quintuplicate concentrations within the working range. The results of the recovery experiments are presented in Appendix B. Regarding accuracy and precision, most compounds meet the validation criteria for concentrations between their LOQ and the highest level studied (40 ng g −1 ). There were some exceptions where recoveries were outside the above range, especially at the lower concentrations. However, these cases are covered, both in the SANTE guideline and in the SWGTOX working document [33,34], which also accepts as a good validation criterion obtaining recoveries between 60% and 140% at some of the concentrations tested, provided that the RSD is less than 15%. Likewise, in some cases, the recoveries were within the established limits with an RSD slightly higher than 15%, a scenario that is also contemplated in the methodological guidelines, provided that the result is reproducible. As a rule, this second exception applies for concentrations equal to or lower than 4 ng g −1 . As SANTE analytical guide recommends, the expanded measurement uncertainty (U') was calculated, from precision and bias, and all analytes presented U' < 50%, that complies with the requirement.
For the calculation of the LOQ, matrix-matched calibration curves were prepared in quintuplicate (0.2-20 ng g −1 ). From these, the lowest concentration level of each analyte that met the criteria for accuracy and precision was considered as the LOQ. As in the original method, the LOQ for the analytes included in this scope extension was calculated from five replicates of fortified blank matrix, within the working range. The lowest non-zero calibrator approximation was used to calculate LOQs. This means that the lowest point on the calibration curve that met the identity, bias, and precision criteria was established as the LOQ for a given compound. The LOQs for the 351 liver analytes are shown in Appendix B. The LOQ was set at 0.4 ng g −1 for 61 compounds, at 0.8 ng g −1 for 82 compounds, at 1.2 ng g −1 for 40 compounds, at 1.6 ng g −1 for 37 compounds, at 2 ng g −1 for 50 compounds, at 4 ng g −1 for 46 compounds, at 8.0 ng g −1 for 24 compounds, at 12 ng g −1 for 5 compounds, at 16 ng g −1 for 4 compounds, and at 20 ng g −1 for 2 compounds. That is, 76.9% of the compounds included in this method can be reliably and accurately quantified at concentrations below 2 ng g −1 , making it suitable not only for food safety or poisoning diagnostic studies, but also for biomonitoring studies.
In the original method from which we started it was observed that there was a strong blood matrix effect on about 40% of the analytes. Presumably, a similar situation would occur with the liver matrix. Nevertheless, we decided to include the study of the matrix effect within the validation strategy of this analytical scope extension to prove it, as recommended in the reference guides. All validation assays involve the addition of known concentrations of analytes to the matrix. For the matrix effect study, we worked with the addition of three known concentrations of all analytes (2 ng g −1 , 10 ng g −1 , and 20 ng g −1 ) on blank liver extract, and the quantification was done against calibration curves prepared in solvent (without matrix). Experiments were performed in quintuplicate for each concentration. One difficulty was that, given the enormous number of substances included in the method, the liver was not completely free of 100% of the chemicals, in particular POPs. Therefore, in these cases, the response of the white matrix sample was subtracted from the calibration standards and QC to calculate the response of the externally added analyte. As we expected, matrix effect (ME) was observed for both, compounds analyzed by LC-MS/MS and GC-MS/MS, especially for compounds analyzed by the latter technique. A strong or medium suppression of the signal was demonstrated for 17.66% of the compounds (n = 62), and signal enhancement was verified for 36.47% of the compounds (n = 128). For the remaining 45.87% (n = 161 contaminants, the ME was considered negligible (−20% < M < 20%). Since for most of the compounds, significant ME was indeed observed, and it was concluded that matrix-matched calibration had to be used to compensate for these interferences. All detailed ME data for individual compounds in liver are shown in Appendix C.
Finally, we also assessed if carryover occurred after injecting a blank matrix fortified at 80 ng g −1 and processed with this method, before a blank matrix extract. We were not able to find a clear response in that blank matrix, so we concluded that in our working range, we had not any carryover effect in any of the analyzed compounds.

Application to Food Safety
In the study of the 34 beef liver samples, the results indicated the presence of a discrete number and concentration of contaminants, which ranged from 0 to 15 residues per sample, with an average of 3.13 residues. Of the 351 contaminants and metabolites included in the method, only 25 were detected in the total of beef liver samples, and of these 19 belong to the group of persistent or semi-persistent contaminants (4,4'-DDE, 4,4'-DDD, Dieldrin, Hexachlorobenzene, beta-hexachlorocyclohexane, BDE 153, PCB congeners #105, 118,126, 138, 155, 156, 157, 180, 189, naphthalene, phenanthrene, and pyrene). It is noteworthy that none of the concentrations in any of the samples exceeded the MRL, or even the value of half the MRL. In general, the concentrations of the contaminants detected were low, with the highest values being those of 4,4'-DDE, which was detected in 65.2% of the samples and with a median value of 92.2 ng g −1 . The relatively high levels of DDT derivatives may seem surprising, as this substance was banned in Spain almost 5 decades ago. However, there is abundant literature that has documented that this pesticide was widely used in the Canary archipelago, and how this translates into the levels of this pesticide detected in food for human consumption produced in this region [39][40][41][42].
The next in frequency and concentration were PCB 153 (26.1%; 35.2 ng g −1 ) and PCB 138 (21.8%; 24.5 ng g −1 ). The other contaminants were detected in frequencies and concentrations much lower than these. Among the non-persistent pesticides detected in this series of consumption livers, very low levels of bifenthrin, fenazaquin, fluquinconazole, flutalonil, flutriafol, and imidacloprid were detected.
If the detection of residues in beef liver was low and of little toxicological relevance, it was even more so in chicken liver. In the 12 samples analyzed, we detected only four contaminants out of the 351 included in the method: fenpropidin, fenpropimorph, levamisole, and 4,4' DDE. The latter was the more relevant, and it was only detected in three of the 12 livers analyzed and at a much lower concentration than that detected in beef liver (mean = 4.3 ng g −1 ).
Although it is not the main objective of this study, we made an estimate of the risk of exposure to these contaminants through liver consumption. The calculations were made according to the standard methodology that has been described previously [43], and in no case were the tolerable daily intake levels for these contaminants exceeded, mainly due to the low consumption of liver by the Spanish population (only 1 g/day for the total offal consumption) [44].

Application to Biomonitoring
Regarding biomonitoring of chemical substances, this method was applied to fresh livers obtained from 151 carcasses of 15 species of wild birds whose causes of death were not related to poisoning (mainly trauma). Table 1 shows the results for each of the species, limited to show only the 52 contaminants that were detected in the series. This represents that 15% of the contaminants included in the method were detected.
The mean value of the number of contaminants per sample was 17. The species with the greatest variety of residues detected was Asio otus (n = 41), followed by Falco tinnunculus (n = 27). In contrast, the species with the lowest number of liver contaminants were Turdus merula (n = 5) and Ciconia ciconia (n = 3). Figure 1 shows the LC-MS/MS and GC-MS/MS chromatograms of one of the birds in the series with the highest number of different contaminants (a long-eared owl).
The most frequently detected contaminant was 4,4'-DDE, which was detected in 138 birds (91.4%), followed by PCB 153, detected in 116 animals (76.8%), brodifacoum in 109 animals (72.2%), bromadiolone in 87 animals (57.6%), and dieldrin in 59 animals (39.1%). With respect to concentrations, the highest concentrations corresponded to enrofloxacin, clindamycin and meloxicam (Table 1). However, these values cannot be considered within the biomonitoring study, since they correspond to drugs used during the hospitalization of many of these animals. Therefore, high concentrations of these substances have been marked with an asterisk. However, other veterinary pharmaceuticals detected in some specimens, such as tetraconazole, metronidazole, or sulfathiazole, are not part of the treatment administered and should be considered contaminants. In general terms, the highest concentrations of contaminants corresponded to 4,4'-DDE in all species. Overall, in quantitative terms, the group of organochlorine pesticides was the most abundant ( Figure 2), and the group of persistent and semi-persistent organic pollutants accounted for more than 92% of the total concentration of pollutants detected in the livers of wild birds sampled in the Canary Islands very recently (September 2020-May 2021). This reflects, once again, that contamination by these compounds, in particular organochlorine insecticides, is still very prevalent in the Canary Islands, as has been reported for wildlife in this region [28,[45][46][47]. As indicated in the previous section, there is a large literature body documenting the high levels of contamination by organochlorine pesticides in this region [48][49][50], which also translates into high levels in the biota that inhabit the archipelago. There is a possibility that the high levels detected could also come from the neighboring African continent [51], but in this biomonitoring study this option is ruled out, since all the birds sampled for this pilot study are residents in the archipelago and not migratory birds.
With respect to non-persistent pollutants, several aspects should be highlighted. First, the high prevalence of second-generation anticoagulant rodenticides in wildlife's liver is noteworthy. It was expected, as it has been described in many parts of the world [52][53][54] and recently in the Canary Islands [45,55,56]. However, the presence of at least one of these compounds in more than 80% of the birds studied is striking, even in non-predatory birds such as the blackbird (Turdus merula) or the common curlew (Burhinus oedicnemus), which would point to the fact that these compounds penetrate the trophic chain by several routes, probably including invertebrates, as suggested by other authors [57,58]. Fludioxonil 16.4 (10) -- 0.4 (14) -- 0.4 (14) 0.6 (14) 7.5 (100) 0.4 (13) Dichlorodiphen yldichloroethane (p,p DDD) 3.6 (20) 3.1   The result for carbofuran is also surprising, given that none of the birds studied had any suspicion of intoxication. However, this potent insecticide, banned in the EU since 2007 [59] was detected in small concentrations in the liver of 10 birds of the series, being higher in the case of canary crows (6/12 positives, median = 94.5 ng/g). In all these crows, the main carbofuran metabolite was also detected. This toxicant has widely affected wildlife worldwide [60][61][62][63] and in the Canary Islands its use has been extensive and also has affected wildlife in the past [64]. From the results of this study, it still is today, and it can be concluded that it even penetrates the trophic chain. With respect to the rest of the non-persistent compounds detected in this series, 2-phenylphenol (PHP) stands out. PHP was detected in eleven birds, including five common curlews (Table 1). PHP is a biocide used as a preservative and surface disinfectant on fibers and other materials in homes, hospitals, and elsewhere, and is recognized as a potential endocrine disruptor [65]. Other authors have also reported that PHP is a highly prevalent contaminant in biota samples, such as river fish of different species, where it is found in up to 100% of samples [66]. pled in the Canary Islands very recently (September 2020-May 2021). This reflects, once again, that contamination by these compounds, in particular organochlorine insecticides, is still very prevalent in the Canary Islands, as has been reported for wildlife in this region [28,[45][46][47]. As indicated in the previous section, there is a large literature body documenting the high levels of contamination by organochlorine pesticides in this region [48][49][50], which also translates into high levels in the biota that inhabit the archipelago. There is a possibility that the high levels detected could also come from the neighboring African continent [51], but in this biomonitoring study this option is ruled out, since all the birds sampled for this pilot study are residents in the archipelago and not migratory birds. With respect to non-persistent pollutants, several aspects should be highlighted. First, the high prevalence of second-generation anticoagulant rodenticides in wildlife's liver is noteworthy. It was expected, as it has been described in many parts of the world [52][53][54] and recently in the Canary Islands [45,55,56]. However, the presence of at least one of these compounds in more than 80% of the birds studied is striking, even in nonpredatory birds such as the blackbird (Turdus merula) or the common curlew (Burhinus oedicnemus), which would point to the fact that these compounds penetrate the trophic chain by several routes, probably including invertebrates, as suggested by other authors [57,58]. Since this was an opportunistic study on carcasses obtained from wildlife recovery centers, we did not have too many quality variables to carry out an in-depth study of the determinants of contamination patterns. Even so, we wanted to explore the influence of the variables inherent to the species studied and found a series of statistically significant differences. Thus, when we compared aquatic versus terrestrial birds, we found that the latter presented significantly higher levels of contamination by the three major chemical groups studied ( Figure 3). The result for carbofuran is also surprising, given that none of the birds studied had any suspicion of intoxication. However, this potent insecticide, banned in the EU since 2007 [59] was detected in small concentrations in the liver of 10 birds of the series, being higher in the case of canary crows (6/12 positives, median = 94.5 ng/g). In all these crows, the main carbofuran metabolite was also detected. This toxicant has widely affected wildlife worldwide [60][61][62][63] and in the Canary Islands its use has been extensive and also has affected wildlife in the past [64]. From the results of this study, it still is today, and it can be concluded that it even penetrates the trophic chain. With respect to the rest of the nonpersistent compounds detected in this series, 2-phenylphenol (PHP) stands out. PHP was detected in eleven birds, including five common curlews (Table 1). PHP is a biocide used as a preservative and surface disinfectant on fibers and other materials in homes, hospitals, and elsewhere, and is recognized as a potential endocrine disruptor [65]. Other authors have also reported that PHP is a highly prevalent contaminant in biota samples, such as river fish of different species, where it is found in up to 100% of samples [66].
Since this was an opportunistic study on carcasses obtained from wildlife recovery centers, we did not have too many quality variables to carry out an in-depth study of the determinants of contamination patterns. Even so, we wanted to explore the influence of the variables inherent to the species studied and found a series of statistically significant differences. Thus, when we compared aquatic versus terrestrial birds, we found that the latter presented significantly higher levels of contamination by the three major chemical groups studied ( Figure 3). This result was expected with regard to rodenticides, since in a previous study by our group focused on these compounds, we had already discarded the group of waterfowl due to their low incidence in these pollutants [56]. Regarding POPs and non-persistent This result was expected with regard to rodenticides, since in a previous study by our group focused on these compounds, we had already discarded the group of waterfowl due to their low incidence in these pollutants [56]. Regarding POPs and non-persistent pesticides, although there is not much literature comparing both types of birds from the same region, the available studies usually indicate results similar to ours, with levels in landbirds usually being higher than in waterbirds [67][68][69].
Another variable that seems to influence the pattern of contamination is the raptor/predator bird status. The raptors in our study presented higher levels of POPs than non-predatory birds (Figure 4), which is logical given that they feed higher in the trophic chain, and has been described in the literature [68,69]. They also presented higher levels of AR, as we expected from having previously observed it in this region [56], and also described by other authors [70]. However, in the case of agricultural pesticides the statistical significance was the opposite, with non-predatory birds presenting the highest levels. There is not much literature to support this finding, but a recent study using the terrestrial pesticide residue exposure (T-REX) model estimated that the highest risk was presented by insectivorous birds, followed by fruit and seed feeders [71].
Toxics 2021, 9, x FOR PEER REVIEW 4 of 52 non-predatory birds (Figure 4), which is logical given that they feed higher in the trophic chain, and has been described in the literature [68,69]. They also presented higher levels of AR, as we expected from having previously observed it in this region [56], and also described by other authors [70]. However, in the case of agricultural pesticides the statistical significance was the opposite, with non-predatory birds presenting the highest levels.
There is not much literature to support this finding, but a recent study using the terrestrial pesticide residue exposure (T-REX) model estimated that the highest risk was presented by insectivorous birds, followed by fruit and seed feeders [71]. Finally, we also studied the influence of the diurnal/nocturnal habits of the birds in the study, and found that diurnal species have higher pesticide levels, but lower POPs and ARs than nocturnal species ( Figure 5). We believe that the pesticide result has to do with the previous variable, in the sense that, in our study, all insectivorous species, and those that feed on fruits and seeds are diurnal, while the nocturnal birds in our series are both raptors that feed mainly on large and small rodents. For this same reason, and as we had already verified in previous studies, the nocturnal birds of the Canary Islands have higher levels of POPs [72] and AR [45,56]. Finally, we also studied the influence of the diurnal/nocturnal habits of the birds in the study, and found that diurnal species have higher pesticide levels, but lower POPs and ARs than nocturnal species ( Figure 5). We believe that the pesticide result has to do with the previous variable, in the sense that, in our study, all insectivorous species, and those that feed on fruits and seeds are diurnal, while the nocturnal birds in our series are both raptors that feed mainly on large and small rodents. For this same reason, and as we had already verified in previous studies, the nocturnal birds of the Canary Islands have higher levels of POPs [72] and AR [45,56]. Toxics 2021, 9, x FOR PEER REVIEW 5 of 52

Conclusions
The validated method allows the simultaneous analysis in liver of 351 substances (POPs, pesticides including rodenticides and drugs), using only 1 gram of sample. This is important, since in veterinary forensic medicine, especially with small animals, the amount of sample available is very limited. The proposed analytical method can detect trace amounts of all chemicals in the liver of multiple species. Therefore, it can be successfully applied and used as a routine method in environmental chemistry and forensic toxicology laboratories. The method we have developed can also be used in residue control studies in food intended for human consumption and for the purpose of food safety assessment.

Conclusions
The validated method allows the simultaneous analysis in liver of 351 substances (POPs, pesticides including rodenticides and drugs), using only 1 gram of sample. This is important, since in veterinary forensic medicine, especially with small animals, the amount of sample available is very limited. The proposed analytical method can detect trace amounts of all chemicals in the liver of multiple species. Therefore, it can be successfully applied and used as a routine method in environmental chemistry and forensic toxicology laboratories. The method we have developed can also be used in residue control studies in food intended for human consumption and for the purpose of food safety assessment.

Conflicts of Interest:
The authors declare no conflict of interest. Table A1. Chromatographic and mass spectrometric conditions of the compounds analyzed in liver.

No.
Compound