Nadeem Muhammadab,
Fenglian Wangb,
Qamar Subhanibc,
Qiming Zhaobe,
Muhammad Abdul Qadird,
Hairong Cuia and
Yan Zhu*b
aDepartment of Environmental Engineering, Wuchang University of Technology, Wuhan, 430223, China
bDepartment of Chemistry, Xixi Campus, Zhejiang University, Hangzhou 310028, China. E-mail: zhuyan@zju.edu.cn; Fax: +86 571 88823446; Tel: +86 571 88273637
cHigher Education Department, Punjab, Lahore, Pakistan
dInstitute of Chemistry, University of the Punjab, Lahore, Pakistan
eCollege of Pharmacy, Zhejiang Chinese Medical University, Hangzhou 310053, China
First published on 2nd March 2018
There are increasing concerns about the dietary risks of neonicotinoids (NNIs); therefore their sensitive and accurate determination in dietary products is indispensable. However, the complex composition of agricultural food matrixes makes their extraction and quantitative determination a challenging task. Realizing this need, we herein report a simple, cost-effective, selective and sensitive fluorescence analytical workflow for analyses of two non-fluorescent neonicotinoids imidacloprid (IMI) and clothianidin (CLT) in six complex food samples (honey, ginger, durian, apple, tomato, cucumber) by online clean-up of sample extracts using two-dimensional ion chromatography (2D-IC) and a subsequent online post column UV induced fluorescence detection system. This online clean-up setup has proven advantageous to improve the limit of detection, potentially diminish matrix effects, and reduce analysis time and labor. The developed method showed excellent analytical figures of merit including linearity, selectivity, repeatability, recovery, and resolution for analysis of IMI and CLT in food samples.
On the other hand, there is rising public concern about the dietary risks of NNIs. A number of research studies in the USA and Europe have manifested that even a trace amount of neonicotinoids may cause disorientation, reduce longevity, impair memory and learning, decrease communication and can disrupt the brood cycles of honey bees.8 Especially, since last decade there is a dramatic weakness of honeybee hives (Apis mellifera L., Hymenoptera: Apidae) and a decrease in production of honey observed due to their contact with the honey bees which lead to their depletion. There is a maximum probability of the existence of their trace residue in honey; and decline of honey production is a threat to global food and medicine due to the growing importance of honey.9,10 The European Commission (EU) has banned the use of imidacloprid, clothianidin and thiamethoxam in crops for two years in order to accentuate the awareness of the potentially harmful effect of these pesticides on food products and in honey.11–13 Later, to preserve global health and food, the amended European Union legislation has set maximum residue limits (MRLs) for neonicotinoid insecticides in a number of agricultural products between 0.1 and 1 mg kg−1 (imidacloprid 0.05 mg kg−1 and clothianidin, 0.01 mg kg−1, respectively).2,14–16
Therefore, to the residual determination of NNIs in agricultural products for better monitoring of food safety, quality and to minimize public health problems; various sensitive, expensive instrumentation being used including gas chromatography coupled with mass spectrometry (GC-MS)17 and liquid chromatography braced with tandem mass spectrometry (LC-MS/MS), but due to thermolability, low volatility and high polarity of neonicotinoids residues later one is preferable technique.
However, NNIs polar nature, cost and matrix effect in complex samples make their access and use limited.18–21 Whereas, conventional analytical approaches, liquid chromatography coupled with a diode array detector (LC-DAD) and electrochemical detector are cost-effective and simple techniques but poor selectivity and sensitivity have discouraged their application for trace determination of NNIs in complex matrices.2,17,22–29 Recently, fluorescence spectrometry appeared as one of the simplest, cost-effective, rapid, sensitive, selective and reliable detection technique, but its application is limited by intrinsic non-fluorescent property of NNIs.30,31 This advantageous detection technique was reinstated by direct online photochemical conversion of non-fluorescent analytes into fluorescent species.32 The inline photo conversion of analytes is advantageous as it allows exploitation of incomplete reaction (photogenerated unstable radical), improve the efficiency of the photochemical process, minimum use of solvents, minimize generation of waste effluents and save derivatization labor and valuable time.33 However, neonicotinoid pesticides only exhibit photochemical fluorescence signal in basic pH, therefore, rare reverse phase chromatographic methods have been reported for determination of IMI only by post-column addition of 0.1 M NaOH.32 This extra post-column addition of pump not only make system complex but also increases the chances of baseline drifting, noise, poor repeatability and LOD. This problem to some extent our lab has overcome by introducing direct ion chromatographic based post-column photochemical induced fluorescence determination of pesticides in basic media.34,35 Despite all of it, HPLC/MS-MS is advantageous for meanwhile identification of target analyte and the qualifier/quantifier ions. Whereas, in the case of photoinduced fluorimetric determination method unable to identify the possible changes in the structure of target non-fluorescent analytes.
However, this and most of the reported chromatographic methods have been developed for the determination of individual IMD in simple samples (water and soil).5,23,30 Therefore, an effective sample preparation is vital not only for analytes pre-concentration but also to remove or minimize other compounds, impurities and matrix interferences without losing analytes of interest. Many off-line sample pretreatment methods have been used to overcome these problems to some extent including solid phase extraction (SPE),17,25,36 dispersive solid-phase extraction, dispersive liquid–liquid microextraction,17,26 diatomaceous earth-assisted extraction37 and ionic liquid phase microextraction.2 Most of the above stated sample preparation methods have some major limitations; huge consumption of chemicals and solvent, time-consuming, multiple operation steps required which make these expensive and laborious; whereas the inability to remove matrix interferences of complex samples is a serious limitation. Recently, to compensate the first problem to some extent, there is increasing trend of the development of miniaturized, more efficient, economical and green extraction methods that could significantly reduce the toxic organic solvents and chemical consumption.38–40 Particularly, QuEChERS sample preparation method got attentions for extraction of pesticides from complex food and other matrices because of its simplicity, miniaturization, inexpensiveness, efficiency, ruggedness and amenability to high throughput. However, it consists of multiple extraction and cleaning steps; despite that, it is unable to complete elimination of matrix interferences and risk of analyte contamination and loss are maximum during the multiple extraction, evaporation and reconstitution steps; especially, from complex food samples like honey, fruit and vegetables. The natural honey is a sticky and viscous solution having a complex composition with contents of carbohydrates, protein, ash, amino acids, vitamins, enzymes as well as phenolic antioxidants.10 Similarly, agricultural products (fruits and vegetables) have high molecular weight and complex interferences like lipids (waxes, triacylglycerols and phospholipids etc.), various pigments (carotenoids, chlorophylls, melanoidin etc.) and plant resins, which have adversely affected the analysis result by peak masking, false identification of impurity as the analyte peak, impaired detectability, deposition of salt on GC inlet, elution of volatile impurity along with analyte and signal enhancement and suppression. Besides these, QuEChERS method every time needs to modify according to the nature of matrix and analytes of interest.19,38,39,41–44 The matrix interferences and analyte loss problems can be resolved by making sample preparation method online with the 2D-IC system just like SPE methods after making online proved more efficient, versatile, reproducible, selective, reduce labor, human error, costs, solvent and time.45–48
Since last two decades, an online column-switching technique for heart-cutting is extensively utilized for complex samples clean-up and concentration to reduce matrix interferences.49–52 However, most of the reported online extraction techniques use the mini hydrophobic column as pre-treatment column which unable to complete elimination of matrix interferences from analytes of interest and matrices co-elute and transfer along the analytes of interest in second dimension column. This drawback is severe in the case of polar analytes.32,53 Whereas, full C18 column as a pretreatment column most of the time help to clean and transfer only single peak to the second dimension.50 Therefore, it is need of time to develop an efficient, cheap, selective and sensitive online analytical approach for the determination of more than one neonicotinoids at trace level in complex food samples to access pesticides exposure to human health.54 Recently, ion chromatography appeared an efficient and rapid alternative technique for separation and determination of polar analytes apart from ions.31,43,55 The electrostatic forces of attraction between polar NNIs and ion exchange columns can be exploited for matric inferences elimination in one dimension and separation of two similar non-fluorescent NNIs in second dimensions with minimum consumption of solvent and time. As far as we know, there are no published reports regarding direct simultaneous separation and determination of two non-fluorescent NNIs (IMI and CLT) by online post column photoinduced fluorescence detection in basic media in complex food samples.
Hence, in order to minimize analyte loss, multiple labor steps, matrix interferences, and to avoid the need of sorbent modification, hereby we report for the first time a systematic online clean-up setup by using 2D-IC column switching technique for the residue analysis of two neonicotinoid pesticides in six complex food samples followed by their sensitive and selective fluorescence determination by exploiting simple post column UV irradiation.
A home-made photo chemical reactor was designed by coiling the knitted PTFE tube (17 cm long × 1.5 mm o.d. × 0.5 mm i.d.) around a 20 cm long low-pressure mercury lamp (18 W, 254 nm). It was further fitted into a PVC cylinder and internally covered with aluminum foil to have back the maximum UV light reflection and to minimize light dissipation in surrounding environment. The data analysis was performed with a personal computer equipped with Chromeleon 7.2 software (Thermo-Fisher Scientific, Waltham, MA, USA). A precise pH-meter (Mettler Toledo, FE20 and China) was used for the pH measurements. The SB-5200DT ultrasonic cleaner (Scientz Biotechnology Co. Ltd., Ningbo, China) was used for sonication of sample and eluent. The agricultural samples ginger, durian, apple, tomato and cucumber were bought from a local market in Hangzhou (China) and stored at 4 °C. The third, multifloral pure honey sample was cultivated under our surveillance in Punjab, Pakistan, which was stored in an amber vial at 4 °C in darkness.
Firstly, the sample was manually loaded into the 25 μL sample loop of valve 1 (Fig. 2A), and it is carried to the pretreatment column Ionpac AS11-HC by pumping mobile phase A (55 mM NaOH + 20% ACN) by HPLC pump A. Most of the high molecular weight hydrophobic matrixes interferences were quickly discarded into waste, while both polar NNIs were retained into the AS11-HC column and after 4.50 min the value 2 having 1.1 mL loop turn into load position and start collecting target analytes within 1.0 min time span (Fig. 2B and C). After 1.0 min of analyte collection, the NNIs from the collection loop were delivered to the analytical column Ionpac AS12A protected by a guard column Ionpac AG12A with mobile phase B (1.5 mM Na2CO3 + 15 mM NaOH + 18% ACN) for their isocratic separation (Fig. 2D). Meanwhile, the pretreatment column was washed and regenerated by high strength mobile A. After 9.0 min the valve 1 switch back to load position and valve 2 switches to inject position to bring the 2D-IC system come back to the initial state (Fig. 2A) and to achieve the equilibrium for the subsequent sample analysis. The complete analytical procedure comprised of 31 min, in which 15.5 min for offline solid–liquid extraction of two analytes from complex samples and 15.5 min for their online clean-up and separation followed by photoinduced fluorescence detection. The complete operation procedure of 2D-IC setup is given in Table 1.
Position | Programming | Time (min) | Mobile phase of pre-treatment column | Mobile phase of analytical column | Valve 1 (V-1) | Valve 2 (V-2) |
---|---|---|---|---|---|---|
a Manual 25 μL sample injection.b Inject for 60 s. | ||||||
1 | Samplinga | 0.0 | 55 mM NaOH + 20% ACN | 1.5 mM Na2CO3 +15 mM NaOH + 18% ACN | Load | Inject |
2 | Sample extract pre-treatment | 0–4.50 | — | — | Inject | Inject |
3 | Analytes collection | 4.50–5.50 | — | — | Load | Load |
4 | Analyses, regeneration, and re-equilibrium | 5.50–15.50 | — | — | Load | Inject |
The linearity of the detector response for both NNIs was assessed under the optimal conditions from the standard and matrices matched calibrations standard at six concentration level in the range 12–5000 μg kg−1 for IMI and CLT. The selectivity of this method was observed by duplicate injection of blank samples extracts to observe the influence of matrix effect at a retention time of both analytes.56 The method limit of detections and quantification were estimated as the lowest concentration resulting in a signal-to-noise ratio equal to 3 and 10-fold, respectively. The LOQs were used as the lowest concentration at the calibration curve according to SANTE/11945/2015 guidelines.57 The precision of the method was determined at matrices matched three concentration levels LOQ, 2× LOQ and 10× LOQ and their relative standard deviations (%RSD) were computed as a measurement of the precision. In addition, due to the complexity of the matrices the intra- and inter-day precisions were also calculated in quintet within one day and over three days, respectively. The accuracy of the method was calculated as mean recoveries by spiking μ-QuEChERS extract at three level (LOQ, 2× LOQ and 10× LOQ) in the quintet.
All recoveries (%) were calculated by following equation.
Recovery (%) = (C1 – C2)/C3 × 100. |
Matrix effect (ME) was determined by comparing the slope ratios of the two sets; matrix matched (honey, durian, ginger, apple, tomato and cucumber) and pure solvent calibration curve obtained at six level of concentration for IMI and CLO. According to guidelines; matrix effect value is within the range of ±20, the results are acceptable.
ME (%) = (slope of calibration curve in matrix/slope of calibration curve in solvent − 1) × 100. |
There are various factors including nature of the basic medium, pH, solvent, and irradiation time were optimized in order to establish the optimum condition for the sensitive determination of NNIs by exploiting post column photo induced fluorescence determination in complex samples. All studies were performed by varying each respective variable in return keeping the others constant. Firstly, it is important to know optimum excitation and emission wavelength in order to have their maximum fluorescence intensity. Both analytes turn by turn scanned by passing online through the homemade photochemical reactor. It was observed IMI and CLT showed big peaks at λex/λem = 332/368 nm and λex/λem = 353/401 nm, respectively as shown in Fig. S1(a).†
A flow injection analysis was conducted to investigate the effect of basic media on the photoinduced fluoresce (PIF) intensity of NNIs. To observe this phenomenon various bases and their different IC column compatible combinations with each other (NaOH, KOH, Na2CO3, NaHCO3, Na2CO3 + NaHCO3, NaOH + NaHCO3 and NaOH + Na2CO3) were prepared with the constant addition of 18% ACN. It was observed IMI exhibited much high PIF intensity as compare to CLT in all bases, this probably due to the presence of the aromatic ring in IMI structure. The ionization of alpha proton in basic media extends the resonance path up to two rings which helped in to have their high PIF intensity. Both NNIs showed high PIF intensity in the strongly basic mobile phase of KOH and NaOH as compared to weakly basic Na2CO3 and NaHCO3. Whereas, CLT showed significantly higher PIF signal in the weakly basic medium in contrary to IMI. However, both CLT and IMI simultaneously showed high fluorescence intensity in a combination of NaOH and Na2CO3 or NaOH and NaHCO3 as depicts in Fig. S1(b).† Especially, CLT PIF signal increased two times in medium composed of strong and weak base (NaOH, Na2CO3) as compared to only strong base (NaOH, KOH). Therefore, weak and strong base combination (NaOH + Na2CO3) was used as a mobile phase for the separation and analysis of both NNIs in order to have their optimum PIF intensity.
The influence of pH on both NNIs PIF intensity was investigated by changing the composition of mobile phase with NaOH and Na2CO3 to have pH in the range of 7–14. A significant increase in PIF intensity of both NNIs was observed as pH increased from 6 to 14 and became maximum at pH 11–12. This effect was more dramatic in the case of IMI, while CLT PIF intensity slowly increased and became constant as shown in Fig. S2(a).† Therefore, the pH 11–12 was selected as optimum pH for their sensitive determination. This effect probably due to ease of NNIs degradation in basic media under UV light and changed into fluorescent species.
The PTFE capillaries have been used for the UV induce fluorescence. Therefore, the residence time of analytes of interest under UV light is very critical to have their optimum fluorescence intensity. The residence time of NNIs was controlled by variation of PTFE coil length and it observed IMI PIF intensity increases as its residence time elongated from 18 to 30 s after that it slowly decreased and became constant, while CLT displayed maximum PIF intensity for 18 s and then it remains almost constant as shown in Fig. S2(b)†.
Finally, IonPacs AS12A was successfully utilized for their isocratic separation by using only 18% ACN to improve peak shape and to finish separation within 10 min. Therefore, 1.5 mM Na2CO3 + 15 mM NaOH + 18% ACN was used as a mobile phase to gradual overcome all type of these interactions and to separate them isocratically as shown in Fig. 3.
Analyte | Retention time (min) | Peak asymmetry | Peak resolution | Linear range (μg kg−1) | Correlation (r2) |
---|---|---|---|---|---|
IMI | 4.90 | 1.32 | 4.29 | 0.12–1250 | 0.9999 |
CLT | 8.90 | 1.31 | — | 0.51–5000 | 0.9991 |
The LODs and LOQs of both NNIs in three samples were retrieved in the range from 0.035–0.154 μg kg−1 and 0.120–0.512 μg kg−1, respectively as given in Table 3. These low LODs are the result of high PIF intensity and maximum recoveries of analytes of interest due to the online clean-up of matrix interferences from sample extracts by using the 2D-IC system.
Analyte | Honey (μg kg−1) | Ginger (μg kg−1) | Durian (μg kg−1) | Apple (μg kg−1) | Tomato (μg kg−1) | Cucumber (μg kg−1) | ||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
LOD | LOQ | LOD | LOQ | LOD | LOQ | LOD | LOQ | LOD | LOQ | LOD | LOQ | |
IMI | 0.036 | 0.119 | 0.035 | 0.119 | 0.037 | 0.123 | 0.034 | 0.113 | 0.35 | 0.116 | 0.037 | 0.123 |
CLT | 0.153 | 0.510 | 0.151 | 0.502 | 0.154 | 0.512 | 0.142 | 0.472 | 0.143 | 0.476 | 0.144 | 0.479 |
The intra- and inter-day precisions of the method at three levels are expressed as the %RSD of the relative peak areas of the day controlled fortified samples. Both intra- and inter-day precisions values were within the range 0.12–12.69% and 2.16–12.36%, respectively as given in Table 4. The high precision of the method reflects the reproducibility of the method. The accuracy of the method was gauged in fortified blank samples at three concentration levels LOQ, 2× LOQ and 10× LOQ in the quintet and good recoveries were obtained in the range 82.0–106% as given in Table 5.
Analyte | Conc. level (μg kg−1) | Honey (% RSDs) | Ginger (% RSDs) | Durian (% RSDs) | Apple (% RSDs) | Tomato (% RSDs) | Cucumber (% RSDs) | ||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Intra | Inter | Intra | inter | Intra | Inter | Intra | Inter | Intra | Inter | Intra | Inter | ||
IMI | LOQ (0.12) | 0.12 | 2.60 | 5.18 | 8.96 | 1.16 | 8.68 | 1.25 | 3.68 | 5.6 | 7.52 | 0.78 | 3.07 |
2× LOQ (0.24) | 10.31 | 3.69 | 2.16 | 4.70 | 12.69 | 9.73 | 5.81 | 8.54 | 8.89 | 7.23 | 0.21 | 5.39 | |
10× LOQ (1.2) | 6.14 | 5.98 | 7.96 | 5.23 | 8.98 | 8.90 | 3.69 | 9.27 | 0.86 | 2.37 | 5.86 | 2.17 | |
CLT | LOQ (0.51) | 7.96 | 10.23 | 10.44 | 12.36 | 0.12 | 5.67 | 6.18 | 3.84 | 1.24 | 4.21 | 3.57 | 8.42 |
2× LOQ (1.02) | 9.38 | 11.40 | 4.16 | 7.23 | 3.19 | 7.65 | 7.23 | 4.87 | 5.98 | 8.71 | 2.79 | 8.12 | |
10× LOQ (5.1) | 0.68 | 9.51 | 3.24 | 5.16 | 5.26 | 11.89 | 8.95 | 4.20 | 0.95 | 3.58 | 0.91 | 4.81 |
Analytes | Con. (μg kg−1) | Honey | Ginger | Durian | Apple | Tomato | Cucumber |
---|---|---|---|---|---|---|---|
a Replicate of five readings. | |||||||
IMI | LOQ (0.12) | 77.41 ± 0.12 | 97.22 ± 5.18 | 93.23 ± 1.16 | 84.14 ± 7.13 | 101.24 ± 8.74 | 99.73 ± 3.16 |
2× LOQ (0.24) | 79.80 ± 10.31 | 100.46 ± 2.16 | 105.6 ± 12.69 | 86.10 ± 9.53 | 100.13 ± 7.12 | 101.2 ± 10.79 | |
10× LOQ (1.2) | 89.36 ± 6.14 | 84.50 ± 7.96 | 83.18 ± 8.98 | 82.43 ± 5.24 | 94.50 ± 7.06 | 93.12 ± 7.38 | |
CLT | LOQ (0.51) | 81.10 ± 7.96 | 83.18 ± 10.44 | 90.06 ± 0.12 | 81.80 ± 5.06 | 93.12 ± 11.04 | 98.03 ± 0.18 |
2× LOQ (1.02) | 95.23 ± 9.38 | 88.32 ± 4.16 | 106.8 ± 3.19 | 89.28 ± 10.18 | 88.72 ± 7.16 | 102.8 ± 8.13 | |
10× LOQ (5.1) | 90.75 ± 0.68 | 74.60 ± 3.24 | 73.71 ± 5.26 | 90.05 ± 1.73 | 84.80 ± 4.94 | 83.21 ± 3.29 |
One of the most important parameter of method validation and this method development is the matrix effect. It was the fundamental purpose of this method development to minimize the matrix effect in order to avoid the camouflage determination of NNIs in complex samples analysis. The blank samples extracts obtained by μ-QuEChERS method were fortified and matrix effect (ME) were calculated by using the matrix effect equation as detailed in Section 2.6. It was observed matrix effect was comprehensively reduced as all ME values are within the range of ±20% by integration of μ-QuEChERS method to the 2D-IC system as given in Table 6.
Analyte | Honey | Ginger | Durian | Apple | Tomato | Cucumber | ||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
ME | r2 | ME | r2 | ME | r2 | ME | r2 | ME | r2 | ME | r2 | |
a IMI: imidacloprid and CLT: clothadine. | ||||||||||||
IMI | −0.22 | 0.9993 | 10.02 | 1 | −14.21 | 0.9999 | 3.52 | 0.9980 | 5.78 | 0.9934 | −0.30 | 1 |
CLT | +11.90 | 0.9963 | +5.30 | 0.9990 | +11.06 | 1 | 1.23 | 0.9974 | 2.98 | 0.9945 | 0.40 | 0.9999 |
This developed method also displayed excellent selectivity after integration of μ-QuEChERS sample preparation method to the 2D-IC system as shown in Fig. 4a–c. The good sensitivity of the method can be observed from the all excellent determination of coefficient (r2) values ≥0.9999 and the low LOD values of both NNIs as given in Table 2. In addition, clean separation and analyses of both NNIs in complex samples are shown in Fig. 4 and comparison of this developed method with other is also given in Table 7.
Sample | Sample preparation method | No. of analytes | Solvent uses (mL) | Analysis time (min) | Detection | LOD | References |
---|---|---|---|---|---|---|---|
a DLLME: dispersive liquid–liquid microextraction; MEKC: micellar electrokinetic capillary chromatography. | |||||||
Honey | DLLME and QuEChERS | 7 | >2.5, >10 | 66 | HPLC-UV | 1.5–2.5 μg kg−1 | 54 |
Cucumber | QuEChERS method | 4 | >10 | 20 | HPLC-DAD | 0.04 mg kg−1 | 1 |
Eggplant, cucumber | Ionic liquid extraction | 7 | >33 | 30 | HPLC-DAD | 0.002–0.003 mg kg−1 | 25 |
Honey | Ionic liquid extraction | 4 | >50 μL | 14 | HPLC-DAD | 50 μg kg−1 | 58 |
Honey | SPE | 7 | >19 | No separation | LC/LC-MS | 0.03 μg kg−1 | 59 |
Potato | SPE | 4 | >80 | 9.10 | LC-DAD | 3.2–15 μg L−1 | 60 |
Cucumber soil | QuEChERS method | 3 | >10 | 14 | LC-DAD | 0.01–0.08 mg kg−1 | 61 |
Commercial formulation | Direct | 2 | N.A | N.A | Spectrophotometric | 0.17–0.32 mg L−1 | 62 |
Honey bee | SPE | 2 | >85 mL | 7 | HPLC-FLD | 0.09–1.5 μg kg−1 | 32 |
Honey bee, maize leave | Sample clean-up | 2 | >50–150 mL | 25 | LC-hν-ED | 2.4 mg L−1 | 53 |
SPE | 2 | 100 mL | 7 | MEKC-UV | 0.71–1.18 mg L−1 | 63 | |
Honey, ginger, durian | Online-sample clean-up | 2 | 900 μL | 9.5 | IC-hv-FLD | 0.036–0.15 μg kg−1 | This work |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c7ra12555k |
This journal is © The Royal Society of Chemistry 2018 |