Bioaccessibility of trace elements and Fe and Al endogenic nanoparticles in farmed insects: Pursuing quality sustainable food.

This study investigated the in vitro bioaccessibility of aluminum, copper, iron, manganese, lead, selenium, and zinc in three important species of farmed insects: the yellow mealworm (Tenebrio molitor), the house cricket (Acheta domesticus) and the migratory locust (Locusta migratoria). Results show that all three insect species constitute excellent sources of essential elements (Fe, Cu and Zn) for the human diet, contributing to the recommended dietary allowance, i.e., 10%, 50%, and 92%, respectively. A higher accumulation of Se (≥1.4 mg Se/kg) was observed with increasing exposure concentration in A. domesticus, showing the possibility of using insects as a supplements for this element. The presence of Al and Fe nanoparticles was confirmed in all three species using single particle-inductively coupled plasma-mass spectrometry and transmission electron microscopy. The results also indicate that Fe bioaccessibility declines with increasing Fe-nanoparticle concentration. These findings contribute to increase the nutritional and toxicological insights of farmed insects.

A comparison of screw and suture button fixation in the management of adolescent ankle syndesmotic injuries.

Purpose: Ankle injuries involving the tibiofibular syndesmosis often necessitate operative fixation to restore stability to the ankle. Recent literature in the adult population has suggested that suture button fixation may be superior to screw fixation. There is little evidence as to which construct is preferable in the pediatric and adolescent population. This study investigates outcomes of suture button and screw fixation in adolescent ankle syndesmotic injuries. Methods: A retrospective matched cohort study over 10 years of pediatric patients who underwent ankle syndesmotic fixation at a large Level 1 Trauma Center was conducted. Both isolated syndesmotic injuries and ankle fractures with syndesmotic disruption were included. Preoperative variables collected include basic patient demographics, body mass index, and fracture type. Suture button and screw cohorts were matched based on age, race, sex, and open fracture utilizing propensity scores. Outcomes assessed include reoperation and implant failure. Results: A total of 44 cases of operative fixation of the ankle syndesmosis were identified with a mean age of 16 years. After matching cohorts based on age, sex, race, and open fracture status, there were 17 patients in the suture button and screw cohorts, respectively. Patients undergoing screw fixation had a six times greater risk of reoperation (p = 0.043) and 13 times greater risk of implant failure (p < 0.001). Out of six cases of reoperation in the screw cohort, five were unplanned. Conclusion: Our findings favor suture button fixation in operative management of adolescent tibiofibular syndesmotic injuries. Compared with screws, suture buttons are associated with lower risk of both reoperation and implant failure. Level of evidence: level III therapeutic.

Peripheral blood mononuclear cell respiratory function is associated with progressive glaucomatous vision loss.

Intraocular pressure (IOP) is currently the only modifiable risk factor for glaucoma and all licensed treatments lower IOP. However, many patients continue to lose vision despite IOP-lowering treatment. Identifying biomarkers for progressive vision loss would have considerable clinical utility. We demonstrate that lower peripheral blood mononuclear cell (PBMC) oxygen consumption rate (OCR) is strongly associated with faster visual field (VF) progression in patients treated by lowering IOP (P < 0.001, 229 eyes of 139 participants), explaining 13% of variance in the rate of progression. In a separate reference cohort of untreated patients with glaucoma (213 eyes of 213 participants), IOP explained 16% of VF progression variance. OCR is lower in patients with glaucoma (n = 168) than in controls (n = 50; P < 0.001) and is lower in patients with low baseline IOP (n = 99) than those with high baseline IOP (n = 69; P < 0.01). PBMC nicotinamide adenine dinucleotide (NAD) levels are lower in patients with glaucoma (n = 29) compared to controls (n = 25; P < 0.001) and strongly associated with OCR (P < 0.001). Our results support PBMC OCR and NAD levels as new biomarkers for progressive glaucoma.

Estimating the Distribution of True Rates of Visual Field Progression in Glaucoma.

Purpose: The purpose of this study was to estimate the distribution of the true rates of progression (RoP) of visual field (VF) loss. Methods: We analyzed the progression of mean deviation over time in series of ≥ 10 tests from 3352 eyes (one per patient) from 5 glaucoma clinics, using a novel Bayesian hierarchical Linear Mixed Model (LMM); this modeled the random-effect distribution of RoPs as the sum of 2 independent processes following, respectively, a negative exponential distribution (the "true" distribution of RoPs) and a Gaussian distribution (the "noise"), resulting in a skewed exGaussian distribution. The exGaussian-LMM was compared to a standard Gaussian-LMM using the Watanabe-Akaike Information Criterion (WAIC). The random-effect distributions were compared to the empirical cumulative distribution function (eCDF) of linear regression RoPs using a Kolmogorov-Smirnov test. Results: The WAIC indicated a better fit with the exGaussian-LMM (estimate [standard error]: 192174.4 [721.2]) than with the Gaussian-LMM (192595 [697.4], with a difference of 157.2 [22.6]). There was a significant difference between the eCDF and the Gaussian-LMM distribution (P < 0.0001), but not with the exGaussian-LMM distribution (P = 0.108). The estimated mean (95% credible intervals, CIs) "true" RoP (-0.377, 95% CI = -0.396 to -0.359 dB/year) was more negative than the observed mean RoP (-0.283, 95% CI = -0.299 to -0.268 dB/year), indicating a bias likely due to learning in standard LMMs. Conclusions: The distribution of "true" RoPs can be estimated with an exGaussian-LMM, improving model accuracy. Translational Relevance: We used these results to develop a fast and accurate analytical approximation for sample-size calculations in clinical trials using standard LMMs, which was integrated in a freely available web application.

Correction: Won et al. Ex Vivo Perfusion Using a Mathematical Modeled, Controlled Gas Exchange Self-Contained Bioreactor Can Maintain a Mouse Kidney for Seven Days. Cells 2022, 11, 1822.

In the original publication [...].

Relationship between Intraocular Pressure Fluctuation and Visual Field Progression Rates in the United Kingdom Glaucoma Treatment Study.

PURPOSE: To investigate whether intraocular pressure (IOP) fluctuation is associated independently with the rate of visual field (VF) progression in the United Kingdom Glaucoma Treatment Study. DESIGN: Randomized, double-masked, placebo-controlled multicenter trial. PARTICIPANTS: Participants with ≥5 VFs (213 placebo, 217 treatment). METHODS: Associations between IOP metrics and VF progression rates (mean deviation [MD] and five fastest locations) were assessed with linear mixed models. Fluctuation variables were mean Pascal ocular pulse amplitude (OPA), standard deviation (SD) of diurnal Goldmann IOP (diurnal fluctuation), and SD of Goldmann IOP at all visits (long-term fluctuation). Fluctuation values were normalized for mean IOP to make them independent from the mean IOP. Correlated nonfluctuation IOP metrics (baseline, peak, mean, supine, and peak phasing IOP) were combined with principal component analysis, and principal component 1 (PC1) was included as a covariate. Interactions between covariates and time from baseline modeled the effect of the variables on VF rates. Analyses were conducted separately in the two treatment arms. MAIN OUTCOME MEASURES: Associations between IOP fluctuation metrics and rates of MD and the five fastest test locations. RESULTS: In the placebo arm, only PC1 was associated significantly with the MD rate (estimate, -0.19 dB/year [standard error (SE), 0.04 dB/year]; P < 0.001), whereas normalized IOP fluctuation metrics were not. No variable was associated significantly with MD rates in the treatment arm. For the fastest five locations in the placebo group, PC1 (estimate, -0.58 dB/year [SE, 0.16 dB/year]; P < 0.001), central corneal thickness (estimate, 0.26 dB/year [SE, 0.10 dB/year] for 10 µm thicker; P = 0.01) and normalized OPA (estimate, -3.50 dB/year [SE, 1.04 dB/year]; P = 0.001) were associated with rates of progression; normalized diurnal and long-term IOP fluctuations were not. In the treatment group, only PC1 (estimate, -0.27 dB/year [SE, 0.12 dB/year]; P = 0.028) was associated with the rates of progression. CONCLUSIONS: No evidence supports that either diurnal or long-term IOP fluctuation, as measured in clinical practice, are independent factors for glaucoma progression; other aspects of IOP, including mean IOP and peak IOP, may be more informative. Ocular pulse amplitude may be an independent factor for faster glaucoma progression. FINANCIAL DISCLOSURE(S): Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Evaluating Primary Treatment for People with Advanced Glaucoma: Five-Year Results of the Treatment of Advanced Glaucoma Study.

PURPOSE: To determine whether primary trabeculectomy or medical treatment produces better outcomes in terms of quality of life (QoL), clinical effectiveness, and safety in patients with advanced glaucoma. DESIGN: Multicenter randomized controlled trial. PARTICIPANTS: Between June 3, 2014, and May 31, 2017, 453 adults with newly diagnosed advanced open-angle glaucoma in at least 1 eye (Hodapp classification) were recruited from 27 secondary care glaucoma departments in the United Kingdom. Two hundred twenty-seven were allocated to trabeculectomy, and 226 were allocated medical management. METHODS: Participants were randomized on a 1:1 basis to have either mitomycin C-augmented trabeculectomy or escalating medical management with intraocular pressure (IOP)-reducing drops as the primary intervention and were followed up for 5 years. MAIN OUTCOME MEASURES: The primary outcome was vision-specific QoL measured with the 25-item Visual Function Questionnaire (VFQ-25) at 5 years. Secondary outcomes were general health status, glaucoma-related QoL, clinical effectiveness (IOP, visual field, and visual acuity), and safety. RESULTS: At 5 years, the mean ± standard deviation VFQ-25 scores in the trabeculectomy and medication arms were 83.3 ± 15.5 and 81.3 ± 17.5, respectively, and the mean difference was 1.01 (95% confidence interval [CI], -1.99 to 4.00; P = 0.51). The mean IOPs were 12.07 ± 5.18 mmHg and 14.76 ± 4.14 mmHg, respectively, and the mean difference was -2.56 (95% CI, -3.80 to -1.32; P < 0.001). Glaucoma severity measured with visual field mean deviation were -14.30 ± 7.14 dB and -16.74 ± 6.78 dB, respectively, with a mean difference of 1.87 (95% CI, 0.87-2.87 dB; P < 0.001). Safety events occurred in 115 (52.2%) of patients in the trabeculectomy arm and 124 (57.9%) of patients in the medication arm (relative risk, 0.92; 95% CI, 0.72-1.19; P = 0.54). Serious adverse events were rare. CONCLUSIONS: At 5 years, the Treatment of Advanced Glaucoma Study demonstrated that primary trabeculectomy surgery is more effective in lowering IOP and preventing disease progression than primary medical treatment in patients with advanced disease and has a similar safety profile. FINANCIAL DISCLOSURE(S): Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Stress fractures of the hand and wrist in athletes.

Stress fractures of the upper extremity are reported less often than their lower extremity counterpart. This review aims to provide a comprehensive overview of an important and often missed diagnosis in pediatric athletes: hand and wrist stress fractures.

European Glaucoma Society research priorities for glaucoma care.

BACKGROUND/AIMS: The goal of health research is to improve patients care and outcomes. Thus, it is essential that research addresses questions that are important to patients and clinicians. The aim of this study was to develop a list of priorities for glaucoma research involving stakeholders from different countries in Europe. METHODS: We used a three-phase method, including a two-round electronic Delphi survey and a workshop. The clinician and patient electronic surveys were conducted in parallel and independently. For phase I, the survey was distributed to patients from 27 European countries in 6 different languages, and to European Glaucoma Society members, ophthalmologists with expertise in glaucoma care, asking to name up to five research priorities. During phase II, participants were asked to rank the questions identified in phase I using a Likert scale. Phase III was a 1 day workshop with patients and clinicians. The purpose was to make decisions about the 10 most important research priorities using the top 20 priorities identified by patients and clinicians. RESULTS: In phase I, 308 patients and 150 clinicians were involved. In phase II, the highest-ranking priority for both patients and clinicians was 'treatments to restore vision'. In phase III, eight patients and four clinicians were involved. The top three priorities were 'treatments to stop sight loss', 'treatments to restore vision' and 'improved detection of worsening glaucoma'. CONCLUSION: We have developed a list of priorities for glaucoma research involving clinicians and patients from different European countries that will help guide research efforts and investment.

Spatial Summation in the Glaucomatous Macula: A Link With Retinal Ganglion Cell Damage.

Purpose: The purpose of this study was to test whether functional loss in the glaucomatous macula is characterized by an enlargement of Ricco's area (RA) through the application of a computational model linking retinal ganglion cell (RGC) damage to perimetric sensitivity. Methods: One eye from each of 29 visually healthy subjects <40 years old, 30 patients with glaucoma, and 20 age-similar controls was tested with a 10-2 grid with stimuli of 5 different area sizes. Structural estimates of point-wise RGC density were obtained from optical coherence tomography (OCT) scans. Structural and functional data from the young healthy cohort were used to estimate the parameters of a computational spatial summation model to generate a template. The template was fitted with a Bayesian hierarchical model to estimate the latent RGC density in patients with glaucoma and age-matched controls. We tested two alternative hypotheses: fitting the data by translating the template horizontally (H1: change in RA) or vertically (H2: loss of sensitivity without a change in RA). Root mean squared error (RMSE) of the model fits to perimetric sensitivity were compared. Ninety-five percent confidence intervals were bootstrapped. The dynamic range of the functional and structural RGC density estimates was denoted by their 1st and 99th percentiles. Results: The RMSE was 2.09 (95% CI = 1.92-2.26) under H1 and 2.49 (95% CI = 2.24-2.72) under H2 (P < 0.001). The average dynamic range for the structural RGC density estimates was only 11% that of the functional estimates. Conclusions: Macular sensitivity loss in glaucoma is better described by a model in which RA changes with RGC loss. Structural measurements have limited dynamic range.

Investigating the Spatiotemporal Summation of Perimetric Stimuli in Dry Age-Related Macular Degeneration.

Purpose: To measure achromatic spatial, temporal, and spatiotemporal summation in dry age-related macular degeneration (AMD) compared to healthy controls under conditions of photopic gaze-contingent perimetry. Methods: Twenty participants with dry AMD (mean age, 74.6 years) and 20 healthy controls (mean age, 67.8 years) performed custom, gaze-contingent perimetry tests. An area-modulation test generated localized estimates of Ricco's area (RA) at 2.5° and 5° eccentricities along the 0°, 90°, 180°, and 270° meridians. Contrast thresholds were measured at the same test locations for stimuli of six durations (3.7-190.4 ms) with a Goldmann III stimulus (GIII, 0.43°) and RA-scaled stimuli. The upper limit (critical duration) of complete temporal summation (using the GIII stimulus) and spatiotemporal summation (using the RA stimuli) was estimated using iterative two-phase regression analysis. Results: Median (interquartile range [IQR]) RA estimates were significantly larger in AMD participants (2.5°: 0.21 [0.09-0.41] deg2; 5°: 0.32 [0.15-0.65 deg2]) compared to healthy controls (2.5°: 0.08 [0.05-0.13] deg2; 5°: 0.15 [0.08-0.22] deg2) at all test locations (all P < 0.05). No significant difference in median critical duration was found in AMD participants with the GIII stimulus (19.6 [9.9-30.4] ms) and RA-scaled stimuli (22.9 [13.9-40.3] ms) compared to healthy controls (GIII: 17.0 [11.3-24.0] ms; RA-scaled: 22.4 [14.3-33.1] ms) at all test locations (all P > 0.05). Conclusions: Spatial summation is altered in dry AMD, without commensurate changes in temporal summation. Translational Relevance: The sensitivity of perimetry to AMD may be improved by utilizing stimuli that probe alterations in spatial summation in the disease.

Validating Trend-Based End Points for Neuroprotection Trials in Glaucoma.

Purpose: The purpose of this study was to evaluate the power of trend-based visual field (VF) progression end points against long-term development of event-based end points accepted by the US Food and Drug Administration (FDA). Methods: One eye from 3352 patients with ≥10 24-2 VFs (median = 11 years) follow-up were analyzed. Two FDA-compatible criteria were applied to these series to label "true-progressed" eyes: ≥5 locations changing from baseline by more than 7 dB (FDA-7) or by more than the expected test-retest variability (GPA-like) in 2 consecutive tests. Observed rates of progression (RoP) were used to simulate trial-like series (2 years) randomly assigned (1000 times) to a "placebo" or a "treatment" arm. We simulated neuroprotective "treatment" effects by changing the proportion of "true progressed" eyes in the two arms. Two trend-based methods for mean deviation (MD) were assessed: (1) linear mixed model (LMM), testing average difference in RoP between the two arms, and (2) time-to-progression (TTP), calculated by linear regression as time needed for MD to decline by predefined cutoffs from baseline. Power curves with 95% confidence intervals were calculated for trend and event-based methods on the simulated series. Results: The FDA-7 and GPA-like progression was achieved by 45% and 55% of the eyes in the clinical database. LMM and TTP had similar power, significantly superior to the event-based methods, none of which reached 80% power. All methods had a 5% false-positive rate. Conclusions: The trend-based methods can efficiently detect treatment effects defined by long-term FDA-compatible progression. Translational Relevance: The assessment of the power of trend-based methods to detect clinically relevant progression end points.

Improving the Accuracy and Speed of Visual Field Testing in Glaucoma With Structural Information and Deep Learning.

Purpose: To assess the performance of a perimetric strategy using structure-function predictions from a deep learning (DL) model. Methods: Visual field test-retest data from 146 eyes (75 patients) with glaucoma with (median [5th-95th percentile]) 10 [7, 10] tests per eye were used. Structure-function predictions were generated with a previously described DL model using cicumpapillary optical coherence tomography (OCT) scans. Structurally informed prior distributions were built grouping the observed measured sensitivities for each predicted value and recalculated for each subject with a leave-one-out approach. A zippy estimation by sequential testing (ZEST) strategy was used for the simulations (1000 per eye). Ground-truth sensitivities for each eye were the medians of the test-retest values. Two variations of ZEST were compared in terms of speed (average total number of presentations [NP] per eye) and accuracy (average mean absolute error [MAE] per eye), using either a combination of normal and abnormal thresholds (ZEST) or the calculated structural distributions (S-ZEST) as prior information. Two additional versions of these strategies employing spatial correlations were tested. Results: S-ZEST was significantly faster, with a mean average NP of 213.87 (SD = 28.18), than ZEST, with a mean average NP of 255.65 (SD = 50.27) (P < 0.001). The average MAE was smaller for S-ZEST (1.98; SD = 2.37) than ZEST (2.43; SD = 2.69) (P < 0.001). Spatial correlations further improved both strategies (P < 0.001), but the differences between ZEST and S-ZEST remained significant (P < 0.001). Conclusions: DL structure-function predictions can significantly improve perimetric tests. Translational Relevance: DL structure-function predictions from clinically available OCT scans can improve perimetry in glaucoma patients.

Screening for new psychoactive substances in wastewater from educational institutions.

Drug (ab)use among young people is a serious issue, negatively impacting their well-being and prospects. The emergence of new psychoactive substances (NPS) further complicates the situation as they are easily accessible (e.g., online), but users are at high risk of intoxication as their chemical identity is often unknown and toxicity poorly understood. While surveys and drug testing are traditionally used in educational institutions to comprehend drug use trends and establish effective prevention programs, they are not without their limitations. Accordingly, we investigated the occurrence of NPS in educational institutions through wastewater analysis and critically evaluated the viability of the approach. The study included eight wastewater samples from primary schools (ages 6-15 years), six from secondary schools (ages 15-19 years), three from institutions for both secondary and higher education (ages 15+), and six from higher educational institutions (ages 19+). Samples were obtained mid-week and evaluated in two Slovenian municipalities; the capital Ljubljana and a smaller one (M1). Samples were screened using liquid chromatography-ion mobility-high-resolution mass spectrometry (LC-IMS-HRMS), and NPS identified at three levels of confidence (Level 1: unequivocal, Level 2: probable, Level 3: tentative) from a suspect list containing over 5600 entries. NPS were identified in all types of educational institutions. Most were synthetic stimulants, with 3-MMC, ephedrine, 4-chloro-α-PPP, and ethcathinone being unequivocally identified. Also, NPS were present in wastewater from all educational institution types revealing potential spatial but no inter-institutional trends. Although specific groups cannot be targeted, the study, as a proof-of-concept, demonstrates that a suspect screening of wastewater employing LC-IMS-HRMS can be used as a radar for NPS in educational institutions and potentially replace invasive drug testing.

Retrospective spatiotemporal study of antidepressants in Slovenian wastewaters.

This study utilizes wastewater-based epidemiology (WBE) to evaluate spatiotemporal changes in the consumption of antidepressants before and during the COVID-19 pandemic in Slovenia. Composite 24-h influent wastewater samples (n = 210) were collected from six wastewater treatment plants between summer 2019 and spring 2021. The samples were extracted using 96-well solid-phase extraction and analysed by liquid chromatography-tandem mass spectrometry. The measured concentrations of target antidepressant biomarkers were then converted to population-normalised mass loads (PNMLs), taking into account flow rate and catchment population. Ten biomarkers, including amitriptyline, bupropion, bupropion-OH, citalopram, norcitalopram, normirtazapine, venlafaxine, O-desmethylvenlafaxine, trazodone, and moclobemide, were above the lower limit of quantification and were included in the spatiotemporal temporal assessment. The highest PNMLs were detected for O-desmethylvenlafaxine (mean ± SD: 82.1 ± 21.2 mg/day/1000 inhabitants) and venlafaxine (38.0 ± 10.6 mg/day/1000 inhabitants), followed by citalopram (27.0 ± 10.7 mg/day/1000 inhabitants). In addition, the mean metabolite/parent compound ratios were comparable with other WBE studies indicating consumption rather than direct disposal. Overall, the results indicated significant spatiotemporal variations depending on the location, and the PNMLs of most biomarkers increased during the first wave of the COVID-19 pandemic (spring of 2020). However, no clear spatial patterns were revealed related to the pandemic.

The Role of Mitophagy in Glaucomatous Neurodegeneration.

This review aims to provide a better understanding of the emerging role of mitophagy in glaucomatous neurodegeneration, which is the primary cause of irreversible blindness worldwide. Increasing evidence from genetic and other experimental studies suggests that mitophagy-related genes are implicated in the pathogenesis of glaucoma in various populations. The association between polymorphisms in these genes and increased risk of glaucoma is presented. Reduction in intraocular pressure (IOP) is currently the only modifiable risk factor for glaucoma, while clinical trials highlight the inadequacy of IOP-lowering therapeutic approaches to prevent sight loss in many glaucoma patients. Mitochondrial dysfunction is thought to increase the susceptibility of retinal ganglion cells (RGCs) to other risk factors and is implicated in glaucomatous degeneration. Mitophagy holds a vital role in mitochondrial quality control processes, and the current review explores the mitophagy-related pathways which may be linked to glaucoma and their therapeutic potential.

Socioeconomic Status Affects Postoperative Time to Union in Pediatric Patients with a Surgically Treated Fracture.

Fractures account for 10% to 25% of all pediatric injuries, and surgical treatment is common. In such cases, postoperative healing can be affected by a number of factors, including those related to socioeconomic status (SES). The purpose of this study was to investigate the relationship between time to fracture union and SES, which was measured with use of the median household income (MHI) and Child Opportunity Index (COI). Methods: A retrospective review was conducted of pediatric patients with a long-bone fracture that had been surgically treated at a Level-I pediatric trauma center between January 2010 and June 2020. Demographic and relevant medical data were collected. Patients were sorted into union and nonunion groups. The ZIP code of each patient was collected and the MHI and COI of that ZIP code were identified. Income brackets were created in increments of $10,000 ranging from $20,000 to $100,000, with an additional category of >$100,000, and patients were sorted into these groups according to MHI. Comparisons among the income groups and among the union status groups were conducted for each of the collected variables. A multiple regression analysis was utilized to determine the independent effect of each variable on time to union. Results: A total of 395 patients were included in the final sample, of whom 51% identified as Hispanic. Patients in the union group had a higher mean COI and MHI. Nonunion occurred in only 8 patients. Patients who achieved fracture union in ≤4 months had a significantly higher mean COI and MHI. When controlling for other demographic variables, the time to union increased by a mean of 9.6 days for every $10,000 decrease in MHI and increased by a mean of 6.8 days for every 10-unit decrease in the COI. Conclusions: The present study is the first, to our knowledge, to investigate the relationship between SES and time to fracture union in pediatric patients. When controlling for other demographic factors, we found a significant relationship between SES and time to union in pediatric patients with a surgically treated fracture. Further investigations of the relationship between SES and time to union in pediatric patients are needed to determine potential mechanisms for this relationship. Level of Evidence: Therapeutic Level III. See Instructions for Authors for a complete description of levels of evidence.

Occurrence and sources of residues of drugs of abuse in an urban aquifer: Chemical analysis and solute transport modelling.

This study investigated the occurrence and potential sources of residues of drugs of abuse in an urban aquifer beneath the City of Ljubljana using water analysis and a solute transport model designed to predict nitrogen distribution. Samples were collected from three sources: 28 wastewater samples (24-h composites), 4 aquifer-recharging river samples (grab), and 22 groundwater samples. The samples were analysed for residues of commonly (ab)used licit drugs (nicotine and alcohol), medications of abuse (morphine, methadone, codeine, and ketamine), and illicit drugs (tetrahydrocannabinol - THC, cocaine, amphetamines, and heroin) using liquid-liquid (alcohol residue) and solid-phase extraction, followed by liquid chromatography-tandem mass spectrometry (LC-MS/MS). Additionally, we used solute transport modelling to predict the spatial distribution of drug residues in the aquifer and their potential sources. Nicotine (up to 45,7 ng/L), cotinine (up to 5.86 ng/L), trans-3'-hydroxycotinine (up to 0.528 ng/L) and benzoylecgonine (up to 0.572 ng/L) were the most commonly detected drug residues in groundwater, followed by cocaine (